1. Introduction
Mesoscale eddies are energetic circulations with spatial scales on the order of 100 km and lifetimes spanning from weeks to months (Chelton et al. 2007; Ferrari and Wunsch 2009). They are ubiquitous in global oceans and are mainly generated by instabilities of large-scale ocean circulation (Gill et al. 1974; Chelton et al. 2011). Near-inertial waves (NIWs) are internal gravity waves with a frequency close to the local inertial frequency f and are energetic motions at high frequencies (Garrett and Munk 1979; Ferrari and Wunsch 2009). NIWs are characterized by circularly polarized velocities and strong shears. The lateral spatial scales of NIWs range from tens to hundreds of kilometers, and are comparable to mesoscale eddies, but their vertical wavelengths are relatively short, typically 100–400 m (Alford et al. 2016). They are usually inspired by the resonant wind forcing and therefore commonly appear in the upper ocean (D’Asaro et al. 1995; Chen et al. 2013; Pallàs-Sanz et al. 2016). Once generated in the surface layer, NIWs propagate both equatorward and downward (Gill 1984; Garrett 2001). During propagation, NIWs transform energy to waves with smaller vertical scales and ultimately break, triggering turbulent mixing in the ocean interior (Alford 2003; Whalen et al. 2020). Therefore, NIWs are supposed to be one of the main drivers of the diapycnal mixing in the interior ocean (Munk and Wunsch 1998; MacKinnon et al. 2017).
Relatively slow group velocities allow NIWs likely to interact with mesoscale eddies (Weller 1982). Mesoscale vorticity could modify the propagation pathways of NIWs by changing their intrinsic frequencies. Using the ray tracing approach based on the WKB approximation, Kunze (1985) found that the relative vorticity ζ shifts the lower bound of internal wave band from f to an effective Coriolis frequency, feff = f + ζ/2, by which the negative vorticity [e.g., the anticyclonic warm-core eddies (WEs)] traps the incident NIWs while the positive vorticity [e.g., the cyclonic cold-core eddies (CEs)] expels them. The trapped NIWs are generally reflected off the turning points in the horizontal direction and stall at the critical layers in the vertical direction, where their amplitudes increase and vertical wavelengths shrink. The critical layer behaviors of NIWs lead to strong vertical shears and give rise to intense turbulence owing to shear instabilities (Kunze 1986). A lot of the trapped NIW energy was found to have been lost to turbulence at the base of the eddy core in the WE (Kunze et al. 1995). However, the WKB approximation assumes the NIW scales are smaller than the mesoscale eddies, which makes it inapplicable to the synoptic-scale NIWs initially generated by winds (D’Asaro et al. 1995; D’Asaro 1995). To overcome this limitation, Young and Ben Jelloul (1997) developed a framework that circumvented the scale separation assumption as a supplement to Kunze’s (1985) theory, which predicts that mesoscale vorticity effectively distorts the near-inertial motions at eddy scales and subsequently accelerate their downward dispersion.
Following the above theoretical studies, several numerical studies have focused on the modification of NIW propagation by mesoscale eddies (e.g., Lee and Niiler 1998; Zhai et al. 2005), which emphasize the trapping of NIWs by the WEs. The WEs efficiently drain near-inertial energy from the surface to the deep layer below the thermocline, a phenomenon referred to as the “inertial chimney” effect (Lee and Niiler 1998). This process is significant for the down-taking of wind-generated near-inertial energy to facilitate diapycnal mixing in the deep ocean (Wunsch and Ferrari 2004; Jing and Wu 2014). The inertial chimney effect of WEs has drawn a lot of attention and has been studied extensively using numerical models (Zhai et al. 2007; Whitt and Thomas 2013; Lelong et al. 2020) and observations (Alford and Whitmont 2007; Joyce et al. 2013; Martínez‐Marrero et al. 2019). The elevated turbulent dissipation driven by the trapped NIWs in the WEs was also observed in different regions across the global ocean (Lueck and Osborn 1986; Cuypers et al. 2012; Kawaguchi et al. 2016; Fer et al. 2018; Cyriac et al. 2021).
Previous studies have significantly elucidated the inertial chimney effect and the accompanying diapycnal mixing in the WEs. According to Young and Ben Jelloul’s (1997) theory, the presence of CEs could also enhance the downward propagation of NIWs. However, the modification of NIW propagation by the CEs has been rarely reported. In numerical models, downward near-inertial energy propagation appeared in the CEs but with low strength and shallow depth (Lee and Niiler 1998; Danioux et al. 2008). Jaimes and Shay (2010) observed that NIWs stalled in the upper layers of CEs and strengthened the vertical shears and turbulence. Kawaguchi et al. (2021) observed that weak NIWs generated by the local wind forcing were trapped in the upper 300 m of a CE. Based on the limited studies mentioned above, the CEs aid in accelerating the downward NIW propagation despite lower strength and shallower residence compared to the inertial chimney effect of the WEs, which eventually leads to elevated diapycnal mixing in the interior ocean. Furthermore, the down-taking of wind-generated near-inertial energy through the CEs has a significant impact on diapycnal mixing in the global ocean (Whalen et al. 2018). This necessitates additional efforts to elucidate the NIW propagation modified by the CEs.
The Kuroshio Extension (KE) in the northwestern Pacific is one of the most dynamically complex regions across the global ocean (Donohue et al. 2008). There are abundant energetic mesoscale eddies pinched off by the eastward jet and NIWs forced by the frequent storms (Qiu and Chen 2010; Rimac et al. 2013). In this study, we compared the NIWs modified by a CE and WE as well as the concurrent turbulent dissipation based on shipboard observations in the KE region, with particular focus on the downward near-inertial energy propagation in the CE. This study begins with an introduction section comprising field observations, data collection, and a brief description of the eddy characteristics in section 2. Section 3 presents a method for separating the balanced motions and waves from the shipboard acoustic Doppler current profiler (SADCP) measurement. Results of the NIW characteristics and downward near-inertial energy propagation modified by the eddies are shown in section 4, followed by an analysis of locations of turbulent dissipation in the eddies in section 5. The energy exchange between the NIWs and eddies, and the implications for the parameterizations of turbulent mixing are discussed in section 6. Section 7 presents the conclusions drawn from the study findings.
2. Data
a. Observations of two mesoscale eddies
A spring cruise in 2018 was conducted to maintain the Kuroshio Extension Mooring System (KEMS), constructed by Ocean University of China, in the KE region in the northwestern Pacific [see Zhu et al. (2021) for details about KEMS]. To explore the dynamics in the mesoscale eddies, we designed a comparative experiment in a CE and WE shed from the KE during 18–29 May (Fig. 1a). Locations of the two eddies were tracked by the near-real-time sea level anomaly (SLA) fields distributed by AVISO (https://www.aviso.altimetry.fr/). Finally, 2 (6) transects and 10 (50) temperature–salinity (T–S) stations were deployed in the CE (WE) (Figs. 1b,c).

(a) Ship tracks (black) and stations (red) during the observations on 18–29 May 2018, overlapped on the map of time-averaged SSH distributed by AVISO. (b) Measurements of velocity using ADCP and TKE dissipation rate ε and T–S by MSS in the warm-core eddy (WE). The background map is the time-averaged SLA during the observation of the eddy. The 50 stations (red dots) and their station numbers are labeled. The six transects T1–T6 are marked by colored segments in the subgraph. The velocity vectors at a depth of 100 m are denoted by gray arrows. (c) As in (b), but for the cold-core eddy (CE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

(a) Ship tracks (black) and stations (red) during the observations on 18–29 May 2018, overlapped on the map of time-averaged SSH distributed by AVISO. (b) Measurements of velocity using ADCP and TKE dissipation rate ε and T–S by MSS in the warm-core eddy (WE). The background map is the time-averaged SLA during the observation of the eddy. The 50 stations (red dots) and their station numbers are labeled. The six transects T1–T6 are marked by colored segments in the subgraph. The velocity vectors at a depth of 100 m are denoted by gray arrows. (c) As in (b), but for the cold-core eddy (CE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
(a) Ship tracks (black) and stations (red) during the observations on 18–29 May 2018, overlapped on the map of time-averaged SSH distributed by AVISO. (b) Measurements of velocity using ADCP and TKE dissipation rate ε and T–S by MSS in the warm-core eddy (WE). The background map is the time-averaged SLA during the observation of the eddy. The 50 stations (red dots) and their station numbers are labeled. The six transects T1–T6 are marked by colored segments in the subgraph. The velocity vectors at a depth of 100 m are denoted by gray arrows. (c) As in (b), but for the cold-core eddy (CE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
The research vessel (R/V DFH2) was equipped with a 75-kHz SADCP to record the underway velocity in the upper 650 m, with a typical ship speed of ∼10 kt (1 kt ≈ 0.51 m s−1). It stopped after every several tens of kilometers to conduct an 800-m cast of turbulent kinetic energy (TKE) dissipation rate ε and T–S using the Microstructure profiler (MSS90, Sea and Sun Technology GmbH). The single ping SADCP velocity was first processed by CODAS (https://currents.soest.hawaii.edu/docs/adcp_doc/index.html) and then averaged into 10-min ensembles. Velocity profiles acquired during the ship stopping at T–S stations (usually 1–2 h) and under slow navigation were discarded because of unacceptable errors for the SADCP measurement. The velocity profiles range from 20 to 650 m and comprise 80 bins, each with a bin size of 8 m. Velocity in the first bin was disregarded owing to numerous disqualified values flagged in CODAS. Finally, we obtained a set of velocity profiles with a spacing of approximately 3 km (navigating speed of ∼10 kt and average ensembles of 10 min). The ε and T–S profiles have an original vertical interval of 1 m. To match the vertical resolution of the SADCP velocity, they were averaged into 8-m bins.
b. Characteristics of the two eddies
As shown in Fig. 1, the two eddies are close to deformed ellipses with scales of ∼200 km in the horizontal direction. Their three-dimensional structures are visible from the velocities along transects in different directions (Fig. 2). The meridional squeezing and zonal stretching of the WE lead to a greater eastward velocity u than the northward velocity υ (Figs. 1b, 2b, and 2d), while the CE is meridionally stretched and zonally squeezed so that u is smaller than υ (Figs. 1c, 2a, and 2c). According to the vertical structure of the velocity, both eddies are surface intensified with the maximum velocity exceeding 1 m s−1, are attenuated with increasing depth, and remain energetic at the maximum observation depth of 650 m. Usually, mesoscale eddies in the KE region can penetrate a 1000-m depth (Dong et al. 2017). The eddy center can be clearly distinguished by the zero line of velocity. Unlike mesoscale eddies that have distinct vertical tilts of their axes in the marginal seas of complicated topography (Zhang et al. 2016), the two eddies’ axes show nearly no tilt in the open ocean. The up-heaving (down-heaving) isopycnals of the CE (WE) reach their maximum heaves around the eddy center, thereby showing that the T–S structure is highly consistent with the velocity field.

(a) Eastward velocity component u and (c) northward velocity component υ along the space–time varying ship tracks in the cold-core eddy (CE). (b),(d) As in (a) and (c), but for the warm-core eddy (WE). The transects and T–S stations in different directions corresponding to those in Fig. 1 are marked at the top of the horizontal axis. The T–S stations near the eddy center are highlighted with their station numbers labeled. Contours of the potential density anomaly σθ are plotted to aid in identifying the eddy structure.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

(a) Eastward velocity component u and (c) northward velocity component υ along the space–time varying ship tracks in the cold-core eddy (CE). (b),(d) As in (a) and (c), but for the warm-core eddy (WE). The transects and T–S stations in different directions corresponding to those in Fig. 1 are marked at the top of the horizontal axis. The T–S stations near the eddy center are highlighted with their station numbers labeled. Contours of the potential density anomaly σθ are plotted to aid in identifying the eddy structure.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
(a) Eastward velocity component u and (c) northward velocity component υ along the space–time varying ship tracks in the cold-core eddy (CE). (b),(d) As in (a) and (c), but for the warm-core eddy (WE). The transects and T–S stations in different directions corresponding to those in Fig. 1 are marked at the top of the horizontal axis. The T–S stations near the eddy center are highlighted with their station numbers labeled. Contours of the potential density anomaly σθ are plotted to aid in identifying the eddy structure.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
Furthermore, we calculated the stratification N2 [

As in Fig. 2, but for the (top) stratification N2 and (bottom) normalized relative vorticity ζ/f. The local inertial frequency f is 8.0 × 10−5 s−1 (8.9 × 10−5 s−1) for the CE (WE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

As in Fig. 2, but for the (top) stratification N2 and (bottom) normalized relative vorticity ζ/f. The local inertial frequency f is 8.0 × 10−5 s−1 (8.9 × 10−5 s−1) for the CE (WE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
As in Fig. 2, but for the (top) stratification N2 and (bottom) normalized relative vorticity ζ/f. The local inertial frequency f is 8.0 × 10−5 s−1 (8.9 × 10−5 s−1) for the CE (WE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
To compare the features of NIWs and turbulent dissipation inside and outside the eddy cores in the following sections, we estimated the radius of each eddy to distinguish the two regions. The location of the depth-averaged eddy center was first determined by the local minimum of the gridded velocity. As shown in Figs. 4a and 4b, the estimated eddy center of the CE (WE) is located at 33.14°N, 152.00°E (37.70°N, 147.00°E), close to that determined by the local minimum (maximum) of SLA as shown in Fig. 1. Usually, the eddy velocity increases away from the eddy center, reaches its maximum at the edge of the eddy core, and then decreases. The linear increase along with the exponential decay and the Gaussian profiles were commonly applied to the radial distribution of the velocity in mesoscale eddies (Olson 1980). Hence, we made a radial composition of the velocity with a 10-km average in the radial direction. The velocity increases first and then decreases as the radius increases, which is quite remarkable in the upper 200 m. The depth-averaged velocity above 200 m is therefore used to show the radial eddy structure, which displays an apparent Gaussian distribution for both eddies (Figs. 4c,d). Fitting the depth-averaged velocity with the Gaussian function,

(a) Map of the depth-averaged (above 200 m) velocity on the 3 km × 3 km horizontal grids by linear interpolation for the cold-core eddy (CE). The depth-averaged eddy center estimated by the local minimum velocity is marked by the red dot. The regions of missing values are masked by gray dots. The ship tracks are overlaid by the gray dashed lines. (c) The depth-averaged velocity above 200 m (black dots) is fitted with the Gaussian profile (red line) to estimate the eddy radius r0 of the CE, using the location of the maximum velocity V0. Outliers (gray dots) at the eddy periphery are excluded from the fitting. (b),(d) As in (a) and (c), but for the warm-core eddy (WE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

(a) Map of the depth-averaged (above 200 m) velocity on the 3 km × 3 km horizontal grids by linear interpolation for the cold-core eddy (CE). The depth-averaged eddy center estimated by the local minimum velocity is marked by the red dot. The regions of missing values are masked by gray dots. The ship tracks are overlaid by the gray dashed lines. (c) The depth-averaged velocity above 200 m (black dots) is fitted with the Gaussian profile (red line) to estimate the eddy radius r0 of the CE, using the location of the maximum velocity V0. Outliers (gray dots) at the eddy periphery are excluded from the fitting. (b),(d) As in (a) and (c), but for the warm-core eddy (WE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
(a) Map of the depth-averaged (above 200 m) velocity on the 3 km × 3 km horizontal grids by linear interpolation for the cold-core eddy (CE). The depth-averaged eddy center estimated by the local minimum velocity is marked by the red dot. The regions of missing values are masked by gray dots. The ship tracks are overlaid by the gray dashed lines. (c) The depth-averaged velocity above 200 m (black dots) is fitted with the Gaussian profile (red line) to estimate the eddy radius r0 of the CE, using the location of the maximum velocity V0. Outliers (gray dots) at the eddy periphery are excluded from the fitting. (b),(d) As in (a) and (c), but for the warm-core eddy (WE).
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
3. Separating balanced motions and waves from shipboard ADCP measurement
a. Method
The SADCP velocity contains multifrequency signals contributed by different motions, including the geostrophic and cyclogeostrophic balanced flows at low frequencies and the unbalanced part dominated by internal gravity waves at high frequencies (Torres et al. 2018). Conventionally, a temporal or spatial filter is used to separate motions of different scales (Thomson and Emery 2014). However, the space–time aliasing of the SADCP velocity precludes the separation of the balanced motions and waves by the frequency-domain filter (Alford et al. 2013). Considering the mismatch between the vertical wavelengths of the balanced motions and waves, applying the vertical wavenumber filter to a single velocity profile to realize the separation is feasible. Wavenumber-domain filter (Karstensen et al. 2017; Martínez‐Marrero et al. 2019) and fit (Kunze 1986; Joyce et al. 2013) were commonly used for a velocity profile to separate the mesoscale balanced flow and NIWs by first obtaining the trend of the balanced flow and then subtracting the trend from the total velocity to derive the NIWs. However, the selection of vertical cutoff wavenumber and fitting functions mainly relies on subjective judgments, requiring sensitivity experiments to identify optimal practice.

Schematic of extracting the balanced velocity (blue) from the SADCP velocity profiles (gray) by MRM. Positions of the T–S station pair and velocity profiles are denoted by red and black dots, respectively. The
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

Schematic of extracting the balanced velocity (blue) from the SADCP velocity profiles (gray) by MRM. Positions of the T–S station pair and velocity profiles are denoted by red and black dots, respectively. The
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
Schematic of extracting the balanced velocity (blue) from the SADCP velocity profiles (gray) by MRM. Positions of the T–S station pair and velocity profiles are denoted by red and black dots, respectively. The
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
b. Validation of mode regression method

(a) T–S stations (black dots) and the station pairs selected for the calculation of geostrophic velocity (red dots) in the cold-core eddy (CE). (b) Comparison of the SADCP (gray), MRM-derived (blue), and geostrophic/cyclogeostrophic balanced (red) velocity profiles in the CE. GB and CGB are the abbreviations of geostrophic balance and cyclogeostrophic balance, respectively. Root-mean-square (RMS) bias between the diagnosed balanced velocity and the MRM-derived balanced velocity is marked in black text. (c),(d) As in (a) and (b), but for the warm-core eddy (WE). The cyclogeostrophic (geostrophic) balanced velocity of the CE (WE) is regarded as the reference for examining the MRM-derived balanced velocity.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

(a) T–S stations (black dots) and the station pairs selected for the calculation of geostrophic velocity (red dots) in the cold-core eddy (CE). (b) Comparison of the SADCP (gray), MRM-derived (blue), and geostrophic/cyclogeostrophic balanced (red) velocity profiles in the CE. GB and CGB are the abbreviations of geostrophic balance and cyclogeostrophic balance, respectively. Root-mean-square (RMS) bias between the diagnosed balanced velocity and the MRM-derived balanced velocity is marked in black text. (c),(d) As in (a) and (b), but for the warm-core eddy (WE). The cyclogeostrophic (geostrophic) balanced velocity of the CE (WE) is regarded as the reference for examining the MRM-derived balanced velocity.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
(a) T–S stations (black dots) and the station pairs selected for the calculation of geostrophic velocity (red dots) in the cold-core eddy (CE). (b) Comparison of the SADCP (gray), MRM-derived (blue), and geostrophic/cyclogeostrophic balanced (red) velocity profiles in the CE. GB and CGB are the abbreviations of geostrophic balance and cyclogeostrophic balance, respectively. Root-mean-square (RMS) bias between the diagnosed balanced velocity and the MRM-derived balanced velocity is marked in black text. (c),(d) As in (a) and (b), but for the warm-core eddy (WE). The cyclogeostrophic (geostrophic) balanced velocity of the CE (WE) is regarded as the reference for examining the MRM-derived balanced velocity.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
As shown in Fig. 6, although the diagnosed balanced velocity shows a trend that is similar to that of the SADCP velocity, the two velocities exhibit large shifts at certain station pairs, which probably results from the fast-varying flow compared to the slow-varying T–S owing to the eddy deformation and/or movement, like station pairs 1–2 and 19–20 (1–2 and 26–27) at the eddy periphery of the CE (WE), and station pairs 3–4 and 16–17 (21–22 and 22–23) at the eddy center of the CE (WE). At other station pairs, the diagnosed balanced velocity is consistent with the SADCP velocity trend and can therefore be regarded as a reference for examining the accuracy of the balanced velocity derived by MRM. The root-mean-square bias between the diagnosed balanced velocity and the MRM-derived balanced velocity does not exceed 5 (6) cm s−1 in the CE (WE), significantly lower than the balanced velocity and also lower than the wave amplitude (i.e., residual velocity), accounting for approximately one-third (one-fifth) of the typical wave amplitude (shown in the subsequent section). Therefore, MRM was confirmed to be valid for separating the balanced motions and waves from the SADCP velocity.
MRM is a fit that is similar to the smooth spline fit used by Joyce et al. (2013). We found that the residual velocities derived by MRM and the smooth spline fit differed slightly (figures not shown). The smooth spline fit was performed using the curve fitting toolbox in MATLAB, with the smoothing parameters artificially selected as 0.1 to better fit the trend of the SADCP velocity. In addition, the 350-m lowpass Butterworth filter used by Martínez‐Marrero et al. (2019) was also applied to the SADCP velocity to separate the balanced motions and waves, but it failed to completely extract the waves owing to the relatively long vertical wavelengths we observed. Therefore, MRM has the advantage of objectivity compared to the vertical wavenumber filter and fit. Theoretically, MRM is applicable for regions dominated by geostrophic/cyclogeostrophic dynamics, such as in the mesoscale eddies as well as western boundary currents and their jet-like extensions.
4. Enhanced NIWs observed in the two eddies
The waves separated from the SADCP velocity are mainly composed of internal waves, like NIWs and internal tides (Torres et al. 2018). Internal tides are weak in the upper 1000 m in this region (Zhao et al. 2016). Therefore, the waves are speculated to be dominated by NIWs. We estimated the vertical wavenumber spectra of the velocity and shear of the waves, and compared them with the GM spectra. The WKB normalization was applied to the waves to remove the effect of stratification. As shown in Fig. 7, the observed spectra are consistent with the GM spectra, and the spectral peak is located at a vertical wavelength ranging between 200 and 400 m, indicative of the dominance of NIWs.

Clockwise (CW) and counterclockwise (CCW) vertical wavenumber spectra of the (a) NIW velocity and (b) NIW shear in the cold-core eddy (CE). (c),(d) As in (a) and (b), but for the warm-core eddy (WE). The displayed spectra are the average spectra inside (r < r0/2, “In”) and outside (r > r0, “Out”) the eddy core, and r0 is the radius of the eddy core. The 95% confidence interval is marked by the error bar. The vertical wavelengths are shown at the top of the horizontal axis. The NIW shear is estimated using the 32-m vertical gradient of the NIW velocity smoothed over a 32-m interval to remove the noise at high wavenumbers. The GM spectra are calculated using the MATLAB toolbox downloaded from http://jklymak.github.io/GarrettMunkMatlab/.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

Clockwise (CW) and counterclockwise (CCW) vertical wavenumber spectra of the (a) NIW velocity and (b) NIW shear in the cold-core eddy (CE). (c),(d) As in (a) and (b), but for the warm-core eddy (WE). The displayed spectra are the average spectra inside (r < r0/2, “In”) and outside (r > r0, “Out”) the eddy core, and r0 is the radius of the eddy core. The 95% confidence interval is marked by the error bar. The vertical wavelengths are shown at the top of the horizontal axis. The NIW shear is estimated using the 32-m vertical gradient of the NIW velocity smoothed over a 32-m interval to remove the noise at high wavenumbers. The GM spectra are calculated using the MATLAB toolbox downloaded from http://jklymak.github.io/GarrettMunkMatlab/.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
Clockwise (CW) and counterclockwise (CCW) vertical wavenumber spectra of the (a) NIW velocity and (b) NIW shear in the cold-core eddy (CE). (c),(d) As in (a) and (b), but for the warm-core eddy (WE). The displayed spectra are the average spectra inside (r < r0/2, “In”) and outside (r > r0, “Out”) the eddy core, and r0 is the radius of the eddy core. The 95% confidence interval is marked by the error bar. The vertical wavelengths are shown at the top of the horizontal axis. The NIW shear is estimated using the 32-m vertical gradient of the NIW velocity smoothed over a 32-m interval to remove the noise at high wavenumbers. The GM spectra are calculated using the MATLAB toolbox downloaded from http://jklymak.github.io/GarrettMunkMatlab/.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
a. Eddy-modified NIW characteristics
Figure 8 shows the NIW structure along the transects in different directions. NIWs are enhanced inside the eddy cores of both eddies. The typical NIW amplitude in the CE is approximately 0.15 m s−1, half of that of the WE, which is 0.3 m s−1. In the CE, the strong NIWs are mainly concentrated around the eddy center at a depth exceeding 350 m. On the contrary, the NIWs are amplified at the base of the eddy core at depths below 350 m in the WE, which is the feature of the trapped NIWs encountering the critical layers (Kunze 1986; Martínez‐Marrero et al. 2019). The NIW structure in the CE has been rarely observed and compared with that in the WE. Kawaguchi et al. (2021) observed attenuated NIWs across the inner core of a CE, but the most energetic waves were packed at a depth above 150 m, probably because the CE was weaker (<0.3 m s−1) and shallower (above 200-m depth) compared to the CE shed from the KE in this study.

As in Fig. 2, but for the NIW velocity; u′ and υ ′ are the eastward and northward velocity components.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

As in Fig. 2, but for the NIW velocity; u′ and υ ′ are the eastward and northward velocity components.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
As in Fig. 2, but for the NIW velocity; u′ and υ ′ are the eastward and northward velocity components.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
The NIWs have constant phase lines almost slanting along the heaving isopycnals, which is consistent with the NIW characteristics shown in previous studies (Alford 2010; Joyce et al. 2013). The heaving NIW phase along with the isopycnals was found to be the manifestation of water parcels oscillating along the isopycnals under the strong baroclinicity of mesoscale eddies (Whitt and Thomas 2013). NIW phase θ can be identified by combining the signs of the velocity components, u′ and υ ′, of the NIW velocity (i.e., the direction of the velocity) in Fig. 8 or calculated by the equation,

As in Fig. 2, but for NIW phase θ (top) before and (bottom) after inertial back-rotation. The θ is calculated with the NIW velocity components horizontally smoothed over 30 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

As in Fig. 2, but for NIW phase θ (top) before and (bottom) after inertial back-rotation. The θ is calculated with the NIW velocity components horizontally smoothed over 30 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
As in Fig. 2, but for NIW phase θ (top) before and (bottom) after inertial back-rotation. The θ is calculated with the NIW velocity components horizontally smoothed over 30 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
After inertial back-rotation, most of the discontinued NIW phase along the ship tracks resulting from the stops at T–S stations and slow navigation is removed (Figs. 8–10; Shcherbina et al. 2003). Upon examining the phase change in transect 3 of the WE, discontinuities are observed to vanish. This transect has opposing NIW phases on the two sides owing to the intrinsic long and slow navigation that is sufficient for the NIWs to rotate half a circle (Figs. 8b, 8d, and 9b), whereas the phases become consistent across the transect after inertial back-rotation (Figs. 9d, 10b, and 10d). Generally, the NIW phases in both eddies become more coherent along the isopycnals across the entire eddy after inertial back-rotation (Fig. 9). For example, most of the NIW phase in 25.5–26.0 kg m−3 σθ (26.0–26.1 kg m−3 σθ) in the CE (WE) is positive along the transects in different directions, indicating that the observed NIWs in each eddy are all part of the same wave packet.

As in Fig. 2, but for the NIW velocity after inertial back-rotation; u′ and υ ′ are the eastward and northward velocity components.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

As in Fig. 2, but for the NIW velocity after inertial back-rotation; u′ and υ ′ are the eastward and northward velocity components.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
As in Fig. 2, but for the NIW velocity after inertial back-rotation; u′ and υ ′ are the eastward and northward velocity components.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
According to the vertical and horizontal variations of the NIW phase after inertial back-rotation (Figs. 9c,d), the wavelengths can be determined by the NIW phase spanning a complete cycle (changing from −180° to 180°). The vertical wavelengths of the NIWs are approximately 300 m in the CE and 400 m in the WE. In the horizontal direction, the NIW phase is almost consistent throughout the eddy for both eddies; therefore, the horizontal NIW scales are comparable to the eddy size.
b. Downward propagation of NIW energy
The vertical NIW propagation is visible from the changing of the NIW phase versus depth (Fig. 9). The NIWs are of clockwise rotation by the dominance of clockwise energy (recall Fig. 7) and a 90° phase lead for u′ versus υ′ through a downward view in Figs. 8 and 10. Generally, clockwise rotating NIWs (in the Northern Hemisphere) have a downward group speed (Leaman and Sanford 1975). Such energetic NIWs of downward propagation are usually generated by the wind. Because NIW energy is dominated by the near-inertial kinetic energy (NIKE), we calculated the NIKE with the NIW velocity that was first processed by WKB normalization to remove the stratification effect and then smoothed over a 32-m interval to eliminate noise at high vertical wavenumbers. In addition, we made a radial composition to the NIKE to show its distribution in the eddy. As shown in Figs. 11 and 12, the enhanced NIKE more than the canonical value of 3 J m−3 (Munk 1981) appears inside the inner core (r < r0/2) of the WE and CE. The enhancement of NIKE inside the inner core is thought to result from the inward and downward energy propagation (Fig. 9). For this reason, we considered the average of the NIW profiles inside the inner cores when comparing the NIW spectra inside and outside the eddy cores in Fig. 7. Then, we would separately analyze the three-dimensional structures of the NIKE in the two eddies.

Three-dimensional structure of the near-inertial kinetic energy (NIKE) in the warm-core eddy (WE). (a) Radial section of the NIKE averaged every 10 km in the radial direction. The contours of the potential density anomaly σθ are overlaid with dashed lines. (b) Radial distribution of the depth-averaged NIKE. Average profiles of the NIKE (full line) and relative vorticity (dashed line) inside (r < r0 /2) and outside (r > r0) the eddy core in the (c) σθ and (d) depth coordinate. (e) Clockwise (CW) and counterclockwise (CCW) NIKE inside the eddy core. The canonical value of 3 J m−3 (gray dashed line) in (b)–(e) is marked for reference. The NIKE is calculated by NIW velocity that is first WKB normalized and vertically smoothed over a 32-m interval. r is the distance from the eddy center, and r0 is the eddy radius of 56 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

Three-dimensional structure of the near-inertial kinetic energy (NIKE) in the warm-core eddy (WE). (a) Radial section of the NIKE averaged every 10 km in the radial direction. The contours of the potential density anomaly σθ are overlaid with dashed lines. (b) Radial distribution of the depth-averaged NIKE. Average profiles of the NIKE (full line) and relative vorticity (dashed line) inside (r < r0 /2) and outside (r > r0) the eddy core in the (c) σθ and (d) depth coordinate. (e) Clockwise (CW) and counterclockwise (CCW) NIKE inside the eddy core. The canonical value of 3 J m−3 (gray dashed line) in (b)–(e) is marked for reference. The NIKE is calculated by NIW velocity that is first WKB normalized and vertically smoothed over a 32-m interval. r is the distance from the eddy center, and r0 is the eddy radius of 56 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
Three-dimensional structure of the near-inertial kinetic energy (NIKE) in the warm-core eddy (WE). (a) Radial section of the NIKE averaged every 10 km in the radial direction. The contours of the potential density anomaly σθ are overlaid with dashed lines. (b) Radial distribution of the depth-averaged NIKE. Average profiles of the NIKE (full line) and relative vorticity (dashed line) inside (r < r0 /2) and outside (r > r0) the eddy core in the (c) σθ and (d) depth coordinate. (e) Clockwise (CW) and counterclockwise (CCW) NIKE inside the eddy core. The canonical value of 3 J m−3 (gray dashed line) in (b)–(e) is marked for reference. The NIKE is calculated by NIW velocity that is first WKB normalized and vertically smoothed over a 32-m interval. r is the distance from the eddy center, and r0 is the eddy radius of 56 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

As in Fig. 11, but for the three-dimensional structure of the NIKE in the CE. The eddy radius r0 is 77 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

As in Fig. 11, but for the three-dimensional structure of the NIKE in the CE. The eddy radius r0 is 77 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
As in Fig. 11, but for the three-dimensional structure of the NIKE in the CE. The eddy radius r0 is 77 km.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
In the WE, the enhanced NIKE appears in regions of negative vorticity, which is known as the inertial chimney effect (Lee and Niiler 1998; Fig. 11). Generally, the NIKE is amplified at the base of the eddy core at depths ranging from 350 to 650 m with the maxima exceeding 40 J m−3 (Fig. 11a), which results from the stall of the NIWs that have their vertical wavelengths shrunk and amplitudes enlarged at the critical layers. Radially, the depth-averaged NIKE reaches 19 J m−3 at the eddy center and decreases away from the center (Fig. 11b). This radial structure of the NIKE is consistent with the azimuthal-mode-one eddy-trapped NIWs (Kunze and Boss 1998). Vertically, approximately 67% of the NIKE propagates downward via the inertial chimney effect. The average NIKE inside the inner core peaks at ∼18 J m−3 and a depth of ∼450 m, which is approximately 6 times larger than that at the same depth outside the eddy core (Figs. 11c–e).
Similar to the WE, the enhanced NIKE appears in regions of positive vorticity inside the inner core of the CE, but with lower strength and shallower residence (Fig. 12). There is an energetic wave packet between 150 and 350 m with the maxima exceeding 10 J m−3 (Fig. 12a), only one-quarter of that in the WE. The depth-averaged NIKE shows a decreasing radial structure, similar to that in the WE (Fig. 12b). The average NIKE inside the inner core peaks at ∼6 J m−3 at a depth of ∼250 m, which is approximately 3 times larger than that at the same depth outside the eddy core (Fig. 12d). Jaimes and Shay (2010) have observed the vertical structure of NIW amplitude in a CE that is similar to Fig. 12d, which is intensified in the upper 250 m and reaches the maximum at ∼125-m depth, shallower than the 250-m depth in our observation. The difference is the NIWs are dominated by the upward propagation component in their study (Jaimes and Shay 2010).
The enhanced NIKE in the inner core of the CE is thought to result from the vertical dispersion of the NIWs accelerated by the existence of relative vorticity in mesoscale eddy fields (Young and Ben Jelloul 1997). Unlike the negative vorticity in the WE, which can trap the NIWs and drain the NIW energy deep into the eddy base through the inertial chimney effect, the positive vorticity seems to have less influence on the downward NIW propagation so that the enhancement of the NIKE is relatively weak and shallow (Figs. 12d,e).
5. Enhanced turbulent dissipation observed in the two eddies
As shown in Figs. 13c and 13f, the TKE dissipation rate ε is as high as 10−7 W kg−1 in both eddies, accompanied by an elevated diapycnal mixing with the vertical diffusivity Kρ (Kρ = 0.2ε/N2, Osborn 1980) up to 5 × 10−4 m2 s−1, which is 50 times larger than the average value in the upper 1000 m (10−5 m2 s−1, Waterhouse et al. 2014). To explain the mechanisms of the turbulence generation, the vertical shears of the NIWs and mesoscale balanced flow were calculated (Figs. 13a,b,d, e). The NIW shear is one order of magnitude larger than the mesoscale shear and is mainly responsible for turbulence production. The gradient Richardson number, Ri, was estimated as follows: Ri = N2/S2. Similar to the estimation of N2, the total shear S2 was calculated using the 32-m vertical gradient of the total velocity that was first smoothed over a 32-m interval, and the values below 10−6 s−1 were replaced by 10−6 to remove noise at high wavenumbers (Fer et al. 2018). As shown in Figs. 14c and 14h, the average Ri significantly exceeds its critical value of 0.25 for the production of turbulence because of the 32-m average and the difficulty in detecting the patchy and intermittent turbulence by a velocity of 8-m vertical space and 10-min average. The average Ri is below 15 with the increasing ε, which might suggest the production of turbulence by shear instability.

Radial sections of the (a) NIW shear
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

Radial sections of the (a) NIW shear
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
Radial sections of the (a) NIW shear
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

Vertical structures of the (a) stratification N2, (b) total shear S2, (c) Richardson number Ri, (d) TKE dissipation rate ε, and (e) vertical diffusivity Kρ inside (red) and outside (blue) the core of the CE. (f)–(j) As in (a)–(e), but for the WE. The N2 and S2 are estimated by the 32-m vertical gradient of the density and total velocity that are smoothed over a 32-m interval. The S2 are calculated with velocity profiles at the T–S stations to match the observation of the TKE dissipation rate ε.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

Vertical structures of the (a) stratification N2, (b) total shear S2, (c) Richardson number Ri, (d) TKE dissipation rate ε, and (e) vertical diffusivity Kρ inside (red) and outside (blue) the core of the CE. (f)–(j) As in (a)–(e), but for the WE. The N2 and S2 are estimated by the 32-m vertical gradient of the density and total velocity that are smoothed over a 32-m interval. The S2 are calculated with velocity profiles at the T–S stations to match the observation of the TKE dissipation rate ε.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
Vertical structures of the (a) stratification N2, (b) total shear S2, (c) Richardson number Ri, (d) TKE dissipation rate ε, and (e) vertical diffusivity Kρ inside (red) and outside (blue) the core of the CE. (f)–(j) As in (a)–(e), but for the WE. The N2 and S2 are estimated by the 32-m vertical gradient of the density and total velocity that are smoothed over a 32-m interval. The S2 are calculated with velocity profiles at the T–S stations to match the observation of the TKE dissipation rate ε.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
In the WE (Figs. 13d,f), the elevated turbulent dissipation shows good consistency with the strong NIW shear, implying the efficient turbulence production driven by the NIWs. The turbulent dissipation is enhanced in two regions. One is in the surface layer of high stratification, which is in funnel shape following the down-heaving isopycnals that extend to depths of 200 and 100 m in the eddy center and periphery, respectively. Another is at the base of the eddy core at a depth ranging from 350 to 650 m, which results from the shrinking of the vertical NIW wavelengths at the critical layers. Between the two regions, a reduced turbulent dissipation following the suppression of the NIW shear is a common feature of these structures because of the increase in vertical wavelength caused by reduced stratification and negative vorticity (Fernández-Castro et al. 2020). In other words, the NIWs generated in the surface layer are radially trapped in the WE, propagate downward through the region of negative vorticity inside the eddy core, and finally break to drive the elevated turbulent dissipation and vertical diffusivity inside the eddy core, which are 3–4 times larger than those on the outside (Figs. 14g–j). The high–low–high sandwiched structure of the turbulent dissipation associated with trapped NIWs in the WE is generally consistent with observations made in other oceans (e.g., Lueck and Osborn 1986; Cuypers et al. 2012; Fer et al. 2018).
Turbulent dissipation following the NIW shear is also expected in the CE. Figures 13a and 13c show that the elevated turbulent dissipation mainly occurs in the sloping region above the 26.0 kg m−3 σθ, which is similar to the enhanced thermocline mixing due to eddy-induced geostrophic shear in a permanent CE, the Mindanao Eddy (Liu et al. 2017). The eddy shear is smaller than the NIW shear in our observation (Figs. 13a,b), and thus is unlikely to induce such strong turbulent dissipation. Therefore, we attempt to partly explain the turbulence production by the shear instability of NIWs. Locally, the turbulent dissipation enhanced in the high-stratified surface layer above a depth of 100 m is consistent with the strong NIW shear there, and that sporadically appears in regions below the surface layer inside the eddy core might result from the strong NIW shear as well. However, the turbulent dissipation is significantly higher and more concentrated in the eddy periphery (25.0–26.0 kg m−3 σθ, deep to 400-m depth), where the NIW shear is weaker than that inside the eddy core but the eddy kinetic energy is higher. We speculate there are other mechanisms of producing turbulence by cascading the eddy energy forward. For example, the submesoscale motions commonly generated in the eddy periphery (Yang et al. 2017) can penetrate far beneath the mixed layer through the presence of a weakly stratified transitional layer (e.g., the layer in 25.0–25.5 kg m−3 σθ in the CE in Figs. 3a) between the mixed layer and the thermocline, and eventually transfer the eddy energy to the turbulence through symmetric instability (Zhang et al. 2021).
In this study, we focus on the turbulent dissipation and diapycnal mixing driven by NIWs. Therefore, the enhancement of turbulent dissipation and diapycnal mixing inside the core of CE is relative to the global-averaged values in the upper ocean, rather than those outside the CE. In fact, the average turbulent dissipation and diapycnal mixing inside the CE are comparable with those driven by the trapped NIWs inside the WE (Figs. 14d,e,i,j). This indicates enhanced turbulent dissipation and diapycnal mixing driven by the downward propagating NIWs inside the CE.
6. Discussion
a. NIWs-eddy energy exchange
As shown in Figs. 15a and 15c, the NIWs both extract energy from (positive Fh) and transfer energy to (negative Fh) the mesoscale balanced flow, indicating the lack of a unique direction for the NIWs–balanced flow energy exchange (Thomas and Daniel 2020). The instantaneous energy transform rate is as high as 10−7 W kg−1, the same as that calculated in numerical model for the NIWs trapped by the WE (Lelong et al. 2020), but is two to three orders larger than the time-averaged values (Jing et al. 2018; Cusack et al. 2020).

(a) Horizontal energy transform rate Fh between NIWs and mesoscale eddies and (b) the average Fh inside the core of the CE vs the balance-to-NIWs energy. (c),(d) As in (a) and (b), but for the WE. The positive and negative Fh indicate the increase and decrease of the NIW energy, respectively.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1

(a) Horizontal energy transform rate Fh between NIWs and mesoscale eddies and (b) the average Fh inside the core of the CE vs the balance-to-NIWs energy. (c),(d) As in (a) and (b), but for the WE. The positive and negative Fh indicate the increase and decrease of the NIW energy, respectively.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
(a) Horizontal energy transform rate Fh between NIWs and mesoscale eddies and (b) the average Fh inside the core of the CE vs the balance-to-NIWs energy. (c),(d) As in (a) and (b), but for the WE. The positive and negative Fh indicate the increase and decrease of the NIW energy, respectively.
Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0160.1
By comparing the average energy transform rate and the balance-to-NIWs energy inside the eddy, there appears to be a positive relation between them (Figs. 15b,d). In addition, there appears to be a critical value for the balance-to-NIWs energy to separate the positive and negative energy transform rate. When the balance-to-NIWs energy is high, the NIWs mainly draw energy from the balanced flow; however, when this ratio drops below the critical value, the NIWs transform energy to the balanced flow. The finding likely supports the numerical result that the NIWs transfer energy to balanced flows in regimes where balance-to-wave energy is low (Thomas and Daniel 2020).
b. Implication for parameterizing the wind-driven diapycnal mixing
Theoretical and numerical studies showed downward NIW energy propagation enhanced by the presence of mesoscale eddies (Kunze 1985; Young and Ben Jelloul (1997); Zhai et al. 2005; Danioux et al. 2008). In this study, the enhanced downward NIW energy propagation was observed inside the eddy cores of a CE and WE, which led to an elevated turbulent dissipation and diapycnal mixing in the interior ocean. The down-taking of wind-generated near-inertial energy through the inertial chimney effect of the WEs plays an important role in furnishing turbulent diapycnal mixing in the deep ocean (Kunze et al. 1995; Jing and Wu 2014; Fer et al. 2018). However, the downward NIW energy propagation enhanced by the CEs is somewhat debatable and lacks observation evidence. In our observations, the NIKE was enhanced inside the eddy core of the CE and facilitated an elevated turbulent dissipation below the surface layer, albeit with lower strength and shallower residence than those in the WE. From a global perspective, the wind-generated near-inertial energy propagating downward through the CEs might be nonnegligible in the eddying ocean, which has significant implications for the parameterization of wind-driven diapycnal mixing in climate models (Wunsch and Ferrari 2004; Whalen et al. 2020).
7. Conclusions
Shipboard observations of the upper-ocean current, temperature–salinity, and turbulent dissipation rate were conducted in a cold-core eddy (CE) and warm-core eddy (WE) in the Kuroshio Extension (KE) region in May 2018. The maximum velocities of the two eddies shed by the KE exceeded 1 m s−1 and the relative vorticity magnitude was as high as 0.6f. They were close to deformed ellipses with a scale of ∼200 km in the horizontal direction, and eddy core radii of 77 and 56 km for the CE and WE, respectively. Both eddies were surface intensified and remained energetic at the maximum observation depth of 650 m. The mode regression method (MRM) was proposed to separate the balanced motions and waves from the space–time aliasing shipboard ADCP velocities, and was confirmed to be valid and objective. In this study, the waves extracted by MRM were dominated by near-inertial waves (NIWs).
The enhanced NIWs were found inside the eddy cores of both eddies, with typical NIW amplitudes reaching up to 0.15 and 0.3 m s−1 in the CE and WE, respectively. The NIWs had their constant phase lines slanting along the heaving isopycnals under the modification of the strong baroclinic eddy flow, with the vertical wavelengths being approximately 300–400 m and the horizontal scale being comparable to the eddy size. They were dominated by clockwise rotation and efficiently propagated downward in regions of strong relative vorticity inside the eddy cores, leading to enhancement of the near-inertial kinetic energy (NIKE) far beneath the surface layer. In the WE, the NIKE was enhanced at the base of the eddy core at depths ranging between 350 and 650 m owing to the amplification of the downward propagating NIWs at the critical layers. The maximum NIKE was as high as 40 J m−3 and its average value inside the eddy core reached up to ∼18 J m−3 at a depth of ∼450 m, which is approximately 6 times larger than that at the same depth outside the eddy. In the CE, the enhanced NIKE was located at depths ranging between 150 and 350 inside the inner core of the eddy, with the maxima exceeding 10 J m−3 and an average of ∼6 J m−3 at a depth of ∼250 m, which is approximately 3 times larger than that at the same depth outside the eddy core. Compared to the deep propagation of the NIKE through the inertial chimney effect in the WE, the enhanced NIKE with lower strength and shallower residence in the CE was speculated to result from the enhanced downward dispersion of the NIWs by the existence of strong relative vorticity.
The downward propagating NIWs had strong shears and drove an elevated turbulent dissipation rate ε (10−7 W kg−1) and a vertical diffusivity Kρ (5 × 10−4 m2 s−1), below the surface layer in both eddies. In the WE, the elevated turbulent dissipation followed the strong NIW shear, which was enhanced in the high-stratified surface layer above the upper 200 m, weak in the low-stratified midlayer, and strongest at the base of the eddy core of the WE because of the variation of NIW vertical wavelengths with the changing stratification and relative vorticity versus depth. In the CE, the strong NIW shear also induced an enhanced turbulent dissipation that was concentrated in the surface layer above 100 m and patchily appeared in regions of strong NIW shear inside the eddy core.
The results in the present study show clear images of the eddy-modified NIW structures and turbulent dissipation in a CE and WE in the KE region and provide observational evidence for the downward NIW propagation accelerated by the presence of mesoscale eddies, which has significant implications for parameterizing the wind-driven diapycnal mixing in the eddying ocean. In the future, it is necessary to estimate the global diapycnal mixing based on high-resolution ocean circulation and climate models that take the downward NIW propagation modified by the CEs into account.
Acknowledgments.
We thank the three reviewers for their detailed and constructive comments that helped us greatly improve our manuscript. We also thank all crew of the R/V DFH2 for their great efforts. This study was funded by the National Natural Science Foundation of China (42076009, 41822601, 41806008 and 41776006), Qingdao Pilot National Laboratory for Marine Science and Technology (2017ASTCP-ES05), and Fundamental Research Funds for the Central Universities (201762013 and 202072001). Z. C. is partially supported by Taishan Scholar Funds (tsqn201812022).
REFERENCES
Alford, M. H., 2003: Redistribution of energy available for ocean mixing by long-range propagation of internal waves. Nature, 423, 159–162, https://doi.org/10.1038/nature01628.
Alford, M. H., 2010: Sustained, full-water-column observations of internal waves and mixing near Mendocino Escarpment. J. Phys. Oceanogr., 40, 2643–2660, https://doi.org/10.1175/2010JPO4502.1.
Alford, M. H., and M. Whitmont, 2007: Seasonal and spatial variability of near-inertial kinetic energy from historical moored velocity records. J. Phys. Oceanogr., 37, 2022–2037, https://doi.org/10.1175/JPO3106.1.
Alford, M. H., A. Y. Shcherbina, and M. C. Gregg, 2013: Observations of near-inertial internal gravity waves radiating from a frontal jet. J. Phys. Oceanogr., 43, 1225–1239, https://doi.org/10.1175/JPO-D-12-0146.1.
Alford, M. H., J. A. MacKinnon, H. L. Simmons, and J. D. Nash, 2016: Near-inertial internal gravity waves in the ocean. Annu. Rev. Mar. Sci., 8, 95–123, https://doi.org/10.1146/annurev-marine-010814-015746.
Bosse, A., and Coauthors, 2017: A submesoscale coherent vortex in the Ligurian Sea: From dynamical barriers to biological implications. J. Geophys. Res. Oceans, 122, 6196–6217, https://doi.org/10.1002/2016JC012634.
Bühler, O., and M. E. McIntyre, 2005: Wave capture and wave–vortex duality. J. Fluid Mech., 534, 67–95, https://doi.org/10.1017/S0022112005004374.
Chelton, D. B., M. G. Schlax, R. M. Samelson, and R. A. de Szoeke, 2007: Global observations of large oceanic eddies. Geophys. Res. Lett., 34, L15606, https://doi.org/10.1029/2007GL030812.
Chelton, D. B., M. G. Schlax, and R. M. Samelson, 2011: Global observations of nonlinear mesoscale eddies. Prog. Oceanogr., 91, 167–216, https://doi.org/10.1016/j.pocean.2011.01.002.
Chen, G., H. Xue, D. Wang, and Q. Xie, 2013: Observed near-inertial kinetic energy in the northwestern South China Sea. J. Geophys. Res. Oceans, 118, 4965–4977, https://doi.org/10.1002/jgrc.20371.
Cusack, J. M., J. A. Brearley, A. C. N. Garabato, D. A. Smeed, K. L. Polzin, N. Velzeboer, and C. J. Shakespeare, 2020: Observed eddy–internal wave interactions in the Southern Ocean. J. Phys. Oceanogr., 50, 3043–3062, https://doi.org/10.1175/JPO-D-20-0001.1.
Cuypers, Y., P. Bouruet-Aubertot, C. Marec, and J.-L. Fuda, 2012: Characterization of turbulence from a fine-scale parameterization and microstructure measurements in the Mediterranean Sea during the BOUM experiment. Biogeosciences, 9, 3131–3149, https://doi.org/10.5194/bg-9-3131-2012.
Cyriac, A., H. E. Phillips, N. L. Bindoff, H. Mao, and M. Feng, 2021: Observational estimates of turbulent mixing in the southeast Indian Ocean. J. Phys. Oceanogr., 51, 2103–2128, https://doi.org/10.1175/JPO-D-20-0036.1.
D’Asaro, E. A., 1995: Upper-ocean inertial currents forced by a strong storm. Part III: Interaction of inertial currents and mesoscale eddies. J. Phys. Oceanogr., 25, 2953–2958, https://doi.org/10.1175/1520-0485(1995)025<2953:UOICFB>2.0.CO;2.
D’Asaro, E. A., C. C. Eriksen, M. D. Levine, C. A. Paulson, P. Niiler, and P. Van Meurs, 1995: Upper-ocean inertial currents forced by a strong storm. Part I: Data and comparisons with linear theory. J. Phys. Oceanogr., 25, 2909–2936, https://doi.org/10.1175/1520-0485(1995)025<2909:UOICFB>2.0.CO;2.
Danioux, E., P. Klein, and P. Rivière, 2008: Propagation of wind energy into the deep ocean through a fully turbulent mesoscale eddy field. J. Phys. Oceanogr., 38, 2224–2241, https://doi.org/10.1175/2008JPO3821.1.
Dong, D., P. Brandt, P. Chang, F. Schütte, X. Yang, J. Yan, and J. Zeng, 2017: Mesoscale eddies in the northwestern Pacific Ocean: Three-dimensional eddy structures and heat/salt transports. J. Geophys. Res. Oceans, 122, 9795–9813, https://doi.org/10.1002/2017JC013303.
Donohue, K., and Coauthors, 2008: Program studies the Kuroshio extension. Eos, Trans. Amer. Geophys. Union, 89, 161–162, https://doi.org/10.1029/2008EO170002.
Fer, I., A. Bosse, B. Ferron, and P. Bouruet-Aubertot, 2018: The dissipation of kinetic energy in the Lofoten Basin eddy. J. Phys. Oceanogr., 48, 1299–1316, https://doi.org/10.1175/JPO-D-17-0244.1.
Fernández-Castro, B., D. G. Evans, E. Frajka-Williams, C. Vic, and A. C. Naveira-Garabato, 2020: Breaking of internal waves and turbulent dissipation in an anticyclonic mode water eddy. J. Phys. Oceanogr., 50, 1893–1914, https://doi.org/10.1175/JPO-D-19-0168.1.
Ferrari, R., and C. Wunsch, 2009: Ocean circulation kinetic energy: Reservoirs, sources, and sinks. Annu. Rev. Fluid Mech., 41, 253–282, https://doi.org/10.1146/annurev.fluid.40.111406.102139.
Garrett, C., 2001: What is the “near-inertial” band and why is it different from the rest of the internal wave spectrum? J. Phys. Oceanogr., 31, 962–971, https://doi.org/10.1175/1520-0485(2001)031<0962:WITNIB>2.0.CO;2.
Garrett, C., and W. Munk, 1979: Internal waves in the ocean. Annu. Rev. Fluid Mech., 11, 339–369, https://doi.org/10.1146/annurev.fl.11.010179.002011.
Gill, A. E., 1984: On the behavior of internal waves in the wakes of storms. J. Phys. Oceanogr., 14, 1129–1151, https://doi.org/10.1175/1520-0485(1984)014<1129:OTBOIW>2.0.CO;2.
Gill, A. E., J. S. A. Green, and A. J. Simmons, 1974: Energy partition in the large-scale ocean circulation and the production of mid-ocean eddies. Deep-Sea Res. Oceanogr. Abstr., 21, 499–528, https://doi.org/10.1016/0011-7471(74)90010-2.
Jaimes, B., and L. K. Shay, 2010: Near-inertial wave wake of Hurricanes Katrina and Rita over mesoscale oceanic eddies. J. Phys. Oceanogr., 40, 1320–1337, https://doi.org/10.1175/2010JPO4309.1.
Jing, Z., P. Chang, S. F. DiMarco, and L. Wu, 2018: Observed energy exchange between low-frequency flows and internal waves in the Gulf of Mexico. J. Phys. Oceanogr., 48, 995–1008, https://doi.org/10.1175/JPO-D-17-0263.1.
Jing, Z., and L. Wu, 2014: Intensified diapycnal mixing in the midlatitude western boundary currents. Sci. Rep., 4, 7412, https://doi.org/10.1038/srep07412.
Jing, Z., L. Wu, and X. Ma, 2017: Energy exchange between the mesoscale oceanic eddies and wind-forced near-inertial oscillations. J. Phys. Oceanogr., 47, 721–733, https://doi.org/10.1175/JPO-D-16-0214.1.
Joyce, T. M., J. M. Toole, P. Klein, and L. N. Thomas, 2013: A near-inertial mode observed within a Gulf Stream warm-core ring. J. Geophys. Res. Oceans, 118, 1797–1806, https://doi.org/10.1002/jgrc.20141.
Karstensen, J., and Coauthors, 2017: Upwelling and isolation in oxygen-depleted anticyclonic modewater eddies and implications for nitrate cycling. Biogeosciences, 14, 2167–2181, https://doi.org/10.5194/bg-14-2167-2017.
Kawaguchi, Y., S. Nishino, J. Inoue, K. Maeno, H. Takeda, and K. Oshima, 2016: Enhanced diapycnal mixing due to near-inertial internal waves propagating through an anticyclonic eddy in the ice-free Chukchi Plateau. J. Phys. Oceanogr., 46, 2457–2481, https://doi.org/10.1175/JPO-D-15-0150.1.
Kawaguchi, Y., T. Wagawa, I. Yabe, D. Ito, T. Senjyu, S. Itoh, and Y. Igeta, 2021: Mesoscale-dependent near-inertial internal waves and microscale turbulence in the Tsushima warm current. J. Oceanogr., 77, 155–171, https://doi.org/10.1007/s10872-020-00583-1.
Kunze, E., 1985: Near-inertial wave propagation in geostrophic shear. J. Phys. Oceanogr., 15, 544–565, https://doi.org/10.1175/1520-0485(1985)015<0544:NIWPIG>2.0.CO;2.
Kunze, E., 1986: The mean and near-inertial velocity fields in a warm-core ring. J. Phys. Oceanogr., 16, 1444–1461, https://doi.org/10.1175/1520-0485(1986)016<1444:TMANIV>2.0.CO;2.
Kunze, E., 1998: A model for vortex-trapped internal waves. J. Phys. Oceanogr., 28, 2104–2115, https://doi.org/10.1175/1520-0485(1998)028<2104:AMFVTI>2.0.CO;2.
Kunze, E., R. W. Schmitt, and J. M. Toole, 1995: The energy balance in a warm-core ring’s near-inertial critical layer. J. Phys. Oceanogr., 25, 942–957, https://doi.org/10.1175/1520-0485(1995)025<0942:TEBIAW>2.0.CO;2.
Leaman, K. D., and T. B. Sanford, 1975: Vertical energy propagation of inertial waves: A vector spectral analysis of velocity profiles. J. Geophys. Res., 80, 1975–1978, https://doi.org/10.1029/JC080i015p01975.
Lee, D.-K., and P. P. Niiler, 1998: The inertial chimney: The near-inertial energy drainage from the ocean surface to the deep layer. J. Geophys. Res., 103, 7579–7591, https://doi.org/10.1029/97JC03200.
Lelong, M.-P., Y. Cuypers, and P. Bouruet-Aubertot, 2020: Near-inertial energy propagation inside a Mediterranean anticyclonic eddy. J. Phys. Oceanogr., 50, 2271–2288, https://doi.org/10.1175/JPO-D-19-0211.1.
Liu, Z., Q. Lian, F. Zhang, L. Wang, M. Li, X. Bai, J. Wang, and F. Wang, 2017: Weak thermocline mixing in the North Pacific low-latitude western boundary current system. Geophys. Res. Lett., 44, 10 530–10 539, https://doi.org/10.1002/2017GL075210.
Lueck, R., and T. Osborn, 1986: The dissipation of kinetic energy in a warm-core ring. J. Geophys. Res., 91, 803–818, https://doi.org/10.1029/JC091iC01p00803.
MacKinnon, J. A., and Coauthors, 2017: Climate process team on internal wave–driven ocean mixing. Bull. Amer. Meteor. Soc., 98, 2429–2454, https://doi.org/10.1175/BAMS-D-16-0030.1.
Martínez‐Marrero, A., and Coauthors, 2019: Near‐inertial wave trapping near the base of an anticyclonic mesoscale eddy under normal atmospheric conditions. J. Geophys. Res. Oceans, 124, 8455–8467, https://doi.org/10.1029/2019JC015168.
Munk, W., 1981: Internal waves and small-scale process. Evolution of Physical Oceanography, N. A. Warren and C. Wunsch, Eds., The MIT Press, 264–291.
Munk, W., and C. Wunsch, 1998: Energetics of tidal and wind mixing. Deep-Sea Res. I, 45, 1977–2010, https://doi.org/10.1016/S0967-0637(98)00070-3.
Olson, D. B., 1980: The physical oceanography of two rings observed by the cyclonic ring experiment. Part II: Dynamics. J. Phys. Oceanogr., 10, 514–528, https://doi.org/10.1175/1520-0485(1980)010<0514:TPOOTR>2.0.CO;2.
Osborn, T. R., 1980: Estimates of the local rate of vertical diffusion from dissipation measurements. J. Phys. Oceanogr., 10, 83–89, https://doi.org/10.1175/1520-0485(1980)010<0083:EOTLRO>2.0.CO;2.
Pallàs‐Sanz, E., J. Candela, J. Sheinbaum, J. Ochoa, and J. Jouanno, 2016: Trapping of the near-inertial wave wakes of two consecutive hurricanes in the Loop Current. J. Geophys. Res. Oceans, 121, 7431–7454, https://doi.org/10.1002/2015JC011592.
Polzin, K. L., 2010: Mesoscale eddy–internal wave coupling. Part II: Energetics and results from polymode. J. Phys. Oceanogr., 40, 789–801, https://doi.org/10.1175/2009JPO4039.1.
Qiu, B., and S. Chen, 2010: Eddy-mean flow interaction in the decadally modulating Kuroshio Extension system. Deep-Sea Res. II, 57, 1098–1110, https://doi.org/10.1016/j.dsr2.2008.11.036.
Rimac, A., J.-S. von Storch, C. Eden, and H. Haak, 2013: The influence of high-resolution wind stress field on the power input to near-inertial motions in the ocean. Geophys. Res. Lett., 40, 4882–4886, https://doi.org/10.1002/grl.50929.
Shcherbina, A. Y., L. D. Talley, E. Firing, and P. Hacker, 2003: Near-surface frontal zone trapping and deep upward propagation of internal wave energy in the Japan/East Sea. J. Phys. Oceanogr., 33, 900–912, https://doi.org/10.1175/1520-0485(2003)33<900:NFZTAD>2.0.CO;2.
Thomas, J., and S. Arun, 2020: Near-inertial waves and geostrophic turbulence. Phys. Rev. Fluids, 5, 014801, https://doi.org/10.1103/PhysRevFluids.5.014801.
Thomas, J., and D. Daniel, 2020: Turbulent exchanges between near-inertial waves and balanced flows. J. Fluid Mech., 902, A7, https://doi.org/10.1017/jfm.2020.510.
Thomson, R. E., and W. J. Emery, 2014: Data Analysis Methods in Physical Oceanography. 3rd ed. Elsevier, 716 pp.
Torres, H. S., P. Klein, D. Menemenlis, B. Qiu, Z. Su, J. Wang, S. Chen, and L.-L. Fu, 2018: Partitioning ocean motions into balanced motions and internal gravity waves: A modeling study in anticipation of future space missions. J. Geophys. Res. Oceans, 123, 8084–8105, https://doi.org/10.1029/2018JC014438.
Waterhouse, A. F., and Coauthors, 2014: Global patterns of diapycnal mixing from measurements of the turbulent dissipation rate. J. Phys. Oceanogr., 44, 1854–1872, https://doi.org/10.1175/JPO-D-13-0104.1.
Weller, R. A., 1982: The relation of near-inertial motions observed in the mixed layer during the JASIN (1978) experiment to the local wind stress and to the quasi-geostrophic flow field. J. Phys. Oceanogr., 12, 1122–1136, https://doi.org/10.1175/1520-0485(1982)012<1122:TRONIM>2.0.CO;2.
Whalen, C. B., J. A. MacKinnon, and L. D. Talley, 2018: Large-scale impacts of the mesoscale environment on mixing from wind-driven internal waves. Nat. Geosci., 11, 842–847, https://doi.org/10.1038/s41561-018-0213-6.
Whalen, C. B., C. de Lavergne, A. C. Naveira Garabato, J. M. Klymak, J. A. MacKinnon, and K. L. Sheen, 2020: Internal wave-driven mixing: Governing processes and consequences for climate. Nat. Rev. Earth Environ., 1, 606–621, https://doi.org/10.1038/s43017-020-0097-z.
Whitt, D. B., and L. N. Thomas, 2013: Near-inertial waves in strongly baroclinic currents. J. Phys. Oceanogr., 43, 706–725, https://doi.org/10.1175/JPO-D-12-0132.1.
Wunsch, C., and R. Ferrari, 2004: Vertical mixing, energy, and the general circulation of the oceans. Annu. Rev. Fluid Mech., 36, 281–314, https://doi.org/10.1146/annurev.fluid.36.050802.122121.
Yang, Q., W. Zhao, X. Liang, J. Dong, and J. Tian, 2017: Elevated mixing in the periphery of mesoscale eddies in the South China Sea. J. Phys. Oceanogr., 47, 895–907, https://doi.org/10.1175/JPO-D-16-0256.1.
Young, W. R., and M. Ben Jelloul, 1997: Propagation of near-inertial oscillations through a geostrophic flow. J. Mar. Res., 55, 735–766, https://doi.org/10.1357/0022240973224283.
Zhai, X., R. J. Greatbatch, and J. Zhao, 2005: Enhanced vertical propagation of storm-induced near-inertial energy in an eddying ocean channel model. Geophys. Res. Lett., 32, L18602, https://doi.org/10.1029/2005GL023643.
Zhai, X., R. J. Greatbatch, and C. Eden, 2007: Spreading of near-inertial energy in a 1/12° model of the North Atlantic Ocean. Geophys. Res. Lett., 34, L10609, https://doi.org/10.1029/2007GL029895.
Zhang, Z., Y. Zhang, W. Wang, and R. X. Huang, 2013: Universal structure of mesoscale eddies in the ocean. Geophys. Res. Lett., 40, 3677–3681, https://doi.org/10.1002/grl.50736.
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., X. Zhang, B. Qiu, W. Zhao, C. Zhou, X. Huang, and J. Tian, 2021: Submesoscale currents in the subtropical upper ocean observed by long-term high-resolution mooring arrays. J. Phys. Oceanogr., 51, 187–206, https://doi.org/10.1175/JPO-D-20-0100.1.
Zhao, Z., M. H. Alford, J. B. Girton, L. Rainville, and H. L. Simmons, 2016: Global observations of open-ocean mode-1 M2 internal tides. J. Phys. Oceanogr., 46, 1657–1684, https://doi.org/10.1175/JPO-D-15-0105.1.
Zhu, R., Z. Chen, Z. Zhang, H. Yang, and L. Wu, 2021: Subthermocline eddies in the Kuroshio Extension region observed by mooring arrays. J. Phys. Oceanogr., 51, 439–455, https://doi.org/10.1175/JPO-D-20-0047.1.