Geology Reports⌕ Search

Geology topics

John R. Nimmo

Publications and source records attributed to John R. Nimmo.

At least 19 recordsLinked to original sources
USGS2025

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS2025

Identifying preferential flow from soil moisture time series: Review of methodologies

Identifying and quantifying preferential flow (PF) through soil—the rapid movement of water through spatially-distinct pathways in the subsurface—is vital to understanding how the hydrologic cycle responds to climate, land cover, and anthropogenic changes. In recent decades, methods have been developed that use measured soil moisture time series to identify PF. Because they allow for continuous monitoring and are relatively easy to implement, these methods have become an important tool for recognizing when, where, and under what conditions PF occurs. The methods seek to identify a pattern or quantification that indicates the occurrence of PF. Most commonly, the chosen signature is either (1) a nonsequential response to infiltrated water, in which soil moisture responses do not occur in order of shallowest to deepest, or (2) a velocity criterion, in which newly infiltrated water is detected at depth earlier than is possible by nonpreferential flow processes. Alternative signatures have also been developed that have certain advantages but are less commonly utilized. Choosing among these possible signatures requires attention to their pertinent characteristics, including susceptibility to errors, possible bias toward false negatives or false positives, reliance on subjective judgments, and possible requirements for additional types of data. We review 77 studies that have applied such methods, to highlight important information for readers who want to identify PF from soil moisture data, and to inform those who aim to develop new methods or improve existing ones.

Vadose Zone Journal↗
USGS2024

New insights on the origin of the Richardson-Richards equation

The Richardson-Richards equation (RRE), despite known shortcomings especially in regard to preferential flow, provides the basis of the vast majority of unsaturated flow models in use today. L.F. Richardson published this equation in 1922, nine years before L.A. Richards. Whereas Richards approached this problem directly from the groundbreaking developments of Edgar Buckingham, Richardson, surprisingly, cited as his starting point only the earlier work of L.J. Briggs. Collectively, these four scientists’ published and unpublished work reveals that: (1) Briggs’ work, though qualitative, captured the essential physical principles needed for quantifying unsaturated flow; (2) Buckingham came very close to deriving the RRE and explained why he stopped short of doing so; (3) derivation of the RRE from the work of either Briggs or Buckingham required only modest developmental work; and (4) besides deriving the RRE, Richards carried through much of the experimental agenda that Buckingham considered a necessary precursor to mathematical treatment.

Hydrological Sciences Journal↗
USGS2023

Detection and quantification of preferential flow using artificial rainfall with multiple experimental approaches

Preferential flow in the unsaturated zone strongly influences important hydrologic processes, such as infiltration, contaminant transport, and aquifer recharge. Because it entails various combinations of physical processes arising from the interactions of water, air, and solid particles in a porous medium, preferential flow is highly complex. Major research is needed to improve the ability to understand, quantify, model, and predict preferential flow. Toward a solution, a combination of diverse experimental measurements at multiple scales, from laboratory scale to mesoscale, has been implemented to detect and quantify preferential paths in carbonate and karstic unsaturated zones. This involves integration of information from (1) core samples, by means of mercury intrusion porosimeter, evaporation, quasi-steady centrifuge and dewpoint potentiometer laboratory methods, to investigate the effect of pore-size distribution on hydraulic characteristics and the potential activation of preferential flow, (2) field plot experiments with artificial sprinkling, to visualize preferential pathways related to secondary porosity, through use of geophysical measurements, and (3) mesoscale evaluation of field data through episodic master recession modeling of episodic recharge. This study demonstrates that preferential flow processes operate from core scale to two different field scales and impact on the qualitative and quantitative groundwater status, by entailing fast flow with subsequent effects on recharge rate and contaminant mobilizing. The presented results represent a rare example of preferential flow detection and numerical modeling by reducing underestimation of the recharge and contamination risks.

Bari↗
USGS2023

Improved calculation of hydraulic conductivity for small-disk tension infiltrometers

