SEARCH · Search NASA

Results for “groundwater infiltration”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

70 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2020

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Quantifying Subsurface Flow and Solute Transport in a Snowmelt‐Recharged Hillslope With Multiyear Water Balance

Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Insights into the thermo-hydraulic properties of compacted MX80 bentonite during hydration under elevated temperature

In high-level radioactive waste geological repositories, compacted bentonite undergoes coupled thermo-hydraulic processes due to heat released from a central waste canister and groundwater imbibition from the surrounding host rock. An understanding of these processes is essential for long-term simulations of radionuclide migration and canister corrosion, which requires an understanding of temperature effects on the coupled thermo-hydraulic properties governing these processes. In this study, a tank-scale radial infiltration test was used to investigate water imbibition processes in compacted bentonite under a central heater temperature of 200 °C that simulates high thermal gradients in a repository. Interpretation of this test focuses on evaluation of the liquid water wetting front during hydration and the interpretation of the temperature-dependent transient soil water retention curve (SWRC), thermal conductivity function (TCF), and hydraulic conductivity function (HCF). The SWRC during imbibition follows a temperature-dependent wetting path. Temperature effects on the hydraulic conductivity of bentonite in saturated conditions had the greatest effect on the shape of the HCF, with minimal temperature effects at higher suctions. The transient thermal conductivity data matched well with a new TCF linked with the SWRC shape. Although the bentonite layer was restrained, local deformations during hydration may have affected the shapes of the TCF and HCF.

Engineering↗

DOE OSTI2020

Phase II Geochemical and Isotopic Evaluation of Groundwater Flow in the Pahute Mesa–Oasis Valley Flow System, Nevada

This report presents the results of the Phase II geochemical evaluation for the Pahute Mesa CAUs. The evaluation was performed by multiple organizations including Los Alamos National Laboratory (LANL), Desert Research Institute (DRI), Lawrence Livermore National Laboratory (LLNL), and Navarro. Section 1.0 presents the introduction and describes the objectives for the Phase II geochemical evaluation. Section 2.0 presents the hydrologic setting, including average annual precipitation and net infiltration, hydraulic-head gradients, geology, and hydrostratigraphy of the PM-OV flow system.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Evaluating GCM land surface hydrology parameterizations by computing river discharges using a runoff routing model: Application to the Mississippi basin

To relate general circulation model (GCM) hydrologic output to readily available river hydrographic data, a runoff routing scheme that routes gridded runoffs through regional- or continental-scale river drainage basins is developed. By following the basin overland flow paths, the routing model generates river discharge hydrographs that can be compared to observed river discharges, thus allowing an analysis of the GCM representation of monthly, seasonal, and annual water balances over large regions. The runoff routing model consists of two linear reservoirs, a surface reservoir and a groundwater reservoir, which store and transport water. The water transport mechanisms operating within these two reservoirs are differentiated by their time scales; the groundwater reservoir transports water much more slowly than the surface reservior. The groundwater reservior feeds the corresponding surface store, and the surface stores are connected via the river network. The routing model is implemented over the Global Energy and Water Cycle Experiment (GEWEX) Continental-Scale International Project Mississippi River basin on a rectangular grid of 2 deg X 2.5 deg. Two land surface hydrology parameterizations provide the gridded runoff data required to run the runoff routing scheme: the variable infiltration capacity model, and the soil moisture component of the simple biosphere model. These parameterizations are driven with 4 deg X 5 deg gridded climatological potential evapotranspiration and 1979 First Global Atmospheric Research Program (GARP) Global Experiment precipitation. These investigations have quantified the importance of physically realistic soil moisture holding capacities, evaporation parameters, and runoff mechanisms in land surface hydrology formulations.

Liston, G. E.↗

DOE OSTI2023

Vadose Zone Soil Flushing for Chromium Remediation: A Laboratory Investigation to Support Field‐scale Application

