Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface–subsurface interactions

Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES↗
DOE OSTI2024

Calibration of V-Notch and Compound Weirs for Subsurface Drainage Water Level Control Structures

Highlights Accurate discharge estimation is important when evaluating edge-of-field conservation practices. V-notch weir equations were developed for three sizes of subsurface drainage water level control structures. Compound weir equations were developed for subsurface drainage water level control structures. The compound weir equation accurately estimates discharge for flows within and overtopping the V-notch. Abstract.Numerous edge-of-field conservation practices use subsurface drainage water level control structures to monitor water levels and estimate discharge. In a control structure, procedures for calculating discharge when flow depth (head) exceeds the V-notch depth and overflows in the rectangular portion of the compound weir (CW) are ambiguous. In this study, we developed calibration equations for V-notch weirs in Agri Drain inline water level control structures of different sizes for flows within the V-notch and overtopping flow events. The discharge equation for overtopping events (Q CW , L s -1 ) was determined as: Q CW = a 1 (h b1 -h 1 b1 )+a 2 (W e -W v )h 1 b2 , where h and h 1 are heads above vertex/bottom and top of V-notch (cm), respectively, W is the effective crest width of rectangular weir (cm), W v is the top width of V-notch (cm), a 1 and b 1 are parameters for V-notch weir obtained by calibration, and a 2 and b 2 are calibration parameters for rectangular weir obtained from literature. Results were compared with a weir equation available in the literature (Q V+R ), which combines a V-notch equation with a head equal to V-depth and a rectangular weir equation for flow above V-depth. Discharge at overflow was estimated with high accuracy with Q CW, whereas Q V+R underestimated discharge (e.g., PBIAS of 0.67% vs. 17.82% for a 15.2 cm structure). An example using Q V+R resulted in a 14% lower annual estimation of nitrate-N load diverted to a saturated buffer than Q CW due to underestimation of drainage discharge during overflow events. Results suggest that the developed equation (Q CW ) accurately estimates discharge and will thus improve the estimated N load compared to Q V+R . Keywords: Compound weir, Flow monitoring, Subsurface drainage, V-notch weir, Water level control structure, Weir calibration.

Agriculture↗
DOE OSTI2020

GBCGE Subsurface Database Explorer and APIs

This submission defines a DOI for the Great Basin Center for Geothermal Energy's (GBCGE) Subsurface Database Explorer web application and underlying data services, and acknowledges the INGENIOUS project as a major source of funding for data compilation and quality assurance. The GBCGE Subsurface Database Explorer is an interactive web mapping application that provides public access to the GBCGE Subsurface Database, and its collection of datasets pertinent to geothermal exploration, oil and gas exploration, critical mineral exploration, and other subsurface characterization for the Great Basin Region, western US. This is a living database, and will be continuously updated with new data and datasets as funding and motivations allow. The underlying database views that populate the web application are on an automated refresh schedule. Data sources and acknowledgements: We thank our partners with the Nevada Division of Minerals (NDOM), the Southern Methodist University (SMU), and Great Basin State Geological Surveys for their active efforts in data curation, schema design, and quality assurance. We also thank contributors among the USGS, Oregon Institute of Technology, State Divisions of Water Resources, State Divisions of Oil, Gas, and Minerals, and State Geological Surveys for open data availability and direct contributions made under the National Geothermal Data System (NGDS).

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

H2 Sensing with an Optical Fiber Sensor in the Subsurface H2 Storage Conditions

The concentration of hydrogen is subject to change due to biological reactions in the underground hydrogen storage reservoirs. Thus, monitoring hydrogen concentration in the subsurface gas deposits is vital to ensure the integrity and safety of the storage facilities. In this study, the optical fiber hydrogen sensor was developed and validated in relevant subsurface storage environments. Specifically, the sensor has demonstrated successful sensing performance at high temperatures (~80 °C) and high pressures (~1,000 psi) under very humid conditions (~100% RH). In addition, the sensor was exposed to real subsurface microbial samples 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. The biotic subsurface sample with hydrogen gas resulted in higher transmission intensity change than the abiotic sample due to the possible hydrogen consumption with microbes.

Kim, Daejin↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

The Influence of Shallow Subsurface Properties on Particle Motion in Acoustic-Seismic Coupling

Atmospheric acoustic waves transmit energy into the solid Earth through air-to-ground coupling. These waves are recorded by seismic sensors and provide insight into both atmospheric phenomena and subsurface properties. Interpreting these signals is often challenging because they are modulated by subsurface structure and the incidence angle of the acoustic wave. This study examines acoustic--seismic coupling generated by the 2012 Camp Minden Explosion, which was recorded by hundreds of seismoacoustic stations. We apply a novel technique to quantify the seismic particle motion, model coupled waves with a propagator matrix approach, and apply a Bayesian inversion to infer properties of the shallow subsurface. Our analysis reveals that prograde motion is widespread and focused in low shear-wave velocity regions, such as the Mississippi Embayment, and retrograde motion is more common in higher shear-velocity areas. Inversion results at some stations produce plausible subsurface models with strong waveform fits, while inversions at other sites are less successful. These results indicate prograde particle motion in air-to-ground coupled waves is more prevalent than previously recognized and may serve as a diagnostic for shallow velocity structure. Our comprehensive modeling and inversion framework provides a potential method to extract layered near-surface properties from acoustic-seismic coupling observations.

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 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 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↗