SEARCH · Search NASA

Results for “Subsurface characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Subsurface imaging of grain microstructure using picosecond ultrasonics

We report on imaging subsurface grain microstructure using picosecond ultrasonics. This approach relies on elastic anisotropy of crystalline materials where ultrasonic velocity depends on propagation direction relative to the crystal axes. Picosecond duration ultrasonic pulses are generated and detected using ultrashort light pulses. In materials that are transparent or semitransparent to the probe wavelength, the probe monitors GHz Brillouin oscillations. The frequency of these oscillations is related to the ultrasonic velocity and the optical index of refraction. Ultrasonic waves propagating across a grain boundary experience a change in velocity due to a change in crystallographic orientation relative to the ultrasonic propagation direction. This change in velocity is manifested as a change in the Brillouin oscillation frequency. Using the ultrasonic propagation velocity, the depth of the interface can be determined from the location in time of the transition in oscillation frequency. An image of the grain boundary is obtained by scanning the beam along the surface. We demonstrate this volumetric imaging capability using a polycrystalline UO 2 sample. As a result, cross section liftout analysis of the grain boundaries using electron microscopy were used to verify our imaging results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2004

Automated Guided-Wave Scanning Developed to Characterize Materials and Detect Defects

The Nondestructive Evaluation (NDE) Group of the Optical Instrumentation Technology Branch at the NASA Glenn Research Center has developed a scanning system that uses guided waves to characterize materials and detect defects. The technique uses two ultrasonic transducers to interrogate the condition of a material. The sending transducer introduces an ultrasonic pulse at a point on the surface of the specimen, and the receiving transducer detects the signal after it has passed through the material. The aim of the method is to correlate certain parameters in both the time and frequency domains of the detected waveform to characteristics of the material between the two transducers. The scanning system is shown. The waveform parameters of interest include the attenuation due to internal damping, waveform shape parameters, and frequency shifts due to material changes. For the most part, guided waves are used to gauge the damage state and defect growth of materials subjected to various mechanical or environmental loads. The technique has been applied to polymer matrix composites, ceramic matrix composites, and metal matrix composites as well as metallic alloys. Historically, guided wave analysis has been a point-by-point, manual technique with waveforms collected at discrete locations and postprocessed. Data collection and analysis of this type limits the amount of detail that can be obtained. Also, the manual movement of the sensors is prone to user error and is time consuming. The development of an automated guided-wave scanning system has allowed the method to be applied to a wide variety of materials in a consistent, repeatable manner. Experimental studies have been conducted to determine the repeatability of the system as well as compare the results obtained using more traditional NDE methods. The following screen capture shows guided-wave scan results for a ceramic matrix composite plate, including images for each of nine calculated parameters. The system can display up to 18 different wave parameters. Multiple scans of the test specimen demonstrated excellent repeatability in the measurement of all the guided-wave parameters, far exceeding the traditional point-by-point technique. In addition, the scan was able to detect a subsurface defect that was confirmed using flash thermography This technology is being further refined to provide a more robust and efficient software environment. Future hardware upgrades will allow for multiple receiving transducers and the ability to scan more complex surfaces. This work supports composite materials development and testing under the Ultra-Efficient Engine Technology (UEET) Project, but it also will be applied to other material systems under development for a wide range of applications.

Martin, Richard E.↗

DOE OSTI2025

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Technical note: Using long short-term memory models to fill data gaps in hydrological monitoring networks

