Search

Deep water in the South China Sea is renewed by the cold and dense Luzon Strait overflow. However, from where and how the deep water upwells is poorly understood yet. Based on the Hybrid Coordinate Ocean Model reanalysis data, vertical velocity is derived to answer these questions. Domain-integrated vertical velocity is of two maxima, one in the shallow water and the other at depth, and separated by a layer of minimum at the bottom of the thermocline. Further analysis shows that this two-segmented vertical transport is attributed to the vertical compensation of subsurface water to the excessive outflow of shallow water and upward push of the dense Luzon Strait overflow, respectively. In the abyssal basin, the vertical transport increases upward from zero at the depth of 3 500–4 000 m and reaches a maximum of 1.5×106 m3/s at about 1 500 m. Deep water upwells mainly from the northeastern and southwestern ends of the abyssal basin and off the continental slopes. To explain the upward velocity arising from slope breaks, a possible mechanism is proposed that an onshore velocity component can be derived from the deep western boundary current above steep slopes under bottom friction.

This study aims to investigate variability of the deep South China Sea (SCS) circulation using the Hybrid Coordinate Ocean Model (HYCOM) global reanalysis product. The results reveal that annual cycle is a dominant component in the deep SCS circulation. Meanwhile, the boundary circulation strength is the weakest in January and peaks between June and September. The eastern and southern boundary currents strengthen and weaken one to three months earlier than that of the western and northern boundaries. Vector Empirical Orthogonal Functions (VEOF) analysis results reveal that semiannual and intraseasonal fluctuations are significant components, of which the spatial patterns are mainly confined in the northern and western boundary areas as well as the southwestern sub-basin. Wavelet analysis results show the strength of significant fluctuation varies year to year. Trend analysis results indicate a decadal weakening in the deep SCS circulation. An anomalous anticyclonic circulation, 50–70 km apart from the slope break, tends to weaken the cyclonic boundary circulation in the western and northern boundaries as well as the southwestern sub-basin. This trend is similar to the observed decadal weakening in the Northern Atlantic deep circulation. Thus, the findings of this study reveal that the variability of the deep SCS circulation has a remarkable response to the climate change. The mechanisms responsible for the variability are worth pursuing if more observations are available.

In this study, the Indian Ocean Subtropical Underwater (IOSTUW) was investigated as a subsurface salinity maximum using Argo floats (2000–20) for the first time. It has mean salinity, potential temperature, and potential density values of 35.54 ± 0.29 psu, 17.91° ± 1.66°C, and 25.56 ± 0.35 kg m−3, respectively, and mainly extends between 10° and 30°S along the isopycnal surface in the subtropical south Indian Ocean. The annual subduction rate of the IOSTUW during the period of 2004–19 was investigated based on a gridded Argo dataset. The results revealed a mean value of 4.39 Sv (1 Sv ≡ 106 m3 s−1) with an interannual variability that is closely related to the southern annular mode (SAM). The variation in the annual subduction rate of the IOSTUW is dominated by the lateral induction term, which largely depends on the winter mixed layer depth (MLD) in the sea surface salinity (SSS) maximum region. The anomalies of winter MLD are primarily determined by SAM-related air–sea heat flux and zonal wind anomalies through modulation of the buoyancy. As a result, the annual subduction rate of the IOSTUW generally increased when the SAM index showed negative anomalies and decreased when the SAM index showed positive anomalies. Exceptional cases occurred when the wind anomaly within the SSS maximum region was weak or was dominated by its meridional component. Significance Statement The Indian Ocean Subtropical Underwater (IOSTUW) is characterized as a subsurface salinity maximum in the subtropical south Indian Ocean, and it has potential contribution to regional sea level change and interoceanic salinity exchange. The purpose of this work is to give a more precise description of the IOSTUW and to better understand the formation process of this water mass. Our results shine new light on hydrological characteristics of the IOSTUW based on an observational dataset collected during the past two decades. In addition, the formation of the IOSTUW was found to be closely related to the leading mode of variability in the atmospheric circulation of the Southern Hemisphere though air–sea heat flux and surface wind anomaly.

