ALBERT

Description: The temporal response of the length of a partially mixed estuary to changes in freshwater discharge Qf and tidal amplitude UT is studied using a 108-day time series collected along the length of the Hudson River estuary in the spring and summer of 2004 and a long-term (13.4 yr) record of Qf , UT, and near-surface salinity. When Qf was moderately high, the tidally averaged length of the estuary L5, here defined as the distance from the mouth to the up-estuary location where the vertically averaged salinity is 5 psu, fluctuated by more than 47 km over the spring–neap cycle, ranging from 28 to 〉75 km. During low flow periods, L5 varied very little over the spring–neap cycle and approached a steady length. The response is quantified and compared to predictions of a linearized model derived from the global estuarine salt balance. The model is forced by fluctuations in Qf and UT relative to average discharge Qo and tidal amplitude UTo and predicts the linear response time scale τ and the steady-state length Lo for average forcing. Two vertical mixing schemes are considered, in which 1) mixing is proportional to UT and 2) dependence of mixing on stratification is also parameterized. Based on least squares fits between L5 and estuary length predicted by the model, estimated τ varied by an order of magnitude from a period of high average discharge (Qo = 750 m3 s−1, τ = 4.2 days) to a period of low discharge (Qo = 170 m3 s−1, τ = 40.4 days). Over the range of observed discharge, Lo ∝ Qo−0.30±0.03, consistent with the theoretical scaling for an estuary whose landward salt flux is driven by vertical estuarine exchange circulation. Estimated τ was proportional to the discharge advection time scale (LoA/Qo, where A is the cross-sectional area of the estuary). However, τ was 3–4 times larger than the theoretical prediction. The model with stratification-dependent mixing predicted variations in L5 with higher skill than the model with mixing proportional to UT. This model provides insight into the time-dependent response of a partially stratified estuary to changes in forcing and explains the strong dependence of the amplitude of the spring–neap response on freshwater discharge. However, the utility of the linear model is limited because it assumes a uniform channel, and because the underlying dynamics are nonlinear, and the forcing Qf and UT can undergo large amplitude variations. River discharge, in particular, can vary by over an order of magnitude over time scales comparable to or shorter than the response time scale of the estuary.

Description: In most estuarine systems it is assumed that the dominant along-channel momentum balance is between the integrated pressure gradient and bed stress. Scaling the amplitude of the estuarine circulation based on this balance has been shown to have predictive skill. However, a number of authors recently highlighted important nonlinear processes that contribute to the subtidal dynamics at leading order. In this study, a previously validated numerical model of the Hudson River estuary is used to examine the forces driving the residual estuarine circulation and to test the predictive skill of two linear scaling relationships. Results demonstrate that the nonlinear advective acceleration terms contribute to the subtidal along-channel momentum balance at leading order. The contribution of these nonlinear terms is driven largely by secondary lateral flows. Under a range of forcing conditions in the model runs, the advective acceleration terms nearly always act in concert with the baroclinic pressure gradient, reinforcing the residual circulation. Despite the strong contribution of the nonlinear advective terms to the subtidal dynamical balance, a linear scaling accurately predicts the strength of the observed residual circulation in the model. However, this result is largely fortuitous, as this scaling does not account for two processes that are fundamental to the estuarine circulation. The skill of this scaling results because of the compensatory relationship between the contribution of the advective acceleration terms and the suppression of turbulence due to density stratification. Both of these processes, neither of which is accounted for in the linear scaling, increase the residual estuarine circulation but have an opposite dependence on tidal amplitude and, consequently, strength of stratification.

Description: The subtidal salt balance and the mechanisms driving the downgradient salt flux in the Hudson River estuary are investigated using measurements from a cross-channel mooring array of current meters, temperature and conductivity sensors, and cross-channel and along-estuary shipboard surveys obtained during the spring of 2002. Steady (subtidal) vertical shear dispersion, resulting from the estuarine exchange flow, was the dominant mechanism driving the downgradient salt flux, and varied by over an order of magnitude over the spring–neap cycle, with maximum values during neap tides and minimum values during spring tides. Corresponding longitudinal dispersion rates were as big as 2500 m2 s−1 during neap tides. The salinity intrusion was not in a steady balance during the study period. During spring tides, the oceanward advective salt flux resulting from the net outflow balanced the time rate of change of salt content landward of the study site, and salt was flushed out of the estuary. During neap tides, the landward steady shear dispersion salt flux exceeded the oceanward advective salt flux, and salt entered the estuary. Factor-of-4 variations in the salt content occurred at the spring–neap time scale and at the time scale of variations in the net outflow. On average, the salt flux resulting from tidal correlations between currents and salinity (tidal oscillatory salt flux) was an order of magnitude smaller than that resulting from steady shear dispersion. During neap tides, this flux was minimal (or slightly countergradient) and was due to correlations between tidal currents and vertical excursions of the halocline. During spring tides, the tidal oscillatory salt flux was driven primarily by oscillatory shear dispersion, with an associated longitudinal dispersion rate of about 130 m2 s−1.

Description: Observations at the mouth of the Fraser River (British Columbia, Canada) indicate an abrupt frontal transition between unstratified river outflow and a highly stratified river plume with differences in salinity greater than 25 psu across a few meters in the vertical direction and several hundred meters in the horizontal direction. The front roughly follows a natural break in the bathymetry, crossing the channel at an angle of approximately 45°, and is essentially stationary for a period of approximately 3.5 h centered on the low tide following the larger of two daily ebbs. The location of the front is coincident with observations of significantly supercritical internal Froude numbers at the front, based on velocities in the along-flow direction. This observation contradicts the one-dimensional theory, which indicates that the Froude number should be 1. However, because the front is oriented obliquely to the outflow, a coordinate system can be selected that is normal to the front and for which a critical Froude number of 1 is obtained. This indicates that a Froude angle, similar in application to a Mach angle for transonic flows, can be used to determine critical conditions when the front is oblique to the principal flow direction.

Description: A tidally and cross-sectionally averaged model based on the temporal evolution of the quasi-steady Hansen and Rattray equations is applied to simulate the salinity distribution and vertical exchange flow along the Hudson River estuary. The model achieves high skill at hindcasting salinity and residual velocity variation during a 110-day period in 2004 covering a wide range of river discharges and tidal forcing. The approach is based on an existing model framework that has been modified to improve model skill relative to observations. The external forcing has been modified to capture meteorological time-scale variability in salinity, stratification, and residual velocity due to sea level fluctuations at the open boundary and along-estuary wind stress. To reflect changes in vertical mixing due to stratification, the vertical mixing coefficients have been modified to use the bottom boundary layer height rather than the water depth as an effective mixing length scale. The boundary layer parameterization depends on the tidal amplitude and the local baroclinic pressure gradient through the longitudinal Richardson number, and improves the model response to spring–neap variability in tidal amplitude during periods of high river discharge. Finally, steady-state model solutions are evaluated for both the Hudson River and northern San Francisco Bay over a range of forcing conditions. Agreement between the model and scaling of equilibrium salinity intrusions lends confidence that the approach is transferable to other estuaries, despite significant differences in bathymetry. Discrepancies between the model results and observations at high river discharge are indicative of limits at which the formulation begins to fail, and where an alternative approach that captures two-layer dynamics would be more appropriate.