The balance of salinity variance in a partially stratified estuary: implications for exchange flow, mixing, and stratification

Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in he balance of salinity variance in a partially stratified estuary: Implications for exchange flow, mixing, and stratification. Journal of Physical Oceanography, 48(12), (2018) 2887-2899., doi: 10.1175/JPO-D-18-0032.1.Salinity variance dissipation is related to exchange flow through the salinity variance balance equation, and meanwhile its magnitude is also proportional to the turbulence production and stratification inside the estuary. As river flow increases, estuarine volume-integrated salinity variance dissipation increases owing to more variance input from the open boundaries driven by exchange flow and river flow. This corresponds to the increased efficient conversion of turbulence production to salinity variance dissipation due to the intensified stratification with higher river flow. Through the spring–neap cycle, the temporal variation of salinity variance dissipation is more dependent on stratification than turbulence production, so it reaches its maximum during the transition from neap to spring tides. During most of the transition time from spring to neap tides, the advective input of salinity variance from the open boundaries is larger than dissipation, resulting in the net increase of variance, which is mainly expressed as vertical variance, that is, stratification. The intensified stratification in turn increases salinity variance dissipation. During neap tides, a large amount of enhanced salinity variance dissipation is induced by the internal shear stress near the halocline. During most of the transition time from neap to spring tides, dissipation becomes larger than the advective input, so salinity variance decreases and the stratification is destroyed.TW was supported by the National Key R&D Program of China (Grant 2017YFA0604104), National Natural Science Foundation of China (Grant 41706002), Natural Science Foundation of Jiangsu Province (Grant BK20170864), and MEL Visiting Fellowship (MELRS1617). WRG was supported by NSF Grant OCE 1736539. Part of this work is finished during TW’s visit in MEL and WHOI. We would like to acknowledge John Warner for providing the codes of the Hudson estuary model, and Parker MacCready, the editor, and two reviewers for their insightful suggestions on improving the manuscript.2019-06-0

Estuarine frontogenesis

Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 546–561, doi:10.1175/JPO-D-14-0082.1.Model studies and observations in the Hudson River estuary indicate that frontogenesis occurs as a result of topographic forcing. Bottom fronts form just downstream of lateral constrictions, where the width of the estuary increases in the down-estuary (i.e., seaward) direction. The front forms during the last several hours of the ebb, when the combination of adverse pressure gradient in the expansion and baroclinicity cause a stagnation of near-bottom velocity. Frontogenesis is observed in two dynamical regimes: one in which the front develops at a transition from subcritical to supercritical flow and the other in which the flow is everywhere supercritical. The supercritical front formation appears to be associated with lateral flow separation. Both types of fronts are three-dimensional, with strong lateral gradients along the flanks of the channel. During spring tide conditions, the fronts dissipate during the flood, whereas during neap tides the fronts are advected landward during the flood. The zone of enhanced density gradient initiates frontogenesis at multiple constrictions along the estuary as it propagates landward more than 60 km during several days of neap tides. Frontogenesis and frontal propagation may thus be essential elements of the spring-to-neap transition to stratified conditions in partially mixed estuaries.Support for this research was provided by NSF Grant OCE 0926427.2015-08-0

Sediment transport time scales and trapping efficiency in a tidal river

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Earth Surface 122 (2017): 2042–2063, doi:10.1002/2017JF004337.Observations and a numerical model are used to characterize sediment transport in the tidal Hudson River. A sediment budget over 11 years including major discharge events indicates the tidal fresh region traps about 40% of the sediment input from the watershed. Sediment input scales with the river discharge cubed, while seaward transport in the tidal river scales linearly, so the tidal river accumulates sediment during the highest discharge events. Sediment pulses associated with discharge events dissipate moving seaward and lag the advection speed of the river by a factor of 1.5 to 3. Idealized model simulations with a range of discharge and settling velocity were used to evaluate the trapping efficiency, transport rate, and mean age of sediment input from the watershed. The seaward transport of suspended sediment scales linearly with discharge but lags the river velocity by a factor that is linear with settling velocity. The lag factor is 30–40 times the settling velocity (mm s−1), so transport speeds vary by orders of magnitude from clay (0.01 mm s−1) to coarse silt (1 mm s−1). Deposition along the tidal river depends strongly on settling velocity, and a simple advection-reaction equation represents the loss due to settling on depositional shoals. The long-term discharge record is used to represent statistically the distribution of transport times, and time scales for settling velocities of 0.1 mm s−1 and 1 mm s−1 range from several months to several years for transport through the tidal river and several years to several decades through the estuary.Hudson River Foundation Grant Number: 004/13A; National Science Foundation Grant Number: 13251362018-05-0

