Eddy-Induced Stratified Turbulence

• Eddy-induced stratified turbulence is defined by anisotropic flows in stably stratified fluids, where background density gradients, eddy instabilities, and wave interactions coexist.
• The phenomenon is characterized by key instabilities such as Kelvin–Helmholtz, Holmboe, and inertial instabilities, analyzed using the Boussinesq framework and dimensionless parameters like Fr, Ro, and Ri.
• This turbulence underpins diapycnal mixing and energy redistribution in geophysical systems, offering improved parameterizations for oceanic, atmospheric, and planetary interior dynamics.

Eddy-induced stratified turbulence comprises turbulent flows in stably stratified fluids where stratification—set by a background potential density gradient—suppresses vertical mixing, enforcing anisotropy, while eddy structures, instabilities, and wave interactions inject and redistribute energy. This regime is central to geophysical and astrophysical systems, such as oceans, atmospheres, and planetary interiors, where the competition among inertia, buoyancy, and rotation produces a diverse array of flow structures including layerwise jets, internal waves, coherent vortices, and intermittent turbulence. Recent research has elucidated key dynamical regimes, spectral features, scaling laws, and parameterizations that enable quantitative predictions and improved modeling of turbulent mixing and transport in stratified environments.

1. Governing Equations and Control Parameters

The canonical framework for stratified turbulence is the Boussinesq system, extended where relevant to include rotation:

$$\begin{aligned} &\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{u} + f\mathbf{e}_z \times \mathbf{u} = - \frac{1}{\rho_0}\nabla p' + \nu \nabla^2\mathbf{u} + b\mathbf{e}_z \ &\frac{\partial b}{\partial t} + (\mathbf{u} \cdot \nabla)b = \kappa \nabla^2 b, \quad \nabla\cdot\mathbf{u}=0 \end{aligned}$$

where the velocity $\mathbf{u}$, buoyancy $b = -(g/\rho_0)\rho'$, kinematic viscosity $\nu$, and diffusivity $\kappa$ characterize the fluid dynamics; stratification is set by the Brunt–Väisälä frequency $N = \sqrt{-(g/\rho_0)\partial \rho/\partial z}$, and rotation by the Coriolis parameter $f$.

Dimensionless groups governing the behavior include:

• Froude number $\mathrm{Fr}=U/(NL)$: ratio of inertial to buoyancy forces.
• Rossby number $\mathrm{Ro}=U/(fL)$: ratio of inertial to rotational effects.
• Richardson number $\mathrm{Ri}=N^2/S^2$: ratio of stable stratification to shear.
• Buoyancy Reynolds number $\mathrm{Re}_b=Re\,Fr^2=U^2/(\nu N^2)$: measures turbulence intensity relative to stratification.

These parameters delineate regimes, e.g., strong stratification ($\mathrm{Fr}\ll 1$), intense rotation ($\mathrm{Ro}\ll 1$), and high turbulence ($\mathrm{Re}_b\gg 1$).

2. Instabilities and Turbulence Generation Mechanisms

Eddy-induced turbulence arises from the destabilization and nonlinear interaction of coherent structures in stratified environments.

• Kelvin–Helmholtz Instability (KHI): Vertical shear in horizontal velocity, common in boundary layers and wake flows, triggers turbulence when the local gradient Richardson number $Ri_g=N^2/S_v^2<1/4$. Strong stratification ($\mathrm{Fr}\ll1$) thins vortex layers and enhances local shear, promoting KHI unless suppressed by rapid rotation, which organizes flow into near-vertical Taylor columns (Liu et al., 10 Feb 2025).
• Centrifugal/Inertial Instability (CI): Occurs in regions with negative absolute vorticity ($f+\omega_z<0$), such as anticyclonic wakes. CI is active at moderate stratification ($\mathrm{Fr}\sim0.3-0.4$) and Rossby numbers near unity; strong rotation ($\mathrm{Ro}\ll1$) or elongation of vertical scales by stratification suppresses CI.
• Holmboe Instability: Layered density and velocity profiles can yield interfacial instabilities (e.g., Holmboe waves) that contribute to anisotropic turbulence in exchange flows (Lefauve et al., 2021).
• Internal Wave Breaking: Nonlinear interactions among waves or between waves and eddies lead to local overturns, particularly when internal gravity wave amplitudes become sufficient for convective instability ($Ri<1/4$ locally) (Yokoyama et al., 2019, Marino et al., 2015).

