How Water Chooses Its Path Through Soil express

Original: "The Statistical physics of unsaturated soil water: kinetic theory and non commutative pore water dynamics"
· Riccardo Rigon
arXiv:2607.09416 · 2026-07-10 · CC BY 4.0 · 2 min
Physicists have derived an equation that explains why water in the ground sometimes flows slowly and sometimes seeps through in swift rivulets.

📄 Showing the "Simple" version — "Popular" is not ready yet. Add it to favorites to help prioritize it.

Soil is not a solid sponge but a labyrinth of pores of various sizes. After rain, water does not fill them all at once: first the narrowest ones, then the wider ones. Physicists built a theory of this process based on a simple image — a set of communicating vessels. The flows between pores obey the laws of viscous fluid flow, and the driving force is the difference in chemical potentials. Roughly speaking, it is water's "desire" to occupy a particular pore, taking into account capillary forces and gravity.

Water in a pore behaves like a capricious visitor: if a neighboring pore is slightly more spacious and energetically favorable, it moves there, but not instantly — with a delay depending on the channel radius.

The key turned out to be the dimensionless Damköhler number for each pore. When it is small, water manages to distribute evenly: this is the classic slow absorption described by the Richards equation. When the number is large, water darts through selected paths without having time to "look around." This gives rise to fast preferential flows, familiar to anyone who has seen a puddle disappear in minutes rather than hours after a downpour. The theory shows that these are not two different phenomena but a smooth transition governed by a single parameter.

Hysteresis, when the drying curve does not coincide with the wetting curve, here is not a property of individual pores but a geometric phase. It arises because the system traverses the same path differently in the forward and reverse directions. The authors derive a simple rule for this: the area of the loop is proportional to the square of the accumulated flux.

This view makes dozens of ad hoc models unnecessary. All the necessary information about the soil is the geometry of the pore network, which can be measured with a microtomograph. There is no need to fit parameters to field data. The equation itself decides where water should flow slowly and where quickly.

The Damköhler number, which governs water behavior in soil, originally came from chemical engineering and described how fast a chemical reaction proceeds relative to the flow of matter. Now it also "works" right under our feet.

\partial_t g + \nabla \cdot \mathbf{F} = C[g] - E - T
The change in the fraction of filled pores over time (∂_t g) plus the fluxes between spatial points (∇·F) equals the difference between sources and sinks: C is the inflow from neighboring pore classes, E is evaporation, T is root uptake.
Scientists
Carl Friedrich GaussGeorge Gabriel StokesHenri PoincaréDaniel BernoulliArchimedesAndrey Markov
Concepts
soil moisture fluid dynamics in porous media porosity Hooke's law Navier-Stokes equations reynolds number hydrodynamics virial theorem Poiseuille's equation Water liquid water Goldman-Hodgkin-Katz equation Starling's law precipitation ergodic theorem free energy random walk ground state drag porous material