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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Mars Sulfate Formation Sourced in Sulfide-Enriched Subsurface Fluids: The Rio Tinto Model

The extensive evidence for sulfate deposits on Mars provided by analyses of MER and Mars Express data shows that the sulfur played an essential role in the geochemical cycles of the planet, including reservoirs in the atmosphere, hydro-sphere and geosphere. Overall the data are consistent with a fluvial/lacustrine-evaporative origin of at least some of the sulfate deposits, with mineral precipitation through oversaturation of salty acidic fluids enriched in sulfates. This scenario requires reservoirs of sulfur and associated cations, as well as an acidic and oxidizing hydrochemistry which could be provided by surface and subsurface catching of meteoric waters resulting in the presence of sulfur-bearing gases and steam photochemistry. In this work we suggest a new scenario for the extensive generation of sulfates in Mars based on the observation of seasonal changes in the redox and pH of subsurface waters enriched in sulfur that supply the acidic Mars process analog of Rio Tinto. This model considers the long-term subsurface storage of sulfur during most of Noachian and its release from the late Noachian to Hesperian time through weathering by meteoric fluids that would acidify and oxidize the sulfur bearing compounds stored in the subsurface.

Fernandez-Remolar, D. C.↗
NASA NTRS2010

Drilling on the Moon and Mars: Developing the Science Approach for Subsurface Exploration with Human Crews

DOMEX (Drilling on the Moon and Mars in Human Exploration) is using analog missions to develop the approach for using human crews to perform science activities on the Moon and Mars involving exploration and sampling of the subsurface. Subsurface science is an important activity that may be uniquely enabled by human crews. DOMEX provides an opportunity to plan and execute planetary mission science activities without the expense and overhead of a planetary mission. Objectives: The objective of this first in a series of DOMEX missions were to 1) explore the regional area to understand the geologic context and determine stratigraphy and geologic history of various geologic units in the area. 2) Explore for and characterize sites for deploying a deep (10 m depth) drilling system in a subsequent field season. 3) Perform GPR on candidate drill sites. 4) Select sites that represent different geological units deposited in different epochs and collect soil cores using sterile procedures for mineralogical, organic and biological analysis. 5) Operate the MUM in 3 different sites representing different geological units and soil characteristics. 6) Collect rock and soil samples of sites visited and analyze them at the habitat. Results: At mission start the crew performed a regional survey to identify major geologic units that were correlated to recognized stratigraphy and regional geologic maps. Several candidate drill sites were identified. During the rest of the mission, successful GPR surveys were conducted in four locations. Soil cores were collected in 5 locations representing soils from 4 different geologic units, to depths up to 1m. Soil cores from two locations were analyzed with PCR in the laboratory. The remainder were reserved for subsequent analysis. XRD analysis was performed in the habitat and in the field on 39 samples, to assist with sample characterization, conservation, and archiving. MUM was deployed at 3 field locations and 1 test location (outside the habitat) where it operated autonomously for 2-4 hours at each site. Depths achieved ranged from 15 to 70 cm depending on the soil compressive strength and the presence and depth of subsurface indurated layers. Subsurface samples weighing 0.5 to 1 g were collected at the deepest depth encountered at each of the sites using the MUM automated sample collection system, and subsequently analyzed with XRD. Downhole inspection of holes produced by MUM with the Raman spectrometer was acquired on two of the holes and spectral features associated with selenite were identified in specific soil layers. Previously unreported fossilized remains of vertebrate fauna from the Jurassic era were discovered during our mission. Analysis of mineral biomarkers associated with this discovery are underway.

Stoker, C. R.↗
DOE OSTI2018

Distributed Fiber Optic Arrays: Integrated Temperature and Seismic Sensing for Detection of CO 2 Flow, Leakage and Subsurface Distribution (Final Report)