In aggregate, the interplay between stratification, rotation, shear, and boundary geometries (e.g., topography or inclined ducts) sets the spatial and temporal intermittency, spectral character, and anisotropy of turbulence.

3. Anisotropy, Energy Partition, and Spectral Transitions

Stratification fundamentally breaks vertical isotropy: vertical velocity fluctuations and scales become strongly suppressed compared to horizontal components. Several research efforts quantify this onset and degree of anisotropy:

• Global Anisotropy Parameter ($A$): Defined as the ratio of (mean-square) horizontal to vertical velocity fluctuations, $$A = \langle v_\perp^2\rangle / (2\langle v_3^2\rangle)$$. For small Richardson number, $A-1$ increases linearly with $\mathrm{Ri}$, as analytically predicted and numerically confirmed (Bhattacharjee et al., 2019).
• Spectral Regimes: Rotating stratified turbulence divides into a large-scale wave–eddy interaction regime and a small-scale inertia–gravity wave regime. A critical wavenumber, $k_R\sim\mathrm{Fr}^{-1}$, marks the transition: for $k<k_R$, slow (vortical) modes dominate; for $k>k_R$ inertia–gravity waves become prevalent, and energy approaches equipartition between kinetic and potential forms (Marino et al., 2015).
• Anisotropic Spectral Decomposition: Wave versus eddy dominance in spectral space may be diagnosed via indices such as the wave-to-total energy ratio, wave versus potential balance, and polarization anomaly. The transition is best captured using two-dimensional $(k_h, k_z)$ spectra and a criterion $\chi(k_h,k_z)=\omega\,T_{\rm eddy}<1/3$ (Yokoyama et al., 2019); this demarcation is substantially sharper than isotropic partition at the Ozmidov scale $k_O$.

Table: Key Scales and Regimes in Stratified Turbulence

Regime	Transition/Diagnostic	Scaling Law
Eddy-dominated	$k<k_R \sim Fr^{-1}$	$E_k\sim k^{-11/5}$
Wave-dominated	$k>k_R$	$E_k/E_p \sim$ const (equipartition)
Buoyancy/Vertical Layer	$\ell_B \sim Fr L_0$

4. Scaling Laws and Multiscale Structure

Recent asymptotic theory and direct numerical simulations identify two central regimes for vertical velocity and length scales in strongly stratified turbulence (Garaud et al., 2024):

• Single-Scale Model (SSA):
  • Assumes strictly anisotropic flow with vertical aspect ratio $\alpha = L_v/L_h \ll 1$.
  • For advection-dominated ($Pe_b \gg 1$): $L_v \sim Fr\,L_*$, $w \sim Fr\,U_*$.
  • For diffusion-dominated ($Pe_b \ll 1$): $L_v \sim (Fr^2/Pe)^{1/4}\,L_*$, $w \sim (Fr^2/Pe)^{1/4}\,U_*$.
• Multiscale Model (MSA):
  • Accommodates coexisting large-scale anisotropic and small-scale isotropic motions (patches).
  • For advection-dominated: $w_{rms}^{\rm turb} \sim Fr^{1/2}U_*$ in turbulent patches.
  • For diffusion-dominated: $w_{rms}^{\rm turb} \sim (Fr^2/Pe)^{1/6}U_*$.
  • Volume fraction of isotropic patches decreases with increasing stratification ($Fr^{-1}$); as $Re_b \to O(1)$, turbulence collapses to the SSA regime.

