Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geological 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 739 records · Page 41

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

Rapid Analyses of Sparse Seismoacoustic Data Reveals the Timing and Size of the Accurate Energetic Systems Explosion

On 10 October 2025 an explosion occurred at a facility operated by Accurate Energetic Systems in Humphreys County, Tennessee. The incident resulted in 16 fatalities and created a debris field over several square kilometers. To address remaining questions about explosion timing and size, we collected about 20 seismic and 19 acoustic records of the blast from sensors up to hundreds of kilometers away. We then deployed 10 distinct physics-based, reduced order models (ROMs) that used validated geological structure and atmospheric conditions from the time of the event, along with observations of body- and surface-wave energy, as well as acoustic overpressure and phase duration. Each ROM predicted either timing, yield estimates, or both. We binned these estimates and their uncertainties according to each ROMs’ assumptions about confinement (aboveground, buried fully coupled, and buried partially coupled) and combined these estimates with other forensic data to conclude that the event occurred as a single, aboveground explosion on 10 December 2025 12:47:50.8 ±0.1 s with a yield equivalent to 11.8 [2.3,16.5] tons of Trinitrotoluene. Our estimates align with the Bureau of Alcohol, Tobacco, Firearms and Explosives inventory reports of 11–13 tons. This multimethod approach demonstrates the use of remotely observed geophysical data to rapidly aid conventional forensic investigations of accidental explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Application of Generic Disposal System Models

This report describes specific GDSA activities in fiscal year 2015 (FY2015) toward the development of the enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. The GDSA framework employs the PFLOTRAN thermal-hydrologic-chemical multi-physics code (Hammond et al., 2011) and the Dakota uncertainty sampling and propagation code (Adams et al., 2013). Each code is designed for massively-parallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through the engineered barriers and natural geologic barriers to a well location in an overlying or underlying aquifer. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Rock Valley Accelerated Weight Drop Preliminary P-wave Tomographic Model

An active source experiment using an accelerated weight drop was conducted in Rock Valley, Nevada National Security Site, during the spring of 2021 in order to characterize the shallow seismic structure of the region. P-wave first arrival travel times picked from this experiment were used to construct a preliminary 3-D compressional wave speed model over an area that is roughly 4 km wide east-west and 8 km north-south to a depth of about 500-600 m below the surface, but with primary data concentration along the transects of the experimental lines. The preliminary model shows good correlation with basic geology and surface features, but geological interpretation is not the focus of this report. We describe the methods used in the tomographic inversion of the data and show results from this preliminary P-wave model.

42 ENGINEERING↗

Hierarchical Fractional Advection-Dispersion Equation (FADE) to Quantify Anomalous Transport in River Corridor over a Broad Spectrum of Scales: Theory and Applications

Fractional calculus-based differential equations were found by previous studies to be promising tools in simulating local-scale anomalous diffusion for pollutants transport in natural geological media (geomedia), but efficient models are still needed for simulating anomalous transport over a broad spectrum of scales. This study proposed a hierarchical framework of fractional advection-dispersion equations (FADEs) for modeling pollutants moving in the river corridor at a full spectrum of scales. Applications showed that the fixed-index FADE could model bed sediment and manganese transport in streams at the geomorphologic unit scale, whereas the variable-index FADE well fitted bedload snapshots at the reach scale with spatially varying indices. Further analyses revealed that the selection of the FADEs depended on the scale, type of the geomedium (i.e., riverbed, aquifer, or soil), and the type of available observation dataset (i.e., the tracer snapshot or breakthrough curve (BTC)). When the pollutant BTC was used, a single-index FADE with scale-dependent parameters could fit the data by upscaling anomalous transport without mapping the sub-grid, intermediate multi-index anomalous diffusion. Pollutant transport in geomedia, therefore, may exhibit complex anomalous scaling in space (and/or time), and the identification of the FADE’s index for the reach-scale anomalous transport, which links the geomorphologic unit and watershed scales, is the core for reliable applications of fractional calculus in hydrology.

97 MATHEMATICS AND COMPUTING↗

Retrospective on Recent DOE-Funded Studies Concerning the Extraction of Rare Earth Elements & Lithium from Geothermal Brines (Final Report)

