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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 19

The topics include: 1) The abundances of Iron-60 in Pyroxene Chondrules from Unequilibrated Ordinary Chondrites; 2) LL-Ordinary Chondrite Impact on the Moon: Results from the 3.9 Ga Impact Melt at the Landing Site of Appolo 17; 3) Evaluation of Chemical Methods for Projectile Identification in Terrestrial and Lunar Impactites; 4) Impact Cratering Experiments in Microgravity Environment; 5) New Achondrites with High-Calcium Pyroxene and Its implication for Igneous Differentiation of Asteroids; 6) Climate History of the Polar Regions of Mars Deduced form Geologic Mapping Results; 7) The crater Production Function for Mars: A-2 Cumulative Power-Law Slope for Pristine Craters Greater than 5 km in Diameter Based on Crater Distribution for Northern Plains Materials; 8) High Resolution Al-26 Chronology: Resolved Time Interval Between Rim and Interior of a Highly Fractionated Compact Type a CAI from Efremovka; 9) Assessing Aqueous Alteration on Mars Using Global Distributions of K and Th; 10) FeNi Metal Grains in LaPaz Mare Basalt Meteorites and Appolo 12 Basalts; 11) Unique Properties of Lunar Soil for In Situ Resource Utilization on the Moon; 12) U-Pb Systematics of Phosphates in Nakhlites; 13) Measurements of Sound Speed in Granular Materials Simulated Regolith; 14) The Effects of Oxygen, Sulphur and Silicon on the Dihedral Angles Between Fe-rich Liquid Metal and Olivine, Ringwoodite and Silicate Perovskite: Implications for Planetary Core Formation; 15) Seismic Shaking Removal of Craters 0.2-0.5 km in Diameter on Asteroid 433 Eros; 16) Focused Ion Beam Microscoopy of ALH84001 Carbonate Disks; 17) Simulating Micro-Gravity in the Laboratory; 18) Mars Atmospheric Sample Return Instrument Development; 19) Combined Remote LIBS and Raman Spectroscopy Measurements; 20) Unusual Radar Backscatter Properties Along the Northern Rim of Imbrium Basin; 21) The Mars Express/NASAS Project at JPL; 22) The Geology of the Viking 2 Lander Site Revisited; 23) An Impact Genesis for Loki Patera? 24) Mars Polar Cap Edges Tracked over 3 Full Mars Years; 25) Elemental Abundance in Presolar SiC: Comparing Grains Separated by Acid Residue and Gently Separation Procedures; 26) First Results from the Descent Imager/Spectral Radiometer (DISR) Experiment on the Huygens Entry Probe of Titan; 27) Minor Element Behavior of Pallasite Olivine: Understanding Pallasite Thermal History and Chronology; 28) Canonical Anorthite in a Grosnaja Forsterite-bearing CAI; 29) Experimental Evidence for Condensation of 'Astrophysical' Carbonate; 30) Distribution and Classification of Multiple Coronae on Venus; 31) Recognition of Rayed Craters on Mars in THEMIS Thermal Infrared Imagery: Implications for Martian Meteorite Source Regions; 32) Geochemical Modeling of Evaporites on Mars: Insight from Meridiani Planum; 33) Hadean Crustal Processes Revealed from Oxygen Isotopes and U-Th-Pb Depth Profiling of Pre-4.0 Ga Detrital Zircons from Western Australia; 34) On Modeling the Seepage of Water into the Martian Subsurface; 35) Martial Gullies and Groundwater: A Series of Unfortunate Exceptions; 36) Olivine and Carbonate Globules in ALH84001: A Terrestrial Analog, and Implications for Water on Mars; 37) A Reevaluation of Mass Movements Within the Valles Marineris Region of Mars Using MOLA and MOC Data; 38) Evidence of Hydrated 109P/Swift-Tuttle Meteoroids from Meteor Spectroscopy; 39) Cr-54 Anomalies in the Solar System: Their Extent and Origin; 40) Reevaluation of the Mn-53-Cr-53 Systematic in the Basaltic Achondrites; 41) Effective Liquid Metal-Silicate Mixing Upon Shock by Power-Law Droplet Size Scaling in Richtmyer-Meshkov Like Perturbations; 42) Post-Impact Deformation of Impact Craters: Towards a Better Understanding Through the Study of Mjolnir Crater; 43) Cutting Silica Aerogel for Particle Extraction; 44) Liquid Hydrocarbons on Titan's Surface? How Cassini ISS Observations Fit into the Story (So Far); and 45) Mesoscale Simulations of Polar Circulations: Late Spring to Late Summe