Abstract. Quantifying the spatiotemporal dynamics in subsurface hydrological flows over a long time window usually employs a network of monitoring wells. However, such observations are often spatially sparse with potential temporal gaps due to poor quality or instrument failure. In this study, we explore the ability of recurrent neural networks to fill gaps in a spatially distributed time-series dataset. We use a well network that monitors the dynamic and heterogeneous hydrologic exchanges between the Columbia River and its adjacent groundwater aquifer at the U.S. Department of Energy's Hanford site. This 10-year-long dataset contains hourly temperature, specific conductance, and groundwater table elevation measurements from 42 wells with gaps of various lengths. We employ a long short-term memory (LSTM) model to capture the temporal variations in the observed system behaviors needed for gap filling. The performance of the LSTM-based gap-filling method was evaluated against a traditional autoregressive integrated moving average (ARIMA) method in terms of error statistics and accuracy in capturing the temporal patterns of river corridor wells with various dynamics signatures. Our study demonstrates that the ARIMA models yield better average error statistics, although they tend to have larger errors during time windows with abrupt changes or high-frequency (daily and subdaily) variations. The LSTM-based models excel in capturing both high-frequency and low-frequency (monthly and seasonal) dynamics. However, the inclusion of high-frequency fluctuations may also lead to overly dynamic predictions in time windows that lack such fluctuations. The LSTM can take advantage of the spatial information from neighboring wells to improve the gap-filling accuracy, especially for long gaps in system states that vary at subdaily scales. While LSTM models require substantial training data and have limited extrapolation power beyond the conditions represented in the training data, they afford great flexibility to account for the spatial correlations, temporal correlations, and nonlinearity in data without a priori assumptions. Thus, LSTMs provide effective alternatives to fill in data gaps in spatially distributed time-series observations characterized by multiple dominant frequencies of variability, which are essential for advancing our understanding of dynamic complex systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

For many problems in computational science and engineering, the so-called forward problem—solution of the underlying mathematical model to yield output quantities of interest, given input parameters—is difficult enough for frontier complex models, which are often characterized by coupled multiphysics and possibly stochastic behavior over a wide range of length and time scales. However, we face the need to move beyond the forward problem, to address what is often the ultimate goal of CS&E: decision making. This requires us to tackle a spectrum of mathematical problems that subsume, and are thus even more difficult than, the forward problem. First, given experimental data, we wish to estimate unknown parameters characterizing a model of a natural or engineered system by solving an inverse problem. Second, we seek the optimal configuration of the system (or experiments) by solving an optimal design problem. Third, the optimal operation of the system must be determined by solving an optimal control problem. And fourth, we must quantify uncertainties as they propagate through all of the preceding problems, from data to model inference to prediction and finally to optimal design and control. Optimal mathematical methods and algorithms for such end-to-end, data-to-decisions modeling and simulation of complex problems—requiring integrated solution of forward, inverse, optimization, and UQ problems for large, multiphysics, multiscale models—entails challenges of the highest order, while also presenting great opportunities for applied mathematics research. Unfortunately, research on forward, inverse, optimization, and UQ problems has all-too-often progressed in isolation. This has led to mathematical methods that perform well when considering the forward simulation in isolation, but are prohibitive or suboptimal or unstable when combined with other methods within the framework of inversion, optimization, or UQ. The reverse is also true: general-purpose methods developed within the inversion, optimization, and UQ fields often become prohibitive when applied to complex CS&E problems since they do not exploit their structure.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Imaging Science Requirements for a Uranus Flagship Mission

