Search for: All records

This work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

 

Abstract Direct numerical simulations are performed to compare the evolution of turbulent stratified shear layers with different density gradient profiles at a high Reynolds number. The density profiles include uniform stratification, two-layer hyperbolic tangent profile and a composite of these two profiles. All profiles have the same initial bulk Richardson number ( $$Ri_{b,0}$$ R i b , 0 ); however, the minimum gradient Richardson number and the distribution of density gradient across the shear layer are varied among the cases. The objective of the study is to provide a comparative analysis of the evolution of the shear layers in term of shear layer growth, turbulent kinetic energy as well as the mixing efficiency and its parameterization. The evolution of the shear layers in all cases shows the development of Kelvin–Helmholtz billows, the transition to turbulence by secondary instabilities followed by the decay of turbulence. Comparison among the cases reveals that the amount of turbulent mixing varies with the density gradient distribution inside the shear layer. The minimum gradient Richardson number and the initial bulk Richardson number do not correlate well with the integrated TKE production, dissipation and buoyancy flux. The bulk mixing efficiency for fixed $$Ri_{b,0}$$ R i b , 0 is found to be highest in the case with two-layer density profile and lowest in the case with uniform stratification. However, the parameterizations of the flux coefficient based on buoyancy Reynolds number and the ratio of Ozmidov and Ellison scales show similar scaling in all cases. 

Observations in the Pacific Equatorial Undercurrents (EUC) show that the nighttime deep-cycle turbulence (DCT) in the marginal-instability (MI) layer of the EUC exhibits seasonal variability that can modulate heat transport and sea surface temperature. Large-eddy simulations (LES), spanning a wide range of control parameters, are performed to identify the key processes that influence the turbulent heat flux at multiple time scales ranging from turbulent (minutes to hours) to daily to seasonal. The control parameters include wind stress, convective surface heat flux, shear magnitude, and thickness of the MI layer. In the LES, DCT occurs in discrete bursts during the night, exhibits high temporal variability within a burst, and modulates the mixed layer depth. At the daily time scale, turbulent heat flux generally increases with increasing wind stress, MI-layer shear, or nighttime convection. Convection is found to be important to mixing under weak wind, weak shear conditions. A parameterization for the daily averaged turbulent heat flux is developed from the LES suite to infer the variability of heat flux at the seasonal time scale. The LES-based parameterized heat flux, which takes into account the effects of all control parameters, exhibits a seasonal variability that is similar to the observed heat flux from theχ-pods.

 

Direct numerical simulations are performed to investigate a stratified shear layer at high Reynolds number ( $Re$ ) in a study where the Richardson number ( $Ri$ ) is varied among cases. Unlike previous work on a two-layer configuration in which the shear layer resides between two layers with constant density, an unbounded fluid with uniform stratification is considered here. The evolution of the shear layer includes a primary Kelvin–Helmholtz shear instability followed by a wide range of secondary shear and convective instabilities, similar to the two-layer configuration. During transition to turbulence, the shear layers at low $Ri$ exhibit a period of thickness contraction (not observed at lower $Re$ ) when the momentum and buoyancy fluxes are counter-gradient. The behaviour in the turbulent regime is significantly different from the case with a two-layer density profile. The transition layers, which are zones with elevated shear and stratification that form at the shear-layer edges, are stronger and also able to support a significant internal wave flux. After the shear layer becomes turbulent, mixing in the transition layers is shown to be more efficient than that which develops in the centre of the shear layer. Overall, the cumulative mixing efficiency ( $E^C$ ) is larger than the often assumed value of 1/6. Also, $E^C$ is found to be smaller than that in the two-layer configuration at moderate Ri . It is relatively less sensitive to background stratification, exhibiting little variation for $0.08 \leqslant Ri \leqslant 0.2$ . The dependence of mixing efficiency on buoyancy Reynolds number during the turbulence phase is qualitatively similar to homogeneous sheared turbulence.