Source record↗
NASA NTRS2004

Mars Reconnaissance Orbiter Navigation

Mars Reconnaissance Orbiter will launch in August 2005 at Cape Canaveral Air Force Station. The heavyweight spacecraft will use a Lockheed-Martin Atlas V-401 launch vehicle. It will be the first mission in a low Mars Orbit to characterize the surface, subsurface, and atmospheric properties. The intensive science operation imposes a great challenge for Navigation to satisfy the stringent requirements. This paper describes navigation key requirements, major challenges, and the sophisticated dynamic modeling. It also details navigation strategy and processes for various mission phases. Mars Reconnaissance Orbiter will return significant amount of scientific data in support of the objectives set by the Mars Exploration Program. A robust and precise navigation is the key to the success of this mission.

navigation↗
DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗
DOE OSTI2025

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗
DOE OSTI2023

Machine learning and deep learning for mineralogy interpretation and CO 2 saturation estimation in geological carbon Storage: A case study in the Illinois Basin

Carbon capture and storage (CCS) is a promising approach to simultaneously maintaining energy security and reducing carbon dioxide (CO 2 ) emissions under the current energy portfolio that is dominated by fossil fuel energy. Pre-injection formation characterization and post-injection CO 2 monitoring are two critical tasks to guarantee storage efficiency in CCS. The CCS projects in the Illinois Basin, the first large-scale CO 2 injection into saline aquifers in the United States, employed conventional and the latest pulsed neutron logging (PNL) tools for mineralogy interpretation and CO 2 saturation estimation, which provide valuable references for future CCS projects. Because of the inherent fuzziness of petrophysical measurements and complex subsurface heterogeneity, interpreting well-logging data is time-consuming, and its accuracy can be user-biased. In recent years, data-driven methods have been widely used to capture the non-linear patterns between input features and interpretation results. This work applied and evaluated four commonly used machine learning (ML) models, including ridge regression (RR), random forest (RF), gradient boosting regression (GBR), support vector regression (SVR), and one deep learning (DL) model, the artificial neural network (ANN). We optimized the hyperparameters of the four ML models and the DL model using the simulated annealing algorithm and the grid search strategy, respectively. The input features of the mineralogy interpretation models were eleven conventional well-logging parameters, and the label data (i.e., ground truth) were the porosity and volumetric fractions of six minerals, including quartz, feldspar, dolomite, calcite, clay, and iron minerals. The results demonstrated that the GBR and RF models were superior in predicting volumetric fractions of minerals and porosity; label data with low coefficient of variation (CV) values tended to yield better performance. For CO 2 saturation estimation, the RF was the best-performing model, followed by SVR, ANN, GBR, and RR. Furthermore, we conducted feature importance ranking using the permutation importance algorithm and found that the formation sigma and well pressure were the most important features in this study. In conclusion, the study of CCS projects in the Illinois Basin bridges the gap between the limited knowledge and understanding of geological carbon storage and the increasing demand for reliable, cost-effective, and sustainable energy solutions.

58 GEOSCIENCES↗
NASA NTRS2001

The Low-Frequency Variability of the Tropical Atlantic Ocean

Upper ocean temperature variability in the tropical Atlantic is examined from the Comprehensive Ocean Atmosphere Data Set (COADS) as well as from an ocean model simulation forced by COADS anomalies appended to a monthly climatology. Our findings are as follows: Only the sea surface temperatures (SST) in the northern tropics are driven by heat fluxes, while the southern tropical variability arises from wind driven ocean circulation changes. The subsurface temperatures in the northern and southern tropics are found to have a strong linkage to buoyancy forcing changes in the northern North Atlantic. Evidence for Kelvin-like boundary wave propagation from the high latitudes is presented from the model simulation. This extratropical influence is associated with wintertime North Atlantic Oscillation (NAO) forcing and manifests itself in the northern and southern tropical temperature anomalies of the same sign at depth of 100-200 meters as result of a Rossby wave propagation away from the eastern boundary in the wake of the boundary wave passage. The most apparent association of the southern tropical sea surface temperature anomalies (STA) arises with the anomalous cross-equatorial winds which can be related to both NAO and the remote influence from the Pacific equatorial region. These teleconnections are seasonal so that the NAO impact on the tropical SST is the largest it mid-winter but in spring and early summer the Pacific remote influence competes with NAO. However, NAO appears to have a more substantial role than the Pacific influence at low frequencies during the last 50 years. The dynamic origin of STA is indirectly confirmed from the SST-heat flux relationship using ocean model experiments which remove either anomalous wind stress forcing or atmospheric forcing anomalies contributing to heat exchange.