The mid-depth ocean circulation is critically linked to actual changes in the long-term global climate system. However, in the past few decades, predictions based on ocean circulation models highlight the lack of data, knowledge, and long-term implications in climate change assessment. Here, using 842,421 observations produced by Argo floats from 2001-2020, and Lagrangian simulations, we show that only 3.8% of the mid-depth oceans, including part of the equatorial Pacific Ocean and the Antarctic Circumpolar Current, can be regarded as accurately modelled, while other regions exhibit significant underestimations in mean current velocity. Knowledge of ocean circulation is generally more complete in the low-latitude oceans but is especially poor in high latitude regions. Accordingly, we propose improvements in forecasting, model representation of stochasticity, and enhancement of observations of ocean currents. The study demonstrates that knowledge and model representations of global circulation are substantially compromised by inaccuracies of significant magnitude and direction, with important implications for modelled predictions of currents, temperature, carbon dioxide sequestration, and sea-level rise trends.

The cold wake caused by a tropical cyclone (TC) extends for hundreds of kilometers and persists for several weeks, thus influencing the surface response for any subsequent TCs that might pass over it. It is commonly accepted that sea-surface temperature (SST) cooling, as produced by a single TC, occurs primarily through vertical mixing. However, when there are sequential TCs, the earlier TC can dramatically change the thermal structure of the upper ocean, which may influence the subsequent development of a latter-occurring TC (LTC). Therefore, the contribution of horizontal advection and vertical mixing to SST-cooling during the passage of LTCs is of great interest. Using a 19-year-long observational dataset and the heat budget analysis of an idealized numerical simulation, the SST change during the passage of sequential TCs is investigated. The results demonstrate that, on average, the SST cooling caused by the LTC shows an overall decreasing trend with enhanced lingering wakes. Budget analysis of the model simulations suggests that an earlier TC can suppress the vertical mixing induced by an LTC mainly through an alteration of dynamics within the deepened mixed layer and that the contribution of vertical mixing to the SST cooling is weaker due to the intensification of the earlier TC. The weakened vertical mixing dominates the decreased SST cooling induced by an LTC. In contrast, the cold wake generated by an earlier TC can produce more cold water on the right side of the TC’s track, which contributes to stronger horizontal advection upon the arrival of the LTC. In general, the effects of the earlier TC can suppress the sea-surface thermal response to an LTC. If the contribution of the horizontal advection to SST cooling is neglected, the SST cooling induced by an LTC could be reduced by about 40%. As for the response of the sub-surface water to the passage of an LTC, the weakened warm anomaly induced by vertical mixing and the enhanced cooling anomaly caused by the vertical advection explain the reduced tendency for the mixed layer to deepen. As a result, the tendency for the mixed layer depth (MLD) to increase is suppressed during the passage of an LTC. These results highlight the importance of optimally depicting cold wakes in numerical simulations to improve the prediction of the upper ocean’s response to sequential TCs.

Exact knowledge on the seasonal variations of main tidal constituents is beneficial for improving tidal prediction. The semi-annual cycles in K1 and S2 tides are abnormally exaggerated by astronomical P1 and K2 tides, which interferes with our understanding on tidal seasonality. The widely-used tidal inference method in previous studies cannot fully separate astronomical P1 and K2 tides from seasonal P1 and K2 tides due to inaccurate inference relationship. In this study, on the basis of the ‘credo of smoothness’ which indicates that tidal admittances are smooth functions of tidal frequencies, we develop a novel but simple method to address this intractable issue and applied this method to explore the seasonality of tidal currents observed in the deep Timor Passage at the depth of 1800m. We find that the timing and range of seasonal modulations of M2, S2, K1, and O1 tides are distinct. Annual variations in tidal currents are much stronger than semi-annual variations in tidal currents. The annual and semi-annual ranges of M2 tide can reach 2.69 cm/s and 1.51 cm/s, which are largest among main constituents. Although the annual range of K1 tide is only 1.85 cm/s, considering the relatively small amplitude of time-averaged K1 tide (2.87cm/s), K1 the most affected tide by the annual cycle. The seasonal cycles of semi-diurnal tides (M2 and S2) are basically synchronous while those of diurnal tides (K1 and O1) are generally out-of-phase. As a general method, the proposed method can be widely applied to other sea areas to explore local tidal seasonality.

