Engineering PapersSearch

SEARCH · Engineering Papers

Results for “seismic waves”

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

Analysis of Seismic Wave Amplification in Sedimentary Basins Using 3D Wavefield Simulations: Nenana Basin, Central Alaska

We investigate seismic wave amplification in Nenana basin, central Alaska, using 3D seismic wavefield simulations. We quantify frequency-dependent amplification by comparing synthetic seismograms among four different seismic velocity models: (a) a tomographic model, (b) a tomographic model with the uppermost 6.5 km replaced by a 1D basin profile, (c) a tomographic model with an embedded ellipsoidal basin, and (d) a tomographic model with an embedded realistic basin. For each model we perform wavefield simulations accurate up to 2 Hz for a set of 10 different earthquakes, which provide a range of frequencies and incident angles for waves interacting with the basin. Average amplification ratios are four on the horizontal components and seven on the vertical component. The dominant influence of the amplification is the 3D geometry of the basin, not the slow-velocity profile near the surface. Our synthetic amplification ratios are in general agreement with amplification estimates obtained from 14 stations that recorded the same set of earthquakes. Our approach offers a general strategy for documenting frequency-dependent basin amplification for a region with realistic basin structures and local earthquakes.

Geosciences

The propagation of seismic waves, misinformation, and disinformation from the 2024-10-05 M 4.5 Iran earthquake

The 2024-10-05 Iran M 4.5 earthquake took place at a time of heightened tensions in the Middle East. We perform a discrimination and moment tensor analysis and identify a shallow-dipping, reverse fault source commensurate with the compressional setting of the Iranian interior. Nonetheless, the event's aftermath saw widespread dissemination of misinformation, and potentially active disinformation, concluding that it was in fact a test of an Iranian nuclear weapon. The 'evidence' for many of these claims was based on inaccurate interpretation of seismic data. In this paper, we analyze how geophysical 'fake news' propagated through social media (mainly Twitter/X) following this event, eventually gaining traction in mainstream, earned media. This event is an illustrative warning of how seismic data can be misinterpreted and/or manipulated in public discourse.