Cr(VI) flushing from the vadose zone to the groundwater (with subsequent Cr(VI) removal in groundwater by pump-and-treat system) is a promising remedial technique that has recently been used at field scale. This laboratory study was conducted to provide the technical basis to design a field soil flushing strategy. The objectives were to (1) quantify the relationship between sediment Cr(VI) and Cr(III) mass and release rates and subsequent Cr(VI) leaching; (2) investigate different methodologies to maximize Cr(VI) leaching, and (3) investigate methods to minimize leaching of remaining residual Cr. Characterization of Cr-contaminated sediments (Hanford Site, WA) exhibited Cr(VI) showed that leach rates that were correlated to different Cr surface phases. Sediments with low leachable Cr(VI) (<2 μg/g) leached Cr rapidly, so slow infiltration of water in a single pulse was sufficient to leach most Cr. In contrast, sediments with high Cr (2 to 200 μg/g) released some Cr(VI) quickly but 10 to 50% Cr(VI) slowly (tens to hundreds of hours). Efficient unsaturated leaching of these sediments required a different infiltration strategy that includes: multiple slow leach pulses with time between flushing cycles; the use of a surfactant to increase Cr leaching from low-permeability zones, and the use of a reductant (Na-dithionite or Ca-polysulfide) in the final leach water was highly effective at decreasing residual Cr leaching. This study clearly demonstrated that the methodology of basing laboratory Cr flushing on parameters such as Cr release mass and rates could be used to improve the efficiency of soil flushing at field scale.

Accelerated Cr Flushing↗

NASA NTRS2022

Analyzing Recharge Dynamics and Storage in a Thick, Karstic Vadose Zone

Karst vadose zone heterogeneity creates complex transmission and storage dynamics that affect the timing and magnitude of aquifer recharge. Young, high-matrix permeability eogenetic karst aquifers may have significantly higher matrix storage than older, lower matrix permeability counterparts. In vulnerable and water-limited karst regions, the timescales of storage may be important for seasonal and sub-seasonal water resource management. We create a framework to quantify storage dynamics in karst aquifers using high-resolution precipitation and groundwater levels from the Northern Guam Lens Aquifer in the US territory of Guam. We estimate recharge using the Water Table Fluctuation method, and then develop transfer functions between precipitation and recharge to quantify storage and release of water from the vadose zone. The transfer functions are partitioned into different flow pathways including conduit, conduit/matrix, and slowly draining matrix. Probability distributions are fit to each pathway to determine the average travel times of infiltrated waters. The results show that aquifer recharge through secondary porosity features typically occurs within a few hours of a rainfall event, but this rapid recharge accounts for only 12%–28% of total recharge. The majority of aquifer recharge (>70%) occurs within a month of a contributing storm event. An additional 10% of recharge, on average, took longer than a month to reach the water table. The framework established can help improve the hydrological modeling and freshwater management for karst aquifers.

Patricia Spellman↗

DOE OSTI2021

Seismic Imaging of a Shale Landscape Under Compression Shows Limited Influence of Topography–Induced Fracturing

We used seismic refraction to image the P-wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill-valley-hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. Here, this contradiction could be explained by the greater importance of infiltration-driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.

58 GEOSCIENCES↗

DOE OSTI2023

Estimating Vadose Zone Flow Properties at the 100 K-East Soil Flushing Site Using ERT Monitoring Data: 2023 Interim Report - 100 KE Soil Property Estimation

In situ soil flushing is being using at the Hanford 100 K-East (100 KE) area to transport mobile chromium contamination in the vadose zone to the water table, where it can be collected and treat through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. In practice, it is infeasible to comprehensively monitor which regions of the vadose zone are being infiltrated through direct sampling of pore water. Consequently, there can be significant uncertainty about which regions of the subsurface have been treated, especially if hydrogeologic conditions are favorable for the development of unstable flows and preferred flow pathways through the vadose zone (Jarvis, Koestel, and Larsbo 2016). Current approaches for quantitative monitoring of soil flushing performance rely on contaminant concentration measurements collected from extractions wells. There is no quantitative information on the volume of flush water delivered to targeted regions of the vadose zone at the Hanford Site, leading to significant uncertainty regarding source term removal and long-term impacts to groundwater. If the subsurface hydrogeologic properties at the 100 KE Area were adequately known, qualitative metrics of soil flushing performance could be simulated, thereby negating expenses required to obtain quantitative performance information through borehole drilling/sampling. However, estimating in situ hydrogeologic properties has long proven elusive, due primarily to a lack of sufficient information to constrain heterogeneous property estimates to a useful degree of certainty. Estimating vadose zone hydrogeologic properties is particularly challenging due the dependence of hydraulic conductivity on saturation. This report describes progress toward a first-of-its-kind demonstration using surface time-lapse 3D electrical resistivity tomography (ERT) monitoring data to estimate the hydrogeologic properties that control flush water transport at the 100 KE soil flushing site. The ultimate objective is (1) to verify sufficient information exists in the ERT monitoring data to adequately resolve vadose zone hydraulic properties, and (2) generate a “digital twin” (i.e., a numerical simulator) that can be used to simulate the amount of flush water that has been delivered to each targeted region of the vadose zone, and thereby assess the efficacy of flush water delivery. Resulting performance estimates can be used in leu of comprehensive borehole drilling and direct sampling (or wellbore logging) that would otherwise be required to obtain the same information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1987

