SEARCH · Engineering Papers

Results for “hydrogeology”

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.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Improved geophysical monitoring of carbon sequestration through parameter linkage to reservoir modeling

Predictive reservoir modeling, even if present in the form of only basic hydrogeological model assumptions, is expected to accompany the majority of carbon capture and sequestration monitoring activities. It thus represents a source of prior information about the migration of injected fluids that can benefit geophysical survey planning and ensuing monitoring. Constraining the imaging of geophysical monitoring data with reservoir modeling is preferable over standalone geophysical imaging because of additional complementary hydrogeological information. However, fully coupled hydrogeophysical data inversion for flow-modeling parameters that control saturation predictions is an involved process. Within the context of three-dimensional electromagnetic (EM) inversion of data from borehole-to-surface layouts, we employ a "poor people's" alternative. The approach constrains geophysical inversion parameters through saturation predictions. The coupling is realized through spatially variable lower and upper parameter bounds that scale with gas saturation magnitudes, the latter provided by reservoir modeling. In this work, enhancement of three-dimensional time-lapse plume EM imaging is demonstrated for simulated sequestration into a depleted gas reservoir.

58 GEOSCIENCES↗

DOE OSTI2023

Continental Scale Hydrostratigraphy: Basin‐Scale Testing of Alternative Data‐Driven Approaches

Abstract Integrated hydrological modeling is an effective method for understanding interactions between parts of the hydrologic cycle, quantifying water resources, and furthering knowledge of hydrologic processes. However, these models are dependent on robust and accurate datasets that physically represent spatial characteristics as model inputs. This study evaluates multiple data‐driven approaches for estimating hydraulic conductivity and subsurface properties at the continental‐scale, constructed from existing subsurface dataset components. Each subsurface configuration represents upper (unconfined) hydrogeology, lower (confined) hydrogeology, and the presence of a vertical flow barrier. Configurations are tested in two large‐scale U.S. watersheds using an integrated model. Model results are compared to observed streamflow and steady state water table depth (WTD). We provide model results for a range of configurations and show that both WTD and surface water partitioning are important indicators of performance. We also show that geology data source, total subsurface depth, anisotropy, and inclusion of a vertical flow barrier are the most important considerations for subsurface configurations. While a range of configurations proved viable, we provide a recommended Selected National Configuration 1 km resolution subsurface dataset for use in distributed large‐and continental‐scale hydrologic modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Natural Systems Data Package for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

This document summarizes the available information regarding the natural system (vadose zone and saturated zone) at the active trenches of the 200 West Area and 200 East Area Low-Level Burial Grounds (LLBGs). This information is relevant for developing the natural system fate and transport model for the update of the active trenches performance assessment. This document presents a description of the major features, events, and processes of the natural system at the active trenches of the 200 West Area and 200 East Area LLBGs that involves the geologic features or framework, and the vadose zone and groundwater flow and transport system. The description of the geologic features includes the following items: geologic setting, major structural features, major hydrostratigraphic units, and geologic framework. The discussion concerning vadose zone flow and transport at the active trenches of the 200 West Area and 200 East Area LLBGs includes the following topics: vadose zone features, events, and processes relevant to the development of the vadose zone conceptual model; hydrogeologic units and conditions; factors affecting moisture movement; and flow and transport properties of major units. The discussion concerning the groundwater flow and transport system in the vicinity of the active trenches of the 200 West Area and 200 East Area LLBGs includes the following topics: saturated zone features, events, and processes relevant to the saturated zone conceptual model; hydrogeologic units; groundwater flow conditions; and flow and transport properties.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Alpha Spectrometry Results for Groundwater Samples Collected in Northern Iraq and a Summary of the Environmental Setting of the Adaya Burial Site

