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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Deep Learning for Full Waveform Inversion of Elastic Active-Source Seismic Data to Estimate P-Wave Velocity Models

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2‐CISM2 Coupled Climate–Ice Sheet Model

Abstract During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Model complexity and choice of model approaches for practical simulations of CO 2 injection, migration, leakage and long-term fate

This report documents the accomplishments achieved during the project titled “Model complexity and choice of model approaches for practical simulations of CO 2 injection,migration, leakage and long-term fate” funded by the US Department of Energy, Office of Fossil Energy. The objective of the project was to investigate modeling approaches of various levels of complexity relevant to geologic carbon storage (GCS) modeling with the goal to establish guidelines on choice of modeling approach.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

An Inversion of Gravity and Topography for Mantle and Crustal Structure on Mars

Analysis of the gravity and topography of Mars presently provides our primary quantitative constraints on the internal structure of Mars. We present an inversion of the long-wavelength (harmonic degree less than or equal to 10) gravity and topography of Mars for lateral variations of mantle temperature and crustal thickness. Our formulation incorporates both viscous mantle flow (which most prior studies have neglected) and isostatically compensated density anomalies in the crust and lithosphere. Our nominal model has a 150-km-thick high-viscosity surface layer over an isoviscous mantle, with a core radius of 1840 km. It predicts lateral temperature variations of up to a few hundred degrees Kelvin relative to the mean mantle temperature, with high temperature under Tharsis and to a lesser extent under Elysium and cool temperatures elsewhere. Surprisingly, the model predicts crustal thinning beneath Tharsis. If correct, this implies that thinning of the crust by mantle shear stresses dominates over thickening of the crust by volcanism. The major impact basins (Hellas, Argyre, Isidis, Chryse, and Utopia) are regions of crustal thinning, as expected. Utopia is also predicted to be a region of hot mantle, which is hard to reconcile with the surface geology. An alternative model for Utopia treats it as a mascon basin. The Utopia gravity anomaly is consistent with the presence of a 1.2 to 1.6 km thick layer of uncompensated basalt, in good agreement with geologic arguments about the amount of volcanic fill in this area. The mantle thermal structure is the dominant contributor to the observed geoid in our inversion. The mantle also dominates the topography at the longest wavelengths, but shorter wavelengths (harmonic degrees greater than or equal to 4) are dominated by the crustal structure. Because of the uncertainty about the appropriate numerical values for some of the model's input parameters, we have examined the sensitivity of the model results to the planetary structural model (core radius and core and mantle densities), the mantle's viscosity stratification, and the mean crustal thickness. The model results are insensitive to the specific thickness or viscosity contrast of the high-viscosity surface layer and to the mean crustal thickness in the range 25 to 100 km. Models with a large core radius or with an upper mantle low-viscosity zone require implausibly large lateral variations in mantle temperature.

Kiefer, Walter S.↗
DOE OSTI2020

Center for the Geologic Storage of CO 2 (GSCO2) (Final Report)

The Center for Geologic Storage of CO2 (GSCO2) focuses on basic-science approaches to recognized industry and technology challenges for commercial deployment of carbon dioxide (CO2) storage in clastic reservoirs. With this focus, the GSCO2 strives to generate new conceptual, mathematical, and numerical models applicable to geologic storage systems in specific and strategically identified research areas. Over the four-year research period (plus two-year no-cost extension), the GSCO2 developed its focus on one central research question: What is the mechanism of injection-induced microseismicity, and can we control and predict its occurrence?

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1975

Evolution of the moon - The 1974 model

Based on the geological discoveries of the Apollo and Luna missions, the evolution of the moon is divided into seven major stages, beginning at some time near the end of the formation of the solar system. These stages are: (1) the beginning, 4.6 billion years ago (b.y.a.); (2) melting of the outer shell, 4.6 to 4.4 b.y.a.; (3) cratering of the highlands, 4.4 to 4.1 b.y.a.; (4) formation of the large basins, 4.1 to 3.9 b.y.a.; (5) formation of the light-colored plains, 3.9 to 3.8 b.y.a.; (6) formation of the basaltic maria, 3.8 to about 3.0 b.y.a.; and (7) the quiet-crust stage, about 3.0 b.y.a. to the present. Each Apollo and Luna mission is reviewed, emphasizing the most important discoveries and experiments performed. The limits on this interpretative understanding of the moon and the outlined sequence of events are pointed out.

Schmitt, H. H.↗
DOE OSTI2015

Elastic-Waveform Inversion with Compressive Sensing for Sparse Seismic Data

Accurate velocity models of compressional- and shear-waves are essential for geothermal reservoir characterization and microseismic imaging. Elastic-waveform inversion of multi-component seismic data can provide high-resolution inversion results of subsurface geophysical properties. However, the method requires seismic data acquired using dense source and receiver arrays. In practice, seismic sources and/or geophones are often sparsely distributed on the surface and/or in a borehole, such as 3D vertical seismic profiling (VSP) surveys. We develop a novel elastic-waveform inversion method with compressive sensing for inversion of sparse seismic data. We employ an alternating-minimization algorithm to solve the optimization problem of our new waveform inversion method. We validate our new method using synthetic VSP data for a geophysical model built using geologic features found at the Raft River enhanced-geothermal-system (EGS) field. We apply our method to synthetic VSP data with a sparse source array and compare the results with those obtained with a dense source array. Our numerical results demonstrate that the velocity models produced with our new method using a sparse source array are almost as accurate as those obtained using a dense source array.

