Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface data models”

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

Thermodynamics of Minerals Stable near Earth's Surface

Project “Thermodynamics of Minerals Stable Near the Earth’s Surface” uses specialized calorimetric techniques to obtain new thermochemical data for solid phases of interest in the Earth’s “critical zone”, the near surface regions of the planet most directly affecting and affected by human activity. Major activity was in three areas: manganese oxide materials, aluminum clusters which are precursors to aluminum oxyhydroxide precipitates, and zeolites. All these classes of minerals are potential carriers of heavy metals and radionuclides in the environment, both within their crystal structures and as adsorbed layers or colloids. Their thermodynamic stability, as well as the sequence of stability of competing phases, governs the aqueous solubility of contaminants (and nutrients) in soil, groundwater, and sediments. The thermodynamics of these phases is strongly affected by particle size and crystallinity, since they often occur as very fine particles which do not coarsen in the near-surface environment. Many of these materials exhibit extensive polymorphism, multicomponent solid solution, and variable hydration states. Thus, one must define and systematize important thermodynamic variables, since one cannot make measurements on each and every possible phase, particle size, and composition. A major finding of the work on manganese oxides was labile water-induced oxidation-reduction reactions. These will be studied further by a combination of thermodynamic and synchrotron-based spectroscopic methods, the latter in collaboration with Peggy O’Day at UC Merced. The adsorption of other species, including simple organics, on such redox-active surfaces will also be investigated by a combination of thermochemical and spectroscopic studies. Studies of the energetics of water adsorption on mineral surfaces will be expanded to emphasize the adsorption of CO2 and small organic molecules, relevant to soil, groundwater, and the subsurface environment. Gas adsorption calorimetry, immersion calorimetry, and solution calorimetry will be combined with appropriate analytical and spectroscopic techniques to obtain a comprehensive picture of the interactions and their energetics. The evolution of aqueous aluminum clusters to form aluminum oxyhydroxide precipitates will be studied in detail by a combination of calorimetric and structural methods applied to aqueous samples aged at room temperature. The energetics of new classes of clusters will be explored. The project addresses several common themes relevant to the broad mission of the Department of Energy. Thermochemical data for minerals stable near the Earth’s surface, though still incomplete, are essential for understanding both natural and anthropogenic processes that affect our very existence on this planet. Such data contribute to the fundamental underpinning for mitigation and management strategies for the extraction and use of energy. Fundamental thermodynamic and mechanistic understanding is required to model processes ranging from heavy metal pollution to nuclear waste disposal to climate change.

58 GEOSCIENCES↗
DOE OSTI2023

Assessment of Carbonated Brine Injection as Low-risk Strategy for Geologic Carbon Storage

The main objective of this early-phase research was to evaluate the techno-economic feasibility and risk associated with combined brine and CO 2 storage in SWD wells using brine dissolution in the North Dakota portion of the Williston Basin. Three simulation studies were conducted to investigate: (1) CO 2 phases at different conditions, (2) wellbore compatibility, and (3) long-term storage fate in reservoir. (1) A simple reservoir model and injection simulations were created using data to represent the BEST (brine extraction and storage test) site, an operational SWD facility located near Watford City, North Dakota. The pressure evolution caused by CO 2 comingled in produced water injectate in a layer cake reservoir was then modeled while tracking aqueous CO 2 throughout the project. The salinity of the injection water, the salinity of the reservoir brine, and the amount of dissolved CO 2 comingled in the injection water were varied. (2) A wellbore corrosion model was performed using the CO 2 concentrations selected based on the reservoir modeling to examine the carbonated produced water impact on wellbore. (3) Reactive transport modeling was conducted with the optimal CO 2 concentration for this injection site to study the rock-fluid interactions and CO 2 fate in the reservoir. Results suggest that CO 2 dissolved in produced water can be injected without appreciably increasing subsurface pressure or leakage risks. Pressure buildup was found to vary with salinity but not with CO 2 mass fraction. Simulation results show that lower CO 2 percent mass fraction leads to a higher amount of CO 2 that can be dissolved at a higher injection salinity. Furthermore, the long-term goal of dissolution trapping in a traditional carbon storage project is accomplished from the start, mitigating risks associated with potential migration of buoyant CO 2 , so long as the reservoir pressure and temperature are used to determine the maximum mass fraction of the dissolved CO 2 .

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Thermal properties of Rhea's Poles: Evidence for a Meter-Deep Unconsolidated Subsurface Layer

