Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A New Decade in Seismoacoustics (2010–2022)

Several sources of interest often generate both low‐frequency acoustic and seismic signals due to energy propagation through the atmosphere and the solid Earth. Seismic and acoustic observations are associated with a wide range of sources, including earthquakes, volcanoes, bolides, chemical and nuclear explosions, ocean noise, and others. The fusion of seismic and acoustic observations contributes to a better understanding of the source, both in terms of constraining source location and physics, as well as the seismic to acoustic coupling of energy. Here in this review, we summarize progress in seismoacoustic data processing, including recent developments in open‐source data availability, low‐cost seismic and acoustic sensors, and large‐scale deployments of collocated sensors from 2010 to 2022. Similarly, we outline the recent advancements in modeling efforts for both source characteristics and propagation dynamics. Finally, we highlight the advantages of fusing multiphenomenological signals, focusing on current and future techniques to improve source detection, localization, and characterization efforts. This review aims to serve as a reference for seismologists, acousticians, and others within the growing field of seismoacoustics and multiphenomenology research.

58 GEOSCIENCES↗
DOE OSTI2023

Testing and Design of Discriminants for Local Seismic Events Recorded during the Redmond Salt Mine Monitoring Experiment

The Redmond Salt Mine (RSM) Monitoring Experiment in Utah was designed to record seismoacoustic data at distances less than 50 km for algorithm testing and development. During the experiment from October 2017 to July 2019, six broadband seismic stations were operating at a time, with three of them having fixed locations for the duration, whereas the three other stations were moved to different locations every one-and-half to two-and-half months. RSM operations consist of nighttime underground blasting several times per week. The RSM is located in proximity to a belt of active seismicity, allowing direct comparison of natural and anthropogenic sources. Using the recorded data set, we built 1373 events with local magnitude (M L ) of -2.4 and lower to 3.3. For 75 blasts (RMEs) from the Redmond Salt Mine and 206 tectonic earthquakes (EQs), both M L and the coda duration magnitude (MC) are well constrained. We used these events to test and design discriminants that separate the RMEs from the EQs and are effective at local distances. The discriminants consist of M L -M C , low-frequency Sg to high-frequency Sg, Pg/Sg phase-amplitude ratios, and Rg/Sg spectral amplitude ratios, as well as different combinations of two or more of these classifiers. The areas under the receiver operating characteristic curves (AUCs) of 0.92–1.0 for M L -M C , low-frequency Sg to high-frequency Sg, and Rg/Sg indicate that these discriminants are very effective. Conversely, the AUC of only 0.57 for Pg/Sg suggests that this discriminant is only slightly better than a random classifier. Among the effective classifiers, Rg/Sg, shows the lowest likelihood of misclassification (4.3%) for the populations. In conclusion, results of joint discriminant analyses suggest that even the arguably ineffective single classifier, like Pg/Sg in this case, can provide some value when used in combination with others.

58 GEOSCIENCES↗
DOE OSTI2024

Seismic Tomography 2024

Seismic tomography is the most abundant source of information about the internal structure of the Earth at scales ranging from a few meters to thousands of kilometers. It constrains the properties of active volcanoes, earthquake fault zones, deep reservoirs and storage sites, glaciers and ice sheets, or the entire globe. It contributes to outstanding societal problems related to natural hazards, resource exploration, underground storage, and many more. The recent advances in seismic tomography are being translated to nondestructive testing, medical ultrasound, and helioseismology. Nearly 50 yr after its first successful applications, this article offers a snapshot of modern seismic tomography. Focused on major challenges and particularly promising research directions, it is intended to guide both Earth science professionals and early-career scientists. The individual contributions by the coauthors provide diverse perspectives on topics that may at first seem disconnected but are closely tied together by a few coherent threads: multiparameter inversion for properties related to dynamic processes, data quality, and geographic coverage, uncertainty quantification that is useful for geologic interpretation, new formulations of tomographic inverse problems that address concrete geologic questions more directly, and the presentation and quantitative comparison of tomographic models. Finally, it remains to be seen which of these problems will be considered solved, solved to some extent, or practically unsolvable over the next decade.

