Search NASA⌕ Search

SEARCH · Search NASA

Results for “Earth Model”

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 307 records · Page 17

A substantial role of soil erosion in the land carbon sink and its future changes

Realistic representation of land carbon sink in climate models is vital for predicting carbon climate feedbacks in a changing world. Although soil erosion that removes land organic carbon has increased substantially since the onset of agriculture, it is rarely included in the current generation of climate models. Using an Earth system model with soil erosion represented, we estimated that on average soil erosion displaces 5% of newly fixed land organic carbon downslope annually in the continental US. In the lower Mississippi river basin and the Cascades, the fraction can be as large as 40%. About 12% of the eroded organic carbon is eventually exported to inland waters, which is equal to 14% of the simulated net carbon gain by terrestrial ecosystems. By comparing the eroded organic carbon export to rivers with the particulate organic carbon export to oceans, we demonstrated that a large fraction of the carbon export to rivers could have been mineralized in inland waters. Importantly, with a direct comparison of eroded and exported soil organic carbon and land net carbon uptake, we found that Earth system models that ignore soil erosion likely offset the erosional carbon loss by increasing heterotrophic respiration implicitly. But as soil erosion and heterotrophic respiration respond differently to a warming climate, this unrealistic compensation would lead to biased predictions of future land carbon sink.

Tan, Zeli↗

Discrete global grid system-based flow routing datasets in the Amazon and Yukon basins

Abstract. Discrete global grid systems (DGGS) are emerging spatial data structures widely used to organize geospatial datasets across scales. While DGGS have found applications in various scientific disciplines, including atmospheric science and ecology, their integration into physically based hydrological models and Earth system models (ESMs) has been hindered by the lack of flow routing datasets based on DGGS. In response to this gap, this study pioneers the development of new flow routing datasets using icosahedral Snyder equal-area (ISEA) DGGS and a novel mesh-independent flow direction model. We present flow routing datasets for two large basins, the tropical Amazon River basin and the Arctic Yukon River basin. These datasets (1) facilitate the adoption of DGGS for hydrological models and (2) provide flow routing inputs for evaluation of DGGS-based flow routing in the Amazon and Yukon river basins. The data are available at https://doi.org/10.5281/zenodo.8377765 (Liao, 2023).

54 ENVIRONMENTAL SCIENCES↗

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

Onslow Beach Washover Fan Morphology from 2012 to 2015

Critical military training and testing on lands along the nation’s coastal and estuarine shorelines are increasingly placed at risk because of encroachment pressures in surrounding areas, impairments due to other anthropogenic disturbances, and changes in climate and sea level. The U.S. Department of Defense (DoD) intends to enhance and sustain its training and testing assets and also optimize its stewardship of natural resources through the development and application of an ecosystem-based management approach on DoD installations. To accomplish this goal, particularly for installations in estuarine/coastal environments, the Strategic Environmental Research and Development Program (SERDP) launched the Defense Coastal/Estuarine Research Program (DCERP) as a 10-year effort at Marine Corps Base Camp Lejeune (MCBCL) in North Carolina. The results of the second 5 years of the program (DCERP2) are presented in the DCERP2 Final Report.There were four overarching objectives of DCERP2. The first objective was to understand the effects of climate change impacts, including warming temperatures, variability in the hydrological cycle, storm events, and sea level rise on the coastal ecosystems at MCBCL from observations and measurements made over the 10-year program. The second objective was to understand the carbon cycle of the coastal and terrestrial ecosystems at MCBCL through a highly integrated sampling program. The third objective was to develop models, tools, and indicators to evaluate current and projected future ecosystem state changes and translate scientific findings into actionable information for installation managers. The last objective was to recommend adaptive management strategies to sustain ecosystem natural resources within the context of an active military installation.These grid files were created from data collected using a Riegl 3-D laser scanner. Millions of x, y, and z points from the F2 Washover Fan on Onslow Beach were processed using the MARS software package to define a bare-earth model. Grid files were imported into Golden Software Surfer version 19 for editing and exported from that application as .grd files. Files can be opened in Surfer or ArcMap, but the coordinate system information will be lost if opened in ArcMap. These data were collected to quantify the morphological changes that occur on a washover fan through time.

54 ENVIRONMENTAL SCIENCES↗

Spatial heterogeneity effects on land surface modeling of water and energy partitioning