Introduction: Our presentation will summarize the requirements for UV/Visible/Near-IR remote-sensing imaging science instruments (including mapping spectrometers) envisioned for the Uranus Orbiter and Probe (UOP) mission recommended by the recent planetary decadal survey. The Uranus system offers many targets for imaging investigations. Imaging science investigations are crucial to answering fundamental questions about various physical processes that shape the Uranian system and drive its evolution. Each science objective addressed through imaging has requirements for observation and illumination geometries as well as coverage and resolution in the spatial, spectral and temporal dimensions. By combining the requirements for the diverse science targets in the Uranian system, we present a set of overarching imaging science requirements for UOP. Atmospheric Science: The study of Uranus’ atmosphere is crucial in understanding the planet's energy balance and interior heat flux. Imaging observations map how the atmosphere scatters incident solar radiation. Imaging is key to characterizing planetary-scale energy/material transport though global atmospheric dynamics and local meteorology. Hyperspectral mapping combined with radiative transfer analysis reveals the vertical distribution of clouds and hazes. Chronicling temporal variabilities of atmospheric features is key to understanding the various processes that operate within the dynamic Uranian atmosphere. A comprehensive understanding of Uranus' atmosphere is essential in deciphering the planet's origins and evolution as well as providing insight into the broader field of planetary and exoplanetary science. Magnetospheric Interactions: Imaging science is critical to understand the elusive Uranian aurorae. Mapping the aurorae and their temporal dynamics will reveal the interaction of Uranus’ complex magnetosphere with the solar wind. The apparent decades-long cooling of Uranus’ upper atmosphere is thought to be due to changes in energy input from the overlying magnetosphere, and so mapping and quantifying the auroral evidence for this energy flux is particularly relevant for upper atmospheric science In addition, Uranus’ large moons may generate auroral footprints, which can be used to constrain the planet’s unusual magnetic field structure. Satellite Interiors: Imaging science can address the potential habitability of the larger moons of Uranus by measuring their shapes and rotational states (including librations), which combined with gravity data are crucial to understanding their internal density structure. Satellite Origins: Origin markers such as organics and volatiles may be preserved on the satellite surfaces, which may be sensed by UV/IR mapping. Imaging combined with the moons’ densities may reveal whether the moons formed in Uranus’ circumplanetary disk, rings, or from ejecta of a giant impact on Uranus. Satellite Geology: Imaging investigations, including stereo imaging to obtain topography, will enable constraining each satellite’s geologic history by investigating the distributions and morphologies of impact craters, tectonic and potential cryovolcanic features, and regolith and mass wasting deposits. These records will help unravel the past complex history of orbital resonances that occurred between Uranian satellites, which may have driven geologic activity. Analyses of these features will also point to resurfacing styles and modification processes important for understanding each moon’s thermal evolution, and any present-day geological activity, including potential subsurface deep oceans or possible current activity. Additionally, hyperspectral imaging will enable mapping the surface composition and identifying the potential presence of material that may point to a subsurface ocean, such as ammonia-bearing species. Rings and Small Moons: Imaging observations elucidate the physical processes that contribute to the complex ring-moon system’s inferred unstable state. Imaging surveys will seek smaller moons and other debris that could provide evidence of past collisions, while precise astrometric data of moons will determine the current dynamical state of the system. Many of Uranus' rings are exceptionally narrow, and images of the rings covering a broad range of longitudes, and high-resolution images of features within the rings will reveal how these rings are confined and maintained, as well as illustrating the rings' internal structure and dynamics. Images of the rings and moons obtained at a variety of wavelengths and lighting conditions are also needed to quantify the composition and the size distribution of the small satellites and ring particles.

K M Sayanagi↗

DOE OSTI2012

An Assessment of Geological Carbon Storage Options in the Illinois Basin: Validation Phase

The Midwest Geological Sequestration Consortium (MGSC) assessed the options for geological carbon dioxide (CO 2 ) storage in the 155,400 km 2 (60,000 mi 2 ) Illinois Basin, which underlies most of Illinois, western Indiana, and western Kentucky. The region has annual CO 2 emissions of about 265 million metric tonnes (292 million tons), primarily from 122 coal-fired electric generation facilities, some of which burn almost 4.5 million tonnes (5 million tons) of coal per year (U.S. Department of Energy, 2010). Validation Phase (Phase II) field tests gathered pilot data to update the Characterization Phase (Phase I) assessment of options for capture, transportation, and storage of CO 2 emissions in three geological sink types: coal seams, oil fields, and saline reservoirs. Four small-scale field tests were conducted to determine the properties of rock units that control injectivity of CO 2 , assess the total storage resources, examine the security of the overlying rock units that act as seals for the reservoirs, and develop ways to control and measure the safety of injection and storage processes. The MGSC designed field test operational plans for pilot sites based on the site screening process, MVA program needs, the selection of equipment related to CO 2 injection, and design of a data acquisition system. Reservoir modeling, computational simulations, and statistical methods assessed and interpreted data gathered from the field tests. Monitoring, Verification, and Accounting (MVA) programs were established to detect leakage of injected CO 2 and ensure public safety. Public outreach and education remained an important part of the project; meetings and presentations informed public and private regional stakeholders of the results and findings. A miscible (liquid) CO 2 flood pilot project was conducted in the Clore Formation sandstone (Mississippian System, Chesterian Series) at Mumford Hills Field in Posey County, southwestern Indiana, and an immiscible CO 2 flood pilot was conducted in the Jackson sandstone (Mississippian System Big Clifty Sandstone Member) at the Sugar Creek Field in Hopkins County, western Kentucky. Up to 12% incremental oil recovery was estimated based on these pilots. A CO 2 huff puff (HNP) pilot project was conducted in the Cypress Sandstone in the Loudon Field. This pilot was designed to measure and record data that could be used to calibrate a reservoir simulation model. A pilot project at the Tanquary Farms site in Wabash County, southeastern Illinois, tested the potential storage of CO 2 in the Springfield Coal Member of the Carbondale Formation (Pennsylvanian System), in order to gauge the potential for large-scale CO 2 storage and/or enhanced coal bed methane recovery from Illinois Basin coal beds. The pilot results from all four sites showed that CO 2 could be injected into the subsurface without adversely affecting groundwater. Additionally, hydrocarbon production was enhanced, giving further evidence that CO 2 storage in oil reservoirs and coal beds offers an economic advantage. Results from the MVA program at each site indicated that injected CO 2 did not leave the injection zone. Topical reports were completed on the Middle and Late Devonian New Albany Shale and Basin CO 2 emissions. The efficacy of the New Albany Shale as a storage sink could be substantial if low injectivity concerns can be alleviated. CO 2 emissions in the Illinois Basin were projected to be dominated by coal-fired power plants.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2019