The Karimata Strait (KS) throughflow between the South China Sea (SCS) and Java Sea plays an essential role in heat and freshwater budget in the SCS and dual roles in strengthening/reducing the primary Indonesian throughflow (ITF) in the Makassar Strait. A sustained long-term monitoring of the ITF is logistically challenging and expensive; therefore, proxies are needed. Here, we use a combination of in situ measurement of the KS throughflow and satellite-derived sea surface height (SSH) and sea surface wind (SSW) to determine the interannual and decadal modulations in seasonal amplitude of the KS throughflow associated with El Niño-Southern Oscillation (ENSO), Indian Ocean dipole (IOD), Pacific Decadal Oscillation (PDO). Linear regression, correlation, harmonic and power spectrum analyses are used. The results manifest that there are significant interannual to decadal modulations in the seasonal amplitude of the KS throughflow. The modulations of the seasonal amplitude in the volume and heat transports range 1.36-1.92 Sv (1 Sv = 106 m3 s-1) and 126.41-173.36 TW (1 TW = 1012 W), respectively, with a significant cycle of ~9 years. From 1994 to 2020, the seasonal amplitude of volume transport through the KS shows an increasing trend of 37.75 ± 15.69 mSv decade-1 (1 mSv = 103 m3 s-1). The seasonal amplitude of the heat transport also increases, at a rate of 4.78 ± 1.52 TW decade-1. The KS volume transport is positively correlated with PDO and ENSO indices (r2 = 0.69 and r2 = 0.58), with a lag of 12 and 10 months, respectively. The results of composite analysis suggest that the interannual variability of the KS transport is related to the interannual anomalies of the SSH gradient and the local SSW fields in boreal winter.

The Japan/East Sea (JES) is known as a mid-latitude “Miniature Ocean” that features multiscale oceanic dynamical processes. Using principal component analysis (PCA), we investigate the variability of the sea surface chlorophyll-a concentration (SSC) and its bloom timing in the JES based on satellite remote sensing products spanning 1998–2019. The JES SSC exhibits strong seasonal variability and blooms twice annually. The spring bloom is induced under combined factors of increased photosynthetically active radiation (PAR), weakened wind speeds and sea ice melting, and terminated by the enhanced stratification. The fall bloom is induced by destratification and active dynamic processes (such as upwelling and front), and terminated by decreased PAR. The interannual variability of spring and fall bloom occur along the northwestern coast of the JES and in the deep Japan Basin, respectively. The positive SSC anomalies along the northwestern coast of the JES in spring is associated with more sea ice in the previous winter, weaker wind speed, and stronger stratification induced by the El Niño events. No significant relationship has been found between the fall bloom and the El Niño events. The bloom timing is controlled by the critical depth hypothesis. The initiation/termination timing of spring blooms has shifted earlier by 0.37/0.45 days, and the counterpart of fall blooms has shifted 0.49/1.28 days earlier per year. The duration and magnitude are independent with each other for spring bloom at interannual time scale. In contrast, they are positively correlated for fall bloom, because of both bloom timing and magnitude are dominated by active oceanic dynamical processes in fall.

Informer model and obsevation data

Informer model and obsevation data

Using in situ observations collected by a drifting air–sea interface buoy (DrIB) in the northern South China Sea from August 30 to September 13, 2018, the extreme air–sea turbulent fluxes that occurred from September 8 to 13 during tropical cyclone (TC) Barijat were investigated. The most striking features were substantial increases in momentum and heat fluxes, with maximum increases of 10.8 m s−1 in the wind speed (WS), 0.73 N m−2 in the wind stress, 68.1 W m−2 in the sensible heat fluxes (SH) and 258.8 W m−2 in the latent heat fluxes (LH). The maximum WS, wind stress, SH and LH values amounted to 15.3 m s−1, 0.8 N m−2, 70.9 W m−2 and 329.9 W m−2, respectively. Using these new DrIB observations, the performance of two state-of-the-art, high-resolution reanalysis products, ERA5 and MERRA2, was assessed. The consistency of the observed values with ERA5 was slightly better than with MERRA2, reflected in higher correlations but both products underestimated the WS during TC conditions. In calm weather conditions, the turbulent heat fluxes were overestimated, because they simulated a too dry and cold atmospheric state, enhancing the air–sea differences in temperature and humidity. Considering that an accurate representation of the air–sea turbulent and momentum fluxes is essential for understanding and predicting ocean and atmospheric variability, our findings indicate that more high-quality temperature and relative humidity observations are required to evaluate and improve existing reanalysis products.

The South Atlantic subtropical underwater (STUW) is a high-salinity water mass formed by subduction within the subtropical gyre. It is a major component of the subtropical cell and affects stratification in the downstream direction due to its high salinity characteristics. Understanding the interannual variability in STUW subduction is essential for quantifying the impact of subtropical variability on the tropical Atlantic. Using the output from the ocean state estimate of the Consortium for Estimating the Circulation and Climate of the Ocean (ECCO), this study investigates the interannual variability in STUW subduction from 1992 to 2016. We find that heat fluxes, wind stress, and wind stress curl cause interannual variability in the subduction rate. Heat fluxes over the subduction area modulate the sea surface buoyancy and regulate the mixed layer depth (MLD) during its deepening and shoaling phases. Additionally, the wind stress curl and zonal wind stress can modulate the size of the subduction area by regulating the probability of particles entrained into the mixed layer within 1 year of tracing. This analysis evaluates the influence of subtropical wind patterns on the South Atlantic subsurface high-salinity water mass, highlighting the impact of heat and wind on the interannual changes in the oceanic component of the hydrological cycle.