Abstract. Understanding the influence of land surface heterogeneity on surface water and energy fluxes is crucial for modeling earth system variability and change. This study investigates the effects of four dominant heterogeneity sources on land surface modeling, including atmospheric forcing (ATM), soil properties (SOIL), land use and land cover (LULC), and topography (TOPO). Our analysis focused on their impacts on the partitioning of precipitation (P) into evapotranspiration (ET) and runoff (R), partitioning of net radiation into sensible heat and latent heat, and corresponding water and energy fluxes. An initial set of 16 experiments were performed over the continental US (CONUS) using the E3SM land model (ELMv1) with different combinations of heterogeneous and homogeneous datasets. The Sobol' total and first-order sensitivity indices were utilized to quantify the relative importance of the four heterogeneity sources. Sobol' total sensitivity index measures the total heterogeneity effects induced by a given heterogeneity source, consisting of the contribution from its own heterogeneity (i.e., the first-order index) and its interactions with other heterogeneity sources. ATM and LULC are the most dominant heterogeneity sources in determining spatial variability of water and energy partitioning, mainly contributed by their own heterogeneity and slightly contributed by their interactions with other heterogeneity sources. Their heterogeneity effects are complementary, both spatially and temporally. The overall impacts of SOIL and TOPO are negligible, except TOPO dominates the spatial variability of R/P across the transitional climate zone between the arid western and humid eastern CONUS. Accounting for more heterogeneity sources improves the simulated spatial variability of water and energy fluxes when compared with ERA5-Land reanalysis dataset. An additional set of 13 experiments identified the most critical components within each heterogeneity source, which are precipitation, temperature, and longwave radiation for ATM, soil texture, and soil color for SOIL and maximum fractional saturated area parameter for TOPO.

54 ENVIRONMENTAL SCIENCES↗

Discussion on “Saving Storage in Climate Ensembles: A Model-Based Stochastic Approach”

Abstract We thank the authors for this interesting paper that highlights important ideas and concepts for the future of climate model ensembles and their storage, as well as future uses of stochastic emulators. Stochastic emulators are particularly relevant because of the statistical nature of climate model ensembles, as discussed in previous work of the authors (Castruccio et al. in J Clim 32:8511–8522, 2019; Hu and Castruccio in J Clim 34:8409–8418, 2021). We thank the authors for sharing of some of their data with us in order to illustrate this discussion. In the following, in Sect. 1 we discuss alternative techniques currently used and studied, namely lossy compression and ideas emerging from the climate modeling community, that could feed the discussion on ensemble and storage. In that section, we also present numerical results of compression performed on the data shared by the authors. In Sect. 2, we discuss the current statistical model proposed by the authors and its context. We discuss other potential uses of stochastic emulators in climate and Earth modeling.

97 MATHEMATICS AND COMPUTING↗

Deep Learning-Assisted Real-Time Forward Modeling of Electromagnetic Logging in Complex Formations

Higher dimensional (i.e., 2-D and 3-D) modeling is indispensable to correctly evaluate the responses of electromagnetic logging tools in complex formation environments. However, limited by the high computational cost of rigorous modeling such as finite difference method and the finite element method, the real-time applications in the well logging industry primarily rely on the 1-D forward solver, which would result in erroneous formation evaluation for complex scenarios. As a result, aiming at realizing fast modeling for electromagnetic logging tools in complex formations, this paper proposes a general framework assisted by deep neural networks (DNNs). The framework consists of three modules: earth model classification, parameter extraction, and surrogate construction. Separate DNNs are trained and tested for different modules. The accuracy and efficiency of the DNN assisted fast modeling are validated by several experiments. Here, this study finds that the fast modeling assisted by DNNs is able to calculate the tool responses and reconstruct the subsurface formations in real-time.

97 MATHEMATICS AND COMPUTING↗

Task 5 Field Scale Stress Modeling Scripts and Data

Project data and Matlab scripts for Task 5 of the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage”. Project includes the static earth models, grids and other data used for simulation and calibration of the two case studies addressed in the project, referred to herein as the “FutureGen2” model (Illinois) and the “Perch” model (Otsego county, Michigan). Large Dataset, please contact EDXSupport@netl.doe.gov

FutureGen2↗

Comprehensive geoneutrino analysis with Borexino

