Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface geologic modeling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine Learning Analysis of Hydrologic Exchange Flows and Transit Time Distributions in a Large Regulated River

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. High-resolution numerical models were often used to resolve the spatial and temporal variations of exchange flows, which are computationally expensive. In this study, we adopt Random Forest (RF) and Extreme Gradient Boosting (XGB) approaches for deriving reduced order models of hydrologic exchange flows and associated transit time distributions, with integrated field observations (e.g., bathymetry) and hydrodynamic simulation data (e.g., river velocity, depth). The setup allows an improved understanding of the influences of various physical, spatial, and temporal factors on the hydrologic exchange flows and transit times. The predictors also contain those derived using hybrid clustering, leveraging our previous work on river corridor system hydromorphic classification. The machine learning-based predictive models are developed and validated along the Columbia River Corridor, and the results show that the top parameters are the thickness of the top geological formation layer, the flow regime, river velocity, and river depth; the RF and XGB models can achieve 70% to 80% accuracy and therefore are effective alternatives to the computational demanding numerical models of exchange flows and transit time distributions. Each machine learning model with its favorable configuration and setup have been evaluated. The transferability of the models to other river reaches and larger scales, which mostly depends on data availability, is also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
DOE OSTI2020

Mars Extant Life: What's Next? Conference Report

On November 5–8, 2019, the ‘‘Mars Extant Life: What’s Next?’’ conference was convened in Carlsbad, New Mexico. The conference gathered a community of actively publishing experts in disciplines related to habitability and astrobiology. Primary conclusions are as follows: A significant subset of conference attendees concluded that there is a realistic possibility that Mars hosts indigenous microbial life. A powerful theme that permeated the conference is that the key to the search for martian extant life lies in identifying and exploring refugia (‘‘oases’’), where conditions are either permanently or episodically significantly more hospitable than average. Based on our existing knowledge of Mars, conference participants highlighted four potential martian refugium (not listed in priority order): Caves, Deep Subsurface, Ices, and Salts. The conference group did not attempt to reach a consensus prioritization of these candidate environments, but instead felt that a defensible prioritization would require a future competitive process. Within the context of these candidate environments, we identified a variety of geological search strategies that could narrow the search space. Additionally, we summarized a number of measurement techniques that could be used to detect evidence of extant life (if present). Again, it was not within the scope of the conference to prioritize these measurement techniques—that is best left for the competitive process. We specifically note that the number and sensitivity of detection methods that could be implemented if samples were returned to Earth greatly exceed the methodologies that could be used at Mars. Finally, important lessons to guide extant life search processes can be derived both from experiments carried out in terrestrial laboratories and analog field sites and from theoretical modeling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the Breakdown of ENSO's Relationship With Thermocline Depth in the Central-Equatorial Pacific

Recent classification of El Nino-Southern Oscillation (ENSO) into two types, Eastern (EP) and Central (CP) events, has highlighted the importance of the central Pacific. We show here that the local correlation between ENSO subsurface temperatures (T sub ) in the upper 100-m and thermocline depth anomalies breaks down in the central equatorial Pacific, whereas T sub remains well correlated with sea surface height anomalies. This observed difference in the central equatorial Pacific is simulated by almost all climate models. It arises from a vertically slanted Tsub anomaly structure that is unique to the central equatorial Pacific. We show that this feature is an adiabatic response to wind-driving that is even present in a linear dynamic model, as long as the model has enough baroclinic modes to adequately represent the observed vertical complexity. Our findings have implications for better understanding the different importance of thermocline feedback in CP and EP events.

58 GEOSCIENCES↗
DOE OSTI2021

Thin Water Films Enable Low-Temperature Magnesite Growth Under Conditions Relevant to Geologic Carbon Sequestration

Injecting supercritical CO2 (scCO2) into basalt formations for long-term storage is a promising strategy for mitigating CO2 emissions. Mineral carbonation can result in permanent entrapment of CO2; however, carbonation kinetics in thin H2O films in humidified scCO2 is not well understood. We investigated forsterite (Mg2SiO4) carbonation to magnesite (MgCO3) via amorphous magnesium carbonate (AMC; MgCO3·xH2O, 0.5 < x < 1), with the goal to establish the fundamental controls on magnesite growth rates at low H2O activity and temperature. Experiments were conducted at 25, 40, and 50 °C in 90 bar CO2 with a H2O film thickness on forsterite that averaged 1.78 ± 0.05 monolayers. In situ infrared spectroscopy was used to monitor forsterite dissolution and the growth of AMC, magnesite, and amorphous SiO2 as a function of time. Geochemical kinetic modeling showed that magnesite was supersaturated by two to three orders of magnitude and grew according to a zero-order rate law. The results indicate that the main drivers for magnesite growth are sustained high supersaturation coupled with low H2O activity, a combination of thermodynamic conditions not attainable in bulk aqueous solution. This improved understanding of reaction kinetics can inform subsurface reactive transport models for better predictions of CO2 fate and transport.

