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

Soil water content, matric potential, carbon dioxide and oxygen concentrations, Oct 2018-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of soil sensor data (temperature, water content, bulk electrical conductivity, porewater dissolved oxygen and porewater dissolved carbon dioxide) in a vertical profile from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package includes: (1) soil temperature, volumetric water content and electrical conductivity at 40, 60 and 82.5 cm depth; (2) soil matric potential at 40, 60, 79 and 100 cm depth; (3) soil CO2 concentrations at 40, 60 and 82.5 cm depth; and (4) soil oxygen concentrations at 60, 82.5, 100, 135, 170 and 182 cm depth. Both the carbon dioxide and oxygen sensors are optical sensors that can measure the partial pressure of oxygen in both saturated and unsaturated conditions. Unfortunately, soil CO2 in the profile is unexpectedly high and above the sensor calibration range (0-25,000 ppm). In addition, soil CO2 sensors failed within a year of deployment, so we only report CO2 data from 2019-2020.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 4 "*data.csv" files, one for each calendar year in the dataset. Each "*data.csv" file has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding for each measured parameter. Because sensors have been added over time, not every sensor has data dating back to Oct 2018. Note that there is a data gap over winter 2019-2020 due to a power outage. While this repository currently only contains data through December 2021, the dataset will be updated as additional years are collected and processed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Short-Term Groundwater Level Fluctuations Drive Subsurface Redox Variability

As global change processes modify the extent and functions of terrestrial-aquatic interfaces, the variability of critical and dynamic transitional zones between wetlands and uplands increases. However, it is still unclear how fluctuating water levels at these dynamic boundaries alter groundwater biogeochemical cycling. Here, we used high-temporal resolution data along gradients from wetlands to uplands and during fluctuating water levels at freshwater coastal areas to capture spatiotemporal patterns of groundwater redox potential (E h ). We observed that topography influences groundwater E h that is higher in uplands than wetlands; however, the high variability within TAI zones challenged the establishment of distinct redox zonation. Declining water levels generally decreased E h , but most locations exhibited significant E h variability, which is associated with rare instances of short-term water level fluctuations introducing oxygen. The E h -oxygen relationship showed distinct hysteresis patterns, reflecting redox poising capacity at higher E h , maintaining more oxidizing states longer than dissolved oxygen presence. Surprisingly, we observed more frequent oxidizing states in transitional areas and wetlands than in uplands. We infer that occasional oxygen entering specific wetland-upland boundaries acts as critical biogeochemical control points. Finally, high-resolution data can capture such rare yet significant biogeochemical instances, supporting redox-informed models and advancing the predictability of climate change feedback.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

3D-Laser-Scanning Technique Applied to Bulk Density Measurements of Apollo Lunar Samples

In order to better interpret gravimetric data from orbiters such as GRAIL and LRO to understand the subsurface composition and structure of the lunar crust, it is import to have a reliable database of the density and porosity of lunar materials. To this end, we have been surveying these physical properties in both lunar meteorites and Apollo lunar samples. To measure porosity, both grain density and bulk density are required. For bulk density, our group has historically utilized sub-mm bead immersion techniques extensively, though several factors have made this technique problematic for our work with Apollo samples. Samples allocated for measurement are often smaller than optimal for the technique, leading to large error bars. Also, for some samples we were required to use pure alumina beads instead of our usual glass beads. The alumina beads were subject to undesirable static effects, producing unreliable results. Other investigators have tested the use of 3d laser scanners on meteorites for measuring bulk volumes. Early work, though promising, was plagued with difficulties including poor response on dark or reflective surfaces, difficulty reproducing sharp edges, and large processing time for producing shape models. Due to progress in technology, however, laser scanners have improved considerably in recent years. We tested this technique on 27 lunar samples in the Apollo collection using a scanner at NASA Johnson Space Center. We found it to be reliable and more precise than beads, with the added benefit that it involves no direct contact with the sample, enabling the study of particularly friable samples for which bead immersion is not possible

