Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6

Horizontal stress in planetary lithospheres from vertical processes

Understanding the stress states in a lithosphere is of fundamental importance for planetary geophysics. It is closely linked to the processes which form and modify tectonic features on the surface and reflects the behavior of the planet's interior, providing a constraint for the difficult problem of determining interior structure and processes. The tectonics on many extraterrestrial bodies (Moon, Mars, and most of the outer planet satellites) appears to be mostly vertical, and the horizontal stresses induced by vertical motions and loads are expected to dominate the deformation of their lithospheres. Herein, only changes are examined in the state of stress induced by processes such as sedimentary and volcanic deposition, erosional denudation, and changes in the thermal gradient that induce uplift or subsidence. This analysis is important both for evaluating stresses for specific regions in which the vertical stress history can be estimated, as well as for applying the proper loading conditions to global stress models. All references to lithosphere herein should be understood to refer to the elastic lithosphere, that layer which deforms elastically or brittlely when subjected to geologically scaled stresses.

Banerdt, W. B.↗

Constraints on the rheological structure of the mantle

Rheological models of the mantle are at present limited to radial symmetry, usually with homogeneous linearly viscous or viscoelastic incompressible layers. While such models are probably overly simple, they readily allow calculation of geophysical effects, such as post glacial rebound and related changes in the Earth's shape. They are also directly applicable to problems of global mantle flow and plate motions.

Oconnell, R. J.↗

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗

Horizontal stresses induced by vertical processes in planetary lithospheres

Understanding the state of stress in the elastic lithosphere is of fundamental importance for planetary geophysics, as it is the link between the observed geologic structures on the surface and the processes which form and modify these structures. As such, it can provide valuable constraints for the difficult problem of determining interior structure and processes. On the Earth, most large scale, organized deformation can be related to lateral tectonics associated with plate dynamics; however, the tectonics on many extraterrestrial bodies (such as the Moon, Mars, and most of the outer-planet satellites) appears to be primarily vertical in nature, and the horizontal stresses induced by vertical motions and loads are expected to dominate the deformation of their lithospheres. The largest stress contributions from vertical loading come from the flexure of the lithosphere, which induces both bending moments and membrane stresses. We are concerned here only with nonflexural changes in the state of stress induced by processes such as sedimentary and volcanic deposition, erosional denudation, and changes in the thermal gradient that induce uplift or subsidence. This analysis is important both for evaluating stresses for specific regions in which the vertical stress history can be estimated, as well as for applying the proper loading conditions to global stress models. It is also of interest for providing a reference state of stress for interpreting stress measurements in the crust of the Earth.

Banerdt, W. B.↗

Exploring the Clouds of Venus: Science Driven Aerobot Missions to Our Sister Planet

Venus is essential to our understanding of the evolution and habitability of Earth-size planets throughout the galaxy. The selection of the VERITAS, EnVision, and DAVINCI missions by NASA and ESA in June 2021 is an important step in advancing the science. However, addressing many of the most challenging questions will require in situ platforms that can operate in the Venus environment for extended periods in order to capture the full complexity of our sister planet. Aerobots are aerial vehicles that exploit buoyancy to achieve long-duration operation in the Venus cloud layer where environmental conditions are comparatively benign. Buoyancy control, explained in more detail in a companion paper at this conference, allows aerobots to change altitude with little energy expenditure enabling new scientific measurement opportunities. These include atmospheric chemistry, dynamics, geophysical measurements of the crust and interior and geological investigations enabled by high resolution surface imaging. One aspect to our approach to defining missions that fit within the resource constraints of competitive missions is keeping the scale small. Today’s science-driven appetite for sophisticated measurements and large volumes of data is driving size upwards but advances in technology can enable aerobots that can be delivered to Venus at manageable costs. The other aspect is supporting the aerobot at Venus with orbiters providing data relay, localization and synergistic science. The recently selected orbiters, equipped with low-cost proximity relay systems routinely used at Mars may obviate the need for dedicated orbiters thereby enabling Discovery mission candidates. Four aerobot mission concepts have been defined which fit comfortably within the current New Frontiers (NF) cost cap ($900M in $FY22). One of these concepts would also be a candidate for a Discovery mission if that cost cap ($500M in $FY19) were raised. Raising the NF cost cap would enable more capable aerobot missions combining both altitude control with synergistic orbital observations. Investigations of surface geology at high resolution with subcloud NIR nightside imaging and dropsondes on the dayside of Venus could also benefit from collaborations with foreign contributions.

"O'Rourke, Joseph"↗

The satellite altimeter as a platform for observation of the oceanic mesoscale

