Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMIC WAVE”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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

Determining Stress Orientation in Rock Valley, Nevada, Using Ambient Seismic Noise

The stress field and the mechanical properties of rocks are important to consider for nuclear explosion monitoring due to their effect on seismic wave radiation from earthquakes and explosions. At the Rock Valley Direct Comparison site, the regional orientation of the maximum horizontal compressive stress (SH max ) is well constrained, but it is unknown whether there are local heterogeneities. Here, I show that stress-induced anisotropy in nonlinear elasticity can be used to estimate the orientation of SH max . Rocks have compliant internal contacts, such as fractures and mineral grain boundaries, that respond to applied strains more strongly than individual mineral crystals. This strain response is asymmetric between compressive and dilatational strains and is affected by anisotropy in the ambient stress field. Traditional seismic velocity measurements are of strain-averaged velocity, which is less sensitive than nonlinear elasticity to the behavior of compliant internal contacts and therefore the stress field and fracture behavior. My results show that the orientation of SH max measured using ambient seismic noise aligns with regional estimates, with some heterogeneity that can also be explained by the limitations of the seismic array. I demonstrate the potential of a passive technique for monitoring the stress field in places that currently lack local measurements.

58 GEOSCIENCES↗
DOE OSTI2025

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

Migration and accumulation of hydrous mantle incipient melt in the Earth's asthenosphere: Constraints from in-situ falling sphere viscometry measurements

The Earth's oceanic lithosphere-asthenosphere boundary (LAB) is marked by a notable decrease in seismic wave velocities and an increase in attenuation. This phenomenon is likely attributed to the accumulation of partial melt at the top of the asthenosphere. Nevertheless, the process involving the upward migration and aggregation of low-degree partial melts, highly likely to be mantle incipient melts in the asthenosphere, remains underexplored. Viscosity is a key factor controlling the flow of melts, thus in this study, we used In-situ X-ray falling sphere viscometry experiments to determine the viscosity of the mantle incipient melt containing 1.8–4.0 wt.% water at 1.5–6 GPa and 2100 K. We found that water and pressure can effectively decrease the viscosity of the mantle incipient melt. Therefore, the mantle incipient melts are highly mobile with a high segregation velocity in the deeper part of the asthenosphere. However, during ascent, their mobility diminishes due to an increase in viscosity and a decrease in water content. In conclusion, these mechanisms potentially contribute to the accumulation of partial melt at the top of the asthenosphere.

58 GEOSCIENCES↗
DOE OSTI2025

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

Shallow Soil Response to a Buried Chemical Explosion With Geophones and Distributed Acoustic Sensing

Abstract Shallow sediments can respond non‐linearly to large dynamic strains and undergo a subsequent healing phase as the material gradually recovers following the passing of seismic waves. This study focuses on the physical changes in the subsurface caused by the shaking from a buried chemical explosion detonated in a borehole in Nevada, USA, as a part of the Source Physics Experiment Phase II. The explosion damaged the shallow subsurface and modified the frequency content recorded by 491 geophones and 2240 Distributed Acoustic Sensing (DAS) channels within 2.5 km from surface ground zero. We observe a gradual shift of resonance frequencies in the 10–25 Hz frequency band in the hours following the explosion and develop a method to characterize the related logarithm‐type healing process of the shallow (i.e., upper ∼25 m) subsurface. We find that stronger levels of ground motion increase the relative degree of damage and duration of the subsurface healing; with the spall region exhibiting the largest degree of damage and longest healing recovery time. We observe coherent spatial patterns of damage with the region located to the southeast of the explosion exhibiting more damage than the southwest region. This study demonstrates that both DAS and co‐located geophones capture similar temporal changes associated with the physical processes occurring in the subsurface, with the high‐density sampling of DAS measurements enabling a new capability to monitor the fine‐scale changes of the Earth's shallow subsurface following the detonation of a buried explosion.

Viens, Loïc↗
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

High-Pressure Melting Experiments of Fe 3 C and a Thermodynamic Model of Fe-C Liquids for the Earth's Core

