Salt wedge estuaries

From Coastal Wiki
Jump to: navigation, search
Definition of Salt wedge:
Seawater intrusion in an estuary as a wedge-shaped bottom layer which hardly mixes with the overlying fresh water layer. Salt wedges occur in estuaries where tidal motion is very weak or absent.
This is the common definition for Salt wedge, other definitions can be discussed in the article

Introduction

In estuaries where tidal motion is very weak or absent, seawater intrudes in the form of a wedge along the bottom, as a consequence of the higher density of seawater compared to fresh water. The intrusion length depends essentially on two parameters: river flow velocity [math]v[/math] and water depth [math]h[/math]. This article deals with the hydrodynamic balances that determine the salt wedge, illustrated by a model for a prismatic estuary (uniform depth and width). We follow the derivation first given by Schijf and Schönfeld in 1953 [1].

In most estuaries, seawater intrusion is strongly influenced by tidal motion (see Estuarine circulation). Even in microtidal estuaries this influence is often significant [2]. In some cases, especially in the Mediterranean Sea, tides are so small that their influence on seawater intrusion is minor. This article focuses on such situations; we assume that the salt wedge is not significally altered during the slight tidal motion up and down the estuary.

Seawater intrusion process

Figure 1: Longitudinal section of river with a fictitious corrugated sloping bottom. The sea boundary is at [math]x=0[/math]. The river discharge is very small; the upper boundary of the salt wedge is almost horizontal. The salt intrusion length [math]L[/math] corresponds to the location where the level of the river bed equals the sea level.

First consider the case of a very small river discharge such that the stress it exerts on the salt wedge is not significant. The absence of forces leads to an equilibrium situation where the salt wedge slope is zero. The level of the river bed rises inland; at some point inland, the level is equal to the sea level. Seawater can enter the river that far. The salt intrusion length is determined by the intersection of the sea level with the river bed, see Fig. 1.

In general, the shear stress exerted by the river on the intruding seawater cannot be neglected. We consider a situation without tide or a very small tide (tidal flow velocities much less than the flow velocity due to river discharge) and constant river discharge. We then have a stationary salt wedge, also called 'arrested salt wedge'. Except for extreme river discharges, the rate at which river water flows over the seawater wedge is not fast enough to produce complete mixing of the water column. However, the river flow velocity is high enough to produce a thin mixing layer. The river water exerts a shear stress on the seawater wedge via this mixing layer. The surface of the seawater wedge is therefore inclined (the pycnocline), such that the hydrostatic force (due to the slopes of the salt wedge and the water surface) balances the shear stress, see Fig. 2. The thin mixing layer is generally turbulent; the shear stress is therefore proportional to the square of the velocity difference over the mixing layer and it is assumed that it can be parameterized as the product of the river flow velocity squared and a shear stress coefficient [math] C_i [/math]. In theory, this coefficient should have a certain fixed value, but experimental observations yield various values.

Figure 2: Longitudinal section of a salt wedge in a prismatic estuary. The sea boundary is situated at [math]x[/math]=0. The shape of the salt wedge shown in the figure corresponds to the theoretical curve (solution of Eq. (1)) for 8 m water depth, river flow velocity 0.3 m/s, salt layer density 1020 kg/m3 (internal river Froude number [math]F_i[/math]=0.24) and interfacial shear stress coefficient [math]C_i[/math]=0.00025.

The reason is that the ideal situation of a stationary salt wedge hardly occurs in practice. The river discharge varies, the tide is not negligibly small, extra turbulence is caused by a non-uniform river geometry, by waves, wind and shipping. Numerical models can partly quantify these influences, but it is difficult to accurately predict the salt intrusion length.


Seawater intrusion length

If the tide has a significant influence on the seawater intrusion process, an estuarine circulation is generated; this situation is discussed in the article Estuarine circulation. Here we assume the simpler case that the influence of the tide can be neglected. In the special case that the river geometry can be approached by a prismatic channel (uniform depth and width, no bottom slope). The relative depth of the upper fresh water layer [math]y=h_1/h[/math] then satisfies the differential equation (see Appendix)

[math] \large\frac{d}{dx} \normalsize [ -4 y^5 + 5y^4 - 10 F_i^2 y(2-y) ] = 20 \large \frac{C_i F_i^2}{h} \normalsize , \qquad (1) [/math]

where the internal Froude number [math]F_i[/math] is given by

