XU Anqi , YU Fei , , NAN Feng , REN Qiang

1 Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China

2 University of Chinese Academy of Sciences, Beijing 100049, China

3 Key Laboratory of Ocean Circulation and Waves, Chinese Academy of Sciences, Qingdao 266071, China 4 Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China

5 Marine Dynamic Process and Climate Function Laboratory, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

Abstract Subsurface eddies (SSEs) are common features of the ocean interior. They are particularly abundant in oceanic basins and the vicinity of major intermediate water outfl ows. They are responsible for subsurface transport of mass, heat, and salt. Analysis of high-resolution general circulation model data has revealed the existence of subsurface anticyclonic eddies (SSAEs) and subsurface cyclonic eddies (SSCEs) in the northwestern tropical Pacifi c Ocean. SSEs are abundant east of the Philippines (0°-22°N, 120°E-140°E) and in latitude bands between 9°N-17°N east of 140°E. The composite structure of SSEs was investigated. SSEs had a core at about 400-m water depth and their maximum meridional velocity exceeded 10 cm/s. They exhibited two cores with diff erent salinity polarities in the surface and subsurface. Additionally, spatial distributions of heat transport induced by SSEs in the northwestern tropical Pacifi c were presented for the fi rst time. A net equatorward heat fl ux toward a temperature up-gradient was observed. The analysis of eddy-mean fl ow interactions revealed that the circulation is baroclinically and barotropically unstable at diff erent depths and to diff ering degrees. The energy conversions suggest that both barotropic and baroclinic instabilities are responsible for SSE generation east of the Philippines, whereas baroclinic instability caused by a horizontal density gradient and vertical eddy heat fl ux are important between 9°N and 17°N east of 140°E. Meridional movement of the north equatorial current and the north equatorial undercurrent can contribute to SSE generation in our study region.

Keyword: the northwestern tropical Pacifi c Ocean; subsurface eddies; composite structures; heat transport; energy source

1 INTRODUCTION

The higher resolution of the merged dataset reveals that more than 50% of the variability over much of the World Ocean is accounted for by eddies with amplitudes of 5-25 cm and diameters of 100-200 km (Chelton et al., 2007). Mesoscale eddies can be simply classifi ed into two distinct categories depending on their vertical position (surface or subsurface). There is evidence that surface-intensifi ed eddies, as a ubiquitous phenomenon in the ocean, play an essential role in transferring physical, biological, and chemical properties across large spatial and temporal scales (Shapiro and Meschanov, 1991; Chaigneau et al., 2011; Chelton et al., 2011; Gaube et al., 2013; Zhang et al., 2014; Thomsen et al., 2016). Using satellite data and Argo fl oats, the characteristics and generation mechanisms of surface-intensifi ed eddies in northwestern tropical oceans have been well studied (Kashino et al., 2013; Chen et al., 2015b). Another type of eddy that is largely invisible at the sea surface is known as a subsurface-intensifi ed eddy (SSE). SSEs typically have their core in a subsurface layer and exhibits, by nature, a completely distinct vertical structure from surface-intensifi ed eddies. SSEs are characterized by oceanographic anomalies indicative of foreign water of distant origin. Isopleths of subsurface anticyclonic (cyclonic) eddies tend to dome (bowl) in shallow water depths and bowl (dome) in deeper water (McGillicuddy, 2015).

Fig.1 Depth sections of zonal geostrophic velocity (color shading, m/s) along 140°E (a); topography of the study area in the north western Pacifi c Ocean (b); the spatial distribution of SSAE occurrence (c); the spatial distribution of SSCE occurrence (d)

Unlike surface eddies, which can be characterized from satellite data, the discovery of SSEs is primarily through hydrographic data. Despite these diffi culties in observation, SSEs have been found in many parts of the global ocean, such as Mediterranean water eddies (Meddies) in the North Atlantic Ocean (Isern-Fontanet et al., 2004; Filyushkin et al., 2014; Damien et al., 2017; Radko and Sisti, 2017), Red Sea outfl ow water eddies (Reddies) in the Indian Ocean (Shapiro and Meschanov, 1991), California undercurrent eddies (Cuddies) in the eastern North Pacifi c Ocean (Collins et al., 2013; Pelland et al., 2013), 13°C eddies in the eastern South Pacifi c Ocean (Johnson and McTaggart, 2010; Combes et al., 2015). SSEs were also observed in the Arctic Ocean (Brannigan et al., 2017), the Bay of Bengal (Gordon et al., 2017) and in Drake Passage (Brearley at al., 2014).

