Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geophysical”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Salton Sea Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to Salton Sea Geothermal Field. It includes all input and output files used with the Geothermal Exploration Artificial Intelligence. Input and output files are sorted into three categories: raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs which are titled Radar, SWIR, Thermal, Geophysics, Geology, and Wells. The files are used with the Geothermal Exploration Artificial Intelligence for the Salton Sea Geothermal Site to identify indicators of blind geothermal systems. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Salton Sea Geothermal Site.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Anatomy of Continuous Mars SEIS and Pressure Data from Unsupervised Learning

The seismic noise recorded by the Interior Exploration using Seismic Investigations, Geodesy, and Heat Transport (InSight) seismometer (Seismic Experiment for Interior Structure [SEIS]) has a strong daily quasi-periodicity and numerous transient microevents, associated mostly with an active Martian environment with wind bursts, pressure drops, in addition to thermally induced lander and instrument cracks. That noise is far from the Earth’s microseismic noise. Quantifying the importance of nonstochasticity and identifying these microevents is mandatory for improving continuous data quality and noise analysis techniques, including autocorrelation. Cataloging these events has so far been made with specific algorithms and operator’s visual inspection. We investigate here the continuous data with an unsupervised deep-learning approach built on a deep scattering network. This leads to the successful detection and clustering of these microevents as well as better determination of daily cycles associated with changes in the intensity and color of the background noise. We first provide a description of our approach, and then present the learned clusters followed by a study of their origin and associated physical phenomena. We show that the clustering is robust over several Martian days, showing distinct types of glitches that repeat at a rate of several tens per sol with stable time differences. We show that the clustering and detection efficiency for pressure drops and glitches is comparable to or better than manual or targeted detection techniques proposed to date, noticeably with an unsupervised approach. Finally, here we discuss the origin of other clusters found, especially glitch sequences with stable time offsets that might generate artifacts in autocorrelation analyses. We conclude with presenting the potential of unsupervised learning for long-term space mission operations, in particular, for geophysical and environmental observatories.

58 GEOSCIENCES↗
DOE OSTI2020

Predicting the proximity to macroscopic failure using local strain populations from dynamic in situ X-ray tomography triaxial compression experiments on rocks

Predicting the proximity of large-scale dynamic failure is a critical concern in the engineering and geophysical sciences. Here we use evolving contractive, dilatational, and shear strain deformation preceding failure in dynamic X-ray tomography experiments to examine which strain components best predict the proximity to failure. We develop machine learning models to predict the proximity to failure using time series of three-dimensional local incremental strain tensor fields acquired in rock deformation experiments under stress conditions of the upper crust. Three-dimensional scans acquired in situ throughout triaxial compression experiments provide a distribution of density contrasts from which we estimate the three-dimensional incremental strain that accumulates between each scan acquisition. Training machine learning models on multiple experiments of six rock types provides suites of feature importance that indicate the predictive power of each feature. Comparing the average importance of groups of features that include information about each strain component quantifies the ability of the contractive, dilatational and shear strain to predict the proximity of macroscopic failure. A total of 24 models of four machine learning algorithms with six rock types indicate that 1) the dilatational strain provides the best predictive power of the strain components, and 2) the intermediate values (25th-75th percentile) of the strain population provide the best predictive power of the statistics of the strain populations. In addition, the success of the predictions of models trained on one rock type and tested on other rock types quantifies the similarities and differences of the precursory strain accumulation process in the six rock types. These similarities suggest the potential existence of a unified theory of brittle rock deformation for a range of rock types.

58 GEOSCIENCES↗
DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Illumination of Damage in Intact Rocks by Ultrasonic Transmission-Reflection and Digital Image Correlation

aboratory-scale experiments on intact rocks are critical to the development of physics-based fundamental understanding of various geophysical phenomena. In this work, the capability of ultrasonic wave transmission (T-mode) and reflection (R-mode) to monitor damage progression in uniaxially loaded prismatic intact rock specimens has been analyzed, as it is imperative to study the observations and document the capabilities of these techniques in a controlled environment. This study is novel in the sense that the R-mode linear ultrasonic testing (LUT) has been rarely employed in studying intact rock damage processes in a laboratory setting, with most of the studies utilizing a direct transmission (T-mode) approach or focusing on macroscopically fractured material. The two-dimensional digital image correlation (2-D DIC) full-field strain measurement approach was also used in-sync with the LUT monitoring to explicitly correlate the stress-induced damage in the specimens with the changes observed in the ultrasonic (T-mode and R-mode) signals. The results show that both the T-mode and R-mode LUT approaches are sensitive to detect the evolution of tensile and shear damage in the specimens, with the R-mode ultrasonic signals showing higher degree of sensitivity to the damage in the rocks, immediately following the initiation of damage in the rock volume.

