An official website of the United States government
Both large and small earthquakes rupture in complex ways. However, microearthquakes are often simplified as point sources and their rupture properties are challenging to resolve. We leverage seismic wavefields recorded by a dense array in Oklahoma to image microearthquake rupture processes. We construct machine-learning enabled catalogs and identify four spatially disconnected seismic clusters. These clusters likely delineate near-vertical strike-slip faults. We develop a new approach to use the maximum absolute SH-wave amplitude distributions (S-wave wavefields) to compare microearthquake rupture processes. We focus on one cluster with earthquakes located beneath the dense array and have a local magnitude range of -1.3 to 2.3. The S-wave wavefields of single earthquakes are generally coherent but differ slightly between the low-frequency (<12 Hz) and high-frequency (>12 Hz) bands. The S-wave wavefields are coherent between different earthquakes at low frequencies with average correlation coefficients greater than 0.95. However, the wavefield coherence decreases with increasing frequency for different earthquakes. This reduced coherence is likely due to the rupture differences among individual earthquakes. Our results suggest that earthquake slip of the microearthquakes dominates the radiated S-wave wavefields at higher frequencies. Our method suggests a new direction in resolving small earthquake source attributes using dense seismic arrays without assuming a rupture model.
more »
« less
What makes low-frequency earthquakes low frequency
Low-frequency earthquakes, atypical seismic events distinct from regular earthquakes, occur downdip of the seismogenic megathrust where an aseismic rheology dominates the subduction plate boundary. Well situated to provide clues on the slip regime of this unique faulting environment, their distinctive waveforms reflect either an unusual rupture process or unusually strong attenuation in their source zone. We take advantage of the unique geometry of seismicity in the Nankai Trough to isolate the spectral signature of low-frequency earthquakes after correcting for empirically derived attenuation. We observe that low-frequency earthquake spectra are consistent with the classical earthquake model, yet their rupture duration and stress drop are orders of magnitude different from ordinary earthquakes. We conclude their low-frequency nature primarily results from an atypical seismic rupture process rather than near-source attenuation.
more »
« less
- PAR ID:
- 10502083
- Publisher / Repository:
- American Association for the Advancement of Science
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 9
- Issue:
- 32
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Low‐frequency earthquakes are a seismic manifestation of slow fault slip. Their emergent onsets, low amplitudes, and unique frequency characteristics make these events difficult to detect in continuous seismic data. Here, we train a convolutional neural network to detect low‐frequency earthquakes near Parkfield, CA using the catalog of Shelly (2017),https://doi.org/10.1002/2017jb014047as training data. We explore how varying model size and targets influence the performance of the resulting network. Our preferred network has a peak accuracy of 85% and can reliably pick low‐frequency earthquake (LFE) S‐wave arrival times on single station records. We demonstrate the abilities of the network using data from permanent and temporary stations near Parkfield, and show that it detects new LFEs that are not part of the Shelly (2017),https://doi.org/10.1002/2017jb014047catalog. Overall, machine‐learning approaches show great promise for identifying additional low‐frequency earthquake sources. The technique is fast, generalizable, and does not require sources to repeat.more » « less
-
ABSTRACT Although the Brune source model describes earthquake moment release as a single pulse, it is widely used in studies of complex earthquakes with multiple episodes of high moment release (i.e., multiple subevents). In this study, we investigate how corner frequency estimates of earthquakes with multiple subevents are biased if they are based on the Brune source model. By assuming complex sources as a sum of multiple Brune sources, we analyze 1640 source time functions of Mw 5.5–8.0 earthquakes in the seismic source characteristic retrieved from deconvolving teleseismic body waves catalog to estimate the corner frequencies, onset times, and seismic moments of subevents. We identify more subevents for strike-slip earthquakes than dip-slip earthquakes, and the number of resolvable subevents increases with magnitude. We find that earthquake corner frequency correlates best with the corner frequency of the subevent with the highest moment release (i.e., the largest subsevent). This suggests that, when the Brune model is used, the estimated corner frequency and, therefore, the stress drop of a complex earthquake is determined primarily by the largest subevent rather than the total rupture area. Our results imply that, in addition to the simplified assumption of a radial rupture area with a constant rupture velocity, the stress variation of asperities, rather than the average stress change of the whole fault, contributes to the large variance of stress-drop estimates.more » « less
-
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
-
Abstract Megathrust earthquakes exhibit a ubiquitous seismic radiation style: low‐frequency (LF) seismic energy is efficiently emitted from the shallowest portion of the fault, whereas high‐frequency (HF) seismic energy is efficiently emitted from the deepest part of the fault. Although this is observed in many case‐specific studies, we show that it is ubiquitous in global megathrust earthquakes between 1995 and 2021. Previous studies have interpreted this as an effect of systematic depth variation in either the plate interface frictional properties (Lay et al., 2012) or the P wavespeeds (Sallarès & Ranero, 2019). This work suggests an alternative hypothesis: the interaction between waves and ruptures due to the Earth's free surface is the leading mechanism that generates this behavior. Two‐dimensional dynamic rupture simulations of subduction zone earthquakes support this hypothesis. Our simulations show that the interaction between the seismic waves reflected at the Earth's free surface and the updip propagating rupture results in LF radiation at the source. In contrast, the downdip propagation of rupture is less affected by the free surface and is thus dominated by HF radiation typical of buried faults. To a second degree, the presence of a realistic Earth structure derived from P‐wave velocity (VP) tomographic images and realistic VP/VSratio estimated in boreholes further enhances the contrast in source radiation. We conclude that the Earth's free surface is necessary to explain the observed megathrust earthquake radiation style, and the realistic structure of subduction zone is necessary to better predict earthquake ground motion and tsunami potential.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1126/sciadv.adh3688