Kerisit, Sebastien N.↗
DOE OSTI2023

Iterative multi-task learning and inference from seismic images

Seismic interpretation aims to extract quantitative and interpretable attributes from a seismic image produced using some migration method to inform characteristics of a subsurface reservoir or target of interest. Current paradigms for computing seismic attributes mostly rely on single-task algorithms. We develop an iterative, multi-task machine learning method to learn and infer multiple attributes from a seismic image. This method is composed of two stages: a multi-task inference stage and a multi-modal, multi-task refinement stage. The basic mechanism of this method is that we train a multi-task inference neural network (NN) to estimate a set of attributes, including a relative geological time (RGT), a denoised higher-resolution (DHR) seismic image, and multiple fault attributes (including probability, dip, and strike), from a low-resolution, noisy seismic image; then we input the inferred attributes to a multi-task refinement NN to enhance the raw inference results iteratively. The two multi-task NNs are trained separately based on synthetic seismic images and associated attributes generated by a geological modeling algorithm. The software we intend to release is a PyTorch implementation of this multi-task learning method for both 2D and 3D cases along with scripts to run the training/validation. The algorithm and software can be a useful tool for automatic seismic interpretation.

Gao, Kai↗
NASA NTRS2001

Geology of Lofn Crater, Callisto

Lofn crater is a 180-km-diameter impact structure in the southern cratered plains of Callisto and is among the youngest features seen on the surface. The Lofn area was imaged by the Galileo spacecraft at regional-scale resolutions (875 m/pixel), which enable the general geology to be investigated. The morphology of Lofn crater suggests that (1) it is a class of impact structure intermediate between complex craters and palimpsests or (2) it formed by the impact of a projectile which fragmented before reaching the surface, resulting in a shallow crater (even for Callisto). The asymmetric pattern of the rim and ejecta deposits suggests that the impactor entered at a low angle from the northwest. The albedo and other characteristics of the ejecta deposits from Lofn also provide insight into the properties of the icy lithosphere and subsurface configuration at the time of impact. The "target" for the Lofn impact is inferred to have included layered materials associated with the Adlinda multiring structure northwest of Loh and ejecta deposits from the Heimdall crater area to the southeast. The Lofn impact might have penetrated through these materials into a viscous substrate of ductile ice or possibly liquid water. This interpretation is consistent with models of the current interior of Callisto based on geophysical information obtained from the Galileo spacecraft.

Greeley, Ronald↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
NASA NTRS2022

Abundant phosphorus expected for possible life in Enceladus’s ocean

Saturn’s moon Enceladus has a potentially habitable subsurface water ocean that contains canonical building blocks of life (organic and inorganic carbon, ammonia, possibly hydro-gen sulfide) and chemical energy (disequilibria for methanogenesis). However, its habit-ability could be strongly affected by the unknown availability of phosphorus (P). Here, we perform thermodynamic and kinetic modeling that simulates P geochemistry based on recent insights into the geochemistry of the ocean–seafloor system on Enceladus. We find that aqueous P should predominantly exist as orthophosphate (e.g., HPO422), and total dissolved inorganic P could reach 1027to 1022mol/kg H2O, generally increasing with lower pH and higher dissolved CO2, but also depending upon dissolved ammonia and silica. Levels are much higher than<10210mol/kg H2Ofrompreviousestimatesand close to or higher than∼1026mol/kg H2O in modern Earth seawater. The high P concentration is primarily ascribed to a high (bi)carbonate concentration, which decreases the concentrations of multivalent cations via carbonate mineral formation, allowing phosphate to accumulate. Kinetic modeling of phosphate mineral dissolution suggests that geologically rapid release of P from seafloor weathering of a chondritic rocky core could supply millimoles of total dissolved P per kilogram of H2Owithin105y, much less than the likely age of Enceladus’s ocean (108to 109y). These results provide further evidence of habitable ocean conditions and show that any oceanic life would not be inhibited by low P availability.

Enceladus ocean↗
DOE OSTI2024

Machine Learning-Driven Quantification of CO2 Plume Dynamics at Illinois Basin Decatur Project Sites Using Microseismic Data