Macke, R. J.↗
DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
DOE OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

2020 Innovation Fair Project: Is Midlatitude Ground Ice on Mars Habitable for Life?

The objective of this project was to perform a modeling study to evaluate whether midlatitude ground ice on Mars that could support life in geologically recent times. Recent studies from Mars orbital imaging, ground penetrating radar data, and infrared thermal data have revealed that large areas of Mars at latitudes higher than 40o N host near surface ground ice. Past climate cycles driven by orbital oscillations are expected to result in the distribution of ground ice changing over geologically recent history. The insolation at any location on Mars varies over time due to quasi-periodic changes in its orbital parameters (obliquity, eccentricity, perihelion longitude). The last 2 MY Mars obliquity varied between 15o and 35o but 5 MY ago varied between 24o and 48o with the average obliquity being 36o. The current low obliquity Mars atmosphere is dryer than average conditions. Recent work (Mellon and Sizemore Icarus (2022) 371;114667) has shown that repeatedly over the last 2 MY ground ice was present even at low Northern latitudes. The ice table depth is controlled by diffusive equilibrium with water vapor in the atmosphere, which also varies with orbital parameters. Determining the surface temperature and water vapor abundance requires a sophisticated climate model that was developed at NASA Ames, and results of high obliquity climate simulations will be used in the study. Work is underway to evaluate the subsurface thermal conditions in ice-rich areas of Mars during high obliquity periods to determine whether, where, and when the subsurface ice becomes warm enough to support metabolic processes or growth of terrestrial organisms.

Carol R Stoker↗
DOE OSTI2025

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
DOE OSTI2021

Stability of Pt Skin Intermetallic Core Catalysts and Adsorption Properties for the Oxygen Reduction Reaction

Density functional theory calculations were used to determine the stability of metal slabs consisting of a Pt surface monolayer and intermetallic supporting layers made from combinations of six transition metal elements (Pt, Fe, Co, Ni, Cu, and Ag), as a model system for Pt skin intermetallic core nanoparticle catalysts. The stability of the slabs is largely determined by strain at the interface of the Pt skin and the subsurface intermetallic, which was described by a lattice matching parameter (r). The surface charge on the Pt skin was found to be correlated with the average electronegativity (EN) of the intermetallic core, so this average EN was used as a descriptor for how the electronic coupling (or ligand effect) affects adsorption energies. A total of 46 slabs were investigated in terms of their stability, from which 10 stable slabs were selected for further studies of adsorbate binding (OOH*, O*, and OH*) that are intermediates in the oxygen reduction reaction (ORR). The correlation between all three adsorption energies and descriptors r and EN was found. Using a linear fit between our descriptors and the calculated adsorption energies, the overpotential for the ORR was obtained as a function of r and EN, from which a volcano plot was produced. The volcano peak was found at r = 0.96 or at EN = 2.025. Interestingly, neither r nor EN was a sufficient single reactivity descriptor as the data points were well off the general trend in both linear fits; this implies that both the strain effect and the ligand effect influence the adsorption energies, although they are partly correlated. The (r, EN) target peak parameters were used to screen over 241 intermetallic combinations of transition metal elements as active ORR activity. In conclusion, this analysis identified 11 intermetallic compounds which can support a Pt skin to have a high predicted ORR activity.

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

NNSA-IAEC Science Area V Environmental ISR: Waste Management and Subsurface Science (Final Report FY18-FY20)

