An official website of the United States government
Abstract Dynamic rupture simulations generate synthetic waveforms that account for nonlinear source and path complexity. Here, we analyze millions of spatially dense waveforms from 3D dynamic rupture simulations in a novel way to illuminate the spectral fingerprints of earthquake physics. We define a Brune-type equivalent near-field corner frequency (fc) to analyze the spatial variability of ground-motion spectra and unravel their link to source complexity. We first investigate a simple 3D strike-slip setup, including an asperity and a barrier, and illustrate basic relations between source properties and fc variations. Next, we analyze >13,000,000 synthetic near-field strong-motion waveforms generated in three high-resolution dynamic rupture simulations of real earthquakes, the 2019 Mw 7.1 Ridgecrest mainshock, the Mw 6.4 Searles Valley foreshock, and the 1992 Mw 7.3 Landers earthquake. All scenarios consider 3D fault geometries, topography, off-fault plasticity, viscoelastic attenuation, and 3D velocity structure and resolve frequencies up to 1–2 Hz. Our analysis reveals pronounced and localized patterns of elevated fc, specifically in the vertical components. We validate such fc variability with observed near-fault spectra. Using isochrone analysis, we identify the complex dynamic mechanisms that explain rays of elevated fc and cause unexpectedly impulsive, localized, vertical ground motions. Although the high vertical frequencies are also associated with path effects, rupture directivity, and coalescence of multiple rupture fronts, we show that they are dominantly caused by rake-rotated surface-breaking rupture fronts that decelerate due to fault heterogeneities or geometric complexity. Our findings highlight the potential of spatially dense ground-motion observations to further our understanding of earthquake physics directly from near-field data. Observed near-field fc variability may inform on directivity, surface rupture, and slip segmentation. Physics-based models can identify “what to look for,” for example, in the potentially vast amount of near-field large array or distributed acoustic sensing data.
more »
« less
Resolving Differences in the Rupture Properties of M5 Earthquakes in California Using Bayesian Source Spectral Analysis
Abstract 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 observedS‐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
- Award ID(s):
- 2121666
- 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
-
-
ABSTRACT Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground-motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well-recorded small-to-mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high-quality data from stations in the direction of rupture, leading to average stress-drop increases of 50% if directivity is not accounted for.more » « less
-
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
-
One of most universal statistical properties of earthquakes is the tendency to cluster in space and time. Yet while clustering is pervasive, individual earthquake sequences can vary markedly in duration, spatial extent, and time evolution. In July 2014, a prolific earthquake sequence initiated within the Sheldon Wildlife Refuge in northwest Nevada, USA. The sequence produced 26 M4 earthquakes and several hundred M3s, with no clear mainshock or obvious driving force. Here we combine a suite of seismological analysis techniques to better characterize this unusual earthquake sequence. High-precision relocations reveal a clear, east dipping normal fault as the dominant structure that intersects with a secondary, subvertical cross fault. Seismicity occurs in burst of activity along these two structures before eventually transitioning to shallower structures to the east. Inversion of hundreds of moment tensors constrain the overall normal faulting stress regime. Source spectral analysis suggests that the stress drops and rupture properties of these events are typical for tectonic earthquakes in the western US. While regional station coverage is sparse in this remote study region, the timely installation of a temporary seismometer allows us to detect nearly 70,000 earthquakes over a 40-month time period when the seismic activity is highest. Such immense productivity is difficult to reconcile with current understanding of crustal deformation in the region and may be facilitated by local hydrothermal processes and earthquake triggering at the transitional intersection of subparallel fault systems.more » « less
-
The 2019 Ridgecrest, CA earthquake sequence has provided a unique opportunity and a rich dataset to understand earthquake source properties and near-fault structure. Using the high-quality seismic data provided by the SCEC Stress Drop Validation group, we first estimate the corner frequency of M2.0-4.5 earthquakes by applying the spectral ratio method based on empirical Green’s function (Liu et al., 2020). We relate corner frequency estimates to stress drops assuming the Brune source model and circular cracks. Our preliminary results show increasing median stress drops with magnitude for both P and S waves, from 1 MPa for M2.0 events to 10 MPa for M4.0 events, though the limited frequency bandwidth may cause underestimation for small events. The estimated moment magnitude is proportional to the catalog magnitude by a factor of 0.72, which is close to 0.74 estimated by Trugman (2020) for the Ridgecrest earthquake sequence. In the second part of the study, we examine the impact of fault zone structure on the azimuthal variation of the source spectra. Using kinematic simulations and observations of the 2003 Big Bear earthquake sequence, Huang et al. (2016) showed that fault damage zones can act as an effective wave guide and cause high-frequency wave amplification along directions close to fault strike. We use clusters of M1.5-3 earthquakes in the Ridgecrest region to further examine the azimuthal variation of the stacked source spectra and investigate if the near-source structure can affect our corner frequency estimates. We aim to develop robust methods that utilize high-quality seismic data to illuminate earthquake source processes and fault zone properties.more » « less