[math]F_i = \large \frac{Q}{hb} \normalsize [ \epsilon g h]^{-1/2} , \qquad (2)[/math]

[math]Q[/math]= river discharge, [math]h[/math]= river depth, [math]b[/math]= river width and [math]\epsilon=\large \frac{\Delta \rho}{\rho} \normalsize [/math]= relative density difference of seawater and river water. The differential equation (1) is easily integrated if conditions at both boundaries of the salt wedge are known.

The conditions at the tip of the salt wedge is obvious: [math]y=1[/math] at [math]x=L [/math].

The condition at the sea boundary is more subtle. The upper layer depth cannot be zero because no fresh water could be discharged to the sea. When during the process of seawater intrusion the depth of the upper layer at the sea boundary decreases, the outflow velocity [math]v_1[/math] increases. However, the intrusion process stops when the outflow velocity equals the critical velocity [math]v_c[/math]. This is because a further upward shift of the interface boundary at sea cannot propagate upstream into the river. The critical velocity for upstream propagation of a perturbation at the interface is given by

[math]v_c=[\large \frac{\Delta \rho}{\rho} \normalsize g h_1]^{1/2} . \qquad (3) [/math]

Substitution in Eq. (2) gives [math]F_i = \large \frac{h_1 v_c}{h} \normalsize [\large \frac{\Delta \rho}{\rho} \normalsize g h]^{-1/2} = \large(\frac{h_1}{h})^{3/2}\normalsize = y(0)^{3/2}[/math],

from which it follows that the boundary conditions at the sea boundary [math]x=0[/math] is given by [math]y(0) = F_i^{2/3} . \qquad (4)[/math]

Integration of the differential equation with the two boundary conditions gives the length [math]L[/math] at constant river discharge:

[math]L = \Large \frac{h}{20 C_i} \normalsize (-10+15F_i^{2/3}-6F_i^{4/3}+F_i^{-2}) . \qquad (5)[/math]

In this expression, the salt intrusion length is inversely proportional to the interfacial shear stress coefficient [math]C_i[/math]. The value of this coefficient cannot be determined via the model and must therefore be determined independently.

The dimensionless salt wedge length [math]C_i L / h[/math] as a function of the internal river Froude number [math]F_i[/math] is shown in Fig. 3. The condition at the sea boundary plays a crucial role for the salt wedge length. The depth of the river outflow layer [math]h_1(0)[/math] is proportional to the river discharge at the power 2/3, according to the equations (2) and (4). Therefore the depth of the salt wedge layer at the mouth, [math]h_2(0)[/math], becomes very small for high river discharge.

Integration of Eq. (1) also gives the shape of the seawater wedge, shown in Fig. 2 for particular values of river depth, river discharge and shear stress coefficient [math] C_i [/math].

Generally, the analytical model cannot be used for predicting salt intrusion at low river discharge. This is an important limitation, because precisely this information is often most relevant in practice. The problem is that there is no analytical solution to the equations in the case of a sloping river bed. At low river runoff, the seawater intrudes thus far that the influence of the bottom slope cannot be neglected. In that case, the equations describing the seawater intrusion must be solved numerically. For rivers with a strongly sloping river bed, a first estimate of the salt intrusion length at low river discharge is obtained by the intersection of the seawater level with the river bottom, as shown in Fig. 1.

Additional mixing of river water and sea water can occur as a consequence of local river contractions, sharp river bends and high bottom sills. In these cases too, the analytical solution is less reliable and it is better to use a numerical model to determine the seawater intrusion.


Comparison with observations

Figure 3: Dimensionless salt wedge length [math]C_i L / h[/math] according to Eq. (5) for a prismatic estuary and observed salt wedge lengths in the Rhône, Ebro [3], the Mississippi [4], Magdalena[5] and Neretva[6] estuaries at different discharges (corresponding to different internal river Froude numbers [math]F_i[/math]). For all estuaries the same value is used for the interfacial shear stress coefficient, [math]C_i[/math]=0.00025. Other parameters are indicated in Table 1.

The Rhône (France) and the Ebro (Spain) are examples of estuaries with almost no tidal motion (both situated on the Mediterranean coast, tidal amplitude in the order of 0.1 m on average). The tidal amplitude for the Mississippi at the US Gulf coast is in the order of 0.25 m on average; the tidal current velocities are small, but not much smaller than the river flow velocities for the lowest discharges considered in the comparison. The Magdalena (Colombia) discharges into the Caribbean Sea with a predominantly diurnal, mixed microtidal regime. The tidal range of 0.13 m to 0.40 m produces small tidal velocities that shift the salt water front over a distance of about 300 m[7].