58 GEOSCIENCES↗
DOE OSTI2024

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗
DOE OSTI2025

Effect of Time Window and Spectral Measurement Options on Empirical Green’s Function Analysis Using DAS Array and Seismic Stations

The recorded seismic waveform is a convolution of event source term, path term, and station term. Removing high-frequency attenuation due to path effect is a challenging problem. Empirical Green’s function (EGF) method uses nearly collocated small earthquakes to correct the path and station terms for larger events recorded at the same station. However, this method is subject to variability due to many factors. Here, we focus on three events that were well recorded by the seismic network and a rapid response distributed acoustic sensing (DAS) array. Using a suite of high-quality EGF events, we assess the influence of time window, spectral measurement options, and types of data on the spectral ratio and relative source time function (RSTF) results. Increased number of tapers (from 2 to 16) tends to increase the measured corner frequency and reduce the source complexity. Extended long time window (e.g., 30 s) tends to produce larger variability of corner frequency. The multitaper algorithm that simultaneously optimizes both target and EGF spectra produces the most stable corner-frequency measurements. The stacked spectral ratio and RSTF from the DAS array are more stable than two nearby seismic stations, and are comparable to stacked results from the seismic network, suggesting that DAS array has strong potential in source characterization.

58 GEOSCIENCES↗
DOE OSTI2025

Denoising Seismograms in the Time Domain Using a Deep Learning Model

Deep learning has emerged as a transformative tool for enhancing the extraction of reliable information from seismograms, addressing the increasing demand for precise and efficient seismic data analysis. We introduce an innovative encoder–decoder deep learning model, named WaveDenoiser, designed for noise reduction in the time domain, thereby eliminating the need for spectrogram computations that have been used for existing deep learning tools and significantly improving processing speed. Utilizing the benchmark dataset that is Stanford Earthquake Dataset, we developed three models of varying sizes: base, medium, and large. Notably, the large (referred to as WaveDenoiser) model demonstrated superior performance, achieving a median signal‐to‐noise ratio improvement of 8.8 dB on in‐distribution unseen data (in the same geographic region) and 7.7 dB on out‐distribution unseen data (in a new geographic region), outpacing both the base and medium models. Further evaluation of the WaveDenoiser model revealed a reduction in median arrival‐time errors by 0.02 s for P waves and 0.01 s for S waves when processing waveforms prior to phase picking using PhaseNet on in‐distribution unseen data. When tested on out‐distribution unseen data, the model also effectively reduced the P‐wave median arrival‐time error by 0.02 and 0.01 s in median arrival‐time error for S waves. Importantly, the application of WaveDenoiser resulted in a significant reduction of phase picking outliers by 1.1% to 3.6% for both P and S waves. In addition, we achieved over five times acceleration in processing speed compared with the seisBench implementation of DeepDenoiser. Our findings underscore the potential of WaveDenoiser as a powerful tool for improving seismic data analysis and processing efficiency.

P-waves↗
DOE OSTI2020

Effects of Shallow-Velocity Reductions on 3D Propagation of Seismic Waves

Abstract We perform 3D simulations of seismic wavefields to clarify effects of strong reductions of shallow velocities on long-period seismic waves. The simulations use a reference Community Velocity Model of southern California and a modified version with strong velocity reductions in the top 500 m of the Los Angeles basin. Differences between wavefields generated by 10 earthquakes in the reference and perturbed models are analyzed. Velocity changes are estimated by measuring relative time shifts between reference and perturbed seismograms using wavelet cross-correlation spectra. The results indicate that strong localized temporal velocity drops near the surface, such as those observed during strong ground motions, may generate regional perturbations of wavefields at periods up to 20 s. These perturbations may be misinterpreted as generated by temporal changes at seismogenic depths. The results also have important implications for waveform tomography studies.

Geochemistry & Geophysics↗
DOE OSTI2021

Analysis of Seismic Signals Generated by Vehicle Traffic with Application to Derivation of Subsurface Q-Values

