Data Sheet 1_An energy-constrained profile parameterization of shear-driven turbulence in the interior ocean.docx

This paper proposes an energy-constrained profile parameterization of both turbulent kinetic energy dissipation rate (ϵ) and vertical diffusivity (κ), for shear instability-induced turbulence that is initiated in an initial unstable layer (IUL) where the gradient Richardson number Ri ∈ (0, 0.25). Large-eddy simulation (LES) experiments provide the data of turbulent processes originating from Kelvin-Helmholtz instability of varied initial shear conditions. The energy-constrained framework posits ϵ and κ as proportional to  Ka and τ–1, where  Ka represents available kinetic energy, measuring the released kinetic energy, τ denotes turbulence evolution timescale. Both are determinable by the thickness of IUL (h0), buoyancy frequency (N0), vertical shear (S0), and Richardson number (Ri0) of the IUL. Notably, unlike conventional schemes that parameterize turbulent mixing for single model grid point layer by layer, the present scheme parameterizes the turbulent mixing not only for the grid point(s) of IUL, but also for all the model grid points that are within a determined vertical turbulent penetration layer, by providing a profile of diffusivity. Therefore, the scheme is termed the energy-constrained profile parameterization (EPP). EPP aligns well with the LES results and direct microstructure measurements, outperforming existing parameterizations.