skip to main content


Title: New Estimates of Magnitude‐Frequency Distribution and b ‐Value Using Relative Magnitudes for the 2011 Prague, Oklahoma Earthquake Sequence
Abstract

The magnitude‐frequency distribution (MFD) describes the relative proportion of earthquake magnitudes and provides vital information for seismic hazard assessment. Theb‐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 theM5.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.

 
more » « less
NSF-PAR ID:
10483831
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
129
Issue:
1
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We conduct a detailed study of the foreshock sequence preceding the 2010Mw6.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 aMw4.6 foreshock approximately 2 hr before theMw6.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 theMw4.6 earthquake. We further invert for the focal mechanisms and rupture directions for both the largest foreshock and the mainshock. TheMw4.6 foreshock likely occurred on a NE‐SW trending fault conjugating to the NW‐SE trending fault of the mainshock. Coulomb stress analysis suggests theMw4.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.

     
    more » « less
  2. How structural segment boundaries modulate earthquake behavior is an important scientific and societal question, especially for the Wasatch fault zone (WFZ) where urban areas lie along multiple fault segments. The extent to which segment boundaries arrest ruptures, host moderate magnitude earthquakes, or transmit ruptures to adjacent fault segments is critical for understanding seismic hazard. To help address this outstanding issue, we conducted a paleoseismic investigation at the Traverse Ridge paleoseismic site (TR site) along the ∼7-km-long Fort Canyon segment boundary, which links the Provo (59 km) and Salt Lake City (40 km) segments of the WFZ. At the TR site, we logged two trenches which were cut across sub-parallel traces of the fault, separated by ∼175 m. Evidence from these exposures leads us to infer that at least 3 to 4 earthquakes have ruptured across the segment boundary in the Holocene. Radiocarbon dating of soil material developed below and above fault scarp colluvial packages and within a filled fissure constrains the age of the events. The most recent event ruptured the southern fault trace between 0.2 and 0.4 ka, the penultimate event ruptured the northern fault trace between 0.6 and 3.4 ka, and two prior events occurred between 1.4 and 6.2 ka (on the southern fault trace) and 7.2 and 8.1 ka (northern fault trace). Colluvial wedge heights of these events ranged from 0.7 to 1.2 m, indicating the segment boundary experiences surface ruptures with more than 1 m of vertical displacement. Given these estimates, we infer that these events were greater than Mw 6.7, with rupture extending across the entire segment boundary and portions of one or both adjacent fault segments. The Holocene recurrence of events at the TR site is lower than the closest paleoseismic sites at the adjacent fault segment endpoints. The contrasts in recurrence rates observed within 15 km of the Fort Canyon fault segment boundary may be explained conceptually by a leaky segment boundary model which permits spillover events, ruptures centered on the segment boundary, and segmented ruptures. The TR site demonstrates the utility of paleoseismology within segment boundaries which, through corroboration of displacement data, can demonstrate rupture connectivity between fault segments and test the validity of rupture models. 
    more » « less
  3. 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.

     
    more » « less
  4. Abstract

    Dike swarms are the fossil remains of regions of the crust that have undergone repeated magma injections. Volcanic earthquake swarms and geodetic measurements are, at least in part, a record of active injection of fluids (water, gas, or magma) into fractures. Here, we link these two ways of observing magmatic systems by noting that dike thicknesses and earthquake magnitudes share similar scaling parameters. In the Jurassic Independence dike swarm of eastern California median dike thickness is ∼1 m, similar to other swarms worldwide, but glacially polished exposures reveal that a typical dike comprises a number of dikelets that are lognormally distributed in thickness with a mean of ∼0.1 m. Assuming that dikes fill penny‐shaped cracks of a given aspect ratio, the geodetic moment and earthquake magnitude of a diking event can be estimated. A Monte Carlo simulation of dike‐induced earthquakes based on observed dike thickness variations yields a frequency‐magnitude distribution remarkably like observed volcanic earthquake swarms, with ab‐value of ∼1.7. We suggest that swarms of dikes composed of dikelets, as well as plutons built incrementally by sheet intrusions, are physical complements to volcanic seismic swarms, and that at least some earthquake swarms are a palpable expression of incremental magma emplacement.

     
    more » « less
  5. Abstract

    Understanding earthquake foreshocks is essential for deciphering earthquake rupture physics and can aid seismic hazard mitigation. With regional dense seismic arrays, we identify immediate foreshocks of 527 0.9M5.4 events of the 2019 Ridgecrest earthquake sequence, including 48 earthquakes with series of immediate foreshocks. These immediate foreshocks are adjacent to the mainshocks occurring within 100 s of the mainshocks, and their P waves share high resemblances with the mainshock P waves. However, attributes of the immediate‐foreshock P waves, including the amplitudes and preceding times, do not clearly scale with the mainshock magnitudes. Our observations suggest that earthquake rupture may initiate in a universal fashion but evolves stochastically. This indicates that earthquake rupture development is likely controlled by fine‐scale fault heterogeneities in the Ridgecrest fault system, and the final magnitude is the only difference between small and large earthquakes.

     
    more » « less