Because tension infiltrometers apply water through a disk of finite size, the infiltrated water moves laterally as well as downward. Only the vertical component of this flow is indicative of the hydraulic conductivity K , so the algorithm for computing K must include a way of isolating that component from the total flow. Some commonly used formulas correct for the multidimensional effects by subtracting an estimate of the laterally spreading flow. For disks smaller than about 200 mm in diameter, however, lateral spreading constitutes so much of the total flow that these subtractive formulas lose considerable accuracy, and sometimes overcorrect so severely as to produce a negative number for K . Other methods rely on empiricisms that are not completely consistent with unsaturated-flow theory and that require prior knowledge of certain soil properties. We developed a new formula that uses a multiplicative factor instead of a subtracted term to achieve the needed correction. For testing we conducted numerical experiments with synthetic data produced by solving the Richardson-Richards equation using the code VS2DRTI, for diverse media and a range of disk sizes, including the widely used 45-mm diameter. We compared K values calculated from our formula to the actual K used to generate the simulated data, as well as to results from other published formulas. This comparison shows that our method provides an algorithm based in unsaturated-flow theory that produces more reliable values for small disks without requiring prior knowledge of soil properties.

Water Resources Research↗
USGS2023

Porosity and pore-size distribution

Porosity, the fraction of soil volume not occupied by solids, is relatively easy to conceptualize and measure. Pore-size distribution is a complex topic, in part from the lack of a clear and unique concept of a soil pore as a discrete object. Available tools for evaluating pore-size distribution involve traditional conventions and operational definitions applied to hydraulic property characterizations, as well as more recently developed imaging techniques. Application of these tools can provide effective characterizations of the pores and pore space that lead to an understanding of the relation between pore characteristics and soil functions.

Book chapter↗
USGS2021

Rapid-response unsaturated zone hydrology: Small-scale data, small-scale theory, big problems

The unsaturated zone (UZ) extends across the Earth’s terrestrial surface and is central to many problems related to land and water resource management. Flow of water through the UZ is typically thought to be slow and diffusive, such that it could attenuate fluxes and dampen variability between atmospheric inputs and underlying aquifer systems. This would reduce water resource vulnerability to contaminants and water-related hazards. Reducing or negating that effect, however, spatially concentrated and rapid flow and transport through the unsaturated zone is surprisingly common and becoming more so with the increasing frequency and magnitude of extreme hydroclimatic events. Arising from the wide range in the rates and complex modes of nonlinear flow processes, these effects are among the most poorly characterized hydrologic phenomena. Issues of scale present additional difficulties. Equations representing unsaturated processes have been developed and tested on the basis of field and laboratory measurements typically made at scales from pore size to plot size. In contrast, related problems of significant interest to society, including floods, aquifer recharge, landslides, and groundwater contamination, range from watershed to regional scales. The disparity between the scale of our understanding and the scale of interest for societal problems has spurred application of these model equations at increasingly coarse resolutions over larger areas than can be justified by existing measurements or theory. This mismatch in scales requires an assumption that spatially averaging slow diffusive flow and rapid preferential flow can effectively represent the influence of both processes across vast areas. Given the currently inadequate recognition and quantitative characterization of focused and rapid processes in unsaturated flow, these phenomena are critically in need of expanded attention and effort.

Frontiers in Earth Science↗
USGS2021

The processes of preferential flow in the unsaturated zone

Preferential flow, a major influence in unsaturated soil and rock almost everywhere, occurs by multiple phenomenologically distinct hydraulic processes. For the mode known as funneled flow, concentrated in particularly conductive portions of the medium, the surface-tension/viscous-flow processes of traditional unsaturated flow theory predominate. Fingered flow, through conductive paths of higher water content than surrounding material, requires amendments to traditional theory concerning instabilities and dynamic flow-regime boundaries. Macropore flow, the most recognized preferential flow mode, poses unanswered questions and major difficulties in practice. Accumulated evidence shows that water flows preferentially mostly through macropores that are (a) only partially filled with water, and (b) surrounded by matrix material that is drier, sometimes much drier, than saturation. With partial filling, geometric characteristics such as aperture have much less influence than was previously thought, and the intra-macropore configuration of the flowing water phase, about which little is conclusively known, is then a dominant controlling influence. With unsaturated surroundings, macropore/matrix exchange interactions control, for given input and medium, the initiating circumstances, conveyed flux, and duration of macropore flow. The multiple processes in play during such interactions have different sensitivities to the matrix water state and different directions of influence. The net influence of matrix water content on macropore flow is thus highly complex and a major research need. Additional high-priority topics are: flowpath connectivity, for watersheds as well as small scales; intra-macropore processes, to discern their importance and possible means of quantification; and the identification of measurable soil and rock properties that can be utilized predictively.

