In November 2011, a M5.0 earthquake occurred less than a day before a M5.7 earthquake near Prague, Oklahoma, which may have promoted failure of the mainshock and thousands of aftershocks along the Wilzetta fault, including a M5.0 aftershock. The M5.0 foreshock occurred in close proximity to active fluid injection wells; fluid injection can cause a buildup of pore fluid pressure, decrease the fault strength, and may induce earthquakes. Keranen et al. [2013] links the M5.0 foreshock with fluid injection, but the relationship between the foreshock and successive events has not been investigated. Here we examine the role of coseismic Coulomb stress transfer on earthquakes that follow the M5.0 foreshock, including the M5.7 mainshock. We resolve the static Coulomb stress change onto the focal mechanism nodal plane that is most consistent with the rupture geometry of the three M ≥ 5.0 earthquakes, as well as specified receiver fault planes that reflect the regional stress orientation. We find that Coulomb stress is increased, e.g., fault failure is promoted, on the nodal planes of~60% of the events that have focal mechanism solutions, and more specifically, that the M5.0 foreshock promoted failure on the rupture plane of the M5.7 mainshock. We test our results over a range of effective coefficient of friction values. Hence, we argue that the M5.0 foreshock, induced by fluid injection, potentially triggered a cascading failure of earthquakes along the complex Wilzetta fault system.

In November 2011, three Mw ≥ 4.8 earthquakes and thousands of aftershocks occurred along the structurally complex Wilzetta fault system near Prague, Oklahoma. Previous studies suggest that wastewater injection induced a Mw 4.8 foreshock, which subsequently triggered a Mw 5.7 mainshock. We examine source properties of aftershocks with a standard Brune‐type spectral model and jointly solve for seismic moment (M0), corner frequency (f0), and kappa (κ) with an iterative Gauss‐Newton global downhill optimization method. We examine 934 earthquakes with initial moment magnitudes (Mw) between 0.33 and 4.99 based on the pseudospectral acceleration and recover reasonable M0, f0, and κ for 87 earthquakes with Mw 1.83–3.51 determined by spectral fit. We use M0 and f0 to estimate the Brune‐type stress drop, assuming a circular fault and shear‐wave velocity at the hypocentral depth of the event. Our observations suggest that stress drops range between 0.005 and 4.8 MPa with a median of 0.2 MPa (0.03–26.4 MPa with a median of 1.1 MPa for Madariaga‐type), which is significantly lower than typical eastern United States intraplate events (>10 MPa). We find that stress drops correlate weakly with hypocentral depth and magnitude. Additionally, we find the stress drops increase with time after the mainshock, although temporal variation in stress drop is difficult to separate from spatial heterogeneity and changing event locations. The overall low median stress drop suggests that the fault segments may have been primed to fail as a result of high pore fluid pressures, likely related to nearby wastewater injection.

[1] We detect and precisely locate over 9500 aftershocks that occurred in the Yuha Desert region during a 2 month period following the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake. Events are relocated using a series of absolute and relative relocation procedures that include Hypoinverse, Velest, and hypoDD. Location errors are reduced to 40 m horizontally and 120 m vertically. Aftershock locations reveal a complex pattern of faulting with en echelon fault segments trending toward the northwest, approximately parallel to the North American-Pacific plate boundary and en echelon, conjugate features trending to the northeast. The relocated seismicity is highly correlated with published surface mapping of faults that experienced triggered surface slip in response to the EMC main shock. Aftershocks occurred between 2 km and 11 km depths, consistent with previous studies of seismogenic thickness in the region. Three-dimensional analysis reveals individual and intersecting fault planes that are limited in their along-strike length. These fault planes remain distinct structures at depth, indicative of conjugate faulting, and do not appear to coalesce onto a throughgoing fault segment. We observe a complex spatiotemporal migration of aftershocks, with seismicity that jumps between individual fault segments that are active for only a few days to weeks. Aftershock rates are roughly consistent with the expected earthquake production rates of Dieterich (1994). The conjugate pattern of faulting and nonuniform aftershock migration patterns suggest that strain in the Yuha Desert is being accommodated in a complex manner.