SEARCH · Search NASA
Results for “hydraulic conductivity”
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Impact of soil surface and subsurface properties on soil saturated hydraulic conductivity in the semi-arid Walnut Gulch Experimental Watershed, Arizona, USA
Explore the source record for details and available documents.
Direct Observation of the Depth of Active Groundwater Circulation in an Alpine Watershed
Abstract The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3 H/ 3 He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization.
Effect of initial water saturation on the performance of fracturing fluids with and without polyallylamine under simulated EGS conditions
Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO 2 - reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO 2 , can be used as a potentially less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO 2 fracturing fluids outperformed other fluids such as water, CO 2 , and CO 2 /water in generating large fractures with less fluid consumed. The objective of this work is to investigate the effect of initial water saturation on the performance of StimuFrac by conducting hydraulic fracturing tests with ½ foot cubic rock samples held under representative EGS stress/temperature conditions and by using cyclic injection strategies (under constant injection rate). The resulting fracture hydraulic conductivities, breakdown pressures, and volumes of fluids required are compared. Methods/Procedures/Process: To simulate geothermal reservoir conditions, in all tests, the rock sample was held under triaxial confinement and at 200 °C, and different volumes of water were initially injected into the rock sample before any fracturing processes were initiated. For the single-cycle PAA (or water) alternating CO 2 (PAG or WAG) injection fracturing experiments, one complete cycle consisted of two steps: (1) injecting a PAA slug (or water slug) followed by (2) injecting CO 2 to initiate and propagate the fracture. For experiments involving multiple injection cycles, the CO 2 injection pressure is increased until it peaks and begins to decline (indicating fracture initiation at this moment), and then continued being injected for another 30 seconds to propagate the fracture. Then, these two-step cycles [injection of PAA (or water) followed by CO 2 injection (up to 2-4 mL/min)] are repeated. Applications/Significance/Novelty: The results of this study suggest that water saturation significantly affects the fracturing fluid transmission into the rock pore space, thus affecting the fracture initiation and propagation. In this study, fracturing tests via a single injection cycle or multiple injection cycles were performed. Splitting the rock samples in half after testing reveals that fracture propagation is significantly limited under high water saturation conditions (three-day initial water injection) compared to stimulation experiments performed in hot dry rock. The fractures propagate less than 1/3 of the distance from the wellbore to the outer rock surface, and in some cases, no fracture is generated. This may be caused by leak-off dominating the fracturing process and the fluid injection rate is insufficient to overcome leak-off, even under high injection rate conditions. Additionally, CO 2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate sufficiently high-pressure gradients away from the near-wellbore region. Under low water saturation conditions (dry rock or after 1-day initial water injection), PAA/CO 2 consistently generated significantly larger fractures compared with the other fluids. CO 2 generated large fractures only in the hot dry rock and only when using high injection rates, though data variability is high.
Lianas and Trees From a Seasonally Dry and a Wet Tropical Forest Did Not Differ in Embolism Resistance but Did Differ in Xylem Anatomical Traits in the Dry Forest
One of the most prominent changes in neotropical forests has been the increase in abundance and size of lianas. Studies suggest that lianas have more acquisitive strategies than trees, which could allow them to take advantage of water more effectively when it is available in water-limited forests, but few studies compared across growth form (i.e., lianas vs. trees) and forest type (i.e., wet vs. seasonally dry). We measured hydraulic and anatomical traits of co-occurring lianas and trees that convey drought resistance (xylem embolism resistance and intervessel pit membranes) and water transport capacity (xylem vessel diameter and density) in a seasonally dry and a wet evergreen tropical forest to address: (1) Are there differences between vulnerability to embolisms (P 50 —water potential at 50% loss of hydraulic conductivity) and hydraulic safety margins (HSM) across growth form and forest type? (2) How do vessel diameter and density vary across growth form and forest type? (3) Are there differences in xylem intervessel pit membrane thickness across growth form and forest type and does it predict xylem embolism vulnerability in trees and lianas? We examined hydraulic and xylem anatomical traits of 32 species—eight lianas and eight trees in each forest type. We found no difference in P 50 and HSMs between lianas and trees and between the wetter and drier forest. Dry forest lianas had 81% greater maximum vessel diameter and 125% greater range in vessel diameter sizes than dry forest trees but, there was no significant difference between life forms in the wet forest. Dry forest species had 50% greater vessel density and 30% greater maximum pit membrane thickness than wet forest ones. Maximum pit membrane thickness was correlated to P 50 and HSMs. The main difference between lianas and trees occurred in the dry forest, where lianas had larger maximum xylem vessel size than trees, implying that they have proportionally greater hydraulic conductive capacity than the trees in seasonal forests.
Physics-Informed Neural Network Method for Forward and Backward Advection-Dispersion Equations
Advection-dispersion equations (ADEs) are commonly used to describe transport phenomena in porous media. Even though mature discretization-based numerical methods for ADEs exist, some challenges still remain, especially when it comes to solving advection-dominated forward ADEs and diffusion-dominated backward ADEs. The latter problem usually arises in the source identification context and leads to numerically unstable grid-based solutions that require a form of regularization or should be treated as an inverse problem that is computationally more expensive because it requires solving the forward problem multiple times. In this study, we propose a discretization-free approach based on the physics-informed neural network (PINN) method for solving coupled ADE and Darcy flow equations with space-dependent hydraulic conductivity. In this approach, the hydraulic conductivity, hydraulic head, and concentration fields are approximated with deep neural networks (DNNs). We assume that the conductivity field is given by its values on a grid, and we use these values to train the conductivity DNN. The head and concentration DNNs are trained by minimizing the residuals of the flow equation and ADE and using the initial and boundary conditions as additional constraints. The PINN method is applied to one- and two-dimensional forward ADE problems, where its performance for various P\'{e}clet numbers ($Pe$) is compared with the analytical and numerical solutions. We find that the PINN method is accurate with errors of less than 1\% and outperforms some conventional discretization-based methods for $Pe$ larger that 100. Next, we demonstrate that the PINN method remains accurate for the backward ADEs, with the relative errors in most cases staying under 5\% compared to the reference concentration field. Finally, we show that when available, the concentration measurements can be easily incorporated in the PINN method and significantly improve (by more than 50\% in the considered cases) the accuracy of the PINN solution of the backward ADE.
Leaves as bottlenecks: The contribution of tree leaves to hydraulic resistance within the soil–plant–atmosphere continuum
Within vascular plants, the partitioning of hydraulic resistance along the soil-to-leaf continuum affects transpiration and its response to environmental conditions. In trees, the fractional contribution of leaf hydraulic resistance (R leaf ) to total soil-to-leaf hydraulic resistance (R total ), or fR leaf (= R leaf /R total ), is thought to be large, but this has not been tested comprehensively. Here, we compiled a multi-biome dataset of fR leaf using new and previously published measurements of pressure differences within trees in situ. Across 80 samples, fR leaf averaged 0.51 (95% confidence interval [CI] = 0.46-0.57) and it declined with tree height. We also used the allometric relationship between field-based measurements of soil-to-leaf hydraulic conductance and laboratory-based measurements of leaf hydraulic conductance to compute the average fR leaf for 19 tree samples, which was 0.40 (95% CI = 0.29-0.56). The in-situ technique produces a more accurate descriptor of fR leaf because it accounts for dynamic leaf hydraulic conductance. Both approaches demonstrate the outsized role of leaves in controlling tree hydrodynamics. A larger fR leaf may play a role in protecting stems from loss of hydraulic conductance. Thus, the decline in fR leaf with tree height would contribute to greater drought vulnerability in taller trees and potentially to their observed disproportionate drought mortality.
Characterization of an electrically conductive proppant for fracture diagnostics
Fracture diagnosis with electromagnetic (EM) and electrical tools requires proppants with high electrical conductivity and mechanical strength. Lab measurements of the electrical and hydraulic conductivity of proppants are critical for selecting the best candidates. Such measurements greatly benefit simulations, field tests, and the ultimate application of such proppants in the field. To that end, a new lab protocol is developed for measuring the electrical and hydraulic conductivity of proppants. The lab setup, which mainly includes a resistivity core holder and a Hassler sleeve core holder, allows for simulation of realistic pressure and temperature conditions when making measurements. Petroleum coke (PC) is proposed as a candidate proppant because of its widespread availability and low cost. Lab measurements show that the effective electrical conductivity of pure PC in a model fracture is approximately 5000 S/m, under a closure stress greater than [Formula: see text] (4000 psi). When PC is mixed with sand, the effective electrical conductivity of the mixture decreases with an increasing weight percentage of sand. Although sand degrades the contact between PC particles, the electrical conductivity stays reasonably high (approximately 1700 S/m) when 50% sand is added. Hydraulic conductivity measurements show that when a fracture is propped with pure PC, the measured fracture conductivity is greater than [Formula: see text] ([Formula: see text]) (dimensionless fracture conductivity greater than 100 for a shale with [Formula: see text] or 100 nD permeability) under a confining pressure of [Formula: see text] (6000 psi). This means that a fracture propped with PC is infinitely conductive in a typical shale formation. When sand is added, the fracture’s hydraulic conductivity becomes even higher, which clearly shows PC’s ability of sustaining high stresses. The proposed protocol provides a robust and effective method that can be generalized for lab testing for other candidate proppants. The data presented clearly show that PC has the potential for field-scale applications in EM hydraulic fracture diagnostics.
Response of the water status of soybean to changes in soil water potentials controlled by the water pressure in microporous tubes
Water transport through a microporous tube-soil-plant system was investigated by measuring the response of soil and plant water status to step change reductions in the water pressure within the tubes. Soybeans were germinated and grown in a porous ceramic 'soil' at a porous tube water pressure of -0.5 kpa for 28 d. During this time, the soil matric potential was nearly in equilibrium with tube water pressure. Water pressure in the porous tubes was then reduced to either -1.0, -1.5 or -2.0 kPa. Sap flow rates, leaf conductance and soil, root and leaf water potentials were measured before and after this change. A reduction in porous tube water pressure from -0.5 to -1.0 or -1.5 kPa did not result in any significant change in soil or plant water status. A reduction in porous tube water pressure to -2.0 kPa resulted in significant reductions in sap flow, leaf conductance, and soil, root and leaf water potentials. Hydraulic conductance, calculated as the transpiration rate/delta psi between two points in the water transport pathway, was used to analyse water transport through the tube-soil-plant continuum. At porous tube water pressures of -0.5 to-1.5 kPa soil moisture was readily available and hydraulic conductance of the plant limited water transport. At -2.0 kPa, hydraulic conductance of the bulk soil was the dominant factor in water movement.
CMP Pits and D-Area Watershed Infiltrometer Testing
Savannah River National Laboratory (SRNL) conducted saturated hydraulic conductivity field testing at the Chemical, Metals, and Pesticide Pits (CMP Pits) and the D-Area watershed to develop infiltration rates for use in future updates to the CMP Pits groundwater flow model. Approximately 18 testing locations were selected at each site based on the soil series that were identified in the Soil Survey of Savannah River Plant Area, Parts of Aiken, Barnwell, and Allendale Counties, South Carolina.
CMP Pits and D-Area Watershed Infiltrometer Testing
Savannah River National Laboratory (SRNL) conducted saturated hydraulic conductivity field testing at the Chemical, Metals, and Pesticide Pits (CMP Pits) and the D-Area watershed to develop infiltration rates for use in future updates to the CMP Pits groundwater flow model. Approximately 18 testing locations were selected at each site based on the soil series that were identified in the Soil Survey of Savannah River Plant Area, Parts of Aiken, Barnwell, and Allendale Counties, South Carolina.
Index tests for geosynthetic clay liners containing bentonite–polymer composites
The hydraulic conductivity of geosynthetic clay liners (GCLs) containing bentonite–polymer composites (BPCs) is controlled by a pore clogging mechanism related to the swelling of sodium bentonite (NaB) granules and formation of polymer gel in the pore space. Three index testing procedures that are influenced by bentonite swelling and rheology of the polymer gel were evaluated in this study for screening BPC GCLs for hydraulic conductivity: viscosity, liquid limit by fall-cone penetration and active water content. Experiments were conducted using polyacrylamide (PAM) polymer, dry-blended mixtures of PAM and NaB and a BPC from a commercially available GCL. These materials were contacted with inorganic solutions representing leachates from a range of waste containment scenarios. For the BPC that was evaluated, hydraulic conductivities less than 10 -10 m/s were obtained when the viscosity exceeded 220 mPa s and the active water content exceeded 300%. Hydraulic conductivities less than 10 -11 m/s were achieved for viscosity exceeding 350 mPa s, fall-cone liquid limit exceeding 190 and active water content exceeding 300%. Additional testing is needed to assess the generality of these index tests and the identified thresholds when applied to other BPC GCLs and other permeant liquids.
3D hydrogeophysical characterization of managed aquifer recharge basins
Aquifers are increasingly stressed. Managed aquifer recharge provides a potential solution to mitigate this stress and provide sustainable groundwater resources. Subsurface properties are known to have a strong control on the infiltration rates that can be achieved. However, these properties are often highly heterogeneous and difficult to assess with conventional probing techniques. Here, we show the application of 3D geophysical imaging to assess the recharge potential and its variation across several basins used for managed aquifer recharge. We link in-situ measurements of saturated hydraulic conductivity with the electrical resistivity of the subsurface to establish petrophysical relationships and use those relationships to estimate the distribution of hydraulic conductivity throughout the five recharge basins. Our results show a considerable variability in the hydraulic properties, i.e., soil texture and saturated hydraulic conductivity, that have a direct impact on potential infiltration rates. We use the 3D hydraulic property distributions to model groundwater recharge and provide estimates for infiltration rates and volumes, and use this approach to assess the impact of management activities on groundwater recharge performance. Having such data not only enables us to predict infiltration rates, but also provides means for optimizing such water infrastructure.
New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles
Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.
Physical and Hydraulic Properties of RCRA Borehole Samples from the Hanford Site : Final Report, Fiscal Years 2023-2024
Sediment from 44 core samples collected from Resource Conservation and Recovery Act (RCRA) boreholes drilled in the 200 East and 200 West areas of the Hanford Site were characterized for physical and hydraulic properties (Table S.1). Characterization data included gravimetric water contents and matric potentials, grain-size distributions, saturated hydraulic conductivity, water retention characteristics, and unsaturated hydraulic conductivity. These properties provide site-specific data and parameters that can be used in subsurface flow and contaminant transport models to assess the transport and fate of contaminants in the vadose zone and underlying aquifer systems. The analyzed core samples come from specific areas and depth intervals at the Hanford Site that were targeted for sampling to address data gaps identified by site contractors (Khaleel 2020). X-ray computed tomography (XCT) was used to evaluate the general textural characteristics of the samples and to determine which samples to use for further physical and hydraulic property characterization. Subsequent sample selection was determined by consensus after review of the XCT images by Pacific Northwest National Laboratory, Central Plateau Cleanup Company, and INTERA staff.
Morphological controls on flow conductivity and viscosity of bentonite–polymer composites
In this study, the authors investigate the morphological organisation of polymers in bentonite–polymer composites (BPCs) used in geosynthetic clay liners and controls on hydraulic conductivity and viscosity. A limited understanding of the relationship between the microstructure of BPCs and their macroscopic engineering behaviour motivates this study. The polymer resides in the inter-granular pores of bentonite based on synchrotron-based small-angle X-ray scattering measurements, scanning electron microscopy imaging and pore size analyses. On wetting, the polymer swells to produce a swollen hydrogel residing in the inter-granular pores. Enhanced tortuosity arising from hydrogel formation reduces the hydraulic conductivity of BPCs. The higher viscosity of the BPC slurry is attributed to the higher shear resistance emerging from swelling of the polymer fraction on wetting. These findings link the morphological features of BPCs to the observed hydraulic conductivity and viscosity of these materials.
Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau
The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.