Difference between revisions of "Estuarine circulation"
Dronkers J (talk  contribs) 
Dronkers J (talk  contribs) 

Line 1:  Line 1:  
{{Definitiontitle=Estuarine circulation  {{Definitiontitle=Estuarine circulation  
definition= Residual flow pattern in an estuary induced by the density difference between seawater and river water. }}  definition= Residual flow pattern in an estuary induced by the density difference between seawater and river water. }}  
+  
+  
+  Seawater has a higher density than river water, due to its high salinity. The density difference is of the order of 3%. In estuaries, where seawater and river water meet, seawater will tend to intrude beneath the river flow. However, mixing will also occur between seawater and river water, whereby seawater is returned to the sea by the river flow via the upper water layer. This allows 'new' seawater to enter the estuary via the lower water layers. In this way, a circulation of seawater is created, which is directed inland along the bottom and seaward along the surface. This circulation is called estuarine circulation or gravitational circulation. It is driven by the longitudinal density gradient along the estuary. In addition to the density differences along the estuary, the circulation depends on the mixing of seawater with river water. Mixing is the result of turbulence, which in estuaries is mainly caused by strong tidal flows. Most of the mixing takes place vertically in the water column, due to turbulence generated by tidal flow along the estuarine bottom. The different processes that drive estuarine circulation are discussed in this article. In estuaries with little tide there is hardly any mixing of seawater with river water. In that case a more or less stationary salt wedge forms along the bottom; this case is discussed in the article [[Salt wedge estuaries]].  
+  
Line 9:  Line 13:  
The classical theory of estuarine circulation was first developed by Hansen and Rattray in 1965 <ref>Hansen, D.V. and Rattray M. (1965) Gravitational circulation in straits and estuaries. J. Mar. Res. 23: 104–22</ref>, after earlier investigations by D. Pritchard in the St. James river <ref>Pritchard, D.W. (1952) Salinity distribution and circulation in the Chesapeake Bay estuaries system. J. Mar. Res. 11: 106–23</ref>. Many other authors have subsequently contributed to further verification and refinement<ref>Geyer, W.R., Trowbridge. J.H. and Bowen, M.M. (2000) The dynamics of a partially mixed estuary. J. Phys. Oceanogr. 30: 2035–48</ref><ref>MacCready, P. (2004) Toward a unified theory of tidallyaveraged estuarine salinity structure. Estuaries 27: 561–70</ref>.  The classical theory of estuarine circulation was first developed by Hansen and Rattray in 1965 <ref>Hansen, D.V. and Rattray M. (1965) Gravitational circulation in straits and estuaries. J. Mar. Res. 23: 104–22</ref>, after earlier investigations by D. Pritchard in the St. James river <ref>Pritchard, D.W. (1952) Salinity distribution and circulation in the Chesapeake Bay estuaries system. J. Mar. Res. 11: 106–23</ref>. Many other authors have subsequently contributed to further verification and refinement<ref>Geyer, W.R., Trowbridge. J.H. and Bowen, M.M. (2000) The dynamics of a partially mixed estuary. J. Phys. Oceanogr. 30: 2035–48</ref><ref>MacCready, P. (2004) Toward a unified theory of tidallyaveraged estuarine salinity structure. Estuaries 27: 561–70</ref>.  
−  In the classical theory of estuarine circulation tidal motion  +  In the classical theory of estuarine circulation, advection of seawater up and down the estuary by the tidal motion and associated mixing processes are not considered explicitly. These models consider a stationary salinity field in a frame that moves cyclically up and down the estuary and a tideaveraged coefficient for the mixing processes. The effect of tidal motion is further restricted to turbulence generation at the bed, which enhances momentum dissipation and mixing in and across the bottom boundary layer. 
==Notation==  ==Notation==  
−  In the following basic notions of the classical theory are presented. We consider a straight prismatic estuary (uniform depth <math>h</math> and width <math>b</math> ); the width is small enough for neglecting the influence of earth rotation (Rossby radius <math>U/fb << 1</math> ). We use the following notations (see also figure 1):  +  In the following, some basic notions of the classical theory are presented. We consider a straight prismatic estuary (uniform depth <math>h</math> and width <math>b</math> ); the width is small enough for neglecting the influence of earth's rotation (Rossby radius <math>U/fb << 1</math> ). We use the following notations (see also figure 1): 
* <math>x</math> is the coordinate along the estuarine axis (<math>x=0</math> at the sea boundary, <math>x=L</math> at the landward boundary of salt intrusion);  * <math>x</math> is the coordinate along the estuarine axis (<math>x=0</math> at the sea boundary, <math>x=L</math> at the landward boundary of salt intrusion);  
−  * <math>u</math> is the velocity in the <math>x</math>direction, <math>u(z) = u_R + u_C(z) + U \cos \omega t,</math> where <math>u_R</math> is the velocity related to river runoff (river discharge divided by estuarine crosssection), <math>u_C(z)</math> is the estuarine circulation, <math>U</math> is the tidal velocity amplitude and <math>\omega = 2 \pi / T</math> is the tidal radial frequency;  +  * <math>u</math> is the velocity in the <math>x</math>direction, <math>u(z) = u_R + u_C(z) + U \cos \omega t,</math> where <math>u_R</math> is the velocity related to river runoff (river discharge divided by estuarine crosssection), <math>u_C(z)</math> is the estuarine circulation, <math>U</math> is the depthaveraged tidal velocity amplitude and <math>\omega = 2 \pi / T</math> is the tidal radial frequency; 
* <math>S=s/s_0</math> is the ratio of the local salinity <math>s</math> and the seawater salinity <math>s_0 :</math> <math>S</math> is related to the density <math>\rho</math> by the relation <math>\rho=\rho_0(1+\beta s_0 S),</math> where <math>\rho_0</math> is the fresh water density and the expansion coefficient <math>\beta=7.7 \times 10^{4};</math>  * <math>S=s/s_0</math> is the ratio of the local salinity <math>s</math> and the seawater salinity <math>s_0 :</math> <math>S</math> is related to the density <math>\rho</math> by the relation <math>\rho=\rho_0(1+\beta s_0 S),</math> where <math>\rho_0</math> is the fresh water density and the expansion coefficient <math>\beta=7.7 \times 10^{4};</math>  
−  * <math> p(z) = g  +  * <math> p(z) = g \int_z^\eta \rho dz</math> is the hydrostatic pressure; <math>\eta</math> is the surface elevation relative to a horizontal plane and <math>z</math> is the depth relative to the surface; 
* <math>K_M</math> is the eddy viscosity (vertical momentum diffusion coefficient <math>[m^2/s] </math> ) and <math>K_S</math> the eddy diffusivity (vertical salinity diffusion coefficient).  * <math>K_M</math> is the eddy viscosity (vertical momentum diffusion coefficient <math>[m^2/s] </math> ) and <math>K_S</math> the eddy diffusivity (vertical salinity diffusion coefficient).  
Line 25:  Line 29:  
The following additional assumptions are made:  The following additional assumptions are made:  
−  
* The diffusion coefficients <math>K_M, K_S</math> do not depend on the spatial and temporal coordinates <math>x,z,t</math> ; this means that the classical model is realistic only for continuously stratified estuaries (no layered salinity structure; such estuaries are often called "partially mixed", see figure 1);  * The diffusion coefficients <math>K_M, K_S</math> do not depend on the spatial and temporal coordinates <math>x,z,t</math> ; this means that the classical model is realistic only for continuously stratified estuaries (no layered salinity structure; such estuaries are often called "partially mixed", see figure 1);  
* The river flow velocity <math>u_R</math> and the tidal velocity amplitude <math>U</math> do not depend on the spatial coordinates <math>x,z</math> ;  * The river flow velocity <math>u_R</math> and the tidal velocity amplitude <math>U</math> do not depend on the spatial coordinates <math>x,z</math> ;  
Line 33:  Line 36:  
* Side wall effects are neglected, implying uniformity of all variables in lateral <math>y</math>direction;  * Side wall effects are neglected, implying uniformity of all variables in lateral <math>y</math>direction;  
* Longitudinal gradients in the velocity <math>u</math> are neglected;  * Longitudinal gradients in the velocity <math>u</math> are neglected;  
−  * The longitudinal salinity gradient is the same at all depths (<math>\  +  * The longitudinal salinity gradient is the same at all depths (<math>\large\frac{\partial S}{\partial x}\normalsize=\large\frac{\partial \overline S}{\partial x}\normalsize ).</math> 
−  
−  
−  
==Model formulation==  ==Model formulation==  
Line 42:  Line 42:  
[[Image:VelSalProfilesNew.jpgthumb300pxrightFigure 2: Vertical profiles of estuarine circulation and salinity according to the classical theory (Eq. 4).]]  [[Image:VelSalProfilesNew.jpgthumb300pxrightFigure 2: Vertical profiles of estuarine circulation and salinity according to the classical theory (Eq. 4).]]  
−  We consider the influence of constant river flow <math>u_R</math> on the salinity distribution in a frame moving with the velocity <math>u_R + U \cos \omega t</math> ; we thus have <math>u=u_C(z)  +  We consider the influence of a constant river flow <math>u_R</math> on the salinity distribution in a frame moving with the velocity <math>u_R + U \cos \omega t</math>. In the moving frame the velocity <math>u</math> is timeindependent and depends only on the vertical coordinate <math>z</math>; we thus have <math>u=u_C(z)</math>. The only driving force acting on the fluid in the moving frame is the pressure gradient <math>\partial p / \partial x </math>, which is a linearly increasing function of depth (<math>z</math>) below mean water level: 
−  <math>\  +  <math>\large\frac{1}{g \rho} \frac{\partial p}{\partial x}\normalsize = \large\frac{\partial \eta}{\partial x}\normalsize + (\eta z) \beta s_0 \large\frac{\partial S}{\partial x}\normalsize . \qquad (1) </math> 
−  The fluid acceleration caused by this pressure gradient is  +  The fluid acceleration caused by this pressure gradient is balanced by momentum dissipation through turbulent diffusion: 
−  <math>\  +  <math>\large\frac{\partial p}{\rho \partial x}\normalsize = \large\frac{\partial }{\partial z}\normalsize ( K_M \large\frac{\partial u_C}{\partial z}\normalsize) . \qquad (2) </math> 
This equation can be solved considering the boundary conditions:  This equation can be solved considering the boundary conditions:  
−  +  No slip at the bottom: <math>u_C(h)=0</math> , no momentum diffusion through the surface: <math>du_C / dz =0</math> at <math>z=\eta \approx 0</math> and no net discharge: <math>\int_{h}^0 u_C dz =0</math>.  
−  
−  
Because the river runoff is constant and the tide cyclic, the salinity distribution should not change with time in the moving frame. However, the estuarine circulation modifies the salinity distribution, because salt advection <math>u_C(z) S(x,z)</math> has a gradient in <math>x</math>direction. The salt balance is restored by vertical diffusion of salinity:  Because the river runoff is constant and the tide cyclic, the salinity distribution should not change with time in the moving frame. However, the estuarine circulation modifies the salinity distribution, because salt advection <math>u_C(z) S(x,z)</math> has a gradient in <math>x</math>direction. The salt balance is restored by vertical diffusion of salinity:  
−  <math>u_C \  +  <math>u_C \large\frac{\partial \overline S}{\partial x}\normalsize = \large\frac{\partial }{\partial z}\normalsize ( K_S \large\frac{\partial S}{\partial z}\normalsize) . \qquad (3) </math> 
−  This equation can be solved with the boundary conditions:  +  This equation can be solved with the boundary conditions: no salt flux through bottom and surface, <math>\partial S / \partial z =0</math> at <math> z=h, z=0 .</math> 
−  
−  The solution of equations Eqs. (13)  +  The solution of equations Eqs. (13) is given by (shown in figure 2): 
<math> u_C(z)=u_E (19 \zeta^28 \zeta^3) , \; S=\overline S + s_C(z), \; s_C(z)=S_E (5+30 \zeta^245 \zeta^424 \zeta^5), \; \zeta=z/h , \qquad (4) </math>  <math> u_C(z)=u_E (19 \zeta^28 \zeta^3) , \; S=\overline S + s_C(z), \; s_C(z)=S_E (5+30 \zeta^245 \zeta^424 \zeta^5), \; \zeta=z/h , \qquad (4) </math>  
−  <math>u_E=\  +  <math>u_E=\large\frac{g \beta s_0 h}{48 K_M} \frac{\partial \overline S}{\partial x}\normalsize , \; S_E=\large\frac{h^2 u_E}{60 K_S} \frac{\partial \overline S}{\partial x}\normalsize . \qquad (5) </math> 
−  A frequently used estimate of the eddy viscosity is <math>K_M=hU_*^2/U,</math> where <math>U_*=\sqrt{c_D}U</math> is the friction velocity. We then find that the estuarine circulation <math>u_E</math> is related to the Simpson stratification number <math> Si</math> :  +  A frequently used estimate of the eddy viscosity in unstratified water is <math>K_M=hU_*^2/U,</math> where <math>U_*=\sqrt{c_D}U</math> is the friction velocity. We then find that the estuarine circulation <math>u_E</math> is related to the Simpson stratification number <math> Si</math> : 
−  <math>u_E=\  +  <math>u_E=\large\frac{U}{48}\normalsize Si , \quad Si=\large\frac{g}{\rho} \frac{\partial \rho}{\partial x} \frac{h^2}{U_*^2}\normalsize . \qquad (6) </math> 
==Salt balance==  ==Salt balance==  
−  The motion in the moving frame is not exactly cyclic; after one tide the moving frame has moved seaward over a distance <math>T u_R</math> . In order to prevent flushing of the salinity a net landward salt transport should exist. Many seawater intrusion mechanisms can contribute to landward salt transport, as discussed in [[Seawater intrusion and mixing in estuaries]]. In this case, we assume that estuarine circulation is the primary mechanism for seawater intrusion. This means that the loss of salt content from the estuary due to the net shift of the moving frame, <math>Thu_rS</math> , should be compensated by salt transport related to estuarine circulation:  +  The motion in the moving frame is not exactly cyclic; after one tide the moving frame has moved seaward over a distance <math>T u_R</math> . In order to prevent complete flushing of the salinity a net landward salt transport should exist. Many seawater intrusion mechanisms can contribute to landward salt transport, as discussed in the article [[Seawater intrusion and mixing in estuaries]]. In this case, we assume that estuarine circulation is the primary mechanism for seawater intrusion. This means that the loss of salt content from the estuary due to the net shift of the moving frame, <math>Thu_rS</math> , should be compensated by salt transport related to estuarine circulation: 
−  <math>hu_RS=  +  <math>hu_RS=\int_{h}^0 u_C(z)S(z) dz . \qquad (7) </math> 
−  Evaluation of the integral gives, using Eqs. (4,5),  +  Evaluation of the integral gives, using Eqs. (4, 5), 
−  <math>S=L^3 (\  +  <math>S=L^3 (\large\frac{\partial \overline S}{\partial x})^3\normalsize , \; L=0.024 h^2 c (\large\frac{c}{u_R K_S K_M^2})^{1/3}\normalsize , \; c=\sqrt{\beta s_0 gh} . \qquad (8) </math> 
−  <math>L</math> is an estimate for the length of the seawater intrusion due to estuarine circulation; <math>c</math> is the internal wave speed in a fresh water layer of depth <math>h</math> overlying seawater. By approximating <math>\partial \overline S / \partial x \approx  1 / L</math> the maximum estuarine circulation velocity can be written as  +  <math>L</math> is an estimate for the length of the seawater intrusion due to estuarine circulation; <math>c</math> is the internal wave speed in a fresh water layer of depth <math>h</math> overlying seawater. By approximating <math>\partial \overline S / \partial x \approx 1 / L </math> the maximum estuarine circulation velocity can be written as 
−  <math>u_E=0.87(u_R c^2 \  +  <math>u_E=0.87(u_R c^2 \large\frac{K_S}{K_M})^{1/3}\normalsize. \qquad (9) </math> 
In general <math>u_E>>u_R .</math>  In general <math>u_E>>u_R .</math>  
Line 94:  Line 91:  
[[Image:EstuarineCirculation.jpgthumb450pxrightFigure 3: Schematic representation of estuarine circulation in a prismatic channel.]]  [[Image:EstuarineCirculation.jpgthumb450pxrightFigure 3: Schematic representation of estuarine circulation in a prismatic channel.]]  
−  +  According to equation (4), estuarine circulation corresponds to a net landward flow in a near bottom layer of thickness <math>h_2 \approx 0.58 h</math> and a net seaward flow of the same order in the upper layer (see figure 3). By calling <math>s_1,Q_1</math> the average salinity and average discharge in the upper layer and <math>s_2,Q_2</math> the average salinity and discharge in the bottom layer, we can determine the average salinity difference <math>\Delta s = s_2s_1</math> ("stratification") from the condition of net zero salt transport:  
<math>s_1Q_1 + s_2Q_2 =0 , \qquad (10) </math>  <math>s_1Q_1 + s_2Q_2 =0 , \qquad (10) </math>  
−  <math>Q_1 \approx Q_2  hbu_R \approx b  +  <math>Q_1 \approx Q_2  hbu_R \approx b \int_{h_2}^0 u_C dz \approx 0.26hu_E . \qquad (11) </math> 
This gives <math> \Delta s / s \approx 4.4 u_R / u_E . \qquad (12) </math>  This gives <math> \Delta s / s \approx 4.4 u_R / u_E . \qquad (12) </math>  
Line 104:  Line 101:  
The "bulk" Richardson number <math>Ri= \Delta \rho gh / U^2</math> is a measure for the degree of stratification in an estuary. For small Richardson numbers (typically smaller than 0.25) the diffusion coefficients are not strongly influenced by vertical density gradients. In this case one may approximate <math>K_M \approx K_S \approx K = c_D hU.</math> The friction coefficient <math>c_D</math> is of the order of 0.002. Substitution in the expressions for <math>L</math> (Eq. 8), <math>u_E</math> (Eq. 9) and <math>\Delta s</math> (Eq. 12) yields  The "bulk" Richardson number <math>Ri= \Delta \rho gh / U^2</math> is a measure for the degree of stratification in an estuary. For small Richardson numbers (typically smaller than 0.25) the diffusion coefficients are not strongly influenced by vertical density gradients. In this case one may approximate <math>K_M \approx K_S \approx K = c_D hU.</math> The friction coefficient <math>c_D</math> is of the order of 0.002. Substitution in the expressions for <math>L</math> (Eq. 8), <math>u_E</math> (Eq. 9) and <math>\Delta s</math> (Eq. 12) yields  
−  <math> u_E=0.87 (u_R c^2)^{1/3} , \; \  +  <math> u_E=0.87 (u_R c^2)^{1/3} , \; \large\frac{\Delta s}{s}\normalsize = 4.4 (\large\frac{u_R}{c})^{2/3}\normalsize , \; L = \large\frac{0.024}{c_D}\normalsize (\large\frac{c}{u_R})^{1/3} \frac{hc}{U}\normalsize . \qquad (13) </math> 
The estuarine circulation and the stratification increase with increasing river runoff; they do not depend on the tidal velocity.  The estuarine circulation and the stratification increase with increasing river runoff; they do not depend on the tidal velocity.  
−  In case of strong (but continuous) stratification the Richardson number is larger than 1. In this case the diffusion coefficients <math>K_M, K_S</math> depend on <math>Ri</math> and are strongly reduced following  +  In case of strong (but continuous) stratification the Richardson number is larger than 1. In this case the diffusion coefficients <math>K_M, K_S</math> depend on <math>Ri</math> and are strongly reduced following relations experimentally established by MunkAnderson <ref>Munk, W.H. and Anderson, E.R. (1948). Notes on the theory of the thermocline. J. Marine Res. 7: 276295</ref> 
−  <math>K_M = K(1+10Ri)^{1/2} \approx 0.316 K Ri^{1/2} , \; K_S = K (1+3.3Ri)^{3/2} \approx 0.167 K Ri^{3/2} , \; Ri=(\large \frac{c}{U} )^2 \frac{ \Delta s}{s} \normalsize . \qquad (14) </math>  +  <math>K_M = K(1+10Ri)^{1/2} \approx 0.316 K Ri^{1/2} , \; K_S = K (1+3.3Ri)^{3/2} \approx 0.167 K Ri^{3/2} , \; Ri=(\large\frac{c}{U})^2 \frac{\Delta s}{s}\normalsize . \qquad (14) </math> 
Substitution of these approximations yields  Substitution of these approximations yields  
−  <math>u_E=0.3 U , \; \  +  <math>u_E=0.3 U , \; \large\frac{\Delta s}{s}\normalsize = 13 \large\frac{u_R}{U}\normalsize , \; L = \large\frac{0.8}{c_D}\normalsize h (\large\frac{c}{U})^3 (\frac{u_R}{U})^{1/2}\normalsize . \qquad (15) </math> 
With an increase of runoff <math>u_R</math> the estuarine circulation <math>u_E</math> remains unchanged, whereas the intrusion length <math>L</math> increases. This counterintuitive result is due to the strong increase of stratification with increasing runoff.  With an increase of runoff <math>u_R</math> the estuarine circulation <math>u_E</math> remains unchanged, whereas the intrusion length <math>L</math> increases. This counterintuitive result is due to the strong increase of stratification with increasing runoff.  
Line 120:  Line 117:  
==Observations==  ==Observations==  
−  The assumption that seawater intrusion is due only to gravitational circulation strongly limits the salt intrusion length <math>L.</math> The previous formulas indicate for realistic values <math>U=1 m/s, u_R=0.05 m/s, c_D=0.0010.003</math> an intrusion length <math>L \approx  +  The assumption that seawater intrusion is due only to gravitational circulation strongly limits the salt intrusion length <math>L.</math> The previous formulas indicate for realistic values <math>U=1 m/s, h=10 m, u_R=0.05 m/s, c_D=0.0010.003</math> an intrusion length <math>L \approx 0.56 km</math>. Observed salt intrusion lengths are typically several times larger, because other processes also contribute to seawater intrusion (see [[Seawater intrusion and mixing in estuaries]]). However, a much larger salt intrusion length (i.e. much smaller longitudinal density gradient) also implies a much smaller estuarine circulation <math>u_E</math>, according to equation (5). This is in contradiction with field observations, that show strong residual circulation (order 10 cm/s), even in the case of small longitudinal and vertical salinity gradients. Reasons for this discrepancy will be discussed is the following sections. 
==Asymmetric Tidal Mixing==  ==Asymmetric Tidal Mixing==  
Line 126:  Line 123:  
[[Image:EbbFloodProfiles.jpgthumb400pxrightFigure 4: Schematic representation of vertical profiles of salinity (right) and alongchannel velocity (left) for flood and ebb.]]  [[Image:EbbFloodProfiles.jpgthumb400pxrightFigure 4: Schematic representation of vertical profiles of salinity (right) and alongchannel velocity (left) for flood and ebb.]]  
−  In 1986, Van Aken <ref>Van Aken, H. M. (1986). The onset of seasonal stratification in shelf seas due to differential advection in the presence of a salinity gradient. Cont. Shelf Res. 5: 475485</ref> explained the onset of seasonal stratification in the North Sea by the process of tidal advection of surface water with lower density, related to the presence of strong salinity gradients. A few years later, Simpson et al. <ref>Simpson, J.H., Brown, J., Matthews, J.and Allen, G. (1990) Tidal straining, density currents, and stirring in the control of estuarine stratification. Estuaries 13: 125132</ref> observed and modeled tidal variation of stratification in Liverpool Bay by this differential advection mechanism, which they called "tidal straining". Numerical model studies by Jay and Musiak <ref>Jay, D.A. and Musiak, J.D. (1996) Internal tidal asymmetry in channel flows: origins and consequences. In Mixing Processes in Estuaries and Coastal Seas, ed. C Pattiaratchi, Am. Geophys. Union, Washington: 211–49</ref> and Scott <ref>Scott, C.F. (1994) A numerical study of the interaction of tidal oscillations and nonlinearities in an estuary. Estuarine Coast. Shelf Sci. 39: 477–96</ref> revealed that tidal straining can strongly influence estuarine circulation. When during flood high salinity seawater enters the estuary, the less dense estuarine surface water is repulsed. Vertical salinity stratification is therefore reduced, a process further strengthened in some cases by convective overturning <ref>Prandle,  +  In 1986, Van Aken <ref>Van Aken, H. M. (1986). The onset of seasonal stratification in shelf seas due to differential advection in the presence of a salinity gradient. Cont. Shelf Res. 5: 475485</ref> explained the onset of seasonal stratification in the North Sea by the process of tidal advection of surface water with lower density, related to the presence of strong salinity gradients. A few years later, Simpson et al. <ref>Simpson, J.H., Brown, J., Matthews, J.and Allen, G. (1990) Tidal straining, density currents, and stirring in the control of estuarine stratification. Estuaries 13: 125132</ref> observed and modeled tidal variation of stratification in Liverpool Bay by this differential advection mechanism, which they called "tidal straining". Numerical model studies by Jay and Musiak <ref>Jay, D.A. and Musiak, J.D. (1996) Internal tidal asymmetry in channel flows: origins and consequences. In Mixing Processes in Estuaries and Coastal Seas, ed. C Pattiaratchi, Am. Geophys. Union, Washington: 211–49</ref> and Scott <ref>Scott, C.F. (1994) A numerical study of the interaction of tidal oscillations and nonlinearities in an estuary. Estuarine Coast. Shelf Sci. 39: 477–96</ref> revealed that tidal straining can strongly influence estuarine circulation. When during flood high salinity seawater enters the estuary, the less dense estuarine surface water is repulsed. Vertical salinity stratification is therefore reduced, a process further strengthened in some cases by convective overturning (seawater being advected by tidal flow on top of less dense estuarine water <ref>Prandle, 2004 Saline intrusion in partially mixed estuaries. Estuarine, Coastal and Shelf Science 59: 385397</ref>). For sufficiently weak stratification, the surface flow exchanges momentum with the nearbottom flow, which is accelerated in upstream direction. During ebb, the tidal flow restores stratification, by advecting less dense estuarine waters seaward over the denser nearbottom waters. The momentum exchange between surface and nearbottom water is now reduced; seaward tidal flow is concentrated near the surface, see figure 4. The net effect is an enhancement of residual landward nearbottom flow, as shown by numerical experiments <ref>Burchard, H. and Baumert, H. (1998) The formation of estuarine turbidity maxima due to density effects in the salt wedge. A hydrodynamic process study. J. Phys. Oceanogr. 28: 209–321</ref>. 
+  
+  [[Image:RdamWaterwayVelSal.jpgthumb600pxrightFigure 5: Velocity and salinity distributions in a longitudinal section of Rotterdam Waterway on 9 September 1967 at different tidal phases: at the end of flood tide (H.S.W.=high slack water), during maximum ebb flow, at the end of ebb tide (L.W.S.low slack water) and during maximum flood flow. The local salinity maximum halfway the estuary at maximum ebb is related to outflow of the Old Meuse tributary. Characteristics of Rotterdam Waterway are: approximately prismatic channel, width 300 m, depth 15 m, tidal range 2 m, maximum depthaveraged tidal velocity 1 m/s. ]]  
Enhancement of estuarine circulation by asymmetric tidal mixing only occurs in situations of weak stratification. In case of strong stratification (weak tides, high river runoff) momentum exchange between surface and nearbottom water hardly occurs (see [[Salt wedge estuaries]]). The major effect of weak tides in the latter case is an increase of eddy viscosity in the bottom layer, which reduces gravitational circulation <ref>Cheng P., De Swart, H.E. and ValleLevinson, A. (2013) Role of asymmetric tidal mixing in the subtidal dynamics of narrow estuaries. J. Geophys. Res. 118: 2623–2639</ref>. Enhancement of estuarine circulation by asymmetric tidal mixing is strongest for small Simpson numbers (typically <math>Ri < 1</math>) and periodic or weak stratification.  Enhancement of estuarine circulation by asymmetric tidal mixing only occurs in situations of weak stratification. In case of strong stratification (weak tides, high river runoff) momentum exchange between surface and nearbottom water hardly occurs (see [[Salt wedge estuaries]]). The major effect of weak tides in the latter case is an increase of eddy viscosity in the bottom layer, which reduces gravitational circulation <ref>Cheng P., De Swart, H.E. and ValleLevinson, A. (2013) Role of asymmetric tidal mixing in the subtidal dynamics of narrow estuaries. J. Geophys. Res. 118: 2623–2639</ref>. Enhancement of estuarine circulation by asymmetric tidal mixing is strongest for small Simpson numbers (typically <math>Ri < 1</math>) and periodic or weak stratification.  
−  
−  
Figure 5 illustrates asymmetric tidal mixing in Rotterdam Waterway. The velocity and salinity profiles were measured on 9 September 1967 at an estimated net discharge of about 1000 m3/s (runoff velocity <math>u_R \approx 0.2 m/s</math>). Characteristic values of the Richardson and Simpson numbers are <math>Ri \approx 1.8, Si \approx 0.5 .</math> The sections at maximum ebb and maximum flood flow display asymmetry between nearbottom flow velocities; nearbottom velocities are substantially higher during flood than during ebb (as schematically depicted in figure 4). The sections at H.S.W. and L.S.W. show asymmetry in stratification: at the end of the flood the water column is less stratified than at the end of the ebb. Another survey in Rotterdam Waterway was conducted by De Nijs et al. in 2006 at higher river discharge (stronger stratification) <ref>De Nijs, M.A.J., Pietrzak, J.D. and Winterwerp, J.C.(2011) Advection of the Salt Wedge and Evolution of the Internal Flow Structure in the Rotterdam Waterway. J. Phys. Oceanogr. 41: 327</ref>. These authors show that in Rotterdam Waterway asymmetry in tidal mixing is primarily related to asymmetry of the tide at the mouth (faster tidal rise than tidal fall, which enhances flood currents over ebb currents), more than to differential advection. This is consistent with the fast thickening of the bottom boundary layer during flood, caused by strong bedgenerated turbulence. A similar behaviour was observed in the Hudson estuary <ref>Nepf, H.M. and Geyer, W.R. (1996) Intratidal variations in stratification and mixing in the Hudson estuary. J. Geophys. Res. 101: 1207912086</ref>.  Figure 5 illustrates asymmetric tidal mixing in Rotterdam Waterway. The velocity and salinity profiles were measured on 9 September 1967 at an estimated net discharge of about 1000 m3/s (runoff velocity <math>u_R \approx 0.2 m/s</math>). Characteristic values of the Richardson and Simpson numbers are <math>Ri \approx 1.8, Si \approx 0.5 .</math> The sections at maximum ebb and maximum flood flow display asymmetry between nearbottom flow velocities; nearbottom velocities are substantially higher during flood than during ebb (as schematically depicted in figure 4). The sections at H.S.W. and L.S.W. show asymmetry in stratification: at the end of the flood the water column is less stratified than at the end of the ebb. Another survey in Rotterdam Waterway was conducted by De Nijs et al. in 2006 at higher river discharge (stronger stratification) <ref>De Nijs, M.A.J., Pietrzak, J.D. and Winterwerp, J.C.(2011) Advection of the Salt Wedge and Evolution of the Internal Flow Structure in the Rotterdam Waterway. J. Phys. Oceanogr. 41: 327</ref>. These authors show that in Rotterdam Waterway asymmetry in tidal mixing is primarily related to asymmetry of the tide at the mouth (faster tidal rise than tidal fall, which enhances flood currents over ebb currents), more than to differential advection. This is consistent with the fast thickening of the bottom boundary layer during flood, caused by strong bedgenerated turbulence. A similar behaviour was observed in the Hudson estuary <ref>Nepf, H.M. and Geyer, W.R. (1996) Intratidal variations in stratification and mixing in the Hudson estuary. J. Geophys. Res. 101: 1207912086</ref>.  
Line 136:  Line 133:  
==Lateral estuarine circulation==  ==Lateral estuarine circulation==  
+  
+  In natural estuaries the depth decreases towards the lateral boundaries (in a meander the depth decreases towards the inner channel bend). Tidal velocities are lower near the channel banks than at the channel centre. The salinity front is therefore not advected at the same speed at different lateral locations. This results in lateral salinity gradients. During flood (especially shortly before the end of the flood phase) the salinity is higher at the channel centre than near the channel banks. During ebb the situation is reversed; near the surface, the salinity at the channel centre is lower than near the channel banks. The lateral salinity gradients related to this differential tidal salinity advection, induce a baroclinic (densitydriven) lateral current. This current has a vertical structure; during flood it is directed near the bottom toward the channel banks and toward the channel centre near the surface, see figure 6. The resulting flow convergence at the channel centre is often visible in the field by a scum line following the channel axis <ref>Nunes, R.A. and Simpson, J.H. (1985) Axial Convergence in a Wellmixed Estuary. Coastal and Shelf Science 20: 637649</ref>. During ebb the lateral exchange flow is divergent at the surface and convergent near the bottom. This lateral baroclinic exchange flow has important consequences for longitudinal dispersion in natural estuaries, as first pointed out by R. Smith <ref>Smith, R. (1976) Longitudinal dispersion of a buoyant contaminant in a shallow channel. J. Fluid Mech. 78: 677–88</ref>, as explained in the article [[Seawater intrusion and mixing in estuaries]].  
−  [[Image:LateralEstCirc.jpgthumb  +  [[Image:LateralEstCirc.jpgthumb900pxcenterFigure 6: Lateral densitydriven exchange flow related to crosschannel depth variation. Panel a: Crosssectional salinity distribution and flow pattern shortly before the end of the flood phase. Higher salinity is indicated by darker blue. Panel b: The same in the last phase of ebb. The strength of alongchannel flood flow in different parts of the crosssection is indicated by crossed circles; the strength of ebb flow by dotted circles. The lateral exchange flow influences the strength and the distribution of longitudinal estuarine circulation by altering the vertical distribution of alongchannel flood and ebb flow. Panel c: Residual flow pattern (estuarine circulation). The crossed circle indicates upestuary residual flow and the dotted circles downestuary residual flow.]] 
−  
Lateral exchange flows also influence estuarine circulation. Downwelling flow at the channel axis during flood transports momentum down the vertical; surface flow is decelerated and near bottom flow accelerated in upestuary direction. During ebb, vertical momentum transport reduces turbulent momentum dissipation by reducing the vertical shear of the tidal current. This also favours the baroclinic flow component <ref>Burchard, H. and Schuttelaars, H.H. (2012) Analysis of Tidal Straining as Driver for Estuarine Circulation in WellMixed Estuaries. J. Phys. Oceanograph. 42: 261271</ref>. Lateral baroclinic exchange flow therefore enhances estuarine circulation. It also induces an a net downstream flow component in the shallower parts of the crosssection (lateral inflow of downestuary momentum from the channel centre during ebb, outflow of upestuary momentum toward the channel centre during flood, see figure 6); this flow component enhances the downestuary flow component related to river runoff <ref>Burchard, H., Hetland, R.D., Schulz, E. and Schuttelaars, H.H. (2011)Drivers of Residual Estuarine Circulation in Tidally Energetic Estuaries: Straight and Irrotational Channels with Parabolic Cross Section. J.Phys.Oceanogr. 41: 548570</ref>.  Lateral exchange flows also influence estuarine circulation. Downwelling flow at the channel axis during flood transports momentum down the vertical; surface flow is decelerated and near bottom flow accelerated in upestuary direction. During ebb, vertical momentum transport reduces turbulent momentum dissipation by reducing the vertical shear of the tidal current. This also favours the baroclinic flow component <ref>Burchard, H. and Schuttelaars, H.H. (2012) Analysis of Tidal Straining as Driver for Estuarine Circulation in WellMixed Estuaries. J. Phys. Oceanograph. 42: 261271</ref>. Lateral baroclinic exchange flow therefore enhances estuarine circulation. It also induces an a net downstream flow component in the shallower parts of the crosssection (lateral inflow of downestuary momentum from the channel centre during ebb, outflow of upestuary momentum toward the channel centre during flood, see figure 6); this flow component enhances the downestuary flow component related to river runoff <ref>Burchard, H., Hetland, R.D., Schulz, E. and Schuttelaars, H.H. (2011)Drivers of Residual Estuarine Circulation in Tidally Energetic Estuaries: Straight and Irrotational Channels with Parabolic Cross Section. J.Phys.Oceanogr. 41: 548570</ref>.  
Line 146:  Line 144:  
==Other lateral circulation patterns==  ==Other lateral circulation patterns==  
−  [[Image:CoriolisEstCirc.jpgthumb300pxrightFigure 7: Influence of earth rotation on estuarine circulation. The upestuary residual flow (crossed circles) is shifted toward the right channel bank (looking in upestuary direction, northern hemisphere) and the downestuary residual flow (dotted circles) toward the left bank. ]]  +  [[Image:CoriolisEstCirc.jpgthumb300pxrightFigure 7: Influence of earth's rotation on estuarine circulation. The upestuary residual flow (crossed circles) is shifted toward the right channel bank (looking in upestuary direction, northern hemisphere) and the downestuary residual flow (dotted circles) toward the left bank. ]] 
−  Several mechanisms can alter the symmetric flow pattern of figure 6. In wide estuaries the symmetry is broken by earth rotation. Currents are deviated to the right on the northern hemisphere and to the left on the southern hemisphere (see [[Coriolis acceleration]]). Looking in upestuary direction, flood currents are concentrated at the right side of a straight tidal channel and ebb currents at the left side (northern hemisphere). Estuarine circulation is influenced by earth rotation in a similar way <ref>Huijts, K.M.H., Schuttelaars, B.H.M., De Swart, H.E. and Friedrichs, C.T. (2009) Analytical study of the transverse distribution of alongchannel and transverse residual flows in tidal estuaries. Continental Shelf Res. 29: 89100</ref>: the residual upestuary flow is shifted to the right bank and the downestuary flow to the left bank, see figure 7. The mean crosssectional salinity distribution is tilted accordingly.  +  Several mechanisms can alter the symmetric flow pattern of figure 6. In wide estuaries the symmetry is broken by earth's rotation. Currents are deviated to the right on the northern hemisphere and to the left on the southern hemisphere (see [[Coriolis acceleration]]). Looking in upestuary direction, flood currents are concentrated at the right side of a straight tidal channel and ebb currents at the left side (northern hemisphere). Estuarine circulation is influenced by earth's rotation in a similar way <ref>Huijts, K.M.H., Schuttelaars, B.H.M., De Swart, H.E. and Friedrichs, C.T. (2009) Analytical study of the transverse distribution of alongchannel and transverse residual flows in tidal estuaries. Continental Shelf Res. 29: 89100</ref>: the residual upestuary flow is shifted to the right bank and the downestuary flow to the left bank, see figure 7. The mean crosssectional salinity distribution is tilted accordingly. 
−  Channel meandering induces transverse circulation, directed to the inner channel bend in the lower part of the vertical and to the outer bend near the surface. The spiraling tidal flow modifies the salinity distribution and the densitydriven longitudinal residual flow in a similar way as earth rotation when the tidal flow follows a channel meander with the inner channel bend at the right bank <ref>Basdurak, N.B. and ValleLevinson, A. (2013) Tidal variability of lateral advection in a coastal plain estuary. Continental Shelf Res. 61–62: 85–97</ref>. For an inner channel bend at the left bank the picture is mirrored; in this case the effects of channel meandering and earth rotation are opposite. Channel bends and other topographic structures (bars, shoals) generate longitudinal circulating flow patterns in estuaries also in the absence of density differences. These horizontal topographically induced flow circulations are often stronger than horizontal estuarine circulation generated by density gradients. However, residual flow circulation in the longitudinalvertical plane is mainly due to the density difference between seawater and river water.  +  Channel meandering induces transverse circulation, directed to the inner channel bend in the lower part of the vertical and to the outer bend near the surface. The spiraling tidal flow modifies the salinity distribution and the densitydriven longitudinal residual flow in a similar way as earth's rotation when the tidal flow follows a channel meander with the inner channel bend at the right bank <ref>Basdurak, N.B. and ValleLevinson, A. (2013) Tidal variability of lateral advection in a coastal plain estuary. Continental Shelf Res. 61–62: 85–97</ref>. For an inner channel bend at the left bank the picture is mirrored; in this case the effects of channel meandering and earth's rotation are opposite. Channel bends and other topographic structures (bars, shoals) generate longitudinal circulating flow patterns in estuaries also in the absence of density differences. These horizontal topographically induced flow circulations are often stronger than horizontal estuarine circulation generated by density gradients. However, residual flow circulation in the longitudinalvertical plane is mainly due to the density difference between seawater and river water. 
==Consequence of estuarine circulation==  ==Consequence of estuarine circulation==  
−  Estuarine circulation is highly relevant for transport of fine sediments and sedimentation (see [[Dynamics of mud transport]]). The upestuary nearbottom flow captures a substantial part of the downstream flux of fluvial sediment, and moves this sediment back toward the tip of the salt intrusion zone. The upestuary nearbottom flow contributes also to import of marine sediment into the estuary. This sediment settles on the bed during periods of neap tide and slack water. It is resuspended in periods of high tidal currents, leading to a strong increase in turbidity. Decomposition of organic material in the high turbidity zone can cause oxygen depletion, especially in the lower part of the vertical where mixing with oxygenrich surface water is inhibited by stratification. However, estuarine circulation is not the only cause of upestuary sediment transport. Many field and model investigations show that tidal asymmetry due to nonlinear tidal propagation in shallow water (leading to stronger flood currents than ebb currents and longer highwater slack periods than lowwater slack periods) is often a more important mechanism.  +  Estuarine circulation is highly relevant for transport of fine sediments and sedimentation (see [[Dynamics of mud transport]]). The upestuary nearbottom flow captures a substantial part of the downstream flux of fluvial sediment, and moves this sediment back toward the tip of the salt intrusion zone. The upestuary nearbottom flow contributes also to import of marine sediment into the estuary. This sediment settles on the bed during periods of neap tide and slack water. It is resuspended in periods of high tidal currents, leading to a strong increase in turbidity. Decomposition of organic material in the high turbidity zone can cause oxygen depletion, especially in the lower part of the vertical where mixing with oxygenrich surface water is inhibited by stratification. However, estuarine circulation is not the only cause of upestuary sediment transport. Many field and model investigations show that tidal asymmetry due to nonlinear tidal propagation in shallow water (leading to stronger flood currents than ebb currents and longer highwater slack periods than lowwater slack periods) is often a more important mechanism, see the article [[Tidal asymmetry and tidal basin morphodynamics]]. 
==Related articles==  ==Related articles==  
−  +  :[[Salt wedge estuaries]]  
−  [[Salt wedge estuaries]]  +  :[[Seawater intrusion and mixing in estuaries]] 
−  +  :[[Coriolis acceleration]]  
−  [[Seawater intrusion and mixing in estuaries]]  +  :[[Dynamics of mud transport]] 
−  
−  [[Coriolis acceleration]]  
−  
−  [[Dynamics of mud transport]]  
==References==  ==References==  
Line 176:  Line 170:  
AuthorName=Dronkers J}}  AuthorName=Dronkers J}}  
−  
[[Category:Estuaries and tidal rivers]]  [[Category:Estuaries and tidal rivers]]  
[[Category:Hydrodynamics]]  [[Category:Hydrodynamics]]  
+  [[Category:Land and ocean interactions]] 
Latest revision as of 21:06, 16 October 2019
Definition of Estuarine circulation:
Residual flow pattern in an estuary induced by the density difference between seawater and river water.
This is the common definition for Estuarine circulation, other definitions can be discussed in the article

Seawater has a higher density than river water, due to its high salinity. The density difference is of the order of 3%. In estuaries, where seawater and river water meet, seawater will tend to intrude beneath the river flow. However, mixing will also occur between seawater and river water, whereby seawater is returned to the sea by the river flow via the upper water layer. This allows 'new' seawater to enter the estuary via the lower water layers. In this way, a circulation of seawater is created, which is directed inland along the bottom and seaward along the surface. This circulation is called estuarine circulation or gravitational circulation. It is driven by the longitudinal density gradient along the estuary. In addition to the density differences along the estuary, the circulation depends on the mixing of seawater with river water. Mixing is the result of turbulence, which in estuaries is mainly caused by strong tidal flows. Most of the mixing takes place vertically in the water column, due to turbulence generated by tidal flow along the estuarine bottom. The different processes that drive estuarine circulation are discussed in this article. In estuaries with little tide there is hardly any mixing of seawater with river water. In that case a more or less stationary salt wedge forms along the bottom; this case is discussed in the article Salt wedge estuaries.
Contents
Introduction
The classical theory of estuarine circulation was first developed by Hansen and Rattray in 1965 ^{[1]}, after earlier investigations by D. Pritchard in the St. James river ^{[2]}. Many other authors have subsequently contributed to further verification and refinement^{[3]}^{[4]}.
In the classical theory of estuarine circulation, advection of seawater up and down the estuary by the tidal motion and associated mixing processes are not considered explicitly. These models consider a stationary salinity field in a frame that moves cyclically up and down the estuary and a tideaveraged coefficient for the mixing processes. The effect of tidal motion is further restricted to turbulence generation at the bed, which enhances momentum dissipation and mixing in and across the bottom boundary layer.
Notation
In the following, some basic notions of the classical theory are presented. We consider a straight prismatic estuary (uniform depth [math]h[/math] and width [math]b[/math] ); the width is small enough for neglecting the influence of earth's rotation (Rossby radius [math]U/fb \lt \lt 1[/math] ). We use the following notations (see also figure 1):
 [math]x[/math] is the coordinate along the estuarine axis ([math]x=0[/math] at the sea boundary, [math]x=L[/math] at the landward boundary of salt intrusion);
 [math]u[/math] is the velocity in the [math]x[/math]direction, [math]u(z) = u_R + u_C(z) + U \cos \omega t,[/math] where [math]u_R[/math] is the velocity related to river runoff (river discharge divided by estuarine crosssection), [math]u_C(z)[/math] is the estuarine circulation, [math]U[/math] is the depthaveraged tidal velocity amplitude and [math]\omega = 2 \pi / T[/math] is the tidal radial frequency;
 [math]S=s/s_0[/math] is the ratio of the local salinity [math]s[/math] and the seawater salinity [math]s_0 :[/math] [math]S[/math] is related to the density [math]\rho[/math] by the relation [math]\rho=\rho_0(1+\beta s_0 S),[/math] where [math]\rho_0[/math] is the fresh water density and the expansion coefficient [math]\beta=7.7 \times 10^{4};[/math]
 [math] p(z) = g \int_z^\eta \rho dz[/math] is the hydrostatic pressure; [math]\eta[/math] is the surface elevation relative to a horizontal plane and [math]z[/math] is the depth relative to the surface;
 [math]K_M[/math] is the eddy viscosity (vertical momentum diffusion coefficient [math][m^2/s] [/math] ) and [math]K_S[/math] the eddy diffusivity (vertical salinity diffusion coefficient).
Depthaveraging is denoted by an overbar.
Assumptions
The following additional assumptions are made:
 The diffusion coefficients [math]K_M, K_S[/math] do not depend on the spatial and temporal coordinates [math]x,z,t[/math] ; this means that the classical model is realistic only for continuously stratified estuaries (no layered salinity structure; such estuaries are often called "partially mixed", see figure 1);
 The river flow velocity [math]u_R[/math] and the tidal velocity amplitude [math]U[/math] do not depend on the spatial coordinates [math]x,z[/math] ;
 The river flow velocity is much smaller than the tidal velocity, [math]u_R \lt \lt U[/math] ;
 The water surface variation is much smaller than depth, [math]\eta \lt \lt h[/math] ; Stokes drift is neglected;
 The vertical flow velocity [math] w[/math] is neglected;
 Side wall effects are neglected, implying uniformity of all variables in lateral [math]y[/math]direction;
 Longitudinal gradients in the velocity [math]u[/math] are neglected;
 The longitudinal salinity gradient is the same at all depths ([math]\large\frac{\partial S}{\partial x}\normalsize=\large\frac{\partial \overline S}{\partial x}\normalsize ).[/math]
Model formulation
We consider the influence of a constant river flow [math]u_R[/math] on the salinity distribution in a frame moving with the velocity [math]u_R + U \cos \omega t[/math]. In the moving frame the velocity [math]u[/math] is timeindependent and depends only on the vertical coordinate [math]z[/math]; we thus have [math]u=u_C(z)[/math]. The only driving force acting on the fluid in the moving frame is the pressure gradient [math]\partial p / \partial x [/math], which is a linearly increasing function of depth ([math]z[/math]) below mean water level:
[math]\large\frac{1}{g \rho} \frac{\partial p}{\partial x}\normalsize = \large\frac{\partial \eta}{\partial x}\normalsize + (\eta z) \beta s_0 \large\frac{\partial S}{\partial x}\normalsize . \qquad (1) [/math]
The fluid acceleration caused by this pressure gradient is balanced by momentum dissipation through turbulent diffusion:
[math]\large\frac{\partial p}{\rho \partial x}\normalsize = \large\frac{\partial }{\partial z}\normalsize ( K_M \large\frac{\partial u_C}{\partial z}\normalsize) . \qquad (2) [/math]
This equation can be solved considering the boundary conditions:
No slip at the bottom: [math]u_C(h)=0[/math] , no momentum diffusion through the surface: [math]du_C / dz =0[/math] at [math]z=\eta \approx 0[/math] and no net discharge: [math]\int_{h}^0 u_C dz =0[/math].
Because the river runoff is constant and the tide cyclic, the salinity distribution should not change with time in the moving frame. However, the estuarine circulation modifies the salinity distribution, because salt advection [math]u_C(z) S(x,z)[/math] has a gradient in [math]x[/math]direction. The salt balance is restored by vertical diffusion of salinity:
[math]u_C \large\frac{\partial \overline S}{\partial x}\normalsize = \large\frac{\partial }{\partial z}\normalsize ( K_S \large\frac{\partial S}{\partial z}\normalsize) . \qquad (3) [/math]
This equation can be solved with the boundary conditions: no salt flux through bottom and surface, [math]\partial S / \partial z =0[/math] at [math] z=h, z=0 .[/math]
The solution of equations Eqs. (13) is given by (shown in figure 2):
[math] u_C(z)=u_E (19 \zeta^28 \zeta^3) , \; S=\overline S + s_C(z), \; s_C(z)=S_E (5+30 \zeta^245 \zeta^424 \zeta^5), \; \zeta=z/h , \qquad (4) [/math]
[math]u_E=\large\frac{g \beta s_0 h}{48 K_M} \frac{\partial \overline S}{\partial x}\normalsize , \; S_E=\large\frac{h^2 u_E}{60 K_S} \frac{\partial \overline S}{\partial x}\normalsize . \qquad (5) [/math]
A frequently used estimate of the eddy viscosity in unstratified water is [math]K_M=hU_*^2/U,[/math] where [math]U_*=\sqrt{c_D}U[/math] is the friction velocity. We then find that the estuarine circulation [math]u_E[/math] is related to the Simpson stratification number [math] Si[/math] :
[math]u_E=\large\frac{U}{48}\normalsize Si , \quad Si=\large\frac{g}{\rho} \frac{\partial \rho}{\partial x} \frac{h^2}{U_*^2}\normalsize . \qquad (6) [/math]
Salt balance
The motion in the moving frame is not exactly cyclic; after one tide the moving frame has moved seaward over a distance [math]T u_R[/math] . In order to prevent complete flushing of the salinity a net landward salt transport should exist. Many seawater intrusion mechanisms can contribute to landward salt transport, as discussed in the article Seawater intrusion and mixing in estuaries. In this case, we assume that estuarine circulation is the primary mechanism for seawater intrusion. This means that the loss of salt content from the estuary due to the net shift of the moving frame, [math]Thu_rS[/math] , should be compensated by salt transport related to estuarine circulation:
[math]hu_RS=\int_{h}^0 u_C(z)S(z) dz . \qquad (7) [/math]
Evaluation of the integral gives, using Eqs. (4, 5),
[math]S=L^3 (\large\frac{\partial \overline S}{\partial x})^3\normalsize , \; L=0.024 h^2 c (\large\frac{c}{u_R K_S K_M^2})^{1/3}\normalsize , \; c=\sqrt{\beta s_0 gh} . \qquad (8) [/math]
[math]L[/math] is an estimate for the length of the seawater intrusion due to estuarine circulation; [math]c[/math] is the internal wave speed in a fresh water layer of depth [math]h[/math] overlying seawater. By approximating [math]\partial \overline S / \partial x \approx 1 / L [/math] the maximum estuarine circulation velocity can be written as
[math]u_E=0.87(u_R c^2 \large\frac{K_S}{K_M})^{1/3}\normalsize. \qquad (9) [/math]
In general [math]u_E\gt \gt u_R .[/math]
Stratification
According to equation (4), estuarine circulation corresponds to a net landward flow in a near bottom layer of thickness [math]h_2 \approx 0.58 h[/math] and a net seaward flow of the same order in the upper layer (see figure 3). By calling [math]s_1,Q_1[/math] the average salinity and average discharge in the upper layer and [math]s_2,Q_2[/math] the average salinity and discharge in the bottom layer, we can determine the average salinity difference [math]\Delta s = s_2s_1[/math] ("stratification") from the condition of net zero salt transport:
[math]s_1Q_1 + s_2Q_2 =0 , \qquad (10) [/math]
[math]Q_1 \approx Q_2  hbu_R \approx b \int_{h_2}^0 u_C dz \approx 0.26hu_E . \qquad (11) [/math]
This gives [math] \Delta s / s \approx 4.4 u_R / u_E . \qquad (12) [/math]
The "bulk" Richardson number [math]Ri= \Delta \rho gh / U^2[/math] is a measure for the degree of stratification in an estuary. For small Richardson numbers (typically smaller than 0.25) the diffusion coefficients are not strongly influenced by vertical density gradients. In this case one may approximate [math]K_M \approx K_S \approx K = c_D hU.[/math] The friction coefficient [math]c_D[/math] is of the order of 0.002. Substitution in the expressions for [math]L[/math] (Eq. 8), [math]u_E[/math] (Eq. 9) and [math]\Delta s[/math] (Eq. 12) yields
[math] u_E=0.87 (u_R c^2)^{1/3} , \; \large\frac{\Delta s}{s}\normalsize = 4.4 (\large\frac{u_R}{c})^{2/3}\normalsize , \; L = \large\frac{0.024}{c_D}\normalsize (\large\frac{c}{u_R})^{1/3} \frac{hc}{U}\normalsize . \qquad (13) [/math]
The estuarine circulation and the stratification increase with increasing river runoff; they do not depend on the tidal velocity.
In case of strong (but continuous) stratification the Richardson number is larger than 1. In this case the diffusion coefficients [math]K_M, K_S[/math] depend on [math]Ri[/math] and are strongly reduced following relations experimentally established by MunkAnderson ^{[5]}
[math]K_M = K(1+10Ri)^{1/2} \approx 0.316 K Ri^{1/2} , \; K_S = K (1+3.3Ri)^{3/2} \approx 0.167 K Ri^{3/2} , \; Ri=(\large\frac{c}{U})^2 \frac{\Delta s}{s}\normalsize . \qquad (14) [/math]
Substitution of these approximations yields
[math]u_E=0.3 U , \; \large\frac{\Delta s}{s}\normalsize = 13 \large\frac{u_R}{U}\normalsize , \; L = \large\frac{0.8}{c_D}\normalsize h (\large\frac{c}{U})^3 (\frac{u_R}{U})^{1/2}\normalsize . \qquad (15) [/math]
With an increase of runoff [math]u_R[/math] the estuarine circulation [math]u_E[/math] remains unchanged, whereas the intrusion length [math]L[/math] increases. This counterintuitive result is due to the strong increase of stratification with increasing runoff.
Observations
The assumption that seawater intrusion is due only to gravitational circulation strongly limits the salt intrusion length [math]L.[/math] The previous formulas indicate for realistic values [math]U=1 m/s, h=10 m, u_R=0.05 m/s, c_D=0.0010.003[/math] an intrusion length [math]L \approx 0.56 km[/math]. Observed salt intrusion lengths are typically several times larger, because other processes also contribute to seawater intrusion (see Seawater intrusion and mixing in estuaries). However, a much larger salt intrusion length (i.e. much smaller longitudinal density gradient) also implies a much smaller estuarine circulation [math]u_E[/math], according to equation (5). This is in contradiction with field observations, that show strong residual circulation (order 10 cm/s), even in the case of small longitudinal and vertical salinity gradients. Reasons for this discrepancy will be discussed is the following sections.
Asymmetric Tidal Mixing
In 1986, Van Aken ^{[6]} explained the onset of seasonal stratification in the North Sea by the process of tidal advection of surface water with lower density, related to the presence of strong salinity gradients. A few years later, Simpson et al. ^{[7]} observed and modeled tidal variation of stratification in Liverpool Bay by this differential advection mechanism, which they called "tidal straining". Numerical model studies by Jay and Musiak ^{[8]} and Scott ^{[9]} revealed that tidal straining can strongly influence estuarine circulation. When during flood high salinity seawater enters the estuary, the less dense estuarine surface water is repulsed. Vertical salinity stratification is therefore reduced, a process further strengthened in some cases by convective overturning (seawater being advected by tidal flow on top of less dense estuarine water ^{[10]}). For sufficiently weak stratification, the surface flow exchanges momentum with the nearbottom flow, which is accelerated in upstream direction. During ebb, the tidal flow restores stratification, by advecting less dense estuarine waters seaward over the denser nearbottom waters. The momentum exchange between surface and nearbottom water is now reduced; seaward tidal flow is concentrated near the surface, see figure 4. The net effect is an enhancement of residual landward nearbottom flow, as shown by numerical experiments ^{[11]}.
Enhancement of estuarine circulation by asymmetric tidal mixing only occurs in situations of weak stratification. In case of strong stratification (weak tides, high river runoff) momentum exchange between surface and nearbottom water hardly occurs (see Salt wedge estuaries). The major effect of weak tides in the latter case is an increase of eddy viscosity in the bottom layer, which reduces gravitational circulation ^{[12]}. Enhancement of estuarine circulation by asymmetric tidal mixing is strongest for small Simpson numbers (typically [math]Ri \lt 1[/math]) and periodic or weak stratification.
Figure 5 illustrates asymmetric tidal mixing in Rotterdam Waterway. The velocity and salinity profiles were measured on 9 September 1967 at an estimated net discharge of about 1000 m3/s (runoff velocity [math]u_R \approx 0.2 m/s[/math]). Characteristic values of the Richardson and Simpson numbers are [math]Ri \approx 1.8, Si \approx 0.5 .[/math] The sections at maximum ebb and maximum flood flow display asymmetry between nearbottom flow velocities; nearbottom velocities are substantially higher during flood than during ebb (as schematically depicted in figure 4). The sections at H.S.W. and L.S.W. show asymmetry in stratification: at the end of the flood the water column is less stratified than at the end of the ebb. Another survey in Rotterdam Waterway was conducted by De Nijs et al. in 2006 at higher river discharge (stronger stratification) ^{[13]}. These authors show that in Rotterdam Waterway asymmetry in tidal mixing is primarily related to asymmetry of the tide at the mouth (faster tidal rise than tidal fall, which enhances flood currents over ebb currents), more than to differential advection. This is consistent with the fast thickening of the bottom boundary layer during flood, caused by strong bedgenerated turbulence. A similar behaviour was observed in the Hudson estuary ^{[14]}.
Lateral estuarine circulation
In natural estuaries the depth decreases towards the lateral boundaries (in a meander the depth decreases towards the inner channel bend). Tidal velocities are lower near the channel banks than at the channel centre. The salinity front is therefore not advected at the same speed at different lateral locations. This results in lateral salinity gradients. During flood (especially shortly before the end of the flood phase) the salinity is higher at the channel centre than near the channel banks. During ebb the situation is reversed; near the surface, the salinity at the channel centre is lower than near the channel banks. The lateral salinity gradients related to this differential tidal salinity advection, induce a baroclinic (densitydriven) lateral current. This current has a vertical structure; during flood it is directed near the bottom toward the channel banks and toward the channel centre near the surface, see figure 6. The resulting flow convergence at the channel centre is often visible in the field by a scum line following the channel axis ^{[15]}. During ebb the lateral exchange flow is divergent at the surface and convergent near the bottom. This lateral baroclinic exchange flow has important consequences for longitudinal dispersion in natural estuaries, as first pointed out by R. Smith ^{[16]}, as explained in the article Seawater intrusion and mixing in estuaries.
Lateral exchange flows also influence estuarine circulation. Downwelling flow at the channel axis during flood transports momentum down the vertical; surface flow is decelerated and near bottom flow accelerated in upestuary direction. During ebb, vertical momentum transport reduces turbulent momentum dissipation by reducing the vertical shear of the tidal current. This also favours the baroclinic flow component ^{[17]}. Lateral baroclinic exchange flow therefore enhances estuarine circulation. It also induces an a net downstream flow component in the shallower parts of the crosssection (lateral inflow of downestuary momentum from the channel centre during ebb, outflow of upestuary momentum toward the channel centre during flood, see figure 6); this flow component enhances the downestuary flow component related to river runoff ^{[18]}.
Other lateral circulation patterns
Several mechanisms can alter the symmetric flow pattern of figure 6. In wide estuaries the symmetry is broken by earth's rotation. Currents are deviated to the right on the northern hemisphere and to the left on the southern hemisphere (see Coriolis acceleration). Looking in upestuary direction, flood currents are concentrated at the right side of a straight tidal channel and ebb currents at the left side (northern hemisphere). Estuarine circulation is influenced by earth's rotation in a similar way ^{[19]}: the residual upestuary flow is shifted to the right bank and the downestuary flow to the left bank, see figure 7. The mean crosssectional salinity distribution is tilted accordingly.
Channel meandering induces transverse circulation, directed to the inner channel bend in the lower part of the vertical and to the outer bend near the surface. The spiraling tidal flow modifies the salinity distribution and the densitydriven longitudinal residual flow in a similar way as earth's rotation when the tidal flow follows a channel meander with the inner channel bend at the right bank ^{[20]}. For an inner channel bend at the left bank the picture is mirrored; in this case the effects of channel meandering and earth's rotation are opposite. Channel bends and other topographic structures (bars, shoals) generate longitudinal circulating flow patterns in estuaries also in the absence of density differences. These horizontal topographically induced flow circulations are often stronger than horizontal estuarine circulation generated by density gradients. However, residual flow circulation in the longitudinalvertical plane is mainly due to the density difference between seawater and river water.
Consequence of estuarine circulation
Estuarine circulation is highly relevant for transport of fine sediments and sedimentation (see Dynamics of mud transport). The upestuary nearbottom flow captures a substantial part of the downstream flux of fluvial sediment, and moves this sediment back toward the tip of the salt intrusion zone. The upestuary nearbottom flow contributes also to import of marine sediment into the estuary. This sediment settles on the bed during periods of neap tide and slack water. It is resuspended in periods of high tidal currents, leading to a strong increase in turbidity. Decomposition of organic material in the high turbidity zone can cause oxygen depletion, especially in the lower part of the vertical where mixing with oxygenrich surface water is inhibited by stratification. However, estuarine circulation is not the only cause of upestuary sediment transport. Many field and model investigations show that tidal asymmetry due to nonlinear tidal propagation in shallow water (leading to stronger flood currents than ebb currents and longer highwater slack periods than lowwater slack periods) is often a more important mechanism, see the article Tidal asymmetry and tidal basin morphodynamics.
Related articles
 Salt wedge estuaries
 Seawater intrusion and mixing in estuaries
 Coriolis acceleration
 Dynamics of mud transport
References
 ↑ Hansen, D.V. and Rattray M. (1965) Gravitational circulation in straits and estuaries. J. Mar. Res. 23: 104–22
 ↑ Pritchard, D.W. (1952) Salinity distribution and circulation in the Chesapeake Bay estuaries system. J. Mar. Res. 11: 106–23
 ↑ Geyer, W.R., Trowbridge. J.H. and Bowen, M.M. (2000) The dynamics of a partially mixed estuary. J. Phys. Oceanogr. 30: 2035–48
 ↑ MacCready, P. (2004) Toward a unified theory of tidallyaveraged estuarine salinity structure. Estuaries 27: 561–70
 ↑ Munk, W.H. and Anderson, E.R. (1948). Notes on the theory of the thermocline. J. Marine Res. 7: 276295
 ↑ Van Aken, H. M. (1986). The onset of seasonal stratification in shelf seas due to differential advection in the presence of a salinity gradient. Cont. Shelf Res. 5: 475485
 ↑ Simpson, J.H., Brown, J., Matthews, J.and Allen, G. (1990) Tidal straining, density currents, and stirring in the control of estuarine stratification. Estuaries 13: 125132
 ↑ Jay, D.A. and Musiak, J.D. (1996) Internal tidal asymmetry in channel flows: origins and consequences. In Mixing Processes in Estuaries and Coastal Seas, ed. C Pattiaratchi, Am. Geophys. Union, Washington: 211–49
 ↑ Scott, C.F. (1994) A numerical study of the interaction of tidal oscillations and nonlinearities in an estuary. Estuarine Coast. Shelf Sci. 39: 477–96
 ↑ Prandle, 2004 Saline intrusion in partially mixed estuaries. Estuarine, Coastal and Shelf Science 59: 385397
 ↑ Burchard, H. and Baumert, H. (1998) The formation of estuarine turbidity maxima due to density effects in the salt wedge. A hydrodynamic process study. J. Phys. Oceanogr. 28: 209–321
 ↑ Cheng P., De Swart, H.E. and ValleLevinson, A. (2013) Role of asymmetric tidal mixing in the subtidal dynamics of narrow estuaries. J. Geophys. Res. 118: 2623–2639
 ↑ De Nijs, M.A.J., Pietrzak, J.D. and Winterwerp, J.C.(2011) Advection of the Salt Wedge and Evolution of the Internal Flow Structure in the Rotterdam Waterway. J. Phys. Oceanogr. 41: 327
 ↑ Nepf, H.M. and Geyer, W.R. (1996) Intratidal variations in stratification and mixing in the Hudson estuary. J. Geophys. Res. 101: 1207912086
 ↑ Nunes, R.A. and Simpson, J.H. (1985) Axial Convergence in a Wellmixed Estuary. Coastal and Shelf Science 20: 637649
 ↑ Smith, R. (1976) Longitudinal dispersion of a buoyant contaminant in a shallow channel. J. Fluid Mech. 78: 677–88
 ↑ Burchard, H. and Schuttelaars, H.H. (2012) Analysis of Tidal Straining as Driver for Estuarine Circulation in WellMixed Estuaries. J. Phys. Oceanograph. 42: 261271
 ↑ Burchard, H., Hetland, R.D., Schulz, E. and Schuttelaars, H.H. (2011)Drivers of Residual Estuarine Circulation in Tidally Energetic Estuaries: Straight and Irrotational Channels with Parabolic Cross Section. J.Phys.Oceanogr. 41: 548570
 ↑ Huijts, K.M.H., Schuttelaars, B.H.M., De Swart, H.E. and Friedrichs, C.T. (2009) Analytical study of the transverse distribution of alongchannel and transverse residual flows in tidal estuaries. Continental Shelf Res. 29: 89100
 ↑ Basdurak, N.B. and ValleLevinson, A. (2013) Tidal variability of lateral advection in a coastal plain estuary. Continental Shelf Res. 61–62: 85–97
Please note that others may also have edited the contents of this article.
