Elizabeth S. Cochran

Potentially induced earthquakes in Oklahoma, USA: Links between wastewater injection and the 2011 Mw 5.7 earthquake sequence

Significant earthquakes are increasingly occurring within the continental interior of the United States, including five of moment magnitude (Mw) ≥ 5.0 in 2011 alone. Concurrently, the volume of fluid 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 Mw 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 Mw 5.0 and a prolific 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 fluid was injected into effectively sealed compartments, and we interpret that a net fluid volume increase after 18 yr of injection lowered effective stress on reservoir-bounding faults. Significantly, this case indicates that decades-long lags between the commencement of fluid injection and the onset of induced earthquakes are possible, and modifies our common criteria for fluid-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 first.

Seismic and geodetic evidence for extensive, long-lived fault damage zones

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.

Seismic Evidence for Rock Damage and Healing on the San Andreas Fault Associated with the 2004 M 6.0 Parkfield Earthquake

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 exper- iment 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 cross- correlations 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 promi- nent at depths above 7 km.

Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence

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.

Postseismic Fault Healing on the Rupture Zone of the 1999 M 7.1 Hector Mine, California, Earthquake

We probed the rupture zone of the October 1999 M 7.1 Hector Mine earthquake using repeated near-surface explosions in October 2000 and November 2001. Three dense linear seismic arrays were deployed across the north and south Lavic Lake faults (LLFs) that broke to the surface in the mainshock and across the Bullion fault (BF) that experienced minor slip in that event. Two explosions each year were detonated in the rupture zone, one on the middle and one on the south LLF. We found that P and S velocities of fault-zone rocks increased by 0.7%-1.4% and 0.5%-1.0% between 2000 and 2001, respectively. In contrast, the velocities for P and S waves in surrounding rocks increased much less. This trend indicates that the Hector Mine rupture zone has been healing by strengthening after the main- shock, most likely due to the closure of cracks that opened during the 1999 earth- quake. The observed fault-zone strength recovery is consistent with an apparent crack density decrease of 1.5% within the rupture zone. The ratio of travel-time decrease for P to S waves was 0.72, consistent with partially fluid-filled cracks near the fault zone. This restrengthening is similar to that observed after the 1992 M 7.4 Landers earthquake, which occurred 25 km to the west (Li and Vidale, 2001). We also find that the velocity increase with time varies from one fault segment to another at the Hector Mine rupture zone. We see greater changes on the LLFs than on the BF, and the greatest change is on the middle LLF at shallow depth. We tentatively conclude that greater damage was inflicted, and thus greater healing is observed, in regions with larger slip in the mainshock.

Anisotropy in the Shallow Crust Observed around the San Andreas Fault Before and After the 2004 M 6.0 Parkfield Earthquake

Local seismic arrays were deployed at two locations along the San An- dreas fault (SAF) near Parkfield, California, before and after the 2004 M 6.0 Parkfield earthquake. Using local earthquakes we determine the anisotropic field within 1- 2 km of the main trace of the SAF at the two array locations separated by 12 km. The initial array, near the SAFOD site, was deployed for six weeks in October and November 2003, and the second array, located near the town of Parkfield, was de- ployed for 3 months following the 28 September 2004 M 6.0 Parkfield earthquake. We find the fast shear-wave polarization direction nearly fault-parallel (N40W) for stations on the main fault trace and within 100 m to the southwest of the SAF at both array locations. These fault-parallel measurements span the 100- to 150-m-wide zone of pervasive cracking and damage interpreted from fault-zone-trapped waves associated with the main fault core (Li et al., 2004, 2006). Outside of this zone, the fast orientations are scattered with some preference for orientations near N10E, roughly parallel to the regional maximum horizontal compressive stress direction (rh). In addition, fast directions are preferentially oriented parallel to a northern branch of the SAF recorded on stations in the 2004 Parkfield deployment. The measured anisotropy is likely due to a combination of stress-aligned micro- cracks away from the fault and shear fabric within the highly evolved fault core. The majority of our measurements are taken outside of the main fault core, and we es- timate the density of microcracks from the measured delay times. Apparent crack densities are approximately 3%, with large scatter. The data suggest weak depth dependence to the measured delay times for source depths between 2 and 7 km. Below 7-km source depth, the delay times do not correlate with depth suggesting higher confining pressure is forcing the microcracks to close. No coseismic variation in the anisotropic parameters is observed, suggesting little to no influence on measured splitting due to the 2004 M 6.0 Parkfield earthquake. However, the premainshock and postmainshock data presented here are from arrays separated by 12 km, limiting our sensitivity to small temporal changes in anisotropy.

A novel strong-motion seismic network for community participation in earthquake monitoring