In August 2017, the U.S. National Nuclear Security Administration (NNSA) and the Israel Atomic Energy Commission (IAEC) launched a joint waste management program to evaluate the feasibility of siting a radioactive waste repository at intermediate depths in Israel. The bilateral collaboration involves three National Laboratories (LLNL, LANL, and Sandia), and the Nuclear Research Center, Negev, Geological Survey of Israel, and Ben Gurion University in Israel. This research is designed to assist with the evaluation of locations for an intermediate nuclear waste repository in Israel. This final report is a high-level summary of the results from the past three years by the team that worked on the WM-2: Radionuclide colloid-facilitated transport in fractured carbonate rock project for Science Area V, Environmental ISR: Subsurface Science and Waste Management. We have completed all of our tasks and met all of our milestones. MW-2 had three integrated tasks. Task 1 was laboratory experiments investigating colloid transport of radionuclides in natural fractured carbonate rocks from the Avdat formation, Israel. Task 2 was field experiments injecting radionuclide analogues pre-sorbed to clay colloids into fractured carbonate rock, and task 3 was numerical modeling of laboratory and field experimental results to assess the importance of colloids in the migration of relevant radionuclides. This task was incorporated into tasks 1 and 2. A 4 th task was initiated in FY20, to investigate the role of organics in facilitating radionuclide transport under the same conditions explored in task 1 due to the high organic content of the rocks being investigated for the immediate borehole locations. Preliminary results will be summarized here and the work will continue in FY21-23. The overall objective of this research is to evaluate the role of colloids (naturally occurring < 1 micron particles) in facilitating the transport of long-lived radionuclides from a nuclear waste repository situated in fractured carbonate rocks. Currently little data exists on radionuclide transport in carbonate rocks, but this is the main rock type available for siting a repository in Israel. Laboratory experiments were carried out in both Israel and the U.S., field experiments have taken place in Israel, and reactive transport modeling involved LLNL, LANL and Israel partner institutions. LLNL hosted a graduate student, Emily Tran, from Ben Gurion University in FY18 and FY19 and much of the work presented here was part of her PhD research.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2002

Planetary Geophysics and Tectonics

We have carried out several studies that explore explanations for the role of chemical density variations in Moon s evolution. Meaningful models for the evolution of the Moon must explain a number of important magmatic characteristics. Volcanic activity subsequent to the formation of its anorthositic crust was dominated by the eruption of mare basalt. 1) The main phase of mare volcanism began approx. 500 Myr after the crystallization of the anorthositic crust and continued for approx. l Gyr. 2) The picitic glasses, considered to be representative of mare basalt least affected by low pressure, near-surface fractionation, were generated by melting, at 400-600 km depth, of a source containing components that, on the basis of the magma ocean hypothesis, should have crystallized at much shallower depth during fractionation of the anorthositic crust. 3) Mare basalts occur primarily in one region of the Moon. Recent topographic data demonstrate that the earlier idea that mare basalt flooded areas of low elevation is not correct. Large areas of very low elevation do not contain mare basalt. The hemispheric asymmetry of mare basalt distribution on the lunar surface must be explained in some other way. 4) A region of the surface roughly correlating with that containing mare basalts also is thought to contain high subsurface concentrations of KREEP which was excavated during the formation of large impact basins. This so-called Procellarum KREEP Terrane (PKT) is responsible for the Imbrium basin-centered thorium anomaly mapped by Lunar Prospector.

Parmentier, E. Marc↗
DOE OSTI2022

Analysis of Defects in Metal Additive Manufacturing with Augmented Data Generation

Laser powder bed fusion (LPBF) is a method of additive manufacturing (AM) that selectively melts and fuses together microscopic metallic powder. LPBF offers the benefit of producing custom structures out of high strength metals that can be difficult to fabricate with conventional methods. The challenge of LPBF is that 3D printed structures often have internal pores due to process flaws. Pulsed thermal tomography (PTT) is a method for reconstructing the depth profile of materials, allowing the visualization internal voids in solids. In prior work, we developed a convolutional neural network (CNN) which, having been trained on simulated 2D PTT images of subsurface elliptical defects, was able to classify the semi-major radii, semi-minor radii, and angular orientation of the best-fit ellipses in previously unseen PTT images. The unseen PTT images contained subsurface irregular defect shapes imported from scanning electron microscopy (SEM) images of metallic LPBF-printed specimens. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Characterization of Regolith Volatile Transport and Storage Properties by The MECA MSP 2001 Lander Payload