Correct identification and modeling of anthropogenic sources of ground motion are of considerable importance for many studies, including detection of small earthquakes and imaging seismic properties below the surface. Here, to understand signals generated by common vehicle traffic, we use seismic data recorded by closely spaced geophones normal to roads at two sites on San Jacinto fault zone. To quantify the spatiotemporal and frequency variations of the recorded ground motions, we develop a simple analytical solution accounting for propagation and attenuation of surface waves. The model reproduces well-observed bell-shaped spectrograms of car signals recorded by geophones close to roads, and it can be used to estimate frequency-dependent Q-values of the subsurface materials. The data analysis indicates Q-values of 3–40, for frequencies up to 150 Hz for road-receiver paths at the two examined sites. The derived Q-values are consistent with attenuation factors of surface waves previously obtained with other methods. The analytical results and analysis procedure provide a highly efficient method for deriving Q-values of shallow subsurface materials.

58 GEOSCIENCES↗
DOE OSTI2021

Wastewater Disposal Has Not Significantly Altered the Regional Stress State in Southern Kansas

Wastewater disposal is primarily responsible for the increased seismicity rate since ~2013 in southern Kansas. Previous work that used shear-wave splitting (SWS) in southern Kansas interpreted an ~90° temporal rotation in the fast polarization direction and attributed it to increased pore pressures resulting from fluid injection. However, this interpreted rotation coincided with a change in the stations used to make the SWS measurements. We investigate the temporal variability of fast azimuths in southern Kansas by making SWS measurements on earthquake families with similar source–receiver paths recorded on a stable local seismic network. We select high-quality SWS measurements by investigating the stability of results across 65 different frequency bands between 0.5 and 15 Hz. We find that the fast polarization direction in southern Kansas is relatively constant with an average east-northeast (~N79°E) orientation between 2014 and 2017. Our fast polarization measurements are primarily a reflection of the maximum principal horizontal stress direction (SHmax). We observe a slight spatial change in SHmax to the northeast (~N55°E) near the Nemaha ridge in Oklahoma. However, we do not observe any significant temporal rotation of SHmax or variation in delay time (i.e., crack density) in southern Kansas, contrary to the earlier study. The previously interpreted ~90° rotation may either be a reflection of a very local stress change or a misinterpretation of SWS results potentially due to the use of inconsistent source–receiver paths. Our SWS measurements cover the period of peak wastewater disposal and seismicity rates and suggest an absence of significant temporal rotations in the local anisotropy and stress orientations associated with wastewater disposal.

Geochemistry & Geophysics↗
DOE OSTI2021

Exploring the Effects of Emplacement Conditions on Explosion P/S Ratios across Local to Regional Distances

High-frequency (~> 2 Hz) seismic P/S amplitude ratios are well-established as a discriminant to distinguish between natural earthquakes and underground explosions at regional distances (~200–1500 km). As research shifts toward identifying lower-yield events, work has begun to investigate the potential of this discriminant for use at local distances (<200 km), in which initial results raise questions about its effectiveness. Here, we utilize data from several chemical explosion experiment series at the Nevada National Security Site in southern Nevada in the United States to study explosion Pg/Lg ratios across the range of local to regional distances. Additionall, the experiments are conducted over differing emplacement conditions, with contrasting geologies and a variety of yields and depths of burial, including surface explosions. We first establish the similarities of Pg/Lg ratios from chemical explosions to those from historic nuclear tests and conclude that, as previous data have suggested, chemical explosion ratios are good proxies for nuclear tests. We then examine Pg/Lg ratios from the new experiment series as functions of distance, yield, depth of burial, and scaled depth of burial (SDOB). At far-local and regional distances, we observe consistently higher ratios from hard-rock explosions compared to ones in a weaker dry alluvium medium, consistent with prior regional distance results. No other trends with yield, depth of burial, or SDOB are strongly evident. Scatter in the observed ratios is very high, particularly at the shortest event-to-station distances, suggesting that small-scale path effects play a significant role. On average, the local distance explosion Pg/Lg ratios show remarkable consistency across all the variations in emplacement. Explosion source models will need to reproduce these results.

58 GEOSCIENCES↗
DOE OSTI2022

Automatic Waveform Quality Control for Surface Waves Using Machine Learning