Cassini's Composite Infrared Spectrometer (CIRS) observed both of Rhea's polar regions during a close (2000 km) flyby on 9th March 2013 during orbit 183. Rhea's southern pole was again observed during a more distant (51,000 km) flyby on 10th February 2015 during orbit 212. The results show Rhea's southern winter pole is one of the coldest places directly observed in our Solar System: surface temperatures of 25.4 +/-7.4 K and 24.7 +/-6.8 K are inferred from orbit 183 and 212 data, respectively. The surface temperature of the northern summer pole inferred from orbit 183 data is warmer: 66.6 +/-0.6 K. Assuming the surface thermophysical properties of the two polar regions are comparable then these temperatures can be considered a summer and winter seasonal temperature constraint for the polar region. Orbit 183 will provide solar longitude ( LS ) coverage at 133 deg and 313 deg for the summer and winter poles respectively, while orbit 212 provides an additional winter temperature constraint at LS 337 deg. Seasonal models with bolometric albedo values between 0.70 and 0.74 and thermal inertia values between 1 and 46 J m( exp −2) K( exp −1) s (exp −1/2) (otherwise known as MKS units) can provide adequate fits to these temperature constraints (assuming the winter temperature is an upper limit). Both these albedo and thermal inertia values agree within the uncertainties with those previously observed on both Rhea's leading and trailing hemispheres. Investigating the seasonal temperature change of Rhea's surface is particularly important, as the seasonal wave is sensitive to deeper surface temperatures ( approximately tens of centimeters to meter depths) than the more commonly reported diurnal wave (typically less than a centimeter), the exact depth difference dependent upon the assumed surface properties. For example, if a surface porosity of 0.5 and thermal inertia of 25 MKS is assumed then the depth of the seasonal thermal wave is 76 cm, which is much deeper than the approximately 0.5 cm probed by diurnal studies of Rhea ( Howett et al., 2010 ). The low thermal inertia derived here implies that Rhea's polar surfaces are highly porous even at great depths. Analysis of a CIRS focal plane 1 (10-600 cm(exp −1 ) stare observation, taken during the orbit 183 encounter between 16:22:33 and 16:23:26 UT centered on 71.7 deg W, 58.7 deg S provides the first analysis of a thermal emissivity spectrum on Rhea. The results show a flat emissivity spectrum with negligible emissivity features. A few possible explanations exist for this flat emissivity spectrum, but the most likely for Rhea is that the surface is both highly porous and composed of small particles ( less than approximately 50 micrometers).

Howett, C. J. A.↗
DOE OSTI2024

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine Learning to Identify Geologic Factors Associated with Production in Geothermal Fields: A Case-Study Using 3D Geologic Data from Brady Geothermal Field and NMFk

In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fuid-fow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fuid-fow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes. This submission includes the published journal article detailing this work, the published 3D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using the GeoThermalCloud framework, which is part of SmartTensors (both are linked below).

15 GEOTHERMAL ENERGY↗
NASA NTRS1999

MOLA-derived Roughness Data Used to Predict Surface Scattering for Mars Subsurface Radar Sounding

