Simulated Indonesian Throughflow in Makassar Strait across the SODA3 products

The Indonesian Throughflow(ITF), which connects the tropical Pacific and Indian oceans, plays important roles in the inter-ocean water exchange and regional or even global climate variability. The Makassar Strait is the main inflow passage of the ITF, carrying about 77% of the total ITF volume transport. In this study, we analyze the simulated ITF in the Makassar Strait in the Simple Ocean Data Assimilation version 3(SODA3) datasets. A total of nine ensemble members of the SODA3 datasets, of which are driven by different surface forcings and bulk formulas, and with or without data assimilation, are used in this study. The annual mean water transports(i.e.,volume, heat and freshwater) are related to the combination of surface forcing and bulk formula, as well as whether data assimilation is employed. The phases of the seasonal and interannual variability in water transports cross the Makassar Strait, are basically consistent with each other among the SODA3 ensemble members. The interannual variability in Makassar Strait volume and heat transports are significantly correlated with El Ni?oSouthern Oscillation(ENSO) at time lags of-6 to 7 months. There is no statistically significant correlation between the freshwater transport and the ENSO. The Makassar Strait water transports are not significantly correlated with the Indian Ocean Dipole(IOD), which may attribute to model deficiency in simulating the propagation of semiannual Kelvin waves from the Indian Ocean to the Makassar Strait.

Subthermocline eddies carrying the Indonesian Throughflow water observed in the southeastern tropical Indian Ocean

We observed a subthermocline eddy(STE)with a cold and fresh core during an observation cruise along a transect of 10°S in the southeastern tropical Indian Ocean(SETIO)in December 2017.The vertical scale,speed radius,and maximum swirl velocity of the STE were about 200 m,55 km,and 0.5 m/s,respectively.The mean Rossby number and Burger number of the STE were then estimated to be about−0.7 and 2.4,indicating the STE was a submesoscale coherent vortex.The STE core water had characteristics of the Indonesian Throughflow(ITF)water and was distinct from that of surrounding areas.By examining Argo float data,another STE was well captured by five successive profiles of the same Argo float.Both STEs showed significant temperature and salinity anomalies at theσ0=26.0-26.5 kg/m3 surfaces.With the assumption that the low-salinity ITF water parcels could be carried only by surface eddies and the STEs,the Argo profiles,which detected low-salinity ITF water and were located outside a surface eddy,were believed to be inside an STE and were used to analyze the distribution,origin,and generation mechanism of the STE.The results suggested that the STEs carrying ITF water may be generated under topography-current interaction at the eastern coastal waters or under front-induced subduction in the area away from coastal waters.Those STEs may be widely distributed in the SETIO and may play a role in ITF water parcel transport.

An overview of 10-year observation of the South China Sea branch of the Pacific to Indian Ocean throughflow at the Karimata Strait

Besides the Indonesian throughflow(ITF), the South China Sea throughflow(SCSTF) also contributes to the water transport from the Pacific to the Indian Ocean. However, this South China Sea(SCS) branch at the Karimata Strait is poorly observed until 2007, even though its importance has been suggested by numerical studies for decades. In this paper, we review the nearly 10-year field measurement in the Karimata Strait by the execution of the projects of "SCS-Indonesian Seas Transport/Exchange(SITE) and Impacts on Seasonal Fish Migration" and "The Transport, Internal Waves and Mixing in the Indonesian Throughflow regions(TIMIT) and Impacts on Marine Ecosystem", which extend the observations from the western Indonesian seas to the east to include the main channels of the ITF, is introduced. Some major achievements from these projects are summarized.

Climatology and Variability of the Indonesian Throughflow in an Eddy-permitting Oceanic GCM

A quasi-global eddy permitting oceanic GCM, LICOM1.0, is run with the forcing of ERA40 daily wind stress from 1958 to 2001. The modelled Indonesian Throughflow (ITF) is reasonable in the aspects of both its water source and major pathways. Compared with the observation, the simulated annual mean and seasonal cycle of the ITF transport are fairly realistic. The interannual variation of the tropical Pacific Ocean plays a more important role in the interannual variability of the ITF transport. The relationship between the ITF and the Indian Ocean Dipole (IOD) also reflects the influence of ENSO. However, the relationship between the ITF transport and the interannual anomalies in the Pacific and Indian Oceans vary with time. During some years, (e.g., 1994), the effect of a strong IOD on the ITF transport is more than that from ENSO.