58 GEOSCIENCES↗
DOE OSTI2023

Hydrogen and Silicon Effects on Hexagonal Close Packed Fe Alloys at High Pressures: Implications for the Composition of Earth's Inner Core

Hexagonal close-packed (hcp) structured Fe-Ni alloy is believed to be the dominant phase in the Earth's inner core. This phase is expected to contain 4%–5% light elements, such as Si and H. While the effects of individual light element candidates on the equation of state (EoS) of the hcp Fe metal have been studied, their combined effects remain largely unexplored. Here, in this study, we report the equations of state for two hcp-structured Fe-Si-H alloys, namely Fe 0.83 Si 0.17 H 0.07 and Fe 0.83 Si 0.17 H 0.46 , using synchrotron X-ray diffraction measurements up to 125 GPa at 300 K. These alloys were synthesized by cold compression of Fe-9wt%Si in either pure H 2 or Ar-H 2 mixture medium in diamond-anvil cells. The volume increase caused by a H atom in hcp Fe-Si-H alloys is approximately eight times greater than that by a Si atom. We used the improved data set to develop a composition-dependent EoS that covers a wide range of compositions. Our calculated density and bulk sound velocity of hcp Fe-Si-H alloys suggest a large trade-off between Si and H contents in fitting the seismic properties of the inner core. Combining our new EoS with geophysical and geochemical constraints, we propose 1.6–3 wt% Si and 0.15–0.6 wt% H in the Earth's inner core.

58 GEOSCIENCES↗
DOE OSTI2024

Thermal Equations of State of Magnesite: Implication for the Complex Mid‐Lower Mantle Seismic Scatterers

Abstract Magnesite (MgCO 3 ) entering the lower mantle together with the subducted oceanic crust is an important carbon carrier. The reaction between magnesite and mantle minerals has been documented, but its influence on the density and velocity profiles of lower mantle remains unexplored. To decipher the deep carbon transportation and its associated effect, here we determined the thermal equations of state of magnesite up to 120 GPa and 2600 K using X‐ray diffraction in laser‐heated diamond anvil cells. The obtained thermal elastic parameters of magnesite facilitated a comprehensive understanding on the influence of magnesite‐SiO 2 reaction, variation of carbon and SiO 2 content, and temperature on the origin of lower‐mantle scatterers at 1,000–1,800 km depth. Our modeling revealed that the depth of the lower‐mantleV S scatterers is mainly controlled by the Al 2 O 3 content in SiO 2 , while its magnitude depends on the SiO 2 content. Along normal geotherm, the magnesite‐SiO 2 reaction would occur before the post‐stishovite transition, consuming substantial SiO 2 in the subducted oceanic crust. Depending on the amount of residual SiO 2 , the post‐stishovite transition can produce a 2.5–5.2 (2)%V S reduction, compatible with the observed seismic scatterers in Izu‐Bonin and Mariana subduction zones. Along slab geotherm, this reaction occurs after the post‐stishovite transition, generating a greaterV S reduction of 4.4–6.4 (4)%. We thus propose that the reaction between sinking MgCO 3 and SiO 2 in the slab is one of the potential factors influencing the magnitude of the lower‐Vs scatterers at 1,000–1,900 km depth. Our results provide new insights into the deep‐mantle carbonate transportation influencing regional geophysics.

Geochemistry & Geophysics↗
DOE OSTI2025

Rising concerns of climate extremes and land subsidence impacts

A recent article in Reviews of Geophysics explores land subsidence drivers, rates, and impacts across the globe. It also discusses the need for improved process representations and the inclusion of the interplay among land subsidence and climatic extremes, including their effects in models and risk assessments. Here, we asked the lead author to explain the concept of land subsidence, its impacts, and future directions needed for improved mitigation.

Earth science↗
DOE OSTI2025