The Japan/East Sea (JES) is known as a mid-latitude “Miniature Ocean” that features multiscale oceanic dynamics processes. We investigate the variability of the sea surface chlorophyll-a concentration (SSC) and bloom timing in the JES based on satellite remote sensing products spanning 1998–2019. The JES SSC exhibits strong seasonal variability and blooms twice annually, which are mainly governed by the physical environmental conditions. However, the influences of local oceanic dynamic processes (e.g., upwelling, oceanic fronts, mesoscale eddies, and near-inertial oscillations) on the bloom magnitude and timing of the entire JES are not critical, compared with the PAT and stratification. In addition, significant interannual variabilities of spring bloom magnitude occur along the JES's northwestern coast, and that of fall bloom magnitude occur in the deep Japan Basin. For spring bloom, the interannual variability of the bloom timing (initiation timing, termination timing and duration), which significantly affect the interannual bloom magnitude anomalies, are correlated with climate modes such as AO and ENSO. For fall bloom, on the interannual time scale, the bloom duration is mainly affected by the initiation timing. Both of them have a significant influence on the bloom magnitude. The initiation/termination timing of spring blooms has shifted earlier by 0.37/0.45 days annually along the JES's northwestern coast; the counterpart of fall blooms has shifted 0.49/1.28 days earlier annually in the deep Japan Basin.

Biogenic secondary organic aerosols (BSOA) are important components of the remote marine atmosphere. However, the response of BSOA changes to sea ice reduction over the Arctic Ocean remains unclear. Here we investigated isoprene and monoterpenes secondary organic aerosol (SOAI and SOAM) tracers in three years of summer aerosol samples collected from the Arctic Ocean atmosphere. The results indicated that methyltetrols were the most abundant SOAI tracers, while the main oxidation products of monoterpenes varied over the years owing to different aerosol aging. The results of the principal component analysis (PCA)-generalized additive model (GAM) combined with correlation analysis suggested that SOAI tracers were mainly generated by the oxidation of isoprene from marine emissions, while SOAM tracers were probably more influenced by terrestrial transport. Estimation of secondary organic carbon (SOC) indicated that monoterpenes oxidation contributed more than isoprene and that sea ice changes had a relatively small effect on biogenic SOC concentration levels. Our study quantified the contribution of influencing factors to the atmospheric concentration of BSOA tracers in the Arctic Ocean, and showed that there were differences in the sources of precursors for different BSOA. Hence, our findings have contributed to a better understanding of the characteristics, sources and formation of SOA in the atmosphere of the Arctic Ocean.

Using observational temperature profiles, this study reveals that eddy-induced oceanic vertical temperature anomalies are geographically varying in the South China Sea (SCS). The maximum anomalies are more intensified in the central and south SCS. Due to ocean stratification, the eddy available potential energy (EAPE) can more comprehensively reflect the eddy effect on the ocean vertical structure. The EAPE within a composite eddy is the highest in the northern region and lowest in the southeastern region. The geographical variation of vertical anomalies induced by eddies is because of the combination effect of the ocean stratification, the background energy (available potential energy and barotropic kinetic energy), and the energy conversion rate. The anomalies induced by cyclonic and anticyclonic eddies might be asymmetrical in the same region. For example, in the regions to the west of Luzon Island and in the central west SCS, the temperature anomalies and EAPE within a composite cyclonic eddy are significantly larger than those within an anticyclonic eddy. This variation is owing to the different oceanic stratification along the tracks of cyclonic and anticyclonic eddies.

Marine heatwaves (MHWs) and low sea surface salinity (SSS) events can significantly impact marine ecosystems and dynamic systems, respectively. Compound marine extreme events can cause more significant damage than individual extreme events. However, the spatiotemporal patterns of compound MHW-low SSS extremes are not well understood. Daily reanalysis data were used to identify the basic patterns of compound extreme events and their drivers. These events mainly occur over the central tropical Pacific Ocean during record-breaking El Niño events. This analysis revealed that extreme sea surface warming associated with El Niño drives increased convection, which subsequently leads to increased rainfall. It ultimately causes extreme sea surface freshening. This analysis highlights the significance of air-sea interactions and low-frequency climate variability in shaping compound extreme events.

