Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “earthquake”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An open-access simulated earthquake ground-motion database for an M7 Hayward Fault earthquake in the San Francisco Bay Region

Comprehensive understanding of earthquake ground motions, particularly in the near-fault region of large-magnitude events, is limited by gaps in strong-motion data. This challenge is prominent in areas with high seismic hazard but infrequent large earthquakes where data is sparse and difficult to interpret. These data limitations lead to uncertainties in the development of site-specific ground motions, which are crucial for engineering risk assessments. To address these challenges, physics-based regional-scale ground-motion simulations have been developed. With the emergence of exaflop-scale computing ecosystems, it is now possible to simulate regional earthquake processes at unprecedented fidelity and generate the large number of fault rupture realizations necessary to characterize both intra- and inter-event ground-motion variability. This article introduces a new database of simulated earthquake ground motions, created for applications in earthquake engineering, earthquake planning, and emergency response. The inaugural version of the database features simulated ground motions for a magnitude 7 Hayward Fault earthquake in the San Francisco Bay Region (SFBR), using the EarthQuake SIMulation (EQSIM) simulation framework and the Graves–Pitarka kinematic rupture model. The aim is to provide high-fidelity, spatially dense, three-component motions generated on the Department of Energy’s (DOE) newest generation of graphics processing unit (GPU)-accelerated supercomputers. These motions are being made openly available to the engineering, scientific, and disaster planning communities. In addition, this work develops protocols for the efficient dissemination of these large data sets and emphasizes community engagement to build confidence in their application. This article discusses the methodology behind the data, underlying software verification and validation, scalable data management, and a user interface for data access. The goal is to facilitate widespread use and elicit expert feedback to maximize the utility and exploitation of simulated motions. While the initial focus is on the San Francisco Region, simulations for additional regions will be added as the DOE program progresses.

Simulated ground-motion database↗
DOE OSTI2026

Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes

This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.

EQSIM↗
DOE OSTI2025

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗
DOE OSTI2025

Exploring Uncertainty in Moment Estimation for Small Earthquakes in Southern Nevada Using the Coda Envelope Method

Compiling source parameter estimates for small earthquakes is important both for our understanding of earthquake physics and for accurately assessing earthquake hazard. Reliable source parameter estimates are difficult to achieve for small earthquakes, in part due to our inability to accurately model the relevant physical processes at high frequencies. The coda envelope methodology developed by Mayeda and Walter (1996) and Mayeda et al. (2003) can mitigate this concern and estimate the moment of small earthquakes by determining the parameters that control the shape of the S-wave coda envelope while eliminating path effects by minimizing the scatter between seismic stations. Here, we use an open-source implementation of this technique called the Coda Calibration Tool (CCT; Barno, 2017) to calculate CCT-based moment magnitude estimates of small earthquakes (M L 0–3) in the Rock Valley, Nevada, region within the Nevada National Security Site. The Rock Valley data set is of particular interest because it allows us to explore the changes in uncertainties of the coda calibration method with earthquake size and depth. We found that a consistent linear relationship exists between the local magnitude M L and our coda-derived M w estimates for earthquakes as small as M L 0–3, but that current CCT workflows do not accurately characterize very shallow events. We also demonstrate that the epistemic uncertainty in the apparent stress value assumed by the CCT algorithm can influence magnitude estimates of small earthquakes. In conclusion, these results provide valuable insight into the seismicity of this region, and inform future analysis and modeling efforts for nuclear monitoring and seismic hazard.

58 GEOSCIENCES↗
DOE OSTI2025

Fiber-Optic Sensing for Earthquake Hazards Research, Monitoring, and Early Warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

58 GEOSCIENCES↗
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 OSTI2025

Fault Network Geometry Modulates Earthquake Source Spectra Across Scales

Earthquake source spectra provide unique insights into the earthquake rupture process. Motivated by previous research suggesting that complex fault geometries enhance high‐frequency seismic radiation, we study the influence of fault network geometry on earthquake source spectra using multiple independent observations. At regional scales, we examine correlations of stress drop measurements with surface fault trace misalignment in Southern California, Japan, and Central Italy. At a global scale, we examine correlations of moment‐rate function complexity of large earthquakes with focal mechanism variability, a proxy for local fault complexity. Despite significant scatter in the observations, we find overall consistent positive correlations. The concept that elastic interactions of discrete fault structures during the earthquake rupture process generates high‐frequency ground motions offers a coherent framework for interpreting our observations. These findings suggest that variations in fault complexity explain why some earthquakes produce stronger high‐frequency ground motions than others.

