skip to main content


Title: Scaling the Mixing Efficiency of Sediment‐Stratified Turbulence
Abstract

The flux Richardson numberRf, also called the mixing efficiency of stratified turbulence, is important in determining geophysical flow phenomena such as ocean circulation and air‐sea transports. MeasuringRfin the field is usually difficult, thus parameterization ofRfbased on readily observed properties is essential. Here, estimates ofRfin a strongly turbulent, sediment‐stratified estuarine flow are obtained from measurements of covariance‐derived turbulent buoyancy fluxes (B) and spectrally fitted values of the dissipation rate of turbulent kinetic energy (ε). We test scalings forRfin terms of the buoyancy Reynolds number (Reb), the gradient Richardson number (Ri), and turbulent Froude number (Frt). Neither theReb‐based nor theRi‐based scheme is able to describe the observed variations inRf, but theFrt‐based parameterization works well. These findings support further use of theFrt‐ based parameterization in turbulent oceanic and estuarine environments.

 
more » « less
Award ID(s):
1830071
NSF-PAR ID:
10370191
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
13
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Turbulence parameters in the lower troposphere (up to ~4.5 km) are estimated from measurements of high-resolution and fast-response cold-wire temperature and Pitot tube velocity from sensors onboard DataHawk Unmanned Aerial Vehicles (UAVs) operated at the Shigaraki Middle and Upper atmosphere (MU) Observatory during two ShUREX (Shigaraki UAV Radar Experiment) campaigns in 2016 and 2017. The practical processing methods used for estimating turbulence kinetic energy dissipation rate ε and temperature structure function parameter C T 2 from one-dimensional wind and temperature frequency spectra are first described in detail. Both are based on the identification of inertial (−5/3) subranges in respective spectra. Using a formulation relating ε and C T 2 valid for Kolmogorov turbulence in steady state, the flux Richardson number R f and the mixing efficiency χ m are then estimated. The statistical analysis confirms the variability of R f and χ m around ~ 0.13 − 0.14 and ~ 0.16 − 0.17 , respectively, values close to the canonical values found from some earlier experimental and theoretical studies of both the atmosphere and the oceans. The relevance of the interpretation of the inertial subranges in terms of Kolmogorov turbulence is confirmed by assessing the consistency of additional parameters, the Ozmidov length scale L O , the buoyancy Reynolds number R e b , and the gradient Richardson number Ri. Finally, a case study is presented showing altitude differences between the peaks of N 2 , C T 2 and ε , suggesting turbulent stirring at the margin of a stable temperature gradient sheet. The possible contribution of this sheet and layer structure on clear air radar backscattering mechanisms is examined. 
    more » « less
  2. Abstract

    We performed laboratory experiments to investigate the influence of sand content on the dynamics of wave‐supported gravity flows in mud‐dominant environments. The experiments were carried out in an oscillatory water tunnel with a sediment bed of either 1% or 13% sand. Low and high energy regimes are differentiated based on a Stokes Reynolds numberReΔ ≈ 500. In the low energy regime, the sand fraction influences flow dynamics primarily through ripple formation; no ripples form in the 1% sand experiments, whereas ripples form in the 13% experiments that increase turbulence and the wave boundary layer thickness,δm. In the high energy regime, small ripples form in both the 1% and 13% sand experiments and we observe high near‐bed suspended sediment concentrations. The influence of stratification on the boundary layer flow is characterized in terms of the gradient Richardson numberRig. The flow is weakly stratified inside the boundary layer for all runs and critically stratified at or above the top of the boundary layer. In the lower energy regime, the sand content reduces the relative influence of stratification in the boundary layer, shifting the elevation of critical stratification,LB, from approximately 1.3δmto 2.5δmin the 1% and 13% experiments, respectively. In both sets of experimentsLB ≈ δmat the strongest wave energy, indicating a transition to strongly stratified dynamics.

     
    more » « less
  3. Abstract

    We formulate an expression for the turbulent kinetic energy dissipation rate,ϵ, associated with shear‐generated turbulence in terms of quantities in the ocean or atmosphere that, depending on the situation, may be measurable or resolved in models. The expression depends on the turbulent vertical length scale,v, the inverse time scaleN, and the Richardson numberRi = N2/S2, whereSis the vertical shear, withvscaled in a way consistent with theories and observations of stratified turbulence. Unlike previous studies, the focus is not so much on the functional form ofRi, but the vertical variation of the length scalev. Using data from two ∼7‐day time series in the western equatorial Pacific, the scaling is compared with the observedϵ. The scaling works well with the estimatedϵcapturing the differences in amplitude and vertical distribution of the observedϵbetween the two times series. Much of those differences are attributable to changes in the vertical distribution of the length scalev, and in particular the associated turbulent velocity scale,ut. We relateutto a measure of the fine‐scale variations in velocity,. Our study highlights the need to consider the length scale and its estimation in environmental flows. The implications for the vertical variation of the associated turbulent diffusivity are discussed.

     
    more » « less
  4. 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. 
    more » « less
  5. Abstract

    Anisotropic turbulence is ubiquitous in atmospheric and oceanic boundary layers due to differences in energy injection mechanisms. Unlike mechanical production that injects energy in the streamwise velocity component, buoyancy affects only the vertical velocity component. This anisotropy in energy sources, quantified by the flux Richardson numberRif, is compensated by a “return to isotropy” (RTI) tendency of turbulent flows. Describing RTI in Reynolds‐averaged models and across scales continues to be a challenge in stratified turbulent flows. Using phenomenological models for spectral energy transfers, the necessary conditions for which the widely‐used Rotta model captures RTI across variousRifand eddy sizes are discussed for the first time. This work unravels adjustments to the Rotta constant, withRifand scale, necessary to obtain consistency between RTI models and the measured properties of the atmospheric surface layer for planar‐homogeneous and stationary flows in the absence of subsidence. A range ofRifand eddy sizes where the usage of a conventional Rotta model is prohibited is also found. Those adjustments lay the groundwork for new closure schemes.

     
    more » « less