Exploring the Use of Non‐Invasive Drone‐Based Ground‐Penetrating Radar (GPR) to Characterize Biogenic Gas Dynamics in Subtropical Peat Soils

Peat soils are a critical component of the global carbon cycle as natural producers of biogenic greenhouse gases (e.g., methane and carbon dioxide) that accumulate within the soil and are released to the atmosphere. Previous studies have showed the ability of ground-based minimally-invasive geophysical methods such as ground-penetrating radar (GPR) to characterize carbon dynamics in peat soils. However, ground-based GPR is limited by scale of measurement and soil disturbance potentially altering gas releases during deployment. Here, we explore the potential of drone-based GPR for identification of hot spots and hot moments of gas accumulation and release in subtropical soils. Here, we collected drone-based GPR data sets across two grids (∼17,500 m 2 ) in the Everglades during January (dry season), September, and November (wet season) of 2023 to characterize peat thickness and seasonal variability of gas content. Results show that drone-based GPR is effective and efficient for: (a) capturing the temporal variation of in situ biogenic gas content in peat soils with changes between 1% and 25 % volumetric gas content over repeatable grids; (b) inferring a total peat thickness between 0.8 and 1.2 m; and (c) estimating flux releases of 63 and 135 mg CH 4 m −2 day −1 for specific locations and periods that are strikingly consistent with our coincident gas trap measurements. This work also indicates that (a) spatial distribution of gas content in the Everglades is strongly controlled by landscape morphology such as ridges and sloughs and (b) the temporal variation of gas content is seasonal with increased gas production during the wet season.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Trans-crustal structural control of CO 2 -rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

58 GEOSCIENCES↗
DOE OSTI2022

Characterization of a 1D-imaging high-energy x-ray backlighter driven by the National Ignition Facility Advanced Radiographic Capability laser

Plastic deformation of samples compressed to Mbar pressures at high strain rates at the National Ignition Facility (NIF) forms the basis of ongoing material strength experiments in conditions relevant to meteor impacts, geophysics, armor development, and inertial confinement fusion. Hard x-ray radiography is the primary means of measuring the evolution of these samples, typically employing a slit-collimated high-Z microdot driven by the NIF laser to generate >40 keV x rays [E. Gumbrell et al., Rev. Sci. Instrum. 89, 10G118 (2018) and C. M. Huntington et al., Rev. Sci. Instrum. 89, 10G121 (2018)]. Alternatively, a dysprosium “micro-flag” target driven by the Advanced Radiographic Capability laser (~2 kJ, 10 ps) can deliver significantly higher spatiotemporal resolution [M. P. Hill et al., Rev. Sci. Instrum. 92, 033535 (2021)], especially in high-opacity samples. Here, initial experiments revealed problematic brightness and spectral gradients from this source, but by radiographing a set of diamond-turned, 105 µm-thick Pb test objects and supported by simulations using the 3D Monte Carlo code GEANT4, these geometry-dependent gradients across the field of view are quantified and mitigation strategies are assessed. In addition to significantly enhancing the modulation transfer function compared to the existing system, image stacking from multiple layers of image plate is shown to almost double the signal to noise ratio that will reduce uncertainties in future dynamic strength experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A generative artificial intelligence framework for long-time plasma turbulence simulations

Generative deep learning techniques are employed in a novel framework for the construction of surrogate models capturing the spatiotemporal dynamics of 2D plasma turbulence. The proposed Generative Artificial Intelligence Turbulence (GAIT) framework enables the acceleration of turbulence simulations for long-time transport studies. GAIT leverages a convolutional variational auto-encoder and a recurrent neural network to generate new turbulence data from existing simulations, extending the time horizon of transport studies with minimal computational cost. The application of the GAIT framework to plasma turbulence using the Hasegawa–Wakatani (HW) model is presented, evaluating its performance via various analyses. Very good agreement is found between the GAIT and the HW models in the spatiotemporal Fourier and Proper Orthogonal Decomposition spectra, the flow topology characterized by the Okubo–Weiss parameter, and the time autocorrelation function of turbulent fluctuations. Excellent agreement has also been obtained in the probability distribution function of particle displacements and the effective turbulent diffusivity. In-depth analyses of the latent space of turbulent states, choice of hyperparameters and alternative deep learning models for the time prediction are presented. Our results highlight the potential of Artificial Intelligence-based surrogate models to overcome the computational challenges in turbulence simulation, which can be extended to other situations such as geophysical fluid dynamics.