The Mars Express orbiter, to be launched by the European Space Agency in 2003, will carry a low-frequency radar sounding instrument, MARSIS (Mars Advanced Radar for Subsurface and Ionospheric Sounding). The primary goal of MARSIS is to map the distribution of water, both solid and (if present) liquid, in the upper several km of the martian crust. Detecting discontinuities in the crust, such as an ice-water transition, presents many challenges for a Mars orbital radar sounder. One challenge that must be overcome is the presence of radar scattering (echoes) from the surface of Mars, expected to be detected by the sounder antennas at the same time as any echoes aris-ing from subsurface interfaces. As the transmitted spherical wavefront spreads within the crust of Mars,. it also interacts with surface topography at off-nadir positions, creating a "clutter" signal that can mask the subsurface echoes. The MARSIS instrument will utilize Doppler filtering to limit the off-nadir clutter in the along-track direction, and a nadir-null secondary antenna to identify strong off-nadir clutter from the cross-track direction. To evaluate the effects of off-nadir surface clutter and the capability of these schemes to reduce the clutter, it is necessary to predict the range of scattering behavior that may be expected from martian surface topography. In this paper, we utilize Mars Orbital Laser Altimeter (MOLA) data from the current Mars Global Surveyor mission to characterize the topographic roughness of a variety of martian terrain types, at scales relevant to the MARSIS clutter problem. Segments of MOLA altimetry profiles are reduced to the topographic parameters rms slope and fractal dimension, which then are used as inputs to a near-nadir radar scattering model to predict the strength of the clutter signal. Additional information is contained in the original extended abstract.

Plaut, J. J.↗
DOE OSTI2021

Evaluation of Energy Storage Potential of Unconventional Shale Reservoirs Using Numerical Simulation of Cyclic Gas Injection

Compressed air energy storage (CAES) stores energy as compressed air in underground formations, typically salt dome caverns. When electricity demand grows, the compressed air is released through a turbine to produce electricity. CAES in the US is limited to one plant built in 1991, due in part to the inherent risk and uncertainty of developing subsurface storage reservoirs. As an alternative to CAES, we propose using some of the hundreds of thousands of hydraulically fractured horizontal wells to store energy as compressed natural gas in unconventional shale reservoirs. To store energy, produced or “sales” natural gas is injected back into the formation using excess electricity and is later produced through an expander to generate electricity. To evaluate this concept, we performed numerical simulations of cyclic natural gas injection into unconventional shale reservoirs using cmg-gem commercial reservoir modeling software. We tested short-term (diurnal) and long-term (seasonal) energy storage potential by modeling well injection and production gas flowrates as a function of bottom-hole pressure. First, we developed a conceptual model of a single fracture stage in an unconventional shale reservoir to characterize reservoir behavior during cyclic injection and production. Next, we modeled cyclic injection in the Marcellus shale gas play using published data. Results indicate that Marcellus unconventional shale reservoirs could support both short- and long-term energy storage at capacities of 100–1000 kWe per well. The results indicate that energy storage in unconventional shale gas wells may be feasible and warrants further investigation.

25 ENERGY STORAGE↗
NASA NTRS2002

Evaluation of Experimental Data from the Gains Balloon GPS Surface Reflection Instrument

The GPS Surface Reflection Instrument was integrated as an experiment on the GAINS (Global Airocean IN-situ System) 48-hour balloon mission flown in June 2002. The data collected by similar instruments in the past has been used to measure sea state from which ocean surface winds can be accurately estimated. The GPS signal has also been shown to be reflected from wetland areas and even from subsurface moisture. The current version of the instrument has been redesigned to be more compact, use less power, and withstand a greater variation in environmental conditions than previous versions. This instrument has also incorporated a new data collection mode to track 5 direct satellites (providing a continuous navigation solution) and multiplex the remaining 7 channels to track the reflected signal of the satellite tracked in channel 0. The new software mode has been shown to increase the signal to noise ratio of the collected data and enhance the science return of the instrument. During the GAINS balloon flight over the Northwest US, the instrument measured surface reflections as they were detected over the balloon's ground track. Since ground surface elevations in this area vary widely from the WGS-84 ellipsoid altitude, the instrument software has been modified to incorporate a surface altitude correction based on USGS 30-minute Digital Elevation Models. Information presented will include facts about instrument design goals, data collection methodologies and algorithms, and will focus on results of the science data analyses for the mission.

Ganoe, George G.↗
NASA NTRS2004

Evaluation of Experimental Data from the GAINS Balloon GPS Surface Reflection Instrument