The Radiation Protection Center (RPC) of the Iraqi Ministry of Environment continues to evaluate the potential health impacts associated with the Adaya Burial Site, which is located 33 kilometers (20.5 miles) southwest of Mosul. This report documents the radiological analyses of 16 groundwater samples collected from wells located in the vicinity of the Adaya Burial Site and at other sites in northern Iraq. The Adaya Burial Site is a high-risk dump site because a large volume of radioactive material and contaminated soil is located on an unsecure hillside above the village of Tall ar Ragrag. The uranium activities for the 16 water samples in northern Iraq are considered to be naturally occurring and do not indicate artificial (man-made) contamination. With one exception, the alpha spectrometry results for the 16 wells that were sampled in 2019 indicate that the water quality concerning the three uranium isotopes (Uranium-233/234, Uranium-235/236, and Uranium-238) was acceptable for potable purposes (drinking and cooking). However, Well 7 in Mosul had a Uranium-233/234 activity concentration that slightly exceeded the World Health Organization guidance level. Eight of the 16 wells are located in the villages of Tall ar Ragrag and Adaya and had naturally occurring uranium concentrations. Wells in the villages of Tall ar Ragrag and Adaya are located near the Adaya Burial Site and should be sampled on an annual schedule. The list of groundwater analytes should include metals, total uranium, isotopic uranium, gross alpha/beta, gamma spectroscopy, organic compounds, and standard water quality parameters. Our current understanding of the hydrogeologic setting in the vicinity of the Adaya Burial Site is solely based on villager's domestic wells, topographic maps, and satellite imagery. To better understand the hydrogeologic setting, a Groundwater Monitoring Program needs to be developed and should include the installation of twelve groundwater monitoring wells in the vicinity of Tall ar Ragrag and the Adaya Burial Site. Characterization of the limestone aquifer and overlying alluvium is needed. RPC should continue to support health assessments for the villagers in Tall ar Ragrag and Adaya. Collecting samples for surface water (storm water), airborne dust, vegetation, and washway sediment should be conducted on a routine basis. Human access to the Adaya Burial Site needs to be strictly limited. Livestock access on or near the burial site needs to be eliminated. The surface-water exposure pathway is likely a greater threat than the groundwater exposure pathway. Installation of a surface-water diversion or collection system is recommended in order to reduce the potential for humans and livestock to come in contact with contaminated water and sediment. To reduce exposure to villagers, groundwater treatment should be considered if elevated uranium or other contaminants are detected in drinking water. Installing water-treatment systems would likely be quicker to accomplish than remediation and excavation of the Adaya Burial Site. The known potential for human exposure to uranium and metals (such as arsenic, chromium, selenium, and strontium) at the Adaya Burial Site is serious. Additional characterization , mitigation, and remediation efforts should be given a high priority.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Effects of Intrusive Confining Units on Groundwater Flow at Pahute Mesa, Nevada National Security Site

Geological, geophysical, hydrogeological, and hydrological evidence indicates that feeder dikes to rhyolitic lava flows may play a role as confining units in the groundwater flow system underlying the Pahute Mesa (PM) area of underground nuclear testing at the Nevada National Security Site (NNSS) (Figure 1). Of geological significance, the PM testing area is largely situated within the Silent Canyon caldera complex (SCCC), which includes both tuffs and lavas that were fed by (and are connected to) intrusive rock units such as feeder dikes prior to and during volcanic eruptions (Orkild et al., 1968; Byers et al., 1976a). Groundwater modeling results presented in this report support extension of “pool-dam” and “breach scenario” concepts presented in Jackson and Fenelon (2018) and Jackson et al. (2021) to low-permeability intrusive confining units (ICUs) such as feeder dikes that are expected to dissect and ring silicic calderas. As laterally extensive, steeply-dipping, and sub-parallel planar structures, ICUs may further compartmentalize groundwater flow in some areas of the volcanic aquifer system at PM. ICUs may explain unusually high horizontal hydraulic gradient southwest of Cheshire. Breaching of an ICU within which the Cheshire test was emplaced may explain changes in water levels before and after the Cheshire test was detonated in the U-20n emplacement hole on 2/14/76 (USDOE, 2015b). This analysis of the effects of ICUs on groundwater flow at PM addresses a need to evaluate uncertainty in the conceptual model or structural aspects of hydrologic source term (HST) or corrective action unit (CAU) models for groundwater flow and contaminant transport (FFACO, 1996, as amended). At Cheshire and possibly other test locations at PM, downgradient extent of radionuclide contamination may be overpredicted by not including the confining effects of ICUs. On the other hand, contaminant extent may be underpredicted by not considering transient flow effects and/or the re-equilibrated flow field caused by breaching of an ICU by an underground nuclear test. Near Cheshire, the effects of ICUs may be most apparent and impactful to the groundwater flow system at PM because of an unusual combination of thickness, extent, and frequency of dikes transecting the saturated CHLFA4 aquifer and feeding into the overlying unsaturated rhyolite lava of Windy Wash. Chapters 2, 3, and 4 present geological, geophysical, hydrogeological, and hydrological evidence for ICUs at PM. Chapters 5 and 6 examine effects of ICUs on groundwater flow including analysis by groundwater modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗

DOE OSTI2025

Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1

