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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analysis of Fault Zone Resonance Modes Recorded by a Dense Seismic Array Across the San Jacinto Fault Zone at Blackburn Saddle

We present observations and modeling of spatial eigen-functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5-Hz geophones that cross the SJFZ with ~10-30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after the S wave arrival are consistent for more than 50 near-fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio shows three main frequency peaks at ~1.3 Hz, ~2.0 Hz, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen-functions of resonance waves in a low-velocity layer (fault zone) between two quarter-spaces. Using a grid search-based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65% S wave velocity reduction compared to the host rock. The SW edge of the low-velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together.

58 GEOSCIENCES↗
DOE OSTI2023

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗
DOE OSTI2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗
DOE OSTI2024

Seismic monitoring of EGS fracture stimulations at Utah FORGE (Part 1): Time-lapse variations of b-values and Shear-wave Splitting Rates of induced microearthquakes

Characterizing geologic formations near geothermal wells is crucial for efficient extraction of subsurface geothermal energy resources. To monitor effective fracture stimulation of enhanced geothermal systems (EGS), we analyze a catalog of induced microearthquakes to characterize seismicity induced by Phase 2C fracture stimulations of the Utah FORGE project and derive the b-values for 424 induced microearthquakes in three stages of fracture stimulation. Here, we find the time-lapse variations of these b-values. We apply eigenvector rotation and waveform correlation to 3C geophone data to separate fast S 1 and slow S 2 waves, determine their directions, measure their lag time, and calculate the S-wave splitting rate (SSR) for six induced microearthquakes in Stage 1 and 21 in Stage 3. Both SSR values and b-values increase with time. The increase in these values may be related to heavily fractured zones created by fracture stimulations during Stage 2 injections. Large and highly variable SSR values in Stage 3 may reflect heterogeneities in the fracture zones. Time-lapse variations of SSR values and b-values may be one of the most useful approaches to monitoring and characterizing the effectiveness and efficiency of various types of fracture stimulations.

58 GEOSCIENCES↗
DOE OSTI2022

Estimating Explosion Yields Using Moment Tensor Solutions and Seismic Moment

We report seismic moment, a measurable and well-understood quantity of seismic sources, is used to estimate the yield of explosions. Application of such a method in the past, as in the manner of m b -derived yields, has been complicated by the effect of variations in the explosion working point, depth, and secondary source effects (such as spalling and tectonic release) on the observed moment. We start using the full (six-element) moment tensor solution, which can capture the relevant source physics and, at least in theory, better isolate the primary explosion source. The moment-to-yield ratio is then estimated using an explosion source model which, provided with emplacement conditions, can relate the two parameters. We discuss the major sources of uncertainty associated with the method, and calibrate it with chemical and nuclear explosions at the Nevada National Security Site. We then apply the method to published moment tensor solutions for the six declared North Korean nuclear explosions that occurred between 2006 and 2017. The results are mostly consistent with other yield estimates made using a variety of high-frequency methods. This technique is a new approach to estimating explosive yield and simple to implement, as much of the complexity is captured by the source models.

58 GEOSCIENCES↗
DOE OSTI2020

Modulation of seismic noise near the San Jacinto fault in southern California: origin and observations of the cyclical time dependence and associated crustal properties

SUMMARY We examine the cyclic amplitude variation of seismic noise recorded by continuous three-component broad-band seismic data with durations spanning 91–713 d (2008–2011) from three different networks: Anza seismic network, IDA network and the transportable seismic array. These stations surround the San Jacinto Fault Zone (SJFZ) in southern California. We find the seismic noise amplitudes exhibit a cyclical variation between 0.3 and 7.2 Hz. The high-frequency (≥0.9 Hz) noise variations can be linked to human activity and are not a concern. Our primary interest is signals in the low frequencies (0.3–0.9 Hz), where the seismic noise is modulated by semi-diurnal tidal mode M2. These long-period (low-frequency) variations of seismic noise can be attributed to a temporal change of the ocean waves breaking at the shoreline, driven by ocean tidal loading. We focus on the M2 variation of seismic noise at f = 0.6 Hz, travelling distances of ∼92 km through the crust from offshore California to the inland Anza, California, region. Relative to the shoreline station, data from the inland stations show a phase lag of ∼ –12°, which we attribute to the cyclic change in M2 that can alter crustal seismic attenuation. We also find that for mode M2 at 0.6 Hz, the amplitude variations of the seismic quality factor (Q) depend on azimuth and varies from 0.22 per cent (southeast to northwest) to 1.28 per cent (northeast to southwest) with Q = 25 for Rayleigh waves. We propose the direction dependence of the Q variation at 0.6 Hz reflects the preferred orientation of subfaults parallel to the main faulting defined by the primarily N45°W strike of the SJFZ.