Figure 3 shows the observed dimensionless salt wedge lengths [math]C_i L / h[/math] in these estuaries for different internal river Froude numbers [math]F_i[/math]. The observed lengths are compared with the theoretical result for a prismatic estuary (Eq. (5)), where the same interfacial shear stress coefficient is used for all the estuaries. For very small river flow the model lengths are much larger than the observed intrusion lengths. This is due in particular to the presence of high sills in the upstream part of the Rhône estuary (at 35 km from the mouth) and Ebro estuary (at 32 km from the mouth). Other less pronounced sills situated downstream also tend to arrest the salt wedge for a certain range of river discharges; the salt wedge length [math]L[/math] therefore does not increase continuously as a function of [math]F_i[/math], as would be the case for a true prismatic canal. The reasonable correspondence between observed and computed salt intrusion lengths for low-to-mean river discharges does not imply that the value of the interfacial shear stress coefficient, [math]C_i=0.00025[/math], should hold in general. In the present comparison this coefficient includes the effect of several factors by which the analytical model deviates from the real estuaries, such as deviations from the prismatic geometry (bed slope, sills, etc.) and tidal influence.


Appendix 1: Analytical model for a prismatic estuary

The following is a simple two-layer description of seawater intrusion in a microtidal prismatic (uniform width and depth) estuary. The salt wedge is represented by a separate layer of thickness [math]h_2[/math] along the bottom of the estuary, see figure 2. The total water depth [math]h[/math] is the sum of the bottom layer thickness [math]h_2[/math] and the upper surface layer thickness [math]h_1[/math]. Lateral variations of depth, density (salinity) and velocity are ignored. The velocity is indicated by [math]v[/math] and the density by [math]\rho[/math]. Subscripts 1 and 2 indicate the surface layer and the bottom layer respectively.

The pressure gradient [math]d p_2 / d x[/math] in the lower layer is related to the slope of the salt wedge (relative to a horizontal reference plane) and to the pressure gradient exerted by the overlying fresh water,

[math]\Large \frac{d p_2}{d x} \normalsize = g \rho_2 ( \beta + \Large \frac{d h_2}{d x} \normalsize )+ g \rho_1 \Large \frac{d h_1}{d x} \normalsize , \qquad (A1) [/math]

where [math]\rho_1=\rho[/math] is the fresh water density, [math]\rho_2=\rho+ \Delta \rho[/math] the seawater density and [math]\beta[/math] the bottom slope. Other symbols are shown in Fig. 2. Because of the lower density of fresh water compared to seawater, this pressure gradient generally acts in upstream direction. A counteracting pressure gradient is produced by friction at the interface (pycnocline) induced by the overlying river flow. The shear stress [math]\tau_i[/math] exerted by river flow on the salt wedge is parameterized as [math]\tau_i / \rho = C_i v_1^2 [/math], where [math]C_i[/math] is the interfacial shear stress coefficient of order [math]10^{-4}-10^{-3} .[/math][8]

The salt wedge intrusion is also influenced by bottom friction. However, in a stationary situation bottom friction does not play a role because the salt wedge is at rest.

River flow is confined to the upper layer which becomes thinner in downstream direction. The resulting spatial acceleration of river flow creates a gradient in the hydraulic head given by [math]\rho v_1 d v_1 / d x .[/math] The surface layer is further subjected to a pressure gradient related to the water surface slope (relative to a horizontal reference plane),

[math]\Large \frac{d p_1}{d x} \normalsize = g \rho_1 ( \beta + \Large \frac{d h_2}{d x} + \frac{d h_1}{d x} \normalsize) , \qquad (A2) [/math]

and to the interfacial shear stress [math]\tau_i .[/math] The stationary hydrodynamic balance of the two layers is therefore given by the equations

[math]\Large \frac{d p_2}{\rho_2 d x} + \frac{\tau_i}{\rho h_2} \normalsize = 0; \; \Large \frac{d p_1}{\rho_1 d x} \normalsize + v_1 \; \Large \frac{d v_1}{d x} - \frac{\tau_i}{\rho h_1} \normalsize = 0 . \qquad (A3)[/math]

