Published June 2025 | Version Supplemental material
Journal Article Open

On Seismological Stress-Drop Estimates for Earthquake Sources Simulated on Rate-and-State Faults

  • 1. ROR icon National Central University
  • 2. ROR icon California State University, Northridge
  • 3. ROR icon California Institute of Technology

Abstract

Seismologists estimate stress drops of small earthquakes based on specific theoretical source models. We explore the accuracy of the stress-drop estimates for several earthquake source models obtained in dynamic simulations on rate-and-state faults. We consider Madariaga-like symmetric circular sources as well as sources with directivity, elongated shapes, partial ruptures, and complex changes in the slip direction. The energy-based average stress drops computed directly on the fault for all simulated source models range from 1.5 to 5 MPa. We consider a range of focal depths and fault dips that results in 980 scenarios overall with respect to a surface network of 16 stations, where we produce synthetic waveforms assuming a known homogeneous velocity structure, and use them to obtain seismologically inferred stress drops. For the second-moment approach and spectral-fitting approach based on S waves and n = 2, the stress drops for most sources are reproduced well on average but with a significant scatter from nearly 0.01 to 100 MPa, representative of scatter for natural earthquakes, despite the actual stress-drop variation of 1.5–5 MPa. The scatter is smaller by a factor of 2 for the second-moment approach. The spectral-fitting approach based on P waves consistently underestimates the stress drops for noncircular sources. All approaches underestimate stress drops for ring-like sources, which leave part of the seismogenic patch unruptured. The spectral estimates are significantly affected by different averages of corner frequencies over the focal sphere for our sources versus typically assumed simple theoretical sources, as was already pointed out for some of the sources by Kaneko and Shearer (2015) and Lin and Lapusta (2018). For both second-moment and spectral methods, the scatter is amplified by partial coverage of the focal sphere by the assumed station geometry, which can also cause systematic depth-dependent artifacts.

Copyright and License

© 2025 Seismological Society of America.

Acknowledgement

This study was supported by the National Science Foundation (NSF; Grant Numbers 1724686 and 2139331). It was also supported by the Taiwan Earthquake Research Center (TEC), funded through the National Science and Technology Council (NSTC) with Grant Numbers 105‐2917‐I‐564‐015‐, 108‐2116‐M‐008‐025‐MY2, and 113‐2116‐M‐008‐007‐. The TEC number for this article is 00196. This work is financially supported by the “Earthquake‐Disaster & Risk Evaluation and Management Center, E‐DREaM” from The Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education (MOE) in Taiwan. The data used are listed in the references, tables, and supplemental material. The authors are grateful for the thorough reviews and suggestions from Rachel Abercrombie, Bill Ellsworth, and an anonymous reviewer, which served to significantly strengthen the article.

Data Availability

All data used in this article are either from our simulations or from published sources listed in the references. Some calculations and plots were made using MATLAB (https://www.mathworks.com/products/matlab.html). Some plots were made using the Generic Mapping Tools (GMT) version 4.5 (www.soest.hawaii.edu/gmtWessel and Smith, 1998). Some plots were modified by Illustrator CC which is available at https://www.adobe.com/tw/products/illustrator.html. All websites were last accessed in December 2024. In supplemental material, we provide all simulated source models that are analyzed in this article (and listed in Table 3) in terms of their slip‐rate variation in space and time as well as our MATLAB code for reading the slip‐rate functions (STF.zip). Table S1 demonstrates the estimates of source parameters for all models.

Supplemental Material

Files

Files (8.1 MB)

Name Size Download all
md5:6974e95184614abd9ffb4e3c7d6c039d
7.9 MB Download
md5:751490c14307b3777ad60b334abcd0ef
234.9 kB Download

Additional details

Dates

Available
2025-03-04
First online

Caltech Custom Metadata

Caltech groups
Seismological Laboratory, Division of Engineering and Applied Science (EAS), Division of Geological and Planetary Sciences (GPS)
Publication Status
Published