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 109 records · Page 6

DeepShake: Shaking Intensity Prediction Using Deep Spatiotemporal RNNs for Earthquake Early Warning

Abstract We propose a deep spatiotemporal recurrent neural network, DeepShake, to project future shaking intensity directly from current ground-motion observations. DeepShake is a network-based forecasting model, able to predict future shaking intensity at all stations within a network given previously measured ground shaking. The model is not given any a priori knowledge of station locations; instead, it learns wave propagation amplitudes and delays solely from training data. We developed DeepShake with the 35,679 earthquakes from the 2019 Ridgecrest sequence. Tasked with alerting for modified Mercalli intensity (MMI) IV+ shaking on 3568 validation earthquakes at least 5 s in advance, DeepShake achieves an equal error rate of 11.4%. For the Mw 7.1 earthquake that hit Ridgecrest on 5 July 2019, DeepShake was able to provide targeted alerts to all stations inside the network 5 s prior to the arrival of MMI IV+ waveforms. DeepShake demonstrates that deep spatiotemporal neural networks can effectively provide one-step earthquake early warning with reasonable accuracy and latency.

Geochemistry & Geophysics↗

Central Asia Seismic Hazard Assessment (Casha) "Complete" Catalog of Earthquakes

Probabilistic Seismic Hazard Assessments (PSHA) greatly benefit from reliable earthquake catalogs. Central Asia, while rich in seismic monitoring history, reports a fraction of their earthquake data to international data collection centres. Furthermore, since the collapse of Soviet Union until now, monitoring and data collection efforts in Central Asia have not been centralized, making monitoring and reporting uneven and disconnected in this time period. Lawrence Livermore National Laboratory (LLNL) has undertaken a multi-year effort to systematically work with national seismic monitoring and data centres and seismic networks in Kazakhstan, Kyrgyzstan, and Tajikistan to compile a unified seismological bulletin for Central Asia. In tandem with this effort, an earthquake catalog that is processed to be PSHA ready was also generated. The processing for this effort included further duplicate identification and removal, magnitude harmonization and catalog completeness analyses. This report describes the “complete” PSHA-ready earthquake catalog, including the completeness intervals that were applied.

58 GEOSCIENCES↗

EQSIM and RAJA: Enabling Exascale Predictions of Earthquake Effects on Critical Infrastructure

Nearly 120 years ago, the great “San Francisco” earthquake of 1906 provided a stark and sobering view of the havoc that can be caused by the sudden and violent movement of Earth’s tectonic plates. According to USGS, the rupture along the San Andreas fault extended 296 miles (447 kilometers) and shook so violently that the motion could be felt as far north as Oregon and east into Nevada. The estimated 7.9-magnitude quake and subsequent fires decimated the major metropolis and surrounding areas: buildings turned to ruins, hundreds of thousands of people left homeless, and a death toll exceeding 3,000. Today, as evidenced by the catastrophic 7.8-magnitude earthquake that struck Turkey in February of 2023, large earthquakes still present a significant potential danger to life and economic security as researchers work to develop ways to better understand earthquake phenomena and quantify associated risks.

42 ENGINEERING↗

Infrasound Detections of Low-Magnitude Earthquakes: Preliminary Results of the West Texas Acoustic Experiment

Infrasound observations have grown increasingly important for the monitoring of earthquakes. While large earthquakes generate infrasound that can be detected thousands of kilometers away, there are few near-field observations of infrasound generated by low-magnitude events. We describe preliminary results of the West Texas Acoustic Experiment, during which infrasound sensors collected continuous data in the Permian Basin for a six-month period spanning January—June 2023. During this time, more than 1000 earthquakes with magnitudes between 1.2 and 4.2 occurred within 50 km of the network. We used spectral analysis, array processing, and manual inspection of waveforms to evaluate arrivals of infrasound signals following 84 events with magnitudes between 2.5 and 4.2. Here, we describe eight such events and the infrasound signals associated with each. We find detections of seismic-to-acoustic infrasound signals associated with seven events. We also find strong evidence of a laterally-propagating, purely acoustic wave generated by an M2.9 earthquake.