58 GEOSCIENCES

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (

owl

O4921: Open Wave Library for seismic wave modeling and FWI

Gao, Kai [Los Alamos National Laboratory]

Shear-wave detection

Angular accelerometers and transducers compared in detection of seismic waves and shear waves

TRANSDUCER

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE

Limited Dynamic Earthquake Triggering in Nevada

Dynamic triggering occurs when seismic waves from distant large earthquakes temporarily alter stress conditions along faults, potentially triggering new earthquakes hundreds to thousands of kilometers away from the source. Previous studies have linked triggered seismicity to anthropogenic activities such as geothermal, oil, and gas production. Although these activities are present in Nevada, little work has been conducted to explore dynamically triggered seismicity in Nevada. Here, we analyze a newly published, high-resolution earthquake catalog for Nevada to identify local seismicity dynamically triggered by teleseismic events (Mw≥7) from 2008 to 2023. We identify 94 dynamically triggered earthquakes concentrated in four distinct regions, which qualitatively show a modest positive correlation with geothermal well locations. Triggered seismicity in Nevada is predominantly delayed, with some instantaneously triggered by Rayleigh waves. The prevalence of delayed triggering indicates that pore fluid interactions may play a critical role in controlling dynamic triggering susceptibility in Nevada. Our results demonstrate that dynamic triggering can provide valuable insight to help identify critically stressed regions.

58 GEOSCIENCES

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala

FY 2026 Midyear Report: Seismic Monitoring of Underground Vibration Sources Using Distributed Acoustic Sensing and Seismometers

Safeguards-relevant temporal changes in underground facilities can be observed using geophysical monitoring techniques. Seismic waves, in particular, provide valuable insights into subsurface activities and can serve as an important tool for detecting anomalous events that may indicate containment breaches at geological repositories. This midyear report summarizes ongoing efforts to automatically and rapidly detect and locate anomalous vibration signals that could be indicative of potential containment breaches. Previous work during FY25 focused on compiling continuous seismic datasets from two underground sites and developing a database of continuous waveforms and ground-truth event data derived from multiple sensing modalities. Building on this foundation, we are adapting anomaly detection and geolocation algorithms to explore methods for monitoring underground activities using two relatively low-maintenance sensing technologies: a dense surface geophone array deployed at the Pleasant Gap mine in Pennsylvania, and a three-dimensional fiber-optic cable array for distributed acoustic sensing (DAS) installed in the subsurface at the Sanford Underground Research Facility (SURF) in South Dakota. This report summarizes work conducted during the first two quarters of FY26, during which we refined a dynamic power spectral density (PSD)-based detector, applied it independently to each geophone station, and then combined the per‑station detections with density-based spatial clustering of applications with noise (DBSCAN) to cluster events and produce spatial maps over a nine‑day interval. In addition, we outline plans for a field trial at the Waste Isolation Pilot Plant (WIPP) in New Mexico to compare traditional seismic monitoring approaches with DAS techniques and to evaluate the benefits of combined data analysis. Activities during the past two quarters have included the preparation and submission of a Field Test Plan to WIPP for approval, as well as submission to headquarters for review and feedback.

58 GEOSCIENCES

Seismic monitoring of underground vibration: database of seismic data and ground truth

Seismic waves provide valuable insights into underground activities, serving as an essential tool for monitoring anomalies that could signal containment breaches in geological repositories. We aim to test and refine underground detection and geolocation techniques to identify anomalous vibration signals indicative of potential breaches, thereby strengthening georepository safeguards. We evaluate the effectiveness of two distinct, low-maintenance sensing technologies (surface geophones and underground distributed acoustic sensing (DAS) fiber optic cable) leveraging existing datasets. This report details the experimental designs, instrumentation, and data characteristics for both seismic and DAS arrays. Additionally, we provide a ground truth database documenting relevant operational activities for each experiment.

58 GEOSCIENCES

Prediction of Complex Observed Shear Wave Splitting Patterns at Ryukyu Subduction Zone Using a Strong Intra-Slab Anisotropy Model

Complex shear wave splitting (SWS) patterns in subduction zones are often interpreted geodynamically as resulting from complex mantle flow; however, this may not always be necessary. We analyzed 7,093 high-quality SWS measurements from teleseismic S waves recorded by Hi-net stations across the Ryukyu arc in Japan. Our findings show a systematic rotation of the fast S polarization from trench-parallel to trench-perpendicular depending on the earthquake backazimuth. For the same earthquake, the measured splitting patterns also vary spatially across the southwest Japan. Using full-wave seismic modeling, we showed that a dipping slab with ~30% shear anisotropy of the tilted transverse isotropy (TTI) type, with a symmetry axis perpendicular to the slab interface, can predict the observed delay times and polarization rotation. Our results highlight the importance of considering dipping anisotropic slabs in interpreting SWS at subduction zones.

58 GEOSCIENCES

Symmetrization of Strong Hydrogen Bond under High Pressure in Bihydroxide-Ion-Containing NaCu 2 (SO 4 ) 2 ·H 3 O 2 Revealed by Experimental Charge Density, Single-Crystal Electron Diffraction, and Neutron Diffraction Studies

In minerals and inorganic compounds, strong hydrogen bonding can lead to the formation of complex ionic species such as the H 3 O 2 – bihydroxide anion and Zundel cation H 5 O 2 + . We studied [NaCu 2 (SO 4 ) 2 ·H 3 O 2 ] natrochalcite, which contains bihydroxide anions and undergoes hydrogen bond symmetrization at the lowest pressure reported so far among inorganic compounds. Hydrogen bond symmetrization leads to changes in the bulk modulus, seismic wave velocities, and proton mobility and plays a primary role in high-temperature superconductivity, but its characteristics are not well understood due to a lack of systematic studies and limitations of experimental methods sensitive to this subtle change. In this work, we applied experimental charge density analysis based on in situ single-crystal X-ray diffraction data, along with the single-crystal neutron and electron diffraction experiments, to probe the behavior of hydrogen atoms during the hydrogen bond symmetrization process under high-pressure conditions. On the way to the symmetrical H-bonding, natrochalcite undergoes a series of complex redistributions of electron density, which we trace with multipole refinement and detailed analysis of changes in the Laplacian of electron density values. Additionally, we deconvoluted the equation of state (volume of the unit cell vs pressure relation) into the atomic equation of states describing dependencies of atomic charges or volumes vs pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Insights into Slip-Rate Time Functions, Rupture Parameter Correlations, and Ground Motions from Validated Multicycle Earthquake Ruptures

Earthquake strong-motion predictions using kinematic source modeling require knowledge of the slip-rate functions (SRFs) along the rupture and their distinct characteristics in asperities, background (off-asperity) areas and near the surface. Here, in this study we analyzed SRFs from well-validated, self-consistent, and fully dynamic rupture models from earthquake cycles obeying a rate-and-state friction law, from our companion study (Galvez et al., 2021). The shapes of SRFs in asperities are well described by the regularized Yoffe function (RYF), which has only two parameters: rise time T r and smoothing time T s , which control the generation of long- and short-period ground motions, respectively. In background areas, we demonstrate that, in addition to the primary rupture, multiple secondary ruptures may also nucleate from rupture heterogeneities related to asperities, resulting in SRFs with multiple peaks. Because it is impossible to fit a multiple-peak SRF by the single-peak RYF, we describe SRFs in background areas in an effective way by fitting their amplitude spectra with the RYF spectra. Such spectrally effective RYFs capture salient aspects of seismic-wave generation and can be used in rupture generators for strong motion prediction. We found that small T s values correlate with small characteristic weakening distances, large peak slip rates (PSRs), and large rupture velocities. T r values are larger in background areas and smaller in asperities. Within the shallow aseismic zone, T s values approximately quadruple whereas T r values approximately double. Because of this dominant T s increase, PSR values decrease in the near-surface zone. These features indicate that the generation of strong motions by the near-surface portions of the rupture is negligible in the studied scenarios.

Geosciences

S&TR September 2025: Computing Grand Challenge Turns 20

Livermore’s Computing Grand Challenge Program enters its 20th year with more unclassified high-performance computing (HPC) power than ever before. This unique, peer-reviewed competition awards HPC allocations on top supercomputers to multidisciplinary teams with high-impact projects. The Grand Challenge encourages researchers to innovate, pushes scientific discovery to new heights, improves the Laboratory’s HPC capabilities, and extends HPC accessibility to collaborators. Awardees must adapt to successive generations of HPC hardware and learn to run simulations at scale. The feature article spotlights three Grand Challenge teams whose research broke new ground in key scientific pursuits—the essence of dark matter, explosion-generated seismic waves, and protein interactions linked to cancer—while underscoring the importance of academic partnerships and considering the program’s future.

07 ISOTOPE AND RADIATION SOURCES

Iron Bonding with Light Elements: Implications for Planetary Cores Beyond the Binary System

Light element alloying in iron is required to explain density deficit and seismic wave velocities in Earth’s core. However, the light element composition of the Earth’s core seems hard to constrain as nearly all light element alloying would reduce the density and sound velocity (elastic moduli). The alloying light elements include oxidizing elements like oxygen and sulfur and reducing elements like hydrogen and carbon, yet their chemical effects in the alloy system are less discussed. Moreover, Fe-X-ray Absorption Near Edge Structure (Fe-XANES) fingerprints have been studied for silicate materials with ferrous and ferric ions, while not many X-ray absorption spectroscopy (XAS) studies have focused on iron alloys, especially at high pressures. To investigate the bonding nature of iron alloys in planetary interiors, we presented X-ray absorption spectroscopy of iron–nitrogen and iron–carbon alloys at high pressures up to 50 GPa. Together with existing literature on iron–carbon, –hydrogen alloys, we analyzed their edge positions and found no significant difference in the degree of oxidation among these alloys. Pressure effects on edge positions were also found negligible. Our theoretical simulation of the valence state of iron, alloyed with S, C, O, N, and P also showed nearly unchanged behavior under pressures up to 300 GPa. This finding indicates that the high pressure bonding of iron alloyed with light elements closely resembles bonding at the ambient conditions. We suggest that the chemical properties of light elements constrain which ones can coexist within iron alloys.

Advanced Photon Source