Rare earth elements (REE) and lithium are non-toxic metals that are considered critical materials due to their use in electronics, magnets, batteries, and a wide variety of industrial processes important for the economy and military preparedness. Demand for REE and lithium is increasing and these critical materials are imported, so identifying and exploiting domestic sources of REE and lithium is a national priority. The U.S. Department of Energy (DOE) Geothermal Technologies Office (GTO) has been in the forefront of sponsoring research investigating the potential recovery of REE, lithium, and other critical minerals from geothermal brines. It has been proposed that the future of geothermal energy should include “hybrid systems” that combine electricity generation with other revenue-generating activities, such as recovery of valuable and critical minerals, including REE and lithium. Two recent GTO funding opportunities have focused on the recovery of REE and other valuable minerals from geothermal brines. The research supported by the GTO’s mineral recovery program is focused on three areas: resource characterization, technology for the extraction of REE, and technology for the extraction of lithium (Tables 1 and 2). This report is a retrospective study examining the outcome of GTO’s two recent mineral recovery programs (DE-FOA-0001016 in FY 2014 and DE-FOA-0001376 in FY 2016). In this report, the knowledge, technology, and techniques that were developed by researchers funded by GTO are summarized and discussed. Four projects were funded to assess the concentrations and amounts of REE found in geothermal brines and oil field produced waters. The GTO-funded studies compiled publically available data on REE concentrations from brines and produced water from all over the USA. In addition, new samples were collected and characterized from major geothermal and hydrocarbon basins in the Western USA. The studies examined the relationship between lithology and REE concentrations and developed models examining the influence of geology on REE concentrations in produced brines. It was determined that REE are frequently found at higher concentrations in oil field produced water than geothermal brines, but that some geothermal areas had significant REE resources. Significant reservoirs of REE were identified in the Western USA. In some cases, concentrations of REE were more than 1000 times the concentrations found in seawater. Collectively, these studies represent a comprehensive picture of REE resources associated with geothermal and hydrocarbon systems in the USA. The studies did not examine lithium resources, but in some cases, lithium concentration data was collected. Data from these studies are housed in the Geothermal Data Repository (GDR) and represent a significant information resource and it is recommended that these data be further analyzed in a future study. Eight projects were funded to develop new technology for REE extraction from geothermal fluids. These projects investigated sorption as an approach for removal and recovery of REE from geothermal brines. The projects investigated cutting-edge technology for selective sorption of ions from complex solutions, including the application of metal-organic frameworks and biosorbent proteins. The REE sorption studies tested different combinations of metal-binding ligands and solid supports. The most promising metal-binding ligands for REE included phosphonic acid, thiol, and carboxylic acid functional groups. Ligands were attached or incorporated into a wide variety of solid supports. In most cases, attachment was via covalent bonding to organic resins, polymers, or silica-based supports. Most of the REE projects were conducted at a low technology readiness level (TRL) and showed promise, but direct comparison between technologies was not possible based on the available information. It is recommended that testing and reporting be standardized to the extent possible to facilitate comparisons between technologies. Two projects were directed at novel lithium extraction technology. Both projects investigated the use of inorganic sorbents, including manganese oxides. One study also examined the use of metal- ion imprinted polymers as selective ion-exchange resins for the separation of lithium and manganese from brines. Both approaches showed promise for the selective extraction of lithium from brines, including potentially geothermal brines. Results from these GTO studies indicated that selective REE and lithium extraction is possible, but interference from co-occurring solutes, such as calcium, magnesium, or heavy metals, will interfere with process efficiency and negatively impact process economics. Techno-economic analysis conducted as part of the resource and technology studies suggest extraction of REE from geothermal brines is unlikely to be economically viable, especially since non-geothermal produced waters frequently have higher REE concentrations. It is recommended that benchmarks for techno-economic analysis be established to the extent possible for future studies, to facilitate direct comparison of various technologies. Based on the collective results of this program, it appears that hybrid geothermal power would benefit more from recovery of lithium and other metals, rather than REE. It is recommended that future studies be conducted at a higher- TRL and that sorbents be tested against actual geothermal fluid samples. Prior higher-TRL efforts to extract metals from geothermal brines should be further evaluated for lessons learned.

36 MATERIALS SCIENCE↗

The geology of Mars - A brief review of some recent results