This paper presents a comprehensive geoneutrino measurement using the Borexino detector, located at Laboratori Nazionali del Gran Sasso (LNGS) in Italy. The analysis is the result of 3262.74 days of data between December 2007 and April 2019. The paper describes improved analysis techniques and optimized data selection, which includes enlarged fiducial volume and sophisticated cosmogenic veto. The reported exposure of ( 1.29 ± 0.05 ) × 10 32 protons × year represents an increase by a factor of two over a previous Borexino analysis reported in 2015. By observing 52.6 - 8.6 + 9.4 ( stat ) - 2.1 + 2.7 ( sys ) geoneutrinos (68% interval) from U 238 and Th 232 , a geoneutrino signal of 47.0 - 7.7 + 8.4 ( stat ) - 1.9 + 2.4 ( sys ) TNU with - 17.2 + 18.3 % total precision was obtained. This result assumes the same Th/U mass ratio as found in chondritic CI meteorites but compatible results were found when contributions from U 238 and Th 232 were both fit as free parameters. Antineutrino background from reactors is fit unconstrained and found compatible with the expectations. The null-hypothesis of observing a geoneutrino signal from the mantle is excluded at a 99.0% C.L. when exploiting detailed knowledge of the local crust near the experimental site. Measured mantle signal of 21.2 - 9.0 + 9.5 ( stat ) - 0.9 + 1.1 ( sys ) TNU corresponds to the production of a radiogenic heat of 24.6 - 10.4 + 11.1 TW (68% interval) from U 238 and Th 232 in the mantle. Assuming 18% contribution of K 40 in the mantle and 8.1 - 1.4 + 1.9 TW of total radiogenic heat of the lithosphere, the Borexino estimate of the total radiogenic heat of the Earth is 38.2 - 12.7 + 13.6 TW , which corresponds to the convective Urey ratio of 0.78 - 0.28 + 0.41 . These values are compatible with different geological predictions, however there is a ~ 2.4 σ tension with those Earth models which predict the lowest concentration of heat-producing elements in the mantle. In addition, by constraining the number of expected reactor antineutrino events, the existence of a hypothetical georeactor at the center of the Earth having power greater than 2.4 TW is excluded at 95% C.L. Particular attention is given to the description of all analysis details which should be of interest for the next generation of geoneutrino measurements using liquid scintillator detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Integrated Modeling Report for CO 2 Storage with Enhanced Oil Recovery in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. MRCSP has been investigating various reservoir characterization and modeling technologies related to Carbon Capture, Utilization, and Storage (CCUS) in conjunction with carbon dioxide enhanced oil recovery (CO 2 -EOR) operations in multiple Silurian-age (Niagaran), oil-bearing, carbonate pinnacle reefs in northern Michigan, USA. This report provides a comprehensive discussion of reservoir modeling studies that were conducted for tracking oil production, forecasting CO 2 plume migration, and estimating associated storage in a number of reefs that were at different stages of their CO 2 -EOR life cycle. The modeling process for simulating oil production, CO 2 injection, and associated storage in these reefs entails two phases. The first phase, geologic framework modeling, integrates all pertinent geological and geophysical data (from logs, cores and seismic surveys) about reservoir structure, geometry, rock types, and property distributions (porosity, permeability, water saturation) into a 3-D distributed grid-based static earth model (SEM). The second phase, dynamic reservoir modeling, uses the SEM as a platform to simulate the movement of oil, gas, water, and CO 2 within the reservoir during primary hydrocarbon production, as well as during subsequent phases such as CO 2 -injection assisted EOR, plume migration, and associated storage. In addition, an assessment of coupled process effects is also carried out, where the impacts of geochemical and geo-mechanical processed induced by CO 2 injection are studied.

01 COAL, LIGNITE, AND PEAT↗

Sound Velocity of Hematite up to 13 GPa

We performed acoustic measurements on polycrystalline hematite using ultrasonic interferometric techniques in a Kawai-type multi-anvil apparatus. The elastic moduli and their pressure derivatives were determined by fitting the experimental data to third-order finite strain equations, yielding: 𝐾 𝑆,0 = 214.8 (20) GPa, 𝐺 0 = 92.0 (2) GPa, 𝐾$^{′}_{S,0}$ = 5.70 (1), and 𝐺$^{′}_{0}$ = 0.88 (4). Modeling of the velocity profiles of hematite along hot and cold slab geotherms shows that compressional (V P ) and shear (V S ) wave velocities are 22%–24% and 53%–55% lower, respectively, compared to the Preliminary Reference Earth Model at the top of the lower mantle. Furthermore, these findings likely offer a potential explanation for the low-velocity anomalies associated with slab graveyards at the base of the lower mantle.

36 MATERIALS SCIENCE↗

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗

The Hubble Space Telescope's Near-UV and Optical Transmission Spectrum of Earth as an Exoplanet