Covariation of the Indonesian Throughflow and South China Sea Throughflow Associated with the 1976/77 Regime Shift

Changes in the Indonesian Throughflow （ITF） and the South China Sea throughflow-measured by the Luzon Strait Transport （LST）-associated with the 1976/77 regime shift are analyzed using the Island Rule theory and the Simple Ocean Data Assimilation dataset. Results show that LST increased but ITF transport decreased after 1975. Such changes were induced by variations in wind stress associated with the regime shift. The strengthening of the easterly wind anomaly east of the Luzon Strait played an important role in the increase of LST after 1975, while the westerly wind anomaly in the equatorial Pacific contributed significantly to the decrease in ITF transport after 1975; accounting for 53% of the change. After 1975, the Kuroshio Current strengthened and the Mindanao Current weakened in response to a decrease in the total transport of the North Equatorial Current. Both the North Equatorial Countercurrent and the South Equatorial Current weakened after 1975, and an anomalous cyclonic circulation in the western equatorial Pacific prevented the tropical Pacific water from entering the Indian Ocean directly.

Influence of Positive and Negative Indian Ocean Dipoles on ENSO via the Indonesian Throughflow: Results from Sensitivity Experiments

The role of the Indonesian Throughflow （ITF） in the influence of the Indian Ocean Dipole （IOD） on ENSO is investigated using version 2 of the Parallel Ocean Program （POP2） ocean general circulation model. We demonstrate the results through sensitivity experiments on both positive and negative IOD events from observations and coupled general circulation model simulations. By shutting down the atmospheric bridge while maintaining the tropical oceanic channel, the IOD forcing is shown to influence the ENSO event in the following year, and the role of the ITF is emphasized. During positive IOD events, negative sea surface height anomalies （SSHAs） occur in the eastern Indian Ocean, indicating the existence of upwelling. These upwelling anomalies pass through the Indonesian seas and enter the western tropical Pacific, resulting in cold anomalies there. These cold temperature anomalies further propagate to the eastern equatorial Pacific, and ultimately induce a La Nifia- like mode in the following year. In contrast, during negative IOD events, positive SSHAs are established in the eastern Indian Ocean, leading to downwelling anomalies that can also propagate into the subsurface of the western Pacific Ocean and travel further eastward. These downwelling anomalies induce negative ITF transport anomalies, and an E1 Nifio-like mode in the tropical eastern Pacific Ocean that persists into the following year. The effects of negative and positive IOD events on ENSO via the ITF are symmetric. Finally, we also estimate the contribution of IOD forcing in explaining the Pacific variability associated with ENSO via ITE

Impact of the South China Sea Throughflow on the Pacific Low-Latitude Western Boundary Current:A Numerical Study for Seasonal and Interannual Time Scales

Prior studies have revealed that,as a part of the Pacific tropical gyre,the South China Sea throughflow（SCSTF） is strongly influenced by the Pacific low-latitude western boundary current（LLWBC）.In this study,ocean general circulation model（OGCM） experiments with and without connection to the South China Sea（SCS） were performed to investigate the impact of the SCSTF on the Pacific LLWBC.These model experiments show that if the SCS is blocked,seasonal variability of the Kuroshio and Mindanao Current becomes stronger,and the meridional migration of the North Equatorial Current（NEC） bifurcation latitude is enhanced.Both in seasonal and interannual time scales,stronger Luzon Strait transport（LST） induces a stronger Kuroshio transport combined with a southward shift of the NEC bifurcation,which is unfavorable for a further increase of the LST;a weaker LST induces a weaker Kuroshio transport and a northward shifting NEC bifurcation,which is also unfavorable for the continuous decrease of the LST.

Implication of the South China Sea Throughflow for the Interannual Variability of the Regional Upper-Ocean Heat Content

