We calculate rupture directivity and velocity for earthquakes in three well‐recorded repeating sequences (2001–2016) on the San Andreas Fault at Parkfield using
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
- Award ID(s):
- 2121666
- NSF-PAR ID:
- 10369512
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 4
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract P waves from borehole recordings and the empirical Green's function method. The individual events in each sequence all show the same directivity; the largest magnitude sequence (M ~ 2.7, 8 events) ruptures unilaterally NW (at ~0.8Vs ), the second sequence (M ~ 2.3, 9 events) ruptures unilaterally SE, and the smallest magnitude sequence (M ~ 2, 11 events) is less well resolved. The highly repetitive rupture suggests that geometry or material properties might control nucleation of small locked patches. The source spectra of theM ~ 2.7 sequence exhibit no detectable temporal variation. The smallerM sequences both exhibit a decrease in high‐frequency energy following theM 6 earthquake that recovers with time. This could indicate a decrease in stress drop, an increase in attenuation, or a combination of the two, followed by gradual healing. -
Fault damage zones can influence various aspects of the earthquake cycle, such as the recurrence intervals and magnitudes of large earthquakes. The properties and structure of fault damage zones are often characterized using dense arrays of seismic stations located directly above the faults. However, such arrays may not always be available. Hence, our research aims to develop a novel method to image fault damage zones using broadband stations at relatively larger distances. Previous kinematic simulations and a case study of the 2003 Big Bear earthquake sequence demonstrated that fault damage zones can act as effective waveguides, amplifying high-frequency waves along directions close to fault strike via multiple reflections within the fault damage zone. The amplified high-frequency energy can be observed by stacking P-wave spectra of earthquake clusters with highly-similar waveforms (Huang et al., 2016), and the frequency band which is amplified may be used to estimate the width and velocity contrast of the fault damage zone. We attempt to identify the high-frequency peak associated with fault zone waves in stacked spectra by conducting a large-scale study of small earthquakes (M1.5–3). We use high quality broadband data from seismic stations at hypocentral distances of 20-80 km in the 2019 Ridgecrest earthquake regions. First, we group the Ridgecrest earthquakes in clusters by their locations and their waveform similarity, and then stack their velocity spectra to average the source effects of individual earthquakes. Our results show that the stations close to the fault strike record more high-frequency energies around the characteristic frequency of fault zone reflections. We find that the increase in the amount of high-frequencies is consistent across clusters with average magnitudes ranging from 1.6-2.4, which suggests that the azimuthal variation in spectra is caused by fault zone amplification rather than rupture directivity. We will apply our method to other fault zones in California, in order to search for fault damage zone structures and estimate their material properties.more » « less
-
Abstract Analysis of earthquake spectra can aid in understanding source characteristics like stress drop and rupture complexity. There is growing interest in probing the similarities and differences of fault rupture for natural and human-induced seismic events. Here, we analyze waveform data from a shallow, buried geophone array that recorded seismicity during a hydraulic fracturing operation near Fox Creek, Alberta. Starting from a quality-controlled catalog of 4000 events between magnitude 0 and 3.2, we estimate source-spectral corner frequencies using methods that account for the band-limited nature of the sensor response. The stress-drop values are found to be approximately self-similar, but with a slight magnitude dependence in which larger events have higher stress drop (∼10 MPa). Careful analysis of the relative corner frequencies shows that individual fault and fracture segments experienced systematic variations in relative corner frequency over time, indicating a possible change in the stress state. Clustering analysis of source spectra based on the relative proportion of high- and low-frequency content relative to the Brune model further shows that event complexity evolves over time. In addition, the faults produce earthquakes with systematically larger stress-drop values than the fractures. Combined, these results indicate that the features activated by hydraulic fracturing experience observable changes in source behavior over time and exhibit different properties depending on the orientation, scale, and fabric of the structural feature on which they occur.more » « less
-
more » « less
SUMMARY Earthquake ruptures are complex physical processes that may vary with the structure and tectonics of the region in which they occur. Characterizing the factors controlling this variability would provide fundamental constraints on the physics of earthquakes and faults. We investigate this by determining finite source properties from second moments of the stress glut for a global data set of large strike-slip earthquakes. Our approach uses a Bayesian inverse formulation with teleseismic body and surface waves, which yields a low-dimensional probabilistic description of rupture properties including the spatial deviation, directivity and temporal deviation of the source. This technique is useful for comparing events because it makes only minor geometric constraints, avoids bias due to rupture velocity parametrization and yields a full ensemble of possible solutions given the uncertainties of the data. We apply this framework to all great strike-slip earthquakes of the past three decades, and we use the resultant second moments to compare source quantities like directivity ratio, rectilinearity, average moment density and vertical deviation. We find that most strike-slip earthquakes have a large component of unilateral directivity, and many of these earthquakes show a mixture of unilateral and bilateral behaviour. We notice that oceanic intraplate earthquakes usually rupture a much larger width of the seismogenic zone than other strike-slip earthquakes, suggesting these earthquakes may often breach the expected thermal boundary for oceanic ruptures. We also use these second moments to resolve nodal plane ambiguity for the large oceanic intraplate earthquakes and find that the rupture orientation is usually unaligned with encompassing fossil fracture zones.
-
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.
