First Evidence of Coherent Bands of Strong Turbulent Layers Associated with High-Wavenumber Internal-Wave Shear in the Upstream Kuroshio

First Evidence of Coherent Bands of Strong Turbulent Layers Associated with High-Wavenumber... www.nature.com/scientificreports Correction: Author Correction OPEN First Evidence of Coherent Bands of Strong Turbulent Layers Associated with High-Wavenumber Internal- Received: 5 May 2017 Wave Shear in the Upstream Accepted: 23 October 2017 Published online: 06 November 2017 Kuroshio 1 2 2 3 Takeyoshi Nagai , Daisuke Hasegawa , Takahiro Tanaka , Hirohiko Nakamura , Eisuke 4 5 6 Tsutsumi , Ryuichiro Inoue & Toru Yamashiro The upstream Kuroshio flows through Okinawa Trough and the Tokara island chain, the region near the continental shelf of the East China Sea and shallow seamounts, where the Kuroshio can induce strong mixing over the shallow topography. Also, tidal currents over the rough topography may produce internal tides, and associated turbulence. The previous observations show energetic high vertical wavenumber near-inertial wave shear in the Kuroshio thermocline, which implies strong turbulent mixing. However, direct turbulence measurements in this region are very scarce. Using high lateral resolution (1–2 km) direct turbulence measurements, we show here, for the first time, that strong turbulent layers form spatially coherent banded structures with lateral scales of >O(10 km), associated with bands of near-inertial wave/diurnal internal tide shear of high vertical wavenumber in the upstream Kuroshio. The turbulent kinetic energy dissipation rates within these turbulent layers −7 −1 −4 2 −1 are >O(10 W kg ), and estimated vertical eddy diffusivity shows >O(10 m s ) on average. These results suggest that the high vertical wavenumber near-inertial waves propagating in the upstream Kuroshio could have large impacts on the watermass modifications, momentum mixing, nutrient supply, and associated biogeochemical responses in its downstream. The Kuroshio, the western boundary current of the subtropical gyer in the North Pacific, has been known to transport large amounts of heat and salt from the tropical ocean . Because the Kuroshio is inherently an ocean dynamic front, where the associated jet flows along the interface between seawaters of different water properties , mixing processes in the Kuroshio Front play very important roles in watermass formation, transformation, and 3–5 subduction . In addition to this physical or hydrographical important role, the mixing in the Kuroshio has profound implications for biogeochemistry because the Kuroshio carries waters of relatively high nutrient con- 6,7 centrations from the tropics through its subsurface layers, as known as the nutrient stream similar to the Gulf 8–11 8,9 12,13 Stream . Earlier studies in the Gulf Stream and relatively recent studies in the Kuroshio have reported that elevated nutrient concentrations are found along these western boundary currents in their subsurface layers compared to that in the ambient waters of the same density. Because along-isopycnal stirring by mesoscale eddies tends to homogenize any tracer anomaly on the density surface in the ocean interior , the elevated nutrient concentrations at the subsurface of the Kuroshio are most likely formed by across-density/diapycnal fluxes by Tokyo University of Marine Science and Technology, Department of Ocean Sciences, Tokyo, 108-8477, Japan. Tohoku National Fisheries Research Institute, Japan Fisheries Research and Education Agency, Fisheries Oceanography and Resources Department, Shiogama, Miyagi, 985-0001, Japan. Kagoshima University, Faculty of Fisheries, Kagoshima, 890-0056, Japan. Kyushu University, Research Institute for Applied Mechanics, Kasuga, Fukuoka, 816-8580, Japan. Research and Development Center for Global Change, Japan Agency for Marine-Earth Science and Technology, Yokosuka, 237-0061, Japan. Graduate School of Science and Engineering, Kagoshima University, Kagoshima, 890-0065, Japan. Daisuke Hasegawa, Takahiro Tanaka, Hirohiko Nakamura, Eisuke Tsutsumi, Ryuichiro Inoue and Toru Yamashiro contributed equally to this work. Correspondence and requests for materials should be addressed to T.N. (email: tnagai@kaiyodai.ac.jp) ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 1 www.nature.com/scientificreports/ irreversible diabatic turbulent diffusion. A previous study speculated that the turbulent diffusive nutrient flux by geostrophic vertical shear contributes to the concentration of elevated nutrients below the Gulf Stream. However, a limited number of direct turbulence observations at that time did not provide the supporting evidence for this 15–17 hypothesis . On the other hand, the Kuroshio flows through shallow topography in its upstream along Okinawa Trough and near Tokara Strait, where the Kuroshio over the topography can induce strong turbulence and diapycnal mix- 18–20 ing . While a number of previous studies have reported that the beams of the M internal tide emanated from 21,22 the rough topography induced strong turbulence in close proximity to canyons and seamounts , recent ship- board and lowered Acoustic Doppler Current Profiler (ADCP) measurements in the upstream Kuroshio showed widespread large amplitude near-inertial internal wave shear of high vertical wavenumber in the Kuroshio and regions between continental margins and the Kuroshio . A very recent field study in Tokara Strait reported near-inertial waves of even shorter vertical wavelengths ~50 m accompanied by pronounced turbulent dissipa- −7 −6 −1 20 tion rates of O(10 –10 W kg ) . One of the most likely mechanisms to form the high vertical wavenumber shear is the near-inertial internal waves generated by the Kuroshio over shallow topography in Okinawa Trough and Tokara Strait . Because the Kuroshio is forced to turn its direction from northeast to southeast when it approaches Kyushu island, high vertical wavenumber near-inertial internal waves could also be generated spon- 24–28 taneously by the meandering Kuroshio . Furthermore, the upstream Kuroshio passes the so-called critical latitude, 28.9°N, where M internal tide energy can be efficiently converted to that of high vertical wavenumber 29,30 near-inertial waves through parametric subharmonic instability (hereinae ft r PSI) . Near the critical latitude, K or O diurnal internal tides generated over the rough topography are also near-inertial internal waves. In a 1 1 strong baroclinic front, where geostrophic shear allows anomalously low frequency internal waves, a recent the- oretical study showed that PSI of wind-driven near-inertial waves could generate high wavenumber internal waves with a frequency of half the local Coriolis frequency f/2. Regardless of their generation mechanism, these propagating near-inertial internal waves can be trapped on the anticyclonic side of the front where the lowest 32–34 internal wave frequency is decreased by geostrophic lateral and vertical shears . The trapped near-inertial 35–37 waves can break into three-dimensional turbulence and enhance mixing of momentum and tracers . However, direct turbulence measurements in the upstream Kuroshio are very scarce. Also, the lateral resolutions of previous turbulence measurements are not sufficient to resolve the coherent spatial patterns associated with high wave- 15–17,38 number near-inertial internal wave shear . As a result, earlier studies did not find elevated turbulence in 15–17 the thermocline of the Gulf Stream , and concluded that turbulence in the ocean interior is caused by patchy random internal wave breaking, even in the Gulf Stream. In contrast, a number of recent studies have reported 27,39–42 43 enhanced turbulent mixing in the thermocline of the Kuroshio and the Gulf Stream caused by propagating internal waves. However, the lateral resolutions of the turbulence data were still insufficient to provide direct evi- dence that the turbulence is caused coherently by banded internal-wave shear in the thermocline of the western 44,45 boundary currents. Although a recent extensive microstructure time-series survey at low latitude reported that the strong turbulence coincided with the relatively large strain caused by wind-induced near-inertial waves of anomalously low frequency in the anticyclonic vorticity field , it is not known if near-inertial internal waves form the spatially coherent strong turbulent layers in ocean fronts. In this study, a new turbulence and microstructure observation scheme is developed and utilized to measure turbulence directly with a high lateral resolution of 1–2 km in the upstream Kuroshio near Tokara Strait, where a large amount of energy from the Kuroshio can be dissipated. For the first time, the observations in this study suggest that propagating high vertical wavenumber near-inertial internal waves form spatially coherent banded layers of strong turbulence in the upstream Kuroshio, with O(100 m) and >O(10 km) vertical and lateral scales, respectively. Results Observations. e Th in-situ survey was conducted using the R/T/V Kagoshima-maru during November 12–20, 2016, near Tokara Strait where the upstream Kuroshio flows over the shallow seamounts and through the narrow straits near the Tokara island chain off Kyushu island (Fig.  1). Two mooring systems, which included an upward looking ADCP (75 kHz 20° beam angle Workhorse Long Ranger ADCP, Teledyne RDI) were deployed near the Kuroshio and measured horizontal currents and directions during November 13–19, 2016, at Stn. M01 (30.35°N, 129.85°E) and M02 (30.00°N, 129.50°E, Fig. 1a). Two transect observations were carried out to obtain current velocity data using a shipboard ADCP (75 kHz 30° beam angle, Ocean Surveyor, Teledyne RDI) and high-resolu- tion towed microstructure data. The transect observations (hereinaer L ft eg A) were conducted on November 14 from west to east closely along the zonal direction for the first half of the transect, and continued in the northeast direction for the rest (thick blue line in Fig. 1b). The transect surveys were repeated (hereinafter Leg B; red in Fig. 1b) on November 18–19 close to the same ship track as for the second half of the Leg A observations. During these two transect observations, a new tow-yo Underway Vertical Microstructure Profiler (Underway-VMP, here- inaer UVMP ft ; Rockland Scientific International, Victoria, Canada) was used in the Kuroshio Front to investigate high-resolution spatial structures of microscale turbulence. Prior to the transect observations of Legs A and B, the Underway Conductivity Temperature Depth profiler (hereinaer UCTD; T ft eledyne Oceanscience, California, USA) was deployed across front direction repeatedly (thick black line in Fig. 1a). Conductivity-Temperature- Depth (CTD) profiles also were made at 36 stations (Fig.  1a). Because the UVMP used in this study was a pro- totype with no conductivity sensor, simultaneously gathered salinity data are not available. However, because water properties obtained by these UCTD and CTD in the observed regions suggest a close relationship between temperature and density (Supplementary Fig. S1), we estimated density in Legs A and B from temperature alone (see Methods section). In addition to the mooring observations made during the shipboard surveys, the relatively long-term moored ADCP records from September 30, 2000 through July 1, 2001 obtained at Stn. TK1 (30.13°N, 130.19°E, Fig. 1a) in the vicinity of the shipboard observation sites were analyzed to investigate the mechanisms generating observed internal waves. ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 2 www.nature.com/scientificreports/ Figure 1. Observation site in Tokara Strait near the Tokara island chain off Kyushu island. (a) Underway-CTD (UCTD) observation line is shown as a thick black line, and Underway-VMP (UVMP) observation lines are in blue for Leg A and red for Leg B. Mooring stations Stn. M01 and M02 are indicated by black (+). CTD stations are denoted by solid white circles (). Total number of the CTD profiles is 36, which consists of 13 profiles each at Stn. M01 and M02 during time series observations and 10 profiles at four other stations near the Kuroshio. e r Th elatively long-term mooring Stn. TK1 is shown as a black triangle. Nakanoshima, and Tanegashima tidal stations are shown as a red diamond, and a yellow square, respectively. Black contours are AVISO sea surface height (m) averaged from November 14–19, 2016. (b) Enlarged map of the UVMP observation site. UVMP in Leg A is denoted by blue diamonds and Leg B by red triangles. Black crosses (×) indicate six XBT stations during Leg A. Color shadings indicate depth (m). The map is created with MATLB R2016b with ETOPO1 topography data available from NGDC NOAA https://www.ngdc.noaa.gov/mgg/global/. High-vertical-wavenumber internal-wave shear. Measured horizontal current velocity using the ship- board ADCP mostly shows a northeastward flow associated with the Kuroshio, consistent with the sea surface height (Fig. 1). During the Leg A observation, this flow was strongest in the westernmost edge of the transect, where the main stream of the Kuroshio flowed (Fig.  2a). On the other hand, in Leg B, the largest velocity magni- tude was recorded at the northeastern part, the last half of the transect (Fig. 2d), as the continuous UVMP profil- ing in Leg B started from the southwestern region to the northeastern direction, reaching the main stream of the Kuroshio Front in the last quarter of the Leg B (Fig. 1). In both the Legs, current velocity magnitudes are largest near the surface, due to the baroclinic flow of the Kuroshio, with some exceptions where the stronger currents are found subsurface, especially on the anticyclonic side of the Kuroshio (Fig. 2a,d). Back-rotated velocity shear, which is the vertical gradient of the lateral velocity rotated horizontally, assuming the inertial rotation with respect to the reference time (see Methods section), shows bands of alternating positive and negative signs, close to the density surface (Fig. 2b,c,e,f ), which are similar to a study by Rainville and Pinkel . These banded shears are typical characteristics of the propagating inertia-gravity waves in a stratified ocean. The wavelengths of these banded shear structures are O(100 m) in vertical and O(10 km) in lateral direction, similar to Rainville and Pinkel’s study in this region . With λ and λ , horizontal and vertical wavelengths (m), respectively, Nb =∂ /∂z h v −1 −2 −2 buoyancy frequency (rad s ) (where b = −gρ/ρ is buoyancy (m s ) with gravitational acceleration g (m s ), −3 −3 seawater density ρ (kg m ) and reference density ρ = 1025 (kg m ), and vertical coordinate z(m)), and f Coriolis −1 frequency (rad s ), the internal-wave dispersion relation, then implies that the frequency of these waves 22 22 2 −2 −1 −4 −1 ω =+ (/ fN λλ /) λ ~ f , based on the estimated N~O(10 rad s ), and local f~O(10 rad s ). The zh z internal-wave ray paths computed, assuming a quiescent condition show that the angles of the internal wave rays with near-inertial frequencies are more consistent with the observed shear bands than those with the M tidal frequency (Fig.  2b,c,e,f ). This also suggests that the observed high-wavenumber velocity shear is due to near-inertial internal waves. The total shear vertical wavenumber spectra averaged for the entire transect of Leg A and B suggest that shear variance is significantly larger than that of Garrett-Munk (GM) internal-wave equilib- 46,47 rium spectra by several factors, at vertical wavelengths larger than 100 m (Fig. 3c,d), at which the resolution of the shipboard ADCP starts to fail to resolve shear of shorter wavelengths. Computed Richardson numbers 22 2 RiN =+ /(uv ), where u and v are zonal and meridional shear, show relatively small values in the banded zz z z shear layers and near the surface. However, because of the coarse resolution of the ADCP shear, the Ri are mostly larger than the critical value for the Kelvin-Helmholtz instability Ri = 0.25 except near the surface and below 500 m in Leg A and 300 m in Leg B (Supplementary Fig. S2). e Th rotary shear spectra in Leg A show that the var - iance of clockwise rotating shear with depth is larger than that of anticlockwise rotating component (Fig. 3c), suggesting the dominant downward energy propagating internal waves. In Leg B, the clockwise and anticlockwise components of rotary shear spectra are comparable with slightly larger shear variance for the anticlockwise ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 3 www.nature.com/scientificreports/ Figure 2. Vertical sections of the shipboard Acoustic Doppler Current Profiler velocity measurements. (a–c) are for Leg A (blue in Fig. 1b) and (d–f) are for Leg B (red in Fig. 1b). Horizontal absolute current velocity −1 −1 (m s ) are shown in (a,d). Back-rotated shear [s ] is shown for zonal shear u (t ) (b,e), and for meridional z 0 −3 shears v (t ) (c,f) (see Methods section). Black contours are σ (kg m ). Magenta curves in (b,c,e,f) are internal- z 0 θ wave ray paths at frequencies of 1.01 f, 1.1 f, 1.2 f, 1.4 f, and M tidal frequency, where f is Coriolis frequency, as indicated in the panels. Red triangles in upper abscissa are the locations at which hodographs of shear are plotted for each Leg in Fig. 3. rotating component (Fig. 3d). The ratio of the integrated anticlockwise shear spectrum to that for the clockwise uu u u z z zz component (where and are anticlockwise and clockwise components of shear spec- φφ dk/ dk φ φ ∫∫ ACWCW ACW CW trum, respectively, and k is the vertical wavenumber) as a function of longitude suggests that large shear variance found at 129.3–129.575°E in Leg A is associated with the downward energy propagating internal waves, and that found at 129.4–129.55°E in Leg B is caused by upward energy propagating internal waves (Fig. 3e,f ). The hodo- graph of the shear at 129.447°E in Leg A, where shear magnitude is large, shows that shear is rotating clockwise with depth, suggesting downward energy propagating internal waves, while that at 129.472°E in Leg B shows anticlockwise rotation for upward energy propagating internal waves, consistent with the trends found in the −3 −3 shear spectra (Fig. 3a,b,e,f ). The estimated average internal wave energy is 6.4 J m for Leg A, and 14.2 J m for −5 −1 Leg B. The estimated mean group velocity, 3.5 × 10 m s with the mean vertical wavelength of the internal −2 −2 wave, 150 m, results in the mean vertical energy flux of 0.2 mW m for Leg A, and 0.5 mW m for Leg B (see Methods section). These estimates are comparable to the rate of wind energy input to inertial motions computed 48–50 by a slab model (see Methods section, Supplementary Fig. S3). es Th e results indicate that the high vertical wavenumber shear is caused by near-inertial internal waves prop- agating in the upward and downward directions in the Kuroshio near Tokara Strait. While the wind energy flux to the inertial currents is found to be sufficient to generate an observed amplitude of near-inertial waves, it is not clear how the observed high vertical wavenumber near-inertial internal waves are formed. We speculate that the possible mechanisms are (1) PSI of the M internal tide upstream, which was advected to the observation sites by the Kuroshio flow, (2) local PSI with anomalously lower minimum internal-wave frequency than f due to mean velocity shears, (3) K internal tides generated over rough topography, (4) near-inertial internal waves generated ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 4 www.nature.com/scientificreports/ −1 Figure 3. Hodographs of the back-rotated shear [s ] are shown for (a) 129.447°E in Leg A (Fig. 2a–c) and for (b) 129.472°E in Leg B (Fig. 2d–f ). Shading indicates depth (m). Average vertical wavenumber shear spectra −2 −1 [s cp m ] are shown for (c) Leg A and for (d) Leg B. Solid black lines (CW) are the clockwise rotating component with depth and dashed lines (ACW) are the anticlockwise rotating component. Solid black-white- black lines are for total shear spectra. Dotted lines (GM) are Garrett-Munk equilibrium internal wave shear spectra. (e,f) Ratio of integrated anticlockwise shear spectra to that for the clockwise component uu u u zz z z φφ dk/ dk (where φ and φ are anticlockwise and clockwise components of shear spectrum, ∫∫ ACWCW ACW CW respectively and k is the vertical wavenumber) as a function of longitude in Leg A (e) and B (f). ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 5 www.nature.com/scientificreports/ Figure 4. Microscale shear and shear spectra obtained from tow-yo observations. (a) One of the turbulent −1 −2 −1 shear profile [s ] and (b,c) turbulent shear spectra [s cp m ]. The shear spectra in red and blue correspond to the depth ranges indicated by the same color in (a). The Nasmyth empirical shear spectra are shown in (b,c) as −1 black curves for each computed turbulent kinetic energy dissipation rate ε [W kg ]. by the Kuroshio over shallow topography in Okinawa Trough and Tokara Strait, (5) spontaneously generated near inertial internal waves by the Kuroshio, and (6) PSI of wind-induced inertial waves to generate high wavenumber anomalously low frequency internal waves of half the Coriolis frequency, f/2 along isopycnal in a strong baroclinic front . High-resolution tow-yo turbulence measurements. The shipboard ADCP measurement shows the large amplitude banded high vertical wavenumber shear, similar to previous studies . Here, we show, for the first time, that these large amplitude near-inertial internal wave shears generate strong turbulence of spatially coherent banded structures, employing a newly developed high-resolution (1–2 km) tow-yo microstructure profiling tech- nique, using 1.5 mm dyneema rope. The thin rope allows the microstructure profiler (VMP-250) to sink relatively smoothly with small drag down to 300 m depth. The measured turbulent shear spectra by the tow-yo profiling show overall good agreement with the Nasmyth empirical shear spectrum (Fig. 4). This novel technique, there- fore, allows us to resolve spatial structures of intermittent turbulence, which have been unrecognized because of the coarse resolution measurements by the previous microstructure profiling. −7 −1 In the upper 100 m, the measured turbulent kinetic energy dissipation rate, ε is large, O(10 W kg ), because of turbulence within the surface boundary layer, which is forced by momentum and heat fluxes at the surface (Fig. 5). In the stratified layers below 100 m depth, away from the surface boundary layer, the dissipation rates −7 −1 −6 −1 still frequently exceed O(10 W kg ), and reach O(10 W kg ) at some locations (Fig. 5a,b). These dissipation rates are 100 to1000 times greater than those found typically in the open water thermocline. More importantly, these large dissipation rates form coherent along-isopycnal bands of strong turbulent layers. These strong turbu- lent layers found beneath and on the less dense side of the Kuroshio correspond very well to the large amplitude near-inertial wave shear bands, suggesting that the observed near-inertial internal waves of high vertical wave- number break into microscale turbulence. The comparison between our observed turbulent dissipation rates and 52,53 previous internal-wave parameterization for energy dissipation rates show a positive correlation (r = 0.63, p-value = 0.0001) within the same orders of magnitude (Supplementary Fig. S4). For the first time, our high ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 6 www.nature.com/scientificreports/ −1 Figure 5. Vertical section of (a,b) measured turbulent kinetic energy dissipation rates log ε [W kg ], and 2 −1 (c,d) estimated vertical turbulent eddy diffusivities log K [m s ] for (a,c) Leg A and (b,d) Leg B. Contours are σ . resolution microstructure surveys provide evidence that strong turbulent layers with spacial scales of O(100 m) and O(10 km) in the vertical and lateral directions, respectively, are coherently formed along the near-inertial internal-wave shear bands (Figs 2 and 5a,b). To investigate the significance of this turbulence not only for energy dissipation but also for the vertical mixing and diffusion, the vertical turbulent eddy diffusivity is estimated using measured dissipation rates ε and buoyancy frequency N with the previous diffusivity model , namely, K = γ , (1) 54 −2 2 −1 where γ = 0.2 is a mixing efficiency factor . The estimated vertical turbulent eddy diffusivities are > O(10 m s ) within the surface boundary layer, which is not surprising considering continuous forcing at the surface with weak stratification in the surface boundary layer (Fig.  5c,d). The vertical eddy diffusivities in the stratified subsur - −4 2 −1 −3 2 −1 face layers, however, are notably high, O(10 m s ) on average, and >O(10 m s ) at many locations below 100 m depth (Fig. 5c,d). The high diffusivity layers are structured in a banded form, reflecting the spatial struc- tures found in the dissipation rates. Discussion Early microstructure surveys in the 1970s–1980s did not find strong turbulence across the western boundary 15–17 current, i.e., the Gulf Stream, but found intermittent turbulent patches due to random internal wave breaking . Relatively recent microstructure measurements in the upstream Gulf Stream, Florida Current , and in the down- stream Gulf Stream , reported the possibility of the internal wave induced turbulence, but whether turbulent layer is formed coherently with the banded internal wave shear is unclear due to the coarse lateral resolution of the surveys. In contrast, this study found coherent bands of very strong turbulence clearly associated with the large amplitude, high vertical wavenumber internal-wave shear, in the upstream Kuroshio. ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 7 www.nature.com/scientificreports/ e o Th bserved strong turbulence associated with high vertical wavenumber internal wave shear implies active occurrence of PSI near the critical latitude 28.9°N. The moored ADCP observations during our shipboard exper - iment suggest that the velocity field is dominantly modulated at M semi-diurnal tidal frequency (Supplementary Fig. S5a,b). In contrast, the moored ADCP vertical shear modulates at subinertial, near-inertial (or K diur- nal), as well as M semi-diurnal tidal frequencies. However, these M tidal, near-inertial and subinertial peaks 2 2 in shear spectra are not clear, especially in upper layers, probably because of the Doppler smearing by the strong Kuroshio flow (Supplementary Fig. S5). The energy of the observed M semi-diurnal internal tide could be converted to that of high vertical wavenumber near-inertial waves through PSI. To investigate the possibility of PSI and other mechanisms to generate observed high vertical wavenumber internal wave shear, the long-term moored ADCP data obtained at Stn. TK1 from September 2000 through July 2001 were analyzed. Assuming the energy source of the PSI-induced near-inertial waves is M internal tide, the amplitude of the PSI-induced near-inertial wave shear should modulate fortnightly with spring-neap tidal cycle , similarly to that of M inter- nal tides. The harmonic analysis of the tidal elevation at Nakanoshima (29.84°N, 129.85°E, Fig.  1a), suggests that the diurnal tidal amplitude (D = K + O ) is about the half of that estimated for the semi-diurnal com- 1 1 1 ponent (D = M + S + N ) (Supplementary Fig. S6b). The 30-h lowpass near-inertial shear amplitude exhibits 2 2 2 2 subinertial variations with larger amplitude found in the upper 250 m (Supplementary Fig. S6c). Most of these subinertial variations are found to occur fortnightly, with a period of 14 days, throughout the water column (Fig. 6a,b). The fortnightly modulation found in the near-inertial shear amplitude may imply occurrence of the PSI . However, our observation sites are very close to the critical latitude, where the frequency of diurnal tidal flow is at the near-inertial frequency. It is also possible that diurnal near-inertial internal tides are generated directly from the nearby rough topography, which then propagate to the observation sites. To elucidate whether fortnight modulations of near-inertial shear amplitude are correlated more with those of semi-diurnal D or diurnal D tidal amplitude, the correlation analyses between D or D tidal amplitude and near-inertial shear 1 2 1 amplitude are conducted (Supplementary Fig. S7). The correlation between D diurnal tidal amplitude and the 30-h lowpass near-inertial shear amplitude is significantly high with a positive Spearman correlation coefficient of 0.59. On the other hand, the correlation with D semi-diurnal tidal amplitude is low, at 0.05. The results suggest that although the near-inertial shear amplitude modulates at the fortnight frequency, the dominant source of the near-inertial energy is D diurnal tides, not D semi-diurnal internal tides through PSI. Our study estimates 1 2 −2 10–100 mW m for the K internal tide energy generation rate at the bottom near Tokara Strait, based on a regional numerical simulation which could be sufficiently large to generate an observed vertical energy ux fl of −2 the near-inertial internal waves, 0.1–1 mW m . Accordingly, the local PSI of M internal tides, with the anoma- lously low minimum internal wave frequency due to geostrophic shears, is unlikely. However, if the PSI-induced near-inertial internal waves are advected by the Kuroshio from the upstream region south of the critical latitude, the fortnightly modulations of the Doppler-smeared near-inertial shear amplitude should be lagged behind that of the D semi-diurnal tide. The lagged Spearman correlation coefficients between the D semi-diurnal tidal 2 2 amplitude and the 30-h lowpass near-inertial shear amplitude show a maximum correlation coefficient of 0.34 with a 94 h time lag (Supplementary Fig. S7). With a distance from the critical latitude region to the observa- −1 tion sites along the Kuroshio, ~160 km, the time lag of 94 h implies an average advection speed at 0.47 m s . The estimated advection speed is consistent with the velocity magnitude observed in the subsurface layers of −1 the Kuroshio, ~0.5 m s between 200–400 m depth (Fig. 2a,d). Assuming that the estimated near-inertial wave −3 energy, O(10 J m ) (see Methods section) is continuously dissipated by the observed bands of the turbulent −8 −7 −1 −3 dissipation rates, O(10 –10 W kg ), with the reference seawater density ρ = 1025 kg m , the internal waves could be dissipated completely aer ft 1–10 days. er Th efore, it is possible that the PSI-induced near-inertial waves can survive ~4 days after the 94-h advection. Even without the tide or PSI, near-inertial internal waves of high vertical wavenumber could be generated by the Kuroshio over the topography . The globally estimated energy generation rate from geostrophic flow to lee wave by the study of Nikurashin and Ferrari in Okinawa Trough, −2 the region near the observation sites, is O(0.1–1 mW m ), consistent with the vertical energy flux of observed near-inertial internal waves. The local generation of the D diurnal internal tide is the most plausible and likely candidate for the mech- anism generating high vertical wavenumber near-inertial internal waves observed near Tokara Strait. The same harmonic analysis for the diurnal and semi-diurnal tidal elevations during the shipboard observations in November 2016 at Tanegashima tidal station (30.47°N, 130.97°E, Fig. 1a) shows a relatively large diurnal tidal amplitude during the UVMP transect observations in Leg A (November 14–15) and Leg B (November 18–19) (Supplementary Fig. S8c). The near-inertial waves induced by the Kuroshio over the topography cannot be ruled out. Other possibilities, such as the spontaneous generation of the near-inertial waves by the meandering 24–28 31 Kuroshio , and the PSI of wind-induced near-inertial waves are not discussed in the present study due to the limited observation data. More intensive field campaigns, and high-resolution numerical simulations, are necessary to definitively identify the source of these near-inertial waves propagating in the upstream Kuroshio. e b Th reaking of the near-inertial waves is influenced by the Kuroshio. This is because the near-inertial waves can be trapped in the region of negative relative vorticity on the less dense side of the Kuroshio and of strong geo- 32–34 strophic vertical shear below the Kuroshio . The modulation of the lowest internal wave frequency, as known as the effective Coriolis frequency f , by the vertical component of the relative vorticity in the Kuroshio can be eff equated for a meridional two-dimensional front,   ∂v −1 ff = f + f , eff   ∂x (2) where v is the mean meridional velocity and x is the zonal coordinate . Considering the contribution from ver- tical shear of the mean flow to draw down the lowest internal wave frequency further, the equation becomes ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 8 www.nature.com/scientificreports/ −2 Figure 6. Long-term moored ADCP data analysis. The variance preserved spectra [s ] of the 30-h lowpass near-inertial shear (a) as a function of depth [m] and frequency [Hz], and (b) the depth averaged uv + zz spectrum. (c) The time series of the amplitude of 30-h lowpass diurnal tidal elevation, D = K + O [m] at 1 1 1 Nakanoshima tidal station (Fig. 1a) is shown in red, and that for semi-diurnal tidal elevation D = M + S + N 2 2 2 2 is shown in black. The diurnal and semi-diurnal constituents are obtained by harmonic analysis for the hourly tidal elevation record at Nakanoshima tidal station. (d) The time series of 30-h lowpass near-inertial shear amplitude measured at Stn. TK1. uv + zz 21 −       ∂v ∂v ∂b −1    ω = ff  + f −  .    min                ∂x ∂z ∂z (3) es Th e lowest internal wave frequencies are computed with the measured ADCP velocity and buoyancy along the transect averaged over 200 m depth and 7 km for vertical and horizontal direction, respectively, replacing v with ∂v the computed velocity component normal to the ship track, x by the direction parallel to the ship track, and () ∂z ∂u ∂v by total velocity shear square . Because our ship tracks were not normal to the front, the contribu- () () ∂z ∂z ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 9 www.nature.com/scientificreports/ Figure 7. Normalized minimum internal wave frequencies by local Coriolis frequency averaged within the upper 500 m as a function of longitude for (a) Leg A and for (b) Leg B. Blue indicates the minimum internal 33 34 wave frequency for a barotropic front f /f (2) , and red indicates that for a baroclinic front ω /f (3) . Shading eff min represents 95% confidence interval for each average value. Horizontal dashed line indicates unity. tion from the vertical relative vorticity to modulate f and ω is most likely underestimated. Nevertheless, both eff min the average f and ω within the upper 500 m show a similar trend in which the lowest internal wave frequencies eff min are significantly decreased around 129.3–129.5°E in Leg A (Fig.  7a) and 129.5–129.65°E in Leg B (Fig. 7b) from the local Coriolis frequencies f. The reduction of the lowest internal-wave frequency and the slightly larger lowest internal-wave frequency formed to the west, create the convex-down structure in f and ω at several longitude eff min ranges, where the near-inertial waves can be trapped. In these regions, the amplitudes of high vertical wavenum- ber shear seems large, especially at 129.3–129.6°E in Leg A (Figs 2b,c and 7a). The large turbulent dissipation rates are found beneath the Kuroshio at 100–150 m depth and on the less dense side of the Kuroshio below 150 m depth at 129.3–129.6°E in Leg A and at 129.5–129.65°E in Leg B (Fig. 5), where the near-inertial wave trapping favorable structures are coincidently formed (Fig. 7). This implies that the trapping mechanisms by the geostrophic shear may be effective in promoting near-inertial internal-wave breaking and causing enhanced turbulence in the upstream Kuroshio. e in Th duced energy dissipations and vertical mixing could have important implications for watermass modi- 57 −3 fications. A previous study estimated the vertical eddy diffusivity at the salinity minimum layer σ = 26.8 kg m −4 2 −1 of North Pacific Intermediate Water (NPIW) in the Okinawa Trough to be 6.8–21.5 × 10 m s using a steady-state salinity balance with historical current data and salinity data obtained by the profiling floats. Although the direct turbulence data obtained in this study cover only the upper 300 m depth, and did not reach −3 the layer of σ = 26.8 kg m , strong turbulence just above salinity minimum layer of NPIW, which is equivalent to −4 2 −1 57 >O(10 m s ) of average eddy diffusivity, is consistent with a previous estimate . This implies that the observed large amplitude high vertical wavenumber internal waves play very important roles in watermass mixing and modification in the study region. Not only is the observed pronounced turbulent diapycnal mixing important for the watermass modifica- tions of NPIW, it is also very important for nutrient supply from the lower layer to the surface euphotic layer. Recent studies of the nutrient transport in the upstream Kuroshio showed that despite its low-nutrient olig- −1 otrophic surface water, the Kuroshio transports a large amount of nitrate at an average rate of 170.8 kmol s −3 6,7 in the subsurface layer, σ = 24–27 kg m and transports it downstream , similar to the nutrient stream of 8–11 the Gulf Stream . In addition, nitrate concentrations below the Gulf Stream are found to be higher than 8,9 in the ambient water of the same density . This along-isopycnal anomaly of nutrients is most likely formed and maintained by diapycnal diffusion, as opposed to homogenizing mesoscale along-isopycnal stirring. −4 2 −1 Therefore, large average turbulent diffusivity of O(10 m s ) in the subsurface layer at 100–300 m depth, −3 σ = 24–26 kg m , reported by this study, indicates that the observed subsurface turbulence induced by the −3 near-inertial internal waves could effectively diffuse up the subsurface nutrient stream of σ = 24–27 kg m to the shallower layer, generate the elevated nutrient concentrations along the Kuroshio, and increase pri- mary production in its downstream. In this study, for the first time, the spatially coherent banded layers of strong turbulence associated with high vertical wavenumber near-inertial internal waves were directly observed using a new tow-yo microstructure pro- filer. Our analyses suggest that the high vertical wavenumber near-inertial waves are most likely generated by the diurnal tide over the rough topography in Tokara Strait. Because the observation site is around 30°N latitude, where the inertial period is 24 h, the generated diurnal internal tides are inherently near-inertial waves with slow group velocities. The ADCP and microstructure data suggest that these near-inertial internal waves are trapped and broken into turbulence, preferentially on the anticyclonic side of the upstream Kuroshio, suggesting that the trapping of the near-inertial waves by the Kuroshio due to anticyclonic vorticity catalyzes the breaking and dis- sipation of the near-inertial/diurnal internal waves, and promoting watermass modification and nutrient supply. Although the contributions of PSI to form observed high vertical wavenumber near-inertial waves were found to be small compared to the diurnal internal tide in our observations, other mechanisms to generate observed ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 10 www.nature.com/scientificreports/ inertial waves are not discussed in detail in this study. Okinawa Trough and Tokara Strait are very peculiar regions because of their strong baroclinic current, the Kuroshio (which can modify the lowest internal wave frequency), flowing through the shallow seamounts and islands, and through the critical latitude for PSI and the boundary where the diurnal period is identical to the inertial period. More comprehensive field and numerical studies are, therefore, necessary to fully understand the generation, propagation, and dissipation mechanisms of these near-inertial internal waves/diurnal internal tides propagating in the upstream Kuroshio. Methods Back-rotation of the shear. Near-inertial wave velocity and shear rotate clockwise in time in the Northern Hemisphere. When the observation period exceeds or is equivalent to the inertial time scale, this rotation needs to be removed to avoid temporal aliasing of lateral velocity and shear structures. The back-rotated shear Z(t ) = u (t ) + iv (t ) (where u and v are the zonal and meridional shear, respectively, and i is an imaginary unit) 0 z 0 z 0 z z referenced to the time t , is estimated from the observations Z(t) = u (t) + iv (t) at the arbitrary time t as 0 z z it () −tf Zt () = Zt ()e , (4) where f is the mean Coriolis frequency, assuming the rotation at the mean inertial frequency, and back-rotated zonal u (t ) and meridional shear v (t ) are obtained by taking the real and imaginary parts of (4), respectively. z 0 z 0 Tow-yo turbulence survey. We used a new tow-yo Underway Vertical Microstructure Profiler (Underway-VMP: UVMP) to measure high-resolution spatial structures of microscale turbulence at the Kuroshio Front. The UVMP consists of a vertical microstructure profiler, VMP-250 (Rockland Scientific International, Victoria, Canada) and a winch for the Underway-CTD (UCTD, Teledyne Oceanscience, USA). The VMP-250 carries two shear probes, two FP07 thermistors, a pressure sensor, and a vibration sensor to measure microscale velocity shear, microscale temperature gradient, pressure and accelerations, respectively, at a frequency of 512 Hz. −1 e UVMP wa Th s tow-yoed at a ship speed relative to water, at 1.5–2.1 m s . The descending speeds of the UVMP, −1 gradually decreasing with depth, were between 0.3–0.8 m s . Taking 7 minutes to reach 300 m depth and 5 min- utes for recovery, the lateral resolution of the tow-yo survey was 1–2 km, depending on the absolute current veloc- ity. The data were internally recorded and recovered ae ft r each transect observation. However, because the both FP07 thermistors mounted initially on the VMP-250 in Leg A malfunctioned, the expendable bathythermographs (XBT-T7s) were deployed six times to obtain temperature profiles. The XBT temperature data were interpolated by objective mapping with decorrelation lengths of 15 km and 40 m along ship-track and vertical directions, respectively, which were then used to estimate the background potential density ρ and σ and stratification for θ θ Leg A (see below). For Leg B, two functional FP07 thermistors were installed on the VMP-250, so that the poten- tial density ρ and σ could be estimated, based on high-resolution FP07 thermistor data. To cross-calibrate the θ θ FP07 thermistors, the calibrated UCTD sensor and VMP-250 were clamped together and lowered to 200 m depth with the BT-winch of the Kagoshima-maru. To compute turbulent kinetic energy dissipation rates ε, the turbulent shear spectra were computed using the microscale shear data measured by the airfoil shear probe over 8 seconds, at approximately every 5–6 m depth, with an overlap of 4 s (2.5–3 m). The obtained shear spectra were integrated over the wavenumber ranges, where spectra agree with the Nasmyth empirical shear spectrum , avoiding the integrations of electronic noise at high wavenumber ranges and low frequency instrument vibration noise at low wavenumber ranges. For the wave- number ranges where the observed and the empirical spectra did not match, we integrated the fitted Nasmyth empirical shear spectrum (see Fig. 4). Assuming isotropic turbulence, ε was then computed by   15 ∂u   εν = ,     2 ∂z (5) where ν is the kinematic viscosity, and is the turbulent shear variance computed by the integration of the (/ ∂∂ uz ) shear spectrum. Because the UVMP does not carry a conductivity sensor to measure salinity, seawater potential density ρ and −3 σ = ρ − 1000 (kg m ) were estimated using the close relationship between potential temperature and potential θ θ density found at the observation sites (Supplementary Fig. S1). The relationship between potential temperature and potential density was derived from the UCTD data obtained across the front prior to the observations in Legs A and B (Fig. 1a), and from CTD data obtained at 36 stations during the surveys (Fig. 1a). The empirical relation- ship was derived by bin-averaging the observed potential density as a function of potential temperature every 0.5 °C (Supplementary Fig. S1). The potential density ρ and σ were then estimated by linear interpolation with θ θ the derived mean density model as a function of potential temperature. The standard deviation of the mean den- −3 −3 sity in each average bin was as small as <0.03 (kg m ) for σ > 25 (deeper than 150 m depth), 0.04–0.1 (kg m ) −3 for 24 < σ < 25 (100–150 m depth), and ~0.12 (kg m ) for σ < 24 (shallower than 100 m depth); the estimated θ θ density had better accuracy for the denser and deeper layers. Internal wave energy flux and wind power input. e in Th ternal wave energy flux along the vertical axis was estimated based on the wave horizontal kinetic energy HKE, and the vertical group velocity Cg , assuming that the near-inertial waves had low potential energy. Using the dispersion relation of the internal waves, the vertical group velocity for the low frequency internal wave is ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 11 www.nature.com/scientificreports/ f − ω ∣∣ Cg = , ωm (6) −1 −1 where ω is the internal wave frequency (rad s ), and m the internal-wave vertical wavenumber (rad m ). The horizontal kinetic energy of the internal waves is [] uv + zz HKE = ρ , o 2 (7) 2 m −1 where u and v are ADCP zonal and meridional shear (s ), respectively, and [] represents the spatiotempo- z z ral average between 100–250 m depth and over the entire transect for Leg A and Leg B. The internal-wave fre- −1 quency ω was estimated as ω = 1.01 f (rad s ), based on the internal ray path and the measured velocity shear band (Fig. 2) following a previous study . The range of the internal-wave vertical wavelength was estimated as 100–200 m with a mean of 150 m, based on the spectral peaks in the shear spectra (Fig. 3c,d). The corresponding −1 wavenumber range m = 0.063–0.031 and its mean m = 0.042 rad m were used to estimate the group velocity −3 Cg and HKE. e es Th timated horizontal HKE were 2.9–11.5 with a mean of 6.4 J m for Leg A and 6.3–25.2 with −3 −5 −1 a mean of 14.2 J m for Leg B. The range of the vertical group velocity was estimated to be 2.3–4.6 × 10 m s −5 −1 −2 with a mean of 3.5 × 10 m s . The resulting vertical internal wave energy flux for Leg A was 0.07–0.5 mW m −2 −2 −2 with a mean of 0.2 mW m , and that for Leg B was 0.1–1.2 mW m with a mean of 0.5 mW m (Supplementary Fig. S3b). To compare the internal wave energy flux with the wind power input to the near-inertial motions, a previous 48 58 slab model was used with the hourly reanalysis Japan Meteorological Agency GPV-MSM wind data from October 28 through November 30, 2016. Considering that the Kuroshio may have a strong advective ee ff ct and that remote wind-induced inertial motions could propagate to the observation sites, the hourly wind data were first averaged spatially over the region within 26–36°N and 125–131°E. The mixed layer depth was defined by the buoyancy frequency jump from the surface, and found to be 100 m depth where the large amplitude near-inertial shear was observed. Strong winds occurring before field observations on November 8–12, and during observa- tions on November 15 input energy to near-inertial motions (Supplementary Fig. S3). Analyses for the tidal elevation and long term moored ADCP data. Hourly tidal elevation data at Nakanoshima tidal station obtained from the Japan Oceanographic Data Center (JODC) were analyzed from August 31, 2000 through August 14, 2001. Harmonic analysis with 59 tidal constituents was conducted for the tidal elevation record to compare with the 30-h lowpass near-inertial shear amplitude (see below). The results of the harmonic analysis showed that the variance of the dominant three components for semi-diurnal (M , S , and 2 2 N ) and two components for diurnal (K and O ) tidal elevation can explain 99% of the total variance, with the M 2 1 1 2 tidal component accounting for 67%, S for 13%, N for 3%, K for 10%, and O for 6%. To obtain the time series of 2 2 1 1 subinertial diurnal and semi-diurnal tidal elevation by the harmonic analysis, the computed diurnal D = K + O 1 1 1 and semi-diurnal D = M + S + N tidal elevation amplitudes were smoothed with a 30-h lowpass filter and 2 2 2 2 compared with the 30-h lowpass near-inertial shear (Fig. 6). The same harmonic analysis was conducted for the tidal elevation data obtained during the shipboard observations in November 2016, at Tanegashima tidal station (Fig. 1a). The variance of total tidal elevation at Tanegashima station is contributed nearly identically from each dominant tidal constituent as in the Nakanoshima tidal data. To compare the obtained 30-h lowpass semi-diurnal and diurnal tidal amplitudes with the near-inertial shear amplitude, the relatively long-term moored ADCP data were analyzed. The moored ADCP (upward looking 75 kHz 20° beam angle Workhorse Long Ranger ADCP, Teledyne RDI) measured the horizontal current at 14 layers with 32 m vertical bins from 480.5 to 64.5 m depth from September 30, 2000 through July 1, 2001, at Stn. TK1 (Fig. 1a) in the vicinity of the Stn. M01 and M02. The hourly ADCP data were differentiated vertically to obtain the vertical shear. To extract the near-inertial shear, a bandpass Butterworth filter with cutoff frequencies at 0.8 f and 1.1 f (equivalent to the periods of 30 h and 21.7 h, respectively) was used. e Th near-inertial amplitude was computed using near-inertial zonal and meridional shear as uv + . The subinertial modulations of the zz near-inertial shear amplitude was then obtained using a lowpass Butterworth filter at 30 h for the uv + zz bandpass near-inertial shear amplitude. Data Availability. The datasets generated and/or analyzed during the current study are available from the corresponding author upon request. References 1. Qiu, B. Kuroshio and Oyashio Currents In: Encyclopedia of Ocean Sciences, Steele, J. H., Thorpe, (Oceanographer), S. A. & Turekian, Karl K. (Eds) (Academic Press, 2001). 2. Uda, M. Researches on “siome” or current rip in the seas and oceans. Geophys. Mag 11, 307–372 (1938). 3. Rudnick, D. L. Intensive surveys of the Azores Front. Part II: Inferring the geostrophic and vertical velocity fields. J. Geophys. Res. 101, 16,291–16,303 (1996). 4. Rudnick, D. L. & Luyten, J. R. Intensive surveys of the Azores Front part I. Tracers and dynamics. J. Geophys. Res. 101, 923–939 (1996). 5. Nurser, A. J. G. & Zhang, J. W. Eddy-induced mixed layer shallowing and mixed layer/thermocline exchange. J. Geophys. Res. 105, 21,851–21,868 (2000). 6. Guo, X., Zhu, X. H., Wu, Q. S. & Huang, D. The Kuroshio nutrient stream and its temporal variation in the East China Sea. Glob. Biogeochem. Cycles 117, C01026, https://doi.org/10.1029/2011JC007292 (2012). ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 12 www.nature.com/scientificreports/ 7. Guo, X., Zhu, X. H., Long, Y. & Huang, D. Spatial variations in the Kuroshio nutrient transport from the East China Sea to south of Japan. Biogeosciences 10, 6403–6417 (2013). 8. Pelegrí, J. L. & Csanady, G. T. Nutrient transport and mixing in the Gulf Stream. J. Geophys. Res. 96, 2577–2583 (1991). 9. Pelegrí, J. L., Csanady, G. T. & Martins, A. The North Atlantic nutrient stream. J. Ocean 52, 275–299 (1996). 10. Williams, R. G., Roussenov, V. & Follows, M. Induction of nutrients into the mixed layer and maintenance of high latitude productivity. Glob. Biogeochem. Cycles 20, GB1016, https://doi.org/10.1029/2005GB002586 (2006). 11. Williams, R. G. et al. Nutrient streams in the North Atlantic: Advective pathways of inorganic and dissolved organic nutrients. Glob. Biogeochem. Cycles 25, GB4008, https://doi.org/10.1029/2010GB003853 (2011). 12. Nagai, T., Aiba, M. & Clayton, S. Multiscale pathway to supply nutrients to the Kuroshio (In Japanese). Kaiyo to Seibutsu 37, 469–477 (2015). 13. Nagai, T. & Clayton, S. Nutrient interleaving below the mixed layer of the Kuroshio Extension Front. Ocean. Dyn. 67, 1027–1046, https://doi.org/10.1007/s10236-017-1070-3 (2017). 14. Rhines, P. B. & Young, W. R. Homogenization of potential vorticity in planetary gyres. J. Fluid Mech. 122, 347–367 (1982). 15. Oakey, N. S. & Elliot, J. A. Vertical temperature gradient structure across the Gulf Stream. J. Geophys. Res. 82, 1369–1380 (1977). 16. Gregg, M. C. & Sanford, T. B. Signatures of mixing from Bermuda slope, the Sargasso Sea and the Gulf Stream. J. Phys. Ocean. 10, 105–127 (1980). 17. Gargett, A. E. & Osborn, T. R. Small-scale shear measurements during the Fine and Microstructure Experiment (Fame). J. Geophys. Res. Ocean. 86, 1929–1944, https://doi.org/10.1029/JC086iC03p01929 (1981). 18. Hasegawa, D., Yamazaki, H., Lueck, R. G. & Seuront, L. How islands stir and fertilize the upper ocean. Geophys. Res. Lett. 31, n/a– n/a, https://doi.org/10.1029/2004GL020143 (2004). 19. Nikurashin, M. & Ferrari, R. Radiation and dissipation of internal waves generated by geostrophic motions impinging on small-scale topography: Theory. J. Phys. Ocean. 40, 1055–1074 (2010). 20. Tsutsumi, E. et al. Turbulent mixing within the Kuroshio in the Tokara Strait. J. Geophys. Res. Ocean. n/a–n/a, https://doi. org/10.1002/2017JC013049 (2017). 21. Lien, R.-C. & Gregg, M. C. Observations of turbulence in a tidal beam and across a coastal ridge. J. Geophys. Res. Ocean. 106, 4575–4591, https://doi.org/10.1029/2000JC000351 (2001). 22. Lueck, R. R. & Mudge, T. D. Topographically induced mixing around a shallow seamount. Sci. 276, 1831–1833, https://doi. org/10.1126/science.276.5320.1831 (1997). 23. Rainville, L. & Pinkel, R. Observations of energetic high-wavenumber internal waves in the Kuroshio. J. Phys. Oceanogr. 34, 1495–1505 (2004). 24. Plougonven, R. & Snyder, C. Inertia-gravity waves spontaneously generated by jets and fronts. Part I: Different baroclinic life cycles. J. Atmos. Sci. 64, 2502–2520 (2007). 25. Alford, M. H., Shcherbina, A. Y. & Gregg, M. C. Observations of near-inertial gravity waves radiating from a frontal jet. J. Phys. Oceanogr. 43, 1225–1239 (2013). 26. Shakespeare, C. J. & Taylor, J. R. The spontaneous generation of inertia-gravity waves during frontogenesis forced by large strain: theory. J. Fluid. Mech. 757, 817–853 (2014). 27. Nagai, T., Tandon, A., Kunze, E. & Mahadevan, A. Spontaneous generation of near-inertial waves by the Kuroshio Front. J. Phys. Oceanogr. 45, 2381–2406, https://doi.org/10.1175/jpo-d-14-0086.1 (2015a). 28. Shakespeare, C. J. & Hogg, A. M. Spontaneous surface generation and interior amplification of internal waves in a regional-scale ocean model. J. Phys. Oceanogr. 47, 811–826, https://doi.org/10.1175/JPO-D-16-0188.1 (2017). 29. Hibiya, T., Niwa, Y. & Fujiwara, K. Numerical experiments of nonlinear energy transfer within the oceanic internal wave spectrum. J. Geophys. Res. Ocean. 103, 18715–18722, https://doi.org/10.1029/98JC01362 (1998). 30. MacKinnon, J. A. & Winters, K. B. Subtropical catastrophe: Significant loss of low-mode tidal energy at 28.9°. Geophys. Res. Lett. 32, n/a–n/a, https://doi.org/10.1029/2005GL023376 (2005). 31. Thomas, L. N. & Taylor, J. R. Damping of inertial motions by parametric subharmonic instability in baroclinic currents. J. Fluid Mech. 743, 280–294, https://doi.org/10.1017/jfm.2014.29 (2014). 32. Mooers, C. N. Several effects of baroclinic current on the cross-stream propagation of inertial-internal waves. Geophys. Fluid. Dyn. 6, 245–275 (1975). 33. Kunze, E. Near-inertial wave propagation in geostrophic shear. J. Phys. Oceanogr. 15, 544–565 (1985). 34. Whitt, D. & Thomas, L. Near-inertial waves in strongly baroclinic currents. J. Phys. Oceanogr. 43, 706–725 (2013). 35. Kunze, E. & Sanford, B. Observations of near-inertial waves in a front. J. Phys. Oceanogr. 15, 566–581 (1984). 36. Kunze, E. The mean and near-inertial velocity fields in a warm-core ring. J. Phys. Oceanogr. 16, 1444–1461 (1986). 37. Lueck, R. G. & Osborn, T. The dissipation of kinetic energy in a warm-core ring. J. Geophys. Res. 91, 803–818 (1986). 38. Winkel, D. P., Gregg, M. C. & Sanford, T. B. Patterns of shear and turbulence across the Florida Current. J. Phys Oceanogr. 32, 3269–3285 (2002). 39. Nagai, T., Tandon, A., Yamazaki, H. & Doubell, M. J. Evidence of enhanced turbulent dissipation in the frontogenetic Kuroshio Front thermocline. Geophys. Res. Lett. 36, L12609, https://doi.org/10.1029/2009GL038832 (2009). 40. Nagai, T., Tandon, A., Yamazaki, H., Doubell, M. J. & Gallager, S. Direct observations of microscale turbulence and thermohaline structure in the Kuroshio Front. J. Geophys. Res. 117, C08013, https://doi.org/10.1029/2011JC00722 (2012). 41. Kaneko, H., Yasuda, I., Komatsu, K. & Itoh, S. Observations of the structure of turbulent mixing across the Kuroshio. Geophys. Res. Lett. 39, n/a–n/a, https://doi.org/10.1029/2012GL052419 (2012). 42. Kaneko, H., Yasuda, I., Komatsu, K. & Itoh, S. Observations of vertical turbulent nitrate flux across the Kuroshio. Geophys. Res. Lett. 40, 3123–3127, https://doi.org/10.1002/grl.50613 (2013). 43. Inoue, R., Gregg, M. C. & Harcourt, R. R. Mixing rates across the Gulf Stream, Part 1: On the formation of Eighteen Degree Water. J. Mar. Res. 68, 643–671 (2010). 44. Alford, M. H. & Gregg, M. C. Near-inertial mixing: Modulation of shear, strain and microstructure at low latitude. J. Geophys. Res. Ocean. 106, 16947–16968, https://doi.org/10.1029/2000JC000370 (2001). 45. Alford, M. H., MacKinnon, J. A., Simmons, H. L. & Nash, J. D. Near-inertial internal gravity waves in the ocean. Annu. Rev. Mar. Sci. 8, 95–123, https://doi.org/10.1146/annurev-marine-010814-015746, PMID: 26331898 (2016). 46. Garrett, C. J. R. & Munk, W. H. Space-time scales of internal waves. Geophys. Fluid Dyn. 2, 225–264 (1972). 47. Garrett, C. J. R. & Munk, W. H. Space-time scales of internal waves. A progress report. J. Geophys. Res. 80, 291–323 (1975). 48. Pollard, R. T. & Millard, R. C. Comparison between observed and simulated wind-generated inertial oscillations. Deep. Res. 17, 813–821 (1970). 49. D’Asaro, E. Wind-forced internal waves in the North Pacific and Sargasso Sea. J. Phys. Ocean. 14, 781–794 (1984). 50. Alford, M. H. Redistribution of energy available for ocean mixing by long-range propagation of internal waves. Nat. 423, 159–162 (2003). 51. Nasmyth, P. W. Oceanic turbulence. Ph.D. thesis, University of British Columbia (1970). 52. Kunze, E., Firing, E., Hummon, J. M., Chereskin, T. K. & Thurnherr, A. M. Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles. J. Phys. Oceanogr. 36, 1553–1576, https://doi.org/10.1175/JPO2926.1 (2006). 53. Whalen, C. B., Talley, L. D. & MacKinnon, J. A. Spatial and temporal variability of global ocean mixing inferred from Argo profiles. Geophys. Res. Lett. 39, n/a–n/a, https://doi.org/10.1029/2012GL053196 (2012). ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 13 www.nature.com/scientificreports/ 54. Osborn, T. Estimates of the local rate of vertical diffusion from dissipation measurements. J. Phys. Oceanogr. 10, 83–89 (1980). 55. Xie, X.-H., Shang, X.-D., van Haren, H., Chen, G.-Y. & Zhang, Y.-Z. Observations of parametric subharmonic instability-induced near-inertial waves equatorward of the critical diurnal latitude. G e o p hy s . R es . L et t . 38, n/a–n/a, ht t ps://doi . org/10.1029/2010GL046521 (2011). 56. Liu, T., Wang, B., Hirose, N., Yamashiro, T. & Yamada, H. High-resolution modeling of the Kuroshio Current power south of Japan. Submitted to Journal of Ocean Engineering and Marine Energy (2017). 57. Nakamura, H. et al. Intermediate and deep water formation in the Okinawa Trough. J. Geophys. Res. Ocean. 118, 6881–6893, https:// doi.org/10.1002/2013JC009326 (2013). 58. Saito, K. et al. The operational JMA nonhydrostatic mesoscale model. Mon. Weather. Rev. 134, 1266–1298 (2006). Acknowledgements We thank Capt. Uchiyama, C/O Azuma, Prof. Nishina, Prof. Kobari, Prof. Matsuno, Prof. Senju, Dr. Li, crews, and all the participants of the R/T/V Kagoshima-maru cruise KG1616. T.N. thanks Hancyk, Dr. Wolk at RSI, Dr. Li, Seki, and Yoshii at JAC, Dr. Kokubu at NORTEK, Rutka at STS, support from OMIX (MEXT KAKENHI JP16H01590) and SKED (“Study of Kuroshio Ecosystem Dynamics for Sustainable Fisheries” funded to FRA from MEXT). D.H., T.T., and R.I. thank OMIX (MEXT KAKENHI) JP15H05818. H.N., and E.T. thank OMIX (MEXT KAKENHI) JP15H05821. Author Contributions T.N. conceived the tow-yo turbulence experiments, D.H., T.N., T.T. materialized the tow-yo microstructure profiling, H.N. organized the cruise, E.T. analysed moored ADCP data and conducted harmonic analyses, the authors (D.H., T.T., H.N., E.T.) conducted the shipboard experiments, T.Y. conducted the long-term mooring observation, R.I. provided a slab model for wind analysis, and T.N. analysed the results. All authors reviewed and edited the manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-15167-1. Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 14 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