In this study the interannual variability of the upper-ocean heat content in the South China Sea （SCS） was revisited using simple ocean data assimilation （SODA） combined with objective analyzed data sets that included the horizontal and vertical structures. The results confirmed that the upper-ocean heat content in the SCS is lower than normal during the mature phase of E1 Nifio events, and two super E1 Nifio events, 1982/1983 and 1997/1998 were also included. The variability of the heat content was consistent with the variability of the dynamic height anomalies. The SCS throughflow （SCSTF） plays an important role in regulating the interannual variability of the heat content, especially during the mature phase of E1 Nifio events.

Origins and pathways of the subsurface and intermediate water masses of the Indonesian Throughflow derived from historical and Argo data

On the basis of Argo data and historic temperature/salinity data from the World Ocean Database 2001 （ WOD01 ）, origins and spreading pathways of the subsurface and intermediate water masses in the Indonesian Throughflow （ITF） region were discussed by analyzing distributions of salinity on representative isopyenal layers. Results were shown that, subsurface water mostly comes from the North Pacific Ocean while the intermediate water originates from both the North and South Pacific Ocean, even possibly from the Indian Ocean. Spreading through the Sulawesi Sea, the Makassar Strait, and file Flores Sea, the North Pacific subsurface water and the North Pacific Intermediate water dominate the western part of the Indonesian Archipelago. Furthermore as the depth increases, the features of the North Pacific sourced water masses become more obvious. In the eastern part of the waters, high sa- linity South Pacific subsurface water is blocked by a strong salinity front between Halmahera and New Guinea. Intermediate water in the eastern interior region owns salinity higher than the North Pacific intermediate water and the antarctic intermediate water （ AAIW）, possibly coming from the vertical mixing between subsurface water and the AAIW from the Pacific Ocean, and possibly coming from the northward extending of the AAIW from the Indian Ocean as well.

Seasonal and Interannual Variability of Bering Strait Throughflow from AO-FVCOM and Observation

A high-resolution Arctic Ocean-Finite Volume Community Ocean Model(AO-FVCOM) and observational current data from 14 mooring stations in Bering Strait and surrounding regions between 1990 and 2015 were used to study the seasonal and interannual variability of Bering Strait throughflow(BST). AO-FVCOM represented the BST with a climatological northward flux of 1.06 Sv, which was close to the observational mean of 0.94 ± 0.26 Sv. From the model results, the strongest volume flux was in summer, approximately 45% larger than that in winter. Interannual variability of BST was also indicated in the model results, and the maximum and minimum annual mean transports are in 2007 and 2012, respectively. AO-FVCOM showed larger differences from the observations in 2000, 2002, and 2015 than in other years, which may be related to the limitation of atmospheric forcing for the model. According to the driving mechanisms of BST, sea level difference(SLD) across the strait dominates the northward volume transport, and local wind is also important in forcing the seasonal variability of the BST and SLD patterns to change the BST indirectly.

Seasonal and interannual variability of water mass sources of Indonesian throughflow in the Maluku Sea and the Halmahera Sea

So far, large uncertainties of the Indonesian throughflow(ITF) reside in the eastern Indonesian seas, such as the Maluku Sea and the Halmahera Sea. In this study, the water sources of the Maluku Sea and the Halmahera Sea are diagnosed at seasonal and interannual timescales and at different vertical layers, using the state-of-the-art simulations of the Ocean General Circulation Model(OGCM) for Earth Simulator(OFES). Asian monsoon leaves clear seasonal footprints on the eastern Indonesian seas. Consequently, the subsurface waters(around 24.5σ_θ and at ～150 m) in both the Maluku Sea and the Halmahera Sea stem from the South Pacific(SP) during winter monsoon, but during summer monsoon the Maluku Sea is from the North Pacific(NP), and the Halmahera Sea is a mixture of waters originating from the NP and the SP. The monsoon impact decreases with depth, so that in the Maluku Sea, the intermediate water(around 26.8σ_θ and at ～480 m) is always from the northern Banda Sea and the Halmahera Sea water is mainly from the SP in winter and the Banda Sea in summer. The deep waters(around27.2σ_θ and at ～1 040 m) in both seas are from the SP, with weak seasonal variability. At the interannual timescale,the subsurface water in the Maluku Sea originates from the NP/SP during El Ni?o/La Ni?a, while the subsurface water in the Halmahera Sea always originates from the SP. Similar to the seasonal variability, the intermediate water in Maluku Sea mainly comes from the Banda Sea and the Halmahera Sea always originates from the SP. The deep waters in both seas are from the SP. Our findings are helpful for drawing a comprehensive picture of the water properties in the Indonesian seas and will contribute to a better understanding of the ocean-atmosphere interaction over the maritime continent.