Ship type identification is an important basis for ship management and monitoring. The paper proposed a new method of ship type identification by combining characteristic parameters from the energy difference between high and low frequencies and the sensitive IMF variance mean value based on the modal decomposition of the underwater radiated noise signals using the Ensemble Empirical Mode Decomposition (EEMD) method. The comparison shows that the characteristic parameters of different types of ship, underwater radiated noises are different, whereas those of the same types of ship, underwater radiated noises fall in close range. Validation experiments based on randomly selected ship underwater radiated noise samples manifest that the method is of good separability for the four types of ship underwater radiated noises in the Deepship dataset. It has a higher identification rate than other methods within the distance range of ship underwater radiated noise detection in the dataset. The accuracy of this method tends to decrease with distance in the classification experiments of the ship underwater radiated noises at different distances.

Exploring multi-timescale tidal variability is fundamental and necessary for numerous practical purposes, such as flood protection, marine cultivation, and ocean transport. It is well known that tides show significant seasonal, inter-annual, and 18.61-year nodal variability. Less known and less discussed is the subseasonal tidal variability (i.e., ter-annual, quarter-annual, and penta-annual cycles) in the coastal ocean. In this study, we explore subseasonal tidal modulations in the Gulf of Tonkin via the combination of four tide gauges and 27-year multi-satellite altimeter observations. Both tide gauges and satellite altimeters indicate that tidal subseasonality is significant in the Gulf of Tokin, although the amplitudes of subseasonal variations are much smaller than those of seasonal variations. Compared to spatially limited tide gauges, satellite altimeters successfully derive the basin-scale tidal subseasonality in the Gulf of Tonkin. The largest amplitude of subseasonal tidal constituents originated from the subseasonality of main tidal constituents, and can reach as high as 31.8 mm. It is suggested that subseasonal variations in ocean environments (e.g., sea levels and ocean stratification) induce tidal subseasonality through changing tidal propagation and dissipation. Although powerful, satellite altimeters also have some defects. Due to tidal aliasing related to long-period sampling intervals, some subseasonal tidal constituents are indistinguishable in satellite altimeter records.

The North Brazil Undercurrent (NBUC) is a narrow (16 Sv) and plays an important role in the Atlantic Meridional Overturning Circulation and the South Atlantic Subtropical Cell. Strong salinity and temperature fronts occur over the NBUC region. The role of temperature and salinity gradients on the genesis of NBUC variability has never been explored. This study uses three high-resolution (≤0.1°) and one low-resolution (=0.25°) model outputs to explore the linear trend of NBUC transport and its variability on annual and interannual time scales. We find that the linear trend and interannual variability of the geostrophic NBUC transport show large discrepancies among the datasets. Thus, the linear trend and variability of the geostrophic NBUC are associated with model configuration. We also find that the relative contributions of salinity and temperature gradients to the geostrophic shear of the NBUC are not model dependent. Salinity-based and temperature-based geostrophic NBUC transports tend to be opposite-signed on all time scales. Despite the limited salinity and temperature profiles, the model results are consistent with the in-situ observations on the annual cycle and interannual time scales. This study shows the relationship of salinity-based and temperature-based geostrophic NBUC variations in the annual and interannual variability and trend among different models and highlights the equal important roles of temperature and salinity in driving the variability of NBUC transport.

This study focuses on different evolutions of the low-level atmospheric circulations between eastern Pacific (EP) El Niño and central Pacific-II (CP-II) El Niño. The western North Pacific anomalous anticyclone (WNPAC) originates from the northern South China Sea for EP El Niño, and moves to the western North Pacific (WNP) afterward. Compared with EP El Niño, the origin of the WNPAC is farther west during CP-II El Niño, with the center over the Indochina Peninsula. Moreover, the WNPAC shows a weaker eastward shift. Such discrepancies are attributed to different evolutions of the cyclonic response over the WNP, which can suppress the convection in the western flank of the anomalous cyclone. The eastward retreat of the anomalous cyclone is significant for EP El Niño, but less evident for CP-II El Niño. These discrepancies are related to zonal evolutions of the increased precipitation over the equatorial Pacific. Following the southward migration of the intertropical convergence zone (ITCZ), the deep-convection region extends eastward along the equator, reinforcing the atmospheric response to the eastern Pacific warming in EP El Niño. For CP-II El Niño, the atmospheric response is insignificant over the eastern Pacific without warming. Moreover, the meridional migration of the ITCZ can modulate zonal variations of the easterly trade wind and specific humidity as well. Due to the combined effects of the climatological background and atmospheric anomalies, the specific humidity–induced and wind-induced moist enthalpy advection contribute to different shifts of the precipitation center.

