Search NASA⌕ Search

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

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.

Smith-Martin, Chris M.↗
DOE OSTI2021

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.

He, Qizhi↗
DOE OSTI2023

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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.

Geochemistry & Geophysics↗
DOE OSTI2025

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

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.

Engineering↗
DOE OSTI2022

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

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.

Harp, Dylan R. (ORCID:0000000197778000)↗
DOE OSTI2025

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.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

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.

Engineering↗
DOE OSTI2025

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.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Convergent evolution of tree hydraulic traits in Amazonian habitats: implications for community assemblage and vulnerability to drought

Amazonian droughts are increasing in frequency and severity. However, little is known about how this may influence species-specific vulnerability to drought across different ecosystem types. We measured 16 functional traits for 16 congeneric species from six families and eight genera restricted to floodplain, swamp, white-sand or plateau forests of Central Amazonia. We investigated whether habitat distributions can be explained by species hydraulic strategies, and if habitat specialists differ in their vulnerability to embolism that would make water transport difficult during drought periods. We found strong functional differences among species. Nonflooded species had higher wood specific gravity and lower stomatal density, whereas flooded species had wider vessels, and higher leaf and xylem hydraulic conductivity. The P 50 values (water potential at 50% loss of hydraulic conductivity) of nonflooded species were significantly more negative than flooded species. However, we found no differences in hydraulic safety margin among species, suggesting that all trees may be equally likely to experience hydraulic failure during severe droughts. Water availability imposes a strong selection leading to differentiation of plant hydraulic strategies among species and may underlie patterns of adaptive radiation in many tropical tree genera. Our results have important implications for modeling species distribution and resilience under future climate scenarios.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Physics-embedded inverse analysis with algorithmic differentiation for the earth’s subsurface

Abstract Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Unsaturated and Saturated Hydraulic Parameters of Soils at the East River Watershed, Colorado Calculated using the Pedotransfer Function Method.

The purpose of the data package is to present the results of calculations of the following unsaturated-saturated hydraulic parameters of soils of the East River Watershed, Colorado: fully saturated water content (Qs), residual saturation (Qr), van genuchten parameters alpha and n, field capacity (FC), and saturated hydraulic conductivity (Ks). Calculations were conducted for the following soil texture types: sandy_loam, loam, silty_clay_loam, clay_loam, silt loam; the particle size distribution included the percentage of sand, silt, clay, and the bulk density (BD) of soils. Calculations were conducted for the following areas- Lower_Subalpine (6 samples),- Lower_Montane_lower_hillslope (7 samples),- Floodplain_at_base_of_lower_montane_hillslope, (11 samples),- PLM-0 (10 samples),- PLM-1 (10 samples), and - PLM-3 (13 samples)Calculations were conducted using the Pedotransfer Function Method developed by Tóth et al. (2014) with the application of the code ”euptf” in Rstudio. The data are being summarized in the Excel file ”Unsaturated_Saturated_Hydraulic_Parameters_PTF_calculated_EASt_EIver.xls" including two worksheets:“full_data_EAST_RIVER” which includes (a) the input data used for calculations, the soil texture triangle, and the plots of the water retention curves and unsaturated hydraulic conductivity for several types of soil texture, and (b) the summary table with the averaged parameters for each area. The the texture triangle was plotted using the Rstudio "soiltexture." Microsoft Excel software is required to open the Excel file. The metadata and a list of references to the Rstudio libraries used for calculations and data visualization are giving in the Word Readme file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling Monthly Nitrate Concentration in a Karst Spring with and without Discrete Conduit Flow

Understanding the coupled continuum pipe-flow framework for modeling contaminant transport in karst systems is critical for protecting water resources therein. This study simulated point and non-point source nitrate concentration in a karst spring and investigated the results generated from the flow and transport model with and without discrete conduit flow. CFPv2 and CMT3D models were integrated to address the changes in nitrate concentration at a monthly scale in a karst spring, and the results were compared with that from an equivalent porous media (EPM) model with high-hydraulic conductivity (K) zones set in the main karstified area to represent conduits. The results show that the CFPv2+CMT3D model is able to describe well the recession of nitrate concentration in spring discharge, and the relatively larger deviation (slower nitrate recession) from the observed trend for the EPM model is probably a result of the limitation of utilizing high hydraulic conductivity cells to represent conduit. Moreover, simulated hydraulic heads in poorly karstified areas from the two models both show slight differences from the observations (the head RMSE values of calibration/validation for CFPv2 and MODFLOW models are 0.16 m/0.25 m and 0.26 m/0.17 m, respectively), indicating the inclusion of conduits may not affect the simulation considerably, and the lower the proportion of karstic area, the slight effects brought from the inclusion of conduits in the model. For highly karstified areas, the CFPv2+CMT3D model could provide more accurate results (head RMSE of calibration/validation for CFPv2 and MODFLOW are 0.22 m/0.06 and 0.52 m/0.47 m, respectively), showing the coupled continuum pipe-flow framework may be more appropriate for applying to highly and maturely karstified areas where the variations in the behavior of flow and contaminant transport are more affected by turbulent flow regime.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Regional Subsidence Analysis Through a Multi-Scale Modeling Framework Based on Breakage Mechanics

Although poroelastic models are often used to explain the delay between subsurface fluid depletion and ground subsidence, inelastic compaction involving permanent changes of rock microstructure may exacerbate hydro-mechanical coupling, thus influencing the interpretation of measurements and long-term forecasts. Here, a multi-scale modeling approach is discussed, which accounts for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. A constitutive model built within the framework of breakage mechanics is proposed to link the hydraulic conductivity of granular rocks with inelastic deformations and changes in grading caused by injection-depletion cycles at stress levels far from yielding. The proposed model has been incorporated into large-scale simulation frameworks, thus enabling the spatiotemporal mapping of regional subsidence through a hybrid, semi-analytical approach. Numerical results based on this strategy show that the model allows isolating near-field and far-field effects into the computation of land subsidence and can generate forecasts for different modeling scenarios (e.g., elastic and inelastic compaction, constant permeability, and concurrent change of compressibility and permeability). In particular, examples of simulations for the case of the Groningen gas field are discussed, showing the model capabilities to use both field measurements and laboratory tests for the generation of reasonable subsidence maps, without expensive computational costs. Results indicate that ignoring coupled inelastic effects has major consequences on the predicted timescale of subsidence. Specifically, while all the model scenarios produced similar long-term ground settlements, those ignoring breakage-dependent permeability changes result in a variation of the temporal window of residual subsidence of the order of several decades.

58 GEOSCIENCES↗