skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Improved Observations of Deep Earthquake Ruptures Using Machine Learning
Abstract Elevated seismic noise for moderate‐size earthquakes recorded at teleseismic distances has limited our ability to see their complexity. We develop a machine‐learning‐based algorithm to separate noise and earthquake signals that overlap in frequency. The multi‐task encoder‐decoder model is built around a kernel pre‐trained on local (e.g., short distances) earthquake data (Yin et al., 2022,https://doi.org/10.1093/gji/ggac290) and is modified by continued learning with high‐quality teleseismic data. We denoise teleseismic P waves of deep Mw5.0+ earthquakes and use the clean P waves to estimate source characteristics with reduced uncertainties of these understudied earthquakes. We find a scaling of moment and duration to beM0 ≃ τ4, and a resulting strong scaling of stress drop and radiated energy with magnitude ( and ). The median radiation efficiency is 5%, a low value compared to crustal earthquakes. Overall, we show that deep earthquakes have weak rupture directivity and few subevents, suggesting a simple model of a circular crack with radial rupture propagation is appropriate. When accounting for their respective scaling with earthquake size, we find no systematic depth variations of duration, stress drop, or radiated energy within the 100–700 km depth range. Our study supports the findings of Poli and Prieto (2016,https://doi.org/10.1002/2016jb013521) with a doubled amount of earthquakes investigated and with earthquakes of lower magnitudes.  more » « less
Award ID(s):
2124722
PAR ID:
10481204
Author(s) / Creator(s):
;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
128
Issue:
12
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract We analyze nearest‐neighbor proximities of earthquakes in California based on the joint distribution (T,R) of rescaled timeTand rescaled distanceRbetween pairs of earthquakes (Zaliapin & Ben‐Zion, 2013a,https://doi.org/10.1002/jgrb.50179), using seismic catalogs from several regions and several catalogs for the San Jacinto Fault Zone (SJFZ). The study aims to identify informative modes in nearest‐neighbor diagrams beyond the general background and clustered modes, and to assess seismic catalogs derived by different methods. The results show that earthquake clusters with large and small‐to‐medium mainshocks have approximately diagonal and horizontal (T,R) distributions of the clustered mode, respectively, reflecting different triggering distances of mainshocks. Earthquakes in the creeping section of San Andreas Fault have a distinct “repeaters mode” characterized by very large rescaled timesTand very small rescaled distancesR, due to nearly identical locations of repeating events. Induced seismicity in the Geysers and Coso geothermal fields follow mostly the background mode, but with larger rescaled timesTand smaller rescaled distancesRcompared to tectonic background seismicity. We also document differences in (T,R) distributions of catalogs constructed by different techniques (analyst‐picks, template‐matching and deep‐learning) for the SJFZ, and detect a mode with very largeRand smallTin the template‐matching and deep‐learning based catalogs. This mode may reflect dynamic triggering by passing waves and/or catalog artifacts. 
    more » « less
  2. ; ; ; ; ; ; ; (, Journal of Geophysical Research: Space Physics)
    Abstract We compared the performance of DREAM3D simulations in reproducing the long‐term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low‐energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics insideL* = 6, with a better model performance at lowerμthan higherμ, whereμis the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data‐driven outer boundary condition is critical to the model performance, while adding in the data‐driven seed population doesn't further improve the performance. (b) The model shows comparable performance withDLLfrom Brautigam and Albert (2000,https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014,https://doi.org/10.1002/2013ja019204), and Liu et al. (2016,https://doi.org/10.1002/2015gl067398), while withDLLfrom Ali et al. (2016,https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data‐based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higherμusing the Wang et al. (2019,https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under‐reproduced in the current model. 
    more » « less
  3. ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with globalM ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023),https://doi.org/10.1029/2023jb026487to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial‐temporal windows in Southern California tested in DeSalvio and Fan (2023),https://doi.org/10.1029/2023jb026487, thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023,https://doi.org/10.1029/2023jb026487), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor‐Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance. 
    more » « less
  4. ; ; ; ; ; ; ; ; (, Journal of Geophysical Research: Atmospheres)
    Abstract We report on the mountain top observation of three terrestrial gamma‐ray flashes (TGFs) that occurred during the summer storm season of 2021. To our knowledge, these are the first TGFs observed in a mountaintop environment and the first published European TGFs observed from the ground. A gamma‐ray sensitive detector was located at the base of the Säntis Tower in Switzerland and observed three unique TGF events with coincident radio sferic data characteristic of TGFs seen from space. We will show an example of a “slow pulse” radio signature (Cummer et al., 2011,https://doi.org/10.1029/2011GL048099; Lu et al., 2011,https://doi.org/10.1029/2010JA016141; Pu et al., 2019,https://doi.org/10.1029/2019GL082743; Pu et al., 2020,https://doi.org/10.1029/2020GL089427), a −EIP (Lyu et al., 2016,https://doi.org/10.1002/2016GL070154; Lyu et al., 2021,https://doi.org/10.1029/2021GL093627; Wada et al., 2020,https://doi.org/10.1029/2019JD031730), and a double peak TGF associated with an extraordinarily powerful and complicated positive‐polarity sferic, where each TGF peak is possibly preceded by a short burst of stepped leader emission. 
    more » « less
  5. ; ; (, Bulletin of the Seismological Society of America)
    ABSTRACT Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time- and frequency-domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016) source time function (STF) database and analyze how rupture complexity drives stress-drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress-drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email