This study utilizes machine learning to quantify CO2 plume extents by analyzing microseismic data from the Illinois Basin Decatur Project (IBDP). Leveraging a unique dataset of well logs, microseismic records, and CO2 injection metrics, this work aims to predict the temporal evolution of subsurface CO2 saturation plumes. The findings illustrate that machine learning can predict plume dynamics, revealing vertical clustering of microseismic events over distinct time periods within certain proximities to the injection well, consistent with an invasion percolation model. The buoyant CO2 plume partially trapped within sandstone intervals periodically breaches localized barriers or baffles, which act as leaky seals and impede vertical migration until buoyancy overcomes gravity and capillary forces, leading to breakthroughs along vertical zones of weakness. Between different unsupervised clustering techniques, K-Means and DBSCAN were applied and analyzed in detail, where K-means outperformed DBSCAN in this specific study by indicating the combination of the highest Silhouette Score and the lowest Davies–Bouldin Index. The predictive capability of machine learning models in quantifying CO2 saturation plume extension is significant for real-time monitoring and management of CO2 sequestration sites. The models exhibit high accuracy, validated against physical models and injection data from the IBDP, reinforcing the viability of CO2 geological sequestration as a climate change mitigation strategy and enhancing advanced tools for safe management of these operations.

Iyegbekedo, Ikponmwosa↗
DOE OSTI2022

On a Unified Core Characterization Methodology to Support the Systematic Assessment of Rare Earth Elements and Critical Minerals Bearing Unconventional Carbon Ores and Sedimentary Strata

A significant gap exists in our understanding and ability to predict the spatial occurrence and extent of rare earth elements (REE) and certain critical minerals (CM) in sedimentary strata. This is largely due to a lack of existing, systematic, and well-distributed REE and CM samples and analyses in United States sedimentary basins. In addition, the type of sampling and characterization performed to date has generally lacked the resolution and approach required to constrain geologic and geographic heterogeneities typical of subsurface, mineral resources. Here, we describe a robust and systematic method for collecting core scale characterization data that can be applied to studies on the contextual and spatial attributes, the geologic history, and lithostratigraphy of sedimentary basins. The methods were developed using drilled cores from coal bearing sedimentary strata in the Powder River Basin, Wyoming (PRB). The goal of this effort is to create a unified core characterization methodology to guide systematic collection of key data to achieve a foundation of spatially and geologically constrained REEs and CMs. This guidance covers a range of measurement types and methods that are each useful either individually or in combination to support characterization and delineation of REE and CM occurrences. The methods herein, whether used in part or in full, establish a framework to guide consistent acquisition of geological, geochemical, and geospatial datasets that are key to assessing and validating REE and CM occurrences from geologic sources to support future exploration, assessment, and techno-economic related models and analyses.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2023

Data-assimilated time-lapse visco-acoustic full-waveform inversion: Theory and application for injected CO 2 plume monitoring

Continuous seismic monitoring for quantifying CO 2 plume migration and detection of any potential leakages in the subsurface is essential for the security of long-term anthropogenic carbon dioxide geologic storage. Traditional time-lapse full-waveform inversion (TLFWI) methods aim to map the CO 2 distribution by estimating seismic velocity changes, but recent studies find that CO 2 -induced attenuation is an important complement to seismic velocity for tracking the CO 2 plumes and even quantifying the CO 2 saturation. We have developed a novel data-assimilated TLFWI method to construct high-resolution time-lapse velocity and attenuation changes from dense time-lapse monitoring data. This method consists of two theoretical developments: visco-acoustic full-waveform inversion (QFWI) and multiparameter hierarchical matrix-powered extended Kalman filter (mHiEKF). The method is capable of (1) posing temporal constraints to retrieve time-lapse information from dense monitoring data by using mHiEKF, (2) accurately recovering high-spatial-resolution velocity and attenuation perturbations using first-order equation system-based QFWI, and (3) providing the model uncertainty by estimating their model standard deviation. With numerical examples, we first find the effectiveness of the new QFWI on estimating accurate velocity and attenuation models simultaneously. Then, a CO 2 leakage case and a realistic Frio-II CO 2 monitoring case are presented to find the advantages and applicability of our data-assimilated QFWI method for estimating time-lapse changes using dense time-lapse monitoring surveys. Here, by assimilating time-lapse seismic monitoring data over time, our data-assimilated QFWI method can improve the resolution of velocity and attenuation changes and decrease their model uncertainties.

58 GEOSCIENCES↗
NASA NTRS1998

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗
NASA NTRS2002

Impact Cratering Calculations