We observed the 2019 January total lunar eclipse with the Hubble Space Telescope’s STIS spectrograph to obtain the first near-UV (1700–3200 Å) observation of Earth as a transiting exoplanet. The observatories and instruments that will be able to perform transmission spectroscopy of exo-Earths are beginning to be planned, and characterizing the transmission spectrum of Earth is vital to ensuring that key spectral features (e.g., ozone, or O{sub 3}) are appropriately captured in mission concept studies. O{sub 3} is photochemically produced from O{sub 2}, a product of the dominant metabolism on Earth today, and it will be sought in future observations as critical evidence for life on exoplanets. Ground-based observations of lunar eclipses have provided the Earth’s transmission spectrum at optical and near-IR wavelengths, but the strongest O{sub 3} signatures are in the near-UV. We describe the observations and methods used to extract a transmission spectrum from Hubble lunar eclipse spectra, and identify spectral features of O{sub 3} and Rayleigh scattering in the 3000–5500 Å region in Earth’s transmission spectrum by comparing to Earth models that include refraction effects in the terrestrial atmosphere during a lunar eclipse. Our near-UV spectra are featureless, a consequence of missing the narrow time span during the eclipse when near-UV sunlight is not completely attenuated through Earth’s atmosphere due to extremely strong O{sub 3} absorption and when sunlight is transmitted to the lunar surface at altitudes where it passes through the O{sub 3} layer rather than above it.

79 ASTRONOMY AND ASTROPHYSICS↗

An Evaluation of the Large‐Scale Atmospheric Circulation and Its Variability in CESM2 and Other CMIP Models

Abstract The Community Earth System Model 2 (CESM2) is the latest Earth System Model developed by the National Center for Atmospheric Research in collaboration with the university community and is significantly advanced in most components compared to its predecessor (CESM1). Here, CESM2's representation of the large‐scale atmospheric circulation and its variability is assessed. Further context is providedthrough comparison to the CESM1 large ensemble and other models from the Coupled Model Intercomparison Project (CMIP5 and CMIP6). This includes an assessment of the representation of jet streams and storm tracks, stationary waves, the global divergent circulation, the annular modes, the North Atlantic Oscillation, and blocking. Compared to CESM1, CESM2 is substantially improved in the representation of the storm tracks, Northern Hemisphere (NH) stationary waves, NH winter blocking and the global divergent circulation. It ranks within the top 10% of CMIP class models in many of these features. Some features of the Southern Hemisphere (SH) circulation have degraded, such as the SH jet strength, stationary waves, and blocking, although the SH jet stream is placed at approximately the correct location. This analysis also highlights systematic deficiencies in these features across the new CMIP6 archive, such as the continued tendency for the SH jet stream to be placed too far equatorward, the North Atlantic westerlies to be too strong over Europe, the storm tracks as measured by low‐level meridional wind variance to be too weak and a lack of blocking in the North Atlantic sector.

Simpson, Isla R.↗

Data for publication Wan et al. (2020, JAMES), Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling

This dataset contains simple condensation mdoel output and CAM4 model output that were used to produce figures in the following paper: Wan, H., Woodward, C. S., Zhang, S., Vogl, C. J., Stinis, P., Gardner, D. J., Rasch, P. J., Zeng, X., Larson, V. E., and Singh, B. (2020): Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling, Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019MS001982

54 ENVIRONMENTAL SCIENCES↗

Amazon drought resilience – emerging results point to new empirical needs

Will climate change turn tropical forests from a carbon sink to a source (Pan et al., 2011)? Tropical forests cycle more carbon, water, and energy than any other biome (Bonan, 2008). The Amazon is the largest tropical forest, spanning one-third of South America, and any perturbation in its function has repercussions on the global climate. Mild and severe droughts in the Amazon are predicted to double and triple in area, respectively, by 2100, while the area under wet extremes will increase (Duffy et al., 2015). Our understanding of the role of plants' hydrological environments in determining their drought response, however, is limited (Chitra-Tarak et al., 2018, 2021). An important review by Costa et al. (2023; pp. 714–733) published in this issue of New Phytologist highlights that shallow-water table (WT) forests constitute c. 50% of the Amazon and may act as hydrological refugia during droughts, yet they have been neglected in Amazon forest research to date. Emerging results suggest that the shallow-WT forests that dominate the Amazon basin may increase in productivity under moderate droughts due to relief from hypoxia, suggesting a potential buffer from drought relative to deeper WT forests (Fig. 1). During severe droughts, however, shallow-WT forests may be vulnerable to collapse due to drought-intolerant traits. Addressing the underrepresentation of relatively drought-resilient shallow-WT forests in the Amazon's forest inventories may help resolve the much-debated incongruence in recent studies of the Amazon's drought resilience. Data from forest inventories (biased toward deep-WT forests) indicate that the Amazon's carbon sink is declining, partially because of drought impacts (Phillips et al., 2009; Brienen et al., 2015), whereas basin-wide satellite-based measures of gross primary productivity, which include shallow-WT forests, indicate varied regional responses to droughts (Saleska et al., 2007; Brando et al., 2010). Costa et al.'s (2023) descriptive insights of the intricate balance of geology, topography, hydrology, vegetation, and drought on ecosystem function, as well as their conceptual predictive framework, are useful for developing new empirical research in these understudied ecosystems and improving Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