The diffusive and adsorptive properties of the Martian regolith influence the exchange of volatiles between the atmosphere and subsurface. Our quantitative knowledge of these properties is extremely poor -introducing substantial uncertainties in efforts to model long-term evolution of ground ice and diurnal, seasonal, and climatic cycles of CO2 and H20. This situation should significantly improve upon arrival of the 2001 Mars Surveyor Lander in 2002. In support of the Human Exploration and Development of Space (HEDS) enterprise, the 2001 mission will include a suite of instruments to characterize the nature of the Martian environment and assess whether it contains hazards that may threaten future human exploration. A major element of this effort is the Mars Environmental Compatibility Assessment (MECA) payload, which consists an optical microscopy system incorporating electrostatic, magnetic, and scratch-hardness materials testing palets, an atomic force microscope with imaging capabilities comparable to an SEM, a wet chemistry laboratory with four independent test cells, an electrometer on the robotic arm, material test patches, a camera also mounted on the arm, and a soil scoop for excavating down to about 50 cm into the soil. Although conceived to address the needs of HEDS, MECA payload is a sophisticated soil science laboratory that should provide a wealth of new data relevant to the volatile transport and storage properties of the regolith. Additional information os contained in the original.

Clifford, S. M.↗
DOE OSTI2023

Mitigation of spatial nonstationarity with vision transformers

Spatial nonstationarity, the location variance of features’ statistical distributions, is ubiquitous in many natural settings. For example, in geological reservoirs rock matrix porosity varies vertically due to geomechanical compaction trends, in mineral deposits grades vary due to sedimentation and concentration processes, in hydrology rainfall varies due to the atmosphere and topography interactions, and in metallurgy crystalline structures vary due to differential cooling. Conventional geostatistical modeling workflows rely on the assumption of stationarity to be able to model spatial features for geostatistical inference. Nevertheless, this is often not a realistic assumption when dealing with nonstationary spatial data and this has motivated a variety of nonstationary spatial modeling workflows such as trend and residual decomposition, cosimulation with secondary features, and spatial segmentation and independent modeling over stationary subdomains. The advent of deep learning technologies has enabled new workflows for modeling spatial relationships. However, there is a paucity of demonstrated best practice and general guidance on mitigation of spatial nonstationarity with deep learning in the geospatial context. We demonstrate the impact of two common types of geostatistical spatial nonstationarity on deep learning model prediction performance and propose the mitigation of such impacts using self-attention (vision transformer) models. We demonstrate the utility of vision transformers for the mitigation of nonstationarity with relative errors as low as 10%, exceeding the performance of alternative deep learning methods such as convolutional neural networks. We establish best practice by demonstrating the ability of self-attention networks for modeling large-scale spatial relationships in the presence of commonly observed geospatial nonstationarity.

58 GEOSCIENCES↗
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
NASA NTRS1997

A 70th Degree Lunar Gravity Model (GLGM-2) from Clementine and other tracking data