This report describes the seismic and heat-pulse monitoring equipment, surveys, results and accomplishments of the fiber optic (FO) project funded by the U.S. Department of Energy (DOE) under Agreement Number DE-FE0012700. The formal title of the project is “Distributed Fiber Optic Arrays: Integrated Temperature and Seismic Sensing for Detection of CO 2 Flow, Leakage and Subsurface Distribution.” The project focuses on the field acquisition of FO data using Distributed Acoustic Sensing (DAS) arrays for subsurface imaging of geologic structure and to track the CO 2 migration in the injection interval using time-lapse techniques. A second component of the project uses a novel heat-pulse monitoring method that relies on Distributed Temperature Sensing (DTS) to evaluate near-field leakage of CO 2 or brine in and around the wellbore. The Southeast Regional Carbon Sequestration (SECARB) partnership hosted the project at the Anthropogenic Test site near Citronelle, Alabama after the original site host, the Kansas Geological Survey (KGS), experienced significant delays in permitting its CO 2 injection well at the Wellington Oil Field site in Kansas. The KGS project was later terminated in 2017 before an injection permit issued. A second field trial was planned for the Livingston Oil Field, located in Louisiana, but the site host ran into financial troubles and the project was terminated in 2015. A replacement field-site host was not found for the Livingston field trial. Both the Wellington, Kansas and Livingston, Louisiana projects were funded by DOE. A novel sensor platform referred to as the Modular Borehole Monitoring (MBM) system deployed at the Citronelle field site was utilized by the FO project. Two different seismic sensor technologies were incorporated in the Citronelle MBM system design including: (1) a semi-permanent, tubing-deployed, 18-level geophone array with custom hydraulic clamps; and (2) two, single mode optical fibers for Distributed Acoustic Sensing (DAS). Incorporation of the 18-level geophone array in the MBM system provided an ideal opportunity to compare the DAS response to conventional geophones. A heat-pulse monitoring system consisting of two, multi-mode optical fibers for Distributed Temperature Sensing (DTS) and copper heater element were also incorporated into the MBM package. This report provides an overview of the design of the MBM, DAS, and DTS equipment and the surveys performed during the study. The heat-pulse monitoring technique was used to diagnose a completion problem with observation well D-9-8 #2 where the MBM system was installed in 2012, thus successfully demonstrating its application for leak detection. In 2015, CO 2 unexpectedly arrived at the observation well approximately one year after CO 2 injection and heat-pulse monitoring ended. It was determined at this time that saltwater intrusion had likely penetrated the stainless-steel capillary tube containing the fibers and copper heater elements, thus preventing the heater from being used. DTS measurements indicated that the fiber had also degraded probably due to hydrogen darkening. It is recommended that the capillary tube be constructed of a more corrosive resistant material (e.g., Inconel) in the future should it be deployed in similar harsh downhole environments. Geophysical imaging of the subsurface and monitoring of the CO 2 were successful using DAS in the vertical seismic profile (VSP) survey configuration. A high resolution VSP image of the subsurface was obtained in 2014 using DAS, which exceeded our expectations in comparison to a lower resolution image obtained by SECARB in 2012 using an 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained using two DAS VSP surveys from June 2014 and December 2015, thus successfully demonstrating its application. DAS data were also acquired during a crosswell seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the crosswell survey configuration because the noise level was too high, increasing by a factor of ten, over the high frequency output (100-1,200) of the piezoelectric source.

58 GEOSCIENCES↗
DOE OSTI2018

Randomization in Characterizing the Subsurface

Current methods for characterizing Earth’s subsurface, such as standard inverse techniques, are not sufficiently accurate to meet the needs of modern applications in the fields of energy exploration, environmental management, and global security. While increasing the quantity of field measurements and robustness of the applied data-/model-analysis methods can improve accuracy, such approaches can be computationally impractical for large data sets and complex site conditions. Thus, there is a need to develop economically-feasible and robust computational methods while maintaining accuracy. For example, in-field drilling for geothermal operations may yield high failure rates, resulting in unacceptably high costs; errors and/or large uncertainties in the estimated subsurface characteristics are the main impediment to the successful siting of an in-field well. This problem is not uniquely geothermal. Accurate characterization of uncertain subsurface properties is also essential for monitoring storage of carbon dioxide, estimating pathways of subsurface contaminant transport, and supervising ground-based nuclear-explosion tests.