Haekkinen, Sirpa↗
NASA NTRS2011

Sensitivity of Depth-Integrated Satellite Lidar to Subaqueous Scattering

A method is presented for estimating subaqueous integrated backscatter from the CALIOP lidar. The algorithm takes into account specular reflection of laser light, laser scattering by wind-generated foam as well as sun glint and solar scattering from the foam Analyses show that the estimated subaqueous integrated backscatter is most sensitive to the estimate of transmittance used in the atmospheric correction, and is very insensitive to the estimate of wind speed used. As a case study, CALIOP data over Tampa Bay were compared to MODIS 645 nm remote sensing reflectance, which previously has been shown to be nearly linearly related to turbidity. The results indicate good correlation on nearly all CALIOP clear-free dates during the period 2006 through 2007, particularly those with relatively high atmospheric transmittance. When data are composited over the entire period the correlation is reduced but still statistically significant, an indication of variability in the biogeochemical composition in the water. Overall, the favorable results show promise for the application of satellite lidar integrated backscatter in providing information about subsurface backscatter properties, which can be extracted using appropriate models

Barton, Jonathan S.↗
DOE OSTI2021

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗
DOE OSTI2020

EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN

This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.

42 ENGINEERING↗
DOE OSTI2026

3D Deep Learning Joint Inversion of Active Seismic Full Waveform and Passive Seismic Traveltime Data for Reservoir Imaging and Uncertainty Quantification

Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.