A spherical harmonic model of the lunar gravity field complete to degree and order 70 has been developed from S band Doppler tracking data from the Clementine mission, as well as historical tracking data from Lunar Orbiters 1-5 and the Apollo 15 and 16 subsatellites. The model combines 361,000 Doppler observations from Clementine with 347,000 historical observations. The historical data consist of mostly 60-s Doppler with a noise of 0.25 to several mm/s. The Clementine data consist of mostly 10-s Doppler data, with a data noise of 0.25 mm/s for the observations from the Deep Space Network, and 2.5 mm/s for the data from a naval tracking station at Pomonkey, Maryland. Observations provided Clementine, provide the strongest satellite constraint on the Moon's low-degree field. In contrast the historical data, collected by spacecraft that had lower periapsis altitudes, provide distributed regions of high-resolution coverage within +/- 29 deg of the nearside lunar equator. To obtain the solution for a high-degree field in the absence of a uniform distribution of observations, we applied an a priori power law constraint of the form 15 x 10(exp -5)/sq l which had the effect of limiting the gravitational power and noise at short wavelengths. Coefficients through degree and order 18 are not significantly affected by the constraint, and so the model permits geophysical analysis of effects of the major basins at degrees 10-12. The GLGM-2 model confirms major features of the lunar gravity field shown in previous gravitational field models but also reveals significantly more detail, particularly at intermediate wavelengths (10(exp 3) km). Free-air gravity anomaly maps derived from the new model show the nearside and farside highlands to be gravitationally smooth, reflecting a state of isostatic compensation. Mascon basins (including Imbrium, Serenitatis, Crisium, Smythii, and Humorum) are denoted by gravity highs first recognized from Lunar Orbiter tracking. All of the major mascons are bounded by annuli of negative anomalies representing significant subsurface mass deficiencies. Mare Orientale appears as a minor mascon surrounded by a horseshoe-shaped gravity low centered on the Inner and Outer Rook rings that is evidence of significant subsurface structural heterogeneity. Although direct tracking is not available over a significant part of the lunar farside, GLGM-2 resolves negative anomalies that correlate with many farside basins, including South Pole-Aitken, Hertzsprung, Korolev, Moscoviense, Tsiolkovsky, and Freundlich-Sharonov.

Lemonie, Frank G. R.↗
DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Illinois Storage Corridor, CarbonSAFE Phase III: UIC Class VI Permitting Plan

The Illinois Storage Corridor (ISC) project evaluated two distinct sites to determine the feasibility of commercial-scale CO₂ storage at each. The project leveraged the region's exceptional geological characteristics, particularly the well-characterized Cambro-Ordovician Storage Complex, to enable permanent geological storage of more than 50 million tonnes of CO₂ over 30 years. The two storage sites are located at One Earth Energy (OEE) facility in northcentral Illinois and Prairie State Generating Company (PSGC) in southcentral Illinois. Once operational, these facilities will combine to capture and store more than 6.5 million tonnes of CO₂ per year, positioning the ISC among the largest carbon storage regions globally. Preliminary homogeneous dynamic modeling based on regional and site-specific reservoir characteristics indicates promising injection capabilities at both locations. For the OEE site, modeling predicts a maximum allowable injection rate of 3.7 MTPA, with a baseline scenario of 1.7 MTPA over 30 years producing a CO₂ plume radius of 1.4 miles at end of injection. For the PSGC site, incorporating recent well data, modeling indicates a single-well maximum injection rate of 2.1 MTPA, with a plume radius of 4.2 miles at end of injection for the 60 MT over 30 years scenario. The primary objective of this CarbonSAFE Phase III project is to develop and submit Class VI Underground Injection Control (UIC) permit applications to the U.S. Environmental Protection Agency Region 5. The permitting plan outlines comprehensive site characterization, Area of Review delineation, monitoring programs, well construction designs, financial responsibility provisions, and post-injection site care procedures necessary to demonstrate safe, permanent CO₂ storage protective of underground sources of drinking water. Three UIC Class VI permit applications for the OEE site were submitted to EPA in October 2022 and are progressing through technical review, with final permit decision projected by June 2026. Through five rounds of Requests for Additional Information and responses, the applications have been refined to address computational modeling, area of review delineation, well integrity, monitoring protocols, and financial assurance requirements. For the PSGC site, finalized characterization and permitting documentation was delivered directly to the facility in July 2023 due to business constraints precluding formal federal regulatory submission. This comprehensive permitting effort builds upon extensive prior subsurface evaluations and demonstration projects that have confirmed the feasibility of widespread commercial-scale carbon storage in the region.

01 COAL, LIGNITE, AND PEAT↗