Turbulent mixing in a strongly forced salt wedge estuary

Abstract

Turbulent mixing of salt is examined in a shallow salt wedge estuary with strong fluvial and tidal forcing. A numerical model of the Merrimack River estuary is used to quantify turbulent stress, shear production, and buoyancy flux. Little mixing occurs during flood tides despite strong velocities because bottom boundary layer turbulence is dislocated from stratification elevated in the water column. During ebbs, bottom salinity fronts form at a series of bathymetric transitions. At the fronts, near-bottom velocity and shear stress are low, but shear, stress, and buoyancy flux are elevated at the pycnocline. Internal shear layers provide the dominant source of mixing during the early ebb. Later in the ebb, the pycnocline broadens and moves down such that boundary layer turbulence dominates mixing. Mixing occurs primarily during ebbs, with internal shear mixing accounting for about 50% of the total buoyancy flux. Both the relative contribution of internal shear mixing and the mixing efficiency increase with discharge, with bulk mixing efficiencies between 0.02 and 0.07. Buoyancy fluxes in the estuary increase with discharge up to about 400 m3 s−1 above which a majority of the mixing occurs offshore. Observed buoyancy fluxes were more consistent with the k-ɛ turbulence closure than the Mellor-Yamada closure, and more total mixing occurred in the estuary with k-ɛ. Calculated buoyancy fluxes were sensitive to horizontal grid resolution, as a lower resolution grid yielded less integrated buoyancy flux in the estuary and exported lower salinity water but likely had greater numerical mixing.