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.


Title: Stress Drop Variation of Deep‐Focus Earthquakes Based on Empirical Green's Functions
Abstract We analyze source characteristics of global, deep‐focus (>350 km) earthquakes with moment magnitudes (Mw) larger than 6.0–8.2 using teleseismic P‐wave and S‐wave spectra and an empirical Green's functions approach. We estimate the corner frequency assuming Brune's source model and calculate stress drops assuming a circular crack model. Based on P‐wave and S‐wave spectra, the one standard deviation ranges are 3.5–369.8 and 8.2–328.9 MPa, respectively. Based on the P‐wave analysis, the median of our stress drop estimates is about a factor of 10 higher than the median stress drop of shallow earthquakes with the same magnitude estimated by Allmann and Shearer (2009,https://doi.org/10.1029/2008JB005821). This suggests that, on average, the shear stress of deep faults in the mantle transition zone is an order of magnitude higher than the shear stress of faults in the crust. The wide range of stress drops implies coexistence of multiple physical mechanisms.  more » « less
Award ID(s):
2019379
PAR ID:
10376561
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
47
Issue:
9
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract Multiple physical mechanisms have been proposed to explain the cause of intermediate‐depth and deep earthquakes, but they are still under debate. Source parameters such as stress drop, have the potential to provide insight into their physical mechanisms. We develop a modified spectral decomposition method to analyze 1‐year seismic data from temporary land‐based and ocean bottom seismographs in a complex subduction zone. By applying this method to investigate 1,083 intermediate‐depth and deep earthquakes in the Tonga slab, we successfully resolve the source spectra and stress drops of 743MW2.6–6.0 earthquakes. Although the absolute stress drops are subject to the choices of source model parameters, the relative stress drops are more reliably resolved. The median stress drop of Tonga earthquakes does not change with respect to magnitude but decreases with depth by 2–3 times in two separate depth ranges of 70–250 and 400–600 km, corresponding to intermediate‐depth and deep earthquakes, respectively. The median stress drops show spatial variations, with two high‐stress‐drop (five times higher than the ambient value) regions, coinciding with strong local deformation where the Tonga slab bends or tears. In the Tonga double seismic zone at 120–300 km depths, the median stress drop appears smaller in the lower plane than in the upper plane, suggesting a slower rupture velocity or a higher fluid content in the lower‐plane region. Our results suggest that intermediate‐depth and deep earthquakes in the Tonga slab generally follow the earthquake self‐similar model and favor the fluid‐related embrittlement hypothesis for both groups of earthquakes. 
    more » « less
  3. Abstract Liu et al. (2022,https://doi.org/10.1029/2021GL093691) used Rayleigh waves extracted from the cross‐correlation of ambient noise recorded by two stations to monitor the seismic velocity variations associated with the 2018 Kı̄lauea eruption. However, their study ignored the fact that the tremors on the Island of Hawai'i were dominated by a source at the Kı̄lauea summit before the eruption. Close inspection of the waveforms of the station pair PAUD‐STCD shows a simple, mistakenly identified wave traveling direction in Liu et al. (2022,https://doi.org/10.1029/2021GL093691). A correct wave traveling direction agrees with the noise source model, where the dominant tremor source should be at the Kı̄lauea summit. Because of the drastic change in the tremor source after the eruption, the cross‐correlation of the tremor records may reflect predominantly changes in the source rather than in the medium properties between the two stations. 
    more » « less
  4. Abstract This letter compares the predictions of two expressions proposed for the porosity evolution in the context of rate and state friction. One (Segall & Rice, 1995,https://doi.org/10.1029/95jb02403) depends only on the sliding velocity; the other (Sleep, 1995,https://doi.org/10.1029/94jb03340) depends only on the state variable. Simulations of both are similar for velocity stepping and slide‐hold‐slide experiments. They differ significantly for normal effective stress jumps at constant sliding velocity. Segall and Rice (1995,https://doi.org/10.1029/95jb02403) predicts no change in the porosity; Sleep (1995,https://doi.org/10.1029/94jb03340) does. Simulation with a spring‐block model indicates that the magnitude of rapid slip events is essentially the same for the two formulations. Variations of porosity and induced pore pressure near rapid slip events are similar and consistent with experimental observations. Predicted porosity variations during slow slip intervals and the time at which rapid slip events occur are significantly different. The simulation indicates that changes in friction stress due to pore pressure changes exceed those due to rate and state effects. 
    more » « less
  5. 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