Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗
DOE OSTI2024

Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes

Abstract Mobile genetic elements (MGEs) like plasmids, viruses, and transposable elements can provide fitness benefits to their hosts for survival in the presence of environmental stressors. Heavy metal resistance genes (HMRGs) are frequently observed on MGEs, suggesting that MGEs may be an important driver of adaptive evolution in environments contaminated with heavy metals. Here, we report the meta-mobilome of the heavy metal-contaminated regions of the Oak Ridge Reservation subsurface. This meta-mobilome was compared with one derived from samples collected from unimpacted regions of the Oak Ridge Reservation subsurface. We assembled 1615 unique circularized DNA elements that we propose to be MGEs. The circular elements from the highly contaminated subsurface were enriched in HMRG clusters relative to those from the nearby unimpacted regions. Additionally, we found that these HMRGs were associated with Gamma and Betaproteobacteria hosts in the contaminated subsurface and potentially facilitate the persistence and dominance of these taxa in this region. Finally, the HMRGs were associated with conjugative elements, suggesting their potential for future lateral transfer. We demonstrate how our understanding of MGE ecology, evolution, and function can be enhanced through the genomic context provided by completed MGE assemblies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Future Circular Collider Feasibility Study Report

Volume 3 of the FCC Feasibility Report presents studies related to civil engineering, the development of a project implementation scenario, and environmental and sustainability aspects. The report details the iterative improvements made to the civil engineering concepts since 2018, taking into account subsurface conditions, accelerator and experiment requirements, and territorial considerations. It outlines a technically feasible and economically viable civil engineering configuration that serves as the baseline for detailed subsurface investigations, construction design, cost estimation, and project implementation planning. Additionally, the report highlights ongoing subsurface investigations in key areas to support the development of an improved 3D subsurface model of the region. The report describes the development of the project scenario based on the ‘avoid-reduce-compensate’ iterative optimisation approach. The reference scenario balances optimal physics performance with territorial compatibility, implementation risks, and costs. Environmental field investigations covering almost 600 hectares of terrain—including numerous urban, economic, social, and technical aspects—confirmed the project’s technical feasibility and contributed to the preparation of essential input documents for the formal project authorisation phase. The summary also highlights the initiation of public dialogue as part of the authorisation process. The results of a comprehensive socio-economic impact assessment, which included significant environmental effects, are presented. Even under the most conservative and stringent conditions, a positive benefit-cost ratio for the FCC-ee is obtained. Finally, the report provides a summary of the studies conducted to document the current state of the environment.