Many Martian craters are surrounded by ejecta blankets which appear to have been fluidized forming lobate and layered deposits terminated by one or more continuous distal scarps, or ramparts. One of the first hypotheses for the formation of so-called rampart ejecta features was shock-melting of subsurface ice, entrainment of liquid water into the ejecta blanket, and subsequent fluidized flow. Our work quantifies this concept. Rampart ejecta found on all but the youngest volcanic and polar regions, and the different rampart ejecta morphologies are correlated with crater size and terrain. In addition, the minimum diameter of craters with rampart features decreases with increasing latitude indicating that ice laden crust resides closer to the surface as one goes poleward on Mars. Our second goal in was to determine what strength model(s) reproduce the faults and complex features found in large scale gravity driven craters. Collapse features found in large scale craters require that the rock strength weaken as a result of the shock processing of rock and the later cratering shear flows. In addition to the presence of molten silicate in the intensely shocked region, the presence of water, either ambient, or the result of shock melting of ice weakens rock. There are several other mechanisms for the reduction of strength in geologic materials including dynamic tensile and shear induced fracturing. Fracturing is a mechanism for large reductions in strength. We found that by incorporating damage into the models that we could in a single integrated impact calculation, starting in the atmosphere produce final crater profiles having the major features found in the field measurements (central uplifts, inner ring, terracing and faulting). This was accomplished with undamaged surface strengths (0.1 GPa) and in depth strengths (1.0 GPa).

Ahrens, Thomas J.↗
DOE OSTI2025

Geologic hydrogen: From natural occurrences to anthropogenic generation – A review of fundamentals, potential, challenges and prospects

Growing demand for hydrogen is exposing the environmental and economic limits of reforming-based and carbon-managed supply chains, while the scale-up of electrolytic capacity remains capital-constrained. Geologic hydrogen, defined as molecular H₂ generated and stored within the Earth's crust offers a complementary, potentially lower-cost resource, yet exploration is still ad hoc. This review (1) revisits a global inventory of confirmed hydrogen seeps and subsurface occurrences; (2) analyzes the controlling reactions, migration pathways, and trapping conditions governing these occurrences; (3) proposes a process-based geologic hydrogen system concept analogous to, yet distinct from, the petroleum system; and (4) evaluates potential geologic hydrogen systems within the United States as a representative case study. Here, we contrast natural systems powered by serpentinization, mantle degassing or radiolysis with anthropogenic systems that stimulate the same reactions or convert in-situ hydrocarbons. Stable hydrogen accumulations require generation rates that outpace combined physical, chemical and microbial losses; the Bourakébougou field (Mali) exemplifies a self-recharging, free-gas reservoir sustained by meteoric-water serpentinization beneath an efficient caprock. Prospective geologic hydrogen resources are likely to occur in regions where iron-rich lithologies, deep-seated faults, and low-permeability sealing formations coexist. Applying this principle, we highlight three promising hydrogen play types in U.S. geological terrains: ophiolite belts (Appalachian and Californian regions), the Midcontinent Rift and the Lake Superior banded‑iron formations. Multiphysics numerical models and positive-unlabeled machine-learning workflows help to accelerate play screening and de-risk future production; yet, reaction kinetics, stimulation strategies, and full techno-economic and life-cycle assessments remain pivotal knowledge gaps.

Anthropogenic hydrogen generation↗
DOE OSTI2022

Using a Multi‐Institutional Ensemble of Watershed Models to Assess Agricultural Conservation Effectiveness in a Future Climate

This study investigates the combined impacts of climate change and agricultural conservation on the magnitude and uncertainty of nutrient loadings in the Maumee River Watershed, the second-largest watershed of the Laurentian Great Lakes. Two scenarios — baseline agricultural management and increased agricultural conservation — were assessed using an ensemble of five Soil and Water Assessment Tools driven by six climate models. The increased conservation scenario included raising conservation adoption rates from a baseline of existing conservation practices to feasible rates in the near future based on farmer surveys. This increased adoption of winter cover crops on 6%–10% to 60% of cultivated cropland; subsurface placement of phosphorus fertilizers on 35%–60% to 68% of cultivated cropland; and buffer strips intercepting runoff from 29%–34% to 50% of cultivated cropland. Increased conservation resulted in statistically significant (p ≤ 0.05) reductions in annual loads of total phosphorus (41%), dissolved reactive phosphorus (18%), and total nitrogen (14%) under the highest emission climate scenario (RCP 8.5). While nutrient loads decreased with increased conservation relative to baseline management for all watershed models, different conclusions on the true effectiveness of conservation under climate change may be drawn if only one watershed model was used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

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

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

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