Surface-wave seismograms are widely used by researchers to study Earth’s interior and earthquakes. To extract information reliably and robustly from a suite of surface waveforms, the signals require quality control screening to reduce artifacts from signal complexity and noise. This process has usually been completed by human experts labeling each waveform visually, which is time consuming and tedious for large data sets. We explore automated approaches to improve the efficiency of waveform quality control processing by investigating logistic regression, support vector machines, K-nearest neighbors, random forests (RF), and artificial neural networks (ANN) algorithms. To speed up signal quality assessment, we trained these five machine learning (ML) methods using nearly 400,000 human-labeled waveforms. The ANN and RF models outperformed other algorithms and achieved a test accuracy of 92%. We evaluated these two best-performing models using seismic events from geographic regions not used for training. The results show that the two trained models agree with labels from human analysts but required only 0.4% of the time. Although the original (human) quality assignments assessed general waveform signal-to-noise, the ANN or RF labels can help facilitate detailed waveform analysis. Our investigations demonstrate the capability of the automated processing using these two ML models to reduce outliers in surface-wave-related measurements without human quality control screening.

58 GEOSCIENCES↗
DOE OSTI2024

Assessing the Expansion of Ground-Motion Sensing Capability in Smart Cities via Internet Fiber-Optic Infrastructure

Monitoring ground motion in smart cities can improve the public safety by providing critical insights on natural and anthropogenic hazards, for example, earthquakes, landslides, explosions, infrastructure failures, and so forth. Although seismic activity is typically measured using dedicated point sensors (e.g., geophones and accelerometers), techniques such as distributed acoustic sensing have demonstrated the utility of using fiber-optic cable to detect seismic activity over comparable distances. In this article, we present the results of a study that quantifies the expansion in an area monitored for low-amplitude ground-motion events by augmenting existing point sensors with the internet fiber-optic cable infrastructure. Here we begin by describing our methodology, which utilizes geospatial data on point sensors and internet optical fiber deployed in metropolitan statistical areas (MSAs) in the United States. We extend these data to identify the area that can be monitored by (1) considering the observed seismic noise data in target locations, (2) applying the model from Wilson et al. (2021) to understand the potential coverage area gains using optical fiber sensing, and (3) optimizing the selection of fiber segments to maximize coverage and minimize deployment costs. We implement our methodology in ArcGIS to assess the additional area that can be monitored for low-amplitude ground-motion events (i.e., magnitude >0.5) by utilizing internet fiber-optic cables in the 100 most populous MSAs in the United States. We find that the addition of internet fiber-based sensors in MSAs would increase the area monitored on average by over an order of magnitude from 1% to 12%, if the subset of fiber cable segments that maximize coverage and minimize deployment costs is chosen even if only 20% of all fibers are used.

58 GEOSCIENCES↗
DOE OSTI2025

Rapid and Temporary Deployment of Fiber-Optic Arrays to Measure Surface Chemical Explosion Series with Distributed Acoustic Sensing

Rapid and temporary distributed acoustic sensing (DAS) deployments are crucial for accurately capturing the seismic wavefield following major events, such as earthquakes, or in anticipation of known events of interest, such as chemical explosions. We provide an overview of two DAS campaigns conducted in May and October 2024 to record the seismoacoustic waves generated by two series of surface chemical explosions in New Mexico, United States, involving 1- and 10-ton trinitrotoluene-equivalent charges. In both campaigns, we deployed approximately 2 km of fiber-optic cables in a dry riverbed, about 12 km east of the explosion sites. For the October campaign, seven geophones and two anemometers were collocated with the fiber. We describe the field deployments and present preliminary results from the recorded signals. Specifically, we cross analyze the data recorded in May and October, validate the DAS data with geophone measurements, and highlight the benefits of burying fibers to reduce wind noise and improve signal-to-noise ratios. This study demonstrates the potential of DAS to record the seismoacoustic waves generated by surface chemical explosions of varying sizes at distances of approximately 10 km from the source.

58 GEOSCIENCES↗
DOE OSTI2025

Applying Machine Learning and Bayesian Inference to Identify and Locate Moving Anthropogenic Sources Using Distributed Acoustic Sensing Data

