ALBERT

Description: A series of dye releases in the Hudson River estuary elucidated diapycnal mixing rates and temporal variability over tidal and fortnightly time scales. Dye was injected in the bottom boundary layer for each of four releases during different phases of the tide and of the spring–neap cycle. Diapycnal mixing occurs primarily through entrainment that is driven by shear production in the bottom boundary layer. On flood the dye extended vertically through the bottom mixed layer, and its concentration decreased abruptly near the base of the pycnocline, usually at a height corresponding to a velocity maximum. Boundary layer growth is consistent with a one-dimensional, stress-driven entrainment model. A model was developed for the vertical structure of the vertical eddy viscosity in the flood tide boundary layer that is proportional to u2*/N∞, where u* and N∞ are the bottom friction velocity and buoyancy frequency above the boundary layer. The model also predicts that the buoyancy flux averaged over the bottom boundary layer is equal to 0.06N∞u2* or, based on the structure of the boundary layer equal to 0.1NBLu2*, where NBL is the buoyancy frequency across the flood-tide boundary layer. Estimates of shear production and buoyancy flux indicate that the flux Richardson number in the flood-tide boundary layer is 0.1–0.18, consistent with the model indicating that the flux Richardson number is between 0.1 and 0.14. During ebb, the boundary layer was more stratified, and its vertical extent was not as sharply delineated as in the flood. During neap tide the rate of mixing during ebb was significantly weaker than on flood, owing to reduced bottom stress and stabilization by stratification. As tidal amplitude increased ebb mixing increased and more closely resembled the boundary layer entrainment process observed during the flood. Tidal straining modestly increased the entrainment rate during the flood, and it restratified the boundary layer and inhibited mixing during the ebb.

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.