Assessment of options for the development of a stacked storage complex in the Northern Michigan Basin, USA

The Northern Michigan Basin (NMB) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Integrated Pre-Feasibility Project is part of the U.S. Department of Energy’s (DOE) CarbonSAFE initiative. Our objective of the NMB CarbonSAFE pre-feasibility project was to assess the feasibility of an integrated, commercial-scale carbon dioxide (CO 2 ) storage complex in the NMB, which includes one or more deep saline geologic formations with sufficient capacity to permanently store at least 50 million metric tons of CO 2 over a 20- to 30-year period. Subsurface injection of CO 2 into depleted oil reservoirs for enhanced oil recovery (EOR) was also evaluated in the pre-feasibility study as an additional geologic storage mechanism because there are hundreds of depleted pinnacle reef reservoirs in the NMB, along with an established CO 2 -EOR operation, a feature that distinguishes this area from most other locations in the United States. Additionally, these reefs are co-located with the most favorable saline formations in the basin – the St. Peter sandstone and the Bass Islands dolomite – making stacked storage a real possibility. The study demonstrates that there are multiple locations in the northern Michigan study area that could support a commercial-scale storage complex with either saline storage or associated storage via CO 2 -EOR in the depleted pinnacle reefs alone (i.e., each storage mechanism has sufficient capacity alone) or via a combination of the two storage mechanisms. The term stacked storage refers to the concept of using both mechanisms together, in the same area, to achieve the required storage capacity. Key sites were identified in the northern part of the southern Michigan Peninsula that had sufficient geologic storage, nearby CO 2 sources, a favorable social and economic setting, and a strong developed team.

58 GEOSCIENCES↗

DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗

DOE OSTI2023

Robust In-Situ Strain Measurements to Monitor CO 2 Storage