The northwestern tropical Pacifi c is characterized by a complicated ocean circulation structure. The North Equatorial Current (NEC) bifurcates into the northward-fl owing Kuroshio and southward-fl owing Mindanao currents (MC) as it approaches the coast of the Philippines (Qu et al., 2012; Schonau et al., 2015). MC carries the north Pacifi c tropical water and the north Pacifi c intermediate water. The main body of the NEC is confi ned to water shallower than 500 m. The MC veers eastward at the southern tip of Mindanao Island near 5°N to form the North Equatorial Countercurrent (NECC) (Qiu et al., 2015). Below the surface currents, the Luzon Undercurrent (LUC) that fl ows equatorward beneath the Kuroshio Current, the poleward-fl owing Mindanao Undercurrent (MUC) beneath the MC and eastwardfl owing North Equatorial Undercurrent (NEUC) below the NEC form a highly variable subsurface current system (Qu et al., 1997; Qu and Lukas, 2003). Besides, the northwestern tropical Pacifi c is a region where interhemispheric waters converge. To clarify the vertical structures of the currents, we plot in Fig.1a the latitude-depth sections of zonal geostrophic velocity U( y, z) along 140°E based on cruise observations conducted by R/V Kexue in April 2017. A Sea-Bird SBE 911 conductivity-temperature-depth (CTD) with dual temperature and conductivity sensors was used to measure the conductivity and temperature profi les at 24 Hz. The accuracy of conductivity was 0.0003 S/m, and the temperature was 0.001°C for CTD. The horizontal geostrophic velocities can be derived as the following:

where Δ φ is dynamic height that can be calculated from temperature and salinity, A and B denote CTD stations, L is the distance between the stations, g is gravitational acceleration, and f is the Coriolis parameter.

Recently, mooring observations showed that SSEs also occasionally observed in the northwestern tropical Pacifi c Ocean. At least two groups of SSEs near the Philippines coast were detected: one originating from the southeast and the other from the east (Chiang et al., 2015). Using a subsurface mooring (4°N, 143°E), Zhang et al. (2017) reported a subthermocline eddy that enhanced shear instability. An extra-large subsurface anticyclonic eddy (horizontal scale of 470 km) was detected east of Taiwan, China, and was found to carry Subtropical Mode Water westward over 1 500 km (Nan et al., 2017). A deep baroclinic eddy induced by a surface cyclonic eddy is detected in the northwestern South China Sea (Chen et al., 2015a). Later, utilizing model data, SSE characteristics in the northwestern Pacifi c Ocean, such as distribution, size, amplitude, polarity and longevity, have been reported and the results show SSEs are strongly nonlinear during most of their lifetime, which implies that they can trap interior fl uid during translation (Xu et al., 2019).

However, compared with surface-intensifi ed eddies, our understanding of SSEs in the northwestern Pacifi c Ocean remains observationally fragmentary and theoretically defective due to spatial and temporal limitations of in situ subsurface observations. Although the main horizontal distributions and kinematic properties of SSEs have been investigated during the last decade, very little is known about their vertical structure, their impact on heat transport, and their generation mechanisms. The rest of this paper is structured as follows. Section 2 describes the data and methodology used. In Section 3, we analyze the characteristics of the vertical structures of the SSE (Section 3.1), heat transport induced by SSE (Section 3.2) and energy source of SSE (Section 3.3). The paper ends with a summary of the results (Section 4).

2 DATA AND METHOD

2.1 Model description

The Oceanic General Circulation Model for the Earth Simulator (OFES) was used in this study. It is based on the Modular Ocean Model V.3 with a spherical grid in the horizontal and z-coordinates in the vertical. Its domain extends from 75°S to 75°N, excluding the Arctic region, with 0.1° horizontal grid spacing. The vertical spacing varies from 5 m at the surface to 330 m near the bottom. After a 50-year climatological spin-up, the model was forced by daily surface wind stress, heat fl ux and salinity fl ux based on National Centers for Environmental Prediction reanalysis products. The readers are referred to Masumoto et al. (2004) and Sasaki et al. (2008) for more details about the confi guration and evaluation of this model. The OFES products have been extensively analyzed by a number of early studies and the results have demonstrated the promising capability of the model in representing realistic variability of diff erent spatial and temporal scales in the ocean (Dutrieux, 2009; Qu et al., 2012; Kashino et al., 2013; Qiu et al., 2015). The validation of the model data and the eddy identifi cation procedure used in this paper can be found in Xu et al. (2019). Comparison of model data and altimeter observations indicated that the model data and our detection algorithm could satisfactorily reproduce eddy activities.

2.2 Subsurface eddy identifi cation

Our detection method is based on Okubo-Weiss parameter contours (Isern-Fontanet et al., 2006). The Okubo-Weiss parameter ( W) is defi ned as:

Because eddies have a deep vertical structure, a deep surface eddy can have signatures at deeper depths (＞400 m), and conversely a few SSEs can have weak signatures at the surface. To avoid double counting, only those eddies with their entire body deeper than 45 m were counted as SSEs (i.e. those that could not be detected in the upper eight layers of the model data).

