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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Utah FORGE: Geochemical Data for Cold Groundwaters and Produced Geothermal Fluids

Geochemical data for cold groundwaters and produced geothermal fluids around the Utah FORGE site. The data is compiled into four tables in the attached Excel File. Table 1 is a compilation of compositions (anions, cations, weak acids, oxygen, hydrogen, and carbon isotopes) for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 2 is a compilation of noble gas (He, Ne, Ar) and He and Ne isotopic compositions for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 3 provides values for calculated advective and diffusive fluxes of helium. Table 4 provides values of calculated subsurface stored heat between the Opal Mound fault and the Utah FORGE site, which are related to volumes of recently solidified magmatic heat sources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2025

Potential Risks Associated with Short-Term Hydrogen Storage in Methane Reservoirs

Hydrogen (H2) has been identified as a flexible energy carrier with zero emissions. It is possible to utilize H2 by storing Hyblend, or H2 blended with CH4, in existing natural gas infrastructure. However, the compatibility of adapting the current CH4 storage strategies to include H2 injection has not been fully demonstrated. It is essential that we understand the impact of H2 gas on the naturally occurring microbial community of subsurface storage reservoirs before deploying large-scale H2-CH4 storage. We designed a series of experiments that allowed us to identify potential challenges of HyBlend Storage in existing methane reservoirs (Figure 1). First, we characterized two field fluid samples from a methane reservoir located in Western United States. Next, we used these samples to complete a series of short-term reactor experiments at reservoir conditions for a natural gas (100% CH4) and HyBlend (80% CH4/20% H2) storage environment.

geomicrobiology↗
DOE OSTI2026

Wyoming CarbonSAFE Phase III: Site Characterization and Permitting Commercial-Scale Carbon Storage Complex Feasibility Study at Dry Fork Station, Wyoming

This report presents the findings of the technical and non-technical site characterization and permitting activities (“Phase III”) conducted under the Wyoming CarbonSAFE: Accelerating CCUS Commercialization and Deployment at Dry Fork Power Station (DFS) and the Wyoming Integrated Test Center project (“Wyoming CarbonSAFE”). Wyoming CarbonSAFE is part of the U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) Carbon Storage Assurance Facility Enterprise (“CarbonSAFE”) Initiative. The results of Phase III demonstrate that the Wyoming CarbonSAFE project - referred to as the Northern Powder River Basin Carbon Sequestration Hub (NPRB-CSH) - meets the technical, regulatory, and commercial feasibility requirements necessary to advance toward commercial development and construction. The activities completed under this project Phase make the NPRB-CSH one of the region’s most commercially ready carbon storage sites. Completion of this Phase included the finalization of all site surface and subsurface characterization activities, completion and testing of two Class VI standard wells, 10 draft Class VI permits-to construct to address the future needs of a storage complex, finalized NEPA assessments, transportation and capture FEED studies, and a full commercialization strategy with economic modeling, operation and site closure strategies. The NPRB-CSH meets all requirements to progress to a CarbonSAFE Phase IV program or advance to full commercial operations under the development of a business partner.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Carbon Utilization and Storage Partnership of the Western United States

This technical report documents research conducted under DOE Award No. DE-FE0031837 focused on evaluating the feasibility of carbon capture, utilization, and storage (CCUS) systems in the central and western United States. The project integrated geologic characterization, reservoir simulation, infrastructure modeling, and economic analysis to assess CO₂ storage potential near industrial sources and develop strategies for transport and sequestration. The work included subsurface modeling, risk assessment, monitoring and verification (MRV) planning, and evaluation of regulatory pathways such as EPA Underground Injection Control (UIC) Class VI permitting and IRS 45Q tax credit eligibility. Results demonstrate the viability of multiple storage approaches, including saline formations, enhanced coalbed methane recovery, and basalt mineralization, supported by data-driven workflows and regional analyses. The project also produced permitting templates, technology transfer activities, and stakeholder engagement efforts to support deployment readiness. These findings contribute to the development of scalable, economically viable CCUS systems and provide a repeatable framework for future carbon management projects.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Nanoscale Editing of Multi and Single Layer Tungsten Disulfide via Gas‐Assisted Focused Electron Beam Induced Etching for Device Prototyping