The use of the satellite radar altimeter as a platform to provide synoptic monitoring of the oceanic mesoscale is faced with two critical issues: removal of geoid error or contamination and election of optimum space/time sampling strategies. Long wavelength orbit determination errors are not critical problems for altimeter measurements of the basin scale circulation. Both issues are addressed within the constraints provided by orbital mechanics which dictates the laydown pattern of the satellite's groundtracks in space/time. Other issues which must be assessed are: adequate mission duration scales and the problems of geophysical noise sources and instrumental noise which degrade the effective alongtrack spatial resolution of the altimeter.

Mitchell, J. L.↗

Internal Dynamics and Crustal Evolution of Mars

The objective of this work is to improve understanding of the internal structure, crustal evolution, and thermal history of Mars by combining geophysical data analysis of topography, gravity and magnetics with results from analytical and computational modeling. Accomplishments thus far in this investigation include: (1) development of a new crustal thickness model that incorporates constraints from Mars meteorites, corrections for polar cap masses and other surface loads, Pratt isostasy, and core flattening; (2) determination of a refined estimate of crustal thickness of Mars from geoid/topography ratios (GTRs); (3) derivation of a preliminary estimate of the k(sub 2) gravitational Love number and a preliminary estimate of possible dissipation within Mars consistent with this value; and (4) an integrative analysis of the sequence of evolution of early Mars. During the remainder of this investigation we will: (1) extend models of degree-1 mantle convection from 2-D to 3-D; (2) investigate potential causal relationships and effects of major impacts on mantle plume formation, with primary application to Mars; (3) develop exploratory models to assess the convective stability of various Martian core states as relevant to the history of dynamo action; and (4) develop models of long-wavelength relaxation of crustal thickness anomalies to potentially explain the degree-1 structure of the Martian crust.

Zuber, Maria↗

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗

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↗

The use of baseline measurements and geophysical models for the estimation of crustal deformations and the terrestrial reference system

Four possible estimators are investigated for the monitoring of crustal deformations from a combination of repeated baseline length measurements and adopted geophysical models, particularly an absolute motion plate model. The first estimator is an extension of the familiar free adjustment. The next two are Bayesian type estimators, one weak and one strong. Finally, a weighted constraint estimator is presented. The properties of these four estimators are outlined and their physical interpretations discussed. A series of simulations are performed to test the four estimators and to determine whether or not to incorporate a plate model for the monitoring of deformations. The application of these estimations to the maintenance of a new conventional terrestrial reference system is discussed.

Bock, Y.↗

Theory and design of interferometric synthetic aperture radars

A derivation of the signal statistics, an optimal estimator of the interferometric phase, and the expression necessary to calculate the height-error budget are presented. These expressions are used to derive methods of optimizing the parameters of the interferometric synthetic aperture radar system (InSAR), and are then employed in a specific design example for a system to perform high-resolution global topographic mapping with a one-year mission lifetime, subject to current technological constraints. A Monte Carlo simulation of this InSAR system is performed to evaluate its performance for realistic topography. The results indicate that this system has the potential to satisfy the stringent accuracy and resolution requirements for geophysical use of global topographic data.

Rodriguez, E.↗

Rover-Based Instrumentation and Scientific Investigations During the 2012 Analog Field Test on Mauna Kea Volcano, Hawaii

Rover-based 2012 Moon and Mars Analog Mission Activities (MMAMA) were recently completed on Mauna Kea Volcano, Hawaii. Scientific investigations, scientific input, and operational constraints were tested in the context of existing project and protocols for the field activities designed to help NASA achieve the Vision for Space Exploration [1]. Several investigations were conducted by the rover mounted instruments to determine key geophysical and geochemical properties of the site, as well as capture the geological context of the area and the samples investigated. The rover traverse and associated science investigations were conducted over a three day period on the southeast flank of the Mauna Kea Volcano, Hawaii. The test area was at an elevation of ~11,500 feet and is known as "Apollo Valley" (Fig. 1). Here we report the integration and operation of the rover-mounted instruments, as well as the scientific investigations that were conducted.

Graham, L. D.↗

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↗

Leveraging the Gravity Field Spectrum for Icy Satellite Interior Structure Determination: The case of Europa with the Europa Clipper Mission

Understanding the interior structures of icy moons is pivotal for addressing their origins and habitability. We introduce an approach employing the gravity field spectrum as an additional constraint for the inversion of differentiated icy bodies' interior structures. After developing the general methodology, we apply it to Europa, utilizing the predicted measurement capability of NASA's Europa Clipper mission, and we prove its effectiveness in resolving key geophysical parameters. Notably, we show that using the gravity field spectrum in combination with the mass and moment of inertia of the body allows us to estimate, depending on the considered end-member interior structure, the hydrosphere thickness with 4–20 km uncertainty and reliably determine the seafloor maximum topographic range and elastic thickness to within 100–600 m and 5–15 km, respectively, together with the power–degree relationship of the seafloor topography. We also show that the proposed method allows us to determine the density of the silicate mantle and the radius of the core to within 0.25 g cc −1 and 50 km, respectively.

