Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geophysical constraints”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Geophysical Constraints on Decarbonized Systems—Building Spatio-Temporal Uncertainties into Future Electricity Grid Planning

Future electricity grids will be characterized by the high penetration of renewables to support the decarbonization process. Yet, this transition will further expose grids to a broad spectrum of geophysical forces, such as weather and climate, land and space, minerals and materials. Understanding their role is thus important to support the grid decarbonization process. Here, we synthesize the current body of knowledge on the relationship between geophysical constraints and electricity grid planning. Additionally, we show that there have been promising advances in the data, methods, and modelling tools needed to incorporate the effect of geophysical constraints on demand, resource availability, and grid operations. There are, however, multiple research avenues that deserve further attention. More system-specific and finer-scale analyses are necessary to better understand how spatio-temporal variability---and associated uncertainty---in geophysical forces affect grid planning. Moreover, we need a broader focus on the multi-sectoral implications of grid decarbonization efforts, so as to limit the risk of unintended consequences. Importantly, all these efforts are challenged by the computational requirements of existing power system models, which often limit our ability to characterize uncertainty, scale analyses across larger domains, or study the multi-sector implications of grid expansion decisions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES↗
DOE OSTI2024

Feasibility of peak temperature targets in light of institutional constraints

Despite faster-than-expected progress in clean energy technology deployment, global annual CO 2 emissions have increased from 2020 to 2023. The feasibility of limiting warming to 1.5 °C is therefore questioned. Here we present a model intercomparison study that accounts for emissions trends until 2023 and compares cost-effective scenarios to alternative scenarios with institutional, geophysical and technological feasibility constraints and enablers informed by previous literature. Our results show that the most ambitious mitigation trajectories with updated climate information still manage to limit peak warming to below 1.6 °C (‘low overshoot’) with around 50% likelihood. However, feasibility constraints, especially in the institutional dimension, decrease this maximum likelihood considerably to 5–45%. Accelerated energy demand transformation can reduce costs for staying below 2 °C but have only a limited impact on further increasing the likelihood of limiting warming to 1.6 °C. Our study helps to establish a new benchmark of mitigation scenarios that goes beyond the dominant cost-effective scenario design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Inference of Multichannel r -process Element Enrichment in the Milky Way Using Binary Neutron Star Merger Observations

Observations of GW170817 strongly suggest that binary neutron star (BNS) mergers produce rapid neutron-capture nucleosynthesis ( r -process) elements. However, it remains an open question whether these mergers can account for all the r -process element enrichment in the Milky Way’s history. Here, we constrain the contributions of the BNS channel using astrophysical neutron star observations. The rate and mass distributions are constrained by LIGO/Virgo/Kagra through the latest catalog GWTC-3, the neutron star equation of state by gravitational-wave, radio, and X-ray observations, and the delay time distribution by short gamma-ray burst (GRB) host galaxy associations. We present a Bayesian framework to consistently combine these observations with abundance information to quantify the contribution and uncertainties of single and multiple astrophysical enrichment sources, and obtain a distribution of per-event BNS r -process element yields consistent with geophysical and astrophysical abundance constraints. We then adopt a Galactic chemical evolution model assuming an instantaneous and fixed amount of Fe enrichment from core-collapse supernovae, and show that BNS-only enrichment scenarios remain inconsistent with the observed r-process abundance trend of disk stars in the Galaxy even with the uncertainties in BNS merger observations. Using stellar abundance observations instead of the short GRB constraints, we can infer a shorter BNS delay time distribution with power-law index α ≤ −2.0 and minimum delay time ${t}_{{\rm{\min }}}\leqslant 40$ Myr at 90% confidence, consistent with detailed Galactic chemical evolution models. Such delay times are in tension with those predicted by standard BNS formation models. Alternatively, we confirm that a two-channel scenario, in which the second channel tracks the star formation history without significant delay, can account for both Galactic stellar and short GRB observations. We estimate that 45%–90% of the r -process abundance in the Milky Way today would have been produced by this star formation-tracking channel, rather than BNS mergers with significant delay times.

gravitational wave astronomy↗
DOE OSTI2025

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↗
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2025

Basin & Range Investigation for Developing Geothermal Energy

