Structure and intensity of estuarine exchange flow depend significantly on the eddy viscosity A_v profile which is dynamically linked to various forces (e.g., gravitational, tidal, wind-driven). The impact of winds on the exchange flow is complex due to its direct (local and remote changes in shear and density stratification) and indirect (modifications to A_v profiles) contributions. This study aims (i) to include wind entrainment effects in the tidally averaged A_v parameterization; (ii) to develop an analytical one-dimensional model for the wind driven exchange flow by using this novel parameterization and assess the tidally averaged dynamics over a relevant physical parameter-space, subdomains of which have not yet been explored numerically. This one-dimensional model is based on a balance between frictional forces and pressure gradient, calibrated with a tidally-resolving one-dimensional water-column model with second-moment closure. Structure and intensity of the resulting exchange flow profiles are analyzed with respect to three dimensionless parameters (the unsteadiness of boundary layer mixing U_n, scaled-directional wind stress W, and horizontal stratification S_i). While down-estuarine winds enhance the gravitational circulation, up-estuarine winds result in either a two-layer inverted circulation opposing the gravitational circulation, or a three-layer flow (favored by relatively strong S_i, weak W, and moderate U_n) that is up-estuarine at the surface with classical two-layer circulation underneath. Relative thickness of surface and bottom boundary layers affect both the intensity and the inflection depth of the exchange flow layers. Up-estuarine winds with W≳0.5 yield unstable stratification and reduce the exchange flow intensity with increasing W.

A new parameterization for the estuarine turbulent eddy viscosity coefficient is developed considering the influence of wind forcing and feedback between stratification and shear. The emerging tidally averaged eddy viscosity profile A_v^T is parameterized as parabolic under well-mixed conditions, and is composed of a skewed-Gaussian-like form for the upper layer, and a parabolic form for the bottom layer under stratified conditions. The precise shape of the profiles depends parametrically on the bottom boundary layer thickness, the bulk Richardson number, and the Wedderburn number. The parameterized A_v^T profiles show excellent agreement with profiles obtained from numerical models. To explore the importance of vertically varying A_v^T with regard to exchange processes, an analytical model is designed. This one-dimensional model is based on a balance between frictional forces and pressure gradient. The resulting exchange flow is analyzed over the relevant parameter space that is associated with horizontal and vertical stratification through the bulk Richardson number, and the bi-directional wind stress via the Wedderburn number. Down-estuary wind enhances the up-estuary flow near the bottom and down-estuary flow near the surface driving an exchange flow pattern typically associated with gravitational circulation. Up-estuary wind results in either a two-layer inverted circulation opposing the gravitational circulation, or a three-layer flow that is up-estuarine at the surface with classical two-layer circulation underneath. Three-layer flow emerges with a weak wind. With increasing runoff velocity, three-layer flow transitions to a single layer flow under weak stratification conditions.