3 RESULT

3.1 Composite structure of subsurface eddies

Figure 1c & d shows the spatial occurrence and dominance of subsurface anticyclonic eddies (SSAEs) and subsurface cyclonic eddies (SSCEs). Here, eddy frequency was defi ned as the percentage of time that a point was located within a cyclonic or anticyclonic eddy. For example, eddy frequency at a specifi c location was considered about 10% if occupied by eddies for 10 out of 100 days. In our study region, SSEs were mainly distributed between 9°N and 17°N in open ocean and east of the Philippines coast (Fig.1c & d). In addition, SSAEs were observed more frequently than SSCEs to the east of the Philippines, which is consistent with previous observational studies (Combes et al., 2015; Song et al., 2017).

To delineate the vertical structure of SSEs, we generated composite maps of eddy anomalies for all detected SSEs. The eddy-induced anomaly for a fi eld at a given time and depth in the analysis was calculated by subtracting the horizontal average of the fi eld over a box fi eld of 5°×5° bin, centered around the eddy. The anomalies are in the depth surface. All identifi ed SSEs were selected for the reconstruction by averaging the snapshots, regardless of their lifespan. Figure 2 shows the vertical section of the composite SSE along the zonal direction at Δ y=0. The positive temperature anomaly core of SSAEs was located at around 650 m, deeper than the SSCE negative core, which was found at 500 m (Fig.2a & b). Temperature anomalies greater than 0.3°C extended from 600 m to 800 m for SSAEs and 400 m to 800 m for SSCEs. Interestingly for SSAEs was the occurrence of a clear, near zero temperature anomaly that coincided with maximum velocity at 400-m water depth. The water in a SSAE was colder than surrounding water above 400 m and warmer below 400 m. Negative temperature anomalies prevailed throughout a SSCE. Both SSAEs and SSCEs showed two salinity anomaly cores of diff erent polarities (Fig.2c & d). SSAEs displayed positive anomalies in their cores and negative anomalies below their cores. The opposite was true for SSCEs; they exhibited negative and positive anomalies in and below their cores, respectively.

The mean vertical structures of meridional velocity associated with SSAEs and SSCEs along the zonal section at Δ y=0 are shown in Fig.2e & f. As expected, maximum speeds prevailed in subsurface layers. In the subsurface, the meridional velocity varied from a few cm/s near the center to 10-13 cm/s approximately halfway between the composite eddy center and its edge. The maxima velocity in the very surface maybe the result of surface eddies which is right above the subsurface eddies. Note that the SSAE and SSCE surveyed by Nan et al. (2017) in east of the Philippines had a maximum swirl speed of 0.4 m/s at ～400 m depth and were located between 200 m and 800 m. In agreement with Nan et al. (2017), the SSE identif ication method used here revealed that both SSAEs and SSCEs have an obvious maximum velocity at ～400 m depth along the zonal section at Δ y=0. These vertical velocity structures give us conf idence in our identif ication method.

The potential vorticity (PV) and PV anomaly (PVA) are particularly suited to identify and describe SSEs (Thomas, 2008; Zhang et al., 2017). The PV is def ined as:

3.2 Heat transport of subsurface eddies

Fig.2 Vertical sections of potential temperature anomalies (°C) (a, b), salinity anomalies (c, d), and meridional speed (cm/s) (e, f) across composite anticyclonic and cyclonic SSEs

Fig.3 Zonal section of PV anomalies (/(m·s)) through the eddy center for SSAE (a) and SSCE (b)

Fig.4 Zonal section of normalized relative vorticity ( ζ/ f) through the eddy center for SSAE (a) and SSCE (b)

Theoretical and observational analyses indicated that mesoscale eddies can transport heat, salt and sea water constituents through individual movements (Volkov et al., 2008; Dong et al., 2014; Zhang et al., 2014). The amount of water trapped inside an eddy depends on the ratio of its average rotational speed to its translation speed, also known as the nonlinearity parameter. An eddy is nonlinear when this ratio exceeds one and maintains a coherent structure as it propagates (Chelton et al., 2007, 2011). Accordingly, we analyzed the nonlinearity of SSEs in the northwestern Pacifi c Ocean. The modeling results showed that most of the observed SSEs were nonlinear and about 89% (37%) of tracked SSEs were nonlinear for at least half (all) of their lifetime (Xu et al., 2019). The composite vertical SSE structure (see Section 3.1) implied strong and complex temperature anomalies. This type of anomaly tends to move with a SSE owing to the trapping of water parcels caused by nonlinearity.

We estimated heat transport induced by SSE movement. For a single SSE, the heat anomaly ( Q) was calculated as:

