Overview
Water movement through unsaturated soils is one of the most fundamental processes in environmental science, controlling everything from agricultural productivity to contaminant transport and groundwater recharge. Unlike saturated flow where all pores are filled with water, unsaturated flow involves complex interactions between water, air, and the soil matrix.
Our research combines advanced measurement techniques, theoretical analysis, and numerical modeling to understand and predict water flow under diverse conditions—from desert soils experiencing extreme dryness to agricultural fields under irrigation.
The Richardson-Richards Equation
At the heart of unsaturated flow modeling is the Richardson-Richards equation (RRE), a partial differential equation that describes water movement based on gradients in water potential. The RRE accounts for:
- Capillary forces: Water retention in small pores
- Gravity: Downward water movement
- Hydraulic conductivity: Soil's ability to transmit water, which varies dramatically with water content
The nonlinear nature of the RRE makes it challenging to solve, especially for heterogeneous soils with contrasting layers or under dynamic boundary conditions. Our work develops novel solution approaches and evaluates their accuracy across diverse soil conditions.
Key Research Areas
1. Hydraulic Property Characterization
Accurate modeling requires knowing two key soil functions:
- Water retention curve (WRC): Relationship between water content and matric potential
- Hydraulic conductivity function (HCF): How water conductivity changes with water content or potential
We use advanced laboratory methods including:
- Evaporation method (HYPROP) for simultaneous WRC and HCF measurement
- Tempe cells for water retention at multiple pressures
- Infiltration experiments for field-scale hydraulic properties
2. Scale Mismatch Problems
A persistent challenge is that laboratory-measured properties often don't represent field-scale behavior. We address this through:
- Inverse modeling: Estimating field-scale properties from in situ moisture measurements
- Physics-informed machine learning: Learning hydraulic functions directly from field data
- Upscaling methods: Mathematical approaches to translate small-scale measurements to larger scales
3. Infiltration and Redistribution
We investigate how water enters soil (infiltration) and subsequently redistributes after rainfall or irrigation ceases. Key questions include:
- How do soil structure and macropores affect infiltration rates?
- What controls the depth and rate of water penetration?
- How long does water remain available to plants versus draining to groundwater?
4. Evaporation from Bare Soils
Desert soils and agricultural fields between crops lose significant water to evaporation. Our research examines:
- Stage transitions in evaporation (energy-limited vs. diffusion-limited)
- Effects of soil texture and structure on evaporative losses
- Improved hydraulic functions for very dry conditions (beyond pF 3.8)
Improved Hydraulic Functions
Peters-Durner-Iden (PDI) Model
Traditional van Genuchten functions often fail for very dry soils. We've demonstrated that the PDI model significantly improves predictions of:
- Moisture redistribution in desert environments
- Evaporation rates from drying soils
- Water retention between pF 2 and pF 3.8 (field capacity to near wilting point)
Applications
Agricultural Water Management
Understanding infiltration and redistribution helps optimize irrigation scheduling, reducing water waste while maintaining crop productivity. Our models predict:
- How much water infiltrates versus runs off
- Water availability in the root zone over time
- Deep drainage losses below the root zone
Wastewater Irrigation
Treated wastewater is increasingly used for irrigation in water-scarce regions. We've shown that long-term wastewater application:
- Reduces soil hydraulic conductivity and infiltration capacity
- Increases runoff and erosion risk
- Requires adjusted irrigation management strategies
Groundwater Recharge
Predicting deep percolation is essential for groundwater management. Our models help quantify recharge under different:
- Climate scenarios (rainfall intensity and distribution)
- Land use practices (tillage, cover crops)
- Soil types and layering
Numerical Modeling Approaches
We employ multiple modeling platforms depending on the problem:
- HYDRUS-1D/2D: Established finite element solver for the Richards equation
- Physics-informed neural networks: Novel AI approach that learns solutions from sparse data
- Analytical solutions: Simplified cases for validation and insight
Critical to all modeling is validation. We rigorously compare simulations against:
- Analytical solutions where available
- Laboratory experiments with known boundary conditions
- Field measurements from instrumented plots
Recent Publications
- Impact of almond shell biochar properties and application rate on soil physical and hydraulic characteristics.Thao, T., Lopez, V. D., Gonzales, M., Berhe, A. A., Diaz, G., & Ghezzehei, T. A.Sustainable Environment, 11(1), 2485688. 2025.
Abstract
We conducted two 64-day incubation experiments to assess how locally produced almond-shell biochar influences soil physical and hydraulic properties. Biochar was created using slow pyrolysis at different temperatures (350 °C or 700 °C), separated into different particle sizes (<250 μm or 1–2 mm), and applied at 10 ton/ha or 60 ton/ha to a coarse-textured soil. While our analysis shows that biochar yielded greater cation exchange capacity (CEC) and specific surface area (SSA) with increasing pyrolysis temperature and finer particle size, its contributions to improving soil hydraulic properties were marginal. In the first experiment, the addition of biochar at high rates slightly improved water stable aggregate (WSA) (3.8%–5.3% increase) but has no effect on saturated hydraulic conductivity (Ksat). Soil respiration measured throughout the experiment were not significantly different among treatments. In the second experiment, the addition of biochar increased soil infiltration rate at the initial stage (8.18E–4 cm/s), but this effect diminished over time. WSA was lower for biochar amended soil and lowest at high application rates (5%–21% reduction). Cumulative carbon dioxide (CO2) flux varied between biochar particle sizes and rates. Additionally, a significant difference between the two experiments was also observed, with cumulative CO2 (38%–56% greater) and WSA (11%–40%) being inversely correlated. Our findings suggest that almond-shell derived biochar has a limited impact on arable loamy sand soil properties, specifically for water retention under short-term conditions.BibTeX
@article{Thao31122025, author = {Thao, Touyee and Lopez, Vivian D. and Gonzales, Melinda and Berhe, Asmeret A. and Diaz, Gerardo and Ghezzehei, Teamrat A.}, title = {Impact of almond shell biochar properties and application rate on soil physical and hydraulic characteristics}, journal = {Sustainable Environment}, volume = {11}, number = {1}, pages = {2485688}, year = {2025}, publisher = {Taylor \& Francis}, doi = {10.1080/27658511.2025.2485688}, pdf = {https://doi.org/10.1080/27658511.2025.2485688}, research-theme = {water-flow, soil-structure, sustainable-agriculture} } - Biochar Impacts on Soil Moisture Retention and Respiration in a Coarse-Textured Soil under Dry Conditions.Thao, T., Harrison, B., Gonzalez, M., Ryals, R., Dahlquist‐Willard, R., Diaz, G. C., & Ghezzehei, T. A.Soil Sci. Soc. Am. J., (Early View), 1–13. 2024.
Abstract
The growing water scarcity jeopardizes crop production for global food security, a problem poised to worsen under climate change–induced drought. Amending soils with locally derived biochar from pyrolyzed agricultural residues may enhance soil moisture retention and resilience, in addition to climate change mitigation. However, prior studies on the hydrologic benefits of biochar focused on optimal moisture, not water-limited conditions where biochar’s large wettable surface area could aid plants and microbes. We hypothesized that biochars differing in feedstocks would positively augment soil moisture and respiration, with overall impacts most beneficial under drier conditions. Using water vapor sorption isotherms, we used film theory to estimate the specific surface area (SSA) of biochars. We then modeled and tested the moisture retention of a coarse-textured soil amended with biochar. Additionally, a 109-day lab incubation experiment was also conducted to examine biochar effects on respiration across a moisture range spanning optimal to wilting point. Among seven tested biochars, almond shell biochar significantly increased soil moisture and yield the second highest SSA. Despite drying treatments, the amended soil maintained higher respiration than the control, indicating enhanced biological activity. The results demonstrate biochars counter drying effects in coarse soils through physical and biological mechanisms linked to increased sorptive capacity. Our findings contribute to the development of sustainable water and waste management strategies tailored to the needs of California Central Valley, where the potential for biochar application is substantial. Above all, our research fills a crucial gap by providing context-specific insights that can inform the effective utilization of locally produced biochars in the face of increasing water scarcity and excess biomass challenges.BibTeX
@article{p2024-Taho-et-al, title = {Biochar Impacts on Soil Moisture Retention and Respiration in a Coarse-Textured Soil under Dry Conditions}, number = {(Early View)}, year = {2024}, language = {en}, journal = {Soil Sci. Soc. Am. J.}, doi = {10.1002/saj2.20746}, pdf = {https://acsess.onlinelibrary.wiley.com/doi/epdf/10.1002/saj2.20746}, author = {Thao, Touyee and Harrison, Brendan and Gonzalez, Melinda and Ryals, Rebecca and Dahlquist‐Willard, Ruth and Diaz, Gerardo C. and Ghezzehei, Teamrat A.}, pages = {1-13}, research-theme = {water-flow, soil-structure, rhizosphere, sustainable-agriculture} } - Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks.Bandai, T., Ghezzehei, T. A., Jiang, P., Kidger, P., Chen, X., & Steefel, C. I.Water Resources Research, 60(7), e2024WR037318. 2024.
Abstract
Abstract The constitutive relations of the Richardson-Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free-form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson-Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network-based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.BibTeX
@article{p2024_bandai-b, author = {Bandai, Toshiyuki and Ghezzehei, Teamrat A. and Jiang, Peishi and Kidger, Patrick and Chen, Xingyuan and Steefel, Carl I.}, title = {Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks}, journal = {Water Resources Research}, volume = {60}, number = {7}, pages = {e2024WR037318}, keywords = {inverse modeling, soil hydraulic functions, physics-informed machine learning, neural networks, soil moisture}, doi = {10.1029/2024WR037318}, pdf = {https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024WR037318}, note = {e2024WR037318 2024WR037318}, year = {2024}, research-theme = {machine-learning, water-flow} }
- Environmental Variability and Moisture-Temperature Coupling Reveal Continental-Scale Controls on Soil Respiration.Rojas, Y. T. P., & Ghezzehei, T. A.JGR-Biogeosciences (under Review).
Abstract
Understanding how temporal patterns of moisture variability control biogeochemical responses remains a fundamental challenge in Earth system science. Jensen’s inequality provides a mathematical framework for quantifying when episodic environmental events dominate over mean conditions. We applied this framework to continental-scale AmeriFlux data (134.5 million hourly observations from 2,004 soil moisture sensors) to quantify how moisture distribution patterns control soil respiration responses across environmental gradients. Sensors in dry regions show large Jensen’s inequality effects (median temporal averaging difference of -63.6%) because they experience highly skewed moisture distributions where brief wet periods drive disproportionate respiratory responses. Wet regions show minimal effects (median -27.1%) because they have more uniform moisture distributions. Data density analysis reveals that sensors operating at θ≈0.05 exhibit severe temporal averaging effects, while sensors at θ≈0.45 show minimal effects, demonstrating the mechanistic basis for where ecosystems operate on the moisture-respiration relationship. Climate gradient analysis shows systematic transitions from severe effects in arid systems to moderate effects in humid systems. Depth analysis reveals that surface soils experience maximum episodic event importance while deeper soils show reduced effects due to environmental buffering. Moisture-temperature coupling demonstrates systematic negative correlations in water-limited systems, indicating that environmental co-variation modulates biogeochemical responses. Jensen’s inequality emerges as a diagnostic tool for identifying when moisture variability patterns dominate biogeochemical processes, with continental-scale patterns revealing fundamental controls on episodic event importance across ecosystems.BibTeX
@article{2025-PerezRojas, title = {Environmental Variability and Moisture-Temperature Coupling Reveal Continental-Scale Controls on Soil Respiration}, language = {en}, journal = {JGR-Biogeosciences (under review)}, author = {Rojas, Yulissa T. Perez and Ghezzehei, Teamrat A.}, pages = {}, research-theme = {water-flow, soil-structure} } - Two Decades of Conservation Agriculture Enhances Soil Structure, Carbon Sequestration, and Water Retention in Mediterranean Soils.Alvarez-Sagrero, J., Berhe, A. A., Chacon, S. S., Mitchell, J. P., & Ghezzehei, T. A.SOIL (under Review).
Abstract
Conservation agriculture offers a pathway for enhancing soil health with climate co-benefits in Mediterranean agricultural systems. This study examined long-term impacts of combining no-till management with cover cropping over 20 years in California’s Central Valley, providing rare insights into soil system equilibrium under sustained conservation management. We assessed soil physical, chemical, and structural properties comparing reduced tillage with cover crops to standard tillage without cover crops, employing density fractionation and spectroscopic analysis to understand carbon protection mechanisms. After two decades, conservation agriculture achieved dynamic equilibrium characterized by fundamental shifts in carbon stabilization pathways. Water-stable aggregate analysis revealed the most pronounced management effects, with conservation practices exhibiting 136% greater stability, indicating substantial improvements in soil structural integrity. These structural enhancements corresponded with a reorganization of carbon protection mechanisms, demonstrating that physical protection within aggregates becomes a dominant carbon stabilization pathway under long-term conservation management. Mineral-associated organic carbon saturation analysis revealed that both management systems remained well below theoretical maximum capacity, indicating substantial remaining potential for carbon sequestration even after reaching equilibrium. Physical property improvements included 15% lower bulk density and 13% greater water retention at field capacity. Our findings demonstrate that two decades of conservation agriculture fundamentally transforms soil functioning through aggregate-mediated physical protection.BibTeX
@article{2025-AlvarezSagrero, title = {Two Decades of Conservation Agriculture Enhances Soil Structure, Carbon Sequestration, and Water Retention in Mediterranean Soils}, language = {en}, journal = {SOIL (under review)}, author = {Alvarez-Sagrero, Jennifer and Berhe, Asmeret Asefaw and Chacon, Stephany S. and Mitchell, Jeffrey P. and Ghezzehei, Teamrat A.}, pages = {}, research-theme = {sustainable-agriculture, soil-structure, water-flow} }
Related Research
Key Equipment
- HYPROP system for evaporation method
- Tempe pressure cells
- Hood infiltrometer
- Tensiometers and TDR probes
- Soil moisture sensors
Collaborate With Us
Interested in hydraulic property measurement or water flow modeling? Contact us to discuss research collaborations.
Get in Touch