Heat contribution of the Indonesian throughflow to the Indian Ocean

Based on the high-resolution Eulerian fields of an ocean general circulation model simulation, the heat contribution of the Indonesian throughflow(ITF) to the Indian Ocean is estimated by Lagrangian tracing method.The heat transport of each particle of ITF waters is calculated by tracing temperature change along the trajectory until the particle exits the Indian Ocean. The simulation reveals that the ITF waters flow westward and branch near Madagascar, further showing the ITF waters are redistributed in both northern and southern Indian Ocean.Heat budget analysis indicates that the ITF waters gain 0.41 PW(Petawatts, 1015 W) in the northern Indian Ocean and lose 0.56 PW in the southern Indian Ocean, respectively. As a result, the ITF waters warm the whole Indian Ocean basin with only 0.15 PW, which shows an "insignificant" role of ITF on the Indian Ocean because of the heat exchange compensation between northern and southern Indian Ocean. Furthermore, the tracing pathways show that the ITF waters mainly flow out the Indian Ocean at both sides of the basin via Agulhas Current and Leeuwin Current. About 89% of the ITF waters leave along western boundary and the rest 11% along eastern boundary. Compared to seeding section, 0.10 PW and 0.05 PW are released to the Indian Ocean, respectively.

Stable carbon and oxygen isotopes of four planktonic foraminiferal species from core-top sediments of the Indonesian throughflow region and their significance

Horizontal and vertical distributions of δ^18 and δ^13 were investigated in shells of four planktonic foraminiferal species, Globigerinoides ruber, Globigerinoides sacculifer, Pulleniatina obliquiloculata and Neogloboquedrina dutertrei, from a total of 62 core-top sediment samples from the Indonesian throughflow region. Results were compared to modern hydrologic conditions in order to explore potential of proxies in reconstructing fluvial discharge and upper ocean water column characteristics in this region. Our results show that, in the Makassar Strait, both of depleted δ^18 and δ^13 of these four species were linked to freshwater input. In the Bali Sea,however, depleted δ^18 and δ^13 for these species may be due to different reasons. Depleted δ^18 was a result of freshwater input and as well influenced by along-shore currents while depleted δ^13 was more likely due to the Java-Sumatra upwelling. Comparison of shell δ^18 records and hydrographic data of World Ocean Atlas 2005 suggests that G. ruber and G. sacculifer calcify within the mixed-layer, respectively at 0–50 m and 20–75 m water depth, and P. obliquiloculata and N. dutertrei within the upper thermocline, both at 75–125 m water depth. N.dutertrei calcifies at slightly deeper water depth than P. obliquiloculata does. In general, δ^13 values of both G.ruber and G. sacculifer are larger than those of P. obliquiloculata and N. dutertrei at all sites, possibly related to depth habitats of these species and vertical distribution of nutrients in the Indonesian throughflow region.

Correlation analysis of the North Equatorial Current bifurcation and the Indonesian Throughflow

Based on monthly mean Simple Ocean Data Assimilation （SODA） products from 1958 to 2007, this study analyzes the seasonal and interannual variability of the North Equatorial Current （NEC） bifurcation latitude and the Indonesian Throughflow （ITF） volume transport. Further, Empirical Mode Decomposition （EMD） method and lag-correlation analysis are employed to reveal the relationships between the NEC bifurcation location, NEC and ITF volume transport and ENSO events. The analysis results of the seasonal variability show that the annual mean location of NEC bifurcation in upper layer occurs at 14.33°N and ITF volume transport has a maximum value in summer, a minimum value in winter and an annual mean transport of 7.75×10^6 m^3/s. The interannual variability analysis indicates that the variability of NEC bifurcation location can be treated as a precursor of El Nino. The correlation coefficient between the two reaches the maximum of 0.53 with a time lag of 2 months. The ITF volume transport is positively related with E1 Nifio events with a maximum coefficient of 0.60 by 3 months. The NEC bifurcation location is positively correlated with the ITF volume transport with a correlation coefficient of 0.43.