where ρ is the density, Cpis the heat capacity (4 200 J/(kg·°C)), T′ is temperature anomaly, r is the eddy radius within a specifi c layer and z is depth layer thickness within a SSE. The vertical integral was calculated by discrete summation over the interpolated vertical levels. In a SSE, Q is dominated by the deeper anomalies that generally span a larger depth range than the shallower anomalies. We set up an average bin with 2.5° meridional width ( Dlat) and 2.5° zonal length ( Dlon) around each grid point. The heat transport carried by every eddy identifi ed in the OFES was Q· Cxwhere Cxis the zonal propagation speed of individual eddies. We accumulated the total SSEinduced transport within each bin for whole periods from 2003 to 2013 and divided it by Dlon, Dlat, and the number of snapshots ( N). The average eddy-induced zonal heat transport (unit: W/m) is

Similar to Hx, we also computed SSE-induced meridional transport

This method is similar to the one used by Zhang et al. (2014) to study eddy-induced mass transport. However, we replaced trapped volume with heat anomalies in their equations to calculate heat transport.

Fig.5 Distribution of zonal (a) and meridional (b) heat transport (W/m) of heat trapped by SSAEs

Fig.6 Distribution of zonal (a) and meridional (b) heat transport (W/m) of heat trapped by SSCEs

Figures 5 and 6 show the distribution of heat transport by SSAEs and SSCEs in our study region. The heat transport patterns were closely associated with eddy distribution and movement. The heat transport was large between 9°N-17°N where SSEs occur frequently. The distribution of meridional heat transport was on the order of 103W/m, which is much smaller than the zonal heat transport (104W/m). This is because of diff erences in propagation characteristics. Xu et al. (2019) reported that most SSEs move westward with slight equatorward or poleward defl ections, and about 10% and 6% of SSAEs and SSCEs, respectively, propagated purely zonally (±1°). The mean zonal velocity from a SSE was about 6.6 cm/s, while it was only about 0.3 cm/s in the meridional direction. In the zonal direction, both SSAEs and SSCEs had negative zonal velocity (westward propagation). With diff ering thermocline displacement tendencies, the meridional integrated zonal heat transport by SSAEs and SSCEs tended to off set each other to some extent. Hence, the total zonal SSE transport was generally smaller than by either type of SSE alone (Fig.7a). However, the zonally integrated meridional heat transport of SSAEs and SSCEs combined to give larger values than from either type alone (between 8°N-22°N) (Fig.7b). The total meridional heat transport peaked at 5°N (～6×109W) and 20°N (～8×109W).

3.3 Energy source

The energy in the gl obal ocean is mainly contained as mean kinetic energy (MKE), eddy kinetic energy (EKE), mean available potential energy (MPE) and eddy available potential energy (EPE) (von Storch et al., 2012). An EKE source comprises baroclinic instability (BC) (MPE→EPE→EKE) in which the generation of EKE occurs at the expense of EPE and barotropic instability (BT) (MKE→EKE) that occurs at the expense of MKE. To determine why SSEs frequently occur between 9°N and 17°N and east of the Philippines, and to understand the mechanism of SSE generation in the northwestern tropical Pacifi c Ocean, we estimated the energy conversion between various components of the total energy. EKE and MKE, the kinetic energies of time-averaged and timevarying circulation, are well defi ned and can be accurately estimated. However, MPE and EPE, the available potential energies of time-averaged and time-varying circulation, cannot be calculated exactly because they are dependent of the choice of the reference density. Despite these limitations, we can examine the energy transfer rate among diff erent reservoirs.

Given that intrinsic instability processes of ocean currents lead to eddy generation, we fi rst calculated the BT and BC conversion terms following (Oey, 2008).

where N2is the square buoyancy frequency, x and y are the longitudinal and latitudinal coordinates, respectively, g is the gravitational acceleration and ρ0is the reference density 1 025 kg/m3. Additionally, the prime symbol denotes the deviation from the mean value and the overbar denotes time averaging over 60 days. Here 60-days are suffi ciently long to cover the expected periods of instability growths, yet not so long that slow variations of the fl ow fi elds will be lumped as “instability”. To estimate the deviation squares, we used:

Fig.7 Meridionally integrated zonal heat transport induced by SSEs as a function of longitude (a); zonal integrated zonal meridional induced by SSEs as a function of latitude (b)

Figure 8 shows the energy conversion from MKE to EKE associated with BT. A positive BT indicates that horizontal eddy fl ux follows the mean momentum gradient where MKE is converted to EKE. A negative BT represents eddy decay that accelerates mean fl ow where EKE is converted to MKE. The maximal values of BT were south of 7°N and east of the Philippines in the surface layer, which is due to the high horizontal shear energy of western boundary currents. In the subsurface, BT was characterized by unique boundary intensifi cation behavior.