58 GEOSCIENCES↗

Probabilistic Seismic Hazard Analysis for Iraq Based on the Updated Earthquake Catalog (1900-2021) and Ground Motion Characteristics

Onur et al. (2017) compiled the first comprehensive earthquake catalog for Iraq, covering 1900 to 2009 within 26°–40°N latitude and 36°–51°E longitude. This catalog was utilized in a probabilistic seismic hazard assessment (PSHA) by Abdulnaby et al. (2020) to aid in updating Iraq’s building code seismic provisions. Recently, we have updated the earthquake catalog for Iraq by adding earthquakes recorded from 2010 to 2021 and directly calculating moment magnitude (Mw) for about 2,800 earthquakes using the coda envelope methodology and waveform data from the Mesopotamian Seismological Network (MPSN) in Iraq.

58 GEOSCIENCES↗

Preseismic Fault Creep and Elastic Wave Amplitude Precursors Scale With Lab Earthquake Magnitude for the Continuum of Tectonic Failure Modes

Tectonic faults fail in a continuum of modes from slow earthquakes to elastodynamic rupture. Precursory variations in elastic wavespeed and amplitude, interpreted as indicators of imminent failure, have been observed in limited natural settings and lab experiments where they are thought to arise from contact rejuvenation and microcracking within and around the fault zone. However, the physical mechanisms and connections to fault creep are poorly understood. Here we vary loading stiffness during frictional shear to generate a range of slip modes and measure fault zone properties using transmitted elastic waves. We find that elastic wave amplitudes show clear changes before fault failure. The temporal onset of amplitude reduction scales with lab earthquake magnitude and the magnitude of this reduction varies with fault slip. Our data provide clear evidence of precursors to lab earthquakes and suggest that continuous seismic monitoring could be useful for assessing fault state and seismic hazard potential.

58 GEOSCIENCES↗

Generalization of Deep-Learning Models for Classification of Local Distance Earthquakes and Explosions across Various Geologic Settings

Although accurately classifying signals from earthquakes and explosions at local distance (<250 km) remains an important task for seismic network operations, the growing volume of available seismic data presents a challenge for analysts using traditional source discrimination techniques. In recent years, deep-learning models have proven effective at discriminating between low-magnitude earthquakes and explosions measured at local distances, but it is not clear how well these models are capable of generalizing across different geological settings. To address the issue of generalization between regions, we train deep-learning models (convolutional neural networks [CNNs]) on time–frequency representations (scalograms) of three-component earthquake and explosion signals from eight different regions in the continental United States. We explore scenarios where models are trained on data from all regions, individual regions, or all but one region. We find that although CNN models trained on individual regions do not necessarily generalize well across different settings, models trained on multiple regions that include diverse path coverage generalize to new regions, with station-level accuracy of up to 90% or more for data sets from unseen regions. In general, CNN-based discrimination models significantly outperform models based on uncorrected P/S ratio (measured in the 10–18 Hz frequency band), even when CNN models are tested on data from entirely unseen regions.

58 GEOSCIENCES↗

Impacts of a Cascadia Subduction Zone Earthquake on Water Levels and Wetlands of the Lower Columbia River and Estuary

Abstract Subsidence after a subduction zone earthquake can cause major changes in estuarine bathymetry. Here, we quantify the impacts of earthquake‐induced subsidence on hydrodynamics and habitat distributions in a major system, the lower Columbia River Estuary, using a hydrodynamic and habitat model. Model results indicate that coseismic subsidence increases tidal range, with the smallest changes at the coast and a maximum increase of ∼10% in a region of topographic convergence. All modeled scenarios reduce intertidal habitat by 24%–25% and shifts ∼93% of estuarine wetlands to lower‐elevation habitat bands. Incorporating dynamic effects of tidal change from subsidence yields higher estimates of remaining habitat by multiples of 0–3.7, dependent on the habitat type. The persistent tidal change and chronic habitat disturbance after an earthquake poses strong challenges for estuarine management and wetland restoration planning, particularly when coupled with future sea‐level rise effects.

