3 Quantifying Groundwater Flow in The Presence of Density Variations
The description of groundwater flow in variable-density systems is considerably more complex than in constant-density systems. While field measurements of hydraulic head conveniently allow for an assessment of groundwater flow direction and magnitude for constant-density systems, hydraulic head measurements cannot be used for flow calculations without making allowance for buoyancy effects in variable-density systems. Without citing anyone or any literature, it was casually remarked in a Dutch report (DZRD, 1936) that “It is known that the measured head of water with a high chloride content requires a correction” [paraphrased translation from Dutch by the authors]. Nevertheless, the use of hydraulic head measurements in a variable-density groundwater system still leads to confusion or even misinterpretation to this day.
Some basic intuition about the effect of density on the flow field can be gained by considering a confined aquifer with an impermeable base in which fresh water and salt water are separated by a vertical interface at time zero (Figure 5a).
All the system boundaries are closed, so there can be no inflow or outflow of water. Under the prevailing hydrostatic conditions at time t = 0, the groundwater pressure p increases with depth according to Equation 1.
(1) |
where:
p | = | pressure (M/(LT2)), e.g., kg m−1 s−2 |
p0 | = | pressure at the top of the aquifer where z = 0 (M/(LT2)), e.g., kg m−1 s−2 |
ρ | = | groundwater density (M/L3), e.g., kg m−3 |
g | = | gravitational acceleration (L/T2), e.g., m s−2 |
z | = | elevation relative to a datum, which in this case is the top of the aquifer where z = 0 and z is positive upward, thus there is a minus sign in the equation because p increases with depth (L), e.g., m |
Assume that the groundwater density on the freshwater side of the interface is ρf = 1000 kg m−3 and on the saltwater side it is ρs = 1025 kg m−3. The pressure thus increases faster with depth on the saltwater side than on the freshwater side (Figure 5a). The horizontal pressure difference gives rise to a flow, which in this example will be from the right to the left. Because the pressure difference between the freshwater and saltwater domains increases with depth, the initial flow will have the greatest magnitude at the base of the aquifer, and the interface will thus tend to rotate. As a result, the flow of fresh water above the interface will have to be from left to right to compensate for the displacement of the salt water along the aquifer base from right to left. The flow of the fresh water is thus in the opposite direction as the salt water (Figure 5b), resulting in a clockwise rotation of the initially vertical interface.
In reality, the description of the flow in an aquifer with a non-horizontal interface is far more complex (Bakker et al., 2004; Verruijt, 1980), but this simple example demonstrates that the presence of density differences gives rise to rotational flow. If this is the case, the hydraulic head is no longer suitable to describe the groundwater flow. The formal proof for this, which relies on determining whether the so-called curl (a property of a vector field that indicates rotational movement) of the force field is zero, can be found in Hubbert (1957).
The fact that the flow field is not irrotational means that measurements of hydraulic head cannot be directly compared to infer flow directions or magnitudes. Instead, flow calculations must be based on the general form of Darcy’s law, which is shown in Equation 2 (Bear, 1972).
(2) |
where:
= | specific discharge (L/T), e.g., m s–1 | |
k | = | intrinsic permeability (L2), e.g., m2 |
μ | = | dynamic viscosity of the groundwater (M/(LT)), e.g., kg m−1 s−1 |
= | gradient operator and represents the rate of change of a variable (p in this case) per unit of distance in the x, y and z direction | |
= | gravitational acceleration vector (L/T2), e.g., m s−2 |
When ρ is spatially constant, Equation 2 may be written as Equation 3.
(3) |
In the expression after the third equal sign, the parameters before the gradient operator were replaced by K (m s−1) the hydraulic conductivity. The term within the last set of parentheses is the hydraulic head h (m) as shown in Equation 4.
(4) |
where:
hp | = | pressure head (L), e.g., m |
The first term in Equation 4, z, is the elevation head, which determines the groundwater’s potential energy (per unit of weight) due to its position in the Earth’s gravitational field. In practice it is simply the elevation of the piezometer screen relative to a standard datum, usually mean sea level. The pressure head is a measure of the energy of groundwater because of its pressure above a reference pressure. The value of p is gage pressure, for which the atmospheric pressure is taken as the reference pressure. Inserting Equation 4 into Equation 3 gives Equation 5.
(5) |
Equation 5 is the familiar form of Darcy’s law based on the hydraulic head gradient. However, given that the step from Equation 2 to 3 relies on the assumption that the density is spatially constant, the hydraulic head gradient is not suitable to quantify groundwater flow in variable-density groundwater flow systems. The hydraulic head form of Darcy’s law (Equation 5) is strictly applicable to constant density systems.