Some Mariner 9 data on geologic provinces, volcanic processes, geophysical models, Eolian activity, surface water, and atmospheric composition on Mars are presented. Ways in which this data is expected to be augmented by the Viking mission are briefly reviewed. Information given in map form includes permanent ice, layered deposits, etched plains, volcanic constructs (shields, domes, or cones), volcanic plains, moderately cratered plains, cratered plains, hummocky terrain, channel deposit, undivided plains, and grooved terrain.

Mutch, T. A.↗

The evolution of impact basins - Viscous relaxation of topographic relief

A topographic profile of the young large lunar basin, Orientale, is presented in order to examine the effects of viscous relaxation on basin topography. Analytical models for viscous flow are considered, showing a wavelength-dependence of time constants for viscous decay on the decrease in viscosity with depth and on the extent of the isostatic compensation of the initial topography. Lunar rheological models which are developed include a half-space model for uniform Newtonian viscosity, density, and gravitational acceleration, a layer over inviscid half space model with material inviscid over geological time scales, and a layer with isostatic compensation where a uniformly viscous layer overlies an inviscid half space of higher density. Greater roughness is concluded, and has been observed, on the moon's dark side due to continued lower temperatures since the time of heavy bombardment.

Solomon, S. C.↗

Geological implications of impacts of large asteroids and comets on the earth

The present conference discusses such topics as large object fluxes in near-earth space and the probabilities of terrestrial impacts, the geological record of impacts, dynamics modeling for large body impacts on continents and oceans, physical, chemical, and biological models of large impacts' atmospheric effects, dispersed impact ejecta and their signatures, general considerations concerning mass biological extinctions, the Cretaceous/Tertiary boundary event, geochemical signatures in the stratigraphic record, and other phanerozoic events. Attention is given to terrestrial impact rates for long- and short-period comets, estimates of crater size for large body impact, a first-order estimate of shock heating and vaporization in oceanic impacts, atmospheric effects in the first few minutes after an impact, a feasibility test for biogeographic extinction, and the planktonic and dinosaur extinctions.

Silver, L. T.↗

Water Content of Earth's Continental Mantle Is Controlled by the Circulation of Fluids or Melts

A key mission of the ARES Directorate at JSC is to constrain models of the formation and geological history of terrestrial planets. Water is a crucial parameter to be measured with the aim to determine its amount and distribution in the interior of Earth, Mars, and the Moon. Most of that "water" is not liquid water per se, but rather hydrogen dissolved as a trace element in the minerals of the rocks at depth. Even so, the middle layer of differentiated planets, the mantle, occupies such a large volume and mass of each planet that when it is added at the planetary scale, oceans worth of water could be stored in its interior. The mantle is where magmas originate. Moreover, on Earth, the mantle is where the boundary between tectonic plates and the underlying asthenosphere is located. Even if mantle rocks in Earth typically contain less than 200 ppm H2O, such small quantities have tremendous influence on how easily they melt (i.e., the more water there is, the more magma is produced) and deform (the more water there is, the less viscous they are). These two properties alone emphasize that to understand the distribution of volcanism and the mechanism of plate tectonics, the water content of the mantle must be determined - Earth being a template to which all other terrestrial planets can be compared.

Peslier, Anne↗

The Permanently Shadowed Regions of Dwarf Planet Ceres

Ceres has only a small spin axis tilt (4 deg), and craters near its rotational poles can experience permanent shadow and trap volatiles, as is the case on Mercury and on Earth's Moon. Topography derived from stereo imaging by the Dawn spacecraft is used to calculate direct solar irradiance that defines the extent of the permanently shadowed regions (PSRs). In the northern polar region, PSRs cover approximately 1800 sq km or 0.13% of the hemisphere, and most of the PSRs are cold enough to trap water ice over geological time periods. Based on modeling of the water exosphere, water molecules seasonally reside around the winter pole and ultimately an estimated 0.14% of molecules get trapped. Even for the lowest estimates of the amount of available water, this predicts accumulation rates in excess of loss rates, and hence, there should be fresh ice deposits in the cold traps.

CERES↗

Characterization of volcanic tuff pores pre- and post-underground nuclear detonation using ultra-small and small angle neutron scattering

