Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geoscience”

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 235 records · Page 13

Tracking the spatial and temporal links between late Cenozoic extension and magmatism in the Death Valley region (California, USA)

The relationship between late Cenozoic magmatism and extension in the central Basin and Range province (western United States) is complex, necessitating high-precision geochronology to understand its spatiotemporal connections. In the Death Valley region (California), the lack of high-precision U-Pb zircon ages has limited our understanding of the timing of pluton formation and its links to regional extension. We present new high-precision chemical abrasion-isotope dilution-thermal ionization mass spectrometry 206 Pb/ 238 U zircon ages and trace element analyses for eight Death Valley plutons. Our findings reveal three distinct phases of intrusive magmatism: (1) emplacement of shallow rapakivi granites at 13.2 Ma, (2) construction of the mid-crustal Black Mountains intrusive complex at 11.3 Ma, and (3) late emplacement of shallow, compositionally diverse intrusions at 8.2 Ma. A gap in zircon crystallization between 10 Ma and 8.2 Ma coincides with exhumation of the Black Mountains and a transition from sill to dike emplacement. The dominance of rapakivi granites in the Death Valley region, which is rare among Cenozoic granitoids, is likely a result of rapid crustal extension that induces adiabatic decompression. A comparison of the timing of volcanism, plutonism, and tectonic events in Death Valley reveals that intrusive magmatism closely tracks the locus of extension, underscoring the plutonic record as a vital link for understanding regional tectonics and changes in plate boundary dynamics during this period.

58 GEOSCIENCES↗

A More La Niña–Like Response to Radiative Forcing after Flux Adjustment in CESM2

Abstract In response to greenhouse gas forcing, most coupled global climate models project the tropical Pacific SST trend toward an “El Niño–like” state, with a reduced zonal SST gradient and a weakened Walker circulation. However, observations over the last five decades reveal a trend toward a more “La Niña–like” state with a strengthening zonal SST gradient. Recent research indicates that the identified trend differences are unlikely to be entirely due to internal variability and probably result, at least in part, from systematic model biases. In this study, Community Earth System Model, version 2 (CESM2), is used to explore how mean-state biases within the model may influence its forced response to radiative forcing in the tropical Pacific. The results show that using flux adjustment to reduce the mean-state bias in CESM2 over the tropical regions results in a more La Niña–like trend pattern in the tropical Pacific, with a strengthening of the tropical Pacific zonal SST gradient and a relatively enhanced Walker circulation, as hypothesized to occur if the ocean thermostat mechanism is stronger than the atmospheric mechanisms which by themselves would weaken the Walker circulation. We also find that the historical strengthening of the tropical Pacific zonal gradient is transient but persists into the near term in a high-emissions future warming scenario. These results suggest the potential of flux adjustment as a method for developing alternative projections that represent a wider range of possible future tropical Pacific warming scenarios, especially for a better understanding of regional patterns of climate risk in the near term.