58 GEOSCIENCES↗
DOE OSTI2024

JAMES BUTTLE REVIEW: Interflow, subsurface stormflow and throughflow: A synthesis of field work and modelling

Interflow, throughflow and subsurface stormflow are interchangeable terms that refer to the lateral subsurface flow above a restricting layer of lower hydraulic conductivity that occurs during and following storm events. Interflow (used here) is a more dominant process in steeper catchments with high infiltration capacity soils overlying a more impermeable soil or geologic layer. Interflow as a runoff process was first recognised in the early 1900s, yet hydrologists still struggle to predict its occurrence, persistence, importance, interaction with other streamflow generation processes, and potential to connect to valleys and streams during and following storms. We review the history of interflow research and address some of the challenges in understanding its role in runoff production. We argue that characterising the controls on interflow initiation and occurrence relies on detailed field observations of subsurface properties, which exist only in limited experimental settings. This data shortcoming contributes to our inability to predict interflow or determine its contribution to streamflow more broadly. There remain many opportunities to advance our understanding of interflow that include both modelling and experimental or observational approaches in hydrology.

hillslope hydrology↗
DOE OSTI2022

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗
DOE OSTI2018

Characterizing Technetium in Subsurface Sediments for Contaminant Remediation

Technetium-99 (Tc) contamination remains a major environmental problem at legacy nuclear reprocessing sites, including the Hanford Site (Washington State, USA) where ~700 Ci of Tc has been released into the subsurface. The environmental risk associated with Tc is due to its long half-life (213,000 years) and high mobility of the oxidized anionic pertechnetate species, Tc(VII)O 4 - , which warranted the US Department of Energy to identify Tc as a primary contaminant of concern. This risk generates a critical need to resolve scientific and technical issues underpinning the ability to safely, cost-effectively, and efficiently remediate Tc-contaminated vadose zone sediments and groundwater and mitigate risk to human health and the environment. Developing enhanced attenuation and efficient remediation strategies for released Tc requires a complete understanding of retardation processes and Tc mass flux, including the different mechanisms by which Tc is immobilized in the subsurface and the effect of localized subsurface conditions. Selection of over 30 sediments from Hanford waste disposal sites, based on historical information and sediment characterization, for analysis by autoradiography revealed that Tc concentrations were generally below the detection limit of 5 mg Tc/g sediment. When Tc was measurable in vadose zone sediments, it was predominantly present as TcO 4 - in water films associated with fine-grained sediments, with a maximum of 12 % of the total Tc present in the acid-extractable fraction, defined here as the immobile fraction. However, beneath one waste disposal site, where sediments containing minerals with reducing capacity intercepted miscellaneous fission product recovery waste and waste from the bismuth phosphate process, the amount of Tc present in the immobile fraction was 53 % of the total. Characterization of Tc-containing phases present in these field-contaminated sediments for the first time using Tc K-edge X-ray absorption near edge structure spectroscopy revealed that, as well as Tc present as Tc(VII)O 4 - in pore water associated with fine-grained sediments, Tc was also: (i) physically encapsulated within solid phases precipitated from other waste components and in multi-component phosphate minerals; and (ii) present as mixed Tc(VII)/Tc(IV) in localized reducing zones. These results will be used to develop improved long-term Tc remediation strategies optimized for field application, through stimulation of conditions that enhance Tc attenuation.

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

Consensus Equilibrium for Subsurface Delineation

Heterogeneity and insufficient site characterization limit our knowledge of the subsurface. Inversion techniques, which minimize the mismatch between observations and model predictions, have become an essential tool of subsurface characterization. Most optimization-based approaches fail to incorporate various implicit priors and capture the geological complexity. We overcome these limitations by deploying the plug-and-play and consensus equilibrium (CE) strategies, which provide a flexible framework for image reconstruction. Our CE methodology for spatial delineation of geologic formations consists of an image denoiser and a variational auto-encoder (deep learning-based emulator). The former ameliorates the reconstruction noise, yielding well-defined geological structures; its mathematical equivalence with the proximal operator allows the deployment of advanced denoisers (e.g., CNN-based denoiser) that do not correspond to a regularization objective. The latter defines a geology prior that imposes a geological constraint, for example, continuity and shape of geological features, onto the reconstructed image. Here, we conduct a series of numerical experiments dealing with transient two-dimensional flow driven by a pumping well and natural hydraulic head gradient. They demonstrate the CE framework's ability to delineate, both probabilistically and deterministically, complex subsurface environments with sufficient quality.

