Le sol n'est pas une éponge homogène, mais un labyrinthe de pores de différentes tailles. Après une pluie, l'eau ne les remplit pas tous d'un coup : d'abord les plus étroits, puis les plus larges. Les physiciens ont construit une théorie de ce processus en s'appuyant sur une image simple : un ensemble de vases communicants. Les échanges entre les pores obéissent aux lois de l'écoulement d'un fluide visqueux, et la force motrice est la différence de potentiel chimique. En gros, c'est le « désir » de l'eau d'occuper tel ou tel pore, compte tenu des forces capillaires et de la gravité.
Le nombre sans dimension de Damköhler pour chaque pore s'est avéré crucial. Lorsqu'il est petit, l'eau a le temps de se répartir uniformément — c'est l'infiltration lente classique décrite par l'équation de Richards. Lorsque le nombre est grand, l'eau file par des chemins privilégiés, sans avoir le temps de « regarder autour d'elle ». C'est ainsi qu'apparaissent les écoulements rapides, que connaît bien quiconque a vu une flaque disparaître en quelques minutes après une averse, et non en quelques heures. La théorie montre qu'il ne s'agit pas de deux phénomènes distincts, mais d'une transition progressive contrôlée par un seul paramètre.
Cette perspective rend superflues des dizaines de modèles particuliers. Toute l'information nécessaire sur le sol — la géométrie du réseau de pores, mesurable par microtomographie — suffit. Pas besoin d'ajuster des paramètres sur des données de terrain. L'équation décide elle-même où l'eau doit s'écouler lentement et où rapidement.
tel que rédigé par les auteurs, en anglais arXiv:2607.09416
We develop a statistical-mechanical theory of water in unsaturated soil whose outcome is a continuum field equation for the pore-occupancy g(r,x,t), the fraction of pores of radius r that are water-filled at position x and time t. The theory is built across three scales: microscopic inter-pore transfers set by Hagen-Poiseuille rates and a driving potential (the difference of pore-class chemical potentials, taken in capillary-gravitational form but open to adsorptive, osmotic, or thermal refinement); a mesoscale master equation relaxing the occupancy toward the equilibrium step g_eq=H(r*-r); and, on contracting the averaging volume to a point, the continuum balance d_t g + div F = C[g] - E - T, of which everything else is a limit, a moment, or a boundary resolution. The kinetic equation is an Onsager gradient flow descending the Gibbs free energy, with an H-theorem for the isothermal unforced system and mass conservation as its zeroth moment. A single dimensionless group, the pore-resolved Damkohler number Da(r,x), organizes the behavior and unifies phenomenologies long modelled separately. A Chapman-Enskog reduction identifies Richards' equation as the quasi-static (Da->0) limit, with matric potential and hydraulic conductivity K emerging only there and K vanishing below the percolation threshold; capillary-bundle and critical-path models are its diagonal and spectral limits. Hysteresis is the holonomy of a forcing bundle, a geometric phase rather than per-pore bistability, with a falsifiable loop-area law H ~ I^2. Preferential flow is what the same equation does where Da>1, so the Richards/preferential-flow dichotomy becomes a continuous Da-controlled crossover. Out of the quasi-static limit g(r) is the irreducible state variable. All inputs are geometric properties of the pore network, measurable from micro-CT and calibrated against no macroscopic data.
Le nombre de Damköhler, qui contrôle le comportement de l'eau dans le sol, vient à l'origine du génie chimique et décrivait la vitesse d'une réaction chimique par rapport au flux de matière. Aujourd'hui, il « travaille » aussi sous nos pieds.