Hidden geothermal systems represent a potentially prolific energy resource that could support critical U.S. public and government energy priorities. Basin and Range Investigations for Developing Geothermal Energy (BRIDGE) addressed some the challenges associated with hidden system exploration by prioritizing cost-effective exploration early on through strategic workflow and informed decision-making that mitigates early risk and shifts resources to later exploration stages (e.g., drilling). Sandia National Laboratories partnered with U.S. Navy Geothermal Office, Geologic Geothermal Group, and independent consultants, with additional collaboration with U.S. Geological Survey and private industry. The primary tool of the BRIDGE project was to deploy a regional-scale airborne electromagnetic method to investigate the shallow resistivity structure in areas with high prospectivity. This was followed up at several prospects by a multidisciplinary exploration approach, including additional geologic, geophysical and geochemical studies. A central tenet to the BRIDGE methodology is that zones of low resistivity frequently occur over geothermal systems in the Basin and Range, and when paired with other data constraints, imaging these zones can enable discovery of these systems. In addition to exploring greenfield areas (i.e., Grover Point), the BRIDGE project also flew HTEM resistivity surveys over known geothermal systems including those with established power plants (Don A. Campbell and Salt Wells) and prospects that are known to the literature but remain undeveloped, at least in part, due to a lack of understanding on the location of their producible reservoirs. BRIDGE produced a comprehensive set of data from prospects identified in the Nevada Play Fairway Analysis along with conceptual models for top ranking prospects, wherein all of the observations are used to inform an interpreted model of the system. These models present a range of possible system parameters such as temperature and size, and they are further informed by system analogues in the Basin and Range province and elsewhere. The results of this work leave space for further exploration that may now occur at prospects ‘down the list’ rather than distribution exploration resources evenly across all prospects.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models

Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.

17 WIND ENERGY↗
DOE OSTI2026

A ModEx Framework for Watershed Subsurface Investigation With Limited Geophysical Data Using Machine Learning and Hydrologic Modeling

Abstract Subsurface heterogeneity influences watershed hydrology strongly but remains difficult to characterize at catchment scales with sparse and costly field data. Geophysical surveys such as electromagnetic induction (EMI) provide local spatial subsurface images yet scaling them to watershed scales and converting EMI‐derived resistivity into hydraulic properties remains a challenge. We present a Model–Experiment (ModEx) framework that integrates limited EMI data with machine learning (ML) and hydrologic modeling to improve process representation and guide field investigations. Sparse EMI surveys were scaled to the catchment scale using a Random Forest model, and the resulting resistivity fields were combined with nearby borehole constraints to parameterize a hydrologic model. The EMI‐informed hydrological simulations improved predictions of streamflow sustained by subsurface flow and shallow saturation patterns. By combining EMI data and ML with hydrologic modeling, the ModEx framework guides future subsurface surveys, providing a transferable and efficient strategy for data–model integration across diverse watersheds. Plain Language Summary Mapping the underground network of soil and rock that controls water is essential for predicting floods and droughts, but seeing underground is difficult and expensive. We cannot drill everywhere, so scientists use geophysical tools to scan broad areas. There are two key challenges: these geophysical scans are often sparse across the whole watershed, and the geophysical data is hard to translate into water‐related properties. We used artificial intelligence to solve these problems. We taught a computer to find patterns linking the limited geophysical data to the land surface properties. This allowed it to fill in the gaps and create a complete, useful subsurface map for the entire watershed. This new map improves hydrologic simulations, leading to more accurate predictions of water movement in the watershed. It also helps scientists build better models with less data and generates a priority map showing where to measure next, making future investigations more efficient. Key Points Limited EMI scaled with ML improves catchment‐scale subsurface parameterization for hydrologic models The framework integrates hydrologic modeling with limited geophysical data to support subsurface investigation design ModEx framework offers a transferable data–model integration strategy that quantifies and reduces uncertainty guiding watershed studies

Chen, Hang↗
DOE OSTI2025

Physics or Knob‐Tuning? Tropical Anvil Peak Is Captured by GCMs

Tropical anvil clouds peak near 200 hPa and significantly impact Earth's climate, yet its physical realism in coarse‐resolution General Circulation Models (GCMs) remains debated. We examine anvil cloud formation by performing simulations with a GCM with a hierarchy of cloud fraction schemes ranging from a complex prognostic Tiedtke scheme to a simple binary scheme. All schemes consistently reproduce the anvil peak. The robust anvil peak arises because extremely cold temperatures at the upper troposphere facilitate frequent saturation events, producing clouds that disproportionately influence mean cloud fraction. Sensitivity experiments with enhanced evaporation of cloud condensate unexpectedly show increased anvil coverage, highlighting how slight evaporative moistening reinforces local saturation in cold upper‐tropospheric conditions. These results demonstrate that the tropical anvil cloud peak emerges from fundamental thermodynamic constraints, rather than specific cloud fraction parameterization choices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Volatiles and Redox Along the East African Rift