Focused electron beam induced etching (FEBIE) with XeF 2 (xenon difluoride) precursor is conducted on multi-layer exfoliated WS 2 (tungsten disulfide) and monolayer WS 2 grown by chemical vapor deposition (CVD). The films are characterized by atomic force microscopy (AFM) and Raman and photoluminescence (PL) spectroscopy post-etching. The etch rates/efficiencies are reported as a function of electron beam energy, current, dwell time, and XeF 2 pressure. Bulk film Raman spectra are unchanged post-FEBIE, indicating minimal subsurface damage. Monolayer WS 2 shows a decrease in Raman and PL intensity post-FEBIE, with a dose-to-clear of ≈2 nC µm −2 . The study reveals regimes affected by the various mass transport contributions such as refresh time and the ratio of electrons/XeF 2 . Spontaneous etching was discovered during FEBIE of large patterned areas due to the long frame/refresh times. Density functional theory and ab initio molecular dynamics simulations compares desorption of SF x and WF x molecules from pristine WS 2 basal planes and pore edges, revealing the spontaneous etching is consistent with etching of partially etched monolayers during each frame. Single-line etching width of 21 nm, and patterning flakes into 100 nm wide channels are demonstrated. In conclusion, this work demonstrates the possibility of editing WS 2 flakes into electronic devices of arbitrary dimensions for semiconductor applications.

2D materials↗
DOE OSTI2024

Pseudocapacitance Facilitates the Electrocatalytic Reduction of Carbon Dioxide

Electroreduction of CO 2 to value-added products is a promising strategy for CO 2 reuse, where copper has a unique ability to produce oxygenates and C 2+ products. Unfortunately, the electronic factors making copper so unique are unknown, which limits the ability to design improved catalysts. By combining in situ surface-sensitive X-ray resonant photoelectron spectroscopy with density functional theory calculations, the complex electronic structure of copper is uncovered during the electrocatalytic reduction of CO 2 . It is found that the pseudocapacitive charging of copper, produced by the incorporation of protons and electrons into the subsurface, facilitates the activation of CO 2 , while simultaneously increasing the barrier for H-H coupling. The net result is that cathodic pseudocapacitive charge suppresses the hydrogen evolution reaction and promotes the production of hydrocarbons and oxygenated products on copper. These results represent a new paradigm in the understanding of CO 2 reduction, highlighting the key role of pseudocapacitive charge in the reaction.

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

Machine‐Learning‐Driven Exploration of Surface Reconstructions of Reduced Rutile TiO 2

Abstract Titanium dioxide (TiO 2 ) is widely used as a catalyst support due to its stability, tunable electronic properties, and surface oxygen vacancies, which are crucial for catalytic processes such as the reverse water‐gas shift (RWGS) reaction. Reduced TiO 2 surfaces undergo complex surface reconstructions that endow unique properties but are computationally challenging to describe. In this study, we utilize machine‐learning interatomic potentials (MLIPs) integrated with an active‐learning workflow to efficiently explore reduced rutile TiO 2 surfaces. This approach enabled the prediction of a phase diagram as a function of oxygen chemical potential, revealing a variety of reconstructed phases, including a previously unreported subsurface shear plane structure. We further investigate the electronic properties of these surfaces and validate our results by comparing experimental and theoretical high‐resolution transmission electron microscopy (HRTEM). Our findings provide new insights into how extreme surface reductions influence the structural and electronic properties of TiO 2 , with potential implications for catalyst design.