The ability to accurately model the subsurface transport of radionuclides is fundamental to the remote detection and characterization of underground nuclear explosion (UNE) events. Developing more sophisticated transport models presents a significant opportunity to enhance monitoring capabilities, particularly in the reliable prediction of signature migration. Experimentally determined characterization of geologic materials associated with transport properties is the pertinent base information for such robust model development and calibration. Here, we report results from an unprecedented study demonstrating changes to the pore and fracture network structures in geological materials in response to UNEs over nanometer to micrometer length scales. Volcanic tuffs of five different lithological formations from pre- and post-UNE environments were collected from the Nevada National Security Site. Combined ultra-small and small-angle neutron scattering techniques were used to characterize the tuff pore structure. The results demonstrate measurable differences in the specific surface area and porosity of samples pre- and post-shot from texturally similar lithological formations, indicating that pore properties can serve as a direct physical signature of a UNE. The results also provide experimentally determined transport parameters in support of advanced model development through the integration of gas migration, hydrodynamic simulations, and geologic framework models.

54 ENVIRONMENTAL SCIENCES↗

Ellicott City Disasters II: Enhancing a Statistical Flood Risk Model to Continue Improving Early Warning Systems and Public Safety in Ellicott City, Maryland

As flooding events in the United States grow in frequency and intensity, the use of technological advancements and applied science are increasingly necessary for effective flood monitoring and warning systems. The NASA DEVELOP Ellicott City Disasters II project investigated the use of machine learning for applications in flood risk detection to support the improvement of early warning systems. To strengthen the efforts of the Howard County Office of Emergency Management (OEM) in building a more robust flood monitoring system, the project improved the original statistical flood risk model, FLuME (Flood Learning Model Environment), programmed by the first DEVELOP term. The enhancements incorporated an additional six years of precipitation and soil moisture data from the North American Land Data Assimilation System (NLDAS), modeled using Aqua Advanced Microwave Scanning Radiometer for EOS and Tropical Rainfall Measuring Mission (TRMM) Microwave Imager. These Earth observations were supplemented by stream gauge data from the OEM and the US Geological Survey. The resultant flood risk model FLASH (Flood Learning Environment and Severity Assessment Hub) was trained to evaluate input variables and predict stage height in Ellicott City in real time. The addition of an advanced deep learning framework known as long short-term memory improved the model’s ability to capture relationships between variables. To assess the effectiveness of the new model, FLASH produced a model efficiency metric of 0.99, a significant improvement over the 0.85 value produced by the previous model. The project assisted the OEM in pursuing the integration of open data and NASA Earth observations into a threat matrix capable of informing near real-time decision making.

Disasters↗

Ellicott City Disasters II: Enhancing a Statistical Flood Risk Model to Continue Improving Early Warning Systems and Public Safety in Ellicott City, Maryland

As flooding events in the United States grow in frequency and intensity, the use of technological advancements and applied science are increasingly necessary for effective flood monitoring and warning systems. The NASA DEVELOP Ellicott City Disasters II project investigated the use of machine learning for applications in flood risk detection to support the improvement of early warning systems. To strengthen the efforts of the Howard County Office of Emergency Management (OEM) in building a more robust flood monitoring system, the project improved the original statistical flood risk model, FLuME (Flood Learning Model Environment), programmed by the first DEVELOP term The enhancements incorporated an additional six years of precipitation and soil moisture data from the North American Land Data Assimilation System (NLDAS), modeled using Aqua Advanced Microwave Scanning Radiometer for EOS and Tropical Rainfall Measuring Mission TRMM Microwave Imager. These Earth observations were supplemented by stream gauge data from the OEM and the US Geological Survey. The resultant flood risk model FLASH (Flood Learning Environment and Severity Assessment Hub) was trained to evaluate input variables and predict stage height in Ellicott City in real time. The addition of an advanced deep learning framework known as long short-term memory improved the model’s ability to capture relationships between variables. To assess the effectiveness of the new model, FLASH produced a model efficiency metric of 0.99, a significant improvement over the 0.85 value produced by the previous model. The project assisted the OEM in pursuing the integration of open data and NASA Earth observations into a threat matrix capable of informing near real-time decision making.

Disasters↗

CO2 exposure at pressure impacts metabolism and stress responses in the model sulfate-reducing bacterium Desulfovibrio vulgaris strain Hildenborough