42 ENGINEERING↗
DOE OSTI2022

Geological activity shapes the microbiome in deep-subsurface aquifers by advection

Subsurface environments host diverse microorganisms in fluid-filled fractures; however, little is known about how geological and hydrological processes shape the subterranean biosphere. Here, we sampled three flowing boreholes weekly for 10 mo in a 1478-m-deep fractured rock aquifer to study the role of fracture activity (defined as seismically or aseismically induced fracture aperture change) and advection on fluid-associated microbial community composition. We found that despite a largely stable deep-subsurface fluid microbiome, drastic community-level shifts occurred after events signifying physical changes in the permeable fracture network. The community-level shifts include the emergence of microbial families from undetected to over 50% relative abundance, as well as the replacement of the community in one borehole by the earlier community from a different borehole. Null-model analysis indicates that the observed spatial and temporal community turnover was primarily driven by stochastic processes (as opposed to deterministic processes). We, therefore, conclude that the observed community-level shifts resulted from the physical transport of distinct microbial communities from other fracture(s) that outpaced environmental selection. Given that geological activity is a major cause of fracture activity and that geological activity is ubiquitous across space and time on Earth, our findings suggest that advection induced by geological activity is a general mechanism shaping the microbial biogeography and diversity in deep-subsurface habitats across the globe.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Differential Expression of Core Metabolic Functions in Candidatus Altiarchaeum Inhabiting Distinct Subsurface Ecosystems

Candidatus Altiarchaea are widespread across aquatic subsurface ecosystems and possess a highly conserved core genome, yet adaptations of this core genome to different biotic and abiotic factors based on gene expression remain unknown. Here, we investigated the metatranscriptome of two Ca. Altiarchaeum populations that thrive in two substantially different subsurface ecosystems. In Crystal Geyser, a high-CO2 groundwater system in the USA, Ca. Altiarchaeum crystalense co-occurs with the symbiont Ca. Huberiarchaeum crystalense, while in the Muehlbacher sulfidic spring in Germany, an artesian spring high in sulfide concentration, Ca. A. hamiconexum is heavily infected with viruses. We here mapped metatranscriptome reads against their genomes to analyse the in situ expression profile of their core genomes. Out of 537 shared gene clusters, 331 were functionally annotated and 130 differed significantly in expression between the two sites. Main differences were related to genes involved in cell defence like CRISPR-Cas, virus defence, replication, transcription and energy and carbon metabolism. Our results demonstrate that altiarchaeal populations in the subsurface are likely adapted to their environment while influenced by other biological entities that tamper with their core metabolism. We consequently posit that viruses and symbiotic interactions can be major energy sinks for organisms in the deep biosphere.

archaea↗
DOE OSTI2022

Data used in manuscript "Quantifying Subsurface Flow and Solute Transport in a Snowmelt-recharged Hillslope with Multiyear Water Balance", submitted to Water Resources Research on 5-26-2022

Data in this package were used to calculate subsurface flow and solute (specifically dissolved nitrogen) transport in the Pumphouse Lower Montane (PLM) hillslope. These data are included in separate CSV files associated with each figure presented in the manuscript. The files describe site locations and elevations, weather data (air temperatures and precipitation), modeled evapotranspiration, time dependent groundwater elevations, calculated annual subsurface flow, hydraulic conductivity values, time trends in subsurface flow within different zones, pore water concentrations of total dissolve nitrogen (TDN), rates of TDN transport, soil water matric potentials following snowmelt, soil water contents during drainage, and evapotranspiration calculated by soil water mass balance versus by CLM4 model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