Geochemistry & Geophysics↗
DOE OSTI2021

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗
DOE OSTI2021

Seismic Spatial Gradients and Machine Learning-Based Classifiers for Explosion Monitoring (LDRD 218327)

This final report summarizes the work completed under the Laboratory Directed Research and Development (LDRD) project “Seismic Spatial Gradients as a Machine Learning-Based Classifier for Explosion Monitoring.” The overarching goal of the project was to explore the efficacy of using machine learning-based classification algorithms where the input data are the spatial gradient of the seismic wavefield collected at a single point on the Earth’s surface. The methods that I describe here are in direct contrast to conventional methods of seismic discrimination which typically rely on a spatially extended network of instruments and physics-based wavefield attributes such as, for example, the ratio between $\textit{P}$ and $\textit{S}$ waves. Rather, we use the spatial gradient of the seismic wavefield observed at a single point on the Earth’s surface and data processing approaches inspired by the machine learning community. We tested two algorithms, a neural network and a modified version of principal component analysis termed Spectrally Filtered Principal Component Analysis (SFPCA). To test these algorithms, we first conducted a series of numerical tests using synthetic data and then conducted a small-scale controlled field experiment. The tests using synthetic data showed that both algorithms had high success rates on gradiometric data, even when simulated noise was added to the signal. Furthermore, we found that using seismic spatial gradients increased the performance of our discrimination algorithms when compared to using just the traditional translational motion seismic data. The tests with field data also showed a high degree of discriminative success.

58 GEOSCIENCES↗
DOE OSTI2021

Illuminating the North Korean nuclear explosion test in 2017 using remote infrasound observations

North Korea conducted its sixth underground nuclear explosion test (mb 6.3) on 2017 September 3. The underground explosion produced substantial low-frequency atmospheric waves, which were detected by infrasound arrays located up to a distance of 566 km. These infrasound waves are formed by the conversion of seismic energy to acoustic energy across the lithosphere–atmosphere interface. While infrasound records at regional distances produce estimates of ground motion amplitude over spatially extended regions covering about 26 500 km 2 , 3-D full seismo-acoustic simulations within the lithosphere and atmosphere provide quantitative information about seismo-acoustic energy partitioning. Our results demonstrate the capability of remote infrasound observations combined with 3-D propagation modelling to further develop discrimination methods for underground sources. These results contribute to enhance the confidence of source identification and characterization in nuclear test monitoring research, which is essential for the enforcement of the Comprehensive Nuclear-Test-Ban Treaty.

58 GEOSCIENCES↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2025

Utah FORGE: Seismic DAS and Geophone Borehole Data Processing and 3D Imaging of Vp/Vs Ratio in the 2024 Stimulated Reservoir

This dataset includes a final report and a 3D velocity model derived from seismic DAS and geophone borehole data collected during the April 2024 stimulation of the reservoir at Utah FORGE. The report details the processing of over 50,000 P- and S-wave travel times used in a tomographic inversion to estimate the Vp/Vs ratio, revealing anomalies adjacent to wells 16A and 16B that may be associated with injected fracturing fluids. Additionally, Wadati analysis of more than 27,000 pairs of differential P- and S-wave travel times supports this interpretation with a high Vp/Vs estimate of 1.86, averaged over the dimensions of the earthquake cluster. The accompanying data file provides a 3D model of P- and S-wave velocities and Vp/Vs ratios, structured with spatial coordinates and velocity values.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Mineralogy, fabric and deformation domains in D" across the southwestern border of the African LLSVP

Recent advances in seismic anisotropy studies that jointly use reflections and shear wave splitting have proven to place tight constraints on the plausible anisotropic and deformation scenarios in the D" region. We apply this novel methodology to a large area of the D" region beneath the South Atlantic, in proximity to and within the African large low seismic velocity province (LLSVP). This area of the mantle is characterized by a transition from fast to slow seismic velocity anomalies and it is thought to be the location of deep-seated plumes responsible for hotspot volcanism. Attempting to probe mantle composition and deformation along the LLSVP borders may provide key information on mantle dynamics. By analysing seismic phases sampling this region, we detect a D" discontinuity over a large area beneath the South Atlantic, with inferred depth ranges ~170 to ~240 km above the core–mantle boundary. We find evidence for a D" reflector within the area of the LLSVP. Shear wave splitting observations suggest that anisotropy is present in this region of the mantle, in agreement with previous studies that partially sampled this region. We model the observations considering lattice- and shape-preferred orientation of materials expected in the D" region. A regional variation of mineralogy, phase transition boundaries, and deformation direction is required to explain the data. We infer two distinct domains of mineralogy and deformation: aligned post-perovskite outside the LLSVP and aligned bridgmanite within the LLSVP. Finally, the scenario depicted by this study agrees well with the current hypotheses for the composition of the LLSVP and with the prevalence of vertical deformation directions expected to occur along the LLSVPs borders.