Zhuo, Jing-Yi [Lamont-Doherty Earth Observatory, C↗

Discrimination of Small Earthquakes and Buried Single-Fired Chemical Explosions at Local Distances (<150 km) in the Western United States from Comparison of Local Magnitude ( M L ) and Coda Duration Magnitude ( M C )

Seismologists distinguish underground nuclear explosions from more commonly occurring earthquakes using moment tensor inversion, high-frequency P/S amplitude ratios, m b :M s comparisons, and P–pP differential travel times. These methods are generally successful for large seismic events (M >3-4) well recorded at regional-to-teleseismic distances (>150 km); however, it is unclear whether they can be modified to work for small events (M <3) well recorded only at local distances (<150 km). Here, we evaluate a recently proposed, local-distance seismic source discriminant—the difference between local magnitude (M L ) and coda duration magnitude (M C )—using seismograms of earthquakes and buried, single-fired chemical explosions recorded in three regions of the western United States. The quantity M L –M C was previously found to be sensitive to source depth, effectively discriminating mine blasts, induced earthquakes, and very shallow tectonic earthquakes from deeper crustal earthquakes. In this study, we report the first evaluation of M L –M C as a depth discriminant using data from buried, single-fired explosions that, unlike the seismic sources studied earlier, are good analogs for underground nuclear explosions. We find that even when using generic, uncalibrated methods of assigning magnitudes, M L –M C separates single-fired explosions and earthquakes. The area under the receiver operating characteristic curve is 0.92 for 19 explosions and 14 earthquakes in Washington, 0.90 for 22 explosions and 90 earthquakes in Wyoming, and 0.99 for three explosions and 149 earthquakes in Nevada. Finally, M L :M C comparisons have the potential to enhance discrimination based on high-frequency P/S amplitudes ratios—which perform less well at local than regional distances—because the two metrics have complementary sensitivities.

58 GEOSCIENCES↗

Relocation of the 8 September 2023 High Atlas, Morocco, Earthquake Aftershock Sequence

The earthquake that occurred on 8 September 2023, with a magnitude of 6.8, was the most destructive earthquake event in Morocco in the past decade. This earthquake took place in the Al Haouz region, located in the western part of the High Atlas Mountain range. To better understand what caused and triggered this earthquake, the earthquake catalogs including P and S arrival times were collected from the Moroccan seismic network and combined with regional data from the International Seismological Centre. The mainshock and aftershocks were relocated by using iLoc, a state-of-the-art single-event location algorithm, and then by the multiple event location double-difference algorithm, hypoDD. The improved earthquake relocations using iLoc and the double-difference methods provide sharper lineation of seismicity and agree well with tomographic images of the earthquake zone. Finally, the seismicity distribution and the focal mechanism of the mainshock indicate that the earthquake sequence has occurred along the South Atlas fault system.

58 GEOSCIENCES↗

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗

SCEC/USGS Community Stress-Drop Validation Study: How Spectral Fitting Approaches Influence Measured Source Parameters

Spectral source parameters used to estimate an earthquake’s stress drop (⁠Δσ⁠) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress‐Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low‐frequency displacement (⁠Ω 0 ⁠) and corner frequency (⁠ƒ c ⁠) by fitting a smooth function to benchmark displacement spectra. Participants used linear‐ or log‐sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate ⁠Ω 0 ⁠, ƒ c ⁠, and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in and can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade‐off between Ω 0 and ƒ c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω 0 and somewhat higher ƒ c compared to those assuming a Brune model, although resulting Δσ estimates are similar. Finally, when compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ⁠.

58 GEOSCIENCES↗

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗

Iron isotope fractionation between solid and liquid metal in the Fe-P±Ni system: Experimental constraints and implications for meteorites

Iron meteorites record a range of Fe isotope compositions that hold valuable information regarding the evolution of their parent bodies. Interpreting this isotopic variability, however, requires experimental constraints on the equilibrium isotope fractionation between phases. It is thought that the cores of many iron meteorite parent bodies experienced fractional crystallization, during which crystallization of solid iron-nickel occurs from an increasingly non-metal-rich liquid alloy. Phosphorus is one component of this alloy, and this study provides the first constraints on Fe-isotope fractionation between solid and liquid alloys in the Fe-Ni-P system. Experiments comprising Fe and P show a clear enrichment in the light isotopes of Fe in the liquid phase, which increases with the amount of phosphorus. Nickel-bearing samples are offset from the trend defined by Ni-free experiments, which is accounted for by the change in the solid alloy phase from a body-centered cubic to face-centered cubic structure upon the addition of Ni. The increasing light isotope enrichment of the liquid with increasing P content suggests interstitial solution of P, which is known to lengthen Fe-Fe bonds in Fe-P liquids (Waseda and Shiraishi 1977). Results suggest a negligible effect of P on Fe isotope fractionation during planetesimal core crystallization. Iron isotopes may, however, prove useful for identifying the petrogenesis of schreibersite in pallasites and iron meteorites.

58 GEOSCIENCES↗

CO₂ Storage prospeCtive Resource Estimation Excel aNalysis (CO₂-SCREEN) User’s Manual

This user’s manual guides the use of the National Energy Technology Laboratory’s (NETL) CO₂ Storage prospeCtive Resource Estimation Excel aNalysis (CO₂-SCREEN) tool, which was developed to aid users screening geologic formations for prospective CO₂ storage resources. This manual is specific to the CO₂-SCREEN 4.0 version which is based in Python. CO₂-SCREEN applies U.S. Department of Energy (DOE) methods and equations for estimating prospective CO₂ storage resources for saline formations, shale formations, and residual oil zones (ROZ). CO₂-SCREEN was developed to be substantive and user-friendly and provide a consistent method for calculating prospective CO₂ storage resources. CO₂-SCREEN uses a Java based graphical user interface for data inputs and uses Python to calculate prospective CO₂ storage resources.

54 ENVIRONMENTAL SCIENCES↗

Report for Department of State (DoS) V-Fund Project "Tying Moment Tensor Solutions to Explosive Yield"

The goal of this study was to use U.S. nuclear explosions with known source parameters (yield, depth, shot point material and/or parameters) to determine moment-derived yield estimates. This work has been accomplished by performing full moment tensor solutions using regional network data from the LLNL network, and other regional broadband stations. As part of this study, we have calculated solutions for 130 U.S. nuclear explosions and 12 additional chemical explosions. We then take several approaches to using moment tensors to estimate yield, considering both the full and isotropic moment tensors, doing straight regression analysis on the whole dataset, then successively refining the calibration with additional information about material and overburden. We have also tried a completely new approach of using the seismic moment to help estimate the radiated seismic energy and tying this to yield through a seismic efficiency. Results appear to be promising, but more work might be required to make this more useful in an operational sense.

58 GEOSCIENCES↗

Interpretation of Mineralogical Diagenesis for Assessment of Radionuclide Transport at Pahute Mesa, Nevada National Security Site

Pahute Mesa (PM) in southern Nevada was one of the four principal underground nuclear testing areas at the Nevada National Security Site (NNSS), known as the “Nevada Test Site” until 2010 (Figure 1). Eighty-two underground nuclear tests were conducted at PM in northwestern corner of NNSS (USDOE, 2015). The PM area of underground nuclear testing is located near the center of the Pahute Mesa and Oasis Valley (PM-OV) groundwater basin, which discharges to the southwest near Oasis Valley (OV) north of Beatty, Nevada (Fenelon et al. 2016; Jackson et al., 2020). Most of the PM tests were located at depths near or below the regional water table and, therefore, present potential sources of groundwater contamination from release of test-related radionuclides.

58 GEOSCIENCES↗

Project "Cropland Abandonment Mapping" (Final Report)

Cropland abandonment is a common type of land-use change across the globe and can be caused by a range of social, economic, and environmental factors. Cropland abandonment has wide-ranging effects on the environment, especially on biodiversity and carbon storage. Despite the importance of cropland monitoring, cropland abandonment is difficult to map and not routinely monitored. In this project, we propose to investigate cropland abandonment using remote sensing technique for the application of carbon storage in Eurasia, where experienced extensive cropland abandonment because of drastic social-economic changes since the 1990s.

54 ENVIRONMENTAL SCIENCES↗

Jaynesian Analysis of Environmental Chemistry: Systems Model Component Integration via the Arctic Aquatic Carbon Cycle

Focal Areas: All three are touched upon from the OBER call. We will consider formal logical inference and entropic information analysis for sparse but intensive chemical measurements, as optimal drivers for managing process complexity. As a centerpiece for the discussion we take soil-to-sea organic structural evolution within the Arctic hydro-geochemosphere, assessing detrital negentropy then algorithmic greed to generate ranked, reduced mechanisms for system integration.

54 ENVIRONMENTAL SCIENCES↗

Marcellus Shale Energy and Environment Laboratory (MSEEL) (Final Report)

The objective of the Marcellus Shale Energy and Environment Laboratory (MSEEL) was to provide a long-term field site to develop and validate new knowledge and technology to improve recovery efficiency and minimize environmental implications of unconventional resource development. MSEEL initiated in October 2014 and completed in September 2021. Total project value was $\$29,765,067$, support from the US Department of energy totaled $\$16,608,355$ with a cost share of $\$13,156,712$ primarily from Northeast Natural Energy. This report in a departure from previous reports summarizes the overarching results and outlines the approach taken. We cover two recent efforts in machine learning and reservoir characterization and simulation. Numerous quarterly reports, public presentations and numerous external publications cover specific results by subtopic and in detail. Publications are listed in the Appendix.

03 NATURAL GAS↗

LLNL 5-2428: Fracture Permeability and Seismic Slip Behavior

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameteruncertainty and better predict and mitigate seismic hazard at Utah FORGE and future EGS sites.We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, herein referred to as “THMC+E” models. This capability will enable improvedengineering decisions at Utah FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. Utah FORGE management and future EGS operators may employ results of THMC+E models for decision making purposes throughout the lifetime of the field operation. For example, before production well installation, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite ofpossible well locations and flow rates, thus enabling evaluation of optimal production well placement. Such efforts will be conducted throughout the lifetime of the project, whereby additional laboratory experiments will constrain key model parameters and machine learning (ML) will reduce the size of the parameter space and the associated uncertainty. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