Throughflow and g-jitter effects on binary fluid saturated porous medium

A non-autonomous complex Ginzburg-Landau equation （CGLE） for the finite amplitude of convection is derived, and a method is presented here to determine the amplitude of this convection with a weakly nonlinear thermal instability for an oscillatory mode under throughflow and gravity modulation. Only infinitesimal disturbances are considered. The disturbances in velocity, temperature, and solutal fields are treated by a perturbation expansion in powers of the amplitude of the applied gravity field. Throughfiow can stabilize or destabilize the system for stress free and isothermal boundary conditions. The Nusselt and Sherwood numbers are obtained numerically to present the results of heat and mass transfer. It is found that throughfiow and gravity modulation can be used alternately to heat and mass transfer. Further, oscillatory flow, rather than stationary flow, enhances heat and mass transfer.

Contributions of the Bering Strait throughflow to oceanic meridional heat transport under modern and Last Glacial Maximum climate conditions

Paleo reconstructions and model simulations have suggested the Bering Strait plays a pivotal role in climate change. However, the contribution of the Bering Strait throughflow to oceanic meridional heat transport (OMHT) is about 100 times smaller than the OMHT at low latitudes in the modern climate and it is generally ignored. Based on model simulations under modern and Last Glacial Maximum (LGM,~21 ka;ka=thousand years ago) climate conditions, this study highlights the importance of the Bering Strait throughflow to OMHT. The interbasin OMHT induced by the Bering Strait throughflow is estimated by interbasin-intrabasin decomposition. Similar to barotropic-baroclinic-horizontal decomposition, we assume the nonzero net mass transport induced by interbasin throughflows is uniform across the entire section, and the interbasin term is separated to force zero net mass transport for the intrabasin term. Based on interbasinintrabasin decomposition, the contribution of the Bering Strait throughflow is determined as ~0.02 PW (1 PW=10 15 W) under the modern climate, and zero under the LGM climate because the closed Bering Strait blocked interbasin throughflows. The contribution of the Bering Strait throughflow to OMHT is rather small, consistent with previous studies. However, comparisons of OMHT under modern and LGM climate conditions indicate the mean absolute changes are typically 0.05 and 0.20 PWin the North Atlantic and North Pacific, respectively. Thus, the contribution of the Bering Strait throughflow should not be ignored when comparing OMHT under diff erent climate conditions.

A METHOD TO ESTIMATE DYNAMIC HEIGHT AND GEOSTROPHIC TRANSPORT RELATIVE TO DEEPER LEVELS FROM 400M EXPENDABLE BATHYTHERMOGRAPH IN THE INDONESIAN THROUGHFLOW

The Expendable Bathythermograph (XBT) Programme used a mix of T4 (450m) and T7(750 m) XBT’s during the pre-TOGA periods. Studies are needed to determine how to use the T4/T7 datatogether, in particular with regard to a reference level for calculation of dynamic height and geostrophiccurrents. Temperature profiles to 750 m collected from 1986 through 1989 on the trackline across theIndonesian throughflow between NW Australia and Java are used to show the relations between dynamicbeight and geostriohic flow using reference levels at 400 db and 750 db. A very high temporalcorrelation between vertically averaged temperture in the upper 400 m and dynamic height at 50 m rela-tive to 750 db was found. The corresponding regression relationships are presented for all one degree lati-tude bins along the section and can be used for dynamical calculation of currents in the upper 400 m rel-ative to 750 db .An attempt is made to estimate volume transport relative to 750 db from 400 m pro-files. Problems which make

The IOD signal in the upper-ocean Indonesian Throughflow since the mid-1970s

