The study, published in the journal Earth Surface Dynamics by Céleste Romon, Eric Lajeunesse, and François Métivier, describes a purpose-built experimental aquifer housed in a rectangular tank measuring 150 by 150 centimeters with a height of 35 centimeters. The researchers filled the upper compartment with a layer of plastic grains between 500 and 1000 micrometers in diameter, forming an erodible aquifer 15 to 20 centimeters thick. Beneath the sand, separated by a thin felt layer held between two metal grids, sits a water reservoir. Rather than sprinkling artificial rain onto the surface, an approach that muddies imaging, splashes grains out of place, and can generate surface runoff, the team injected water from below at a precisely controlled rate. The felt layer, whose hydraulic conductivity is roughly two orders of magnitude lower than that of the sand, acts as a resistor that distributes the incoming water uniformly across the base of the aquifer, mimicking the steady recharge that infiltration provides in nature.

Water leaves the tank through a narrow rectangular outlet cut into one wall, 15 centimeters above the aquifer bottom. As the water table rises to meet this opening, groundwater converges on the outlet and flows out, and if the discharge is strong enough, it drags plastic grains with it. The outflowing mixture, collected and weighed in a beaker at regular intervals, allowed the team to track both water discharge and the much smaller sediment flux, which amounted to only about 20 grams per hour. Within hours of the first run, seepage erosion had carved a channel extending backward from the outlet, its head taking on the amphitheater-shaped form that characterizes groundwater-sapped valleys on Earth and even on Mars. Over longer times, single channels split into two, branches competed for drainage area, and neighboring channels widened and merged, producing genuinely branching drainage networks several channels wide.

The most striking result concerns what happens when growth stops. In every preliminary run, erosion slowed steadily and eventually ceased entirely, with no further activity even after twelve hours of waiting. The only way to restart the network was to increase the recharge of the aquifer, which immediately triggered fresh channel growth until the system once again ground to a halt. The researchers interpret this behavior as evidence that, for a fixed recharge rate, the drainage network relaxes toward a steady-state morphology in which the sediments everywhere sit just at the threshold of motion. Too little groundwater flux, and nothing moves; the network has grown just large enough that the available water is spread thin across its entire perimeter. In this picture, the recharge rate effectively selects the size of the network, a hypothesis the team put to a demanding test.

In their central experiment, which ran continuously for 25 days, the researchers stepped the recharge upward ten times, beginning at roughly 0.1 liters per minute and increasing in increments of about 0.1 liters per minute. After each increase, they waited for the network to stabilize, verified by comparing photographs until no changes could be detected, and then measured the area enclosed by the channel network, tracing its contours manually with a measurement accuracy better than 4 percent. The result was remarkably clean: the steady-state network area increased linearly with the recharge rate, with a fitted proportionality constant of about 3.8 times 10 to the 3 seconds per meter. Because this relationship emerged from a single experiment, the authors caution that definitive conclusions require repetition, and they speculate that the slope of the relationship may depend on aquifer properties such as hydraulic conductivity and grain size. Still, the linear law suggests a simple and testable prediction for natural landscapes: where infiltration dominates over overland flow, the size of a drainage network may encode the intensity of local groundwater recharge.

Throughout the experiment, 22 piezometers, thin plastic tubes threaded along the aquifer bottom with narrowed tips that admit water but exclude grains, monitored the pressure field. Because the aquifer is shallow compared with its lateral extent, the groundwater flow obeys the Dupuit-Boussinesq approximation, a classical result dating to the 1860s and 1870s stating that vertical flow can be neglected to leading order. Under this approximation the pressure is hydrostatic, so the water level in each piezometer directly records the height of the water table above the aquifer bottom. Each step up in recharge produced an almost immediate rise in every piezometer, followed by a gradual relaxation as the network adjusted toward its new equilibrium. In steady state, the water table sloped downward toward the drainage network, reaching its minimum at the channels, exactly as field hydrogeologists observe in real catchments.

Armed with these boundary conditions, the team reconstructed the entire water table numerically. They solved the Dupuit-Boussinesq equation, which in time-averaged form becomes a Poisson equation for the square of the water-table height, using the finite-element method on a mesh refined where the water-table gradient is steep. The impervious tank walls supplied one boundary condition, and the drainage network supplied the other: the water table intersects the channels at the elevation of the streams. Since the experimental channels are only a few millimeters deep with slopes below 3 percent, the researchers approximated the whole network as lying at the elevation of the outlet. An iterative optimization adjusted the hydraulic conductivity to minimize the mismatch between the numerical solution and the 22 piezometric measurements, converging on a value of about 4.1 times 10 to the minus 3 meters per second, slightly higher than the value measured independently in a Darcy column, as expected for the looser packing of the experimental sand bed.