Response to channel deepening of the salinity intrusion, estuarine circulation, and stratification in an urbanized estuary

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research-Oceans 124(7), (2019): 4784-4802, doi: 10.1029/2019JC015006.Modifications for navigation since the late 1800s have increased channel depth (H) in the lower Hudson River estuary by 10–30%, and at the mouth the depth has more than doubled. Observations along the lower estuary show that both salinity and stratification have increased over the past century. Model results comparing predredging bathymetry from the 1860s with modern conditions indicate an increase in the salinity intrusion of about 30%, which is roughly consistent with the H5/3 scaling expected from theory for salt flux dominated by steady exchange. While modifications including a recent deepening project have been concentrated near the mouth, the changes increase salinity and threaten drinking water supplies more than 100 km landward. The deepening has not changed the responses to river discharge (Qr) of the salinity intrusion (~Qr−1/3) or mean stratification (Qr2/3). Surprisingly, the increase in salinity intrusion with channel deepening results in almost no change in the estuarine circulation. This contrasts sharply with local scaling based on local dynamics of an H2 dependence, but it is consistent with a steady state salt balance that allows scaling of the estuarine circulation based on external forcing factors and is independent of depth. In contrast, the observed and modeled increases in stratification are opposite of expectations from the steady state balance, which could be due to reduction in mixing with loss of shallow subtidal regions. Overall, the mean shift in estuarine parameter space due to channel deepening has been modest compared with the monthly‐to‐seasonal variability due to tides and river discharge.Funding was provided by NSF Coastal SEES (OCE 1325136). Data supporting this study are posted to Zenodo (https://doi.org/10.5281/zenodo.2551285) or are available by contacting the author.2019-12-0

The role of advection, straining, and mixing on the tidal variability of estuarine stratification

Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 855–868, doi:10.1175/JPO-D-10-05010.1.Data from the Hudson River estuary demonstrate that the tidal variations in vertical salinity stratification are not consistent with the patterns associated with along-channel tidal straining. These observations result from three additional processes not accounted for in the traditional tidal straining model: 1) along-channel and 2) lateral advection of horizontal gradients in the vertical salinity gradient and 3) tidal asymmetries in the strength of vertical mixing. As a result, cross-sectionally averaged values of the vertical salinity gradient are shown to increase during the flood tide and decrease during the ebb. Only over a limited portion of the cross section does the observed stratification increase during the ebb and decrease during the flood. These observations highlight the three-dimensional nature of estuarine flows and demonstrate that lateral circulation provides an alternate mechanism that allows for the exchange of materials between surface and bottom waters, even when direct turbulent mixing through the pycnocline is prohibited by strong stratification.The funding for this research was obtained from NSF Grant OCE-08-25226.2012-11-0

The role of morphology and wave-current interaction at tidal inlets : an idealized modeling analysis

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 119 (2014): 8818–8837, doi:10.1002/2014JC010191.The outflowing currents from tidal inlets are influenced both by the morphology of the ebb-tide shoal and interaction with incident surface gravity waves. Likewise, the propagation and breaking of incident waves are affected by the morphology and the strength and structure of the outflowing current. The 3-D Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system is applied to numerically analyze the interaction between currents, waves, and bathymetry in idealized inlet configurations. The bathymetry is found to be a dominant controlling variable. In the absence of an ebb shoal and with weak wave forcing, a narrow outflow jet extends seaward with little lateral spreading. The presence of an ebb-tide shoal produces significant pressure gradients in the region of the outflow, resulting in enhanced lateral spreading of the jet. Incident waves cause lateral spreading and limit the seaward extent of the jet, due both to conversion of wave momentum flux and enhanced bottom friction. The interaction between the vorticity of the outflow jet and the wave stokes drift is also an important driving force for the lateral spreading of the plume. For weak outflows, the outflow jet is actually enhanced by strong waves when there is a channel across the bar, due to the “return current” effect. For both strong and weak outflows, waves increase the alongshore transport in both directions from the inlet due to the wave-induced setup over the ebb shoal. Wave breaking is more influenced by the topography of the ebb shoal than by wave-current interaction, although strong outflows show intensified breaking at the head of the main channel.We are grateful to the Career Training Interexchange program that facilitated the training period of Maitane Olabarrieta within the USGS. Maitane Olabarrieta also acknowledges funding from the “Cantabria Campus International Augusto Gonzalez Linares Program.”WRG was supported by ONR grant N00014-13-1–0368.2015-06-2