In a prismatic estuary the depth [math]h[/math] and width [math]b[/math] do not vary along the channel axis. By subtracting the two layer equations the salt wedge slope [math]d h_2 / d x[/math] and the bottom slope [math]\beta[/math] are eliminated. By assuming that the bottom slope is not much larger than the surface slope (water depth [math]h[/math] approximately constant), we can substitute [math]h_2 \approx h-h_1[/math]. The velocity in the upper layer [math]v_1[/math] can be eliminated from the equations by using the condition of continuity of river discharge per unit width,

[math]h_1 v_1 = h v = Q / b, \qquad (A4) [/math]

where [math]Q[/math] is the river discharge. These manipulations lead to the ordinary nonlinear differential equation for the variable [math]y=h_1 / h[/math],

[math] (y^3-F_i^2)(1-y) \large\frac{dy}{dx} \normalsize = \large \frac{C_i F_i^2}{h} \normalsize , \qquad (A5) [/math]

with [math]F_i = \large \frac{Q}{hb} [ \frac{gh \Delta \rho}{ \rho}]^{-1/2} \normalsize .[/math]

This is equivalent to Eq. (1) indicated above.


Appendix 2: Parameters of a few salt wedge estuaries

Table 1. Data of observed salt wedge lengths for the Rhône (France)[3] and Ebro (Spain)[3], the Mississippi (US)[4], Magdalena (Colombia)[5] and Neretva (Croatia)[6] estuaries. The indicated depths are not mean depths, but the smallest depths along the thalweg below the salt wedge taking into account natural sills.
Estuary depth [m] cross section [m2] river discharge [m3/s] observed salt wedge intrusion length [math]L[/math] [km] [math]L[/math] from Eq. (5), [math]C_i=2.5 \, 10^{-4}[/math]
Rhône 8 2500 500 35 53
8 2500 600 28 34
8 3000 850 20 23
8 3000 1200 9 9
Ebro 5.5 1000 80 32 177
5.5 1000 250 18 13
5.5 1000 400 5 3
Mississippi main channel 18 14220 2830 220 296
18 14220 4250 140 120
Mississppi southwest pass 14 6440 2550 32 35
14 6440 2830 23 26
12 5520 6120 1.5 0.7
Magdalena river 12 6000 4000 7 7
12 6000 8000 0 0.1
Neretva 8 800 300 20 15
8 800 400 10 6.5
8 800 500 5 2.9


Related articles

Estuarine circulation
Seawater density
Seawater intrusion and mixing in estuaries


References

  1. Schijf, J.B., Schönfeld, J.C. 1953. Theoretical considerations on the motion of salt and fresh water. Proceedings of the 5th I.A.H.R. Congress, 321–333
  2. Geyer, R.W. and Farmer, D.H. 1989. Tide-induced variation of the dynamics of a salt wedge estuary. J. Phys. Oceanography 19: 1060-1072
  3. 3.0 3.1 3.2 Ibanez, C., Pont, D. and Prat, N. 1997. Characterization of the Ebre and Rhone estuaries: A basis for defining and classifying salt-wedge estuaries. Limnol. Oceanogr. 42: 89–101
  4. 4.0 4.1 Rattray, M. Jr and Mitsuda, E. 1974. Theoretical analysis and conditions in a salt wedge. Estuarine and Coastal Marine Science 2: 375-394
  5. 5.0 5.1 Arevalo, F.M., Alvarez-Silva, O., Caceres-Euse, A. and Cardona, Y. 2022. Mixing mechanisms at the strongly-stratified Magdalena River’s estuary and plume. Estuarine, Coastal and Shelf Science 277, 108077
  6. 6.0 6.1 Krvavica, N., Gotovac, H. and Loncar, G. 2021. Salt-wedge dynamics in microtidal Neretva River estuary. Regional Studies in Marine Science 43, 101713
  7. Ospino, S., Restrepo, J.C., Otero, L., Pierini, J. and Alvarez-Silva, O. 2018. Saltwater intrusion into a river with high fluvial discharge: a microtidal estuary of the Magdalena River, Colombia. J. Coast Res. 34: 1273–1288
  8. Abraham, G., Karelse, M. and Van Os, A.G. 1979. On the magnitude of the interfacial shear stress of subcritical stratified flows in relation with interfacial stability. Hydr. Res. 17: 273-287
The main author of this article is Job Dronkers
Please note that others may also have edited the contents of this article.

Citation: Job Dronkers (2022): Salt wedge estuaries. Available from http://www.coastalwiki.org/wiki/Salt_wedge_estuaries [accessed on 6-12-2024]