58 GEOSCIENCES↗
DOE OSTI2025

Discrete fracture network model benchmarks developed and applied in a DECOVALEX-2023 repository performance assessment study

This study presents newly developed benchmarks for modeling flow and transport within discrete fracture networks (DFNs) and useful methods for analyzing the results. The new benchmarks are designed to test modeling approaches for use in probabilistic performance assessment models of deep geologic repositories in fractured rock. The benchmarks simulate flow and transport through a 1 km 3 block of fractured rock. The first simulates migration of a short pulse of tracer through a simple network of four intersecting fractures. The second adds 1089 stochastically generated fractures. The third changes the pulse to a continuous point source. Evaluation of model performance relies on moment analysis and comparison of the results of different models. The expected nondimensional first moment of the conservative tracer for each benchmark is 1. The benchmarks were simulated by teams from Canada, Czechia, Germany, Korea, Sweden, Taiwan, and the United States as part of a DECOVALEX-2023 study (decovalex.org). The teams used various approaches, including explicit DFN modeling, DFN upscaling to an equivalent continuous porous medium (ECPM), and a combination of both methods. Transport mechanisms are modeled using either the advection-dispersion equation or particle tracking. Results demonstrate strong agreement among the models in breakthrough behavior up to the 75th percentile. Significant deviations in first moments and well-clustered outputs led to the identification of inaccuracies in several models. Such findings exemplify the benefit of exercising these benchmarks and using the presented methods to test DFN flow and transport models.

Benchmark↗
NASA NTRS1990

Satellite laser ranging and geological constraints on plate motion

Satellite laser ranging (SLR) observed baseline rates of change were measured and compared with rates determined from sea floor spreading rates and directions, and earth-quake solutions. With the number of years of observation now over six for many of the baselines, the inaccuracy of determining baseline rates of change has diminished so that in some cases it is less than a few mm per year. Thus, a direct comparison between baseline rates of change and rates of change established using geophysical information (called geological rates) is now feasible. In most cases, there is good agreement between the rates determined from SLR and geological rates, but in some cases there appear to be discrepancies. These discrepancies involve many of the data for which one end of the baseline is either Quincy (California), Huahine (French Polynesia), or Simosato (Japan). A method for looking at the discrepancies for these SLR observatories has been devised which makes it possible to calculate the motion not modeled by the geologic information.

Harrison, C. G. A.↗
DOE OSTI2016

Additive effects of acidification and mineralogy on calcium isotopes in Triassic/Jurassic boundary limestones

The end-Triassic mass extinction coincided with a negative δ 13 C excursion, consistent with release of 13 C-depleted CO 2 from the Central Atlantic Magmatic Province. However, the amount of carbon released and its effects on ocean chemistry are poorly constrained. The co upled nature of the carbon and calcium cycles allows calcium isotopes to be used for constraining carbon cycle dynamics and vice versa. We present a high-resolution calcium isotope (δ 44/40 Ca) record from 100 m of marine limestone spanning the Triassic/Jurassic boundary in two stratigraphic sections from northern Italy. Immediately above the extinction horizon and the associated negative excursion in δ 13 C, δ 44/40 Ca decreases by ca. 0.8‰ in 20 m of section and then recovers to preexcursion values. Coupled numerical models of the geological carbon and calcium cycles demonstrate that this δ 44/40 Ca excursion is too large to be explained by changes to seawater δ 44/40 Ca alone, regardless of CO 2 injection volume and duration. Less than 20% of the δ 44/40 Ca excursion can be attributed to acidification. The remaining 80% likely reflects a higher proportion of aragonite in the original sediment, based largely on high concentrations of Sr in the samples. Our study demonstrates that coupled models of the carbon and calcium cycles have the potential to help distinguish contributions of primary seawater isotopic changes from local or diagenetic effects on the δ 44/40 Ca of carbonate sediments. Finally, differentiating between these effects is critical for constraining the impact of ocean acidification during the end-Triassic mass extinction, as well as for interpreting other environmental events in the geologic past.

58 GEOSCIENCES↗
DOE OSTI2018

Modeling dynamic fracture in granite under in situ conditions at high temperatures and pressures

Reproducing the in situ mechanical response of granite at a depth of a few kilometer under dynamic loading presents significant challenges. We have developed and calibrated a new thermo-mechanical model for granite using extensive experimental data from both quasi-static and dynamic test data. The model includes parameters which can be calibrated for a specific geologic setting. Using the model developed we have performed a series of 2D and 3D simulations to evaluate effects of in situ stress and temperature on the extent of damage surrounding underground explosions. We compare our results with available data for cavity size, plastic zone and peak velocity attenuation with range in granite.

58 GEOSCIENCES↗