Total exchange flow, entrainment, and diffusive salt flux in estuaries

Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 1205-1220, doi:10.1175/JPO-D-16-0258.1.The linkage among total exchange flow, entrainment, and diffusive salt flux in estuaries is derived analytically using salinity coordinates, revealing the simple but important relationship between total exchange flow and mixing. Mixing is defined and quantified in this paper as the dissipation of salinity variance. The method uses the conservation of volume and salt to quantify and distinguish the diahaline transport of volume (i.e., entrainment) and diahaline diffusive salt flux. A numerical model of the Hudson estuary is used as an example of the application of the method in a realistic estuary with a persistent but temporally variable exchange flow. A notable finding of this analysis is that the total exchange flow and diahaline salt flux are out of phase with respect to the spring–neap cycle. Total exchange flow reaches its maximum near minimum neap tide, but diahaline salt transport reaches its maximum during the maximum spring tide. This phase shift explains the strong temporal variation of stratification and estuarine salt content through the spring–neap cycle. In addition to quantifying temporal variation, the method reveals the spatial variation of total exchange flow, entrainment, and diffusive salt flux through the estuary. For instance, the analysis of the Hudson estuary indicates that diffusive salt flux is intensified in the wider cross sections. The method also provides a simple means of quantifying numerical mixing in ocean models because it provides an estimate of the total dissipation of salinity variance, which is the sum of mixing due to the turbulence closure and numerical mixing.T. Wang was supported by the Open Research Fund of State Key Laboratory of Estuarine and Coastal Research (Grant SKLEC-KF201509), the Fundamental Research Funds for the Central Universities (Grant 2017B03514), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant XDA11010203). W. R. Geyer was supported by NSF Grant OCE 0926427 and ONR Grant N00014-16-1-2948. P. MacCready was supported by NSF Grant OCE-1634148.2017-09-1

Estuarine exchange flow is related to mixing through the salinity variance budget

Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1375-1384, doi:10.1175/JPO-D-17-0266.1.The relationship between net mixing and the estuarine exchange flow may be quantified using a salinity variance budget. Here “mixing” is defined as the rate of destruction of volume-integrated salinity variance, and the exchange flow is quantified using the total exchange flow. These concepts are explored using an idealized 3D model estuary. It is shown that in steady state (e.g., averaging over the spring–neap cycle) the volume-integrated mixing is approximately given by Mixing ≅ SinSoutQr, where Sin and Sout are the representative salinities of in- and outflowing layers at the mouth and Qr is the river volume flux. This relationship provides an extension of the familiar Knudsen relation, in which the exchange flow is diagnosed based on knowledge of these same three quantities, quantitatively linking mixing to the exchange flow.The work was supported by the National Science Foundation through Grants OCE-1736242 to PM and OCE-1736539 to WRG and by the German Research Foundation through Grants TRR 181 and GRK 2000 to HB

Tidal and spring-neap variations in horizontal dispersion in a partially mixed estuary

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C07023, doi:10.1029/2007JC004644.A sequence of dye releases in the Hudson River estuary provide a quantitative assessment of horizontal dispersion in a partially mixed estuary. Dye was released in the bottom boundary layer on 4 separate occasions, with varying tidal phase and spring-neap conditions. The three-dimensional distribution of dye was monitored by two vessels with in situ, profiling fluorometers. The three-dimensional spreading of the dye was estimated by calculating the time derivative of the second moment of the dye in the along-estuary, cross-estuary and vertical directions. The average along-estuary dispersion rate was about 100 m2/s, but maximum rates up to 700 m2/s occurred during ebb tides, and minimum rates occurred during flood. Vertical shear dispersion was the principal mechanism during neap tides, but transverse shear dispersion became more important during springs. Suppression of mixing across the pycnocline limited the vertical extent of the patch in all but the maximum spring-tide conditions, with vertical diffusivities in the pycnocline estimated at 4 × 10−5 m2/s during neaps. The limited vertical extent of the dye patch limited the dispersion of the dye relative to the overall estuarine dispersion rate, which was an order of magnitude greater than that of the dye. This study indicates that the effective dispersion of waterborne material in an estuary depends sensitively on its vertical distribution as well as the phase of the spring-neap cycle.This research was supported by National Science Foundation Grant OCE04-52054 (W. Geyer), OCE00-99310 (R. Houghton), and OCE00-95913 (R. Chant)

Broadband acoustic quantification of stratified turbulence

Author Posting. © Acoustical Society of America, 2013. This article is posted here by permission of Acoustical Society of America for personal use, not for redistribution. The definitive version was published in Journal of the Acoustical Society of America 134 (2013): 40-54, doi:10.1121/1.4807780.High-frequency broadband acoustic scattering techniques have enabled the remote, high-resolution imaging and quantification of highly salt-stratified turbulence in an estuary. Turbulent salinity spectra in the stratified shear layer have been measured acoustically and by in situ turbulence sensors. The acoustic frequencies used span 120–600 kHz, which, for the highly stratified and dynamic estuarine environment, correspond to wavenumbers in the viscous-convective subrange (500–2500 m−1). The acoustically measured spectral levels are in close agreement with spectral levels measured with closely co-located micro-conductivity probes. The acoustically measured spectral shapes allow discrimination between scattering dominated by turbulent salinity microstructure and suspended sediments or swim-bladdered fish, the two primary sources of scattering observed in the estuary in addition to turbulent salinity microstructure. The direct comparison of salinity spectra inferred acoustically and by the in situ turbulence sensors provides a test of both the acoustic scattering model and the quantitative skill of acoustical remote sensing of turbulence dissipation in a strongly sheared and salt-stratified estuary.This work was supported by NSF grant OCE- 0824871, ONR grant N00014-0810495, and WHOI internal funds