skip to main content

Title: Turbulence, Sediment‐Induced Stratification, and Mixing Under Macrotidal Estuarine Conditions (Qiantang Estuary, China)

Time series of in situ measured velocity and suspended sediment concentration from Qiantang Estuary (China), and estimates of turbulence and sediment stratification parameters are presented. The data span a period of 9 days, and after phase averaged, they are used to explore spring‐neap tidal variations in flow, turbulence, and sediment stratification. A local balance between shear production, sediment‐induced buoyancy flux, and dissipation is found to hold during ebb for both neap and spring tides. During flood elevated turbulence dissipation rates are observed, attributed to nonlocal turbulence, most likely due to horizontal advection. Our results show that the effect of sediment stratification is successfully parameterized by adding the Monin‐Obukhov length scale to the classical logarithmic layer theory. Flood‐ebb asymmetry in bothRigandRfis observed, with higher values attained during flood due to the higher sediment concentrations and the corresponding weaker velocity shear found at these times.Rfis found to increase withRigwhenRig< 0.25, and it attains a maximum value of ~0.2. Our estimates, consistently with those from previous studies, fall slightly below the commonly used theoretical predictions. Sediment stratification contributes to the decay of turbulence, and weak mixing still exists under highRignumbers.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Oceans
Page Range / eLocation ID:
p. 4058-4077
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The salinity structure in an estuary is controlled by time‐dependent mixing processes. However, the locations and temporal variability of where significant mixing occurs is not well‐understood. Here we utilize a tracer variance approach to demonstrate the spatial and temporal structure of salinity mixing in the Hudson River Estuary. We run a 4‐month hydrodynamic simulation of the tides, currents, and salinity that captures the spring‐neap tidal variability as well as wind‐driven and freshwater flow events. On a spring‐neap time scale, salinity variance dissipation (mixing) occurs predominantly during the transition from neap to spring tides. On a tidal time scale, 60% of the salinity variance dissipation occurs during ebb tides and 40% during flood tides. Spatially, mixing during ebbs occurs primarily where lateral bottom salinity fronts intersect the bed at the transition from the main channel to adjacent shoals. During ebbs, these lateral fronts form seaward of constrictions located at multiple locations along the estuary. During floods, mixing is generated by a shear layer elevated in the water column at the top of the mixed bottom boundary layer, where variations in the along channel density gradients locally enhance the baroclinic pressure gradient leading to stronger vertical shear and more mixing. For both ebb and flood, the mixing occurs at the location of overlap of strong vertical stratification and eddy diffusivity, not at the maximum of either of those quantities. This understanding lends a new insight to the spatial and time dependence of the estuarine salinity structure.

    more » « less
  2. 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
  3. Abstract

    Dye‐release and numerical modeling experiments are conducted to investigate horizontal dispersion in a partially mixed estuary. Longitudinal dispersion of a dye patch shows marked flood‐ebb asymmetry in the first two tidal cycles after a dye release, with most of the dispersion occurring during ebb tides. There are large differences in the dispersion rate between spring and neap tides. Due to increased vertical mixing during spring tides, a dye patch quickly extends from the bottom to the surface and is exposed to the vertical shear in the entire water column, enhancing longitudinal dispersion. In contrast, most of the dye patch is limited to the bottom few meters during the neap tides. Although decreased vertical mixing during neap tides facilitates longitudinal dispersion, the vertical shear across the thin dye patch in the bottom layer is much weaker than the full water column shear in spring tides. Decreased shear leads to reduced longitudinal dispersion. After four tidal cycles, the second moment of the dye patch increases with time in the along‐channel direction at a power of between 2 and 3. The longitudinal dispersion rate varies as the four thirds power of the dye patch size, indicating scale‐dependent diffusion.

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