This study examined the characteristics and vertical propagation of near inertial waves (NIWs) induced by Typhoon Linfa (2015), based on in situ observations conducted southeast of Dongsha Islands in the South China Sea. The results demonstrate that the near inertial currents induced by Linfa had velocities up to 35cm s−1 in the mixed layer and 20cm s−1 in the ocean interior. The near inertial currents were polarized with predominantly clockwise-rotating components, the magnitudes of which were about 10 times larger than the counter-clockwise rotating components. The energy density spectrum showed that the emergence of NIWs resulted in energy redistribution from the diurnal band to the near inertial band. The wavenumber spectrum and the downward/ upward current decomposition demonstrated that the NIWs and energy flux propagated mainly downward. The estimated vertical phase velocity and group velocity are 1.44 and 0.48m h−1, respectively, corresponding to a vertical wavelength of 49.7m. The e-folding time scale was 7.5 d based on the near inertial kinetic energy in the ocean interior. We found no obvious wave-wave interaction during the decay process of the NIWs. The frequency was blue-shifted, being 0.03f0 higher than the local inertial frequency, which was caused by the background vorticity. The normal mode analysis suggests that the higher mode plays a dominant role in the propagation stage of the NIWs.

Numerous oceanic mesoscale eddies occur in the South China Sea (SCS). The present study employs eight automatic eddy detection algorithms to identify these mesoscale eddies and compares the results. Eddy probabilities and areas detected by various algorithms differ substantially. Most regions of the SCS with a high discrepancy of eddy probabilities are those with few mesoscale eddies, except for the area west of the Luzon Strait, the area west of Luzon Island between 12° and 17°N, and the southernmost end of the SCS basin. They are primarily caused by strong interference, noncircular eddy shapes, and gentle sea level anomaly (SLA) gradients, respectively. The SLA, winding angle, and hybrid methods can easily detect the mesoscale eddies with wavelike features. The Okubo–Weiss (OW) and the spatially smoothed OW methods better identify grouping phenomena of mesoscale eddies in the SCS. Suggestions are presented on choosing suitable algorithms for studying mesoscale eddies in the SCS. No single algorithm is perfect for all research purposes. For different studies, the most suitable algorithm is different.

The South China Sea (SCS) is a semi-enclosed ocean basin, with its sole deepwater renewal channel in the northeast, the Luzon Strait, connected to the northwestern Pacific. The circulation in the SCS is endowed with a distinctive multi-layer structure, which is influenced by the seasonally reversing Southeast Asian monsoon at the surface, modulated Kuroshio intrusion in the upper layer, and enhanced diapycnal mixing and persistent Luzon Strait deepwater overflow in the deep layer. After being diapycnally mixed by energetic internal tides, internal waves, and mesoscale eddies, waters entering the SCS in the upper and deep layers through the Luzon Strait are returned to the northwestern Pacific in the intermediate layer. This inflow-outflow-inflow structure in vertical induces net inflows of potential vorticity in the upper and deep layers and a net outflow of potential vorticity in the intermediate layer of the Luzon Strait, which are believed to be closely related with the circulation structure in the SCS. A sandwiched circulation in the SCS, namely, anticyclonic in the intermediate layer and cyclonic in the upper and deeper layers, was recently reported, but controversial issues still remain on its spatial pattern and temporal variability. To ascertain this sandwiched circulation and reveal its driving mechanism is of much scientific interest to the on-going discussion on the SCS circulation and its role in regional climate. Different from previous review articles focusing on eddies or circulation in a specific layer or area, we overview the full-depth circulation of the SCS layer-by-layer and outline its possible forcing mechanisms, and furthermore, provide perspectives for further investigation.

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.

Near-bottom currents play important roles in the formation and dynamics of deep-water sedimentary systems. This study examined the characteristics and temporal variations of near-bottom currents, especially the tidal components, based on two campaigns (2014 and 2016) of in situ observations conducted southeast of the Dongsha Island in the South China Sea. Results demonstrated near-bottom currents are dominated by tidal currents, the variance of which could account for ~70% of the total current variance. Diurnal tidal currents were found stronger than semidiurnal currents for both barotropic and baroclinic components. The diurnal tidal currents were found polarized with predominantly clockwise-rotating constituents, whereas the clockwise and counterclockwise constituents were found comparable for semidiurnal tidal currents. It was established that diurnal tidal currents could induce strong current shear. Baroclinic tidal currents showed pronounced seasonal variation with large magnitude in winter and summer and weak magnitude in spring and autumn in 2014. The coherent components accounted for ~65% and ~50% of the diurnal and semidiurnal tidal current variances, respectively. The proportions of the coherent and incoherent components changed little in different seasons. In addition to tidal currents, it was determined that the passing of mesoscale eddies could induce strong nearbottom currents that have considerable influence on the deep circulation.

