Foreshocks can provide valuable information about possible nucleation process of a mainshock. However, their physical mechanisms are still under debate. In this study, we present a comprehensive analysis of the earthquake sequence preceding the 2010 Mw7.2 El Mayor‐Cucapah mainshock, including waveform detection of missing smaller events, relative relocation, and source parameter analysis. Based on a template matching method, we find a tenfold increase in the number of earthquakes than reported in the Southern California Seismic Network catalog. The entire sequence exhibits nearly continuous episodes of foreshocks that can be loosely separated into two active clusters. Relocated foreshocks show several seismicity streaks at depth, with a consistently active cluster at depths between 14 and 16 km where the mainshock was nucleated. Stress drop measurements from a spectral ratio approach based on empirical Green's functions show a range between 3.8 and 41.7 MPa with a median of 13.0 MPa and no clear temporal variations. The relocation results, together with the source patches estimated from earthquake corner frequencies, revealed a migration front toward the mainshock hypocenter within last 8 hr and a chain of active burst immediately 6 min prior to the mainshock. Our results support combined effects of aseismic slip and cascading failure on the evolution of foreshocks.
- NSF-PAR ID:
- 10447103
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 94
- Issue:
- 1
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 507 to 521
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Volcanic earthquake catalogs are an essential data product used to interpret subsurface volcanic activity and forecast eruptions. Advances in detection techniques (e.g., matched-filtering, machine learning) and relative relocation tools have improved catalog completeness and refined event locations. However, most volcano observatories have yet to incorporate these techniques into their catalog-building workflows. This is due in part to complexities in operationalizing, automating, and calibrating these techniques in a satisfactory way for disparate volcano networks and their varied seismicity. In an effort to streamline the integration of catalog-enhancing tools at the Alaska Volcano Observatory (AVO), we have integrated four popular open-source tools: REDPy, EQcorrscan, HypoDD, and GrowClust. The combination of these tools offers the capability of adding seismic event detections and relocating events in a single workflow. The workflow relies on a combination of standard triggering and cross-correlation clustering (REDPy) to consolidate representative templates used in matched-filtering (EQcorrscan). The templates and their detections are then relocated using the differential time methods provided by HypoDD and/or GrowClust. Our workflow also provides codes to incorporate campaign data at appropriate junctures, and calculate magnitude and frequency index for valid events. We apply this workflow to three datasets: the 2012–2013 seismic swarm sequence at Mammoth Mountain (California), the 2009 eruption of Redoubt Volcano (Alaska), and the 2006 eruption of Augustine Volcano (Alaska); and compare our results with previous studies at each volcano. In general, our workflow provides a significant increase in the number of events and improved locations, and we relate the event clusters and temporal progressions to relevant volcanic activity. We also discuss workflow implementation best practices, particularly in applying these tools to sparse volcano seismic networks. We envision that our workflow and the datasets presented here will be useful for detailed volcano analyses in monitoring and research efforts.
-
null (Ed.)SUMMARY We examine localization processes of low magnitude seismicity in relation to the occurrence of large earthquakes using three complementary analyses: (i) estimated production of rock damage by background events, (ii) evolving occupied fractional area of background seismicity and (iii) progressive coalescence of individual earthquakes into clusters. The different techniques provide information on different time scales and on the spatial extent of weakened damaged regions. Techniques (i) and (ii) use declustered catalogues to avoid the occasional strong fluctuations associated with aftershock sequences, while technique (iii) examines developing clusters in entire catalogue data. We analyse primarily earthquakes around large faults that are locked in the interseismic periods, and examine also as a contrasting example seismicity from the creeping Parkfield section of the San Andreas fault. Results of analysis (i) show that the M > 7 Landers 1992, Hector Mine 1999, El Mayor-Cucapah 2010 and Ridgecrest 2019 main shocks in Southern and Baja California were preceded in the previous decades by generation of rock damage around the eventual rupture zones. Analysis (ii) reveals localization (reduced fractional area) 2–3 yr before these main shocks and before the M > 7 Düzce 1999 earthquake in Turkey. Results with technique (iii) indicate that individual events tend to coalesce rapidly to clusters in the final 1–2 yr before the main shocks. Corresponding analyses of data from the Parkfield region show opposite delocalization patterns and decreasing clustering before the 2004 M6 earthquake. Continuing studies with these techniques, combined with analysis of geodetic data and insights from laboratory experiments and model simulations, might improve the ability to track preparation processes leading to large earthquakes.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. -
Abstract The magnitude‐frequency distribution (MFD) describes the relative proportion of earthquake magnitudes and provides vital information for seismic hazard assessment. The
b ‐value, derived from the MFD, is commonly used to estimate the probability that a future earthquake will exceed a specified magnitude threshold. Improved MFD andb ‐value estimates are of great importance in the central and eastern United States where high volumes of fluid injection have contributed to a significant rise in seismicity over the last decade. In this study, we recalculate the magnitudes of 8,775 events for the 2011 Prague, Oklahoma sequence using a relative magnitude approach that depends only on waveform data to calculate magnitudes. We also compare the distribution of successive magnitude differences to the MFD and show that a combination of the magnitude difference distribution (MDFD) and relative magnitudes yields a reliable estimate ofb ‐value. Using the MDFD and relative magnitudes, we examine the temporal and spatial variations in theb ‐value and show thatb ‐value ranges between ∼0.6 and 0.85 during the aftershock sequence for at least 5 months after theM 5.7 mainshock, though areas surrounding the northeast part of the sequence experience higherb ‐values (0.7–0.85) than the southwestern part of the Meeker‐Prague fault whereb ‐value is the lowest (0.6–0.7). We also identify a cluster of off‐fault events with the highestb ‐values in the catalog (0.85). These new estimates of MFD andb ‐value will contribute to understanding of the relations between induced and tectonic earthquake sequences and promote discussion regarding the use ofb ‐value in induced seismic hazard estimation.