The spectra of earthquake waveforms can provide important insight into rupture processes, but the analysis and interpretation of these spectra is rarely straightforward. Here we develop a Bayesian framework that embraces the inherent data and modeling uncertainties of spectral analysis to infer key source properties. The method uses a spectral ratio approach to correct the observed
We calculate rupture directivity and velocity for earthquakes in three well‐recorded repeating sequences (2001–2016) on the San Andreas Fault at Parkfield using
- NSF-PAR ID:
- 10451675
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 47
- Issue:
- 23
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract S ‐wave spectra of nearby earthquakes for path and site attenuation. The objective then is to solve for a joint posterior probability distribution of three source parameters—seismic moment, corner frequency, and high‐frequency falloff rate—for each earthquake in the sequence, as well as a measure of rupture directivity for select target events with good azimuthal station coverage. While computationally intensive, this technique provides a quantitative understanding of parameter tradeoffs and uncertainties and allows one to impose physical constraints through prior distributions on all source parameters, which guide the inversion when data is limited. We demonstrate the method by analyzing in detail the source properties of 14 different target events of magnitude M5 in southern California that span a wide range of tectonic regimes and fault systems. These prominent earthquakes, while comparable in size, exhibit marked diversity in their source properties and directivity, with clear spatial patterns, depth‐dependent trends, and a preference for unilateral directivity. These coherent spatial variations source properties suggest that regional differences in tectonic setting, hypocentral depth or fault zone characteristics may drive variability in rupture processes, with important implications for our understanding of earthquake physics and its relation to hazard. -
more » « less
Abstract We demonstrate the complex source process of a potentially induced Mw 4.0 earthquake near Guthrie, Oklahoma. Relative source time functions (RSTFs), retrieved using a time domain empirical Green's function (EGF) deconvolution method, show clear evidence of rupture complexity with the presence of multiple pulses and systematic azimuthal variations. Directivity analysis reveals a well‐defined rupture propagation at ~120° toward the southeast. Detailed modeling of the RSTFs results in a four‐subevent model spanning ~1 s in time and ~1.5 km in space, in which the earthquake ruptured four successive slip episodes unilaterally to the southeast. Multiple EGF spectral ratio analysis also shows significant complexity as evidenced by azimuth‐dependent deviations from the omega‐square (Brune) model. The observations of the concentration of early aftershocks at the periphery of the inferred rupture model along the rupture direction and amplified ground motions at stations to the southeast can both be well explained by the directivity effect.
-
more » « less
Abstract We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of
P andS phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage. -
more » « less
Abstract In September and October 2015, three
M 4+ earthquakes occurred as a sequence along a fault northwest of the Cushing city, Oklahoma, followed by anotherM 5 earthquake in November 2016. While previous studies have shown that moderate‐size earthquakes in Oklahoma are likely induced by wastewater injections, it is still not clear what controls the rupture process and spatiotemporal evolutions of seismicity during individual sequences. In this study, we investigated the rupture process of these fourM 4‐5 events in 2015–2016 with finite fault model (FFM) inversions, and computed the static stress changes during this sequence. We found that the rupture processes of fourM 4‐5 earthquakes were very complex, and each of them had several subevents with different rupture directivities. The 2016M 5 earthquake started near the region where threeM 4+ events initiated, but the majority of the slip occurred a few kilometers away in the northeast direction. In comparison, the 2015M 4.3 event mainly ruptured toward the southwest direction. Due to data limitation and inversion uncertainties, we were unable to constrain the rupture directivities for the other twoM 4+ events. The foreshocks 3 days before the firstM 4+ earthquake in 2015 occurred in a region of positive shear stress changes caused by previous earthquakes in 2014–2015 on unmapped faults several kilometers to the south. Our results suggest small‐scale heterogeneity in controlling complex seismicity and rupture patterns in the 2015–2016 Cushing sequence, and possible triggering of this sequence by a small stress perturbation on order of a few kilopascals. -
more » « less
Abstract We conduct a detailed study of the foreshock sequence preceding the 2010
M w6.7 Yushu, Qinghai earthquake in the Tibetan plateau by examining continuous waveforms recorded at a seismic station near the mainshock rupture zone. By using a deep learning phase picker—EQTransformer and a matched‐filter technique, we identify 120 foreshocks with magnitude ranging from −0.7 to 1.6, starting with aM w4.6 foreshock approximately 2 hr before theM w6.7 Yushu mainshock. Our analyses show that the foreshock sequence follows a typical Omori's law decay with ap ‐value of 0.73 and the Gutenberg‐Richer frequency‐magnitudeb ‐value of 0.66. We do not find any evidence of accelerating events leading up to the Yushu mainshock. Hence, they could be considered as aftershocks of theM w4.6 earthquake. We further invert for the focal mechanisms and rupture directions for both the largest foreshock and the mainshock. TheM w4.6 foreshock likely occurred on a NE‐SW trending fault conjugating to the NW‐SE trending fault of the mainshock. Coulomb stress analysis suggests theM w4.6 foreshock induces negative stress on the mainshock source area. These observations do not support either the pre‐slip or the cascade triggering model for foreshock generation. The occurrence of the foreshock, mainshock and large aftershocks appear to be modulated by the Earth's tidal forces, likely reflecting the role of high pore‐fluid pressures. Our observations, together with other recent studies, suggest that extensional step‐overs and conjugate faults along major strike‐slip faults play an important role in generating short‐term foreshock sequences.
