skip to main content

Title: A Regime Diagram for Internal Lee Waves in Coastal Plain Estuaries
Abstract Using an idealized channel representative of a coastal plain estuary, we conducted numerical simulations to investigate the generation of internal lee waves by lateral circulation. It is shown that the lee waves can be generated across all salinity regimes in an estuary. Since the lateral currents are usually subcritical with respect to the lowest mode, mode-2 lee waves are most prevalent but a hydraulic jump may develop during the transition to subcritical flows in the deep channel, producing high energy dissipation and strong mixing. Unlike flows over a sill, stratified water in the deep channel may become stagnant such that a mode-1 depression wave can form higher up in the water column. With the lee wave Froude number above 1 and the intrinsic wave frequency between the inertial and buoyancy frequency, the lee waves generated in coastal plain estuaries are nonlinear waves with the wave amplitude Δ h scaling approximately with , where V is the maximum lateral flow velocity and is the buoyancy frequency. The model results are summarized using the estuarine classification diagram based on the freshwater Froude number Fr f and the mixing parameter M . The Δ h decreases with increasing Fr f as stronger stratification suppresses waves, and no internal waves are generated at large Fr f . The Δ h initially increases with increasing M as the lateral flows become stronger with stronger tidal currents, but decreases or saturates to a certain amplitude as M further increases. This modeling study suggests that lee waves can be generated over a wide range of estuarine conditions.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Physical Oceanography
Page Range / eLocation ID:
3049 to 3064
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Recent mooring observations at a cross-channel section in Chesapeake Bay showed that internal solitary waves regularly appeared during certain phases of a tidal cycle and propagated from the deep channel to the shallow shoal. It was hypothesized that these waves resulted from the nonlinear steepening of internal lee waves generated by lateral currents over channel-shoal topography. In this study numerical modeling is conducted to investigate the interaction between lateral circulation and cross-channel topography and discern the generation mechanism of the internal lee waves. During ebb tides, lateral bottom Ekman forcing drives a counterclockwise (looking into estuary) lateral circulation, with strong currents advecting stratified water over the western flank of the deep channel and producing large isopycnal displacements. When the lateral flow becomes supercritical with respect to mode-2 internal waves, a mode-2 internal lee wave is generated on the flank of the deep channel and subsequently propagates onto the western shoal. When the bottom lateral flow becomes near-critical or supercritical with respect to mode-1 internal waves, the lee wave evolves into an internal hydraulic jump. On the shallow shoal, the lee waves or jumps evolve into internal bores of elevation.

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

    The fate of discharges from small rivers and mountainous streams are little studied relative to the dynamics of large river plumes. Flows from small watersheds are episodic, forming transient low‐salinity surface plumes that vary tidally due to interaction of outflow inertia with buoyancy and ambient tidal currents. An implementation of the Regional Ocean Modeling System (ROMS v3.0), a three‐dimensional, free‐surface, terrain‐following numerical model, is applied to small river outflows (≤10 m3 s−1) entering a tidal ocean, where they form small plumes of scale ~103 m or smaller. Analysis of the momentum balance points to three distinct zones: (i) an inertia‐driven near field, where advection, pressure gradient, and vertical stress divergence control the plume dynamics, (ii) a buoyancy‐driven midfield, where pressure gradient, lateral stress divergence, and rotational accelerations are dominant, and (iii) an advective far field, where local accelerations induced by ambient tidal currents determine the fate of the river/estuarine discharge. The response of these small tidal plumes to different buoyancy forcing and outflow rate is explored. With increasing buoyancy, plumes change from narrow/elongated, bottom‐attached flows to radially expanding, surface‐layer flows. Weak outflow promotes stratification within the plume layer, and with stronger outflow, the plume layer becomes thinner and well‐mixed. When compared with prototypical large‐river plumes in which Coriolis effects are important, (i) in these small plumes, there is no bulge and no coastal buoyancy current, that is, shore contact is mostly absent, and (ii) the plume is strongly influenced by ambient tidal currents, forming a tidal plume that is deflected upcoast/downcoast from the river mouth.

    more » « less
  4. Observations on the lee of a topographic ridge show that the turbulence kinetic energy (TKE) dissipation rate due to shear instabilities is three orders of magnitude higher than the typical value in the open ocean. Laboratory-scale studies at low Reynolds number suggest that high turbulent dissipation occurs primarily within the core region of shear instabilities. However, field-scale studies indicate that high turbulence is mainly populated along the braids of shear instabilities. In this study, a high-resolution, resolving the Ozmidov-scale, non-hydrostatic model with Large Eddy Simulation (LES) turbulent closure is applied to investigate dominant mechanisms that control the spatial and temporal scales of shear instabilities and resulting mixing in stratified shear flow at high Reynolds number. The simulated density variance dissipation rate is elevated in the cusp-like bands of shear instabilities with a specific period, consistent with the acoustic backscatter taken by shipboard echo sounder. The vertical length scale of each cusp-like band is nearly half of the vertical length scale of the internal lee wave. However, it is consistent with instabilities originating from a shear layer based on linear stability theory. The model results indicate that the length scale and/or the period of shear instabilities are the key parameters to the mixing enhancement that increases with lateral Froude number Fr L , i.e. stronger shear and/or steeper ridge. 
    more » « less
  5. 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