FY21 NNSA-IAEC Science Area V, Environmental ISR, Waste Management and Subsurface Science (Final Report)

The FY21 NNSA NA-22 Final Report represents the progress achieved over the past year (FY21) by the research team for Topic Area 3, Waste Management & Subsurface Science (WM3), and Topic Area 7, Radiation and Thermal Effects on Bituminous Rocks (WM7). WM3 and WM7 are part of Science Area V (Subsurface Science and Waste Management) in the NNSA-IAEC Science and Technology Working Group. The project was initiated four years ago (2017) with an MOU agreement between the U.S. National Nuclear Security Administration (NNSA) and the Israel Atomic Energy Commission (IAEC) to evaluate the feasibility of geological subsurface disposal of radioactive waste in Israel. The core WM3 and WM7 teams are comprised of scientists and engineers from Los Alamos National Laboratory, the Geologic Survey of Israel (GSI), and the Nuclear Research Center - Negev (NRCN), with a close collaboration to the teams from Sandia National Laboratory (SNL) and Lawrence Livermore Laboratory (LLNL)

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

Optical Fiber H2 Sensor Operating in Harsh Environments of Subsurface H2 Storage Reservoirs

Monitoring hydrogen concentration in the subsurface storage reservoirs is vital to ensure the integrity and safety of the storage facilities. An optical fiber hydrogen sensor consisting of a palladium-based sensing layer and a protective polymer layer was developed and evaluated in simulated subsurface hydrogen storage conditions. The developed optical fiber hydrogen sensor has demonstrated successful sensing performance at ~80 °C, ~1,000 psi, and ~100% RH. In addition, the sensor was exposed to real subsurface microbial samples in the harsh environments to monitor microbially induced changes in hydrogen concentration. The sensor has shown stable H2 sensing responses in the replicated underground hydrogen storage conditions without deterioration or loss of H2 sensitivity in the presence of biological samples.

filter layer↗
DOE OSTI2020

Culturing of “Unculturable” Subsurface Microbes: Natural Organic Carbon Source Fuels the Growth of Diverse and Distinct Bacteria From Groundwater

Recovery and cultivation of diverse environmentally-relevant microorganisms from the terrestrial subsurface remain a challenge despite recent advances in modern molecular technology. Here, we applied complex carbon (C) sources, i.e., sediment dissolved organic matter (DOM) and bacterial cell lysate, to enrich groundwater microbial communities for 30 days. As comparisons, we also included enrichments amended with simple C sources including glucose, acetate, benzoate, oleic acid, cellulose, and mixed vitamins. Our results demonstrate that complex C is far more effective in enriching diverse and distinct microorganisms from groundwater than simple C. Simple C enrichments yield significantly lower biodiversity, and are dominated by few phyla (e.g., Proteobacteria and Bacteroidetes ), while microcosms enriched with complex C demonstrate significantly higher biodiversity including phyla that are poorly represented in published culture collections (e.g., Verrucomicrobia , Planctomycetes , and Armatimonadetes ). Subsequent isolation from complex C enrichments yielded 228 bacterial isolates representing five phyla, 17 orders, and 56 distinct species, including candidate novel, rarely cultivated, and undescribed organisms. Results from this study will substantially advance cultivation and isolation strategies for recovering diverse and novel subsurface microorganisms. Obtaining axenic representatives of “once-unculturable” microorganisms will enhance our understanding of microbial physiology and function in different biogeochemical niches of terrestrial subsurface ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Temporal Coupling of Subsurface and Surface Soil CO 2 Fluxes: Insights From a Nonsteady State Model and Cross-Wavelet Coherence Analysis