G. Cascioli↗

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↗

Early arrival waveform inversion using data uncertainties and matching filters with application to near-surface seismic refraction data

We develop an early arrival waveform inversion (EAWI) technique for high-resolution near-surface velocity estimation by iteratively updating the P-wave velocity model to minimize the difference between the observed and calculated seismic refraction data. Traditional EAWI uses a least-squares penalty function and an acoustic forward-modeling engine. Conventional least-squares error is sensitive to data with low signal-to-noise ratio (S/N) and iterations of EAWI stop at a local-minimum data misfit or at the preassigned maximum number of iterations. These stopping criteria can result in overfitting the data. In addition, fitting the elastic field data with an acoustic modeling engine can introduce artifacts in velocity estimation, especially in land data with significant elastic effects. To overcome these challenges, we develop a robust EAWI (REAWI) method by (1) incorporating the data uncertainties into the penalty function and (2) mitigating the elastic effects using a matching filter workflow. The data uncertainties are estimated from waveform reciprocal errors. When full-waveform reciprocity is not available, trace interpolation is applied. The proposed method prevents closely fitting data with low S/N, avoids overall overfitting by stopping the iterations when a normalized chi-square ([Formula: see text]) waveform misfit of one is achieved, and is less affected by elastic effects. Numerical examples and application to near-surface refraction data at a groundwater contamination site suggest that the final REAWI models are more accurate than the corresponding EAWI models, at the same level of misfit. This is the first known application of a matching filter workflow to real land data. The final REAWI models satisfy an appropriate misfit between the real data and predicted elastic P-wave data, making this approach in this respect equivalent to elastic waveform inversion. We also develop a method to analyze model constraint by examining the energy of the wavefield Fréchet derivative thereby avoiding the influence of the data residuals in traditional Fréchet kernels.

Geochemistry & Geophysics↗

On the coupling of fluid dynamics and electromagnetism at the top of the earth's core

A kinematic approach to short-term geomagnetism has recently been based upon pre-Maxwell frozen-flux electromagnetism. A complete dynamic theory requires coupling fluid dynamics to electromagnetism. A geophysically plausible simplifying assumption for the vertical vorticity balance, namely that the vertical Lorentz torque is negligible, is introduced and its consequences are developed. The simplified coupled magnetohydrodynamic system is shown to conserve a variety of magnetic and vorticity flux integrals. These provide constraints on eligible models for the geomagnetic main field, its secular variation, and the horizontal fluid motions at the top of the core, and so permit a number of tests of the underlying assumptions.

Benton, E. R.↗

Exploring for Superhot Geothermal Targets in Magmatic Settings: 2022 Field Campaign at Newberry Volcano

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN), supported by the U.S. Department of Energy (DOE) and Geothermica, a joint effort by EU member states and associated countries. The DEEPEN project will develop a customized approach to exploration for supercritical and superhot geothermal plays in magmatic systems, which will be applied to two demonstration sites. This paper summarizes field activities carried out at the U.S. demonstration site, Newberry Volcano in central Oregon. The objective of this work effort is to refine the subsurface model of Newberry Volcano, with special focus on deeper zones including the magmatic plumbing system and other key geologic elements. New data collection included gravity and wideband magnetotelluric (MT) surveys, as well as reinstallation of a seismic network. The National Renewable Energy Laboratory (NREL) and Enthalpion Energy LLC (Enthalpion) worked with the Deschutes National Forest Fort Rock District to use a low ground disturbance method of MT deployment to collect MT data inside the caldera and other restricted areas inside the National Volcanic Monument. This opened these areas to geophysical exploration for the first time in decades. Sites along and adjacent to the south rim of the volcano constituted the primary survey objectives. A team from Lawrence Berkeley National Laboratory (LBNL), the U.S. Geological Survey (USGS), and AltaRock also began the process of reinstalling the seismic network from the AltaRock enhanced geothermal system (EGS) demonstration in anticipation of further development activities at the site. The data ingestion, reduction, and analysis phase of the project is ongoing. We are currently processing the MT and gravity data and are developing a new, highly GPU-accelerated, 3D joint MT and gravity inversion to better localize the south rim/south flank conductive target and better understand its relationship to deep heat, fluid sources, and surface extrusive features. Joint inversions, which have not yet been undertaken at Newberry, will allow us to obtain constraints on the geologic model that cannot be determined from each method in isolation, improving our ability to image key geologic features at depth.

geophysics↗