Signifi cant earthquakes are increasingly occurring within the continental interior of the United States, including fi ve of moment magnitude (M w) ≥ 5.0 in 2011 alone. Concurrently, the volume of fl uid injected into the subsurface related to the production of unconventional resources continues to rise. Here we identify the largest earthquake potentially related to injection, an M w 5.7 earthquake in November 2011 in Oklahoma. The earthquake was felt in at least 17 states and caused damage in the epicentral region. It occurred in a sequence, with 2 earthquakes of M w 5.0 and a prolifi c sequence of aftershocks. We use the aftershocks to illuminate the faults that ruptured in the sequence, and show that the tip of the initial rupture plane is within ~200 m of active injection wells and within ~1 km of the surface; 30% of early aftershocks occur within the sedimentary section. Subsurface data indicate that fl uid was injected into effectively sealed compartments, and we interpret that a net fl uid volume increase after 18 yr of injection lowered effective stress on reservoir-bounding faults. Signifi cantly, this case indicates that decades-long lags between the commencement of fl uid injection and the onset of induced earthquakes are possible, and modifi es our common criteria for fl uid-induced events. The progressive rupture of three fault planes in this sequence suggests that stress changes from the initial rupture triggered the successive earthquakes, including one larger than the fi rst.

We show a correlation between the occurrence of shallow thrust earthquakes and the occurrence of the strongest tides. The rate of earthquakes varies from the background rate by a factor of 3 with the tidal stress. The highest correlation is found when we assume a coefficient of friction of mu = 0.4 for the crust, although we see good correlation for mu between 0.2 and 0.6. Our results quantify the effect of applied stress on earthquake triggering, a key factor in understanding earthquake nucleation and cascades whereby one earthquake triggers others.

During earthquakes, slip is often localized on preexisting faults, but it is not well understood how the structure of crustal faults may contribute to slip localization and energetics. Growing evidence suggests that the crust along active faults undergoes anomalous strain and damage during large earthquakes. Seismic and geodetic data from the Calico fault in the eastern California shear zone reveal a wide zone of reduced seismic velocities and effective elastic moduli. Using seismic traveltimes, trapped waves, and interferometric synthetic aperture radar observations, we document seismic velocities reduced by 40%-50% and shear moduli reduced by 65% compared to wall rock in a 1.5-km-wide zone along the Calico fault. Observed velocity reductions likely represent the cumulative mechanical damage from past earthquake ruptures. No large earthquake has broken the Calico fault in historic time, implying that fault damage persists for hundreds or perhaps thousands of years. These fi ndings indicate that faults can affect rock properties at substantial distances from primary fault slip surfaces, and throughout much of the seismogenic zone, a result with implications for the amount of energy expended during rupture to drive cracking and yielding of rock and development of fault systems.

We deployed a dense linear array of 45 seismometers across and along the San Andreas fault near Parkfield a week after the M 6.0 Parkfield earthquake on 28 September 2004 to record fault-zone seismic waves generated by aftershocks and explosions. Seismic stations and explosions were co-sited with our previous experiment conducted in 2002. The data from repeated shots detonated in the fall of 2002 and 3 months after the 2004 M 6.0 mainshock show ϳ1.0%-1.5% decreases in seismic-wave velocity within an ϳ200-m-wide zone along the fault strike and smaller changes (0.2%-0.5%) beyond this zone, most likely due to the coseismic damage of rocks during dynamic rupture in the 2004 M 6.0 earthquake. The width of the damage zone characterized by larger velocity changes is consistent with the low-velocity waveguide model on the San Andreas fault, near Parkfield, that we derived from fault-zone trapped waves (Li et al., 2004). The damage zone is not symmetric but extends farther on the southwest side of the main fault trace. Waveform crosscorrelations for repeated aftershocks in 21 clusters, with a total of ϳ130 events, located at different depths and distances from the array site show ϳ0.7%-1.1% increases in S-wave velocity within the fault zone in 3 months starting a week after the earthquake. The velocity recovery indicates that the damaged rock has been healing and regaining the strength through rigidity recovery with time, most likely due to the closure of cracks opened during the mainshock. We estimate that the net decrease in seismic velocities within the fault zone was at least ϳ2.5%, caused by the 2004 M 6.0 Parkfield earthquake. The healing rate was largest in the earlier stage of the postmainshock healing process. The magnitude of fault healing varies along the rupture zone, being slightly larger for the healing beneath Middle Mountain, correlating well with an area of large mapped slip. The fault healing is most prominent at depths above ϳ7 km.

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.

The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. We leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.

No abstract

We used dense linear seismic arrays across and along the San Andreas Fault (SAF) at Parkfield, California to record fault zone trapped waves generated by explosions and microearthquakes in 2002. Prominent trapped waves appeared at stations close to the SAF main fault trace while some energy was trapped in the north strand at the array site. Observations and 3‐D finite‐difference simulations of trapped waves at 2–5 Hz show evidence of a damaged core zone on the main SAF. The zone from the surface to seismogenic depths is marked by a low‐velocity waveguide ∼150 m wide, in which Q is 10–50 and shear velocities are reduced by 30–40% from wall‐rock velocities, with the greatest velocity reduction at shallow depth. We interpret that this distinct low‐velocity zone on the main SAF is a remanent of damage due to past large earthquakes on the principal fault plane at Parkfield. A less‐developed low‐velocity zone may be evident on the north strand that experienced minor breaks in the 1966 M6 event.