Geologic carbon dioxide (CO2) sequestration drives physical and geochemical changes in deep subsurface environments that impact indigenous microbial activities. The combined effects of pressurized CO2 on a model sulfate-reducing microorganism, Desulfovibrio vulgaris, have been assessed using a suite of genomic and kinetic measurements. Novel high-pressure NMR time-series measurements using 13C-lactate were used to track D. vulgaris metabolism. We identified cessation of respiration at CO2 pressures of 10 bar, 25 bar, 50 bar, and 80 bar. Concurrent experiments using N2 as the pressurizing phase had no negative effect on microbial respiration, as inferred from reduction of sulfate to sulfide. Complementary pressurized batch incubations and fluorescence microscopy measurements supported NMR observations, and indicated that non-respiring cells were mostly viable at 50 bar CO2 for at least four hours, and at 80 bar CO2 for two hours. The fraction of dead cells increased rapidly after four hours at 80 bar CO2. Transcriptomic (RNA-Seq) measurements on mRNA transcripts from CO2-incubated biomass indicated that cells up-regulated the production of certain amino acids (leucine, isoleucine) following CO2 exposure at elevated pressures, likely as part of a general stress response. Evidence for other poorly understood stress responses were also identified within RNA-Seq data, suggesting that while pressurized CO2 severely limits the growth and respiration of D. vulgaris cells, biomass retains intact cell membranes at pressures up to 80 bar CO2. Together, these data show that geologic sequestration of CO2 may have significant impacts on rates of sulfate reduction in many deep subsurface environments where this metabolism is a key respiratory process.

Desulfovibrio↗

Simplified predictive models for CO 2 sequestration performance assessment

CO2 sequestration in deep saline formations is increasingly being considered as a viable strategy for the mitigation of greenhouse gas emissions from anthropogenic sources. In this context, detailed numerical simulation based models are routinely used to understand key processes and parameters affecting pressure propagation and buoyant plume migration following CO2 injection into the subsurface. As these models are data and computation intensive, the development of computationally-efficient alternatives to conventional numerical simulators has become an active area of research. Such simplified models can be valuable assets during preliminary CO2 injection project screening, serve as a key element of probabilistic system assessment modeling tools, and assist regulators in quickly evaluating geological storage projects. We present three strategies for the development and validation of simplified modeling approaches for CO2 sequestration in deep saline formations: (1) simplified physics-based modeling, (2) statisticallearning based modeling, and (3) reduced-order method based modeling. In the first category, a set of full-physics compositional simulations is used to develop correlations for dimensionless injectivity as a function of the slope of the CO2 fractional-flow curve, variance of layer permeability values, and the nature of vertical permeability arrangement. The same variables, along with a modified gravity number, can be used to develop a correlation for the total storage efficiency within the CO2 plume footprint. Furthermore, the dimensionless average pressure buildup after the onset of boundary effects can be correlated to dimensionless time, CO2 plume footprint, and storativity contrast between the reservoir and caprock. In the second category, statistical “proxy models” are developed using the simulation domain described previously with two approaches: (a) classical Box-Behnken experimental design with a quadratic response surface, and (b) maximin Latin Hypercube sampling (LHS) based design with a multidimensional kriging metamodel fit. For roughly the same number of simulations, the LHS-based metamodel yields a more robust predictive model, as verified by a k-fold cross-validation approach (with data split into training and test sets) as well by validation with an independent dataset. In the third category, a reduced-order modeling procedure is utilized that combines proper orthogonal decomposition (POD) for reducing problem dimensionality with trajectory-piecewise linearization (TPWL) in order to represent system response at new control settings from a limited number of training runs. Significant savings in computational time are observed with reasonable accuracy from the PODTPWL reduced-order model for both vertical and horizontal well problems – which could be important in the context of history matching, uncertainty quantification and optimization problems. The simplified physics and statistical learning based models are also validated using an uncertainty analysis framework. Reference cumulative distribution functions of key model outcomes (i.e., plume radius and reservoir pressure buildup) generated using a 97-run full-physics simulation are successfully validated against the CDF from 10,000 sample probabilistic simulations using the simplified models. The main contribution of this research project is the development and validation of a portfolio of simplified modeling approaches that will enable rapid feasibility and risk assessment for CO2 sequestration in deep saline formations.

54 ENVIRONMENTAL SCIENCES↗