Um, Evan Schankee [Lawrence Berkeley National Labo↗
DOE OSTI2021

A Deep Learning Modeling Framework to Capture Mixing Patterns in Reactive-Transport Systems

Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, in this work we will present in this paper a deep learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework—placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training. The time needed to obtain a forecast using the model is a fraction (≈ O(-6)) of the time needed to obtain the result using a high-fidelity simulation. To achieve an error of 10% (measured using the infinity norm) for capturing local-scale mixing features such as interfacial mixing, only 24% to 32% of the sequence data for model training is required. To achieve the same level of accuracy for capturing global-scale mixing features, the sequence data required for model training is 64% to 70% of the total spatial-temporal data. Hence, the proposed approach—a fast and accurate way to forecast long-time spatial-temporal mixing patterns in heterogeneous and anisotropic media—will be a valuable tool for modeling reactive-transport in a wide range of applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2018

Mission Concept for a Europa Lander

A NASA HQ-directed study team led by the Jet Propulsion Laboratory (JPL) with partners including Applied Physics Lab, Marshall Space Flight Center, Goddard Space Flight Center, Langley Research Center and Sandia National Laboratory has recently presented a mission concept for a Europa Lander that would search for bio-signatures and signs of life in the near-subsurface of the Jovian moon. The mission would follow the Europa Clipper multiple-flyby mission, planned for launch in June of 2022, which would provide reconnaissance imagery and other data to the Lander for use in selecting a scientifically compelling site and certifying it for engineering safety. The Europa Lander concept accommodates the Model Payload identified by the Europa Lander Science Definition Team (SDT) and documented in the Europa Lander Study 2016 Report released in February of 2017. The currently envisioned Europa Lander would launch on an SLS Block 1B as early as October of 2025 into a VEGA trajectory, arriving in the Jovian system as early as July of 2030. The baseline design of the integrated flight system includes a dedicated Carrier and Relay Stage, a Deorbit Vehicle composed of a Deorbit Stage consisting of a solid rocket motor (SRM), an MSL-like sky-crane Descent Stage, and a Lander which accommodates the instrument suite. The Lander would be powered by primary batteries over a 20-day surface mission. The science goals envisioned by the SDT require five samples taken from a depth of 10cm, a depth chosen to ensure minimal radiation processing of the potential biomarkers. Mission challenges include the large launch mass, unknown terrain topography, surface composition and materials properties, the high radiation environment, and complying with stringent planetary protection requirements. The mission concept uses a strategy of early risk reduction and overlapping requirements to provide robustness to harsh and uncertain environments. Early risk reduction efforts are aimed at maturing technologies associated with the sampling system, the intelligent landing system, high specific energy batteries, low mass and power motor controllers, and a thermal sterilization system.

Dooley, Jennifer↗
DOE OSTI2024

Predicting Geologic Behavior in Carbon Storage Projects Using Graph Neural Network

This study was invited to presented at NVIDIA's GTC conference to highlight the potential of Graph Neural Network as a novel and promising methodology for predicting pressure and saturation evolution in carbon storage projects. Carbon capture and storage (CCS) technology plays a pivotal role in mitigating greenhouse gas emissions, facilitating the transition to a low-carbon future. Effective management of subsurface reservoirs is essential to ensure the safe and efficient storage of captured carbon dioxide (CO₂). Accurate predictions of pressure and saturation over time are critical for evaluating the long-term performance and integrity of CCS projects. In recent years, Graph Neural Network (GNN) has emerged as a powerful framework for analyzing complex data in graph-structured domains. This abstract explores the application of GNN to forecast pressure and saturation evolution in carbon storage projects. Traditional numerical simulations of subsurface reservoirs have proven successful in providing pressure and saturation forecasts. However, these simulations involve massive amounts of computational effort and require extensive domain expertise for proper model calibration and validation. Graph Neural Operator offers an alternative approach that harnesses the inherent graph structure of reservoirs, where nodes represent reservoir grid cells and edges represent the geological connectivity between them.

Shih, Chung Yan↗
DOE OSTI2026

Daily, 30 m Resolution NDSI Data for the East River Watershed, CO for 2000-2020

This dataset contains daily Normalized Difference Snow Index (NDSI) values at 30 m spatial resolution for the East River watershed in Colorado, USA. The temporal range of these data includes water years 2001-2020. These data were created using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). This model fuses low spatial and high temporal resolution data from MODIS (500 m, daily) with high spatial and low temporal resolution data from Landsat (30 m, 16 days) to create a 30m synthetic daily snow product. This product allows for the analysis of historical snow covered area trends in the East River Watershed at fine spatiotemporal resolutions where it was not available previously. This research was performed as a part of the Department of Energy’s Subsurface Biogeochemical Research Program with the primary intent of better understanding the timing and spatial patterns of water delivery to the Critical Zone in mountain watersheds. Each .zip file contains one "water year" of data (October 1 - September 30; i.e., water year 2010 starts October 1, 2010 and ends September 30, 2011). Each zip file contains the following: STARFM daily Normalized Difference Snow Index (NDSI) fusion data files in GeoTiff format with one layer for each day between Landsat data acquisition dates (i.e., for dates of Landsat acquisition, the Landsat image is included for that date). The study area is located in an area of Landsat path overlap, so Landsat dates acquisitions are every 7-9 days. Landsat NDSI files containing the high spatial (30m), low temporal (7-9 days due to Landsat path overlap) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Dates for which no Landsat data were obtained are included as NoData layers. MODIS NDSI files containing the high temporal (daily), low spatial (500m) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Please note the MODIS data were resampled to 30m pixels for input into the STARFM model. The data have a scale factor of 10,000 and a no data value of -32767. The projection of all datasets is WGS 84 (EPSG: 4326), which has a latitude/longitude based degree resolution of 0.0002694946 X 0.0002694946, and approximates to the 30 m spatial resolution mentioned above. The Layer Index files in .csv format. They contain information for each layer in the above GeoTiff files regarding the corresponding date for each layer, the fraction of pixels in the image that contain valid data (missing data is due to either cloud cover or poor data quality; these values are not percent snow cover). Dates of Landsat overpass are indicated in these files. If no Landsat data were able to be obtained due to cloud cover or lack of Landsat Tier 1 data available on Google Earth Engine, this is also noted.

EARTH SCIENCE > CRYOSPHERE > SNOW/ICE↗
DOE OSTI2023

Machine Learning Enabled Sensor Fusion for In-Situ Defect Detection in Laser Powder Bed Fusion