Distributed acoustic sensing (DAS) systems, which use existing telecommunication fibers, offer high‐resolution capabilities ideal for recording anthropogenic sources. However, the complexity of urban environments and the large amount of data recorded by DAS require automated methods to efficiently detect and categorize anthropogenic sources. Here, we evaluate how well three machine learning models (k‐nearest neighbor [k‐NN], convolutional neural networks, and recurrent‐convolutional neural networks) can identify various anthropogenic sources recorded by DAS. Our findings reveal that both k‐NN and neural network methods perform well in high signal‐to‐noise ratio (SNR) settings. However, their accuracy decreases at SNRs <4. We also use Kalman filtering, a form of Bayesian inference, on backprojected locations of these sources to recover locations that generally fall within standard smartphone Global Positioning System errors. By combining machine learning and Kalman filter results, we calculate a multidimensional model of moving anthropogenic sources. These results demonstrate the potential of DAS data in urban seismology for accurately identifying and locating such sources. Depending on the research objectives, these sources can be further studied or filtered out to improve the quality of seismic data for earthquake studies. Such methods provide a valuable tool for urban seismology and seismic hazard analysis.

Luckie, Thomas William [Sandia National Laboratori↗
DOE OSTI2022

Multi-Event Explosive Seismic Source for the 2022 M w 6.3 Hunga Tonga Submarine Volcanic Eruption

The eruption of the Hunga Tonga–Hunga Ha’apai submarine volcano on 15 January 2022 produced a variety of geophysical responses, including a significant seismic signal. We study the seismic source process of this event by inverting for moment tensors (MTs) using regional surface waves (Rayleigh, Love). By comparing inversion results for the eruption with eight nearby earthquakes, we show that it is possible to discriminate MT source types. Our inversion yields a shallow explosive source for the eruption and reveals the importance of trade-offs among depth, magnitude, and source type. We illustrate these trade-offs by representing the misfit variations over the eigenvalue lune. Finally, we invert for the source-time function of the sequence of explosions that occurred in the first minutes of the eruption. The multi-event source-time function comprises four subevents spanning ~270 s, with a total magnitude estimate of M w 6.34 ± 0.10.

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 OSTI2024

Making Phase-Picking Neural Networks More Consistent and Interpretable

Improving the interpretability of phase-picking neural networks remains an important task to facilitate their deployment to routine, real-time seismic monitoring. The popular phase-picking neural networks published in the literature lack interpretability because their output prediction scores do not necessarily correspond with the reliability of phase picks and can even be highly inconsistent depending on how we window the waveform data. Here, we show that systematically shifting the waveforms during training and using an antialiasing filter within the neural network architecture can substantially improve the consistency of the output prediction scores and can even make them scale with the signal-to-noise ratios of the waveforms. We demonstrate the improvements by applying these approaches to a commonly used phase-picking neural network architecture and using waveform data from the 2019 Ridgecrest earthquake sequence.

58 GEOSCIENCES↗
DOE OSTI2025

Regional Source-Type Discrimination Using Nonlinear Alignment Algorithms

The discrimination problem in seismology aims to accurately classify different underground source types based on local, regional, and/or teleseismic observations of ground motion. Typical discriminant approaches are rooted in fundamental, physics-based differences in radiation pattern or wave excitation, which can be frequency-dependent and may not make use of the full waveform. In this article, we explore whether phase and amplitude distances derived from dynamic time warping (DTW) and elastic shape analysis (ESA) can inform event discrimination. We demonstrate the ability to distinguish underground point sources using synthetic waveforms calculated for a 1D Earth model and various source mechanisms. We then apply the method to recorded data from events in the Korean Peninsula, which includes declared nuclear explosions, a collapse event, and naturally occurring earthquakes. Phase and amplitude distances derived from DTW and ESA are then used to classify the event types via dendrogram and k-nearest-neighbor clustering analyses. Using information from the full waveform, we show how different underground sources can be distinguished at regional distances. We highlight the potential of these nonlinear alignment algorithms for discrimination and comment on ways we can extend the framework presented here.

58 GEOSCIENCES↗