First Evidence of Coherent Bands of Strong Turbulent Layers Associated with High-Wavenumber Internal-Wave Shear in the Upstream Kuroshio

Free
14 pages
Loading next page...
 
/lp/springer_journal/first-evidence-of-coherent-bands-of-strong-turbulent-layers-associated-kt83mX16QG
Publisher
Nature Publishing Group UK
Copyright
Copyright © 2017 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
D.O.I.
10.1038/s41598-017-15167-1
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports Correction: Author Correction OPEN First Evidence of Coherent Bands of Strong Turbulent Layers Associated with High-Wavenumber Internal- Received: 5 May 2017 Wave Shear in the Upstream Accepted: 23 October 2017 Published online: 06 November 2017 Kuroshio 1 2 2 3 Takeyoshi Nagai , Daisuke Hasegawa , Takahiro Tanaka , Hirohiko Nakamura , Eisuke 4 5 6 Tsutsumi , Ryuichiro Inoue & Toru Yamashiro The upstream Kuroshio flows through Okinawa Trough and the Tokara island chain, the region near the continental shelf of the East China Sea and shallow seamounts, where the Kuroshio can induce strong mixing over the shallow topography. Also, tidal currents over the rough topography may produce internal tides, and associated turbulence. The previous observations show energetic high vertical wavenumber near-inertial wave shear in the Kuroshio thermocline, which implies strong turbulent mixing. However, direct turbulence measurements in this region are very scarce. Using high lateral resolution (1–2 km) direct turbulence measurements, we show here, for the first time, that strong turbulent layers form spatially coherent banded structures with lateral scales of >O(10 km), associated with bands of near-inertial wave/diurnal internal tide shear of high vertical wavenumber in the upstream Kuroshio. The turbulent kinetic energy dissipation rates within these turbulent layers −7 −1 −4 2 −1 are >O(10 W kg ), and estimated vertical eddy diffusivity shows >O(10 m s ) on average. These results suggest that the high vertical wavenumber near-inertial waves propagating in the upstream Kuroshio could have large impacts on the watermass modifications, momentum mixing, nutrient supply, and associated biogeochemical responses in its downstream. The Kuroshio, the western boundary current of the subtropical gyer in the North Pacific, has been known to transport large amounts of heat and salt from the tropical ocean . Because the Kuroshio is inherently an ocean dynamic front, where the associated jet flows along the interface between seawaters of different water properties , mixing processes in the Kuroshio Front play very important roles in watermass formation, transformation, and 3–5 subduction . In addition to this physical or hydrographical important role, the mixing in the Kuroshio has profound implications for biogeochemistry because the Kuroshio carries waters of relatively high nutrient con- 6,7 centrations from the tropics through its subsurface layers, as known as the nutrient stream similar to the Gulf 8–11 8,9 12,13 Stream . Earlier studies in the Gulf Stream and relatively recent studies in the Kuroshio have reported that elevated nutrient concentrations are found along these western boundary currents in their subsurface layers compared to that in the ambient waters of the same density. Because along-isopycnal stirring by mesoscale eddies tends to homogenize any tracer anomaly on the density surface in the ocean interior , the elevated nutrient concentrations at the subsurface of the Kuroshio are most likely formed by across-density/diapycnal fluxes by Tokyo University of Marine Science and Technology, Department of Ocean Sciences, Tokyo, 108-8477, Japan. Tohoku National Fisheries Research Institute, Japan Fisheries Research and Education Agency, Fisheries Oceanography and Resources Department, Shiogama, Miyagi, 985-0001, Japan. Kagoshima University, Faculty of Fisheries, Kagoshima, 890-0056, Japan. Kyushu University, Research Institute for Applied Mechanics, Kasuga, Fukuoka, 816-8580, Japan. Research and Development Center for Global Change, Japan Agency for Marine-Earth Science and Technology, Yokosuka, 237-0061, Japan. Graduate School of Science and Engineering, Kagoshima University, Kagoshima, 890-0065, Japan. Daisuke Hasegawa, Takahiro Tanaka, Hirohiko Nakamura, Eisuke Tsutsumi, Ryuichiro Inoue and Toru Yamashiro contributed equally to this work. Correspondence and requests for materials should be addressed to T.N. (email: tnagai@kaiyodai.ac.jp) ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 1 www.nature.com/scientificreports/ irreversible diabatic turbulent diffusion. A previous study speculated that the turbulent diffusive nutrient flux by geostrophic vertical shear contributes to the concentration of elevated nutrients below the Gulf Stream. However, a limited number of direct turbulence observations at that time did not provide the supporting evidence for this 15–17 hypothesis . On the other hand, the Kuroshio flows through shallow topography in its upstream along Okinawa Trough and near Tokara Strait, where the Kuroshio over the topography can induce strong turbulence and diapycnal mix- 18–20 ing . While a number of previous studies have reported that the beams of the M internal tide emanated from 21,22 the rough topography induced strong turbulence in close proximity to canyons and seamounts , recent ship- board and lowered Acoustic Doppler Current Profiler (ADCP) measurements in the upstream Kuroshio showed widespread large amplitude near-inertial internal wave shear of high vertical wavenumber in the Kuroshio and regions between continental margins and the Kuroshio . A very recent field study in Tokara Strait reported near-inertial waves of even shorter vertical wavelengths ~50 m accompanied by pronounced turbulent dissipa- −7 −6 −1 20 tion rates of O(10 –10 W kg ) . One of the most likely mechanisms to form the high vertical wavenumber shear is the near-inertial internal waves generated by the Kuroshio over shallow topography in Okinawa Trough and Tokara Strait . Because the Kuroshio is forced to turn its direction from northeast to southeast when it approaches Kyushu island, high vertical wavenumber near-inertial internal waves could also be generated spon- 24–28 taneously by the meandering Kuroshio . Furthermore, the upstream Kuroshio passes the so-called critical latitude, 28.9°N, where M internal tide energy can be efficiently converted to that of high vertical wavenumber 29,30 near-inertial waves through parametric subharmonic instability (hereinae ft r PSI) . Near the critical latitude, K or O diurnal internal tides generated over the rough topography are also near-inertial internal waves. In a 1 1 strong baroclinic front, where geostrophic shear allows anomalously low frequency internal waves, a recent the- oretical study showed that PSI of wind-driven near-inertial waves could generate high wavenumber internal waves with a frequency of half the local Coriolis frequency f/2. Regardless of their generation mechanism, these propagating near-inertial internal waves can be trapped on the anticyclonic side of the front where the lowest 32–34 internal wave frequency is decreased by geostrophic lateral and vertical shears . The trapped near-inertial 35–37 waves can break into three-dimensional turbulence and enhance mixing of momentum and tracers . However, direct turbulence measurements in the upstream Kuroshio are very scarce. Also, the lateral resolutions of previous turbulence measurements are not sufficient to resolve the coherent spatial patterns associated with high wave- 15–17,38 number near-inertial internal wave shear . As a result, earlier studies did not find elevated turbulence in 15–17 the thermocline of the Gulf Stream , and concluded that turbulence in the ocean interior is caused by patchy random internal wave breaking, even in the Gulf Stream. In contrast, a number of recent studies have reported 27,39–42 43 enhanced turbulent mixing in the thermocline of the Kuroshio and the Gulf Stream caused by propagating internal waves. However, the lateral resolutions of the turbulence data were still insufficient to provide direct evi- dence that the turbulence is caused coherently by banded internal-wave shear in the thermocline of the western 44,45 boundary currents. Although a recent extensive microstructure time-series survey at low latitude reported that the strong turbulence coincided with the relatively large strain caused by wind-induced near-inertial waves of anomalously low frequency in the anticyclonic vorticity field , it is not known if near-inertial internal waves form the spatially coherent strong turbulent layers in ocean fronts. In this study, a new turbulence and microstructure observation scheme is developed and utilized to measure turbulence directly with a high lateral resolution of 1–2 km in the upstream Kuroshio near Tokara Strait, where a large amount of energy from the Kuroshio can be dissipated. For the first time, the observations in this study suggest that propagating high vertical wavenumber near-inertial internal waves form spatially coherent banded layers of strong turbulence in the upstream Kuroshio, with O(100 m) and >O(10 km) vertical and lateral scales, respectively. Results Observations. e Th in-situ survey was conducted using the R/T/V Kagoshima-maru during November 12–20, 2016, near Tokara Strait where the upstream Kuroshio flows over the shallow seamounts and through the narrow straits near the Tokara island chain off Kyushu island (Fig.  1). Two mooring systems, which included an upward looking ADCP (75 kHz 20° beam angle Workhorse Long Ranger ADCP, Teledyne RDI) were deployed near the Kuroshio and measured horizontal currents and directions during November 13–19, 2016, at Stn. M01 (30.35°N, 129.85°E) and M02 (30.00°N, 129.50°E, Fig. 1a). Two transect observations were carried out to obtain current velocity data using a shipboard ADCP (75 kHz 30° beam angle, Ocean Surveyor, Teledyne RDI) and high-resolu- tion towed microstructure data. The transect observations (hereinaer L ft eg A) were conducted on November 14 from west to east closely along the zonal direction for the first half of the transect, and continued in the northeast direction for the rest (thick blue line in Fig. 1b). The transect surveys were repeated (hereinafter Leg B; red in Fig. 1b) on November 18–19 close to the same ship track as for the second half of the Leg A observations. During these two transect observations, a new tow-yo Underway Vertical Microstructure Profiler (Underway-VMP, here- inaer UVMP ft ; Rockland Scientific International, Victoria, Canada) was used in the Kuroshio Front to investigate high-resolution spatial structures of microscale turbulence. Prior to the transect observations of Legs A and B, the Underway Conductivity Temperature Depth profiler (hereinaer UCTD; T ft eledyne Oceanscience, California, USA) was deployed across front direction repeatedly (thick black line in Fig. 1a). Conductivity-Temperature- Depth (CTD) profiles also were made at 36 stations (Fig.  1a). Because the UVMP used in this study was a pro- totype with no conductivity sensor, simultaneously gathered salinity data are not available. However, because water properties obtained by these UCTD and CTD in the observed regions suggest a close relationship between temperature and density (Supplementary Fig. S1), we estimated density in Legs A and B from temperature alone (see Methods section). In addition to the mooring observations made during the shipboard surveys, the relatively long-term moored ADCP records from September 30, 2000 through July 1, 2001 obtained at Stn. TK1 (30.13°N, 130.19°E, Fig. 1a) in the vicinity of the shipboard observation sites were analyzed to investigate the mechanisms generating observed internal waves. ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 2 www.nature.com/scientificreports/ Figure 1. Observation site in Tokara Strait near the Tokara island chain off Kyushu island. (a) Underway-CTD (UCTD) observation line is shown as a thick black line, and Underway-VMP (UVMP) observation lines are in blue for Leg A and red for Leg B. Mooring stations Stn. M01 and M02 are indicated by black (+). CTD stations are denoted by solid white circles (). Total number of the CTD profiles is 36, which consists of 13 profiles each at Stn. M01 and M02 during time series observations and 10 profiles at four other stations near the Kuroshio. e r Th elatively long-term mooring Stn. TK1 is shown as a black triangle. Nakanoshima, and Tanegashima tidal stations are shown as a red diamond, and a yellow square, respectively. Black contours are AVISO sea surface height (m) averaged from November 14–19, 2016. (b) Enlarged map of the UVMP observation site. UVMP in Leg A is denoted by blue diamonds and Leg B by red triangles. Black crosses (×) indicate six XBT stations during Leg A. Color shadings indicate depth (m). The map is created with MATLB R2016b with ETOPO1 topography data available from NGDC NOAA https://www.ngdc.noaa.gov/mgg/global/. High-vertical-wavenumber internal-wave shear. Measured horizontal current velocity using the ship- board ADCP mostly shows a northeastward flow associated with the Kuroshio, consistent with the sea surface height (Fig. 1). During the Leg A observation, this flow was strongest in the westernmost edge of the transect, where the main stream of the Kuroshio flowed (Fig.  2a). On the other hand, in Leg B, the largest velocity magni- tude was recorded at the northeastern part, the last half of the transect (Fig. 2d), as the continuous UVMP profil- ing in Leg B started from the southwestern region to the northeastern direction, reaching the main stream of the Kuroshio Front in the last quarter of the Leg B (Fig. 1). In both the Legs, current velocity magnitudes are largest near the surface, due to the baroclinic flow of the Kuroshio, with some exceptions where the stronger currents are found subsurface, especially on the anticyclonic side of the Kuroshio (Fig. 2a,d). Back-rotated velocity shear, which is the vertical gradient of the lateral velocity rotated horizontally, assuming the inertial rotation with respect to the reference time (see Methods section), shows bands of alternating positive and negative signs, close to the density surface (Fig. 2b,c,e,f ), which are similar to a study by Rainville and Pinkel . These banded shears are typical characteristics of the propagating inertia-gravity waves in a stratified ocean. The wavelengths of these banded shear structures are O(100 m) in vertical and O(10 km) in lateral direction, similar to Rainville and Pinkel’s study in this region . With λ and λ , horizontal and vertical wavelengths (m), respectively, Nb =∂ /∂z h v −1 −2 −2 buoyancy frequency (rad s ) (where b = −gρ/ρ is buoyancy (m s ) with gravitational acceleration g (m s ), −3 −3 seawater density ρ (kg m ) and reference density ρ = 1025 (kg m ), and vertical coordinate z(m)), and f Coriolis −1 frequency (rad s ), the internal-wave dispersion relation, then implies that the frequency of these waves 22 22 2 −2 −1 −4 −1 ω =+ (/ fN λλ /) λ ~ f , based on the estimated N~O(10 rad s ), and local f~O(10 rad s ). The zh z internal-wave ray paths computed, assuming a quiescent condition show that the angles of the internal wave rays with near-inertial frequencies are more consistent with the observed shear bands than those with the M tidal frequency (Fig.  2b,c,e,f ). This also suggests that the observed high-wavenumber velocity shear is due to near-inertial internal waves. The total shear vertical wavenumber spectra averaged for the entire transect of Leg A and B suggest that shear variance is significantly larger than that of Garrett-Munk (GM) internal-wave equilib- 46,47 rium spectra by several factors, at vertical wavelengths larger than 100 m (Fig. 3c,d), at which the resolution of the shipboard ADCP starts to fail to resolve shear of shorter wavelengths. Computed Richardson numbers 22 2 RiN =+ /(uv ), where u and v are zonal and meridional shear, show relatively small values in the banded zz z z shear layers and near the surface. However, because of the coarse resolution of the ADCP shear, the Ri are mostly larger than the critical value for the Kelvin-Helmholtz instability Ri = 0.25 except near the surface and below 500 m in Leg A and 300 m in Leg B (Supplementary Fig. S2). e Th rotary shear spectra in Leg A show that the var - iance of clockwise rotating shear with depth is larger than that of anticlockwise rotating component (Fig. 3c), suggesting the dominant downward energy propagating internal waves. In Leg B, the clockwise and anticlockwise components of rotary shear spectra are comparable with slightly larger shear variance for the anticlockwise ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 3 www.nature.com/scientificreports/ Figure 2. Vertical sections of the shipboard Acoustic Doppler Current Profiler velocity measurements. (a–c) are for Leg A (blue in Fig. 1b) and (d–f) are for Leg B (red in Fig. 1b). Horizontal absolute current velocity −1 −1 (m s ) are shown in (a,d). Back-rotated shear [s ] is shown for zonal shear u (t ) (b,e), and for meridional z 0 −3 shears v (t ) (c,f) (see Methods section). Black contours are σ (kg m ). Magenta curves in (b,c,e,f) are internal- z 0 θ wave ray paths at frequencies of 1.01 f, 1.1 f, 1.2 f, 1.4 f, and M tidal frequency, where f is Coriolis frequency, as indicated in the panels. Red triangles in upper abscissa are the locations at which hodographs of shear are plotted for each Leg in Fig. 3. rotating component (Fig. 3d). The ratio of the integrated anticlockwise shear spectrum to that for the clockwise uu u u z z zz component (where and are anticlockwise and clockwise components of shear spec- φφ dk/ dk φ φ ∫∫ ACWCW ACW CW trum, respectively, and k is the vertical wavenumber) as a function of longitude suggests that large shear variance found at 129.3–129.575°E in Leg A is associated with the downward energy propagating internal waves, and that found at 129.4–129.55°E in Leg B is caused by upward energy propagating internal waves (Fig. 3e,f ). The hodo- graph of the shear at 129.447°E in Leg A, where shear magnitude is large, shows that shear is rotating clockwise with depth, suggesting downward energy propagating internal waves, while that at 129.472°E in Leg B shows anticlockwise rotation for upward energy propagating internal waves, consistent with the trends found in the −3 −3 shear spectra (Fig. 3a,b,e,f ). The estimated average internal wave energy is 6.4 J m for Leg A, and 14.2 J m for −5 −1 Leg B. The estimated mean group velocity, 3.5 × 10 m s with the mean vertical wavelength of the internal −2 −2 wave, 150 m, results in the mean vertical energy flux of 0.2 mW m for Leg A, and 0.5 mW m for Leg B (see Methods section). These estimates are comparable to the rate of wind energy input to inertial motions computed 48–50 by a slab model (see Methods section, Supplementary Fig. S3). es Th e results indicate that the high vertical wavenumber shear is caused by near-inertial internal waves prop- agating in the upward and downward directions in the Kuroshio near Tokara Strait. While the wind energy flux to the inertial currents is found to be sufficient to generate an observed amplitude of near-inertial waves, it is not clear how the observed high vertical wavenumber near-inertial internal waves are formed. We speculate that the possible mechanisms are (1) PSI of the M internal tide upstream, which was advected to the observation sites by the Kuroshio flow, (2) local PSI with anomalously lower minimum internal-wave frequency than f due to mean velocity shears, (3) K internal tides generated over rough topography, (4) near-inertial internal waves generated ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 4 www.nature.com/scientificreports/ −1 Figure 3. Hodographs of the back-rotated shear [s ] are shown for (a) 129.447°E in Leg A (Fig. 2a–c) and for (b) 129.472°E in Leg B (Fig. 2d–f ). Shading indicates depth (m). Average vertical wavenumber shear spectra −2 −1 [s cp m ] are shown for (c) Leg A and for (d) Leg B. Solid black lines (CW) are the clockwise rotating component with depth and dashed lines (ACW) are the anticlockwise rotating component. Solid black-white- black lines are for total shear spectra. Dotted lines (GM) are Garrett-Munk equilibrium internal wave shear spectra. (e,f) Ratio of integrated anticlockwise shear spectra to that for the clockwise component uu u u zz z z φφ dk/ dk (where φ and φ are anticlockwise and clockwise components of shear spectrum, ∫∫ ACWCW ACW CW respectively and k is the vertical wavenumber) as a function of longitude in Leg A (e) and B (f). ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 5 www.nature.com/scientificreports/ Figure 4. Microscale shear and shear spectra obtained from tow-yo observations. (a) One of the turbulent −1 −2 −1 shear profile [s ] and (b,c) turbulent shear spectra [s cp m ]. The shear spectra in red and blue correspond to the depth ranges indicated by the same color in (a). The Nasmyth empirical shear spectra are shown in (b,c) as −1 black curves for each computed turbulent kinetic energy dissipation rate ε [W kg ]. by the Kuroshio over shallow topography in Okinawa Trough and Tokara Strait, (5) spontaneously generated near inertial internal waves by the Kuroshio, and (6) PSI of wind-induced inertial waves to generate high wavenumber anomalously low frequency internal waves of half the Coriolis frequency, f/2 along isopycnal in a strong baroclinic front . High-resolution tow-yo turbulence measurements. The shipboard ADCP measurement shows the large amplitude banded high vertical wavenumber shear, similar to previous studies . Here, we show, for the first time, that these large amplitude near-inertial internal wave shears generate strong turbulence of spatially coherent banded structures, employing a newly developed high-resolution (1–2 km) tow-yo microstructure profiling tech- nique, using 1.5 mm dyneema rope. The thin rope allows the microstructure profiler (VMP-250) to sink relatively smoothly with small drag down to 300 m depth. The measured turbulent shear spectra by the tow-yo profiling show overall good agreement with the Nasmyth empirical shear spectrum (Fig. 4). This novel technique, there- fore, allows us to resolve spatial structures of intermittent turbulence, which have been unrecognized because of the coarse resolution measurements by the previous microstructure profiling. −7 −1 In the upper 100 m, the measured turbulent kinetic energy dissipation rate, ε is large, O(10 W kg ), because of turbulence within the surface boundary layer, which is forced by momentum and heat fluxes at the surface (Fig. 5). In the stratified layers below 100 m depth, away from the surface boundary layer, the dissipation rates −7 −1 −6 −1 still frequently exceed O(10 W kg ), and reach O(10 W kg ) at some locations (Fig. 5a,b). These dissipation rates are 100 to1000 times greater than those found typically in the open water thermocline. More importantly, these large dissipation rates form coherent along-isopycnal bands of strong turbulent layers. These strong turbu- lent layers found beneath and on the less dense side of the Kuroshio correspond very well to the large amplitude near-inertial wave shear bands, suggesting that the observed near-inertial internal waves of high vertical wave- number break into microscale turbulence. The comparison between our observed turbulent dissipation rates and 52,53 previous internal-wave parameterization for energy dissipation rates show a positive correlation (r = 0.63, p-value = 0.0001) within the same orders of magnitude (Supplementary Fig. S4). For the first time, our high ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 6 www.nature.com/scientificreports/ −1 Figure 5. Vertical section of (a,b) measured turbulent kinetic energy dissipation rates log ε [W kg ], and 2 −1 (c,d) estimated vertical turbulent eddy diffusivities log K [m s ] for (a,c) Leg A and (b,d) Leg B. Contours are σ . resolution microstructure surveys provide evidence that strong turbulent layers with spacial scales of O(100 m) and O(10 km) in the vertical and lateral directions, respectively, are coherently formed along the near-inertial internal-wave shear bands (Figs 2 and 5a,b). To investigate the significance of this turbulence not only for energy dissipation but also for the vertical mixing and diffusion, the vertical turbulent eddy diffusivity is estimated using measured dissipation rates ε and buoyancy frequency N with the previous diffusivity model , namely, K = γ , (1) 54 −2 2 −1 where γ = 0.2 is a mixing efficiency factor . The estimated vertical turbulent eddy diffusivities are > O(10 m s ) within the surface boundary layer, which is not surprising considering continuous forcing at the surface with weak stratification in the surface boundary layer (Fig.  5c,d). The vertical eddy diffusivities in the stratified subsur - −4 2 −1 −3 2 −1 face layers, however, are notably high, O(10 m s ) on average, and >O(10 m s ) at many locations below 100 m depth (Fig. 5c,d). The high diffusivity layers are structured in a banded form, reflecting the spatial struc- tures found in the dissipation rates. Discussion Early microstructure surveys in the 1970s–1980s did not find strong turbulence across the western boundary 15–17 current, i.e., the Gulf Stream, but found intermittent turbulent patches due to random internal wave breaking . Relatively recent microstructure measurements in the upstream Gulf Stream, Florida Current , and in the down- stream Gulf Stream , reported the possibility of the internal wave induced turbulence, but whether turbulent layer is formed coherently with the banded internal wave shear is unclear due to the coarse lateral resolution of the surveys. In contrast, this study found coherent bands of very strong turbulence clearly associated with the large amplitude, high vertical wavenumber internal-wave shear, in the upstream Kuroshio. ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 7 www.nature.com/scientificreports/ e o Th bserved strong turbulence associated with high vertical wavenumber internal wave shear implies active occurrence of PSI near the critical latitude 28.9°N. The moored ADCP observations during our shipboard exper - iment suggest that the velocity field is dominantly modulated at M semi-diurnal tidal frequency (Supplementary Fig. S5a,b). In contrast, the moored ADCP vertical shear modulates at subinertial, near-inertial (or K diur- nal), as well as M semi-diurnal tidal frequencies. However, these M tidal, near-inertial and subinertial peaks 2 2 in shear spectra are not clear, especially in upper layers, probably because of the Doppler smearing by the strong Kuroshio flow (Supplementary Fig. S5). The energy of the observed M semi-diurnal internal tide could be converted to that of high vertical wavenumber near-inertial waves through PSI. To investigate the possibility of PSI and other mechanisms to generate observed high vertical wavenumber internal wave shear, the long-term moored ADCP data obtained at Stn. TK1 from September 2000 through July 2001 were analyzed. Assuming the energy source of the PSI-induced near-inertial waves is M internal tide, the amplitude of the PSI-induced near-inertial wave shear should modulate fortnightly with spring-neap tidal cycle , similarly to that of M inter- nal tides. The harmonic analysis of the tidal elevation at Nakanoshima (29.84°N, 129.85°E, Fig.  1a), suggests that the diurnal tidal amplitude (D = K + O ) is about the half of that estimated for the semi-diurnal com- 1 1 1 ponent (D = M + S + N ) (Supplementary Fig. S6b). The 30-h lowpass near-inertial shear amplitude exhibits 2 2 2 2 subinertial variations with larger amplitude found in the upper 250 m (Supplementary Fig. S6c). Most of these subinertial variations are found to occur fortnightly, with a period of 14 days, throughout the water column (Fig. 6a,b). The fortnightly modulation found in the near-inertial shear amplitude may imply occurrence of the PSI . However, our observation sites are very close to the critical latitude, where the frequency of diurnal tidal flow is at the near-inertial frequency. It is also possible that diurnal near-inertial internal tides are generated directly from the nearby rough topography, which then propagate to the observation sites. To elucidate whether fortnight modulations of near-inertial shear amplitude are correlated more with those of semi-diurnal D or diurnal D tidal amplitude, the correlation analyses between D or D tidal amplitude and near-inertial shear 1 2 1 amplitude are conducted (Supplementary Fig. S7). The correlation between D diurnal tidal amplitude and the 30-h lowpass near-inertial shear amplitude is significantly high with a positive Spearman correlation coefficient of 0.59. On the other hand, the correlation with D semi-diurnal tidal amplitude is low, at 0.05. The results suggest that although the near-inertial shear amplitude modulates at the fortnight frequency, the dominant source of the near-inertial energy is D diurnal tides, not D semi-diurnal internal tides through PSI. Our study estimates 1 2 −2 10–100 mW m for the K internal tide energy generation rate at the bottom near Tokara Strait, based on a regional numerical simulation which could be sufficiently large to generate an observed vertical energy ux fl of −2 the near-inertial internal waves, 0.1–1 mW m . Accordingly, the local PSI of M internal tides, with the anoma- lously low minimum internal wave frequency due to geostrophic shears, is unlikely. However, if the PSI-induced near-inertial internal waves are advected by the Kuroshio from the upstream region south of the critical latitude, the fortnightly modulations of the Doppler-smeared near-inertial shear amplitude should be lagged behind that of the D semi-diurnal tide. The lagged Spearman correlation coefficients between the D semi-diurnal tidal 2 2 amplitude and the 30-h lowpass near-inertial shear amplitude show a maximum correlation coefficient of 0.34 with a 94 h time lag (Supplementary Fig. S7). With a distance from the critical latitude region to the observa- −1 tion sites along the Kuroshio, ~160 km, the time lag of 94 h implies an average advection speed at 0.47 m s . The estimated advection speed is consistent with the velocity magnitude observed in the subsurface layers of −1 the Kuroshio, ~0.5 m s between 200–400 m depth (Fig. 2a,d). Assuming that the estimated near-inertial wave −3 energy, O(10 J m ) (see Methods section) is continuously dissipated by the observed bands of the turbulent −8 −7 −1 −3 dissipation rates, O(10 –10 W kg ), with the reference seawater density ρ = 1025 kg m , the internal waves could be dissipated completely aer ft 1–10 days. er Th efore, it is possible that the PSI-induced near-inertial waves can survive ~4 days after the 94-h advection. Even without the tide or PSI, near-inertial internal waves of high vertical wavenumber could be generated by the Kuroshio over the topography . The globally estimated energy generation rate from geostrophic flow to lee wave by the study of Nikurashin and Ferrari in Okinawa Trough, −2 the region near the observation sites, is O(0.1–1 mW m ), consistent with the vertical energy flux of observed near-inertial internal waves. The local generation of the D diurnal internal tide is the most plausible and likely candidate for the mech- anism generating high vertical wavenumber near-inertial internal waves observed near Tokara Strait. The same harmonic analysis for the diurnal and semi-diurnal tidal elevations during the shipboard observations in November 2016 at Tanegashima tidal station (30.47°N, 130.97°E, Fig. 1a) shows a relatively large diurnal tidal amplitude during the UVMP transect observations in Leg A (November 14–15) and Leg B (November 18–19) (Supplementary Fig. S8c). The near-inertial waves induced by the Kuroshio over the topography cannot be ruled out. Other possibilities, such as the spontaneous generation of the near-inertial waves by the meandering 24–28 31 Kuroshio , and the PSI of wind-induced near-inertial waves are not discussed in the present study due to the limited observation data. More intensive field campaigns, and high-resolution numerical simulations, are necessary to definitively identify the source of these near-inertial waves propagating in the upstream Kuroshio. e b Th reaking of the near-inertial waves is influenced by the Kuroshio. This is because the near-inertial waves can be trapped in the region of negative relative vorticity on the less dense side of the Kuroshio and of strong geo- 32–34 strophic vertical shear below the Kuroshio . The modulation of the lowest internal wave frequency, as known as the effective Coriolis frequency f , by the vertical component of the relative vorticity in the Kuroshio can be eff equated for a meridional two-dimensional front,   ∂v −1 ff = f + f , eff   ∂x (2) where v is the mean meridional velocity and x is the zonal coordinate . Considering the contribution from ver- tical shear of the mean flow to draw down the lowest internal wave frequency further, the equation becomes ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 8 www.nature.com/scientificreports/ −2 Figure 6. Long-term moored ADCP data analysis. The variance preserved spectra [s ] of the 30-h lowpass near-inertial shear (a) as a function of depth [m] and frequency [Hz], and (b) the depth averaged uv + zz spectrum. (c) The time series of the amplitude of 30-h lowpass diurnal tidal elevation, D = K + O [m] at 1 1 1 Nakanoshima tidal station (Fig. 1a) is shown in red, and that for semi-diurnal tidal elevation D = M + S + N 2 2 2 2 is shown in black. The diurnal and semi-diurnal constituents are obtained by harmonic analysis for the hourly tidal elevation record at Nakanoshima tidal station. (d) The time series of 30-h lowpass near-inertial shear amplitude measured at Stn. TK1. uv + zz 21 −       ∂v ∂v ∂b −1    ω = ff  + f −  .    min                ∂x ∂z ∂z (3) es Th e lowest internal wave frequencies are computed with the measured ADCP velocity and buoyancy along the transect averaged over 200 m depth and 7 km for vertical and horizontal direction, respectively, replacing v with ∂v the computed velocity component normal to the ship track, x by the direction parallel to the ship track, and () ∂z ∂u ∂v by total velocity shear square . Because our ship tracks were not normal to the front, the contribu- () () ∂z ∂z ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 9 www.nature.com/scientificreports/ Figure 7. Normalized minimum internal wave frequencies by local Coriolis frequency averaged within the upper 500 m as a function of longitude for (a) Leg A and for (b) Leg B. Blue indicates the minimum internal 33 34 wave frequency for a barotropic front f /f (2) , and red indicates that for a baroclinic front ω /f (3) . Shading eff min represents 95% confidence interval for each average value. Horizontal dashed line indicates unity. tion from the vertical relative vorticity to modulate f and ω is most likely underestimated. Nevertheless, both eff min the average f and ω within the upper 500 m show a similar trend in which the lowest internal wave frequencies eff min are significantly decreased around 129.3–129.5°E in Leg A (Fig.  7a) and 129.5–129.65°E in Leg B (Fig. 7b) from the local Coriolis frequencies f. The reduction of the lowest internal-wave frequency and the slightly larger lowest internal-wave frequency formed to the west, create the convex-down structure in f and ω at several longitude eff min ranges, where the near-inertial waves can be trapped. In these regions, the amplitudes of high vertical wavenum- ber shear seems large, especially at 129.3–129.6°E in Leg A (Figs 2b,c and 7a). The large turbulent dissipation rates are found beneath the Kuroshio at 100–150 m depth and on the less dense side of the Kuroshio below 150 m depth at 129.3–129.6°E in Leg A and at 129.5–129.65°E in Leg B (Fig. 5), where the near-inertial wave trapping favorable structures are coincidently formed (Fig. 7). This implies that the trapping mechanisms by the geostrophic shear may be effective in promoting near-inertial internal-wave breaking and causing enhanced turbulence in the upstream Kuroshio. e in Th duced energy dissipations and vertical mixing could have important implications for watermass modi- 57 −3 fications. A previous study estimated the vertical eddy diffusivity at the salinity minimum layer σ = 26.8 kg m −4 2 −1 of North Pacific Intermediate Water (NPIW) in the Okinawa Trough to be 6.8–21.5 × 10 m s using a steady-state salinity balance with historical current data and salinity data obtained by the profiling floats. Although the direct turbulence data obtained in this study cover only the upper 300 m depth, and did not reach −3 the layer of σ = 26.8 kg m , strong turbulence just above salinity minimum layer of NPIW, which is equivalent to −4 2 −1 57 >O(10 m s ) of average eddy diffusivity, is consistent with a previous estimate . This implies that the observed large amplitude high vertical wavenumber internal waves play very important roles in watermass mixing and modification in the study region. Not only is the observed pronounced turbulent diapycnal mixing important for the watermass modifica- tions of NPIW, it is also very important for nutrient supply from the lower layer to the surface euphotic layer. Recent studies of the nutrient transport in the upstream Kuroshio showed that despite its low-nutrient olig- −1 otrophic surface water, the Kuroshio transports a large amount of nitrate at an average rate of 170.8 kmol s −3 6,7 in the subsurface layer, σ = 24–27 kg m and transports it downstream , similar to the nutrient stream of 8–11 the Gulf Stream . In addition, nitrate concentrations below the Gulf Stream are found to be higher than 8,9 in the ambient water of the same density . This along-isopycnal anomaly of nutrients is most likely formed and maintained by diapycnal diffusion, as opposed to homogenizing mesoscale along-isopycnal stirring. −4 2 −1 Therefore, large average turbulent diffusivity of O(10 m s ) in the subsurface layer at 100–300 m depth, −3 σ = 24–26 kg m , reported by this study, indicates that the observed subsurface turbulence induced by the −3 near-inertial internal waves could effectively diffuse up the subsurface nutrient stream of σ = 24–27 kg m to the shallower layer, generate the elevated nutrient concentrations along the Kuroshio, and increase pri- mary production in its downstream. In this study, for the first time, the spatially coherent banded layers of strong turbulence associated with high vertical wavenumber near-inertial internal waves were directly observed using a new tow-yo microstructure pro- filer. Our analyses suggest that the high vertical wavenumber near-inertial waves are most likely generated by the diurnal tide over the rough topography in Tokara Strait. Because the observation site is around 30°N latitude, where the inertial period is 24 h, the generated diurnal internal tides are inherently near-inertial waves with slow group velocities. The ADCP and microstructure data suggest that these near-inertial internal waves are trapped and broken into turbulence, preferentially on the anticyclonic side of the upstream Kuroshio, suggesting that the trapping of the near-inertial waves by the Kuroshio due to anticyclonic vorticity catalyzes the breaking and dis- sipation of the near-inertial/diurnal internal waves, and promoting watermass modification and nutrient supply. Although the contributions of PSI to form observed high vertical wavenumber near-inertial waves were found to be small compared to the diurnal internal tide in our observations, other mechanisms to generate observed ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 10 www.nature.com/scientificreports/ inertial waves are not discussed in detail in this study. Okinawa Trough and Tokara Strait are very peculiar regions because of their strong baroclinic current, the Kuroshio (which can modify the lowest internal wave frequency), flowing through the shallow seamounts and islands, and through the critical latitude for PSI and the boundary where the diurnal period is identical to the inertial period. More comprehensive field and numerical studies are, therefore, necessary to fully understand the generation, propagation, and dissipation mechanisms of these near-inertial internal waves/diurnal internal tides propagating in the upstream Kuroshio. Methods Back-rotation of the shear. Near-inertial wave velocity and shear rotate clockwise in time in the Northern Hemisphere. When the observation period exceeds or is equivalent to the inertial time scale, this rotation needs to be removed to avoid temporal aliasing of lateral velocity and shear structures. The back-rotated shear Z(t ) = u (t ) + iv (t ) (where u and v are the zonal and meridional shear, respectively, and i is an imaginary unit) 0 z 0 z 0 z z referenced to the time t , is estimated from the observations Z(t) = u (t) + iv (t) at the arbitrary time t as 0 z z it () −tf Zt () = Zt ()e , (4) where f is the mean Coriolis frequency, assuming the rotation at the mean inertial frequency, and back-rotated zonal u (t ) and meridional shear v (t ) are obtained by taking the real and imaginary parts of (4), respectively. z 0 z 0 Tow-yo turbulence survey. We used a new tow-yo Underway Vertical Microstructure Profiler (Underway-VMP: UVMP) to measure high-resolution spatial structures of microscale turbulence at the Kuroshio Front. The UVMP consists of a vertical microstructure profiler, VMP-250 (Rockland Scientific International, Victoria, Canada) and a winch for the Underway-CTD (UCTD, Teledyne Oceanscience, USA). The VMP-250 carries two shear probes, two FP07 thermistors, a pressure sensor, and a vibration sensor to measure microscale velocity shear, microscale temperature gradient, pressure and accelerations, respectively, at a frequency of 512 Hz. −1 e UVMP wa Th s tow-yoed at a ship speed relative to water, at 1.5–2.1 m s . The descending speeds of the UVMP, −1 gradually decreasing with depth, were between 0.3–0.8 m s . Taking 7 minutes to reach 300 m depth and 5 min- utes for recovery, the lateral resolution of the tow-yo survey was 1–2 km, depending on the absolute current veloc- ity. The data were internally recorded and recovered ae ft r each transect observation. However, because the both FP07 thermistors mounted initially on the VMP-250 in Leg A malfunctioned, the expendable bathythermographs (XBT-T7s) were deployed six times to obtain temperature profiles. The XBT temperature data were interpolated by objective mapping with decorrelation lengths of 15 km and 40 m along ship-track and vertical directions, respectively, which were then used to estimate the background potential density ρ and σ and stratification for θ θ Leg A (see below). For Leg B, two functional FP07 thermistors were installed on the VMP-250, so that the poten- tial density ρ and σ could be estimated, based on high-resolution FP07 thermistor data. To cross-calibrate the θ θ FP07 thermistors, the calibrated UCTD sensor and VMP-250 were clamped together and lowered to 200 m depth with the BT-winch of the Kagoshima-maru. To compute turbulent kinetic energy dissipation rates ε, the turbulent shear spectra were computed using the microscale shear data measured by the airfoil shear probe over 8 seconds, at approximately every 5–6 m depth, with an overlap of 4 s (2.5–3 m). The obtained shear spectra were integrated over the wavenumber ranges, where spectra agree with the Nasmyth empirical shear spectrum , avoiding the integrations of electronic noise at high wavenumber ranges and low frequency instrument vibration noise at low wavenumber ranges. For the wave- number ranges where the observed and the empirical spectra did not match, we integrated the fitted Nasmyth empirical shear spectrum (see Fig. 4). Assuming isotropic turbulence, ε was then computed by   15 ∂u   εν = ,     2 ∂z (5) where ν is the kinematic viscosity, and is the turbulent shear variance computed by the integration of the (/ ∂∂ uz ) shear spectrum. Because the UVMP does not carry a conductivity sensor to measure salinity, seawater potential density ρ and −3 σ = ρ − 1000 (kg m ) were estimated using the close relationship between potential temperature and potential θ θ density found at the observation sites (Supplementary Fig. S1). The relationship between potential temperature and potential density was derived from the UCTD data obtained across the front prior to the observations in Legs A and B (Fig. 1a), and from CTD data obtained at 36 stations during the surveys (Fig. 1a). The empirical relation- ship was derived by bin-averaging the observed potential density as a function of potential temperature every 0.5 °C (Supplementary Fig. S1). The potential density ρ and σ were then estimated by linear interpolation with θ θ the derived mean density model as a function of potential temperature. The standard deviation of the mean den- −3 −3 sity in each average bin was as small as <0.03 (kg m ) for σ > 25 (deeper than 150 m depth), 0.04–0.1 (kg m ) −3 for 24 < σ < 25 (100–150 m depth), and ~0.12 (kg m ) for σ < 24 (shallower than 100 m depth); the estimated θ θ density had better accuracy for the denser and deeper layers. Internal wave energy flux and wind power input. e in Th ternal wave energy flux along the vertical axis was estimated based on the wave horizontal kinetic energy HKE, and the vertical group velocity Cg , assuming that the near-inertial waves had low potential energy. Using the dispersion relation of the internal waves, the vertical group velocity for the low frequency internal wave is ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 11 www.nature.com/scientificreports/ f − ω ∣∣ Cg = , ωm (6) −1 −1 where ω is the internal wave frequency (rad s ), and m the internal-wave vertical wavenumber (rad m ). The horizontal kinetic energy of the internal waves is [] uv + zz HKE = ρ , o 2 (7) 2 m −1 where u and v are ADCP zonal and meridional shear (s ), respectively, and [] represents the spatiotempo- z z ral average between 100–250 m depth and over the entire transect for Leg A and Leg B. The internal-wave fre- −1 quency ω was estimated as ω = 1.01 f (rad s ), based on the internal ray path and the measured velocity shear band (Fig. 2) following a previous study . The range of the internal-wave vertical wavelength was estimated as 100–200 m with a mean of 150 m, based on the spectral peaks in the shear spectra (Fig. 3c,d). The corresponding −1 wavenumber range m = 0.063–0.031 and its mean m = 0.042 rad m were used to estimate the group velocity −3 Cg and HKE. e es Th timated horizontal HKE were 2.9–11.5 with a mean of 6.4 J m for Leg A and 6.3–25.2 with −3 −5 −1 a mean of 14.2 J m for Leg B. The range of the vertical group velocity was estimated to be 2.3–4.6 × 10 m s −5 −1 −2 with a mean of 3.5 × 10 m s . The resulting vertical internal wave energy flux for Leg A was 0.07–0.5 mW m −2 −2 −2 with a mean of 0.2 mW m , and that for Leg B was 0.1–1.2 mW m with a mean of 0.5 mW m (Supplementary Fig. S3b). To compare the internal wave energy flux with the wind power input to the near-inertial motions, a previous 48 58 slab model was used with the hourly reanalysis Japan Meteorological Agency GPV-MSM wind data from October 28 through November 30, 2016. Considering that the Kuroshio may have a strong advective ee ff ct and that remote wind-induced inertial motions could propagate to the observation sites, the hourly wind data were first averaged spatially over the region within 26–36°N and 125–131°E. The mixed layer depth was defined by the buoyancy frequency jump from the surface, and found to be 100 m depth where the large amplitude near-inertial shear was observed. Strong winds occurring before field observations on November 8–12, and during observa- tions on November 15 input energy to near-inertial motions (Supplementary Fig. S3). Analyses for the tidal elevation and long term moored ADCP data. Hourly tidal elevation data at Nakanoshima tidal station obtained from the Japan Oceanographic Data Center (JODC) were analyzed from August 31, 2000 through August 14, 2001. Harmonic analysis with 59 tidal constituents was conducted for the tidal elevation record to compare with the 30-h lowpass near-inertial shear amplitude (see below). The results of the harmonic analysis showed that the variance of the dominant three components for semi-diurnal (M , S , and 2 2 N ) and two components for diurnal (K and O ) tidal elevation can explain 99% of the total variance, with the M 2 1 1 2 tidal component accounting for 67%, S for 13%, N for 3%, K for 10%, and O for 6%. To obtain the time series of 2 2 1 1 subinertial diurnal and semi-diurnal tidal elevation by the harmonic analysis, the computed diurnal D = K + O 1 1 1 and semi-diurnal D = M + S + N tidal elevation amplitudes were smoothed with a 30-h lowpass filter and 2 2 2 2 compared with the 30-h lowpass near-inertial shear (Fig. 6). The same harmonic analysis was conducted for the tidal elevation data obtained during the shipboard observations in November 2016, at Tanegashima tidal station (Fig. 1a). The variance of total tidal elevation at Tanegashima station is contributed nearly identically from each dominant tidal constituent as in the Nakanoshima tidal data. To compare the obtained 30-h lowpass semi-diurnal and diurnal tidal amplitudes with the near-inertial shear amplitude, the relatively long-term moored ADCP data were analyzed. The moored ADCP (upward looking 75 kHz 20° beam angle Workhorse Long Ranger ADCP, Teledyne RDI) measured the horizontal current at 14 layers with 32 m vertical bins from 480.5 to 64.5 m depth from September 30, 2000 through July 1, 2001, at Stn. TK1 (Fig. 1a) in the vicinity of the Stn. M01 and M02. The hourly ADCP data were differentiated vertically to obtain the vertical shear. To extract the near-inertial shear, a bandpass Butterworth filter with cutoff frequencies at 0.8 f and 1.1 f (equivalent to the periods of 30 h and 21.7 h, respectively) was used. e Th near-inertial amplitude was computed using near-inertial zonal and meridional shear as uv + . The subinertial modulations of the zz near-inertial shear amplitude was then obtained using a lowpass Butterworth filter at 30 h for the uv + zz bandpass near-inertial shear amplitude. Data Availability. The datasets generated and/or analyzed during the current study are available from the corresponding author upon request. References 1. Qiu, B. Kuroshio and Oyashio Currents In: Encyclopedia of Ocean Sciences, Steele, J. H., Thorpe, (Oceanographer), S. A. & Turekian, Karl K. (Eds) (Academic Press, 2001). 2. Uda, M. Researches on “siome” or current rip in the seas and oceans. Geophys. Mag 11, 307–372 (1938). 3. Rudnick, D. L. Intensive surveys of the Azores Front. Part II: Inferring the geostrophic and vertical velocity fields. J. Geophys. Res. 101, 16,291–16,303 (1996). 4. Rudnick, D. L. & Luyten, J. R. Intensive surveys of the Azores Front part I. Tracers and dynamics. J. Geophys. Res. 101, 923–939 (1996). 5. Nurser, A. J. G. & Zhang, J. W. Eddy-induced mixed layer shallowing and mixed layer/thermocline exchange. J. Geophys. Res. 105, 21,851–21,868 (2000). 6. Guo, X., Zhu, X. H., Wu, Q. S. & Huang, D. The Kuroshio nutrient stream and its temporal variation in the East China Sea. Glob. Biogeochem. Cycles 117, C01026, https://doi.org/10.1029/2011JC007292 (2012). ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 12 www.nature.com/scientificreports/ 7. Guo, X., Zhu, X. H., Long, Y. & Huang, D. Spatial variations in the Kuroshio nutrient transport from the East China Sea to south of Japan. Biogeosciences 10, 6403–6417 (2013). 8. Pelegrí, J. L. & Csanady, G. T. Nutrient transport and mixing in the Gulf Stream. J. Geophys. Res. 96, 2577–2583 (1991). 9. Pelegrí, J. L., Csanady, G. T. & Martins, A. The North Atlantic nutrient stream. J. Ocean 52, 275–299 (1996). 10. Williams, R. G., Roussenov, V. & Follows, M. Induction of nutrients into the mixed layer and maintenance of high latitude productivity. Glob. Biogeochem. Cycles 20, GB1016, https://doi.org/10.1029/2005GB002586 (2006). 11. Williams, R. G. et al. Nutrient streams in the North Atlantic: Advective pathways of inorganic and dissolved organic nutrients. Glob. Biogeochem. Cycles 25, GB4008, https://doi.org/10.1029/2010GB003853 (2011). 12. Nagai, T., Aiba, M. & Clayton, S. Multiscale pathway to supply nutrients to the Kuroshio (In Japanese). Kaiyo to Seibutsu 37, 469–477 (2015). 13. Nagai, T. & Clayton, S. Nutrient interleaving below the mixed layer of the Kuroshio Extension Front. Ocean. Dyn. 67, 1027–1046, https://doi.org/10.1007/s10236-017-1070-3 (2017). 14. Rhines, P. B. & Young, W. R. Homogenization of potential vorticity in planetary gyres. J. Fluid Mech. 122, 347–367 (1982). 15. Oakey, N. S. & Elliot, J. A. Vertical temperature gradient structure across the Gulf Stream. J. Geophys. Res. 82, 1369–1380 (1977). 16. Gregg, M. C. & Sanford, T. B. Signatures of mixing from Bermuda slope, the Sargasso Sea and the Gulf Stream. J. Phys. Ocean. 10, 105–127 (1980). 17. Gargett, A. E. & Osborn, T. R. Small-scale shear measurements during the Fine and Microstructure Experiment (Fame). J. Geophys. Res. Ocean. 86, 1929–1944, https://doi.org/10.1029/JC086iC03p01929 (1981). 18. Hasegawa, D., Yamazaki, H., Lueck, R. G. & Seuront, L. How islands stir and fertilize the upper ocean. Geophys. Res. Lett. 31, n/a– n/a, https://doi.org/10.1029/2004GL020143 (2004). 19. Nikurashin, M. & Ferrari, R. Radiation and dissipation of internal waves generated by geostrophic motions impinging on small-scale topography: Theory. J. Phys. Ocean. 40, 1055–1074 (2010). 20. Tsutsumi, E. et al. Turbulent mixing within the Kuroshio in the Tokara Strait. J. Geophys. Res. Ocean. n/a–n/a, https://doi. org/10.1002/2017JC013049 (2017). 21. Lien, R.-C. & Gregg, M. C. Observations of turbulence in a tidal beam and across a coastal ridge. J. Geophys. Res. Ocean. 106, 4575–4591, https://doi.org/10.1029/2000JC000351 (2001). 22. Lueck, R. R. & Mudge, T. D. Topographically induced mixing around a shallow seamount. Sci. 276, 1831–1833, https://doi. org/10.1126/science.276.5320.1831 (1997). 23. Rainville, L. & Pinkel, R. Observations of energetic high-wavenumber internal waves in the Kuroshio. J. Phys. Oceanogr. 34, 1495–1505 (2004). 24. Plougonven, R. & Snyder, C. Inertia-gravity waves spontaneously generated by jets and fronts. Part I: Different baroclinic life cycles. J. Atmos. Sci. 64, 2502–2520 (2007). 25. Alford, M. H., Shcherbina, A. Y. & Gregg, M. C. Observations of near-inertial gravity waves radiating from a frontal jet. J. Phys. Oceanogr. 43, 1225–1239 (2013). 26. Shakespeare, C. J. & Taylor, J. R. The spontaneous generation of inertia-gravity waves during frontogenesis forced by large strain: theory. J. Fluid. Mech. 757, 817–853 (2014). 27. Nagai, T., Tandon, A., Kunze, E. & Mahadevan, A. Spontaneous generation of near-inertial waves by the Kuroshio Front. J. Phys. Oceanogr. 45, 2381–2406, https://doi.org/10.1175/jpo-d-14-0086.1 (2015a). 28. Shakespeare, C. J. & Hogg, A. M. Spontaneous surface generation and interior amplification of internal waves in a regional-scale ocean model. J. Phys. Oceanogr. 47, 811–826, https://doi.org/10.1175/JPO-D-16-0188.1 (2017). 29. Hibiya, T., Niwa, Y. & Fujiwara, K. Numerical experiments of nonlinear energy transfer within the oceanic internal wave spectrum. J. Geophys. Res. Ocean. 103, 18715–18722, https://doi.org/10.1029/98JC01362 (1998). 30. MacKinnon, J. A. & Winters, K. B. Subtropical catastrophe: Significant loss of low-mode tidal energy at 28.9°. Geophys. Res. Lett. 32, n/a–n/a, https://doi.org/10.1029/2005GL023376 (2005). 31. Thomas, L. N. & Taylor, J. R. Damping of inertial motions by parametric subharmonic instability in baroclinic currents. J. Fluid Mech. 743, 280–294, https://doi.org/10.1017/jfm.2014.29 (2014). 32. Mooers, C. N. Several effects of baroclinic current on the cross-stream propagation of inertial-internal waves. Geophys. Fluid. Dyn. 6, 245–275 (1975). 33. Kunze, E. Near-inertial wave propagation in geostrophic shear. J. Phys. Oceanogr. 15, 544–565 (1985). 34. Whitt, D. & Thomas, L. Near-inertial waves in strongly baroclinic currents. J. Phys. Oceanogr. 43, 706–725 (2013). 35. Kunze, E. & Sanford, B. Observations of near-inertial waves in a front. J. Phys. Oceanogr. 15, 566–581 (1984). 36. Kunze, E. The mean and near-inertial velocity fields in a warm-core ring. J. Phys. Oceanogr. 16, 1444–1461 (1986). 37. Lueck, R. G. & Osborn, T. The dissipation of kinetic energy in a warm-core ring. J. Geophys. Res. 91, 803–818 (1986). 38. Winkel, D. P., Gregg, M. C. & Sanford, T. B. Patterns of shear and turbulence across the Florida Current. J. Phys Oceanogr. 32, 3269–3285 (2002). 39. Nagai, T., Tandon, A., Yamazaki, H. & Doubell, M. J. Evidence of enhanced turbulent dissipation in the frontogenetic Kuroshio Front thermocline. Geophys. Res. Lett. 36, L12609, https://doi.org/10.1029/2009GL038832 (2009). 40. Nagai, T., Tandon, A., Yamazaki, H., Doubell, M. J. & Gallager, S. Direct observations of microscale turbulence and thermohaline structure in the Kuroshio Front. J. Geophys. Res. 117, C08013, https://doi.org/10.1029/2011JC00722 (2012). 41. Kaneko, H., Yasuda, I., Komatsu, K. & Itoh, S. Observations of the structure of turbulent mixing across the Kuroshio. Geophys. Res. Lett. 39, n/a–n/a, https://doi.org/10.1029/2012GL052419 (2012). 42. Kaneko, H., Yasuda, I., Komatsu, K. & Itoh, S. Observations of vertical turbulent nitrate flux across the Kuroshio. Geophys. Res. Lett. 40, 3123–3127, https://doi.org/10.1002/grl.50613 (2013). 43. Inoue, R., Gregg, M. C. & Harcourt, R. R. Mixing rates across the Gulf Stream, Part 1: On the formation of Eighteen Degree Water. J. Mar. Res. 68, 643–671 (2010). 44. Alford, M. H. & Gregg, M. C. Near-inertial mixing: Modulation of shear, strain and microstructure at low latitude. J. Geophys. Res. Ocean. 106, 16947–16968, https://doi.org/10.1029/2000JC000370 (2001). 45. Alford, M. H., MacKinnon, J. A., Simmons, H. L. & Nash, J. D. Near-inertial internal gravity waves in the ocean. Annu. Rev. Mar. Sci. 8, 95–123, https://doi.org/10.1146/annurev-marine-010814-015746, PMID: 26331898 (2016). 46. Garrett, C. J. R. & Munk, W. H. Space-time scales of internal waves. Geophys. Fluid Dyn. 2, 225–264 (1972). 47. Garrett, C. J. R. & Munk, W. H. Space-time scales of internal waves. A progress report. J. Geophys. Res. 80, 291–323 (1975). 48. Pollard, R. T. & Millard, R. C. Comparison between observed and simulated wind-generated inertial oscillations. Deep. Res. 17, 813–821 (1970). 49. D’Asaro, E. Wind-forced internal waves in the North Pacific and Sargasso Sea. J. Phys. Ocean. 14, 781–794 (1984). 50. Alford, M. H. Redistribution of energy available for ocean mixing by long-range propagation of internal waves. Nat. 423, 159–162 (2003). 51. Nasmyth, P. W. Oceanic turbulence. Ph.D. thesis, University of British Columbia (1970). 52. Kunze, E., Firing, E., Hummon, J. M., Chereskin, T. K. & Thurnherr, A. M. Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles. J. Phys. Oceanogr. 36, 1553–1576, https://doi.org/10.1175/JPO2926.1 (2006). 53. Whalen, C. B., Talley, L. D. & MacKinnon, J. A. Spatial and temporal variability of global ocean mixing inferred from Argo profiles. Geophys. Res. Lett. 39, n/a–n/a, https://doi.org/10.1029/2012GL053196 (2012). ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 13 www.nature.com/scientificreports/ 54. Osborn, T. Estimates of the local rate of vertical diffusion from dissipation measurements. J. Phys. Oceanogr. 10, 83–89 (1980). 55. Xie, X.-H., Shang, X.-D., van Haren, H., Chen, G.-Y. & Zhang, Y.-Z. Observations of parametric subharmonic instability-induced near-inertial waves equatorward of the critical diurnal latitude. G e o p hy s . R es . L et t . 38, n/a–n/a, ht t ps://doi . org/10.1029/2010GL046521 (2011). 56. Liu, T., Wang, B., Hirose, N., Yamashiro, T. & Yamada, H. High-resolution modeling of the Kuroshio Current power south of Japan. Submitted to Journal of Ocean Engineering and Marine Energy (2017). 57. Nakamura, H. et al. Intermediate and deep water formation in the Okinawa Trough. J. Geophys. Res. Ocean. 118, 6881–6893, https:// doi.org/10.1002/2013JC009326 (2013). 58. Saito, K. et al. The operational JMA nonhydrostatic mesoscale model. Mon. Weather. Rev. 134, 1266–1298 (2006). Acknowledgements We thank Capt. Uchiyama, C/O Azuma, Prof. Nishina, Prof. Kobari, Prof. Matsuno, Prof. Senju, Dr. Li, crews, and all the participants of the R/T/V Kagoshima-maru cruise KG1616. T.N. thanks Hancyk, Dr. Wolk at RSI, Dr. Li, Seki, and Yoshii at JAC, Dr. Kokubu at NORTEK, Rutka at STS, support from OMIX (MEXT KAKENHI JP16H01590) and SKED (“Study of Kuroshio Ecosystem Dynamics for Sustainable Fisheries” funded to FRA from MEXT). D.H., T.T., and R.I. thank OMIX (MEXT KAKENHI) JP15H05818. H.N., and E.T. thank OMIX (MEXT KAKENHI) JP15H05821. Author Contributions T.N. conceived the tow-yo turbulence experiments, D.H., T.N., T.T. materialized the tow-yo microstructure profiling, H.N. organized the cruise, E.T. analysed moored ADCP data and conducted harmonic analyses, the authors (D.H., T.T., H.N., E.T.) conducted the shipboard experiments, T.Y. conducted the long-term mooring observation, R.I. provided a slab model for wind analysis, and T.N. analysed the results. All authors reviewed and edited the manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-15167-1. Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 ScientiFic REpo R TS | 7:14555 | DOI:10.1038/s41598-017-15167-1 14

Journal

Scientific ReportsSpringer Journals

Published: Nov 6, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off