Laser Powder Bed Fusion (L-PBF) Additive Manufacturing (AM) is among the metal 3D printing technologies most broadly adopted by the manufacturing industry. The current industry qualification paradigm for critical-application L-PBF parts relies heavily on expensive non-destructive inspection techniques such as X-Ray Computed Tomography (XCT), which significantly limits the use-cases of L-PBF. In situ monitoring of the process promises a less expensive alternative to ex situ testing, but existing sensor technologies and data analysis techniques struggle to detect sub-surface flaws (e.g., porosity and cracking) on production-scale L-PBF printers. RTX Technologies Research Center (RTRC) has licensed ORNL’s Peregrine software package – a printer- and camera-agnostic data analytics tool designed specifically for detecting process anomalies using in situ data collected during powder bed printing. The goal of this project was to feed temporally rich, multi-modal sensor data, including visible light, integrated near infrared (NIR), and spatially mapped co-axial melt pool thermal emission data into Peregrine to enable detection of subsurface flaws. XCT data was used as ground truth training data to allow Peregrine’s deep learning algorithms to recognize anomalies in these complex data streams in both test artifacts and industrially relevant geometries. Completion of this program has seen the successful implementation of multi-modal, multi-layer sensor data footprints for training of machine learning models in Peregrine. Flaws detected in XCT data have been successfully detected directly from this in situ data footprint, and initial analyses of the in situ probability-of-detection has been conducted, showing performance levels commensurate with traditional non-destructive evaluation (NDE) methods. The in situ monitoring methodology was then applied to an industrially relevant component that was using post-build NDE, highlighting the utility of the proposed method for hard-to-inspect AM components. As a direct result of this program, two journal manuscripts [1], [2] have been published in Additive Manufacturing, with additional manuscripts planned following program completion.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Fiber Optic Thermal Detection of Composite Delaminations

A recently developed technique is presented for thermographic detection of delaminations in composites by performing temperature measurements with fiber optic Bragg gratings. A single optical fiber with multiple Bragg gratings employed as surface temperature sensors was bonded to the surface of a composite with subsurface defects. The investigated structure was a 10-ply composite specimen with prefabricated delaminations of various sizes and depths. Both during and following the application of a thermal heat flux to the surface, the individual Bragg grating sensors measured the temporal and spatial temperature variations. The data obtained from grating sensors were analyzed with thermal modeling techniques of conventional thermography to reveal particular characteristics of the interested areas. Results were compared and found to be consistent with the calculations using numerical simulation techniques. Also discussed are methods including various heating sources and patterns, and their limitations for performing in-situ structural health monitoring.

Wu, Meng-Chou↗
DOE OSTI2020

Landslide monitoring using seismic refraction tomography – The importance of incorporating topographic variations

Seismic refraction tomography provides images of the elastic properties of subsurface materials in landslide settings. Seismic velocities are sensitive to changes in moisture content, which is a triggering factor in the initiation of many landslides. However, the application of the method to long-term monitoring of landslides is rarely used, given the challenges in undertaking repeat surveys and in handling and minimizing the errors arising from processing time-lapse surveys. This work presents a simple method and workflow for producing a reliable time-series of inverted seismic velocity models. This method is tested using data acquired during a recent, novel, long-term seismic refraction monitoring campaign at an active landslide in the UK. Potential sources of error include those arising from inaccurate and inconsistent determination of first-arrival times, inaccurate receiver positioning, and selection of inappropriate inversion starting models. At our site, a comparative analysis of variations in seismic velocity to real-world variations in topography over time shows that topographic error alone can account for changes in seismic velocity of greater than ±10% in a significant proportion (23%) of the data acquired. The seismic velocity variations arising from real material property changes at the near-surface of the landslide, linked to other sources of environmental data, are demonstrated to be of a similar magnitude. Over the monitoring period we observe subtle variations in the bulk seismic velocity of the sliding layer that are demonstrably related to variations in moisture content. This highlights the need to incorporate accurate topographic information for each time-step in the monitoring time-series. The goal of the proposed workflow is to minimize the sources of potential errors, and to preserve the changes observed by real variations in the subsurface. Following the workflow produces spatially comparable, time-lapse velocity cross-sections formulated from disparate, discretely-acquired datasets. These practical steps aim to aid the use of the seismic refraction tomography method for the long-term monitoring of landslides prone to hydrological destabilization.

Whiteley, J. S.↗
DOE OSTI2023

NNSA-IAEC Area: Workshop on Nuclear Waste Management and Subsurface Science (8-10 November 2022)

The 7 th workshop of the NNSA-IAEC collaboration was conducted on November 8-10 2022. The workshop was held at the geologic survey of Israel, Jerusalem. The IAEC hosted the workshop. The main topic of the workshop evolved around the results of the characterization borehole drilled in the Yamin plain during 2022 and associated experimental work and modeling results. The new data will be integrated for the applicability analysis of intermediate depth borehole (IDB) disposal for spent fuel in the Yamin plain. The first day of the workshop (November 8 th ) was dedicated to a geological field trip in the Judea desert depicting a similar cross section to that found in the Yamin plain. During the 9 th and 10 th of November the characterization BH results were presented, results of new work on the Miocene Hazeva fm. were presented, and this year’s outcome of the different collaboration groups were given.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