Artificial intelligence↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
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 OSTI2024

Ensemble Kalman inversion of induced polarization data

SUMMARY This paper explores the applicability of ensemble Kalman inversion (EKI) with level-set parametrization for solving geophysical inverse problems. In particular, we focus on its extension to induced polarization (IP) data with uncertainty quantification. IP data may provide rich information on characteristics of geological materials due to its sensitivity to characteristics of the pore–grain interface. In many IP studies, different geological units are juxtaposed and the goal is to delineate these units and obtain estimates of unit properties with uncertainty bounds. Conventional inversion of IP data does not resolve well sharp interfaces and tends to reduce and smooth resistivity variations, while not readily providing uncertainty estimates. Recently, it has been shown for DC resistivity that EKI is an efficient solver for inverse problems which provides uncertainty quantification, and its combination with level set parametrization can delineate arbitrary interfaces well. In this contribution, we demonstrate the extension of EKI to IP data using a sequential approach, where the mean field obtained from DC resistivity inversion is used as input for a separate phase angle inversion. We illustrate our workflow using a series of synthetic and field examples. Variations with uncertainty bounds in both DC resistivity and phase angles are recovered by EKI, which provides useful information for hydrogeological site characterization. Although phase angles are less well-resolved than DC resistivity, partly due to their smaller range and higher percentage data errors, it complements DC resistivity for site characterization. Overall, EKI with level set parametrization provides a practical approach forward for efficient hydrogeophysical imaging under uncertainty.

Geochemistry & Geophysics↗
DOE OSTI2021

An Automatic Processing Framework for In Situ Determination of Ecohydrological Root Water Content by Ground-Penetrating Radar

This report concerns Root water content (RWC) as a vital component in water flux in soil-plant-atmosphere continuum. Knowledge of RWC helps to better understand the root function and the soil-root interaction and improves water cycle modeling. However, due to the lack of appropriate methods, field monitoring of RWC is seriously constrained. In this study, we used ground-penetrating radar (GPR), a common geophysical technique, to characterize RWC of coarse roots noninvasively. An automatic GPR data processing framework was proposed to (1) identify hyperbolic root reflections and locate roots in GPR images and (2) extract waveform parameters from the reflected wave of identified roots. These waveform parameters were then used to establish an empirical model and a semiempirical model to determine RWC. We validated the developed models using GPR root data at three antenna center frequencies (500 MHz, 900 MHz, and 2 GHz) that were produced from simulation experiments (with RWC ranging from 70% to 150%) and field experiments in sandy soils (with RWC ranging from 66% to 144%). Our results show that both the empirical and the semiempirical models achieved a good performance in estimating RWC with similar accuracy, i.e., the prediction error [root-mean-square error (RMSE)] was less than 8% for the simulation data and 12% for the field data. For both models, the accuracy of RWC estimation was the highest when applied to 2-GHz data. This study renders a new opportunity to determine RWC under field conditions that enhances the application of GPR for root study and the understanding and modeling of ecohydrology in the rhizosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers part I: Computational models and workflow

Computational simulations have become central to the seismic analysis and design of major infrastructure over the past several decades. Most major structures are now “proof tested” virtually through representative simulations of earthquake-induced response. More recently, with the advancement of high-performance computing (HPC) platforms and the associated massively parallel computational ecosystems, simulation is beginning to play a role in increased understanding and prediction of ground motions for earthquake hazard assessments. However, the computational requirements for regional-scale geophysics-based ground motion simulations are extreme, which has restricted the frequency resolution of direct simulations and limited the ability to perform the large number of simulations required to numerically explore the problem parametric space. In this article, recent developments toward an integrated, multidisciplinary earth science-engineering computational framework for the regional-scale simulation of both ground motions and resulting structural response are described with a particular emphasis on advancing simulations to frequencies relevant to engineered systems. This multidisciplinary computational development is being carried out as part of the US Department of Energy (DOE) Exascale Computing Project with the goal of achieving a computational framework poised to exploit emerging DOE exaflop computer platforms scheduled for the 2022–2023 timeframe.

58 GEOSCIENCES↗