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.