In the past decades, remarkable progresses have been made in developing soil water retention (SWR) functions and relative hydraulic conductivity (
Despite the vast extent of desert soils on the earth's surface, our understanding of the moisture dynamics of near‐surface desert soils (i.e., the top centimeters to few meters of the soil profile) remain limited. The goal of this study was to explore the use of the Peters–Durner–Iden (or PDI) instead of bimodal van Genuchten (or BVG) hydraulic functions to improve water redistribution simulations using HYDRUS‐1D for drier soils in desert environments. The PDI hydraulic functions take capillary and film flow into account, whereas BVG hydraulic functions are limited to capillary flow. By comparing measured with simulated water content data, we found that moisture redistribution simulations were improved by using PDI instead of BVG soil water retention and hydraulic conductivity functions. Compared with the BVG simulations, the PDI simulations particularly improved for drier soil conditions (i.e., volumetric water contents ranging from 6 to 10%; suction heads between pF 2 and pF 3.8, and saturation degrees between 19 and 32%, respectively) for the studied sandy soil of Scaling Environmental Processes in Heterogeneous Arid Soils (SEPHAS) Lysimeter 1. For pF >3, the PDI functions predicted higher hydraulic conductivity than the BVG functions, which confirmed the hypothesis that a hydraulic conductivity function, which can capture film flow, may improve moisture distribution simulations for dry soils. For pF between 2 and 3, however, simulation results improved due to the difference in the water retention rather than the hydraulic conductivity function.
more » « less- NSF-PAR ID:
- 10448359
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Vadose Zone Journal
- Volume:
- 19
- Issue:
- 1
- ISSN:
- 1539-1663
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract k r) models for the full range of matric suction. However, such an existing SWR function and the correspondingk rmodel were found not to have a consistent prediction performance, that is, for a given soil, the SWR function could perform very well but thek rmodel might perform poorly or even fail. Given that both water retention and hydraulic conductivity need to be accurately predicted to simulate land‐atmosphere interactions, vadose‐zone soil–water dynamics, and shallow groundwater recharges, this paper proposed an upgraded SWR function and derived a new closed‐formk rmodel by combining the function with the well‐known Mualem model. The upgradation was done by introducing an exponential function for the pendular (i.e., dry and very dry) range of matric suction. The SWR function andk rmodel keep the conventional simplicity, reflect the capillary versus adsorptive soil–water mechanisms, eliminate fictitious parameters, and are continuously differentiable with closed forms. The assessment using the experimental data for 35 soils with nine different textures indicated that the upgraded SWR function and newk rmodel reproduced the data more accurately than the existing functions and models, as measured by visualization plots and four statistics, namely coefficient of determination, sum squared error, root mean squared error, and Nash–Sutcliffe efficiency. The contribution of this paper is to overcome the inconsistent‐performance issue mentioned above. -
The Caldor fire burned ~222,000 acres of the Eastern Sierra Nevada during summer–fall 2021. We evaluated the effects of this “megafire” on the physical properties of a sandy soil developed from glacial tills to document fire-induced soil modifications in this region. We measured soil water retention and hydraulic conductivity functions as well as the thermal properties of five core samples from control (unburned) areas and eight core samples from burned soil of the same soil unit. Soil water repellency was measured in terms of water drop penetration time (WDPT) in the field and apparent contact angle in the laboratory on control and burned soil as well as ash samples. Soil organic matter (SOM) and particle and aggregate size distributions were determined on control and burned soil samples. Additionally, scanning electron microscopy (SEM) was used to image microaggregates of control and burned soil samples. We found a significant difference in SOM content and sand and silt aggregate size distribution between control and burned samples, which we associated with the disintegration of microaggregates due to the fire. We found no significant difference between soil water retention and hydraulic conductivity functions of control and burned soil but observed greater variation in saturated hydraulic conductivity and systematic shifts in thermal conductivity functions of burned compared to control samples. WDPT and apparent contact angle values were significantly higher for burned soils, indicating the occurrence of fire-induced soil hydrophobicity (FISH). Interestingly, the average apparent contact angle of the control soil was >90°, indicating that even the unburned soil was hydrophobic. However, the ash on top of the burned soil was found to be hydrophilic, having apparent contact angles <10°. Our results indicate that SOM and microaggregates were readily affected by the Caldor fire, even for sandy soil with a weakly developed structure. The fire seemed to have moderated thermal properties, significantly and soil wettability but had only minimal effects on water retention and hydraulic conductivity functions. Our findings demonstrate the complex nature of fire-soil interactions in a natural environment and highlight the need for additional investigation into the causes and processes associated with FISH and structure alterations due to fire to improve our ability to rapidly determine potential problem areas in terms of hazards commonly associated with fire-affected soils.more » « less
-
more » « less
Abstract The capture and use of water are critically important in drylands, which collectively constitute Earth's largest biome. Drylands will likely experience lower and more unreliable rainfall as climatic conditions change over the next century. Dryland soils support a rich community of microphytic organisms (biocrusts), which are critically important because they regulate the delivery and retention of water. Yet despite their hydrological significance, a global synthesis of their effects on hydrology is lacking. We synthesized 2,997 observations from 109 publications to explore how biocrusts affected five hydrological processes (times to ponding and runoff, early [sorptivity] and final [infiltration] stages of water flow into soil, and the rate or volume of runoff) and two hydrological outcomes (moisture storage, sediment production). We found that increasing biocrust cover reduced the time for water to pond on the surface (−40%) and commence runoff (−33%), and reduced infiltration (−34%) and sediment production (−68%). Greater biocrust cover had no significant effect on sorptivity or runoff rate/amount, but increased moisture storage (+14%). Infiltration declined most (−56%) at fine scales, and moisture storage was greatest (+36%) at large scales. Effects of biocrust type (cyanobacteria, lichen, moss, mixed), soil texture (sand, loam, clay), and climatic zone (arid, semiarid, dry subhumid) were nuanced. Our synthesis provides novel insights into the magnitude, processes, and contexts of biocrust effects in drylands. This information is critical to improve our capacity to manage dwindling dryland water supplies as Earth becomes hotter and drier.
-
more » « less
Abstract The heat transfer and water retention in soils, governed by soil thermal conductivity (λ) and soil water retention curve (SWRC), are coupled. Soil water content (θ) significantly affects λ. Several models have been developed to describe λ(θ) relationships for unsaturated soils. Ghanbarian and Daigle presented a percolation‐based effective‐medium approximation (P‐EMA) for λ(θ) with two parameters: scaling exponent (
t s) and critical water content (θc). In this study, we explored the new insights into the correlation between soil thermal conductivity and water retention using the P‐EMA and van Genuchten models. The θcwas strongly correlated to selected soil hydraulic and physical properties, such as water contents at wilting point (θpwp), inflection point (θi), and hydraulic continuity (θhc) determined from measured SWRCs for a 23‐soil calibration dataset. The established relationships were then evaluated on a seven‐soil validation dataset to estimate θc. Results confirmed their robustness with root mean square error ranging from 0.011 to 0.015 cm3cm−3, MAE ranging from 0.008 to 0.013 cm3cm−3, andR 2of 0.98. Further discussion investigated the underlying mechanism for the correlation between θcwith θhcwhich dominate both heat transfer and water flow. More importantly, this study revealed the possibility to further investigate the general relationship between λ(θ) and SWRC data in the future. -
more » « less
Abstract A density‐dependent, variably saturated groundwater flow and solute transport model was used to investigate the influence of swash motions on subsurface flow and moisture dynamics in beach aquifers with heterogeneous distributions of hydraulic conductivity (
K ) and capillarity. The numerical simulations were performed within a Monte Carlo framework using field measurements conducted in the swash zone of a sandy beach. Our results show that heterogeneous capillarity causes spatially variable capillary rise above the groundwater table. In response to swash motions, heterogeneity creates capillary barriers that result in pockets of elevated moisture content beneath the swash zone. These moisture hotspots persist within the unsaturated zone even at ebb tide when the swash motions recede seaward. Heterogeneous capillarity also results in highly tortuous preferential flow paths and alters the flow rates from the sand surface to the water table. HeterogeneousK greatly enhances the seawater infiltration into the swash zone and modulates its spatial distribution along the beach surface. Due to heterogeneousK and capillarity, complex mixing patterns emerge. Both strain‐dominated and vorticity‐dominated flow regions develop and dissipate as tides and waves move across the beach surface. Complex mixing patterns of seawater percolating from the swash zone surface to the water table, with localized areas of high and low mixing intensities, are further demonstrated by analysis of dilution index. Our findings reveal the influence of geologic heterogeneity on swash zone moisture and flow dynamics, which may have important implications for sediment transport and chemical processing in beach aquifers.