The reconstruction matched the piezometric data closely across nearly the entire aquifer, deviating only inside or immediately adjacent to the network, where the simplified boundary condition ignores the roughly one centimeter of actual water depth in the channels. A one-dimensional analysis in the paper’s appendix quantifies the consequence: neglecting the channel water depth overestimates the groundwater flux near the channel tips by a factor of about two. Away from the network, however, the method proved robust, and it revealed the mechanism at the heart of seepage-driven landscape evolution. Computed streamlines converged on the drainage network and concentrated near the channel tips, where the groundwater flux, given by the product of hydraulic conductivity, water-table height, and the gradient of the water table, reached values far higher than anywhere else in the aquifer. Close to the network, the iso-heads, the lines of constant water-table elevation, bent to follow the network’s outline, while at larger distances they smoothed out, indicating that the network’s hydrological influence extends over a region roughly proportional to its own size.

This concentration of flux at the tips closes the feedback loop that builds branching river networks. A larger groundwater discharge at a channel head entrains more sediment and drives faster erosion, while the flux along the channel sides remains too weak to mobilize grains. Erosion therefore undermines the land at the tip, the bank collapses, and the channel head advances into the aquifer, which in turn redirects still more groundwater toward the advancing front. The findings, consistent with earlier field reconstructions in the Florida Panhandle, provide the clearest laboratory demonstration yet that seepage erosion alone, without any surface runoff, can initiate and grow a branching drainage network, and that the network’s ultimate size is set by the recharge it must drain. The authors note that if the same logic applies in nature, many groundwater-dominated networks may operate near steady state today, though caution is warranted: landscapes preserve memory. During a field campaign in the Sanwara catchment in central India in the summer of 2024, the upper network carried no water despite a heavy monsoon, suggesting it was carved under a wetter past regime. The team is now extending the reconstruction method to field settings, aiming to estimate groundwater flow, storage, and river discharge from topographic maps alone, and testing fringe projection photogrammetry to capture channel topography and sharpen the flux estimates at the tips where rivers are born.

Subject of Research: Seepage erosion and the growth of groundwater-driven drainage networks in laboratory aquifer experiments

Article Title: Evolution of seepage driven networks in the lab

Article References: Romon, C., Lajeunesse, E., & Métivier, F. (2026). Evolution of seepage driven networks in the lab. Earth Surface Dynamics, 14(4), 517-525. https://doi.org/10.5194/esurf-14-517-2026

Image Credits: AI Generated

DOI: 10.5194/esurf-14-517-2026

Keywords: seepage erosion, drainage networks, groundwater flow, Dupuit-Boussinesq, laboratory experiments, aquifer recharge, channel heads, water table, geomorphology, river networks, steady state, hydraulic conductivity

The research, published in Ecology and Evolution, set out to determine whether aquatic invertebrate communities across the lower Colorado River Basin are shaped mainly by their position along the river network, by overland dispersal across the surrounding desert and mountains, or by the local environmental conditions of each habitat. To answer this, the team adapted a set of conceptual connectivity models originally developed in population genetics, including the stream hierarchy model, which predicts that communities are structured by hydrological connectivity, and the headwater model, which predicts that terrestrial connections among isolated headwater habitats dominate. These models had been proposed for metacommunity analysis but had rarely been tested against community-level data at such a large spatial scale.

The scale of the survey itself is a technical achievement. Between March and May 2022, researchers collected four one-liter water samples at each of 70 sites across five geographic regions: the Colorado River mainstem, Grand Canyon tributaries, the Verde River, Mogollon Rim headwater streams, and the Salt and Gila Rivers. Environmental DNA was extracted by filtering water onto 0.45-micrometer filters, and a 142-base-pair fragment of the mitochondrial cytochrome oxidase I gene was amplified with the fwhF2/EPTDr2n primer set and sequenced on an Illumina NextSeq 2000 platform. After filtering, reads were clustered into operational taxonomic units with less than 3 percent sequence divergence, yielding 1,445 aquatic-associated invertebrate OTUs, of which 72 percent belonged to the insect order Diptera.