These scalings rationalize observed dichotomies between quiescent, layer-like flows and intermittent, deeply turbulent patches.

5. Turbulence Modeling and Closure Schemes

Turbulence closures for stratified flows must represent the coupled energetics and anisotropy produced by stratification and shear. The EFB (Energy- and Flux-Budget) hierarchy provides a comprehensive framework (Zilitinkevich et al., 2011):

• Budget Equations: Separate TKE and TPE, with explicit conversion by buoyancy flux and shear production.
• Prognostic Models: Ranging from five-equation closures (prognosing $E_K$, $E_P$, stress, flux, and time scale) to minimal energy-only versions.
• Turbulent Prandtl Number: $Pr_t \approx 0.8$ for strong turbulence ($\mathrm{Ri}\ll 1$), but rises as $Pr_t \sim \mathrm{Ri}/R_\infty$ ($R_\infty \approx 0.25$) in strongly stratified ($\mathrm{Ri}\gg 1$) regimes.
• Flux Richardson Number: Tied closely to the mixing efficiency and bounds of shear-supported turbulence.
• Parameter Recommendations: For vigorous geophysical turbulence ($Re_b\gg30, Fr_t\ll1$) expect $Pr_T\approx3$, mixing efficiency $\Gamma\approx R_f\approx0.05$, and flux parameterizations $\kappa_T\sim0.2\,\varepsilon/N^2$ (Lefauve et al., 2021).

6. Spectral Fluxes, Observational Signatures, and Power Laws

Cospectral behavior at high wavenumbers under stratification is critical for flux estimation in models and observations:

• Cospectral Power Laws: Classical dimensional arguments yield $C_{w\theta}(k) \sim k^{-7/3}$, but measurements consistently find a $k^{-2}$ scaling at high $k$ in stably stratified flows (Cheng et al., 2018). Adjusting flux corrections in eddy-covariance and model closures to the –2 law improves missing-flux estimates.
• Measurement Techniques: In the deep Mediterranean, dissipation rates obtained via Thorpe-scale overturns and band-pass–filtered Ellison scales show that sub-mesoscale eddy and wave-induced turbulence exceeds geothermal and open-ocean interior levels by factors of $3$–$10$ (Haren, 7 Jan 2026).
• Internal-Wave–Eddy Interplay: At finite $Ro$ and $Fr$, a significant fraction of energy originating in turbulence can be radiated away by inertial-gravity waves, especially for $Fr \gtrsim Ro$. This partitioning is quantitatively tracked in high-resolution simulations and fits linear wave predictions for the tilt and propagation of energy columns (Li et al., 2023).

7. Geophysical and Practical Implications

Eddy-induced stratified turbulence is a primary driver of diapycnal mixing, energy redistribution, and nutrient/chemical transport in the stratified interiors of oceans, lakes, and the stably stratified atmospheric boundary layer:

• Deep-Ocean/Ecosystem Impact: Observations indicate that episodic submesoscale eddy and internal-wave–driven events elevate turbulence and mixing, dominating nutrient and oxygen supply to abyssal layers over geothermal or convective contributions (Haren, 7 Jan 2026).
• Topographic Wakes and Submesoscale Mixing: The interplay between rotation and stratification modulates turbulence intensity and instability dominance past obstacles, constraining mixing parameterizations for global models (Liu et al., 10 Feb 2025).
• Prediction and Parameterization: Parsimonious relations such as $\varepsilon \sim C(Fr,Ro) U^3/h$ with regime-dependent $C$ provide practical approaches for subgrid mixing in ocean and atmospheric models, with consistent behavior across a range of $Fr,Ro$.

In sum, eddy-induced stratified turbulence embodies multiscale, anisotropic, and regime-dependent dynamics crucial for accurate representation of mixing and transport in natural stratified flows. The interplay of deterministic instabilities, turbulence, and wave dynamics necessitates both high-resolution observational analysis and rigorously constructed dynamical models, supported by recent advances in multiscale theory, turbulence closure, and spectral diagnostics.