Valley development on Hawaiian volcanoes

Work in progress on Hawaiian drainage evolution indicates an important potential for understanding drainage development on Mars. Similar to Mars, the Hawaiian valleys were initiated by surface runoff, subsequently enlarged by groundwater sapping, and eventually stabilized as aquifers were depleted. Quantitative geomorphic measurements were used to evaluate the following factors in Hawaiian drainage evolution: climate, stream processes, and time. In comparing regions of similar climate, drainage density shows a general increase with the age of the volcani island. With age and climate held constant, sapping dominated valleys, in contrast to runoff-dominated valleys, display the following: lower drainage densities, higher ratios of valley floor width to valley height, and more positive profile concavities. Studies of stream junction angles indicate increasing junction angles with time on the drier leeward sides of the major islands. The quantitative geomorphic studies and earlier field work yielded important insights for Martian geomorphology. The importance of ash mantling in controlling infiltration on Hawaii also seems to apply to Mars. The Hawaiian valley also have implications for the valley networks of Martian heavily cratered terrains.

Baker, Victor R.↗

DOE OSTI2024

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗

DOE OSTI2023

Advances in 3D Geologic Modeling of Alluvial Basins with a Focus on Facies and Property Modeling

The unsaturated zone alluvium reference case is one of several geologic systems under consideration by the U.S. Department of Energy Office of Nuclear Energy for hosting repositories for spent nuclear fuel and associated waste (Sevougian et al., 2019). As noted by Mariner et al. (2018), the generic alluvial basin offers positive attributes that merit its consideration as a reference case by the Spent Fuel and Waste Science and Technology (SFWST) campaign. There are hundreds of alluvial basins and sub-basins scattered across the arid western United States (Figure 1-1). Precipitation and infiltration rates are relatively low with high evapotranspiration, resulting in vertical separation between repository and water table and thus longer transport paths to an aquifer. Accumulations of alluvial sediments within these basins are typically on the order of hundreds of meters, and locally may exceed 1,000 m in thickness, as is the case for the Deming sub-basin in southern New Mexico. A thick geologic host medium, which serves as the natural barrier system (NBS) in the conceptual model framework of a geologic disposal system, isolates the waste packages from receptors in the biosphere. Further, alluvial basin fill is typically comprised of stacked playa and lacustrine deposits along the basin axis (Perry et al., 2018). Characterized by low permeability, these layers protect the biosphere above the repository and the groundwater resources below the repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗

DOE OSTI2025

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

Abstract 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. Core Ideas Soil moisture data can be used to identify the occurrence of preferential flow (PF) and its initiating conditions. Various data‐analysis methods to identify PF differ in susceptibility to error, bias, and subjectivity. These methods can utilize vast amounts of data from soil moisture monitoring networks to develop understanding of when, where, and under what conditions PF occurs. Newly developed methods may lead to better accuracy and reliability, and reduce the need for subjective judgments. Plain Language Summary Preferential flow through soil occurs when a large amount of water is suddenly available, as during an intense storm. This type of flow moves rapidly through the soil in distinct narrow pathways rather than moving evenly throughout the body of soil, with major consequences for groundwater resources, ecosystems, spreading of contaminants, and other vital concerns. Methods of detecting preferential flow have been developed that utilize measurements of soil water content made by sensors installed at various depths. This measurement technology has been widely implemented, many locations now having datasets years in length, and various methods have been developed for using these to identify preferential flow. The various methods are based on different features in the soil moisture records and vary in their advantages and shortcomings. In this review, we explain and evaluate these methods, highlighting important information for their implementation to identify preferential flow from soil moisture data and for efforts to develop new methods or improve existing ones.

Nimmo, John R↗

DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗

DOE OSTI2026

Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes

Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099

Li, Bonan↗