How advanced modeling reveals the invisible journeys of contaminants through underground environments
Beneath our feet, an invisible world of constant movement and transformation determines the safety of the water we drink and the health of our ecosystems. Every day, countless contaminantsâfrom industrial solvents and pharmaceuticals to agricultural chemicals and pathogensâbegin hidden journeys through the subsurface, traveling through layers of soil and rock, sometimes transforming along the way, eventually reaching aquifers, rivers, and eventually, our taps.
Understanding these complex journeys is the realm of subsurface fate and transport modeling, an interdisciplinary science that combines hydrology, chemistry, and advanced computation to predict how pollutants move through underground environments. This field stands as our first line of defense in protecting drinking water supplies, guiding environmental cleanup, and preventing public health crises. From the devastating discovery of contaminated wellfields to the silent threat of vapor intrusion into homes, subsurface modeling provides the predictive power needed to make informed decisions about our most precious resource: clean water.
In environmental science, "fate and transport" refers to the combined study of how contaminants change and move through the environment. Fate describes what eventually happens to a contaminantâwhether it persists, breaks down, or transforms into something else. Transport covers the physical movement of contaminants within and between environmental compartments like soil, groundwater, and air 4 . These processes are governed by fundamental principles of chemistry, physics, and biology that determine whether a contaminant plume will spread, stagnate, or dissipate over time.
The science of subsurface modeling has evolved dramatically from simple analog models to sophisticated computational tools. Early approaches relied heavily on analytical solutionsâmathematical equations that could predict contaminant movement under idealized, simplified conditions. These remain valuable for screening-level assessments but struggle to capture the complex reality of heterogeneous subsurface environments.
The advent of powerful computing brought numerical models that can simulate contaminant behavior in complex, real-world conditions by breaking down the subsurface into thousands or millions of discrete elements and solving the governing equations for each one.
Contaminants carried along with flowing groundwater
Spreading and mixing due to velocity variations
Attachment of contaminants to soil particles
Chemical and biological breakdown processes
Movement from water/soil into air
Different transport behaviors for microbes vs chemicals
The theoretical foundations of subsurface modeling rest on solving the advection-dispersion-reaction equation, a mathematical representation that describes how contaminants move through porous media while undergoing various transformations. Early models provided solutions for this equation under simplified conditions, but real-world applications demanded more sophisticated approaches.
Recent research has focused on integrating multiple processes that were previously studied in isolation. For instance, scientists have developed novel frameworks that simultaneously account for advection, rate-limited adsorption on suspended sediments, and first-order degradation in river systems 5 . These integrated approaches are particularly valuable for understanding short-term solute accumulation in riverbedsâa critical yet previously understudied process that significantly affects pollutant fate.
One of the most important theoretical realizations in recent decades is that microbes and dissolved chemicals travel differently through subsurface environments. While this might seem intuitiveâgiven the vast difference in size between a bacterium and a solvent moleculeâthe implications are profound for protecting drinking water supplies.
Unlike dissolved chemicals, microbes behave as colloids (particles suspended in water) and experience size exclusion from smaller pores. This means their transport is limited to larger, typically more conductive pore spaces, resulting in faster breakthrough compared to solutes 2 . This discovery overturned previous assumptions that solute tracers were conservative proxies for microbial transport and explained why groundwater pumped from riverbank filtration systems can contain high microbial loads following heavy rain events, even when chemical indicators suggest safe conditions.
| Characteristic | Solute Transport | Microbial Transport |
|---|---|---|
| Primary Mechanism | Advection-dispersion | Advection with size exclusion |
| Flow Path | All available pore spaces | Only larger, more conductive pores |
| Breakthrough Timing | Later peak concentration | Earlier first arrival |
| Key Retention Processes | Sorption, degradation | Attachment, detachment, straining, inactivation |
| Health Concern | Chronic exposure | Acute exposure |
Riverbank filtrationâthe natural purification of river water as it infiltrates through adjacent sediments to pumping wellsâhas been used for centuries to produce drinking water. However, the detection of pathogens in supposedly protected wellfields prompted researchers to investigate why the conventional models based on solute transport were failing to protect public health.
A pivotal research effort led to the implementation of a dual-permeability, two-site kinetic deposition formulation for microbial transport in the integrated surface-subsurface hydrological model HydroGeoSphere (HGS) 2 . This groundbreaking work explicitly simulated microbial transport in river-groundwater systems, accounting for attachment, detachment, and inactivation of microbes in different permeability regions, while allowing for multispecies transport.
Riverbank filtration systems provide natural purification but require advanced modeling to ensure protection against microbial contamination.
Researchers implemented a dual-permeability formulation with two kinetic deposition sites in HGS, considering both fast, reversible retention in large-pore regions and slow, irreversible retention in small-pore regions.
The implementation was first verified against an established analytical solution for dual-permeability colloid transport to ensure mathematical correctness.
The model's suitability was demonstrated through two illustrative scenarios: an integrated rainfall-runoff and streamflow generation benchmark with added microbial transport from a conceptual manure application, and an idealized alluvial riverbank filtration site simulating parallel transport of reactive microbes, conservative helium-4, and reactive radon-222.
The model enabled direct comparison of breakthrough curves, travel times, and mixing ratios between different tracers, highlighting the divergent behavior of microbes versus solutes.
The simulations revealed several crucial findings with significant implications for drinking water protection:
| Contaminant Type | Example | Travel Time (days) |
|---|---|---|
| Conservative Solute | Chloride, Helium-4 | 45 |
| Reactive Solute | Radon-222 | 52 |
| Viruses | Bacteriophage MS2 | 28 |
| Bacteria | E. coli | 31 |
The advancement of subsurface fate and transport modeling relies on a sophisticated toolkit of computational methods and analytical approaches. These tools range from conceptual models that help frame our understanding to numerical models that provide quantitative predictions of contaminant behavior.
The U.S. Environmental Protection Agency maintains a portfolio of groundwater modeling tools, including MT3D for simulating advection, dispersion, and chemical reactions of dissolved constituents; BIOPLUME III for modeling natural attenuation of organic contaminants; and WhAEM2000 for delineating wellhead protection areas 3 . These publicly available models represent valuable resources for researchers and practitioners alike.
More recent advances include integrated simulators like HydroGeoSphere (HGS) with its microbial transport capabilities 2 , and analytical frameworks implemented in libraries like SWASHES that provide reliable benchmarks for validating numerical models 5 .
| Tool/Model | Primary Application | Key Features |
|---|---|---|
| HydroGeoSphere (HGS) | Integrated surface-subsurface microbial and solute transport | Dual-permeability, two-site kinetic deposition; multispecies support |
| MT3D | 3D solute transport in groundwater | Simulation of advection, dispersion, and chemical reactions |
| BIOPLUME III | Natural attenuation of organics | Models biodegradation processes in groundwater |
| SWASHES Library | Solute transport in rivers | Analytical solutions for model validation |
| Quick_Domenico | Plume concentration calculations | Spreadsheet-based tool for screening assessments |
| 3-Benzylazetidin-3-ol | 1236862-03-5 | C10H13NO |
| Aphidicolin glycinate | C22H37NO5 | |
| N-Hexadecyl-L-alanine | 671247-18-0 | C19H39NO2 |
| N-Isobutylphthalimide | 304-19-8 | C12H13NO2 |
| Iodobenzene diacetate | C10H13IO4 |
A receptor modeling technique used to identify and quantify pollution sources
Artificial neural networks that visualize and cluster high-dimensional data
Computational technique using random sampling to account for uncertainty
Statistical method that reduces dimensionality of complex datasets
As we look to the future of subsurface fate and transport modeling, several critical challenges and research priorities emerge. Perhaps most pressing is the need to address emerging pollutantsâcontaminants not previously recognized or regulated, such as pharmaceuticals, personal care products, and antibiotic resistance genes 1 . The presence of antiretroviral drugs in water systems and their toxicity to cyanobacteria, as highlighted in recent research 1 , illustrates the expanding scope of contaminants requiring investigation.
Future research should adopt a holistic approach that considers the "prevention, generation, monitoring, quantification, interaction, and removal or remediation of pollutants throughout the whole system" 1 . This systems perspective acknowledges that contaminants move between environmental compartments and that effective management requires integrated understanding and solutions.
The future of subsurface modeling will likely be shaped by several innovative approaches:
To raise public awareness about pollutants generated at home and proper disposal methods 1
To accurately detect and quantify pollutants at increasingly lower concentrations 1
Including "effective phytoremediation technologies and biochar-enhanced filtration" 1
Subsurface fate and transport modeling represents one of our most powerful tools for ensuring water security in an era of increasing contamination threats. What was once a science focused primarily on simple dissolved plumes has evolved into a sophisticated discipline that recognizes the complex, interconnected nature of subsurface environments and the diverse behaviors of different contaminant types.
From the critical realization that microbes travel faster than solutes to the development of integrated models that simulate both surface and subsurface processes, the field has made remarkable advances. Yet challenges remain as new contaminants emerge and climate change alters hydrological systems. The continued development of modeling toolsâcoupled with enhanced monitoring techniques and interdisciplinary collaborationâwill be essential for protecting drinking water supplies and ecosystem health for generations to come.
As research continues to illuminate the hidden world beneath our feet, we gain not only scientific knowledge but also practical wisdom for managing our water resources more sustainably. The journey of contaminants through the subsurface may be invisible to our eyes, but through the power of modeling, we can trace their paths, predict their destinations, and intervene to protect what matters mostâthe safety of our water and the health of our communities.
References will be added here in the appropriate format.