Lee, Jaeseok [Brown Univ., Providence, RI (United ↗
DOE OSTI2026

Learning earthquake ground motions via conditional generative modeling

Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.

geophysics↗
DOE OSTI2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

Science & Technology - Other Topics↗
DOE OSTI2026

How earthquakes organize stress

Stress is not uniform in the Earth. Therefore, we must use natural experiments to measure the distribution of stresses and related quantities, rather than single values. For instance, dynamic triggering shows that faults are uniformly distributed over their loading cycles in Southern California. The probability that a fault ruptures across a barrier measures the in situ energy distribution. Fault roughness reflects the distribution of strength. These natural experiments produce observable distributions that are surprisingly consistent and suggest some degree of self-organization in the Earth’s crust. Once established, the functional form of the distributions can be used to track changes in response to earthquakes as well as to distinguish fundamentally different fault systems. Transient fault locking before stress release in laboratory experiments can be interpreted as a consequence of self-organization of fault stress. The robust self-organization of multiple variables in earthquake systems suggests that the most consequential mechanical outcome of earthquakes may be the redistribution of stress and the strain energy associated with it. The low friction on a fault during seismic slip as inferred by temperature measurements of the Tohoku earthquake is consistent with dissipation playing a secondary role to this redistribution process. Through stress redistribution and interaction, subduction zone faults tend to synchronize, perhaps due to their geometric simplicity, while the continental system of Southern California cannot synchronize, perhaps due to the complexity of the fault network. Earthquakes organize stress in the crust and produce a suite of well-defined, consistent distributions.

earthquakes↗
DOE OSTI2025

Earthquake detection in a simulated lunar regolith using distributed acoustic sensing

Current models of inner lunar geology have largely been inferred from the seismic experiments and observations performed during the Apollo missions that comprised a relatively small number of seismic instruments. Refining constraints on fundamental lunar relationships such as crust-mantle and mantle-core boundaries in the future will require seismic arrays spanning larger epicentral distances. A promising technology for installing dense seismic arrays with minimal human effort is distributed acoustic sensing (DAS), an approach that allows a single length of fiber optic cable to act as hundreds or thousands of sensors when coupled with a DAS interrogator. While terrestrial uses of DAS technology for seismic monitoring rely on burying the cable to maximize fidelity of seismic signal transmission to the fiber, digging meters of trench to bury optical fiber on lunar or planetary surfaces is logistically infeasible. To evaluate DAS signal attenuation due to surface deployment of cable in lunar regolith, we completed earthquake detection analyses that evaluated the sensitivity of an optic-fiber DAS system to seismic signals at different burial depths. We deployed a single-mode fiber in a 10-m open-bottom wooden box filled with a lunar regolith simulant (LRS) with fiber buried at different depths within the LRS and recorded signals for four regional and local earthquakes. The results were used to identify and evaluate signal attenuation in surface-deployed fiber compared to buried fiber in the LRS. Burial depth responses to active-source signals were also evaluated similar to previous studies characterizing DAS sensitivity of surface-deployed fiber. Atmospheric noise was minimal as the cable was deployed in an indoor environment; however, where observed, atmospheric and anthropogenic noise was filtered out using the same bandpass filtering used to identify earthquake events. We found that signal attenuation of the surface-deployed fiber compared to buried fiber was relatively high in active-source experiments but was not consistently observed in earthquake signals. That burial depth is not highly correlated to attenuation of the observed earthquake signals indicates that in a noise-limited environment, placing DAS-interrogated fiber directly at the regolith surface may be a promising deployment strategy to consider for sensing remote seismic signals during lunar exploration.

58 GEOSCIENCES↗
DOE OSTI2026

Constructing a High‐Resolution Aftershock Catalog for the 2017 Mw 8.2 Tehuantepec Earthquake Sequence Using a Machine Learning–Based Workflow

The 8 September 2017 Mw 8.2 Tehuantepec earthquake was the largest instrumentally recorded normal‐faulting earthquake in Mexico. The mainshock occurred offshore within the Tehuantepec seismic gap, generating >30,000 aftershocks in the following year. We applied an open‐source, machine learning (ML)–assisted workflow to construct a high‐resolution aftershock catalog using data from temporary and permanent seismic networks in southern Mexico. The workflow integrates PhaseNet for phase detection; GaMMA for phase association; and VELEST, HypoInverse, and HypoDD for velocity modeling and relocation. We processed seven months of continuous waveform data from 29 broadband stations, including a temporary rapid‐response deployment that improved station coverage of the offshore rupture zone. To evaluate performance, we compared our results against analyst‐reviewed picks and event locations from the Servicio Sismológico Nacional catalog. The resulting catalog contains 11,374 relocated earthquakes and represents the most comprehensive published dataset for this sequence, incorporating the first full use of the temporary network. Relocated hypocenters show improved depth control and align well with the Slab2.0 subduction geometry, revealing clearer separation between offshore slab events and onshore crustal seismicity. This study demonstrates that combining ML‐based detection with established methods provides a scalable and reproducible approach for constructing high‐quality earthquake catalogs in tectonically complex environments and offers practical guidance for adapting similar workflows to other earthquake sequences.

Garcia, Marc [The University of Texas at El Paso, ↗
DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗
DOE OSTI2026

Limited Dynamic Earthquake Triggering in Nevada

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

58 GEOSCIENCES↗
DOE OSTI2025

Real-time reconstruction of ground motion during small magnitude earthquakes: A pilot study

This study presents a pilot investigation into a novel method for reconstructing real-time ground motion during small magnitude earthquakes (M < 4.5), removing the need for computationally expensive source characterization and simulation processes to assess ground shaking. Small magnitude earthquakes, which occur frequently and can be modeled as point sources, provide ideal conditions for evaluating real-time reconstruction methods. Utilizing sparse observation data, the method applies the Gappy Auto-Encoder (Gappy AE) algorithm for efficient field data reconstruction. This is the first study to apply the Gappy AE algorithm to earthquake ground motion reconstruction. Numerical experiments conducted with SW4 simulations demonstrate the method’s accuracy and speed across varying seismic scenarios. The reconstruction performance is further validated using real seismic data from the Berkeley area in California, USA, demonstrating the potential for practical application of real-time earthquake data reconstruction using Gappy AE. As a pilot investigation, it lays the groundwork for future applications to larger and more complex seismic events.

58 GEOSCIENCES↗
DOE OSTI2026

A Nonergodic Ground-Motion Model for the San Francisco Bay Area for Small-Magnitude Earthquakes

ABSTRACT Recently, generative models have become a computationally efficient alternative to physics-based numerical simulations of ground motions. Neural networks can learn from existing ground-motion data to generate unobserved ground-motion data at new source and site locations. A key challenge with generative models is ensuring that predicted ground motions remain within a physically realistic range. For this purpose, we developed an empirical, nonergodic ground-motion model (GMM) for small-magnitude earthquakes in the San Francisco Bay area based on about 5000 recordings per component for Mw ≤ 4 earthquakes. The nonergodic GMM predicts spatially varying median source, site, and path effects for both the Fourier amplitude spectrum (FAS) and the Fourier phase derivative (a proxy for duration), as well as the corresponding epistemic uncertainty for each term. For FAS, our model shows above-average source and site effects in the western part of the region and below-average effects in the eastern part, with regional effects exhibiting larger spatial correlation lengths with increasing frequency. For duration, the source term is negligible for small-magnitude earthquakes, and the site term leads to site-specific variations up to 5 s. Path effects for FAS and duration depend on the source–site pair and are extrapolated spatially using recent methods for path-effect modeling. The aleatory variability of the within-site within-path residuals is similar to the variability found in previous studies for other regions. The nonergodic model provides two key contributions: first, median adjustment terms that are transferable to larger magnitude earthquakes, further reducing aleatory variability in probabilistic seismic hazard analysis; second, region-specific criteria for validating machine learning-based ground-motion generators to evaluate whether synthetic ground motions exhibit physically realistic source, site, and path effects.

Lacour, Maxime↗
DOE OSTI2025

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

Automatic speech recognition predicts contemporaneous earthquake fault displacement

Abstract Significant progress has been made in probing the state of an earthquake fault by applying machine learning to continuous seismic waveforms. The breakthroughs were originally obtained from laboratory shear experiments and numerical simulations of fault shear, then successfully extended to slow-slipping faults. Here we apply the Wav2Vec-2.0 self-supervised framework for automatic speech recognition to continuous seismic signals emanating from a sequence of moderate magnitude earthquakes during the 2018 caldera collapse at the Kīlauea volcano on the island of Hawai’i. We pre-train the Wav2Vec-2.0 model using caldera seismic waveforms and augment the model architecture to predict contemporaneous surface displacement during the caldera collapse sequence, a proxy for fault displacement. We find the model displacement predictions to be excellent. The model is adapted for near-future prediction information and found hints of prediction capability, but the results are not robust. The results demonstrate that earthquake faults emit seismic signatures in a similar manner to laboratory and numerical simulation faults, and artificial intelligence models developed for encoding audio of speech may have important applications in studying active fault zones.

58 GEOSCIENCES↗