The GPS Surface Reflection Instrument was integrated as an experiment on the GAINS (Global Airocean IN-situ System) 48-hour balloon mission flown in September 2001. The data collected by similar instruments in the past has been used to measure sea state from which ocean surface winds can be accurately estimated. The GPS signal has also been shown to be reflected from wetland areas and even from subsurface moisture. The current version of the instrument has been redesigned to be more compact, use less power, and withstand a greater variation in environmental conditions than previous versions. This instrument has also incorporated a new data collection mode to track 5 direct satellites (providing a continuous navigation solution) and multiplex the remaining 7 channels to track the reflected signal of the satellite tracked in channel 0. The new software mode has been shown to increase the signal to noise ratio of the collected data and enhance the science return of the instrument. During the 48-hour flight over the Northwest US, the instrument will measure surface reflections that can be detected over the balloon's ground track. Since ground surface elevations in this area vary widely from the WGS-84 ellipsoid altitude, the instrument software has been modified to incorporate a surface altitude correction based on USGS 30-minute Digital Elevation Models. Information presented will include facts about instrument design goals, data collection methodologies and algorithms, and results of the science data analyses for the 48-hour mission.

Gance, George G.↗
DOE OSTI2025

Drought-induced changes in groundwater-surface water exchange at Lake Mead area

This study focuses on the Lake Mead region in the southwestern United States, a key water reservoir serving over 25 million people and agricultural lands across several states. The area has experienced recurring anthropogenic droughts since the early 2000s. We investigate the hydrological response of the coupled surface-groundwater system to the 2020–2022 drought, one of the most severe on record. To this end, we use Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data to quantify vertical land motion caused by the elastic response of the crust to water-mass loss in the lake vicinity. Next, we apply an inverse elastic load modeling framework to quantify water loss. We further assess possible hydraulic connectivity between Lake Mead and adjacent groundwater reservoirs. We detect ground uplift of up to 8 mm/yr near the lake center, likely due to crustal rebound from reduced water-mass loading. We estimated the total water storage loss at 3.03 ± 0.25 km 3 /yr across a 3150 km 2 area surrounding Lake Mead. Groundwater accounts for approximately a third of that, being 0.94 ± 0.32 km 3 /yr. In addition, observed time lags of 6–98 days between lake and groundwater level responses, corresponding to a lateral diffusivity of 3.2–86 m 2 /s, suggest spatially variable connectivity between the lake and aquifers. These findings highlight that drought impacts propagate through the subsurface within interconnected systems, resulting in reduced buffering capacity of groundwater resources following droughts, and emphasizing the need for more integrated surface and groundwater management strategies to enhance resilience under climate and anthropogenic stressors.

58 GEOSCIENCES↗
DOE OSTI2021

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗
DOE OSTI2026

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗
NASA NTRSDate not supplied

Simulating Ordinary Chondrite (Tamdakht) Ablation at the Hypersonic Materials Environmental Test System (HyMETS) Facility.