Soil Science Society of America Journal↗
USGS2020

Does the Darcy-Buckingham Law apply to flow through unsaturated porous rock?

The Darcy–Buckingham (DB) law, critical to the prediction of unsaturated flow, is widely used but has rarely been experimentally tested, and therefore may not be adequate in certain conditions. Failure of this law would imply that the unsaturated hydraulic conductivity is not constant for a given water content, as assumed in nearly all subsurface flow models. This study aims to test the DB law on unsaturated porous rock, complementing the few previous tests, all done on soils. Two lithotypes of calcareous porous rocks were tested. The quasi-steady centrifuge method was used to measure the flux density for different centrifugal driving forces while maintaining essentially constant water content, as required. Any deviations from the direct proportionality of the measured flux and the applied force would indicate a violation of the DB law. Our results show that, for the tested rocks and conditions, no physical phenomena occurred to cause a failure of the DB law.

Water↗
USGS2018

Episodic master recession evaluation of groundwater and streamflow hydrographs for water-resource estimation

Hydrograph analysis tools using a master recession curve (MRC) can produce many types of hydrologically important watershed-response quantifications, including aquifer recharge and stormflow characterization. An MRC is the relation between the value of a measured response R and its rate of change with time, d R /d t , occurring on the falling limb when there is no infiltration or other water input. We have developed MRC and episodic hydrograph-evaluation methods for multiple purposes, utilizing both water table and streamflow data. The determination of a parameterized MRC through a structured procedure provides a basis for quantification of hydrologic variables and characteristics that can be validly compared among different events, sites, and periods of time. Application of the MRC to needed hydrologic quantifications is done with our revised episodic master recession (EMR) method. Expert-guided iterative procedures are used to quantify parameters needed in applying the MRC and EMR methods to a given site. Hydrologic judgments such as the significance threshold for response magnitude, and the time window within which the precipitation is assumed to be the cause of an observed response, inherently involve some elements of subjectivity. Our structured iterative approach, however, affords much flexibility in formulating expert judgments and serves to confine them to statements and procedures that can be quantified and documented. Parallel application to streamflow and water table hydrographs can produce new hydrologic insights and understanding, not least in the role of unsaturated zone processes in controlling exchanges among components of the water cycle.

Vadose Zone Journal↗
USGS2018

UZIG research: Measurement and characterization of unsaturated zone processes under wide-ranging climates and changing conditions

Unsaturated zone properties and processes are central to understanding the interacting effects of land-use change, contamination, and hydroclimate on our ability to grow food, sustain clean water supplies, and minimize loss of life and property. Advances in unsaturated zone science are being achieved through collaborations across traditional boundaries where information from biological, physical, and chemical disciplines is combined for new insights. The Unsaturated Zone Interest Group (UZIG) is an organization that exists principally to promote multidisciplinary collaborations and the sharing of ideas, expertise, and technical assets. Here we summarize key findings from 14 papers, several of which originated from a meeting convened by UZIG in 2017 at the University of Florida in Gainesville titled “Land-Use Change, Climate Change, and Hydrologic Extremes: Unsaturated Zone Responses and Feedbacks.” This special section of Vadose Zone Journal contains multidisciplinary research in three general categories relevant to measuring and understanding unsaturated zone responses to changing land uses and climate: (i) unsaturated zone properties and processes; (ii) soil–plant–atmosphere interactions; and (iii) novel field sampling devices. A strong cross-cutting theme in these papers is the value of continuous monitoring data and ways of utilizing them to discover novel hydrologic, biologic, and pedologic information. As climatic and land-use conditions change and demands for resources and stresses on ecosystems continue to intensify, it is vital to improve our fundamental understanding of the processes at work in the unsaturated zone. Toward that goal, we discuss the need for improved ground-based unsaturated zone monitoring networks.

Vadose Zone Journal↗
USGS2018

Vegetation influences on infiltration in Hawaiian soils