On the basis of simple ocean data assimilation （SODA） reanalysis product, the interannual variability of upper-ocean Indonesian Throughflow （ITF） volume transport since the mid 1970s is examed. The wavelet analysis shows a second prominent interannual oscillation with a period of about 2 ～ 4 a. To reveal any relationship between this band-scale oscillation of upper-ocean ITF and the Indian Ocean dipole （IOD）, the correlation and wavelet analyses are used. The correlation coefficient between the upper-ocean ITF and the IOD reaches - 0.40 with upper-ocean ITF lagging an IOD index by eight months. The wavelet power spectrum of upper-ocean ITF shows similar structure to that of the IOD index. And the evolution of IOD is reproduced by lagged correlation between the upper-ocean ITF and the sea surface temperature anomaly （SSTA） over the Indian Ocean. It suggests that the 2 ～ 4 a band-scale oscillation of upper-ocean ITF is related uniquely to the IOD over the tropical Indian Ocean.

Geochemical characteristics from tests of four modern planktonic foraminiferal species in the Indonesian Throughflow region and their implications

Test geochemistry of planktonic foraminifera is an indispensable tool in reconstructing past ocean hydrological changes. It is essential to investigate region-specific implications of test geochemistry,although those established from other regions can be broadly applied. In this study, characteristics of6180 and Mg/Ca from tests of four planktonic foraminiferal species, Globigerinoides ruber sensu stricto(s.s.), Globigerinoides sacculifer, Pulleniatina obliquiloculata and Neogloboquadrina dutertrei, from 60 coretop sediment samples retrieved from the Indonesian Throughflow(ITF) region were studied. These geochemical data were compared with modern hydrographic profiles in order to assess their relations and to investigate potential implications of test geochemical parameters in reconstructing past oceanographic change in the ITF region. Calcification depths of these four species were first estimated based on comparison between measured test δ180 and predicted calcite δ^(18)O that was calculated from modern temperature and salinity. The results indicate that G. ruber s.s. and G. sacculifer calcify within the mixedlayer at 0-50 m and 20-75 m, respectively, whereas P. obliquiloculata and N. dutertrei calcify within the thermocline at around 75 to 125 m. A combined study of excess Mg/Ca(difference between measured and predicted Mg/Ca) and salinity suggests that salinity exerts a negligible impact on test Mg/Ca of these foraminiferal species in the ITF region. Comparison of test Mg/Ca-derived temperatures with temperature profiles of the upper 200 m of the water column from the seas of the ITF region also indicate calcification depths of these species, which match well with the above estimations using test δ^(18)O. It further indicates that G. sacculifer may be more sensitive in reflecting changes in the depth of the mixedlayer, highlighting a potential use of Mg/Ca temperature difference between G. ruber s.s. and G. sacculifer in reconstructing the depth of the mixed-layer in the ITF region.

The South China Sea throughflow: linkage with local monsoon system and impact on upper thermal structure of the ocean

The tendency of South China Sea throughflow (SCSTF) variation associated with the local monsoon system, and its impact on upper-layer thermal structure, are studied using the Simple Ocean Data Assimilation (SODA) dataset, combined with Ishii reanalysis data. Luzon Strait Transport (LST) is measured and used as an index for studying the SCSTF variation. Results show that LST had an increasing tendency over the last 50 years, mainly in summer and fall. The increasing tendency was 0.017 1 Sv/a in summer and 0.027 4 Sv/a in fall, as estimated by SODA, and 0.018 0Sv/a in summer and 0.018 9 Sv/a in fall, as estimated by "Island Rule" theory. LST increased by 0.53Sv in JJA (June-July-August) and 0.98Sv in SON (September-October-November) after climate shift, as inferred by SODA data. The average LST anomaly in JJA and SON is strongly related to the local monsoon system, especially to variability of the meridional wind stress anomaly after application of a 3-year running mean, with correlation coefficients 0.57 and 0.51, respectively. In addition to the basin-scale wind forcing, the local northeasterly wind stress anomaly in the SCS can push Pacific water entering the SCS more readily in JJA and SON after climate shift, and an SCSTF-associated cooling effect may favor subsurface cooling more frequently after climate shift.