Fig.8 Horizontal distribution of BT (W/m3) at 20 m (a), 100 m (b), 200 m (c), and 400 m (d)

The transfer of MPE to EPE through baroclinic instabilities of mean fl ow is represented by the abbreviation BC (Fig.9). A positive BC indicates horizontal eddy density fl ux against the direction of the mean density gradient where MPE→EPE. Conversely, a negative BC indicates down gradient fl ux where EPE→MPE. The largest conversion from MPE to EPE occurred in the upper ocean (Fig.9a). The widespread positive BC values at 20 m are consistent with the energy fl ow path of EPE along major boundary currents and frontal zones. In the subsurface layers, the most noticeable feature of BC was the prevailing positive value between 9°N and 17°N (Fig.9).

Fig.9 Horizontal distribution of BC (W/m3) at 20 m (a), 100 m (b), 200 m (c), and 400 m (d)

As shown in a depth section at 143°E, the highest values exceeded 2×10-5W/m3and 1.3×10-4W/m3for BC and PKC, respectively. Positive BC and PKC values were located roughly at 9°N-17°N and 9°N-15°N, respectively, and trend to the north with increasing depth (Fig.11). Negative conversion values were below the positive patches and the interface between them deepens northward from 200 m at 9°N to 600 m at 15°N. These depths cover the interface between the surface westward-fl owing NEC and subsurface eastward-fl owing NEUC. The northward shift is associated with meridional movement of NEC bifurcation, which is depth dependent. On average annually, the shift occurs from about 13.3°N near the surface to north of 20°N at depths around 1 000 m (Qu and Lukas, 2003). Note that when the NEC moves north or south, there can be accompanied meridional movement of the NEC and NEUC interface. Thus, meridional movement of the NEC/NEUC can contribute to SSE generation in the study region. To estimate the robustness of our results and study possible SSE energy sources, the vertical profi le of mean energy conversion was calculated for the area east of the Philippines coast (region 1: 0°-22°N, 120°E-140°E) and in latitude bands between 9°N-17°N east of 140°E (region 2), where SSEs are frequently observed.

In region 1, all of the energy conversions showed maximum values in the surface. The BT was dominant in the upper 300 m because of the strong horizontal shear energy from the NECC. With increases in depth, BT decreased drastically. However, it remained an important source of subsurface EKE below 600 m, where it was the largest of the three conversion parameters. The mean PKC reached its maximum value at ～400 m and changed around 600 m with positive values in the upper water column and negative values at deeper depths. The mean energy transfer from MPE to EPE in the area decreased gradually with depth. In region 2, BT was nearly zero from about 100 m (Fig.12), implying that little MKE is converted to EKE and BT is not a stable source of SSEs in this region. Both BC and PKC showed maximum values in the subsurface at about 300 m. In conclusion, both BT and BC were responsible for SSE generation in region 1, whereas BC caused by a horizontal density gradient and vertical eddy heat fl ux were important in region 2.

Fig.10 Horizontal distribution of PKC (W/m 3) at 20 m (a), 100 m (b), 200 m (c), and 400 m (d)

Fig.11 Vertical structure of BC (a) and PKC (b) along 143°E (W/m 3)

4 SUMMARY

Fig.12 Vertical profi le of mean energy conversions BC, BT, and PKC (10 -5 W/m3) averaged over region 1 (0°-22°N, 120°E-140°E) (a) and region 2 (9°N-17°N, 140°E-180°E) (b)

SSEs are abundant in the northwestern tropical Pacifi c Ocean. Based on high resolution OFES data, a method based on the Okubo-Weiss parameter was used to detect SSEs. This study analyzed the mean vertical structure of both SSAEs and SSCEs. We calculated the zonal heat and meridional transport induced by SSEs. Possible energy sources of SSEs were also investigated. The statistical results revealed high occurrences of SSEs in two regions within our study domain, one east of the Philippines and the other between 9°N and 17°N.

This study characterized the mean vertical structure of SSEs in northwestern tropical Pacifi c Ocean, where several SSEs were previously reported from Argo fl oat or hydrographic data. The composite structures indicated that the maximum meridional SSE velocity exceeded 10 cm/s centered at ～400 m. Our results indicated that SSE signals could penetrate to 1 200 m depth. SSAEs showed negative anomalies in the upper 400 m in the eastern part of the study area and strong and extensive positive anomalies from 500 to 1 000 m. However, temperature anomalies were negative throughout the water column of a SSCE.

Previous studies have shown that SSEs have strong nonlinearity, implying they have considerable impacts on heat transport within the subsurface, especially in regions with abundant SSEs. For the fi rst time, we presented detailed spatial distributions of heat transport by SSE movements in the northwestern tropical Pacifi c Ocean. The meridionally integrated zonal heat transport for SSAEs and SSCEs largely compensate for each other, resulting in weak zonal heat transport overall. The zonally integrated meridional heat transport of total SSEs was on the order of 1010W, which is comparable in magnitude to meridional integrated zonal heat transport. A net equatorward heat fl ux toward a temperature upgradient was observed in our study region, owing to preferential equatorward movement of SSAEs (no obvious tendency of defl ection for SSCE) and their associated large vertical trapping depths.