NASA’s Science Mission Directorate (SMD) created the Asteroid Threat Assessment Project (ATAP) to inform decision-makers of risks associated with Potentially Hazardous Objects (PHOs), which may pose an existential threat to human civilization. ATAP assesses risk, in part, by developing analytical physics-based damage models which draw upon data collected during material property characterization efforts, entry simulations, hazard simulations, and ground-based arc-jet testing. Therefore, a pathfinder test campaign was con-ducted at the Hypersonic Materials Environmental Test System (HyMETS) at the NASA Langley Re-search Center to investigate the ablation mechanisms of an ordinary chondrite (Tamdakht H5) and a terrestrial analog (basalt). The HyMETS facility is a 400 kW constricted arc heater used to screen material performance under simulated aerothermal conditions of hypersonic flow.1 Facility conditions were chosen to simulate Earth entry conditions be-tween the upper-mesosphere and the lower thermo-sphere. The ablation mechanisms of Tamdakht and bas-alt are illustrated with still images collected from high-speed video cameras (fig. 1, Tamdakht). Tamdakht produces a relatively stable melt flow, as evidenced by the perpetual growth of a flange at the sidewall and the exiguous detachment of molten material into the flow. Furthermore, an inspection of high-speed video and chemical analysis of the post-test melt layer suggests that the dominant mode of mass loss for Tamdakht, at these test conditions, is expressed through the vaporization of volatiles (e.g., iron, potassium, sodium, phosphorus, and sulfur). The vaporization rates of Tamdakht are greater than basalt under the most extreme test conditions, which leads to an enhanced blowing layer, in-creased thermal shielding, and reduced recession. The ablation mechanisms of basalt are markedly different from Tamdakht and are attributed to the presence of hydrated minerals dispersed at irregular intervals throughout the silicate matrix. Rapid de-composition of secondary minerals and the subsequent formation of water lead to the ejection of the subsurface and overlying melt layer. Furthermore, the presence of water is suspected of lowering the viscosity of the melt layer resulting in increased mass loss near the edge of the test article where drag forces overcome rheological properties (fig. 1, bas-alt). Finally, the results will be discussed in the broader context of meteorite-based material response model development and the potential impact on the Asteroid Threat Assessment Project (ATAP).

Brody K. Bessire↗
DOE OSTI2024

Arctic Ocean Hydroacoustics

As the Arctic warms and loses its perennial ice cover, it is becoming more attractive for a variety of human uses. Hydroacoustic monitoring of this activity will grow in importance over the coming years and decades. Changes to the physical environment affect acoustic propagation and noise, with ramifications for our ability to detect and locate events and activities of interest. In this report, we use two long-term data sets from the Beaufort Sea, western Arctic Ocean, to determine how acoustic propagation conditions and seismic source detections are impacted by the changing environment: 1) oceanographic observations from ice-tethered profilers, and 2) passive acoustic recordings from a hydrophone. We find that changes to Beaufort Sea thermohaline stratification is stabilizing a subsurface duct, leading to more focused acoustic energy arrivals. Detections of catalogued submarine earthquakes show geographic differences in signal strength between seismic and acoustic waves. Signal strength increases with earthquake magnitude, but relationships to other source and path factors are less clear. Ambient noise also has clear seasonal patterns in the Arctic, with relatively low noise in the spring, higher noise near 1 Hz in summer, and higher noise near 10 Hz in winter. Climate change is expected to modify these seasonal noise patterns, impacting event detection. Future work will further investigate the mechanisms of ice effects on sound and couple acoustic modeling to an Earth System Model.

54 ENVIRONMENTAL SCIENCES↗
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↗
NASA NTRS2010

Active Neutron and Gamma Ray Instrumentation for In Situ Planetary Science Applications