Four trawl-resistant bottom mounts, with acoustic Doppler current profilers (ADCPs) embedded, were deployed in the Karimata Strait from November 2008 to June 2015 as part of the South China Sea-Indonesian Seas Transport/Exchange and Impact on Seasonal Fish Migration (SITE) Program, to estimate the volume and property transport between the South China Sea and Indonesian seas via the strait. The observed current data reveal that the volume transport through the Karimata Strait exhibits significant seasonal variation. The winteraveraged (from December to February) transport is–1.99 Sv (1 Sv=1×106 m3/s), while in the boreal summer (from June to August), the average transport is 0.69 Sv. Moreover, the average transport from January 2009 to December 2014 is–0.74 Sv (the positive/negative value indicates northward/southward transport). May and September are the transition period. In May, the currents in the Karimata Strait turn northward, consistent with the local monsoon. In September, the southeasterly trade wind is still present over the strait, driving surface water northward, whereas the bottom flow reverses direction, possibly because of the pressure gradient across the strait from north to south.

The Island Rule, derived from the Sverdrup theory, is widely used to estimate and analyze water transport through a strait. Previous studies presented single- or multi-island rules with either lateral or bottom friction. In this paper, an analytical model of wind-driven circulation is assumed based on linear dynamics. Considering both lateral and bottom friction, the analytic solutions of the transport streamfunction around the islands are derived and the volume transport through the channel is presented. The results are similar to those of Wajsowicz, but the frictional constants represent different values. The analytic solution shows that the relationship between the lateral frictional and bottom frictional dissipation is complex in terms of the frictional constants. To understand the interaction between the two friction types, lateral and bottom friction values were randomly chosen on a barotropic beta plane. The result shows an approximately linear relationship between the lateral and bottom friction in consisting of the combined frictional constants. We studied the effect of the channel width on the transport through the channel. The results show that the friction enhances the flow under some widths, which is similar to the flow behavior when only the lateral friction is considered. We also compared the transport through the channel at different depths and founded that the deeper the water, the smaller the transport reduction ratio when the horizontal eddy viscosity coefficient and the bottom drag coefficient remained constants. To further present the combined role of lateral frictional and bottom frictional dissipation, we compared our model with the model of Wajsowicz for two islands, where only the lateral or bottom friction were considered, with different channel widths. The results showed that the effect of the lateral friction is greater than the bottom friction when the channel is narrow, especially in the Munk boundary layer thickness. When the channel is much wider than the Munk boundary layer thickness, the role of the bottom friction is greater than that of the lateral friction. The model was applied to the Indonesian throughflow and yielded a reduction of approximately 20% in the transport.

Through analysis of the results of a verified high-fidelity numerical model, the intra-seasonal variations (ISVs) in the depth of the 22°C isotherm (D22) in the South China Sea (SCS) basin are investigated. The results show that the ISVs in the D22 exhibit distinct seasonality in the SCS. The ISVs in the D22 are quite significant, especially within a band along the northwestern boundary of the basin and at the southern end of the basin during boreal winter. In these areas, the ratio of the standard deviations (STDs) of intra-seasonal band to the STDs of total data could exceed 0.6. Although the ISVs in the D22 are detectable in the area affected by the Vietnam Offshore Current during boreal summer and autumn, these variations are sometimes overwhelmed by oscillations with other frequencies. An analysis of the causes of the ISVs in the D22 in the SCS indicates that sea surface fluxes and wind stirring are not the dominant external driving mechanisms of the phenomena described above. The ISVs in the D22 are thought to be induced mainly by the thermodynamic adjustment of the ocean itself and the associated instabilities. The energy of the northern and southern bands that display strong ISVs in the D22 may be derived from eddy kinetic energy, rather than eddy available potential energy. The diversity of the propagation of the ISVs in the D22 is very conspicuous within these two bands.