Melting experiments of Fe 3 C were conducted to 85 GPa in laser-heated diamond anvil cells with in situ X-ray diffraction and post-experiment textural observation. From the determined pressure-temperature conditions of the melting curve for Fe 3 C, together with literature data on the melting point of diamond and eutectic point of the system Fe-Fe 3 C/Fe 7 C 3 under high pressures, we established a self-consistent thermodynamic model for high-pressure melting of the system Fe-C including the mixing parameters for liquids. The results show that mixing of Fe and C liquids is negatively nonideal from 1 bar to the pressure at the center of the Earth. The departure from ideal mixing becomes progressively larger with increasing pressure, which leads to greatly stabilized liquids under core pressures. The modeled carbon content in eutectic melts under core pressures is 3.3–4.4 wt%. From the Gibbs free energy, we derived an internally consistent parameters for Fe-C outer cores which included the crystallizing points at their bottoms, isentropic thermal profiles, and densities and longitudinal seismic wave speeds (Vp). While the addition of carbon in excess of the eutectic melt composition effectively reduces the density of iron liquid, the Vp of iron liquid is not greatly changed. Therefore, the low density and high Vp of PREM relative to pure iron cannot be reconciled by an Fe-C liquid. Therefore, the Earth's core cannot be approximated by the system Fe-C and should include another light element.

58 GEOSCIENCES↗
DOE OSTI2024

Utah FORGE: Direct Shear Test Data for Investigating Seismic Precursors to Shear Failure of Fractures

This dataset includes results of direct shear tests to investigate the mechanical and geophysical response of dry and saturated fractures in Indiana limestone and Sierra White granite. Direct shear tests were performed on tensile-induced fractures in Indiana limestone and Sierra White granite in a custom water-pressurized chamber. The provided Excel files include the representative seismic wave signals and the normalized wave amplitudes of ultrasonic wave transducers. A link to the published journal article presenting the data and describing the experiment in detail is provided as well.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Proximal Observations of Epicentral Infrasound Generated by Shallow Low-Magnitude Earthquakes in the Permian Basin, West Texas

Infrasound generated by earthquakes and explosions is generally detected at receivers at epicentral distances of hundreds to thousands of kilometers. However, proximal (<50 km) observations are especially important for low-magnitude earthquakes and low-yield explosions that may not generate signals capable of being detected at great ranges. Here, in this study, we present on the signals detected on an infrasound array 3 km away from two M L 2.9 earthquakes in the Permian Basin of west Texas. Local infrasound (LIS), generated at receivers during the passage of seismic waves, was detected following each earthquake. Epicentral infrasound (EIS), created at or near the epicenter and propagating away as a sound wave, was also detected. Array processing methods show that the EIS signals arrive from the same direction as the earthquake epicenters and at acoustic speeds. To our knowledge, these are the first observations of laterally propagating EIS at proximal ranges following an earthquake of any magnitude.

58 GEOSCIENCES↗
DOE OSTI2025

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

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

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

The Influence of Shallow Subsurface Properties on Particle Motion in Acoustic-Seismic Coupling

Atmospheric acoustic waves transmit energy into the solid Earth through air-to-ground coupling. These waves are recorded by seismic sensors and provide insight into both atmospheric phenomena and subsurface properties. Interpreting these signals is often challenging because they are modulated by subsurface structure and the incidence angle of the acoustic wave. This study examines acoustic--seismic coupling generated by the 2012 Camp Minden Explosion, which was recorded by hundreds of seismoacoustic stations. We apply a novel technique to quantify the seismic particle motion, model coupled waves with a propagator matrix approach, and apply a Bayesian inversion to infer properties of the shallow subsurface. Our analysis reveals that prograde motion is widespread and focused in low shear-wave velocity regions, such as the Mississippi Embayment, and retrograde motion is more common in higher shear-velocity areas. Inversion results at some stations produce plausible subsurface models with strong waveform fits, while inversions at other sites are less successful. These results indicate prograde particle motion in air-to-ground coupled waves is more prevalent than previously recognized and may serve as a diagnostic for shallow velocity structure. Our comprehensive modeling and inversion framework provides a potential method to extract layered near-surface properties from acoustic-seismic coupling observations.

58 GEOSCIENCES↗
DOE OSTI2026

Elastic-wave sensitivity-guided adaptive seismic survey design for cost-effective monitoring of geological carbon storage

Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments

58 GEOSCIENCES↗
DOE OSTI2025

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗
DOE OSTI2025

​​Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities​

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