Abstract The upper mantle under the Afar Depression in the East African Rift displays some of the slowest seismic wave speeds observed globally. Despite the extreme nature of the geophysical anomaly, lavas that erupted along the East African Rift record modest thermal anomalies. We present measurements of major elements, H 2 O, S, and CO 2 , and Fe 3+ /ΣFe and S 6+ /ΣS in submarine glasses from the Gulf of Aden seafloor spreading center and olivine‐, plagioclase‐, and pyroxene‐hosted melt inclusions from Erta Ale volcano in the Afar Depression. We combine these measurements with literature data to place constraints on the temperature, H 2 O, andfO 2 of the mantle sources of these lavas as well as the initial and final pressures of melting. The Afar mantle plume is C/FOZO/PHEM in isotopic composition, and we suggest that this mantle component is damp, with 852 ± 167 ppm H 2 O, not elevated infO 2 compared to the depleted MORB mantle, and has temperatures of ∼1401–1458°C. This is similar infO 2 and H 2 O to the estimates of C/FOZO/PHEM in other locations. Using the moderate H 2 O contents of the mantle together with the moderate thermal anomaly, we find that melting begins at around 93 km depth and ceases at around 63 km depth under the Afar Depression and at around 37 km depth under the Gulf of Aden, and that ∼1%–29% partial melts of the mantle can be generated under these conditions. We speculate that the presence of melt, and not elevated temperatures or high H 2 O contents, are the cause for the prominent geophysical anomaly observed in this region.

Geochemistry & Geophysics↗
DOE OSTI2024

ML-Based Rock Properties and Seismic Volume Enhancement

This project aims to improve field-scale Carbon Capture and Storage (CCS) assessments by enhancing petrophysical and geophysical log predictions through machine learning and neural networks. In our work during EY23, we applied Conditional Variational Autoencoders (CVAEs) to predict compressional velocity (Vp) and assess CO2 saturation levels in geological formations at the Illinois Basin Decatur Project (IBDP). In another task, we improved full-waveform inversion (FWI) methods with machine-learning approaches using lithological constraints. Full-waveform inversion (FWI) of seismic data estimates the elastic properties of subsurface rocks with high spatial resolution.

Nathanail, Athanasios↗
DOE OSTI2026

Passive Microwave Radiometry and Active Radar Sounding as Complementary Tools for Geophysical Investigations of Icy Ocean Worlds

Juno Microwave Radiometer (MWR) observations of Europa and Ganymede offer critical insights into the icy shells of these moons ahead of NASA's Europa Clipper and ESA's JUpiter ICy moons Explorer (JUICE) missions. Both missions are equipped with active radar sounders designed to address key unknowns such as ice shell thickness, thermal state, and composition. In this study, we explore how passive microwave radiometry and active radar sounding can constrain ice shell properties, focusing on Europa. Using modeled microwave brightness temperature observations at 0.6 and 1.2 GHz alongside simulated radar attenuation rate observations, we show that each instrument can independently produce robust ice shell thickness constraints under idealized conditions. We then relax these assumptions, quantifying how uncertainties from non-ideal properties—including convective layers, freezing-point depression, and chloride-doped ice—affect thickness estimates. Finally, we demonstrate how combining observations from these complementary techniques breaks degeneracies between ice shell properties, enabling more robust constraints than either method alone. This approach will maximize the science return of Europa Clipper and JUICE, advancing our understanding of the thermophysical structure and habitability of icy ocean worlds.

58 GEOSCIENCES↗
DOE OSTI2024

High Sodium Solubility in Magnesiowüstite in Iron‐Rich Deep Lower Mantle