Lee, Yonghyuk [Chemistry and Biochemistry Universi↗
DOE OSTI2025

Machine‐Learning‐Driven Exploration of Surface Reconstructions of Reduced Rutile TiO 2

Titanium dioxide (TiO 2 ) is widely used as a catalyst support due to its stability, tunable electronic properties, and surface oxygen vacancies, which are crucial for catalytic processes such as the reverse water-gas shift (RWGS) reaction. Reduced TiO 2 surfaces undergo complex surface reconstructions that endow unique properties but are computationally challenging to describe. In this study, we utilize machine-learning interatomic potentials (MLIPs) integrated with an active-learning workflow to efficiently explore reduced rutile TiO 2 surfaces. This approach enabled the prediction of a phase diagram as a function of oxygen chemical potential, revealing a variety of reconstructed phases, including a previously unreported subsurface shear plane structure. We further investigate the electronic properties of these surfaces and validate our results by comparing experimental and theoretical high-resolution transmission electron microscopy (HRTEM). Our findings provide new insights into how extreme surface reductions influence the structural and electronic properties of TiO 2 , with potential implications for catalyst design.

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

Quantifying the Site Heterogeneities of Non‐Uniform Catalysts Using QuantEXAFS

Abstract We present the Multi‐site (MS) QuantEXAFS approach, designed to model the EXAFS data from samples containing an element in different local bonding geometries. Building upon our QuantEXAFS method, which maps experimental extended X‐ray absorption fine structures (EXAFS) data to DFT‐optimized structures, MS‐QuantEXAFS introduces the key capability to probe fractional contributions of multiple sites that may be present in an experimental sample. Specifically, we demonstrate effectiveness of this technique by investigating mixed samples containing known fractions of site‐isolated subsurface Pt/MgO with Pt nanoparticles uniformly supported on MgO. The ‘site‐fractions’ obtained through MS‐QuantEXAFS closely match (i. e., ±6%) the known fractions of the physically mixed samples. This approach has been generalized to other oxides, and thus represents an important advance in quantifying the speciation of non‐uniform catalyst samples.

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

Temporal Patterns in Soil Redox Potential Vary Across a Freshwater Coastal Delta

Widespread and persistent flooding is submerging coastal ecosystems, particularly along the Louisiana Gulf coast, where flooding results from the compounding effects of ground subsidence and rising sea level. Restoration projects aim to mitigate land loss by diverting sediment loads from rivers into degraded areas to increase ground elevation. To predict how coastal ecosystems will change over time in response to projected changes in relative sea level and restoration, it is necessary to understand how subsurface biogeochemical processes respond to dynamic hydrologic forcings. Here, this study evaluates how environmental parameters that integrate biogeochemical processes vary with water table fluctuations in the freshwater Wax Lake Delta (WLD) in Louisiana, USA, where water diversions have formed one of the only active deltas along the coast. High‐frequency observations of water level, soil redox potential, specific conductance and pH were made for 1 year along elevation transects located on the older, proximal and younger, distal ends of a deltaic island. Redox responded rapidly to changing water tables, with fluctuations occurring primarily in shallow soils (< 20 cm) and at higher elevations. Deeper soils and those at lower elevation remained inundated and reduced. Semi‐diurnal tidal fluctuations were pronounced in younger, distal soils, presumably due to rapid groundwater exchange with the river channel. Tidal signals were muted in older soils that instead exhibited seasonal variability associated with river discharge and evapotranspiration. Although much of the delta sediments are persistently reducing and anoxic, redox fluctuations in the natural levees that border the deltaic islands likely drive high rates of biogeochemical activity. Evaluating how hydrology drives the frequency and duration of redox fluctuations provides a basis for understanding how biogeochemical processes might vary with complex hydrological interactions in coastal systems.

Herndon, Elizabeth [Oak Ridge National Laboratory ↗
DOE OSTI2025

Phosphorus and cover crop management practices affect phosphorus speciation in soils and eroded sediments

Abstract Agricultural runoff often contains P in dissolved and sediment‐bound forms, decreasing surface water quality. No‐till and cover cropping conservation practices have been recommended for reducing erosion and nutrient loss from cropping systems. The overall aims of this study were to characterize and evaluate the effects of fertilizer (placement and source) and cover crop management on P speciation in surface runoff sediments and source soil. In 2014, a field‐scale experiment was established in a no‐till, corn (Zea maysL.)–soybean (Glycine maxL.) cropping system with two cover crop treatments (with and without a winter crop; winter wheat [Triticum aestivumL.], rapeseed [Brassica napusL.], hairy vetch [Vicia villosaRoth], winter triticale [×Triticosecale Wittm.], and cereal rye [Secale cerealeL.]) and three P fertilizer management treatments (no P, fall broadcast diammonium phosphate, and spring subsurface injected ammonium polyphosphate). Phosphorus fractionation in the source soil collected in the fall of 2019 and sediment samples collected throughout 2020 were analyzed using a modified sequential P extraction method to evaluate the cumulative effects of imposing the treatment factors over 5 years. The direct P speciation was done using X‐ray absorption near edge structure spectroscopy. The indirect P speciation (fractionation) results showed that the management practices affected the exchangeable, organic matter‐associated, and Fe‐bound P fractions in sediments and the exchangeable and residual fractions in source soil. Direct P speciation results showed a depletion of Fe‐associated P in soil and sediment from cover crop treatment, suggesting that Fe‐associated P species were affected by cover crops. Changes in soil and runoff sediment P speciation would change the proportions and forms of soluble and particulate P in runoff sediments and may influence P bioavailability in aquatic ecosystems. Developing P fertilizer and cropping system management options with an understanding of soil P transformations helps maintain environmental sustainability.

Environmental Sciences & Ecology↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Identifying preferential flow from soil moisture time series: Review of methodologies

Abstract Identifying and quantifying preferential flow (PF) through soil—the rapid movement of water through spatially distinct pathways in the subsurface—is vital to understanding how the hydrologic cycle responds to climate, land cover, and anthropogenic changes. In recent decades, methods have been developed that use measured soil moisture time series to identify PF. Because they allow for continuous monitoring and are relatively easy to implement, these methods have become an important tool for recognizing when, where, and under what conditions PF occurs. The methods seek to identify a pattern or quantification that indicates the occurrence of PF. Most commonly, the chosen signature is either (1) a nonsequential response to infiltrated water, in which soil moisture responses do not occur in order of shallowest to deepest, or (2) a velocity criterion, in which newly infiltrated water is detected at depth earlier than is possible by nonpreferential flow processes. Alternative signatures have also been developed that have certain advantages but are less commonly utilized. Choosing among these possible signatures requires attention to their pertinent characteristics, including susceptibility to errors, possible bias toward false negatives or false positives, reliance on subjective judgments, and possible requirements for additional types of data. We review 77 studies that have applied such methods to highlight important information for readers who want to identify PF from soil moisture data and to inform those who aim to develop new methods or improve existing ones. Core Ideas Soil moisture data can be used to identify the occurrence of preferential flow (PF) and its initiating conditions. Various data‐analysis methods to identify PF differ in susceptibility to error, bias, and subjectivity. These methods can utilize vast amounts of data from soil moisture monitoring networks to develop understanding of when, where, and under what conditions PF occurs. Newly developed methods may lead to better accuracy and reliability, and reduce the need for subjective judgments. Plain Language Summary Preferential flow through soil occurs when a large amount of water is suddenly available, as during an intense storm. This type of flow moves rapidly through the soil in distinct narrow pathways rather than moving evenly throughout the body of soil, with major consequences for groundwater resources, ecosystems, spreading of contaminants, and other vital concerns. Methods of detecting preferential flow have been developed that utilize measurements of soil water content made by sensors installed at various depths. This measurement technology has been widely implemented, many locations now having datasets years in length, and various methods have been developed for using these to identify preferential flow. The various methods are based on different features in the soil moisture records and vary in their advantages and shortcomings. In this review, we explain and evaluate these methods, highlighting important information for their implementation to identify preferential flow from soil moisture data and for efforts to develop new methods or improve existing ones.

Nimmo, John R↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
DOE OSTI2026

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗
DOE OSTI2026

A comparison of long-term 37 Ar and 39 Ar transport and implications for 39 Ar as a persistent indicator of underground nuclear explosions

Measurements were first made in 2016 of legacy 39 Ar gas from the decades-old Barnwell underground nuclear explosion site within the Nevada National Security Site. This has raised the question of the ubiquity of the radionuclide 39 Ar following underground nuclear explosions in general locations. This article presents a high-level evaluation of the scientific case for using 39 Ar as a long-term indicator of underground nuclear explosions. Assessment is made based on a broad set of assumptions to reasonably estimate the persistence of 39 Ar in the subsurface environment following an explosion to establish potential detectability years to decades after the originating event.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