The goal of this project was to develop and demonstrate robust instrumentation to monitor the in-situ strain tensor in order to improve the reliability and security of CO 2 storage in geologic formations. We met the original goals of the project and the major overarching accomplishment is the advancement of strain tensor monitoring from an intriguing concept to a commercially available technology with a solid foundation of novel instruments supported by theoretical analyses and validation experiments. The main accomplishments of the project are summarized below. We designed, built and evaluated nine new optical fiber strainmeters and tiltmeters using Michelson interferometers to measure deformation with ultra-high resolution at both shallow and deep point locations in the subsurface. These are the robust strainmeters that motivated the title of the project. We designed, built and evaluated a novel method of measuring distributed strain in optical fibers with nanostrain resolution, and cm-scale location, and sampling into the seismic band. The new method is called Coherence-length-gated Microwave Photonics Interfereometry (CMPI). CMPI technology has advantages over existing commercial DAS and DSS methods. We developed and demonstrated capabilities to deploy instruments in the field and used them to measure strain caused by ambient signals like barometric pressure and tides, as well as induced signals like surface loading and pore pressure changes from pumping tests. We deployed a working strainmeter at 1,700 ft depth, slightly above an active reservoir. This is to our knowledge the greatest depth a strainmeter has been deployed and the techniques we used can readily be extended to greater depths. Optical fiber borehole tensor strainmeter techology was advanced from a TRL 4 at the start, to a TRL of 7 at the conclusion of the project. The project included advances in simulations and theoretical analyses. We developed and demonstrated a computational workflow that uses machine learning to reduce the computational requirements and make it practical to use Bayesian inversion to solve large numerical poroelastic analyses needed to interpret strain tensor field data. We evaluated the strain tensor fields and time series that would be caused by leaks of CO 2 or other fluids from reservoirs. These simulations demonstrated that signals from leaks could be measured with instruments developed for the project, opening a potentially new method for ensuring storage security. We showed that strains in caprock can be used to estimate pressure in a reservoir. This avoids the need to drill monitoring wells into the reservoir, and it expands the capabilities of monitoring in the caprock. The project includes a derivation and application of a novel analytical solution to the strains in the vicinity of a pressurized poroelastic inclusion. This solution explains field data measured during injeciton tests at the North Avant Field, and it will simplify future interpretation of strain tensor data. The project included a broad range of experiments, and of the most significant is the characterization of the strain tensor at an array three strainmeters during six injection tests at the North Avant Field, Oklahoma. This demonstrated repeatability of the strain signal measured by the new instruments developed for the project, and it showed similarities between the strain signal at shallow depths and pressure in the underlying reservoir. We also demonstrated that useful strain data can be measured at reservoir depths. This confirms that strain tensor data can be measured throughout the caprock over a reservoir. The project demonstrated the feasibility of using the strain tensor and distributed strain measured in caprock during a variety of different well tests where the pumping rate was constant, sinusoidal and positive, or a periodic square wave with zero net rate. This further strengthens the validity of using strain data to characterize reservoirs and aquifers. We also demonstrated that strain caused be fluctuations of air pressure and water pressure in the vadose zone can be measured and interpreted, suggesting that high resolution distributed strain measurements hold promise for monitoring the vadose zone. The project partially supported nine graduate students in the Environmental Engineering, Hydrogeology, Electrical Engineering programs at Clemson University. The research was described in nine journal papers, 23 talks and conference abstracts. Additional journal papers are in preparation. A new company called Tensora was started to provide strainmeter technology for commercial applications.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2014

Generic Argillite/Shale Disposal Reference Case

Radioactive waste disposal in a deep subsurface repository hosted in clay/shale/argillite is a subject of widespread interest given the desirable isolation properties, geochemically reduced conditions, and widespread geologic occurrence of this rock type (Hansen 2010; Bianchi et al. 2013). Bianchi et al. (2013) provides a description of diffusion in a clay-hosted repository based on single-phase flow and full saturation using parametric data from documented studies in Europe (e.g., ANDRA 2005). The predominance of diffusive transport and sorption phenomena in this clay media are key attributes to impede radionuclide mobility making clay rock formations target sites for disposal of high-level radioactive waste. The reports by Hansen et al. (2010) and those from numerous studies in clay-hosted underground research laboratories (URLs) in Belgium, France and Switzerland outline the extensive scientific knowledge obtained to assess long-term clay/shale/argillite repository isolation performance of nuclear waste. In the past several years under the UFDC, various kinds of models have been developed for argillite repository to demonstrate the model capability, understand the spatial and temporal alteration of the repository, and evaluate different scenarios. These models include the coupled Thermal-Hydrological-Mechanical (THM) and Thermal-Hydrological-Mechanical-Chemical (THMC) models (e.g. Liu et al. 2013; Rutqvist et al. 2014a, Zheng et al. 2014a) that focus on THMC processes in the Engineered Barrier System (EBS) bentonite and argillite host hock, the large scale hydrogeologic model (Bianchi et al. 2014) that investigates the hydraulic connection between an emplacement drift and surrounding hydrogeological units, and Disposal Systems Evaluation Framework (DSEF) models (Greenberg et al. 2013) that evaluate thermal evolution in the host rock approximated as a thermal conduction process to facilitate the analysis of design options. However, the assumptions and the properties (parameters) used in these models are different, which not only make inter-model comparisons difficult, but also compromise the applicability of the lessons learned from one model to another model. The establishment of a reference case would therefore be helpful to set up a baseline for model development. A generic salt repository reference case was developed in Freeze et al. (2013) and the generic argillite repository reference case is presented in this report. The definition of a reference case requires the characterization of the waste inventory, waste form, waste package, repository layout, EBS backfill, host rock, and biosphere. This report mainly documents the processes in EBS bentonite and host rock that are potentially important for performance assessment and properties that are needed to describe these processes, with brief description other components such as waste inventory, waste form, waste package, repository layout, aquifer, and biosphere. A thorough description of the generic argillite repository reference case will be given in Jové Colon et al. (2014).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