Abstract (Mg,Fe)O ferropericlase‐magnesiowüstite has been proposed to host the majority of Earth's sodium, but the mechanism and capacity for incorporating the alkali cation remain unclear. In this work, experiments in the laser‐heated diamond anvil cell and first‐principles calculations determine the solubility of sodium and favorability of sodium incorporation in iron‐rich magnesiowüstite relative to (Mg,Fe)SiO 3 bridgmanite. Reaction of Mg/(Mg + Fe) (Mg#) 55 and 28 olivine with NaCl at 33–128 GPa and 1600–3000 K produces iron‐rich magnesiowüstite containing several percent sodium, while iron‐rich bridgmanite contains little to no detectable sodium. In sodium‐saturated magnesiowüstite, sodium number [Na/(Na + Mg + Fe)] is 2–5 atomic percent at pressures below 60 GPa and drastically increases to 10–20 atomic percent at deep lower mantle pressures. For these two compositions, there is no significant dependence of the results on Mg#. Our calculations not only show consistent results with experiments but further indicate that such an increase in solubility and partitioning of Na into magnesiowüstite is driven by the spin transition in iron. These results provide fundamental constraints on the crystal chemistry of sodium at lower‐mantle conditions. If the sodium capacity of (Mg,Fe)O is not strongly dependent on Mg#, (Mg,Fe)O in the lower mantle may have the capacity to store the entire sodium budget of the Earth.

Geochemistry & Geophysics↗
DOE OSTI2025

Impact of Salinity on Ground Ice Distribution Across an Arctic Coastal Polygonal Tundra Environment

The heterogeneous distribution of ground ice in the Arctic is a key driver of uneven ground subsidence as permafrost thaws, significantly impacting infrastructure and surface/subsurface hydrology. These topographic and hydrological changes contribute to major uncertainties in energy and carbon fluxes and storage in a warming Arctic. This study aims to improve our understanding of the controls on ground ice and organic matter distribution within the top 3 m of permafrost in coastal polygonal tundra near Utqiagvik, Alaska. To this end, we apply a neural network approach to bulk density distributions derived from nondestructive X-ray tomography of soil cores, trained with laboratory analyses, to improve the resolution and spatial coverage of estimates of dry bulk density, ice content, and organic matter content. In addition, we use capacitively coupled geophysical imaging to map soil electrical conductivity and salinity variations. The results show that sedimentary deposits from ocean transgressions, along with subsequent ice wedge polygon geomorphological processes, jointly influence the distribution of ice content at various scales. The impact of the latter decreases with depth, whereas the influence of salinity and sedimentary history increases. Although the controls on the distribution of soil organic matter content (g/cm 3 ) remain unclear, the pronounced heterogeneity in bulk density strongly influences its calculation from laboratory mass fraction measurements (g/g). From a methodological perspective, the interdependencies among soil components and the need for increased data coverage underscore the value of high-resolution density measurements, such as using X-ray tomography. Overall, this study emphasizes the importance of considering salinity constraints on ice content distribution in coastal permafrost regions. The results are expected to aid in the development of data products and process representations in geomorphological and ecosystem models.

Dafflon, Baptiste [Lawrence Berkeley National Labo↗
DOE OSTI2024

Wasatch Fault Structure from Machine Learning Arrival Times and High-Precision Earthquake Locations

Abstract On 18 March 2020, a magnitude 5.7 earthquake hit the Salt Lake valley in the state of Utah, United States. Using a dense geophone deployment and machine learning (ML), an additional several thousand events were detected and located. Currently, both the mainshock and the majority of the aftershocks are suspected to have occurred on or near a deeper portion of the Salt Lake segment of the Wasatch fault—part of a large range-bounding fault system thought to be capable of generating an Mw 7.2 earthquake. However, a small subset of aftershocks may have occurred on a portion of the more steeply, eastward dipping, and poorly understood West Valley fault. Unfortunately, the catalog locations and lack of focal mechanisms for this subset of aftershocks provide only a crude constraint on the true fault structure. To better illuminate fault structure, we relocate the ML-generated catalog with a range of magnitudes from −2 to 4.6, using: (1) NonLinLoc, a nonlinear location algorithm, (2) source-specific station terms, and (3) waveform coherence. We further compute first-motion focal mechanisms for 68 events. Results of the relocation suggest a simpler, minimally listric Wasatch fault geometry, contrary to what has been previously proposed. We also find that analysis of the focal mechanisms and waveform similarity indicates minimal event similarity throughout the Magna sequence, suggesting a highly complex and heterogeneous rupture zone, as opposed to rupture on a single plane. These findings suggest an increased seismic hazard due to the overall shallowness of the earthquake sequence and highly varied rupture mechanisms.

Geochemistry & Geophysics↗