Search for: All records

Abstract The shoaling of high-amplitude internal solitary waves (ISWs) of depression in the South China Sea (SCS) is examined through large-scale parallel turbulence-resolving high-accuracy/resolution simulations. A select, near-isobath-normal, bathymetric transect of the gentle SCS continental slope is employed together with stratification and current profiles obtained by in situ measurements. Three simulations of separate ISWs with initial deep-water amplitudes in the range [136, 150 m] leverage a novel wave-tracking capability for a propagation distance of 80 km and accurately reproduce key features of in situ–observed phenomena with significantly higher spatiotemporal resolution. The interplay between convective and shear instability and the associated turbulence formation and evolution as a function of deep-water ISW amplitude are further studied in part revealing processes previously not observed in the field. Across all three waves, the convective instability develops in a similar fashion. Heavier water entrained from the wave rear plunges into its interior, giving rise to transient, yet distinct, subsurface vortical structures. Ultimately, a gravity current is triggered which horizontally advances through the wave interior and mixes it down to pycnocline’s base. Although the waveform remains distinctly symmetric, Kelvin–Helmholtz billows emerge near the well-mixed ISW trough, disturb the wave’s trailing edge, and give rise to an active wake. The evolution of the kinetic energy associated with fine-scale perturbations to the ISW-induced velocity field shows two different growth regimes, each dominated by either convective or shear instability. The wake’s perturbation kinetic energy is nonlinearly dependent on deep-water wave amplitude and can become a sizable fraction of the kinetic energy of the deep-water ISW. 

Abstract Physical processes behind flow‐topography interactions and turbulent transitions are essential for parameterization in numerical models. We examine how the Kuroshio cascades energy into turbulence upon passing over a seamount, employing a combination of shipboard measurements, tow‐yo microstructure profiling, and high‐resolution mooring. The seamount, spanning 5 km horizontally with two summits, interacts with the Kuroshio, whose flow speed ranges from 1 to 2 m s⁻¹, modulated by tides. The forward energy cascade process is commenced by forming a train of 2–3 nonlinear lee waves behind the summit with a wavelength of 0.5–1 km and an amplitude of 50–100 m. A train of Kelvin‐Helmholtz (KH) billows develops immediately below the lee waves and extends downstream, leading to enhanced turbulence. The turbulent kinetic energy dissipation rate isO(10⁻⁷–10⁻⁴) W kg⁻¹, varying in phase with the upstream flow speed modulated by tides. KH billows occur primarily at the lee wave's trailing edge, where the combined strong downstream shear and low‐stratification recirculation trigger the shear instability,Ri < 1/4. The recirculation also creates an overturn susceptible to gravitational instability. This scenario resembles the rotor, commonly found in atmospheric mountain waves but rarely observed in the ocean. A linear stability analysis further suggests that critical levels, where the KH instability extracts energy from the mean flow, are located predominantly at the strong shear layer of the lee wave's upwelling portion, coinciding with the upper boundary of the rotor. These novel observations may provide insights into flow‐topography interactions and improve physics‐based turbulence parameterization. 

Abstract Generating mechanisms and parameterizations for enhanced turbulence in the wake of a seamount in the path of the Kuroshio are investigated. Full-depth profiles of finescale temperature, salinity, horizontal velocity, and microscale thermal-variance dissipation rate up- and downstream of the ∼10-km-wide seamount were measured with EM-APEX profiling floats and ADCP moorings. Energetic turbulent kinetic energy dissipation ratesand diapycnal diffusivitiesabove the seamount flanks extend at least 20 km downstream. This extended turbulent wake length is inconsistent with isotropic turbulence, which is expected to decay in less than 100 m based on turbulence decay time ofN⁻¹∼ 100 s and the 0.5 m s⁻¹Kuroshio flow speed. Thus, the turbulent wake must be maintained by continuous replenishment which might arise from (i) nonlinear instability of a marginally unstable vortex wake, (ii) anisotropic stratified turbulence with expected downstream decay scales of 10–100 km, and/or (iii) lee-wave critical-layer trapping at the base of the Kuroshio. Three turbulence parameterizations operating on different scales, (i) finescale, (ii) large-eddy, and (iii) reduced-shear, are tested. Averageεvertical profiles are well reproduced by all three parameterizations. Vertical wavenumber spectra for shear and strain are saturated over 10–100 m vertical wavelengths comparable to water depth with spectral levels independent ofεand spectral slopes of −1, indicating that the wake flows are strongly nonlinear. In contrast, vertical divergence spectral levels increase withε. 

