This essential reference for graduate students and researchers provides a unified treatment of earthquakes and faulting as two aspects of brittle tectonics at different timescales. The intimate connection between the two is manifested in their scaling laws and populations, which evolve from fracture growth and interactions between fractures. The connection between faults and the seismicity generated is governed by the rate and state dependent friction laws - producing distinctive seismic styles of faulting and a gamut of earthquake phenomena including aftershocks, afterslip, earthquake triggering, and slow slip events. The third edition of this classic treatise presents a wealth of new topics and new observations. These include slow earthquake phenomena; friction of phyllosilicates, and at high sliding velocities; fault structures; relative roles of strong and seismogenic versus weak and creeping faults; dynamic triggering of earthquakes; oceanic earthquakes; megathrust earthquakes in subduction zones; deep earthquakes; and new observations of earthquake precursory phenomena.

The Mechanics of Earthquakes and Faulting

Direct shear experiments on ground surfaces of a granodiorite from Raymond, California, at normal stresses of ∼6 MPa demonstrate that competing time, displacement, and velocity effects control rock friction. It is proposed that the strength of the population of points of contacts between sliding surfaces determines frictional strength and that the population of contacts changes continuously with displacements. Previous experiments demonstrate that the strength of the contacts increases with the age of the contacts. The present experiments establish that a characteristic displacement, proportional to surface roughness, is required to change the population of contacts. Hence during slip the average age of the points of contact and therefore frictional strength decrease as slip velocity increases. Displacement weakening and consequently the potential for unstable slip occur whenever displacement reduces the average age of the contacts. In addition to this velocity dependency, which arises from displacement dependency and time dependency, the experiments also show a competing but transient increase in friction whenever slip velocity increases. Creep of the sliding surface at stresses below that for steady state slip is also observed. Constitutive relationships are developed that permit quantitative simulation of the friction versus displacement data as a function of surface roughness and for different time and velocity histories. Unstable slip in experiments is controlled by these constitutive effects and by the stiffness of the experimental system. It is argued that analogous properties control earthquake instability.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/JB084iB05p02161%4010.1002/%28ISSN%292169-9356.FAULT1

Modeling of rock friction: 1. Experimental results and constitutive equations

Background Human-induced earthquakes have become an important topic of political and scientific discussion, owing to the concern that these events may be responsible for widespread damage and an overall increase in seismicity. It has long been known that impoundment of reservoirs, surface and underground mining, withdrawal of fluids and gas from the subsurface, and injection of fluids into underground formations are capable of inducing earthquakes. In particular, earthquakes caused by injection have become a focal point, as new drilling and well-completion technologies enable the extraction of oil and gas from previously unproductive formations. Earthquakes with magnitude (M) ≥ 3 in the U.S. midcontinent, 1967–2012. After decades of a steady earthquake rate (average of 21 events/year), activity increased starting in 2001 and peaked at 188 earthquakes in 2011. Human-induced earthquakes are suspected to be partially responsible for the increase. Advances Microearthquakes (that is, those with magnitudes below 2) are routinely produced as part of the hydraulic fracturing (or “fracking”) process used to stimulate the production of oil, but the process as currently practiced appears to pose a low risk of inducing destructive earthquakes. More than 100,000 wells have been subjected to fracking in recent years, and the largest induced earthquake was magnitude 3.6, which is too small to pose a serious risk. Yet, wastewater disposal by injection into deep wells poses a higher risk, because this practice can induce larger earthquakes. For example, several of the largest earthquakes in the U.S. midcontinent in 2011 and 2012 may have been triggered by nearby disposal wells. The largest of these was a magnitude 5.6 event in central Oklahoma that destroyed 14 homes and injured two people. The mechanism responsible for inducing these events appears to be the well-understood process of weakening a preexisting fault by elevating the fluid pressure. However, only a small fraction of the more than 30,000 wastewater disposal wells appears to be problematic—typically those that dispose of very large volumes of water and/or communicate pressure perturbations directly into basement faults. Outlook Injection-induced earthquakes, such as those that struck in 2011, clearly contribute to the seismic hazard. Quantifying their contribution presents difficult challenges that will require new research into the physics of induced earthquakes and the potential for inducing large-magnitude events. The petroleum industry needs clear requirements for operation, regulators must have a solid scientific basis for those requirements, and the public needs assurance that the regulations are sufficient and are being followed. The current regulatory frameworks for wastewater disposal wells were designed to protect potable water sources from contamination and do not address seismic safety. One consequence is that both the quantity and timeliness of information on injection volumes and pressures reported to regulatory agencies are far from ideal for managing earthquake risk from injection activities. In addition, seismic monitoring capabilities in many of the areas in which wastewater injection activities have increased are not capable of detecting small earthquake activity that may presage larger seismic events.