Abstract. Groundwater plays a key role in meeting water demands, supplying over 40 % of irrigation water globally, with this role likely to grow as water demands and surface water variability increase. A better understanding of the future role of groundwater in meeting sectoral demands requires an integrated hydro-economic evaluation of its cost and availability. Yet substantial gaps remain in our knowledge and modeling capabilities related to groundwater availability, recharge, feasible locations for extraction, extractable volumes, and associated extraction costs, which are essential for large-scale analyses of integrated human–water system scenarios, particularly at the global scale. To address these needs, we developed Superwell, a physics-based groundwater extraction and cost accounting model that operates at sub-annual temporal and at the coarsest 0.5° (≈50 km × 50 km) gridded spatial resolution with global coverage. The model produces location-specific groundwater supply–cost curves that provide the levelized cost to access different quantities of available groundwater. The inputs to Superwell include recent high-resolution hydrogeologic datasets of permeability, porosity, aquifer thickness, depth to water table, recharge, and hydrogeological complexity zones. It also accounts for well capital and maintenance costs, as well as the energy costs required to lift water to the surface. The model employs a Theis-based scheme coupled with an amortization-based cost accounting formulation to simulate groundwater extraction and quantify the cost of groundwater pumping. The result is a spatiotemporally flexible, physically realistic, economics-based model that produces groundwater supply–cost curves. We show examples of these supply–cost curves and the insights that can be derived from them across a set of scenarios designed to explore model outcomes. The supply–cost curves produced by the model show that most (90 %) nonrenewable groundwater in storage globally is extractable at costs lower than USD 0.57 m−3, while half of the volume remains extractable at under USD 0.108 m−3. The global unit cost is estimated to range from a minimum of USD 0.004 m−3 to a maximum of USD 3.971 m−3. We also demonstrate and discuss examples of how these cost curves could be used by linking Superwell's outputs with other models to explore coupled human–environmental system challenges, such as water resources planning and management, or broader analyses of multisectoral feedbacks.

Global Change Analysis Model (GCAM)↗

DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗

NASA NTRS1998

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗

NASA NTRS2023

Large-Scale Groundwater Monitoring in Brazil Assisted With Satellite-Based Artificial Intelligence Techniques

Here, we develop and test an artificial intelligence (AI)-based approach to monitor major Brazilian aquifers. The approach combines Gravity Recovery and Climate Experiment (GRACE) data and ground-based hydrogeological measurements from Brazil’s Integrated Groundwater Monitoring Network at hundreds of wells distributed in twelve aquifers across the country. We tested model ensembles based on three AI approaches: Extreme Gradient Boost, Light Gradient Boosting Model and CatBoost, followed by a Linear Regression (LR) step. The approach is further boosted with wavelet and seasonal decomposition processes applied to GRACE data. To determine the AI-based model’s sensitivity to data availability, we propose four experiments combining hydrogeological measurements from different aquifers. Groundwater storage estimates from the Global Land Data Assimilation System (GLDAS) are used as benchmark. A sensitivity analysis shows that the LR-based model ensemble is the best suited and to reproduce groundwater storage change in all studied Brazilian aquifers. Results show that the proposed approach outperforms GLDAS in all experiments, with an RMSE value of 2.68cm for the experiment that covers all monitored wells in Brazil. GLDAS resulted in RMSE=6.76cm. Using our AI model outputs, we quantified the groundwater storage change of two major aquifers, Urucuia and Bauru-Caiuá, over the past two decades: -31km 3 and -6km 3 , respectively. Water loss is driven by a prolonged drought across most of the country and intensification of groundwater pumping for irrigation. This study demonstrates that combining satellite data and AI can be a cost-effective alternative to monitor poorly equipped aquifers at the continental scale, with possible global replicability.

GRACE↗

DOE OSTI2022

Machine Learning-Assisted High-Temperature Reservoir Thermal Energy Storage Optimization: Numerical Modeling and Machine Learning Input and Output Files

This data set includes the numerical modeling input files and output files used to synthesize data, and the reduced-order machine learning models trained from the synthesized data for reservoir thermal energy storage site identification. In this study, a machine-learning-assisted computational framework is presented to identify High-Temperature Reservoir Thermal Energy Storage (HT-RTES) site with optimal performance metrics by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems. All the simulations and the neural network model were done by Idaho National Laboratory. A detailed description of the work was reported in publication linked below.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Community Geothermal: Final Thermal Conductivity Test Report and Data Logs - Carbondale, CO

Provided here are a final thermal conductivity test report, a drilling log, and a heat rejection log from Carbondale, CO. Also attached is a report made before drilling, which contains predictions on the hydrogeologic conditions of the drill site. The forty-eight (48.9) hour in-situ thermal conductivity test was performed on the 28th through the 30th, of November 2023. The test was performed at the borehole drilled on November 16th through the 20th, at the 3rd Street Center at 520 S. 3rd, Street in Carbondale, Colorado. Testing was done with a certified Ewbank portable test unit.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Inverse Analysis with Variational Autoencoders: A Comparison of Shallow and Deep Networks