CHESS 2025: Discrete-return LiDAR point clouds from NEON AOP surveys

This dataset provides Level 1 (L1) discrete-return light detection and ranging (LiDAR) point cloud data collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). These data were acquired to enable characterization of vegetation structure and other three-dimensional features of the land surface, and to evaluate structural changes that may have occurred between a prior LiDAR acquisition in 2018 and the 2025 overflight. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. LiDAR data were acquired using the Optech Galaxy Prime Airborne LiDAR Terrain Mapper onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). These are the primary unclassified discrete-return LiDAR data delivered by NEON and are provided per flightline as LASzip (LAZ) 1.4 Format 6 files. Data were processed following the workflow described in the NEON L0-to-L1 Discrete Return LiDAR Algorithm Theoretical Basis Document (Krause and Goulden 2022). Each record in the unclassified point clouds represents a geolocated laser target/return recorded by the LiDAR system, with values for X, Y, Z position and return intensity. All point coordinates are provided in meters. Horizontal coordinates are referenced in Universal Transverse Mercator (UTM) zone 13N and the World Geodetic System (WGS) 1984 ensemble datum. Elevations are referenced to Geoid12A. Flight metadata describing flightline boundaries and positional uncertainty by point are also included. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

DOE OSTI2020

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2013

Persistence of uranium groundwater plumes: Contrasting mechanisms at two DOE sites in the groundwater-river interaction zone

We examine subsurface uranium (U) plumes at two U.S. Department of Energy sites that are located near large river systems and that are influenced by groundwater-river hydrologic interaction. Following surface excavation of contaminated materials, both sites were projected to naturally flush remnant uranium contamination to levels below regulatory limits (e.g., 30 µg/L or 0.126 µmol/L; U.S. EPA drinking water standard), with 10 years projected for the Hanford 300 Area (Columbia River) and 12 years for the Rifle site (Colorado River). The rate of observed uranium decrease was much lower than expected at both sites. While uncertainty remains, a comparison of current understanding suggests that the two sites have common, but also different mechanisms controlling plume persistence. At the Hanford 300 A, the persistent source is adsorbed U(VI) in the vadose zone that is released to the aquifer during spring water table excursions. The release of U(VI) from the vadose zone and its transport within the oxic, coarse-textured aquifer sediments is dominated by kinetically-limited surface complexation. Modeling implies that annual plume discharge volumes to the Columbia River are small (< one pore volume). At the Rifle site, slow oxidation of naturally reduced, contaminant U(IV) in the saturated zone and a continuous influx of U(VI) from natural, up-gradient sources influences plume persistence. Rate-limited mass transfer and surface complexation also control U(VI) migration velocity in the sub-oxic Rifle groundwater. Flux of U(VI) from the vadose zone at the Rifle site may be locally important, but it is not the dominant process that sustains the plume. A wide range in microbiologic functional diversity exists at both sites. Strains of Geobacter and other metal reducing bacteria are present at low natural abundance that are capable of enzymatic U(VI) reduction in localized zones of accumulated detrital organic carbon or after organic carbon amendment. Major differences between the sites include the geochemical nature of residual, contaminant U; the rates of current kinetic processes (both biotic and abiotic) influencing U(VI) solid-liquid distribution; the presence of detrital organic matter and the resulting spatial heterogeneity in microbially-driven redox properties; and the magnitude of groundwater hydrologic dynamics controlled by river-stage fluctuations, geologic structures, and aquifer hydraulic properties. The comparative analysis of these sites provides important guidance to the characterization, understanding, modeling, and remediation of groundwater contaminant plumes influenced by surface water interaction that are common world-wide.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2020

