Saltwater Intrusion: Pumping a Coastal Aquifer Out of Balance

Two fluids in equilibrium under a coastline

Fresh groundwater is lighter than seawater (densities of about 1.000 vs 1.025 g/cm³), so under a coastline the two don't mix cleanly — fresh water floats in a lens on top of a wedge of denser salt water. In 1888-1901, W. Badon Ghyben and Alexander Herzberg independently worked out where the interface between them has to sit for that stack to be in hydrostatic equilibrium: the depth of the interface below sea level is about 40 times the height of the water table above sea level.

z = (rho_f / (rho_s - rho_f)) * h

rho_f = 1.000  (fresh water)     rho_s = 1.025  (seawater)
=> z ≈ (1.000 / 0.025) * h  ≈  40 * h

// a 1-metre-high freshwater lens implies ~40 m of fresh water
// sitting above the saltwater interface

That 40:1 ratio is the reason coastal aquifers are so fragile: dropping the water table by just half a metre through pumping can pull the interface up by roughly 20 metres — a small, easily overlooked change at the surface produces a disproportionate response underground, purely because of how thin the density difference between the two fluids is.

Upconing: how a single well tips the balance

Pump a well that penetrates the freshwater lens and the drawdown doesn't just lower the water table locally — it locally reduces the freshwater pressure holding the interface down, and the saltwater wedge rises beneath the well in a cone, a process called upconing. Pump too hard or too long and the cone can reach the well screen, at which point the well starts drawing salt water directly and often cannot be flushed clean again, since the interface below a screen rarely resets to its original shape once breached.

Sea-level rise and abstraction: a two-way squeeze

The Ghyben-Herzberg balance is set by the height of the freshwater table relative to sea level, so anything that raises the sea or lowers the water table pushes the interface inland and upward. Groundwater pumping does the second directly; sea-level rise does the first, and also narrows the coastal freshwater lens by simple geometry as the shoreline itself moves inland. Low-lying islands and deltas — where the natural freshwater lens is thin to begin with — are hit by both mechanisms simultaneously.

Where the simple picture breaks down

Ghyben-Herzberg assumes a sharp interface and hydrostatic (non-moving) fluids, which is a good first approximation but not the real story: in practice there's a diffusive mixing zone tens of metres thick where fresh and salt water blend by dispersion, and active recharge, tides and pumping keep the whole system in transient disequilibrium rather than the clean static balance the formula assumes. Numerical models like SEAWAT solve the coupled flow and salt-transport equations directly for a realistic mixing zone; Ghyben-Herzberg remains the fast, back-of-envelope check every one of those models is validated against.

Frequently asked questions

Try it live

Everything above runs in your browser — open Saltwater Intrusion and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.