Abstract Finescale properties of Kelvin–Helmholtz (KH)-like shear instabilities on the trailing edge of a nonlinear lee wave generated by the Kuroshio impinging on a seamount were measured using a towed CTD chain, shipboard ADCP, and echosounder. Lee-wave vertical velocity amplitudes vary in phase with the upstream semidiurnal along-stream current. The instabilities are analogous to atmospheric billows induced by a recirculation on the trailing edge of mountain lee waves. A total of 135 KH billows were identified in a 4-day-long time series roughly 300 m downstream of the center of the lee wave. The KH billows have heightsH= 52 ± 11 m, widthsL= 162 ± 72 m, and aspect ratiosH/L= 0.39 ± 0.18. Positive reduced shear squaredS²− 4N²(whereSis the vertical shear magnitude andNis the buoyancy frequency) in the shear-stratified billows suggests actively growing instabilities, with comparable contributions from across- and along-flow vertical shear. Billow cores are convectively unstable (N²< 0). Large turbulent kinetic energy dissipation ratesare inferred from density overturns. Density, shear, and inferred turbulence properties vary with billow aspect ratios. AsH/Lincreases, density gradients smear out. For 122 billows withH/L< 0.6, dissipation rates increase by one order of magnitude with increasingH/L. These observations of ∼1-m vertical and ∼5-m horizontal resolution billow structures and density overturn dissipation rates can provide a reference for future high-Reynolds-number direct numerical simulations. 

Abstract Horizontal and vertical wavenumbers ( k x , k z ) immediately below the Ozmidov wavenumber ( N 3 / ε ) 1/2 are spectrally distinct from both isotropic turbulence ( k x , k z > 1 cpm) and internal waves as described by the Garrett–Munk (GM) model spectrum ( k z < 0.1 cpm). A towed CTD chain, augmented with concurrent Electromagnetic Autonomous Profiling Explorer (EM-APEX) profiling float microstructure measurements and shipboard ADCP surveys, are used to characterize 2D wavenumber ( k x , k z ) spectra of isopycnal slope, vertical strain, and isopycnal salinity gradient on horizontal wavelengths from 50 m to 250 km and vertical wavelengths of 2–48 m. For k z < 0.1 cpm, 2D spectra of isopycnal slope and vertical strain resemble GM. Integrated over the other wavenumber, the isopycnal slope 1D k x spectrum exhibits a roughly +1/3 slope for k x > 3 × 10 −3 cpm, and the vertical strain 1D k z spectrum a −1 slope for k z > 0.1 cpm, consistent with previous 1D measurements, numerical simulations, and anisotropic stratified turbulence theory. Isopycnal salinity gradient 1D k x spectra have a +1 slope for k x > 2 × 10 −3 cpm, consistent with nonlocal stirring. Turbulent diapycnal diffusivities inferred in the (i) internal wave subrange using a vertical strain-based finescale parameterization are consistent with those inferred from finescale horizonal wavenumber spectra of (ii) isopycnal slope and (iii) isopycnal salinity gradients using Batchelor model spectra. This suggests that horizontal submesoscale and vertical finescale subranges participate in bridging the forward cascade between weakly nonlinear internal waves and isotropic turbulence, as hypothesized by anisotropic turbulence theory. 

Observations on the lee of a topographic ridge show that the turbulence kinetic energy (TKE) dissipation rate due to shear instabilities is three orders of magnitude higher than the typical value in the open ocean. Laboratory-scale studies at low Reynolds number suggest that high turbulent dissipation occurs primarily within the core region of shear instabilities. However, field-scale studies indicate that high turbulence is mainly populated along the braids of shear instabilities. In this study, a high-resolution, resolving the Ozmidov-scale, non-hydrostatic model with Large Eddy Simulation (LES) turbulent closure is applied to investigate dominant mechanisms that control the spatial and temporal scales of shear instabilities and resulting mixing in stratified shear flow at high Reynolds number. The simulated density variance dissipation rate is elevated in the cusp-like bands of shear instabilities with a specific period, consistent with the acoustic backscatter taken by shipboard echo sounder. The vertical length scale of each cusp-like band is nearly half of the vertical length scale of the internal lee wave. However, it is consistent with instabilities originating from a shear layer based on linear stability theory. The model results indicate that the length scale and/or the period of shear instabilities are the key parameters to the mixing enhancement that increases with lateral Froude number Fr L , i.e. stronger shear and/or steeper ridge.