54 ENVIRONMENTAL SCIENCES↗

Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array

SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.

58 GEOSCIENCES↗

Performance-based earthquake early warning for tall buildings

The ShakeAlert Earthquake Early Warning (EEW) system aims to issue an advance warning to residents on the West Coast of the United States seconds before the ground shaking arrives, if the expected ground shaking exceeds a certain threshold. However, residents in tall buildings may experience much greater motion due to the dynamic response of the buildings. Therefore, there is an ongoing effort to extend ShakeAlert to include the contribution of building response to provide a more accurate estimation of the expected shaking intensity for tall buildings. Currently, the supposedly ideal solution of analyzing detailed finite element models of buildings under predicted ground-motion time histories is not theoretically or practically feasible. The authors have recently investigated existing simple methods to estimate peak floor acceleration (PFA) and determined these simple formulas are not practically suitable. Instead, this article explores another approach by extending the Pacific Earthquake Engineering Research Center (PEER) performance-based earthquake engineering (PBEE) to EEW, considering that every component involved in building response prediction is uncertain in the EEW scenario. Additionally, while this idea is not new and has been proposed by other researchers, it has two shortcomings: (1) the simple beam model used for response prediction is prone to modeling uncertainty, which has not been quantified, and (2) the ground motions used for probabilistic demand models are not suitable for EEW applications. In this article, we address these two issues by incorporating modeling errors into the parameters of the beam model and using a new set of ground motions, respectively. We demonstrate how this approach could practically work using data from a 52-story building in downtown Los Angeles. Using the criteria and thresholds employed by previous researchers, we show that if peak ground acceleration (PGA) is accurately estimated, this approach can predict the expected level of human comfort in tall buildings.

58 GEOSCIENCES↗

A Complete Machine-Learning-Based Workflow to Illuminate Earthquake Processes

Under this grant we developed, tested, and made available, machine learning models to improve the tasks in the earthquake monitoring workflow (Figure 1). We implemented these models as a part of an end-to-end workflow for seismic network processing and demonstrated, in a variety of settings, that these methods generalize and that they result in dramatically more comprehensive earthquake catalogs. These catalogs illuminate earthquake processes in detail and to an extent that had previously not been possible, and they do so for both tectonic seismicity and seismicity induced by fluid injection related to unconventional hydrocarbon development. This report summarizes the results from the 15 publications that resulted from this grant. Those contributes are divided into: (1) the development of specific tasks related to monitoring (7 publications), (2) the organization of those tasks into workflows for seismic monitoring (2 publications), and (3) applications of those workflows to data (6 publications).

58 GEOSCIENCES↗

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARY Observational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

, RSQSim↗

Probing Slow Earthquakes With Deep Learning

Abstract Slow earthquakes may trigger failure on neighboring locked faults that are stressed sufficiently to break, and slow slip patterns may evolve before a nearby great earthquake. However, even in the clearest cases such as Cascadia, slow earthquakes and associated tremor have only been observed in intermittent and discrete bursts. By training a convolutional neural network to detect known tremor on a single seismic station in Cascadia, we isolate and identify tremor and slip preceding and following known larger slow events. The deep neural network can be used for the detection of quasi‐continuous tremor, providing a proxy that quantifies the slow slip rate. Furthermore, the model trained in Cascadia recognizes tremor in other subduction zones and also along the San Andreas Fault at Parkfield, suggesting a universality of waveform characteristics and source processes, as posited from experiments and theory.

58 GEOSCIENCES↗

Acoustic Energy Release During the Laboratory Seismic Cycle: Insights on Laboratory Earthquake Precursors and Prediction