The Pulsed Neutron Generator-Gamma Ray And Neutron Detectors (PNG-GRAND) experiment is an innovative application of the active neutron-gamma ray technology so successfully used in oil field well logging and mineral exploration on Earth. The objective of our active neutron-gamma ray technology program at NASA Goddard Space Flight Center (NASA-GSFC) is to bring the PNG-GRAND instrument to the point where it can be flown on a variety of surface lander or rover missions to the Moon, Mars, Menus, asteroids, comets and the satellites of the outer planets. Gamma-Ray Spectrometers (GRS) have been incorporated into numerous orbital planetary science missions and, especially its the case of the Mars Odyssey GRS, have contributed detailed maps of the elemental composition over the entire surface of Mars. However, orbital gamma ray measurements have low spatial sensitivity (100's of km) due to their low surface emission rates from cosmic rays and subsequent need to be averaged over large surface areas. PNG-GRAND overcomes this impediment by incorporating a powerful neutron excitation source that permits high sensitivity surface and subsurface measurements of bulk elemental compositions. PNG-GRAND combines a pulsed neutron generator (PNG) with gamma ray and neutron detectors to produce a landed instrument to determine subsurface elemental composition without needing to drill into a planet's surface a great advantage in mission design. We are currently testing PNG-GRAND prototypes at a unique outdoor neutron instrumentation test facility recently constructed at NASA/GSFC that consists of a 2 m x 2 in x 1 m granite structure placed outdoors in an empty field. Because an independent trace elemental analysis has been performed on the material, this granite sample is a known standard with which to compare both Monte Carlo simulations and our experimentally measured elemental composition data. We will present data from operating PNG-GRAND in various experimental configurations on a known sample in a geometry that is identical to that on a planetary surface. We will also illustrate the use of gamma ray timing techniques to improve sensitivity and will compare the material composition results from our experiments to both an independent laboratory elemental composition analysis and MCNPX computer modeling results.

Parsons, A.↗
DOE OSTI2024

Enhanced Attenuation of Uranium in the Subsurface

Over time, the focus for the Hanford Site has moved from production to decommissioning and site remediation. Several hundred metric tons of near-surface vadose zone sediment have been excavated from the 300 Area in an effort to remove the source of the uranium contamination (DOE/RL, 2011). Despite remediation efforts, uranium continues to be problematic with concentrations greater than the maximum contaminant level (MCL, which is 30 µg/L) for drinking water. Due to the close proximity of the Columbia River, the uranium continues to be a focus of investigation and remediation. After finding that the uranium concentrations were not decreasing as quickly as expected (based on the 1996 ROD) the Tri-Parties set new milestones, including an update of the conceptual model for the 300 Area uranium plume (Williams et al., 2007). In order to more accurately complete a feasibility study and identify how to remediate the uranium in the groundwater, the DOE authorized a remedial investigation/feasibility study (RI/FS). RI/FSs are designed to screen technologies to determine effectiveness, implementability, and cost (see Chapter 2b). These studies are designed to be conducted concurrently as part of the CERCLA process, so the data collected during site characterization (part of the RI) can inform the treatability studies and field tests conducted under the FS (EPA 1988). For the 300 Area, this study included laboratory, field, and modeling components as well as additional sediment and groundwater characterization to determine a suitable remediation strategy to reduce the uranium concentration in the groundwater plume within the expected time frame and cost.

300 area, uranium, polyphosphate, Hanford book pro↗
NASA NTRS2021

The benefit of brightness temperature assimilation for the SMAP Level-4 surface and root-zone soil moisture analysis

The Soil Moisture Active Passive (SMAP) Level-4 (L4) product provides global estimates of surface soil moisture (SSM) and root-zone soil moisture (RZSM) via the assimilation of SMAP brightness temperature (Tb) observations into the NASA Catchment Land Surface Model (CLSM). Here, using in situ measurements from 2474 sites in China, we evaluate the performance of soil moisture estimates from the L4 data assimilation (DA) system and from a baseline “open-loop” (OL) simulation of CLSM without Tb assimilation. Using random forest regression, the efficiency of the L4 DA system (i.e., the performance improvement in DA relative to OL) is attributed to eight control factors related to the CLSM as well as τ–ω radiative transfer model (RTM) components of the L4 system. Results show that the Spearman rank correlation (R) for L4 SSM with in situ measurements increases for 77 % of the in situ measurement locations (relative to that of OL), with an average R increase of approximately 14 % (ΔR=0.056). RZSM skill is improved for about 74 % of the in situ measurement locations, but the average R increase for RZSM is only 7 % (ΔR=0.034). Results further show that the SSM DA skill improvement is most strongly related to the difference between the RTM-simulated Tb and the SMAP Tb observation, followed by the error in precipitation forcing data and estimated microwave soil roughness parameter h. For the RZSM DA skill improvement, these three dominant control factors remain the same, although the importance of soil roughness exceeds that of the Tb simulation error, as the soil roughness strongly affects the ingestion of DA increments and further propagation to the subsurface. For the skill of the L4 and OL estimates themselves, the top two control factors are the precipitation error and the SSM–RZSM coupling strength error, both of which are related to the CLSM component of the L4 system. Finally, we find that the L4 system can effectively filter out errors in precipitation. Therefore, future development of the L4 system should focus on improving the characterization of the SSM–RZSM coupling strength.

Jianxiu Qiu↗