The resulting community dataset was compared against four measures of distance between sites: straight-line Euclidean distance, topographic least-cost distance derived from a slope-weighted digital elevation model, river network distance along the dendritic channel system, and a composite environmental distance built from five significant variables selected by forward selection: drainage area, annual precipitation, elevation, stream type, and 100-year flood magnitude. Generalized dissimilarity models, which accommodate nonlinear relationships, revealed that environmental distance was the strongest single predictor of community dissimilarity, explaining 50.1 percent of the variance. Topographic least-cost distance followed at 42.6 percent, Euclidean distance at 40.4 percent, and river network distance at 38.8 percent, all highly significant.

When all three major distance types were entered into a full model, the picture sharpened further. The combined model explained 59.7 percent of the variance in community dissimilarity, with environmental distance and topographic distance both contributing significantly, but river network distance falling to statistical insignificance with a negligible effect. In other words, once the shape of the land and the character of the environment were accounted for, the length of the watery path between two sites added little explanatory power. For a branching river network long treated as the defining scaffold of aquatic ecology, this is a striking demotion of hydrological distance.

Perhaps the most vivid evidence came from the shape of the distance-decay curves themselves. Community dissimilarity rose with topographic and environmental distance, as expected, but along the river network axis the relationship was hook-shaped: dissimilarity increased over short to moderate network distances and then declined at the greatest separations. That decline matches the theoretical signature of the headwater model, which arises because hydrologically distant headwater streams converge in geographic space. Streams at opposite tips of the network may be separated by enormous lengths of river, yet they perch in similar high-elevation landscapes, and their invertebrate communities reflect that similarity rather than their separation by water.

Hierarchical clustering of the communities reinforced the same story. A two-cluster solution cleanly separated the dam-regulated Colorado River mainstem, including the entire Grand Canyon reach, from all remaining sites, likely reflecting the intense flow modifications that homogenize conditions and suppress invertebrate diversity below major dams. A five-cluster solution, validated by the Dunn index, subdivided the rest into groups that mostly followed basin boundaries, with one conspicuous exception: six Grand Canyon tributaries, including Bright Angel, Clear, Deer, Royal Arch, Shinumo, and Tapeats Creeks, clustered not with their hydrologic neighbors but with streams draining the Mogollon Rim, a highland region beyond a major drainage divide. Composition, in this case, ignored the map of the water.

Variance partitioning using distance-based redundancy analysis with hierarchical partitioning quantified the interplay. The overall model explained 47 percent of community variation, and the single largest fraction was the variance shared jointly by environment, terrestrial geography, and river network connectivity, indicating that these landscape characteristics covary in ecologically meaningful ways. The next largest fraction was uniquely attributable to the environment, followed by shared variance between geography and the network. All three predictors remained individually significant, confirming that local filtering, terrestrial dispersal pathways, and network structure each leave distinct fingerprints, even though more than 70 percent of the explained variance involved the environment in some form, a strong signal of species sorting.

The authors are careful to note what the analysis cannot resolve. Roughly half of the community variation remained unexplained, a common outcome in ecology but one that may point to processes outside the classical metacommunity archetypes, such as floods, droughts, wildfire, and the pervasive imprint of dams and diversions that fragment the Colorado River system. Desert rivers are disproportionately altered by impoundments, and earlier work has shown that hydropeaking operations can extirpate aquatic insects and undermine river food webs. Disentangling these disturbance effects from the connectivity signals documented here remains an open challenge.

Still, the implications reach well beyond the American Southwest. If environmental similarity and terrestrial landscape structure can override river network distance in structuring aquatic communities, then conservation planning that treats rivers as linear corridors alone may miss the true pathways of biodiversity. Headwater habitats, often dismissed as remote and disconnected, emerge as nodes linked across basin boundaries by the topography they share. Managing for overland connectivity, protecting environmentally diverse headwater refugia, and recognizing that two streams far apart by water can be close in life may prove essential for safeguarding freshwater biodiversity in an era when drought, diversion, and dam operations continue to reshape the world’s great river networks.

Subject of Research: Ecological connectivity and metacommunity structure of aquatic invertebrate communities in the lower Colorado River Basin assessed with environmental DNA metabarcoding

Article Title: The Further You Are, the Closer You Get: Environmental Similarity Connects Aquatic Invertebrate Communities Across the Terrestrial Landscape in a Large River Network

Article References: The Further You Are, the Closer You Get: Environmental Similarity Connects Aquatic Invertebrate Communities Across the Terrestrial Landscape in a Large River Network. (n.d.). https://doi.org/10.1002/ece3.74302

Image Credits: AI Generated

DOI: 10.1002/ece3.74302

Keywords: environmental DNA, aquatic invertebrates, river networks, metacommunity ecology, distance decay, ecological connectivity, species sorting, headwater streams, Colorado River, biogeography, generalized dissimilarity modeling, desert rivers