Abstract Measures of foreshock occurrence are systematically examined using earthquake catalogs for eight regions (Italy, southern California, northern California, Costa Rica, Onshore Japan, Alaska, Turkey, and Greece) after imposing a magnitude ≥3.0 completeness level. Foreshocks are identified using three approaches: a magnitude-dependent space + fixed-time windowing method, a nearest-neighbor clustering method, and a modified magnitude-dependent space + variable-time windowing method. The method with fixed-time windows systematically yields higher counts of foreshocks than the other two clustering methods. We find similar counts of foreshocks across the three methods when the magnitude aperture is equalized by including only earthquakes in the magnitude range M*−2≤ M< M*, in which M* is the mainshock magnitude. For most of the catalogs (excluding Italy and southern California), the measured b-values of the foreshocks of all region-specific mainshocks are lower by 0.1–0.2 than b-values of respective aftershocks. Allowing for variable-time windows results in relatively high probabilities of having at least one foreshock in Italy (∼43%–56%), compared to other regional catalogs. Foreshock probabilities decrease to 14%–41% for regions such as Turkey, Greece, and Costa Rica. Similar trends are found when requiring at least five foreshocks in a sequence to be considered. Estimates of foreshock probabilities for each mainshock are method dependent; however, consistent regional trends exist regardless of method, with regions such as Italy and southern California producing more observable foreshocks than Turkey and Greece. Some regions with relatively high background seismicity have comparatively low probabilities of detectable foreshock activity when using methods that account for variable background, possibly due to depletion of near-failure fault conditions by background activity.
more »
« less
This content will become publicly available on June 28, 2024
Global Characteristics of Observable Foreshocks for Large Earthquakes
Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.
more » « less- Award ID(s):
- 1761987
- NSF-PAR ID:
- 10473633
- Publisher / Repository:
- SSA
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 94
- Issue:
- 5
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 2313 to 2325
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
SUMMARY Earthquakes come in clusters formed of mostly aftershock sequences, swarms and occasional foreshock sequences. This clustering is thought to result either from stress transfer among faults, a process referred to as cascading, or from transient loading by aseismic slip (pre-slip, afterslip or slow slip events). The ETAS statistical model is often used to quantify the fraction of clustering due to stress transfer and to assess the eventual need for aseismic slip to explain foreshocks or swarms. Another popular model of clustering relies on the earthquake nucleation model derived from experimental rate-and-state friction. According to this model, earthquakes cluster because they are time-advanced by the stress change imparted by the mainshock. This model ignores stress interactions among aftershocks and cannot explain foreshocks or swarms in the absence of transient loading. Here, we analyse foreshock, swarm and aftershock sequences resulting from cascades in a Discrete Fault Network model governed by rate-and-state friction. We show that the model produces realistic swarms, foreshocks and aftershocks. The Omori law, characterizing the temporal decay of aftershocks, emerges in all simulations independently of the assumed initial condition. In our simulations, the Omori law results from the earthquake nucleation process due to rate and state friction and from the heterogeneous stress changes due to the coseismic stress transfers. By contrast, the inverse Omori law, which characterizes the accelerating rate of foreshocks, emerges only in the simulations with a dense enough fault system. A high-density complex fault zone favours fault interactions and the emergence of an accelerating sequence of foreshocks. Seismicity catalogues generated with our discrete fault network model can generally be fitted with the ETAS model but with some material differences. In the discrete fault network simulations, fault interactions are weaker in aftershock sequences because they occur in a broader zone of lower fault density and because of the depletion of critically stressed faults. The productivity of the cascading process is, therefore, significantly higher in foreshocks than in aftershocks if fault zone complexity is high. This effect is not captured by the ETAS model of fault interactions. It follows that a foreshock acceleration stronger than expected from ETAS statistics does not necessarily require aseismic slip preceding the mainshock (pre-slip). It can be a manifestation of a cascading process enhanced by the topological properties of the fault network. Similarly, earthquake swarms might not always imply transient loading by aseismic slip, as they can emerge from stress interactions.
-
more » « less
Abstract The number of aftershocks increases with mainshock size following a well‐defined scaling law. However, excursions from the average behavior are common. This variability is particularly concerning for large earthquakes where the number of aftershocks varies by factors of 100 for mainshocks of comparable magnitude. Do observable factors lead to differences in aftershock behavior? We examine aftershock productivity relative to the global average for all mainshocks (
) from 1990 to 2019. A global map of earthquake productivity highlights the influence of tectonic regimes. Earthquake depth, lithosphere age, and plate boundary type correspond well with earthquake productivity. We investigate the role of mainshock attributes by compiling source dimensions, radiated seismic energy, stress drop, and a measure of slip heterogeneity based on finite‐fault source inversions for the largest earthquakes from 1990 to 2017. On an individual basis, stress drop, normalized rupture width, and aspect ratio most strongly correlate with aftershock productivity. A multivariate analysis shows that a particular set of parameters (dip, lithospheric age, and normalized rupture area) combines well to improve predictions of aftershock productivity on a cross‐validated data set. Our overall analysis is consistent with a model in which the volumetric abundance of nearby stressed faults controls the aftershock productivity rather than variations in source stress. Thus, we suggest a complementary approach to aftershock forecasts based on geological and rupture properties rather than local calibration alone. -
Abstract Recognizing earthquakes as foreshocks in real time would provide a valuable forecasting capability. In a recent study, Gulia and Wiemer (2019) proposed a traffic-light system that relies on abrupt changes in b-values relative to background values. The approach utilizes high-resolution earthquake catalogs to monitor localized regions around the largest events and distinguish foreshock sequences (reduced b-values) from aftershock sequences (increased b-values). The recent well-recorded earthquake foreshock sequences in Ridgecrest, California, and Maria Antonia, Puerto Rico, provide an opportunity to test the procedure. For Ridgecrest, our b-value time series indicates an elevated risk of a larger impending earthquake during the Mw 6.4 foreshock sequence and provides an ambiguous identification of the onset of the Mw 7.1 aftershock sequence. However, the exact result depends strongly on expert judgment. Monte Carlo sampling across a range of reasonable decisions most often results in ambiguous warning levels. In the case of the Puerto Rico sequence, we record significant drops in b-value prior to and following the largest event (Mw 6.4) in the sequence. The b-value has still not returned to background levels (12 February 2020). The Ridgecrest sequence roughly conforms to expectations; the Puerto Rico sequence will only do so if a larger event occurs in the future with an ensuing b-value increase. Any real-time implementation of this approach will require dense instrumentation, consistent (versioned) low completeness catalogs, well-calibrated maps of regionalized background b-values, systematic real-time catalog production, and robust decision making about the event source volumes to analyze.more » « less
-
more » « less
Abstract We compare source parameter estimates for earthquakes in the 2011 Prague Mw 5.7, Oklahoma, sequence to investigate random uncertainty and systematic bias, and resolve reliable relative variations in stress drop. Source parameters provide insight into the earthquake rupture processes but large variations between studies occur. The Prague earthquake sequence is a prime example of this, with different studies reaching contrasting interpretations of the effects of injection on source parameters. We examine the Prague earthquake sequence using a single coherent catalog for all the events detected by the Oklahoma Geological Survey (OGS) and McMahon et al. (2017). We use three principal approaches to estimate stress drop in order to understand the biases of each: a spectral decomposition method based on stacking, individual event spectral modeling, and a spectral ratio method based on highly correlated events. We also compare our results with previous studies for the Prague sequences aftershocks, as well as past results for the Mw 4.8 foreshock and Mw 4.8 aftershock and Mw 5.7 mainshock. The absolute values of stress drop vary significantly between methods, but the relative patterns remain consistent, except when low quality or low bandwidth data are included. The consistent relative patterns reveal that the stress drops of aftershocks are dependent on the fault orientation and the proximity of the events to the mainshocks slip. These results indicate that fault structure as well as past events play an important role in stress drop patterns.