Benedikt, M. [European Organization for Nuclear Re↗
DOE OSTI2025

Improving Stability of an Optical Fiber pH Sensor with a Calcined Polyethylenimine-Coating at High Pressures and Temperatures

With an increased interest in subsurface gas storage technology for various energy applications, monitoring wellbore structural stability and subsurface geochemistry has become more pressing, and pH is a key parameter to measure. As high pressure and elevated temperatures in subsurface conditions are comparatively harsh relative to that expected for most standard pH sensor designs, any pH monitoring hardware must be designed for extended exposure to high pressures and temperatures. We previously reported that an optical fiber pH sensor functionalized with a calcined polyethylenimine coating had shown some promise as a high temperature and pressure pH sensor but with some drifting when operating for longer than 8 hours. In this paper, we investigated the coating composition and potential cause of the drifting and improved the stability of the prepared coating to minimize sensor drift under simulated wellbore conditions. Scanning electron microscopy (SEM) had previously shown moderate cracking at high pressures over short tests. By applying X-ray photoelectron spectroscopy (XPS) to characterize the sensor coating before and after one week of testing in an H2/CH4 gas blend at 80 °C and 900 psi, a compositional change in the coating was observable, which may indicate susceptibility to alteration by subsurface gas storage conditions. Non-reducing (CH4, N2) environments were also tested, and confirmed that both temperature and pressure were also contributing to the drift.

energy infrastructure↗
DOE OSTI2026

Illinois Storage Corridor - CarbonSAFE Phase III: Policy, Regulatory, Legal and Permitting Characteristics; Subtask 5.5

The Illinois Storage Corridor (ISC) project aims to advance the commercial development and implementation of carbon capture and storage (CCS) technologies within a region in Illinois of suitable geology for carbon dioxide (CO 2 ) storage in deep saline aquifers. The project partners have annual emissions greater than 6.5 million tonnes per year; storage hubs are being explored at two sites—one near the One Earth Energy facility in east-central Illinois, and one near the Prairie State Generating Company campus in southwest Illinois. To implement the technology, legal, policy, and economic considerations must be explored. The United States Environmental Protection Agency (US EPA) administers the Underground Injection Control (UIC) program and is responsible for issuing Class VI permits to construct and operate CO 2 injection wells, i.e., a UIC Class VI well permit is required to inject CO 2 into the subsurface for geologic sequestration. The largest consideration of the Class VI well requirements is to protect underground sources of drinking water (USDWs). Additionally discussed are considerations for Class VI permits relating to public engagement. Besides permitting, property rights to storage sites, subsurface pore spaces, and areas for pipeline transportation must be secured. Pore space rights is still a novel concept being explored and not yet addressed by the Illinois legislature. It is believed that surface property rights are required for the subsurface pore space below, so long as there is not a separated mineral estate in the subsurface. Illinois legislature has addressed securing rights-of-way for CO 2 transportation, allowing for easements and the exercise of eminent domain to secure such rights. Economically, the incentives for CCS are ever-expanding. Recently, the federal government broadened the availability and increased the dollar-amounts for the § 45Q tax credits for geological storage of carbon oxides. Congress has also authorized the Advanced Industrial Facilities Development Program which allocates billions of dollars in funds to installing technology at industrial facilities to reduce greenhouse gas emissions. On the state level, Illinois has had its own incentives for CCS projects since 2009 and expanded its emission-related goals again in 2021. The law enacted in 2021 specifically creates a commission to explore implementing CCS at Prairie State Generating Company, a partner on the Illinois Storage Corridor project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗
DOE OSTI2025

High-Temperature Aquifer Thermal Energy Storage (HT-ATES) Projects in Germany and the Netherlands—Review and Lessons Learned

Aquifer thermal energy storage (ATES) is a concept that can help to address heating and cooling needs through the use of the subsurface as a seasonal thermal energy storage (STES) system. Over 2800 ATES systems have been deployed with storage temperatures typically below 25 °C and only a few with higher temperatures (>40 °C), which would increase the energy density and utility of the stored thermal fluids. Until now, only a few high-temperature aquifer thermal energy storage (HT-ATES) projects have been initiated and are still in operation. These HT-ATES projects have encountered a range of technical and non-technical challenges. This study reviews ten such projects: four in Germany and six in the Netherlands. The non-technical issues include public acceptance, a lack of regulatory framework for these systems, managing overlapping uses of the subsurface, managing changes with the providers and off-takers of thermal energy, and obtaining financing to implement these projects. Common technical issues include geological factors such as incomplete characterization of the subsurface and reservoir heterogeneity; geochemical issues such as mineral scaling, corrosion, and biofouling; lower than expected thermal recovery; and issues with system design and reliability. This review highlights benefits and challenges faced by HT-ATES projects with the goal to use the lessons learned to improve the siting, design, development, and operation of such systems. Recommendations include improved initial subsurface site characterization, use of coupled process models to optimize system design and predict system performance, cascaded uses of stored thermal energy to better utilize the stored heat, monitoring networks to provide feedback on system performance, and expanded system scale to allow for continued operation even when maintenance of some system components is required. Techno-economic modeling and risk analysis could be used to optimize such HT-ATES project design and identify key factors that will affect sustained economic viability. In addition, design flexibility is important for these systems to allow for changing conditions regarding the supply and demand of thermal energy. Adopting these findings should improve the performance and reduce the risks for future HT-ATES projects worldwide.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗
DOE OSTI2025

Improving Stability of an Optical Fiber pH Sensor with a Calcined Polyethylenimine-Coating at High Pressures and Temperatures

With an increased interest in subsurface gas storage technology for various energy applications, monitoring wellbore structural stability and subsurface geochemistry has become more pressing, and pH is a key parameter to measure. As high pressure and elevated temperatures in subsurface conditions are comparatively harsh relative to that expected for most standard pH sensor designs, any pH monitoring hardware must be designed for extended exposure to high pressures and temperatures. We previously reported that an optical fiber pH sensor functionalized with a calcined polyethylenimine coating had shown some promise as a high temperature and pressure pH sensor but with some drifting when operating for longer than 8 hours. In this paper, we investigated the coating composition and potential cause of the drifting and improved the stability of the prepared coating to minimize sensor drift under simulated wellbore conditions. Scanning electron microscopy (SEM) had previously shown moderate cracking at high pressures over short tests. By applying X-ray photoelectron spectroscopy (XPS) to characterize the sensor coating before and after one week of testing in an H2/CH4 gas blend at 80°C and 900 psi, a compositional change in the coating was observable, which may indicate susceptibility to alteration by subsurface gas storage conditions. Non-reducing (CH4, N2) environments were also tested, and confirmed that both temperature and pressure were also contributing to the drift.

energy infrastructure↗
DOE OSTI2024

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Oxide nucleation via threshold-driven volume instability in Ni oxidation

Oxide nucleation dictates the onset of metal oxidation, yet the atomistic pathways bridging initial oxygen adsorption to bulk phase transformation remain elusive. Here, we utilize in-situ atomic-resolution imaging to directly capture the nucleation of nickel oxide, revealing a threshold-driven, cooperative transformation of Ni into NiO. Our observations identify a distinct incubation period during which oxygen progressively accumulates in subsurface layers. Once a critical concentration and penetration depth are reached, a collective lattice reconfiguration is triggered, abruptly converting multiple Ni layers into NiO. Atomistic simulations corroborate this mechanism, identifying a cooperative lattice instability induced by subsurface oxygen saturation. These results establish subsurface oxygen incorporation as the missing mechanistic link between surface adsorption and three-dimensional oxide formation, providing an atomistic framework to understand and control reactive phase transformations in materials synthesis, catalysis, and degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2026

Data & Code from Phoenix CPPP Phase 2 Analysis

This data and code package supports the analysis presented in “Beyond Surface Cooling: Comprehensive Field Assessment of Reflective Pavement Thermal Performance in Phoenix, Arizona” and provides fully reproducible workflows for evaluating the thermal performance of cool pavement treatments in a hot urban environment. The dataset integrates multi-modal field measurements collected across residential and nonresidential settings, including mobile air temperature traverses, stationary air temperature monitoring, residential mean radiant temperature (MRT) measurements, subsurface temperature profiles, and controlled testbed observations. The data package contains raw and processed datasets in comma-separated value (CSV) format, accompanying metadata files describing site characteristics and measurement protocols, and R scripts (.R files) used for data cleaning, time synchronization, spatial and temporal matching, quality control filtering, statistical comparison, and figure generation. All analyses were conducted using R (version ≥ 4.2.0) with commonly available packages (e.g., tidyverse, lubridate, data.table, ggplot2). No proprietary software is required to reproduce results. Field campaigns were designed to quantify the effects of high-reflectance pavement coatings on surface temperature, near-surface air temperature, subsurface heat propagation, and radiative heat exposure. Temporal alignment procedures include standardized timestamp conversion and nearest-neighbor matching of high-frequency sensor measurements to stop-based metadata within defined tolerance windows to ensure comparability across instruments. The workflows generate summary statistics, treatment–control contrasts, depth-dependent thermal gradients, and time-series visualizations used in the associated publication. By integrating mobile, stationary, radiative, and subsurface measurements within a unified and transparent processing framework, this package enables comprehensive evaluation of cool pavement performance across multiple thermal exposure pathways and supports reuse in future urban heat mitigation and climate resilience studies.

AIR TEMPERATURE↗
DOE OSTI2026

Electromagnetic Induction (EMI) Data, 2024, Trail Creek, Colorado

This dataset contains Electromagnetic Induction (EMI) data collected at Trail Creek, Colorado, in 2024. EMI surveys were conducted to investigate the spatial distribution of electrical conductivity in the subsurface, providing insights into soil moisture and subsurface geological features. The surveys were performed along multiple transects to capture variations in conductivity influenced by changes in soil composition, moisture content, and underlying geological structures. This dataset complements other geophysical data collected in the region, including Electrical Resistivity Tomography (ERT) and Terrestrial LiDAR Scanning (TLS), providing a detailed understanding of the subsurface and its impact on surface vegetation and hydrological processes. The data are valuable for environmental geophysics, ecological research, and hydrological modeling in mountainous ecosystems. The files include: - data.zip: the raw EMI data (.csv) - inversion.zip: the inverted resistivity model (.csv and .kml) - kriging.zip: the kriging resistivity model (.csv, .tif, .kmz) - flmd.csv: file level metadata file describing all files within this dataset - dd.csv: data dictionary file describing the column headers within CSV files This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

CMD Mini-Explorer↗
DOE OSTI2025

High-Resolution Computed Tomography Dataset of Mount Simon Sandstone

The Illinois Basin is a critical structure for subsurface energy related activities and their implementation in the United States. The Mount Simon Sandstone has been identified as a storage target for permanent and transient storage of fluids in the basin. Known for its exceptional thickness, depth, porosity, and sealing properties of overlying formations, this saline reservoir is crucial for long-term subsurface energy efforts. We present an extensive Computed Tomography (CT) dataset on a high porosity and permeability zone in the lower Mount Simon Sandstone available on the Energy Data eXchange® (EDX). This publicly accessible database comprises over 500 GB of high-resolution CT scans of six core samples, with resolutions ranging from 14.8 µm to 0.7 µm per pixel. The scans include both dry sandstone samples and those saturated with multiple fluids, allowing for comparative analyses across different conditions and resolutions. Coarser scans capture the bedding structure of the sandstone, while finer resolutions reveal detailed pore infill and throat characteristics. Metadata on location, depth, and saturation state enhance usability, enabling quick identification and cross-sample comparisons. By providing a robust resource for research and collaboration, the database contributes to domestic energy advancement by supporting continued progress in the use of the subsurface for energy solutions.

characterization↗
DOE OSTI2024

Anomalous Role of Carbon in Pd‐Catalyzed Selective Hydrogenation

Carbonaceous species, including subsurface carbidic carbon and surface carbon, play crucial roles in heterogeneous catalysis. Many reports suggested the importance of subsurface carbon in the selective hydrogenation of alkynes over Pd‐based catalysts. However, the role of surface carbon has been largely overlooked. Here, we demonstrate that subsurface carbon in Pd is not responsible for the selectivity in acetylene hydrogenation. In contrast, the structure of surface carbonaceous species plays a decisive role in hydrogenation selectivity. Electron microscopy and spectroscopy evidence, along with theoretical modelling, reveal that partial graphitization of surface carbonaceous species results in unique spatial confinement of surface reaction intermediates, thus altering the reaction energy landscape in favour of ethylene desorption as opposed to over‐hydrogenation. This mechanism for selectivity control is analogous to enzyme catalysis, where the active centers selectively attract reactants and release products. Similar mechanism may be present in CO/CO 2 hydrogenation and alkane dehydrogenation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Explosive Byproduct Gas Transport Through Sorptive Geomedia

Current underground nuclear explosion (UNE) detection strategies rely heavily on atmospheric noble gas sampling of radioxenon. However, discriminating nuclear weapons testing programs from civilian sources is difficult due to highly variable atmospheric radioxenon backgrounds and processes affecting subsurface transport of parent radionuclides. Here, we aim to study the transport of gases produced by subsurface explosions as novel stable signatures for underground nuclear explosion (UNE) monitoring. These gases may be produced in large quantities with distinct molecular ratios, which will be impacted by subsurface transport processes. To demonstrate how ratios of gases produced by explosions can change during transport in geomaterials, we conducted laboratory benchtop experiments on the transport of carbon dioxide (CO 2 ) and hydrogen (H 2 ) gases through variably saturated zeolitic tuff, which is abundant at the historic US testing site. We observed that zeolitic tuff sorbs substantial quantities of CO 2 while allowing H 2 to transport more freely, leading to changes in the molecular ratios of the two gases along the transport pathway. Gas uptake in the dry zeolitic tuff core was 72.3% for CO 2 , compared with 53.4% for xenon and 7.6% for H 2 . The presence of 20% water saturation disrupted the CO 2 sorption process, though to a lesser extent than observed for noble gases, with a 36.7% drop in xenon sorption compared with a 21.9% drop for CO 2 . These results represent the first observations of zeolite sorption altering explosive gas ratios during transport through geomedia relevant to nuclear proliferation monitoring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