We also investigated the sources of energy in our study region. As an indicator of energy transfer between MKE and EKE, BT appeared mainly along the western boundaries, especially within the NECC. When integrating the energy conversion terms in region 1, we determined that BT drew 6.08×109W from MKE to EKE, and the BC converted 3.02×109W to EKE through the MPE→EPE→EKE pathway. The mean of the three conversion terms indicated that BT showed maximum value in surface layers and decreased drastically with increasing depth in region 1. In subsurface layers (300 m to 1 000 m), both BT and BC were very important. In region 2, BT was weak and confi ned to the upper 100 m; thus, it was not a stable source for SSEs. The energy was mainly transferred from mean fl ow to eddy fi elds through the MPE→EPE→EKE pathway. The BC was mainly caused by the westward-fl owing NEC at the surface and eastward-fl owing NEUC in the subsurface. In both regions 1 and 2, with negative PKCs below ～600 m depth, the eddy fi elds were found to transfer energy to the mean fl ow through a baroclinic inverse pathway of EKE→EPE.

The results of this study highlight the importance of SSEs in redistributing heat in the northwestern tropical Pacifi c Ocean, and suggest that SSEs might play an important role in modifying intermediate water mass properties. However, this needs to be confi rmed by more focused observations in the future. In particular, this study off ers a useful basis for future rigorous studies on SSEs in the Pacifi c Ocean.

5 DATA AVAILABILITY STATEMENT

All data generated and analyzed during this study are available from the corresponding author upon reasonable request.

6 ACKNOWLEDGMENT