Inverse problems are applied to determine unknown properties by matching observational data with a physical model that often takes many parameters as input. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented involving regularization with a technique known as a variational autoencoder (VAE), which is trained to map a high-dimensional parameter space with a complex structure to a low-dimensional latent space with a simple structure. We apply this approach to unconditioned realizations of the parameters (heterogeneous hydraulic fields) for a hydrogeological inverse problem. Two types of hydraulic conductivity fields are used to evaluate the characterization for the different levels of heterogeneity complexity of the physical inputs. This approach keeps the computational cost of generating the training data low. The reason is unconditioned realizations neither rely on the observational data used to perform the inverse analysis nor require any groundwater flow model forward runs. In addition, this approach applies regularization on a low-dimensional latent space from the VAE and increases optimization efficiency through automatic differentiation. Furthermore, two different neural network (NN) structures are tested for their utility in using the VAE for inverse analysis. The performance of a deep, convolutional neural network strongly depends on the dimensionality of the latent space, which requires tuning. In contrast, a shallow, dense neural network provides consistently accurate characterization without tuning. Furthermore, our approach evaluates the advantages of the shallow, dense neural network over the deep, convolutional one and enables future application to a wide range of inverse problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Review: “Jacob’s Zoo”—how using Jacob’s method for aquifer testing leads to more intuitive understanding of aquifer characteristics

The interpretation of aquifer responses to pumping tests is an important tool for assessing aquifer geometry and properties, which are critical in the assessment of water resources or in environmental remediation. However, the responses of aquifers, measured by time-drawdown relationships in monitoring wells, are nonunique solutions that are affected by many factors. Jacob's Zoo is a collection of graphical interpretations that allows students and practitioners to develop an intuitive feel for how natural hydrogeological systems work, and develop a set of skills that provide a better understanding of aquifer properties much beyond interpretation of pumping tests. Jacob's Zoo, based on the work of Jacob (1950), fosters a deeper understanding, although few practitioners realize the full utility of the method. Jacob CE (1950) Flow of groundwater, In: Rouse H (ed) Engineering Hydraulics, Wiley, New York. P 321-386.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Using groundwater noble gas measurements to confirm paleorecharge hypotheses at Pahute Mesa, Nevada National Security Site, Nevada, USA

Pahute Mesa (Nevada, USA) was the site of 85 underground nuclear tests between 1965 and 1992 whose residual radio-chemical inventory poses a contaminant threat to local groundwater resources. To assess the fate and transport of these radionuclides in groundwater, several geochemical studies, including this one, have been conducted to understand long-term groundwater flow patterns at, and downgradient of, Pahute Mesa. Groundwater samples from the Pahute Mesa area at the Nevada National Security Site were analyzed for noble gases ( 3 He/ 4 He isotopic ratio and He, Ne, Kr, Ar and Xe concentrations) to infer groundwater recharge temperatures and provide additional evidence in support of groundwater flow patterns and the presence of late-pluvial Pleistocene groundwater. These data confirm recent findings (based on groundwater geochemistry and isotope compositions) that most recharge in the area occurred under temperature conditions much cooler than are found today, except near major drainages where ephemeral run of continues to recharge shallow groundwater. These noble gas data indicate that noble gas recharge temperatures and terrigenic He components have been reset where groundwater has flowed upward over barriers formed by caldera margins, or near the regional discharge area (Oasis Valley) where shallow groundwater has interacted with shallow soil gas. This evidence from groundwater noble gas measurements, in combination with other isotopic and geochemical studies, supports the current understanding of the unusual hydrogeology and recharge history of Pahute Mesa.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Factors influencing water quality in surface water and alluvial groundwaters downgradient of a reclaimed surface coal mine in the Powder River Basin of southeastern Montana, USA

Abstract Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Keeshin, Skye I. (ORCID:0009000986449872)↗

DOE OSTI2024

Temperature effects on adsorption and capillarity water retention mechanisms in constrained unsaturated soils

Abstract This paper focuses on the impact of elevated temperatures on the adsorptive and capillarity water retention mechanisms of unsaturated soils under constrained (constant volume) conditions. This topic is critical for simulating the thermo-hydraulic behavior of soils in hydrogeological or geotechnical applications, including climate change effects on near surface soils, energy piles or soil borehole thermal energy storage systems in unsaturated soil layers, and buffers for geological nuclear waste repositories. A nonisothermal soil water retention curve (SWRC) that separately considers the temperature-dependency of the key parameters governing adsorptive and capillarity water retention mechanisms and soil physical parameters (e.g., surface tension, contact angle, adsorption capacity, cation exchange capacity, mean cavitation suction, air entry value and equilibrium film thickness) was developed to provide insights into the impact of temperature on water retention over the full suction range. The nonisothermal SWRC was validated using experimental data on high plasticity clays, with a good prediction of temperature effects on adsorption and capillarity water retention mechanisms in constrained unsaturated soils.

Lu, Yu↗