Abstract Machine learning can predict the timing and magnitude of laboratory earthquakes using statistics of acoustic emissions. The evolution of acoustic energy is critical for lab earthquake prediction; however, the connections between acoustic energy and fault zone processes leading to failure are poorly understood. Here, we document in detail the temporal evolution of acoustic energy during the laboratory seismic cycle. We report on friction experiments for a range of shearing velocities, normal stresses, and granular particle sizes. Acoustic emission data are recorded continuously throughout shear using broadband piezo‐ceramic sensors. The coseismic acoustic energy release scales directly with stress drop and is consistent with concepts of frictional contact mechanics and time‐dependent fault healing. Experiments conducted with larger grains (10.5 μm) show that the temporal evolution of acoustic energy scales directly with fault slip rate. In particular, the acoustic energy is low when the fault is locked and increases to a maximum during coseismic failure. Data from traditional slide‐hold‐slide friction tests confirm that acoustic energy release is closely linked to fault slip rate. Furthermore, variations in the true contact area of fault zone particles play a key role in the generation of acoustic energy. Our data show that acoustic radiation is related primarily to breaking/sliding of frictional contact junctions, which suggests that machine learning‐based laboratory earthquake prediction derives from frictional weakening processes that begin very early in the seismic cycle and well before macroscopic failure.

58 GEOSCIENCES↗

The Spatiotemporal Evolution of Granular Microslip Precursors to Laboratory Earthquakes

Laboratory earthquake experiments provide important observational constraints for our understanding of earthquake physics. Here we leverage continuous waveform data from a network of piezoceramic sensors to study the spatial and temporal evolution of microslip activity during a shear experiment with synthetic fault gouge. We combine machine learning techniques with ray theoretical seismology to detect, associate, and locate tens of thousands of microslip events within the gouge layer. Microslip activity is concentrated near the center of the system but is highly variable in space and time. While microslip activity rate increases as failure approaches, the spatiotemporal evolution can differ substantially between stick-slip cycles. These results illustrate that even within a single, well-constrained laboratory experiment, the dynamics of earthquake nucleation can be highly complex.

58 GEOSCIENCES↗

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗

The First Detection of an Earthquake From a Balloon Using Its Acoustic Signature

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground-based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon-borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo-acoustic phenomenology demonstrates that seismic activity is detectable from a high-altitude platform on Earth, and that Rayleigh wave-induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus.

58 GEOSCIENCES↗

Seismic Imaging of the Mw 7.1 Ridgecrest Earthquake Rupture Zone From Data Recorded by Dense Linear Arrays

Abstract We analyze seismograms recorded by four arrays (B1–B4) with 100 m station spacing and apertures of 4–8 km that cross the surface rupture of the 2019 Mw 7.1 Ridgecrest earthquake. The arrays extend from B1 in the northwest to B4 in the southeast of the surface rupture. Delay times between P wave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and local P waves travel faster on the northeast than the southwest side of the fault beneath arrays B1 and B4, but the velocity contrast is less reliably resolved at arrays B2 and B3. We identify several 1–2 km wide low‐velocity zones with much slower inner cores that amplify S waveforms, inferred as damage zones, beneath each array. The damage zones at arrays B2 and B4 also generate fault‐zone head and trapped waves. An automated detector, based on peak ground velocities and durations of high‐amplitude waves, identifies candidate fault‐zone trapped waves (FZTWs) in a localized zone for ∼600 earthquakes at array B4. Synthetic waveform modeling of averaged FZTWs, generated by ∼30 events with high‐quality signals, indicates that the trapping structure at array B4 has a width of ∼300 m, depth of 3–5 km, S wave velocity reduction of ∼20% with respect to the surrounding rock, Q ‐value of ∼30, and S wave velocity contrast of ∼4% across the fault (faster on the northeast side). The results show complex fault‐zone internal structures (velocity contrasts and low‐velocity zones) that vary along fault strike.

Qiu, Hongrui↗