Inferences about subsurface CO 2 fluxes often rely on surface soil respiration (R soil ) estimates because directly measuring subsurface microbial and root respiration (collectively, CO 2 production, S Total ), is difficult. To evaluate how well R soil serves as a proxy for S Total , we applied the nonsteady state DEconvolution of Temporally varying Ecosystem Carbon componenTs model (0.01-m vertical resolution), using 6-hourly data from a Wyoming grassland, in six simulations that cross three soil types (clay, sandy loam, and sandy) with two depth distributions of subsurface biota. We used cross-wavelet coherence analysis to examine temporal coherence (localized linear correlation) and offsets (lags) between S Total and R soil and fluxes and drivers (e.g., soil temperature and moisture). Cross-wavelet coherence revealed higher coherence between fluxes and drivers than linear regressions between concurrent variables. Soil texture and moisture exerted the strongest controls over coherence between CO 2 fluxes. Coherence between CO 2 fluxes in all soil types was strong at short (~1 day) and long periods (>8 days), but soil type controlled lags, and rainfall events decoupled the fluxes at periods of 1–8 days for several days in sandy soil, up to 1 week in sandy loam, and for a month or more in clay soil. Concentrating root and microbial biomass nearer the surface decreased lags in all soil types and increased coherence up to 10% in clay soil. The assumption of high temporal coherence between R soil and S Total is likely valid in dry, sandy soil, but may lead to underestimates of short-term STotal in semiarid grasslands with fine-grained and/or wet soil.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗
DOE OSTI2026

Thermo-hydro-mechanical analysis of subsurface ice-based thermal energy storage

Ice-based thermal energy storage systems are widely utilized for cooling and managing peak electrical demand globally, offering daily or weekly storage capabilities for both individual homes and larger office buildings. However, scaling these systems for district-level cooling or integrating them with renewable energy sources presents challenges, especially in accommodating larger volumes and addressing seasonal storage requirements in densely populated urban areas. This paper proposes a novel solution by evaluating subsurface ice-based thermal energy storage, in which the underground is subjected to seasonal freeze/thaw cycles. However, these cycles may influence ground behavior, affecting pore pressure and inducing ground movement. To systematically investigate these challenges, we enhance the TOUGH-FLAC simulator by integrating water/ice phase change capabilities and updating the effective stress–strain constitutive relation. Both modifications are validated against analytical solutions or experimental data. Through numerical simulations spanning a decade with ten seasonal freeze/thaw cycles, we evaluate the performance and long-term stability of a generic subsurface ice-based thermal energy storage system, considering factors such as ground permeability, freezing pipe spacing, freeze/thaw damage, and glycol solution temperature. The simulations indicate that ice formation induces pore pressure variations that drive seasonal surface heave and settlement, controlled by ground permeability, pipe spacing, and glycol solution temperature, along with tensile and localized shear deformation around freeze pipes. This highlights the need for accurate ground property characterization and geomechanical analysis for subsurface ice-based thermal energy storage.

Ground freezing↗
DOE OSTI2022

Weathering of chlorite from grain to watershed: The role and distribution of oxidation reactions in the subsurface

The reaction mechanism of weathering of chlorite, an important rock-forming phyllosilicate, is not well understood in natural settings. Here, in this work we investigated the weathering of Fe-rich chlorite from deep protolith to saprock to soil across a small shale-underlain watershed in the Appalachian Mountains, USA (Shale Hills). We found that oxidation of Fe(II) in chlorite always occurs prior to dissolution of the interlayers of the mineral. The oxidation of pyrite and chlorite commence near the water table across narrow depth intervals under the upper-catchment ridges, but well below the water table across wide depth intervals under the valley. We hypothesize that these patterns can be explained by hydrological and geochemical differences between the ridge and the valley: oxygenated water descends sub-vertically (1D flow) under the ridge, while under the valley, oxygen-depleted water moves upward to the stream and laterally out of the watershed in the subsurface (3D flow). Geochemical and mineralogical characterization indicates that the transformation of Fe-rich chlorite at Shale Hills is initiated by the oxidation of Fe(II). Next, the interlayer hydroxide sheet dissolves to form hydroxy-interlayered vermiculite and then vermiculite. During the transformation, Mg and Fe are released into solution and Fe is reprecipitated as goethite in pore space. Delivery of oxygen to the deep subsurface by infiltration of meteoric water is thought to control the initial transformation of chlorite at Shale Hills. It is possible that weathering of many Fe(II)-rich minerals is initiated by oxidation as mediated by rates of subsurface oxygen delivery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