5.2 Groundwater flow and Darcy's Law

Groundwater flow is governed by Darcy's Law, which links flow rate to hydraulic gradient and conductivity. This principle is foundational for understanding how water moves through aquifers and directly shapes how we manage water resources and assess contamination risks.

Factors like hydraulic gradient, conductivity, and aquifer geometry influence groundwater flow. Together, they determine flow direction, speed, and patterns across different geological settings.

Groundwater Flow Principles and Darcy's Law

Principles of Darcy's Law

Darcy's Law describes how groundwater moves through porous media like soil, sand, and rock. Henry Darcy developed it in 1856 through experiments on water flow in sand filters used for water purification in Dijon, France.

The law establishes a linear relationship between groundwater flow rate and the hydraulic gradient, with hydraulic conductivity as the proportionality constant:

$Q = -KA\frac{dh}{dl}$

• $Q$ = volumetric flow rate (e.g., $m^3/day$ or gallons per minute)
• $K$ = hydraulic conductivity, how easily water moves through the material (e.g., $m/day$ or $ft/day$)
• $A$ = cross-sectional area perpendicular to flow direction (e.g., $m^2$)
• $\frac{dh}{dl}$ = hydraulic gradient, the change in hydraulic head over distance in the flow direction (dimensionless)

The negative sign indicates that flow moves in the direction of decreasing hydraulic head. In practice, you'll often see the equation written with the negative sign absorbed so that $Q$ comes out positive when head drops along the flow path.

Darcy's Law relies on several simplifying assumptions:

• Laminar flow with low Reynolds numbers, meaning smooth, orderly flow without turbulence. This holds true for most natural groundwater systems but breaks down in very coarse gravels or fractured rock with wide apertures.
• Fully saturated porous media, where all pore spaces are filled with water (no air pockets).
• Homogeneous and isotropic aquifer properties, meaning $K$ is constant everywhere and the same in all directions.
• Incompressible fluid and aquifer material, meaning densities stay constant under typical groundwater pressure changes.

Real aquifers rarely satisfy all these assumptions perfectly, but Darcy's Law still provides an excellent approximation for most practical problems.

Calculations with Darcy's Law

Calculating flow rate ($Q$):

To find the volumetric flow rate, plug in your values for $K$, $A$, and $\frac{dh}{dl}$:

$Q = -KA\frac{dh}{dl}$

For example, if $K = 10 \; m/day$, $A = 50 \; m^2$, and the hydraulic head drops 2 m over a 500 m distance (so $\frac{dh}{dl} = -\frac{2}{500} = -0.004$):

$Q = -(10)(50)(-0.004) = 2 \; m^3/day$

Calculating groundwater velocity ($v$):

The Darcy velocity (also called specific discharge) is $q = Q/A$, but this isn't the actual speed water travels through pores. The real seepage velocity accounts for the fact that water can only move through the connected pore spaces:

$v = \frac{Q}{n_e A}$

Effective porosity ($n_e$) represents only the interconnected pore space available for flow. Total porosity includes dead-end and isolated pores that don't contribute to transport. Because $n_e$ is always less than 1, the seepage velocity is always faster than the Darcy velocity.

Hydraulic conductivity in more detail:

$K$ depends on both the aquifer material and the fluid flowing through it:

$K = \frac{k \rho g}{\mu}$

• $k$ = intrinsic permeability of the aquifer material ($m^2$), a property of the pore structure alone
• $\rho$ = fluid density ($kg/m^3$)
• $g$ = gravitational acceleration ($9.81 \; m/s^2$)
• $\mu$ = dynamic viscosity of the fluid ($Pa \cdot s$)

This separation matters because intrinsic permeability ($k$) is purely a rock/sediment property, while $K$ changes if the fluid changes. For most groundwater problems you'll work with $K$ directly, but the distinction becomes important when dealing with fluids other than freshwater (e.g., saltwater intrusion or contaminant plumes with different densities).

Hydraulic Gradient in Groundwater

The hydraulic gradient is defined as:

$i = \frac{dh}{dl}$

where $dh$ is the change in hydraulic head over a distance $dl$ measured along the flow direction.

Hydraulic head itself has two components:

• Elevation head: height of the measurement point above a reference datum
• Pressure head: water pressure at that point divided by the specific weight of water ($\rho g$)

Groundwater always flows from higher hydraulic head toward lower hydraulic head. The steeper the gradient (larger head change over a shorter distance), the faster the flow. Think of it like slope on a hill: a steeper slope means water runs faster.

In the field, you measure hydraulic head using monitoring wells (piezometers) installed at different locations. With three or more wells, you can determine both the magnitude and direction of the gradient using triangulation methods.

Factors Affecting Groundwater Flow

Hydraulic gradient

• Groundwater flows in the direction of decreasing hydraulic head.
• Steeper gradients produce higher flow rates. Near pumping wells, for instance, the gradient steepens dramatically, creating a cone of depression that draws water inward.

Hydraulic conductivity

• Higher $K$ values mean easier flow and higher flow rates.
• $K$ depends on grain size, sorting, and porosity of the aquifer material, plus fluid properties like density and viscosity.
• Well-sorted, coarse-grained sediments like clean gravel can have $K$ values of $10^2$ to $10^3 \; m/day$, while clay may be as low as $10^{-7} \; m/day$. That's a billion-fold difference, which is why aquifer material matters so much.

Aquifer heterogeneity and anisotropy

• Real aquifers vary in composition from place to place (heterogeneity), creating preferential flow paths and non-uniform flow patterns.
• Anisotropy means $K$ differs depending on direction. Layered sedimentary rocks typically have much higher conductivity parallel to bedding planes than perpendicular to them, sometimes by a factor of 10 or more. This can redirect flow away from what the hydraulic gradient alone would predict.

Recharge and discharge areas

• Recharge areas (zones of infiltration from precipitation, streams, or irrigation) add water and increase hydraulic head, driving flow downward and outward.
• Discharge areas (springs, pumping wells, gaining rivers) remove water, lowering hydraulic head and drawing groundwater toward them.

Geologic structures and boundaries

• Faults, folds, and confining layers (aquitards) can act as barriers or conduits depending on their properties. A fault filled with crushed, permeable rock may channel flow, while one sealed with clay gouge blocks it.
• Boundary conditions constrain flow at aquifer edges. An impermeable rock contact creates a no-flow boundary, while a large lake or river can act as a constant-head boundary. These boundaries strongly influence flow patterns near the margins of an aquifer system.