58 GEOSCIENCES↗
DOE OSTI2022

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

Local-Regional Surface Wave Dispersion Analysis using Ambient Noise Data of the Southern Nevada National Security Site Area

This report details seismic ambient noise analysis to improve seismic velocity estimates of the southeastern Nevada National Security Site (NNSS). We compare two different methods for estimating surface wave dispersion curves from ambient noise cross-correlations: frequency time analysis (FTAN) and Aki’s cross-spectral method (XSpec). We find that XSpec performs better for our local-regional dataset and frequency content. Using phase velocity estimates from XSpec, we build a preliminary phase velocity dispersion dataset, which contains data for 1,054 station pairs and 8,905 discrete phase velocity measurements for periods between 0.5 and 13 s. This phase velocity dataset will be incorporated into an updated local-regional P and S wave velocity model of the southeastern NNSS in the future and is expected to improve upon shallow velocity estimates.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Survey of Deep Learning and Physics-Based Approaches in Computational Wave Imaging

Computational wave imaging (CWI) extracts hidden structure and physical properties of a volume of material by analyzing wave signals that traverse that volume. Applications include seismic exploration of the Earth’s subsurface, acoustic imaging and nondestructive testing (NDT) in material science, and ultrasound computed tomography (USCT) in medicine. Current approaches for solving CWI problems can be divided into two categories: those rooted in traditional physics and those based on deep learning. Physics-based methods stand out for their ability to provide high-resolution and quantitatively accurate estimates of acoustic properties within the medium. However, they can be computationally intensive and are susceptible to ill-posedness and nonconvexity typical of CWI problems. Machine learning (ML)-based computational methods have recently emerged, offering a different perspective to address these challenges. Diverse scientific communities have independently pursued the integration of deep learning in CWI. This review discusses how contemporary scientific ML techniques, and deep neural networks in particular, have been developed to enhance and integrate with traditional physics-based methods for solving CWI problems. We present a structured framework that consolidates existing research spanning multiple domains, including computational imaging, wave physics, and data science. This study concludes with important lessons learned from existing ML-based methods and identifies technical hurdles and emerging trends through a systematic analysis of the extensive literature on this topic.

42 ENGINEERING↗
DOE OSTI2026

Comprehensive earthquake catalogue update and spatiotemporal distribution analysis for Iraq and surrounding regions, northeastern Arabian Plate

The updated earthquake catalogue for Iraq covers the period from 1900 to the end of 2021 and includes over 37 000 recorded earthquakes. To create this comprehensive catalogue, five key steps were taken: compiling bulletins, calculating moment magnitudes, harmonizing magnitudes, establishing empirical conversion relations and evaluating the completeness of the catalogue. A notable enhancement in this update is the direct calculation of moment magnitudes for approximately 2800 earthquakes, achieved through the coda envelope technique and waveform data from the Mesopotamian Seismological Network (MPSN) in Iraq. This updated catalogue serves as a valuable resource for examining the spatiotemporal distribution of earthquakes, with respect to earthquake density, maximum moment magnitude and seismogenic depths. Additionally, the Gutenberg–Richter relationship was applied to calculate the a- and b-values specific to Iraq. The findings show that the Zagros Fold-Thrust Belt has a seismogenic layer (source) that ranges from 2 to 33 km deep and experiences high seismic activity. In contrast, the Mesopotamian Foredeep has a seismogenic layer ranging from 1 to 25 km deep and has lower seismic activity. The greatest seismic activity is concentrated around the Mandili-Badra-Teeb fault, which has experienced significant ruptures over time. The Outer Arabian Platform is identified as the main area of seismic activity, while additional activity occurs on the Inner Arabian Platform. Three major tectonic boundaries define the distribution of earthquakes in the northeastern Arabian Plate. These boundaries are defined by the Main Zagros Reverse Fault, the Zagros Foredeep Fault and the Anah Graben and Abu Jir-Euphrates Fault Zone. These boundaries highlight variations in seismicity levels and the spatial distribution of deformation in the region. The updated earthquake catalogue presented in this study is expected to play a vital role in regional seismicity assessments and seismic hazard analyses for Iraq and its surrounding areas.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

Abstract Subsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity ( V S ) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface V S variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred V S variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The V S variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

58 GEOSCIENCES↗