Changes in vegetation communities caused by removing trees, introducing grazing ungulates, and replacing native plants with invasive species have substantially altered soil infiltration processes and rates in Hawaii. These changes directly impact run-off, erosion, plant-available water, and aquifer recharge. We hypothesize that broad vegetation communities can be characterized by distributions of field-saturated hydraulic conductivity (Kfs). We used 290 measurements of Kfs calculated from infiltration tests from 5 of the Hawaiian Islands to show this effect. We classified the data using 3 broad ecosystem categories: grasses, trees and shrubs, and bare soil. The soils of each site have coevolved with past and present ecological communities without significant mechanical disturbance by agriculture or urban development. Geometric mean values Kfs are 203 mm/hr for soils hosting trees and shrubs, 50 mm/hr for grasses, and 13 mm/hr for bare soil. Differences are statistically significant at the 95% confidence level. These examples show that it is feasible to make maps of relative Kfs based on field and ecosystem data. These ecosystem trends can be used to estimate ongoing changes to run-off and recharge from climate and land use change. Greater Kfs for ecosystems with primarily trees and shrubs suggests that management decisions concerning reforestation or other changes of vegetation can have substantial hydrologic impacts.

Hawaii↗
USGS2017

Evaluation of diffuse and preferential flow pathways of infiltratedprecipitation and irrigation using oxygen and hydrogen isotopes

Subsurface-water flow pathways in three different land-use areas (non-irrigated grassland, poplar forest, and irrigated arable land) in the central North China Plain were investigated using oxygen ( 18 O) and hydrogen ( 2 H) isotopes in samples of precipitation, soils, and groundwater. Soil water in the top 10 cm was significantly affected by both evaporation and infiltration. Water at 10–40 cm depth in the grassland and arable land, and 10–60 cm in poplar forest, showed a relatively short residence time, as a substantial proportion of antecedent soil water was mixed with a 92-mm storm infiltration event, whereas below those depths (down to 150 cm), depleted δ 18 O spikes suggested that some storm water bypassed the shallow soil layers. Significant differences, in soil-water content and δ 18 O values, within a small area, suggested that the proportion of immobile soil water and water flowing in subsurface pathways varies depending on local vegetation cover, soil characteristics and irrigation applications. Soil-water δ 18 O values revealed that preferential flow and diffuse flow coexist. Preferential flow was active within the root zone, independent of antecedent soil-water content, in both poplar forest and arable land, whereas diffuse flow was observed in grassland. The depleted δ 18 O spikes at 20–50 cm depth in the arable land suggested the infiltration of irrigation water during the dry season. Temporal isotopic variations in precipitation were subdued in the shallow groundwater, suggesting more complete mixing of different input waters in the unsaturated zone before reaching the shallow groundwater.

Yellow Sea, Bohai Sea↗
USGS2016

Preferential flow, diffuse flow, and perching in an interbedded fractured-rock unsaturated zone

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.

Idaho↗
USGS2016

Quantitative framework for preferential flow initiation and partitioning

A model for preferential flow in macropores is based on the short-range spatial distribution of soil matrix infiltrability. It uses elementary areas at two different scales. One is the traditional representative elementary area (REA), which includes a sufficient heterogeneity to typify larger areas, as for measuring field-scale infiltrability. The other, called an elementary matrix area (EMA), is smaller, but large enough to represent the local infiltrability of soil matrix material, between macropores. When water is applied to the land surface, each EMA absorbs water up to the rate of its matrix infiltrability. Excess water flows into a macropore, becoming preferential flow. The land surface then can be represented by a mesoscale (EMA-scale) distribution of matrix infiltrabilities. Total preferential flow at a given depth is the sum of contributions from all EMAs. Applying the model, one case study with multi-year field measurements of both preferential and diffuse fluxes at a specific depth was used to obtain parameter values by inverse calculation. The results quantify the preferential–diffuse partition of flow from individual storms that differed in rainfall amount, intensity, antecedent soil water, and other factors. Another case study provided measured values of matrix infiltrability to estimate parameter values for comparison and illustrative predictions. These examples give a self-consistent picture from the combination of parameter values, directions of sensitivities, and magnitudes of differences caused by different variables. One major practical use of this model is to calculate the dependence of preferential flow on climate-related factors, such as varying soil wetness and rainfall intensity.

Vadose Zone Journal↗