Geophysical borehole logging data of wells ER-GLS1, ER-GUM1, ER-PLM7, and ER-PLM8 at the East River Watershed, Colorado

The purpose of this data set is to provide a means of characterizing the bedrock variability both within a single borehole location and across the East River Watershed. The East River is part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States, and all boreholes are located within Mancos shale bedrock. The data set contains geophysical borehole logging data acquired in October 2018 in boreholes ER-GUM1 (East River- Gothic Upper Montane Borehole #1) and ER-GLS1 (East River- Gothic Lower Subalpine Borehole #1), and acquired in October 2019 in boreholes ER-PLM7 (East River- Pumphouse Lower Montane Borehole #7) and ER-PLM8 (East River- Pumphouse Lower Montane Borehole #8). It contains depth, natural gamma ray emission, spectral gamma ray emission(Th, K, U), magnetic susceptibility, fluid temperature and electrical conductivity, electromagnetic (EM) induction derived rock electrical resistivity, heat pulse and spinner flow measurements to estimate groundwater flow, and Nuclear Magnetic Resonance (NMR) derived water content and hydraulic conductivity. Data for each borehole are provided as comma delimited text files. For each borehole accurate location information is provided (based on RTK GPS (Real Time Kinematic Global Positioning System) surveying). Borehole logging data were acquired using Mount Sopris borehole tools, and Vista Clara Javelin for NMR data.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2015

Neutron Spectrometer Prospecting in the Mojave Volatiles Project Analog Field Test

We know that volatiles are sequestered at the poles of the Moon. While we have evidence of water ice and a number of other compounds based on remote sensing, the detailed distribution, and physical and chemical form are largely unknown. Additional orbital studies of lunar polar volatiles may yield further insights, but the most important next step is to use landed assets to fully characterize the volatile composition and distribution at scales of tens to hundreds of meters. To achieve this range of scales, mobility is needed. Because of the proximity of the Moon, near real-time operation of the surface assets is possible, with an associated reduction in risk and cost. This concept of operations is very different from that of rovers on Mars, and new operational approaches are required to carry out such real-time robotic exploration. The Mojave Volatiles Project (MVP) was a Moon-Mars Analog Mission Activities (MMAMA) program project aimed at (1) determining effective approaches to operating a real-time but short-duration lunar surface robotic mission, and (2) performing prospecting science in a natural setting, as a test of these approaches. Here we describe some results from the first such test, carried out in the Mojave Desert between 16 and 24 October, 2014. The test site was an alluvial fan just E of the Soda Mountains, SW of Baker, California. This site contains desert pavements, ranging from the late Pleistocene to early-Holocene in age. These pavements are undergoing dissection by the ongoing development of washes. A principal objective was to determine the hydration state of different types of desert pavement and bare ground features. The mobility element of the test was provided by the KREX-2 rover, designed and operated by the Intelligent Robotics Group at NASA Ames Research Center. The rover-borne neutron spectrometer measured the neutron albedo at both thermal and epithermal energies. Assuming uniform geochemistry and material bulk density, hydrogen as either hydroxyl/water in mineral assemblages or as moisture will significantly enhance the return of thermalized neutrons. However, in the Mojave test setting there is little uniformity, especially in bulk material density. We find that lighter toned materials (immature pavements, bar and swale, and wash materials) have lower thermal neutron flux, while mature, darker pavements with the greatest desert varnish development have higher neutron fluxes. Preliminary analysis of samples from the various terrain types in the test area indicates a prevailing moisture content of 2-3 wt% H2O. However, soil mineralogy suggests that the welldeveloped Av1 soil horizon beneath the topmost dark pavement clast layer contains the highest clay content. Structural water (including hydroxyl) in these clays may explain the enhanced neutron albedo over dark pavements. On the other hand, surface and subsurface bulk density can also play a role in neutron albedo - lower density of materials found in washes, for example, can result in a reduction in neutron flux. Analysis is ongoing.

rover-borne neutron spectrometer↗