1. Introduction
The Luzon Strait, which is 300 km wide and located between Taiwan Island and Luzon Island, is a key channel connecting the South China Sea (SCS) and the Pacific Ocean. The warm and saline Kuroshio waters from the Pacific enter the SCS through this gap and are of vital importance to the thermohaline structures and circulation of the SCS. The intrusion of Kuroshio waters into the SCS appears in different patterns and can be classified into direct intrusion branches, loop currents, and detached anticyclonic eddies (Hu et al. 2000; Caruso et al. 2006; Nan et al. 2011a, 2015). The loop currents (or the Kuroshio loops) are not persistent, but occur occasionally (Sun et al. 2020). The presence of a cyclonic eddy often facilitates an anticyclonic eddy shed from these loop currents (Zhang et al. 2017). In rare cases, nonlinear Rossby waves from the Pacific can directly penetrate the SCS (Hu et al. 2012). Winter is the most favorable season for Kuroshio intrusion; in contrast, summer Kuroshio intrusion is rarely reported in the literature (e.g., Yuan et al. 2014; Yang et al. 2019). The seasonality of the Kuroshio intrusion has been attributed to the seasonal variability of winds (Farris and Wimbush 1996; Wu and Hsin 2012) or the strength of the upstream Kuroshio (Sheremet 2001; Sheu et al. 2010; Yuan and Wang 2011).
The west of the Luzon Strait is one of the regions with the most frequent mesoscale eddies in the SCS (Wang et al. 2003; D. Wang et al. 2008; G. Wang et al. 2008; Xiu et al. 2010; Chen et al. 2011; Zhang et al. 2013, 2016, 2017). These eddies can maintain the properties of water in their upper layer and transport them to other regions (Zu et al. 2013; Yang et al. 2019) and have greater impacts than seasonal winds on the thermohaline structure of the water column (Zhao et al. 2017). Strong seasonality is a feature of eddies generated here. For example, in winter, the Luzon Cold Eddy (LCE) occurs northwest of Luzon Island (Yang and Liu 2003), whereas in summer the Luzon Warm Eddy (LWE) is dominant. Their generation is closely related to the wind stress curl (WSC) from orographic wind jets (G. Wang et al. 2008). However, LWEs identified from altimeter maps have only been captured by single hydrographic surveys and a few Argo profiles, with relatively little research attention (Li et al. 1998; Yuan et al. 2007; Chen et al. 2009).
Li et al. (1998) first reported an anticyclonic eddy carrying Kuroshio waters near the Dongsha Islands in the summer of 1994. However, by analyzing merged sea surface heights (SSHs) derived from altimeters, Yuan et al. (2007) argued that this eddy originated northwest of Luzon Island instead of the Kuroshio. They also found that the occurrence of the anticyclonic eddy is seasonal, forced by the southwest monsoon and the negative WSC generated northwest of Luzon Island in summer, and moves toward the continental shelf in fall and winter, leading to current variations. The same authors note, however, that the generation of eddies could not be fully explained by seasonal winds. Chen et al. (2009) demonstrated the seasonal pattern of the LWE from the empirical orthogonal function (EOF) decomposition of merged sea level anomaly data. From the temperature and salinity (T/S) properties of Argo profiles inside an LWE during 2006, they suggested that the LWE incorporated mixed waters from the Kuroshio and the SCS. From a summer hydrographic cruise, Yang et al. (2019) reported an anticyclonic eddy trapping and transporting saline Kuroshio waters into the SCS. Although this eddy was not an LWE, it originated from the central Luzon Strait. This implies that the anticyclonic eddies generated around Luzon Island in summer may be closely related to the Kuroshio intrusion.
To date, studies of LWEs remain limited. Because the LWEs reported by Yuan et al. (2007) and Chen et al. (2009) originated west or northwest of Luzon Island, the waters that they trapped came mostly from the SCS; however, in some cases, LWEs can arrest large amounts of Kuroshio waters when generated in the region north of Luzon Island near the Kuroshio. These eddies are named North Luzon Warm Eddies (NLWEs) to distinguish them from previously reported LWEs. Here, we report two NLWEs that occurred during the summer of 2018. Specifically, we used observational data from an array consisting of 27 current- and pressure-recording inverted echo sounders (CPIESs) to investigate the structures of these NLWEs.
The rest of the paper is organized as follows: section 2 describes the data used and the methods for CPIES data processing and estimation; section 3 describes the observations, including the evolution of the NLWEs, the evidence for Kuroshio waters, and the structure of the NLWEs; section 4 discusses the generation mechanism of the NLWEs; and section 5 provides a summary and conclusions.
2. Data and methods
a. Data
CPIESs are a type of bottom mooring equipped with a mounted instrument measuring acoustic travel time from the sea surface and back (τ), bottom pressure (Pbot), and near-bottom velocities 50 m off the seafloor (v). The body of the CPIES is seated 1 m off the seafloor on a mooring platform with additional iron anchors. One CPIES can estimate the steric and mass-loading components of SSH from τ and Pbot, respectively. Large CPIES arrays are deployed in regions such as the Kuroshio Extension (Donohue et al. 2010), the Japan/East Sea (Mitchell et al. 2005), and the Antarctic Circumpolar Current (Tracey et al. 2013), and they have proved successful in mapping the three-dimensional density and geostrophic current field.
A large CPIES array was deployed in the SCS from June 2018 to August 2019 under the joint efforts of the Second Institute of Oceanography, Ministry of Natural Resources, and Institute of Oceanology, Chinese Academy of Sciences. The array spanning a 230 km × 400 km region west of the Luzon Strait comprised 28 CPIESs (Fig. 1). To the best of our knowledge, the scale of the CPIES mooring array is the largest ever reported for the SCS. One CPIES (C19) failed to recover data for unknown reasons, so data were retrieved from the remaining 27 arrays. CPIESs C21, C23, C29, and C30 were deployed in December 2017 as planned and other deployments were postponed to the following summer because of bad weather conditions. Two tall moorings (M1 and M2; Table 1) consisting of RDI Workhorse acoustic Doppler current profilers (ADCPs; 75 kHz), conductivity–temperature–depth measurers (CTDs; Seabird SBE37SM), current meters (CMs; Nortek, Aquadopp), and temperature chains (Seabird SBE56) were deployed near C21 and C27. Each ADCP observed 74 bins with a size of 8 m. All the mooring data were detided and applied with a fourth-order Butterworth low-pass filter going forward and backward to avoid phase shift with a cutoff period of 72 h, and were subsampled at one-day intervals. A summary of the mooring information is provided in Table 1.
Information of the tall moorings.
We also used observations from the following sources: 6-hourly sea surface wind products from the Cross-Calibrated Multi-Platform (CCMP) with a spatial resolution of 1/4° × 1/4°; daily merged sea surface height anomaly data products from the European Copernicus Marine Environment Monitoring Service (CMEMS) with a spatial resolution of 1/4° × 1/4°; Argo profiles (Liu et al. 2021) provided by the China Argo Real-time Data Center (http://www.argo.org.cn/); and General Digital Environmental Model (GDEM; version 3.0) monthly climatology data with 78 standard depths and a horizontal resolution of 1/4° × 1/4°. The GDEM incorporated temperature and salinity profiles from the Mater Oceanographic Observational Dataset of the Naval Research Laboratory (NRL). The GDEM has been found to perform well in reproducing hydrographic mean fields in the ocean (Qu et al. 2006; Carnes 2009; Wang et al. 2011).
b. Methods
1) Processing of the raw CPIES data
All CPIESs were configured to measure τ, Pbot, and v at intervals of 24 pings per hour, 6 per hour, and hourly, respectively. A sliding event was identified at the C13 and C29 stations, where the Pbot records suffered an abrupt “jump.” This effect was removed by adding an offset (dp) to the Pbot records after the event. The offset was determined by iteration, whereby the detided residuals of the entire time series had a minimum deviation. The τ time series after this jump was also adjusted accordingly by adding the following offset:
2) Estimation of the three-dimensional field
The optimal interpolation (OI) technique adapted from Bretherton et al. (1976) was used to map the two-dimensional τ1000 with a correlation length scale of 120 km. The far-field mean τ1000 was converted from the GDEM climatology. The three-dimensional specific volume anomaly field was determined from the GEM lookup table, and the geopotential height was derived by its integration on the vertical. Then, we calculated the baroclinic velocity shears Vbc(p) from the geopotential height estimates by assuming geostrophy. For each CPIES, we obtained velocities at the reference depth (1000 m) of each station by adding baroclinic velocity shear Vbc(Pref) − Vbc(Pbot) between the reference depth and the ocean bottom to the measured near-bottom velocities V(Pbot). The velocity field of the reference depth was further obtained by combining the above velocities and Pbot anomalies using a multivariate OI scheme (Watts et al. 2001; Donohue et al. 2010) with a correlation length of 100 km.
3) Estimation of the sea surface height
3. Results
a. NLWEs identified from altimeters
Two NLWEs influenced our observation region successively during the summer–fall period of 2018; their centers are marked as black stars and triangles in Fig. 4. The first NLWE (herein “eddy A”) was generated northeast of Luzon Island in early July and was accompanied by a strong negative WSC forced from the atmosphere (Fig. 5). We also observed a weak anticyclonic eddy northwest of Luzon Island on 3 July, which was absorbed by eddy A on 7 July to form a larger eddy. The merged new eddy continued to move northwestward and weakened until a weak anticyclonic SSH anomaly was generated northeast of Luzon Island and merged again on 27 July. Eddy A grew extremely large after the merging event and peaked on 8 August. The eddy covered the region southwest of Taiwan Island, spanning 118.5°–119°E and 19.5°–22.2°N. It then moved westward and quickly dissipated after 12 August.
The second NLWE (“eddy B”) formed in early August and originated northeast of Luzon Island when there was a strong negative WSC in the region (Fig. 5). Eddy B remained north of Luzon Island for almost three weeks, and then moved northwestward following a similar path to eddy A from September. Eddy B was much smaller than eddy A, with one likely contributing factor being that it did not absorb other anticyclonic eddies during its evolution. Once reaching southwest Taiwan, eddy B did not turn westward as in the case of eddy A, but continued its heading and moved out of our observational region after 20 September.
The difference between our observed anticyclonic eddies and that reported by Yuan et al. (2007) is that their LWE originated northwest rather than northeast of Luzon Island. Nevertheless, there are many common characteristics of these anticyclonic eddies that are generated in summer to early fall north of Luzon Island. The eddies we observed seemed to be closely related to forcing by the negative WSC, which is most favorable in summer. Because their origins are located closer to the path of the Kuroshio, it is expected that these NLWEs may carry Kuroshio waters to the SCS during their evolution, as discussed in the following section.
b. Evidence of Kuroshio waters intrusion
1) CPIES-estimated SSH
As demonstrated in section 2b, CPIESs can be used to estimate SSH. Because the barotropic components of SSH are an order of magnitude smaller than the baroclinic components, the estimation of SSH essentially involves the linear mapping of τ based on its relationship with different waters. As shown in Fig. 3c, the estimated SSHs from Kuroshio waters (HKur) are higher than from SCS waters (HSCS) with the same τ (τ > 1.321 s). Therefore, we analyzed the SSH estimations with different hydrocast origins and compared them with the merged SSH products derived from satellite altimetry (Hsat).
The Hsat of each CPIES station was interpolated from the merged absolute dynamic topography (MADT) of the CMEMS. It is evident that the estimated HSCS values are generally in good agreement with Hsat, especially for the stations in the SCS interior (Fig. 6). In comparison, the correlations are much weaker near the Luzon Strait, probably because of the strong mixing of different waters. From July to September when our observational region was impacted successively by the two NLWEs, Hsat northeast of the Luzon Strait increased and appeared as bumps or humps. These humps were significantly underestimated by HSCS but were better reproduced by HKur. This strongly suggests that NLWEs carry large portions of Kuroshio waters, while HSCS failed to capture these SSH signals. The disagreement in SSHs is unlikely attributable to Hsat mapping errors, which tend to be random (Le Traon and Morrow 2001). There is no evidence that the disagreement is attributable to internal tides, which may bring considerable errors to the altimeter SSH (Ray and Zaron 2011; Zhao 2014), because the observed internal tidal signals during July to September were not particularly strong (Fig. S3). Because the SCS hydrographic dataset contains Kuroshio waters incorporated during the winter, and the signals are fully captured by the seasonal correction scheme of SSH estimation, HSCS generally matched Hsat during the winter when the Kuroshio intrusion was the strongest.
2) Argo-observed salinity
One Argo profile (WMO: 2901482) entered our observational region and was captured by eddies A and B (Fig. 4). The Argo float (serial number 6875) was deployed under the U.S. Argo project in May 2018, headed north and turned eastward into July, and then headed north again after 23 July. The float drifted around the rim of eddy A and entered its interior between 4 and 12 August. Although it never reached the center of eddy A, the maximum salinity of the profiles increased to values above 34.7 psu, as is adopted in the literature (Hu et al. 2012) as a criterion for Kuroshio water intrusion. The core of the high-salinity waters (Fig. 7) was located at a depth of 140 m with a thickness of less than 40 m. From the T/S scatters of the Argo profiles, it was evident that this small patch of saline water was closer to the mean Kuroshio waters; however, the waters at other depths had T/S properties closer to those of the SCS. Because the profiles only stayed at the rim of the anticyclonic eddy, saline Kuroshio waters were not captured fully and were expected to be concentrated at the center.
When eddy A quickly moved westward after 16 August, the Argo float exited the eddy and the measured maximum salinity dropped below 34.7 psu. It then stayed near 21°N, 120°E west of the Luzon Strait and encountered eddy B on 5 September. In this case, it entered the center of eddy B and recorded a maximum salinity of 34.95 psu (on 5 September). The saline water parcels (S > 34.7) occupied a depth range of 100–200 m. The T/S scatters further confirmed the existence of Kuroshio waters; the Argo float reported warm and saline waters closer to the Kuroshio origins in the upper layer. For depths with temperatures greater than 14.8°C, where the T/S curves of mean SCS waters and Kuroshio waters intersect (Fig. 3b), the waters maintained their SCS properties. In this case, the Argo profile was located at the center of eddy B. Therefore, we assumed that the profile captured the Kuroshio waters to the greatest extent in the vertical direction. Although the deep waters inside the NLWEs originated from the SCS, the CPIES-estimated SSH from Kuroshio-derived profiles (HKur) still showed better agreement with Hsat because SSH is mostly determined by steric changes in the upper layer.
c. NLWEs from CPIES estimation
1) Validation of CPIES estimates
To examine the performance of the CPIES estimates, data in the upper layer (<1000 m) from the two moorings were used for comparison. We considered the vertical motion of the moorings by interpolating the CPIES estimates to the changing depths of the mooring instruments. Note that the temperature chain data constructed using the SBE56 temperature loggers were not used in this comparison because they lacked pressure records.
As shown in Fig. 8, there was good agreement between the temperature estimates from the CPIES (red) and the observations from the ADCP/CTD. The root-mean-square errors (RMSEs) of the ADCPs of M1 and M3 were 0.33° and 0.16°C, respectively. The RMSEs for the CTDs ranged between 0.14° and 0.19°C, which is small compared to the temperature variations of the observation depths. On 28 October 2018, the CPIES overestimated and underestimated temperature by more than 0.5°C at depths shallower and deeper than 500 m, respectively. This disagreement is likely attributed to mixing events. The performance of the temperature estimates was expected because the GEM RMSEs were less than 0.2°C for depths greater than 400 m. Salinity record comparisons (Fig. S4) are not significant because salinity contributes little to density, and salinity variations are generally small in the SCS below 300 m.
We also used the Argo profiles that were arrested by NLWEs to check the performance of CPIES temperature and salinity estimates in the upper layer. As shown in Figs. 9a and 9c, the CPIES estimated profiles (red lines) match the Argo observations (blue lines) quite well. We also note that the subsurface saline waters (S > 34.7 psu) were well reproduced by the two-GEM method; the estimated subsurface salinity will be much lower if only the GEM of the SCS is used. Throughout the water column, the mean temperature and salinity differences obtained by CPIES estimates and Argo observations are generally less than 1°C and 0.1 psu, respectively. Exceptions to this are depths near 100 and 200 m, where the temperature difference is much lower and higher, and the salinity difference much higher and lower, respectively. These patterns may reflect the mixing process of different waters given that the Argo-observed temperature and salinity profiles have great deviations. Based on these results, we conclude that the CPIES estimates can accurately reproduce the thermohaline structure of NLWEs that trap Kuroshio waters in the subsurface.
Figure 10 shows the zonal and meridional velocities of the upper 250 m from the ADCP observations and CPIES estimates for M1 and M2, respectively. A 20-day running average was applied to all time series. According to ADCP observations, M1 had relatively large velocity variations associated with strong events, such as mesoscale eddies or Kuroshio intrusions; M2 had much lower velocity variations and was only affected by some weaker eddies. The CPIES estimates successfully captured these events and matched the ADCP observations quite well. For the zonal components, the RMSE was less than 0.09 and 0.04 m s−1 for M1 and M2, respectively; for the meridional components, the RMSE was less than 0.09 m s−1 in both cases. The errors in the calculated meridional velocities were much larger because there were no CPIES stations to the east of M1 and M2. Based on these comparisons, we found that the CPIES array accurately estimated the three-dimensional temperature (density) and velocity field. The structures of the NLWEs revealed by the CPIES estimates are discussed in the next section.
2) Three-dimensional structures of NLWEs
The horizontal distributions of temperature, salinity, and velocity fields from the CPIES estimations during the evolution of the NLWEs are shown in Fig. 11. Because eddy-induced variability was mainly concentrated in the upper layers, we only illustrate the estimates at depths of 0, 150, 300, and 500 m. From July to September 2018, two anticyclonic eddies successively entered the CPIES observational region and propagated northwestward. These maps were consistent with those from the MSLA. With the arrival of NLWEs, high-temperature and salinity cores could be identified in all layers except for the surface high-temperature cores. These cores were not significant for eddy A, for which salinity at 500 m was slightly lower (∼0.02 psu). The NLWEs were also poorly distinguishable from satellite-mapped sea surface temperature (SST) products (data not shown), which might be because of the SCS summer stratification feature. Eddy-induced anomalies were most significant in the subsurface and decreased significantly with depth. With the GDEM monthly mean subtracted, the maximum temperature and salinity increase at 150 m was 3.98°C and 0.31 psu, respectively. The Kuroshio waters carried by the NLWEs were wrapped by the 34.7 psu isohaline at 150 m.
To further demonstrate the vertical structures of the NLWEs, we chose a section (the dotted line in Fig. 4) that approximately follows the propagation path. Figure 12 shows the distributions of temperature, salinity, and velocities across this section. It is evident that the NLWEs, characterized by subsurface saline waters (S > 34.7 psu), moved northwestward along this section. When the centers of the NLWEs moved around the section, water cores with salinities higher than 34.8 psu were observed. Anticyclonic velocity structures can also be seen around the center of the eddy. Superimposed on the westward mean flow, the eddy-induced velocity field was not symmetrical but was stronger in the negative components (i.e., southwestward) across the section. This implies that mean flows are important for the evolution of NLWEs.
4. Discussion
Based on section 3a, we identified that the generation of the studied NLWEs was accompanied by negative WSC over the northeast of Luzon Island. Yuan et al. (2007) discussed the relationship between NLWEs and monsoons. In their study, the variation of winds northwest of Luzon Island led to local SSH variation by one to two months, and the seasonal negative WSC was quite weak, such that seasonal winds did not seem to be directly responsible for the generation of NLWEs. As is well known, the winds and SSH in the SCS are dominated by seasonal signals, indicating that a strong correlation is not sufficient to interpret their relationship.
To assess the role of winds on the generation of NLWEs, we first identified historical (1993–2021) anticyclonic eddies generated within the blue box (Fig. 13a) north of Luzon Island from MSLA products. An eddy center was defined as the point with the maximum SSH within the area of the eddy. A total of 123 anticyclonic eddies were identified, of which 71 were locally dissipated and only 52 (Table S1) traveled farther to the northwest or west (Fig. 13a). The mean date of eddy generation was 11 July. Their statistics suggest that these eddies can only be generated between April and October, and summer (June–August) is the most favorable season for eddy generation. Composites of WSC during summer (Fig. 13c) and winter (Fig. 13d) show significant differences; a negative WSC core dominates the region north of Luzon Island in summer, while a positive WSC core dominates in winter. The strong positive WSC is likely to drive seasonal LCE (Yang and Liu 2003) and is responsible for the absence of NLWEs in winter. Although the negative WSC (up to −2.8 × 10−7 N m−3) in summer is much weaker than the positive WSC in winter (up to 1.2 × 10−6 N m−3), it provides the condition under which NLWEs can be generated.
To further investigate the role of WSC in driving the NLWEs, we extracted the signals of winds and SSH by averaging the WSC and MSLA covering the area 18°–19.5°N, 120°–122°E. We applied a fourth-order Butterworth low-pass filter going forward and backward with a cutoff period of 300 days to the spatially averaged WSC and MSLA time series to obtain the seasonal signals. Because MSLA is a smoothed product that lacks high-frequency signals, a bandpass filter with cutoff periods of 20 and 300 days was applied to both the WSC and MSLA.
Seasonally, WSC and SSH have a maximum correlation coefficient of −0.73 when the WSC leads SSH by 20 days (Fig. 14a). The lag time (20 days) was close to the lag (one month) with the maximum cross-correlation between local WSC and SSH northwest of Luzon calculated by Yuan et al. (2007). The mean annual cycle of regional WSC showed that this reached its minimum on 15 July, while regional SSH reached its maximum on 3 August (19 days later). The seasonal signal of the spatially averaged WSC was only negative a few days around 15 July. Such a weak seasonal WSC is unlikely to directly drive NLWEs. The weak seasonal WSC was also noticed by Yuan et al. (2007); although these authors suggested that seasonal WSC played an important role in generating LWEs, they reported that the detailed generation mechanism of LWEs is complicated and may be related to the momentum imbalance of the North Luzon Current.
On intraseasonal time scales, WSC and SSH have a maximum correlation coefficient of −0.49 when WSC leads SSH by 3 days (Fig. 14b). The WSC and SSH time series span 19 years, from 2002 to 2021. The lag time was consistent with the SSH and WSC distributions shown in Figs. 4 and 5, with positive SSH anomalies being generated to the north of the Luzon Island soon after the prevailing negative WSC. During our observational period, the correlation reached −0.68 when WSC led SSH by 3 days. Moreover, the mean generation date for NLWEs (11 July) was closer to the seasonal WSC minimum (15 July) than the seasonal SSH maximum (3 August). This implies that the generation of NLWEs is more closely linked to negative WSC over intraseasonal time scales than seasonally. In other words, in winter, negative WSC is rare to the north of the Luzon Island because of the strong positive seasonal WSC, while in summer, although the negative seasonal WSC is weak, NLWEs can be triggered by intraseasonal negative WSC to the north of Luzon Island.
5. Summary and conclusions
A large CPIES array provided an unprecedented opportunity to investigate the rich multiscale dynamical processes near the Luzon Strait, where phenomena such as strong eddies and Kuroshio intrusion occur. CPIES estimates proved to be excellent at reproducing the full-depth hydrographic fields in this hotspot and captured two anticyclonic eddies (NLWEs) carrying Kuroshio waters during the summer–fall season of 2018.
Similar to previous studies focusing on LWEs generated west of Luzon Island, the NLWEs generated north of Luzon Island were found to be seasonal anticyclonic eddies. During the CPIES observational period, two NLWEs (A and B) were generated north of Luzon Island accompanied by a strong anticyclonic WSC. Eddy A was much larger in size because it absorbed two small anticyclonic eddies nearby. A strong cyclonic eddy separated eddy A and eddy B, which successively moved northwestward across the CPIES array forming an eddy train. These eddies demised northwest of the observational region near 21.5°N and 118.5°E. Because the NLWEs were generated close to the Kuroshio main path, a large amount of Kuroshio water was entrained, as identified by SSH signatures and Argo T/S profiles. Saline water led to significant changes in temperature and salinity in the upper layer of the NLWEs, particularly around the subsurface layer near 150 m. The water column below the 14.8°C isotherm (at a depth of approximately 200 m) inside the NLWEs was still controlled by the SCS waters. The total volume and salt content of Kuroshio waters carried by eddy A was estimated to be 9.74 × 1012 m3 and 3.38 × 1014 kg on 8 August compared to 4.62 × 1012 m3 and 1.59 × 1014 kg for eddy B, respectively. The annual mean transport for historical NLWEs was estimated to be 0.38 Sv, which accounts for 38% of the westward Luzon Strait volume transport during summer (Zhao et al. 2009).
The seasonal monsoon and its associated negative WSC north of Luzon Island provide favorable conditions for the generation of NLWEs in summer, but the mean WSC is not strong enough to drive anticyclonic eddies. The local mean WSC reached its minimum on 15 July, leading to the local highest SSH after 19 days. Historical statistics of anticyclonic eddies generated north of Luzon Island indicate that NLWEs do not occur between November and March, with a mean generation date of 11 July. With the seasonal signals removed, local WSC is also strongly correlated with local SSH (r = −0.48 overall and r = −0.68 for the CPIES observational period). The short lag (3 days) is more consistent with the observation that anticyclonic eddies always appeared after the prevailing local negative WSC, which explains why the mean generation date of NLWEs (11 July) was closer to the date when WSC was minimal (15 July). Thus, we conclude that the local negative intraseasonal WSC is directly responsible for the generation of NLWEs. In other words, the positive seasonal WSC during winter prevents the generation of NLWEs, which are only generated during other seasons under the influence of local strongly negative WSC events. Intraseasonal WSC may also explain the generation of the LWE reported by Yuan et al. (2007), who suggested that the generation mechanism is complicated and could not be fully driven by seasonal winds.
In previous studies, the Kuroshio intrusion into the SCS was considered to be weakest during the summer. The seasonal NLWEs discussed in this paper (Fig. 15) are a new type of Kuroshio water intrusion that occur in summer and have not been previously documented. Because NLWE-induced volume transport likely accounts for almost half of the LST during summer, NLWEs play an important role in stratification west of the Luzon Strait and may also contribute to SCS seasonal circulation. Moreover, as revealed by previous studies, eddies hitting the SCS continental shelf can excite cross-shelf exchange (Wang et al. 2018) and topographic Rossby waves (Wang et al. 2019; Zheng et al. 2021a,b, 2022a,b), and the latter can induce considerable abyssal velocity variation. Further research is now required to better understand the interannual variability of NLWEs and their impacts on deep circulation in the SCS.
Acknowledgments.
The State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography, Ministry of Natural Resources, and the Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) contributed to the work equally and should be regarded as co–first institutional affiliations. This study was supported by the Scientific Research Fund of the Second Institute of Oceanography, MNR (Grants QNYC2102 and JZ2001), the National Natural Science Foundation of China (Grants 41920104006, 41906023, and 41906024), the Project of State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography (SOEDZZ2106 and SOEDZZ2207), the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University (Project No. SL2021MS021), and the Innovation Group Project of the Southern Marine Science and Engineering Guangdong Laboratory, Zhuhai (No. 311020004).
Data availability statement.
Bathymetry data were obtained from ETOPO1 (https://doi.org/10.7289/V5C8276M). Surface geostrophic currents and absolute dynamic topography data were obtained from the CMEMS (http://marine.copernicus.eu/). The moored data analyzed here can be obtained from the corresponding author upon reasonable request. Quality-controlled Argo profiles were retrieved from the Observation and Research Station of the Global Ocean Argo System (Hangzhou), MNR. CTD profiles were obtained from the World Ocean Database (https://www.ncei.noaa.gov/access/world-ocean-database-select/dbsearch.html) and the SCS Open Cruise. The CTD profiles from “SCS Open Cruise” are not publicly available and only the cruise participants are allowed to use the data.
REFERENCES
Bretherton, F., R. E. Davis, and C. B. Fandry, 1976: A technique for objective analysis and design of oceanographic experiments applied to MODE-73. Deep-Sea Res. Oceanogr. Abstr., 23, 559–582, https://doi.org/10.1016/0011-7471(76)90001-2.
Carnes, M. R., 2009: Description and evaluation of GDEM-V3.0. NRL Rep. NRL/MR/7330-09-9165, Naval Research Laboratory, 24 pp.
Caruso, M. J., G. G. Gawarkiewicz, and R. C. Beardsley, 2006: Interannual variability of the Kuroshio intrusion in the South China Sea. J. Oceanogr., 62, 559–575, https://doi.org/10.1007/s10872-006-0076-0.
Chen, G., Y. Hou, X. Chu, P. Qi, and P. Hu, 2009: The variability of eddy kinetic energy in the South China Sea deduced from satellite altimeter data. Chin. J. Oceanol. Limnol., 27, 943, https://doi.org/10.1007/s00343-009-9297-6.
Chen, G., Y. Hou, and X. Chu, 2011: Mesoscale eddies in the South China Sea: Mean properties, spatiotemporal variability, and impact on thermohaline structure. J. Geophys. Res., 116, C06018, https://doi.org/10.1029/2010JC006716.
Donohue, K. A., D. R. Watts, K. Tracey, A. D. Greene, and M. Kennelly, 2010: Mapping circulation in the Kuroshio Extension with an array of current and pressure recording inverted echo sounders. J. Atmos. Oceanic Technol., 27, 507–527, https://doi.org/10.1175/2009JTECHO686.1.
Farris, A., and M. Wimbush, 1996: Wind-induced Kuroshio intrusion into the South China Sea. J. Oceanogr., 52, 771–784, https://doi.org/10.1007/BF02239465.
Hu, J., H. Kawamura, H. Hong, and Y. Qi, 2000: A review on the currents in the South China Sea: Seasonal circulation, South China Sea warm current and Kuroshio intrusion. J. Oceanogr., 56, 607–624, https://doi.org/10.1023/A:1011117531252.
Hu, J., Q. Zheng, Z. Sun, and C.-K. Tai, 2012: Penetration of nonlinear Rossby eddies into South China Sea evidenced by cruise data. J. Geophys. Res., 117, C03010, https://doi.org/10.1029/2011JC007525.
Huang, X., Z. Chen, W. Zhao, Z. Zhang, C. Zhou, Q. Yang, and J. Tian, 2016: An extreme internal solitary wave event observed in the northern South China Sea. Sci. Rep., 6, 30041, https://doi.org/10.1038/srep30041.
Huang, X., Z. Zhang, X. Zhang, H. Qian, W. Zhao, and J. Tian, 2017: Impacts of a mesoscale eddy pair on internal solitary waves in the northern South China Sea revealed by mooring array observations. J. Phys. Oceanogr., 47, 1539–1554, https://doi.org/10.1175/JPO-D-16-0111.1.
Kennelly, M. A., K. L. Tracey, and D. R. Watts, 2007: Inverted echo sounder data processing manual. University of Rhode Island Physical Oceanography Tech. Rep. 2, 87 pp., https://digitalcommons.uri.edu/cgi/viewcontent.cgi?article51001& context5physical_oceanography_techrpts.
Kennelly, M. A., K. A. Donohue, A. Greene, K. L. Tracey, and D. R. Watts, 2008: Inverted echo sounder data report: Kuroshio Extension system study (KESS), April 2004 to July 2006. University of Rhode Island GSO Tech. Rep. 2008-02, 79 pp., https://digitalcommons.uri.edu/cgi/viewcontent.cgi?article=1000&context=physical_oceanography_techrpts.
Le Traon, P.-Y., and R. A. Morrow, 2001: Ocean currents and mesoscale eddies, in satellite altimetry and Earth sciences. A Handbook of Techniques and Applications, L.-L. Fu and A. Cazenave, Eds., Academic Press, 171–215.
Li, L., W. D. Nowlin, and J. Su, 1998: Anticyclonic rings from the Kuroshio in the South China Sea. Deep-Sea Res. I, 45, 1469–1482, https://doi.org/10.1016/S0967-0637(98)00026-0.
Liu, Z. H., Z. Q. Li, S. L. Lu, F. X. Wu, C. H. Sun, and J. P. Xu, 2021: Scattered dataset of global ocean temperature and salinity profiles from the international Argo program. J. Global Change Data Discovery, 5, 312–321, https://doi.org/10.3974/geodp.202104.03.09.
Long, Y., X.-H. Zhu, X. Guo, F. Ji, and Z. Li, 2021: Variations of the Kuroshio in the Luzon Strait revealed by EOF analysis of repeated XBT data and sea-level anomalies. J. Geophys. Res. Oceans, 126, e2020JC016849, https://doi.org/10.1029/2020JC016849.
Mitchell, D. A., W. J. Teague, M. Wimbush, D. R. Watts, and G. G. Sutyrin, 2005: The Dok cold eddy. J. Phys. Oceanogr., 35, 273–288, https://doi.org/10.1175/JPO-2684.1.
Nan, F., H. Xue, F. Chai, L. Shi, M. Shi, and P. Guo, 2011a: Identification of different types of Kuroshio intrusion into the South China Sea. Ocean Dyn., 61, 1291–1304, https://doi.org/10.1007/s10236-011-0426-3.
Nan, F., H. Xue, P. Xiu, F. Chai, M. Shi, and P. Guo, 2011b: Oceanic eddy formation and propagation southwest of Taiwan. J. Geophys. Res., 116, C12045, https://doi.org/10.1029/2011JC007386.
Nan, F., H. Xue, and F. Yu, 2015: Kuroshio intrusion into the South China Sea: A review. Prog. Oceanogr., 137, 314–333, https://doi.org/10.1016/j.pocean.2014.05.012.
Okubo, A., 1970: Horizontal dispersion of floatable particles in the vicinity of velocity singularity such as convergences. Deep-Sea Res. Oceanogr. Abstr., 17, 445–454, https://doi.org/10.1016/0011-7471(70)90059-8.
Park, J.-H., and D. Farmer, 2013: Effects of Kuroshio intrusions on nonlinear internal waves in the South China Sea during winter. J. Geophys. Res. Oceans, 118, 7081–7094, https://doi.org/10.1002/2013JC008983.
Park, J.-H., D. R. Watts, K. L. Tracey, and D. A. Mitchell, 2005: A multi-index GEM technique and its application to the southwestern Japan/East Sea. J. Atmos. Oceanic Technol., 22, 1282–1293, https://doi.org/10.1175/JTECH1778.1.
Park, J.-H., D. R. Watts, K. A. Donohue, and K. L. Tracey, 2012: Comparisons of sea surface height variability observed by pressure-recording inverted echo sounders and satellite altimetry in the Kuroshio Extension. J. Oceanogr., 68, 401–416, https://doi.org/10.1007/s10872-012-0108-x.
Qu, T., J. B. Girton, and J. A. Whitehead, 2006: Deepwater overflow through Luzon Strait. J. Geophys. Res., 111, C01002, https://doi.org/10.1029/2005JC003139.
Ray, R. D., and E. D. Zaron, 2011: Non-stationary internal tides observed with satellite altimetry. Geophys. Res. Lett., 38, L17609, https://doi.org/10.1029/2011GL048617.
Sheremet, V. A., 2001: Hysteresis of a western boundary current leaping across a gap . J. Phys. Oceanogr., 31, 1247–1259, https://doi.org/10.1175/1520-0485(2001)031,1247:HOAWBC.2.0.CO;2.
Sheu, W.-J., C.-R. Wu, and L.-Y. Oey, 2010: Blocking and westward passage of eddies in the Luzon Strait. Deep-Sea Res. II, 57, 1783–1791, https://doi.org/10.1016/j.dsr2.2010.04.004.
Sun, Z., Z. Zhang, B. Qiu, X. Zhang, C. Zhou, X. Huang, W. Zhao, and J. Tian, 2020: Three-dimensional structure and interannual variability of the Kuroshio Loop Current in the northeastern South China Sea. J. Phys. Oceanogr., 50, 2437–2455, https://doi.org/10.1175/JPO-D-20-0058.1.
Tracey, K. L., K. A. Donohue, D. R. Watts, and T. K. Chereskin, 2013: cDrake CPIES data report November 2007 to December 2011. GSO Tech. Rep. 4, 81 pp., https://digitalcommons.uri.edu/cgi/viewcontent.cgi?article=1003&context=physical_oceanography_techrpts.
Wang, D., H. Xu, J. Lin, and J. Hu, 2008: Anticyclonic eddies in the northeastern South China Sea during winter 2003/2004. J. Oceanogr., 64, 925–935, https://doi.org/10.1007/s10872-008-0076-3.
Wang, G., J. Su, and P. C. Chu, 2003: Mesoscale eddies in the South China Sea observed with altimeter data. Geophys. Res. Lett., 30, 2121, https://doi.org/10.1029/2003GL018532.
Wang, G., D. Chen, and J. Su, 2008: Winter eddy genesis in the eastern South China Sea due to orographic wind jets. J. Phys. Oceanogr., 38, 726–732, https://doi.org/10.1175/2007JPO3868.1.
Wang, G., S.-P. Xie, T. Qu, and R. X. Huang, 2011: Deep South China Sea circulation. Geophys. Res. Lett., 38, L05601, https://doi.org/10.1029/2010GL046626.
Wang, Q., L. Zeng, J. Li, J. Chen, Y. He, J. Yao, D. Wang, and W. Zhou, 2018: Observed cross-shelf flow induced by mesoscale eddies in the northern South China Sea. J. Phys. Oceanogr., 48, 1609–1628, https://doi.org/10.1175/JPO-D-17-0180.1.
Wang, Q., and Coauthors, 2019: Energetic topographic Rossby waves in the northern South China Sea. J. Phys. Oceanogr., 49, 2697–2714, https://doi.org/10.1175/JPO-D-18-0247.1.
Watts, D. R., X. Qian, and K. L. Tracey, 2001: Mapping abyssal current and pressure fields under the meandering Gulf Stream. J. Atmos. Oceanic Technol., 18, 1052–1067, https://doi.org/10.1175/1520-0426(2001)018<1052:MACAPF>2.0.CO;2.
Weiss, J., 1991: The dynamics of enstrophy transfer in 2-dimensional hydrodynamics. Physica D, 48, 273–294, https://doi.org/10.1016/0167-2789(91)90088-Q.
Wu, C.-R., and Y.-C. Hsin, 2012: The forcing mechanism leading to the Kuroshio intrusion into the South China Sea. J. Geophys. Res., 117, C07015, https://doi.org/10.1029/2012JC007968.
Xiu, P., F. Chai, L. Shi, H. J. Xue, and Y. Chao, 2010: A census of eddy activities in the South China Sea during 1993–2007. J. Geophys. Res., 115, C03012, https://doi.org/10.1029/2009JC005657.
Yang, H., and Q. Liu, 2003: Forced Rossby wave in the northern South China Sea. Deep-Sea Res. I, 50, 917–926, https://doi.org/10.1016/S0967-0637(03)00074-8.
Yang, Y., and Coauthors, 2019: Eddy-induced transport of saline Kuroshio water into the northern South China Sea. J. Geophys. Res. Oceans, 124, 6673–6687, https://doi.org/10.1029/2018JC014847.
Yuan, D., and Z. Wang, 2011: Hysteresis and dynamics of a western boundary current flowing by a gap forced by impingement of mesoscale eddies. J. Phys. Oceanogr., 41, 878–888, https://doi.org/10.1175/2010JPO4489.1.
Yuan, D., W. Han, and D. Hu, 2007: Anti-cyclonic eddies northwest of Luzon in summer–fall observed by satellite altimeters. Geophys. Res. Lett., 34, L13610, https://doi.org/10.1029/2007GL029401.
Yuan, Y., G. Liao, C. Yang, Z. Liu, H. Chen, and Z.-G. Wang, 2014: Summer Kuroshio Intrusion through the Luzon Strait confirmed from observations and a diagnostic model in summer 2009. Prog. Oceanogr., 121, 44–59, https://doi.org/10.1016/j.pocean.2013.10.003.
Zhang, Z., W. Zhao, J. Tian, and X. Liang, 2013: A mesoscale eddy pair southwest of Taiwan and its influence on deep circulation. J. Geophys. Res. Oceans, 118, 6479–6494, https://doi.org/10.1002/2013JC008994.
Zhang, Z., J. Tian, B. Qiu, W. Zhao, P. Chang, D. Wu, and X. Wan, 2016: Observed 3D structure, generation, and dissipation of oceanic mesoscale eddies in the South China Sea. Sci. Rep., 6, 24349, https://doi.org/10.1038/srep24349.
Zhang, Z., W. Zhao, B. Qiu, and J. Tian, 2017: Anticyclonic eddy sheddings from Kuroshio loop and the accompanying cyclonic eddy in the northeastern South China Sea. J. Phys. Oceanogr., 47, 1243–1259, https://doi.org/10.1175/JPO-D-16-0185.1.
Zhao, R., X.-H. Zhu, and X. Guo, 2017: The impact of monsoon winds and mesoscale eddies on thermohaline structures and circulation patterns in the northern South China Sea. Cont. Shelf Res., 143, 240–256, https://doi.org/10.1016/j.csr.2016.06.009.
Zhao, W., Y.-J. Hou, P. Qi, K.-T. Le, and M.-K. Li, 2009: The effects of monsoons and connectivity of South China Sea on the seasonal variations of water exchange in the Luzon Strait. J. Hydrodyn., 21, 264–270, https://doi.org/10.1016/S1001-6058(08)60144-4.
Zhao, Z., 2014: Internal tide radiation from the Luzon Strait. J. Geophys. Res. Oceans, 119, 5434–5448, https://doi.org/10.1002/2014JC010014.
Zheng, H., C. Zhang, R. Zhao, X.-H. Zhu, Z.-N. Zhu, Z.-J. Liu, and M. Wang, 2021a: Structure and variability of abyssal current in northern South China Sea based on CPIES observations. J. Geophys. Res. Oceans, 126, e2020JC016780, https://doi.org/10.1029/2020JC016780.
Zheng, H., X.-H. Zhu, C. Zhang, R. Zhao, Z.-N. Zhu, and Z.-J. Liu, 2021b: Propagation of topographic Rossby waves in the deep basin of the South China Sea based on abyssal current observations. J. Phys. Oceanogr., 51, 2783–2791, https://doi.org/10.1175/JPO-D-21-0051.1.
Zheng, H., and Coauthors, 2022a: Observation of abyssal circulation to the west of the Luzon Strait, South China Sea. J. Phys. Oceanogr., 52, 2091–2109, https://doi.org/10.1175/JPO-D-21-0284.1.
Zheng, H., and Coauthors, 2022b: Observation of bottom-trapped topographic Rossby waves to the west of the Luzon Strait, South China Sea. J. Phys. Oceanogr., 52, 2853–2872, https://doi.org/10.1175/JPO-D-22-0065.1.
Zu, T. T., D. X. Wang, C. X. Yan, I. Belkin, W. Zhuang, and J. Chen, 2013: Evolution of an anticyclonic eddy southwest of Taiwan. Ocean Dyn., 63, 519–531, https://doi.org/10.1007/s10236-013-0612-6.