The Quake-Catcher Network (QCN) is breaking new ground in seismology by combining new micro-electro-mechanical systems (MEMS) technology with volunteer seismic station distributed computing. Rather than distributing just computations, the QCN allows volunteers to participate in scientific data collection and computation. Using these innovative tools, QCN will increase the number of strong-motion observations for improved earthquake detection and analysis in California, and throughout the world. QCNs increased density of seismic measurements will revolutionize seismology. The QCN, in concert with current seismic networks, may soon provide advanced alerts when earthquakes occur, estimate the response of a building to earthquakes even before they happen, and generate a greater understanding of earthquakes for scientists and the general public alike. Details of how one can join the QCN are outlined. In addition, we have activities on our website that can be used in K-16 classrooms to teach students basic seismology and physics concepts.

Seismic velocity variations on the San Andreas fault caused by the 2004 M6 Parkfield Earthquake and their implications

Repeated earthquakes and explosions recorded at the San Andreas fault (SAF) near Parkfield before and after the 2004 M6 Parkfield earthquake show large seismic velocity variations within an approximately 200- m-wide zone along the fault to depths of approximately 6 km. The seismic arrays were co-sited in the two experiments and located in the middle of a high-slip part of the surface rupture. Waveform cross-correlations of microearthquakes recorded in 2002 and subsequent repeated events recorded a week after the 2004 M6 mainshock show a peak of an approximately 2.5% decrease in seismic velocity at stations within the fault zone, most likely due to the co-seismic damage of fault-zone rocks during dynamic rupture of this earthquake. The damage zone is not symmetric; instead, it extends farther on the southwest side of the main fault trace. Seismic velocities within the fault zone measured for later repeated aftershocks in the following 3–4 months show an approximate 1.2% increase at seismogenic depths, indicating that the rock damaged in the mainshock recovers rigidity—or heals—through time. The healing rate was not constant but was largest in the earliest post-mainshock stage. The magnitude of fault damage and healing varies across and along the rupture zone, indicating that the greater damage was inflicted and thus greater healing is observed in regions with larger slip in the mainshock. Observations of rock damage during the mainshock and healing soon thereafter are consistent with our interpretation of the low-velocity waveguide on the SAF being at least partially softened in the 2004 M6 mainshock, with additional cumulative effects due to recurrent rupture.

Strong-motion observations of the M 7.8 Gorkha, Nepal, earthquake sequence and development of the N-shake strong-motion network

We present and describe strong-motion data observations from the 2015 M 7.8 Gorkha, Nepal, earthquake sequence collected using existing and new Quake-Catcher Network (QCN) and U.S. Geological Survey NetQuakes sensors located in the Kathmandu Valley. A comparison of QCN data with waveforms recorded by a conventional strong-motion (NetQuakes) instrument validates the QCN data. We present preliminary analysis of spectral accelerations, and peak ground acceleration and velocity for earthquakes up to M 7.3 from the QCN stations, as well as preliminary analysis of the mainshock recording from the NetQuakes station. We show that mainshock peak accelerations were lower than expected and conclude the Kathmandu Valley experienced a pervasively nonlinear response during the mainshock. Phase picks from the QCN and NetQuakes data are also used to improve aftershock locations. This study confirms the utility of QCN instruments to contribute to ground-motion investigations and aftershock response in regions where conventional instrumentation and open-access seismic data are limited. Initial pilot installations of QCN instruments in 2014 are now being expanded to create the Nepal–Shaking Hazard Assessment for Kathmandu and its Environment (N-SHAKE) network. Online Material: Figures of Pg arrivals, earthquake locations, epicenter change vectors, and travel-time misfit vector residuals, and tables of QCN and NetQuake stations and relocated hypocenter timing, location, and magnitude.

Performance of Several Low-Cost Accelerometers

Several groups are implementing low-cost host-operated systems of strong-motion accelerographs to support the somewhat divergent needs of seismologists and earthquake engineers. The Advanced National Seismic System Technical Implementation Committee (ANSS TIC, 2002), managed by the U.S. Geological Survey (USGS) in cooperation with other network operators, is exploring the efficacy of such systems if used in ANSS networks. To this end, ANSS convened a working group to explore available Class C strong-motion accelerometers (defined later), and to consider operational and quality control issues, and the means of annotating, storing, and using such data in ANSS networks. The working group members are largely coincident with our author list, and this report informs instrument-performance matters in the working group’s report to ANSS. Present examples of operational networks of such devices are the Community Seismic Network (CSN; csn.caltech.edu), operated by the California Institute of Technology, and Quake-Catcher Network (QCN; Cochran et al., 2009; qcn.stanford.edu; November 2013), jointly operated by Stanford University and the USGS. Several similar efforts are in development at other institutions. The overarching goals of such efforts are to add spatial density to existing Class-A and Class-B (see next paragraph) networks at low cost, and to include many additional people so they become invested in the issues of earthquakes, their measurement, and the damage they cause.