A 72-h fine-resolution atmosphere-wave-ocean coupled forecasting system was developed for the South China Sea and its adjacent seas. The forecasting model domain covers from from 15°S to 45°N in latitude and 99°E to 135°E in longitude including the Bohai Sea, the Yellow Sea, the East China Sea, the South China Sea and the Indonesian seas. To get precise initial conditions for the coupled forecasting model, the forecasting system conducts a 24-h hindcast simulation with data assimilation before forecasting. The Ensemble Adjustment Kalman Filter (EAKF) data assimilation method was adopted for the wave model MASNUM with assimilating Jason-2 significant wave height (SWH) data. The EAKF data assimilation method was also introduced to the ROMS model with assimilating sea surface temperature (SST), mean absolute dynamic topography (MADT) and Argo profiles data. To improve simulation of the structure of temperature and salinity, the vertical mixing scheme of the ocean model was improved by considering the surface wave induced vertical mixing and internal wave induced vertical mixing. The wave and current models were integrated from January 2014 to October 2015 driven by the ECMWF reanalysis 6 hourly mean dataset with data assimilation. Then the coupled atmosphere-wave-ocean forecasting system was carried out 14 months operational running since November 2015. The forecasting outputs include atmospheric forecast products, wave forecast products and ocean forecast products. A series of observation data are used to evaluate the coupled forecasting results, including the wind, SHW, ocean temperature and velocity. The forecasting results are in good agreement with observation data. The prediction practice for more than one year indicates that the coupled forecasting system performs stably and predict relatively accurate, which can support the shipping safety, the fisheries and the oil exploitation.

Understanding nodal tidal characteristics is essential for accurate long-term tidal prediction. Observational nodal evolution of tides is mainly based on tide gauge records in coastal areas which are limited in time and space, thus impeding coherent determinations of basin-wide patterns of tidal variability. In this paper, we indicate the potential of satellite altimeter data to investigate 18.61-year nodal modulations of main constituents in the Gulf of Tonkin. Three tide gauges and multi-source satellite altimeter observations (TOPEX/Poseidon, Jason1, Jason2, and Jason3) revealed that 18.61-year nodal cycles in tidal amplitudes have noticeable deviations from the equilibrium tidal theory in the Gulf of Tonkin. In general, M2 and N2 nodal modulations are anomalously larger than theoretical values while K2, K1, and O1 nodal modulations are noticeably smaller than theoretical values. Compared to point-based tide gauges, satellite altimeter records can provide basin-wide features of nodal modulations of main constituents. Although overlapping geographical blocks are applied to eliminate the effect of tidal alias originated from long-period sampling intervals, the estimation of nodal cycles of minor constituents are still questionable. Nevertheless, the methods described here provide a strong foundation for future research on time-varying tidal dynamics using the combination of tide gauges and satellite altimeter data.

As a ubiquitous movement in the ocean, tides are vital for marine life and numerous marine activities such as fishing and ocean engineering. Tidal dynamics are complicated in the East Asian marginal seas (EAMS) due to changing complex topography and coastlines related to human activities (e.g., land reclamation and channel deepening) and natural variability (e.g., seasonal variations of ocean stratification and river flow). As an important tool, numerical models are widely used because they can provide basin-scale patterns of tidal dynamics compared to point-based tide gauges. This paper aims to overview the development history of the numerical simulation of tides in the EAMS, including the Bohai Sea, the Yellow Sea, the East China Sea, the East/Japan Sea, and the South China Sea, provide comprehensive understanding of tidal dynamics, and address contemporary research challenges. The basic features of major tidal constituents obtained by tidal models are reviewed, and the progress in the inversion of spatially and temporally changing model parameters via the adjoint method are presented. We review numerical research on how a changing ocean environment induces tidal evolution and how tides and tidal mixing influence ocean environment in turn. The generation, propagation, and dissipation of internal tides in the EAMS are also reviewed. Although remarkable progresses in tidal dynamics have been made, nonstationary tidal variations are not fully explained yet, and further efforts are needed. In addition, tidal influences on ocean environment still receive limited attention, which deserves special attention.

The green tides caused by Ulva prolifera have become a recurrent phenomenon in Yellow Sea, China. Investigating the factors governing the biomass of green tides is important for developing management strategies. In this study, an U. prolifera growth model was combined with a hydrodynamic model. This biophysical model can reasonably reproduce the spatiotemporal variation of the green tides in 2012. Among three zones (northern, central, and southern-zones) of Porphyra mariculture region, the northern and central zones were more important in controlling the bloom intensity, and the central zone was the key area in controlling the amount of biomass landed on beaches. Due to the limitation of temperature and nutrients, an earlier or postponed facility recycling might effectively reduce the magnitude of green tides in 2012. This study provides useful information for mitigation of green tides and management of Porphyra mariculture.