The OFES data for this paper are available from the Asian-Pacifi c Data-Research Center (http://apdrc.soest.hawaii.edu/) and we are grateful to H. SASAKI and colleagues from the Earth Simulator for assistance in processing the OFES outputs.

References

Brannigan L, Johnson H, Lique C, Nycander J, Nilsson J. 2017. Generation of subsurface anticyclones at Arctic surface fronts due to a surface stress. Journal of Physical Oceanography, 47(11): 2 653-2 671, https://doi.org/10. 1175/jpo-d-17-0022.1.

Brearley J A, Sheen K L, Naveira Garabato A C, Smeed D A, Speer K G, Thurnherr A M, Meredith M P, Waterman S. 2014. Deep boundary current disintegration in Drake Passage. Geophysical Research Letters, 41(1): 121-127, https://doi.org/10.1002/2013GL058617.

Chaigneau A, Le Texier M, Eldin G, Grados C, Pizarro O. 2011. Vertical structure of mesoscale eddies in the eastern South Pacifi c Ocean: a composite analysis from altimetry and Argo profi ling fl oats. Journal of Geophysical Research, 116(C16): C11025, https://doi.org/10.1029/2011jc007134.

Chelton D B, Schlax M G, Samelson R M. 2011. Global observations of nonlinear mesoscale eddies. Progress in Oceanography, 91(2): 167-216, https://doi.org/10.1016/j.pocean.2011.01.002.

Chelton D B, Schlax M G, Samelson R M, de Szoeke R A. 2007. Global observations of large oceanic eddies. Geophysical Research Letters, 34(15): C15606, https://doi.org/10.1029/2007gl030812.

Chen G X, Wang D X, Dong C M, Zu T T, Xue H J, Shu Y Q, Chu X Q, Qi Y Q, Chen H. 2015a. Observed deep energetic eddies by seamount wake. Scientifi c Reports, 5: 17416, https://doi.org/10.1038/srep17416.

Chen L J, Jia Y L, Liu Q Y. 2015b. Mesoscale eddies in the Mindanao Dome region. Journal of Oceanography, 71(1): 133-140, https://doi.org/10.1007/s10872-014-0255-3.

Chiang T L, Wu C R, Qu T D, Hsin Y C. 2015. Activities of 50-80 day subthermocline eddies near the Philippine coast. Journal of Geophysical Research, 120(5): 3 606-3 623, https://doi.org/10.1002/2013jc009626.

Collins C A, Margolina T, Rago T A, Ivanov L. 2013. Looping RAFOS fl oats in the California current system. Deep Sea Research Part II: Topical Studies in Oceanography, 85: 42-61, https://doi.org/10.1016/j.dsr2.2012.07.027.

Combes V, Hormazabal S, Di Lorenzo E. 2015. Interannual variability of the subsurface eddy fi eld in the Southeast Pacifi c. Journal of Geophysical Research, 120(7): 4 907-4 924, https://doi.org/10.1002/2014jc010265.

Damien P, Bosse A, Testor P, Marsaleix P, Estournel C. 2017. Modeling postconvective submesoscale coherent vortices in the northwestern Mediterranean Sea. Journal of Geophysical Research, 122(12): 9 937-9 961, https://doi.org/10.1002/2016jc012114.

Dong C M, McWilliams J C, Liu Y, Chen D K. 2014. Global heat and salt transports by eddy movement. Nature Communications, 5: 3 294, https://doi.org/10.1038/ncomms4294.

Dutrieux P. 2009. Tropical Western Pacifi c Currents and the Origin of Intraseasonal Variability Below the Thermocline. University of Hawaii at Manoa, Honolulu, USA, 150p.

Filyushkin B N, Sokolovskiy M A, Kozhelupova N G, Vagina I M. 2014. Lagrangian methods for observation of intrathermocline eddies in ocean. Oceanology, 54(6): 688-694, https://doi.org/10.1134/s0001437014050051.

Gaube P, Chelton D B, Strutton P G, Behrenfeld M J. 2013. Satellite observations of chlorophyll, phytoplankton biomass, and Ekman pumping in nonlinear mesoscale eddies. Journal of Geophysical Research, 118(12): 6 349-6 370, https://doi.org/10.1002/2013jc009027.

Gordon A L, Shroyer E, Murty V S N. 2017. An intrathermocline eddy and a tropical cyclone in the Bay of Bengal. Scientifi c Reports, 7: 46218, https://doi.org/10.1038/srep46218.

Isern-Fontanet J, Font J, García-Ladona E, Emelianov M, Millot C, Taupier-Letage I. 2004. Spatial structure of anticyclonic eddies in the Algerian basin (Mediterranean Sea) analyzed using the Okubo-Weiss parameter. Deep Sea Research Part II: Topical Studies in Oceanography, 51(25-26): 3 009-3 028, https://doi.org/10.1016/j.dsr2. 2004.09.013.

Isern-Fontanet J, García-Ladona E, Font J. 2006. Vortices of the Mediterranean Sea: an altimetric perspective. Journal of Physical Oceanography, 36(1): 87-103, https://doi.org/10.1175/jpo2826.1.

Johnson G C, McTaggart K E. 2010. Equatorial pacifi c 13°C water eddies in the eastern subtropical South Pacifi c Ocean. Journal of Physical Oceanography, 40(1): 226-236, https://doi.org/10.1175/2009jpo4287.1.

Kashino Y, Atmadipoera A, Kuroda Y, Lukijanto. 2013. Observed features of the Halmahera and Mindanao Eddies. Journal of Geophysical Research, 118(12): 6 543-6 560, https://doi.org/10.1002/2013jc009207.

Kurian J, Colas F, Capet X, McWilliams J C, Chelton D B. 2011. Eddy properties in the California Current system. Journal of Geophysical Research, 116(C8): C08027, https://doi.org/10.1029/2010jc006895.

Magalh?es F C, Azevedo J L L, Oliveira L R. 2017. Energetics of eddy-mean fl ow interactions in the Brazil current between 20°S and 36°S. Journal of Geophysical Research, 122(8): 6 129-6 146, https://doi.org/10.1002/2016jc012609.

Masumoto Y, Sasaki H, Kagimoto T, Komori N, Ishida A, Sasai Y, Miyama T, Motoi T, Mitsudera H, Takahashi K, Sakuma H, Yamagata T. 2004. A fi fty-year eddy-resolving simulation of the world ocean—preliminary outcomes of OFES (OGCM for the Earth simulator). Journal of the Earth Simulator, 1: 35-56.

McGillicuddy Jr D J. 2015. Formation of intrathermocline lenses by eddy-wind interaction. Journal of Physical Oceanography, 45(2): 606-612, https://doi.org/10.1175/jpo-d-14-0221.1.

Nan F, Yu F, Wei C J, Ren Q, Fan C H. 2017. Observations of an extra-large subsurface anticyclonic eddy in the northwestern Pacifi c subtropical gyre. Journal of Marine Science: Research & Development, 7: 234, https://doi.org/10.4172/2155-9910.1000234.

Oey L Y. 2008. Loop current and deep eddies. Journal of Physical Oceanography, 38(7): 1 426-1 449, https://doi.org/10.1175/2007jpo3818.1.

Pelland N A, Eriksen C C, Lee C M. 2013. Subthermocline eddies over the Washington continental slope as observed by seagliders, 2003-09. Journal of Physical Oceanography, 43(10): 2 025-2 053, https://doi.org/10. 1175/jpo-d-12-086.1.

Qiu B, Rudnick D L, Cerovecki I, Cornuelle B D, Chen S M, Schonau M C, McClean J L, Gopalakrishnan G. 2015. The Pacifi c north equatorial current: new insights from the origins of the Kuroshio and Mindanao Currents (OKMC) Project. Oceanography, 28(4): 24-33, https://doi.org/10.5670/oceanog.2015.78.

Qu T D, Chiang T L, Wu C R, Dutrieux P, Hu D X. 2012. Mindanao current/undercurrent in an eddy-resolving GCM. Journal of Geophysical Research, 117(C6): C06026, https://doi.org/10.1029/2011jc007838.

Qu T D, Kagimoto T, Yamagata T. 1997. A subsurface countercurrent along the east coast of Luzon. Deep Sea Research Part I: Oceanographic Research Papers, 44(3): 413-423, https://doi.org/10.1016/s0967-0637(96)00121-5.

Qu T D, Lukas R. 2003. The bifurcation of the North Equatorial Current in the Pacifi c. Journal of Physical Oceanography, 33(1): 5-18, https://doi.org/10.1175/1520-0485(2003) 033＜0005:tbotne＞2.0.co;2.

Radko T, Sisti C. 2017. Life and demise of intrathermocline mesoscale vortices. Journal of Physical Oceanography, 47(12): 3 087-3 103, https://doi.org/10.1175/jpo-d-17-0044.1.

Sasaki H, Nonaka M, Masumoto Y, Sasai Y, Uehara H, Sakuma H. 2008. An eddy-resolving hindcast simulation of the quasiglobal ocean from 1950 to 2003 on the Earth Simulator. In: High Resolution Numerical Modelling of the Atmosphere and Ocean, Hamilton K and Ohfuchi W eds. Chapter 10, pp.157-185, Springer, New York.

Schonau M C, Rudnick D L, Cerovecki I, Gopalakrishnan G, Cornuelle B D, McClean J L, Qiu B. 2015. The Mindanao current: mean structure and connectivity. Oceanography, 28(4): 34-45, https://doi.org/10.5670/oceanog.2015.79.

Shapiro G I, Meschanov S L. 1991. Distribution and spreading of Red Sea water and salt lens formation in the northwest Indian Ocean. Deep Sea Research Part A. Oceanographic Research Papers, 38(1): 21-34, https://doi.org/10. 1016/0198-0149(91)90052-h.

Song L, Li Y L, Liu C Y, Wang F. 2017. Subthermocline anticyclonic gyre east of Mindanao and its relationship with the Mindanao Undercurrent. Chinese Journal of Oceanology and Limnology, 35(6): 1 303-1 318, https://doi.org/10.1007/s00343-017-6111-8.

Thomas L N. 2008. Formation of intrathermocline eddies at ocean fronts by wind-driven destruction of potential vorticity. Dynamics of Atmospheres and Oceans, 45(3-4): 252-273, https://doi.org/10.1016/j.dynatmoce.2008.02.002.

Thomsen S, Kanzow T, Krahmann G, Greatbatch R J, Dengler M, Lavik G. 2016. The formation of a subsurface anticyclonic eddy in the Peru-Chile Undercurrent and its impact on the near-coastal salinity, oxygen, and nutrient distributions. Journal of Geophysical Research, 121(1): 476-501, https://doi.org/10.1002/2015jc010878.

Volkov D L, Lee T, Fu L L. 2008. Eddy-induced meridional heat transport in the ocean. Geophysical Research Letters, 35(20): L20601, https://doi.org/10.1029/2008gl035490.

von Storch J S, Eden C, Fast I, Haak H, Hernández-Deckers D, Maier-Reimer E, Marotzke J, Stammer D. 2012. An estimate of the Lorenz energy cycle for the world ocean based on the° STORM/NCEP simulation. Journal of Physical Oceanography, 42(12): 2 185-2 205, https://doi.org/10.1175/jpo-d-12-079.1.

Xu A Q, Yu F, Nan F. 2019. Study of subsurface eddy properties in northwestern Pacifi c Ocean based on an eddy-resolving OGCM. Ocean Dynamics, 69(4): 463-474, https://doi.org/10.1007/s10236-019-01255-5.

Yang H Y, Wu L X, Liu H L, Yu Y Q. 2013. Eddy energy sources and sinks in the South China Sea. Journal of Geophysical Research, 118(9): 4 716-4 726, https://doi.org/10.1002/jgrc.20343.

Zhang Z G, Wang W, Qiu B. 2014. Oceanic mass transport by mesoscale eddies. Science, 345(6194): 322-324, https://doi.org/10.1126/science.1252418.

Zhang Z G, Zhang Y, Wang W. 2017. Three-compartment structure of subsurface-intensifi ed mesoscale eddies in the ocean. Journal of Geophysical Research, 122(3): 1 653-1 664, https://doi.org/10.1002/2016jc012376.