http://science.sciencemag.org/content/sci/341/6142/1225942.full.pdf

Injection-Induced Earthquakes

This paper reviews rock friction and the frictional properties of earthquake faults. The basis for rate- and state-dependent friction laws is reviewed. The friction state variable is discussed, including its interpretation as a measure of average asperity contact time and porosity within granular fault gouge. Data are summarized showing that friction evolves even during truly stationary contact, and the connection between modern friction laws and the concept of “static” friction is discussed. Measurements of frictional healing, as evidenced by increasing static friction during quasistationary contact, are reviewed, as are their implications for fault healing. Shear localization in fault gouge is discussed, and the relationship between microstructures and friction is reviewed. These data indicate differences in the behavior of bare rock surfaces as compared to shear within granular fault gouge that can be attributed to dilation within fault gouge. Physical models for the characteristic friction distance are discussed and related to the problem of scaling this parameter to seismic faults. Earthquake afterslip, its relation to laboratory friction data, and the inverse correlation between afterslip and shallow coseismic slip are discussed in the context of a model for afterslip. Recent observations of the absence of afterslip are predicted by the model.

/pdf/laboratory-derived-friction-laws-and-their-application-to-44s35v6won.pdf

Laboratory-derived friction laws and their application to seismic faulting

To date, more than 7300 in situ stress orientations have been compiled as part of the World Stress Map project. Of these, over 4400 are considered reliable tectonic stress indicators, recording horizontal stress orientations to within <±25°. Remarkably good correlation is observed between stress orientations deduced from in situ stress measurements and geologic observations made in the upper 1–2 km, well bore breakouts extending to 4–5 km depth and earthquake focal mechanisms to depths of ∼20 km. Regionally uniform stress orientations and relative magnitudes permit definition of broad-scale regional stress patterns often extending 20–200 times the approximately 20–25 km thickness of the upper brittle lithosphere. The “first-order” midplate stress fields are believed to be largely the result of compressional forces applied at plate boundaries, primarily ridge push and continental collision. The orientation of the intraplate stress field is thus largely controlled by the geometry of the plate boundaries. There is no evidence of large lateral stress gradients (as evidenced by lateral variations in stress regime) which would be expected across large plates if simple resistive or driving basal drag tractions (parallel or antiparallel to absolute motion) controlled the intraplate stress field. Intraplate areas of active extension are generally associated with regions of high topography: western U.S. Cordillera, high Andes, Tibetan plateau, western Indian Ocean plateau. Buoyancy stresses related to crustal thickening and/or lithospheric thinning in these regions dominate the intraplate compressional stress field due to plate-driving forces. These buoyancy forces are just one of several categories of “second-order” stresses, or local perturbations, that can be identified once the first-order stress patterns are recognized. These second-order stress fields can often be associated with specific geologic or tectonic features, for example, lithospheric flexure, lateral strength contrasts, as well as the lateral density contrasts which give rise to buoyancy forces. These second-order stress patterns typically have wavelengths ranging from 5 to 10+ times the thickness of the brittle upper lithosphere. A two-dimensional analysis of the amount of rotation of regional horizontal stress orientations due to a superimposed local stress constrains the ratio of the magnitude of the horizontal regional stress differences to the local uniaxial stress. For a detectable rotation of 15°, the local horizontal uniaxial stress must be at least twice the magnitude of the regional horizontal stress differences. Examples of local rotations of SHmax orientations include a 75°–85° rotation on the northeastern Canadian continental shelf possibly related to margin-normal extension derived from sediment-loading flexural stresses, a 50°–60° rotation within the East African rift relative to western Africa due to extensional buoyancy forces caused by lithospheric thinning, and an approximately 90° rotation along the northern margin of the Paleozoic Amazonas rift in central Brazil. In this final example, this rotation is hypothesized as being due to deviatoric compression oriented normal to the rift axis resulting from local lithospheric support of a dense mass in the lower crust beneath the rift (“rift pillow”). Estimates of the magnitudes of first-order (plate boundary force-derived) regional stress differences computed from modeling the source of observed local stress rotations magnitudes can be compared with regional stress differences based on the frictional strength of the crust (i.e., “Byerlee's law”) assuming hydrostatic pore pressure. The examples given here are too few to provide a definitive evaluation of the direct applicability of Byerlee's law to the upper brittle part of the lithosphere, particularly in view of uncertainties such as pore pressure and relative magnitude of the intermediate principal stresses. Nonetheless, the observed rotations all indicate that the magnitude of the local horizontal uniaxial stresses must be 1–2.5+ times the magnitude of the regional first-order horizontal stress differences and suggest that careful evaluation of such local rotations may be a powerful technique for constraining the in situ magnitude stress differences in the upper, brittle part of the lithosphere.

/pdf/first-and-second-order-patterns-of-stress-in-the-lithosphere-595hs86ljh.pdf

First‐ and second‐order patterns of stress in the lithosphere: The World Stress Map Project

An experiment in an oil field at Rangely, Colorado, has demonstrated the feasibility of earthquake control. Variations in seismicity were produced by controlled variations in the fluid pressure in a seismically active zone. Precise earthquake locations revealed that the earthquakes clustered about a fault trending through a zone of high pore pressure produced by secondary recovery operations. Laboratory measurements of the frictional properties of the reservoir rocks and an in situ stress measurement made near the earthquake zone were used to predict the fluid pressure required to trigger earthquakes on preexisting fractures. Fluid pressure was controlled by alternately injecting and recovering water from wells that penetrated the seismic zone. Fluid pressure was monitored in observation wells, and a computer model of the reservoir was used to infer the fluid pressure distributions in the vicinity of the injection wells. The results of this experiment confirm the predicted effect of fluid pressure on earthquake activity and indicate that earthquakes can be controlled wherever we can control the fluid pressure in a fault zone.

https://science.sciencemag.org/content/sci/191/4233/1230.full.pdf

An Experiment in Earthquake Control at Rangely, Colorado

Faulting and Crustal Stress at Rangely, Colorado

Earthquakes have been recorded from the vicinity of Rangely, Colo., since construction of the Uinta Basin Seismological Observatory in 1962. Waterflooding by injection into the Rangely oil field began in 1958. From subhydrostatic pore pressures in 1958, the reservoir fluid pressure in some peripheral parts of the field has been raised to 0.6 of the overburden pressure. The earthquakes have been precisely located since Sept. 1969, by a network of 14 stations. The principal seismic activity occurs over a depth interval of 1 to 4 km along a vertical fault which cuts across the field. The earthquake epicenters are located on the S. margin and die out toward the interior of the field, disappearing where the pore pressure is at or below the normal hydrostatic gradient. Focal plane solutions idicate right-lateral strike-slip faulting. Numerical modeling of the reservoir pressure history indicates that pumping from 4 wells in the hypocentral region will draw the pressures down to approaching normal hydrostatic in one year in an area larger than the error in hypocentral location.

J. D. Bredehoeft

Papers

An Experiment in Earthquake Control at Rangely, Colorado

Faulting and Crustal Stress at Rangely, Colorado

Earthquakes and waterflooding in the Rangely oil field