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

A global plate motion model, named NUVEL-1, which describes current plate motions between 12 rigid plates is described, with special attention given to the method, data, and assumptions used Tectonic implications of the patterns that emerged from the results are discussed It is shown that wide plate boundary zones can form not only within the continental lithosphere but also within the oceanic lithosphere; eg, between the Indian and Australian plates and between the North American and South American plates Results of the model also suggest small but significant diffuse deformation of the oceanic lithosphere, which may be confined to small awkwardly shaped salients of major plates

Current plate motions

SUMMARY 
We determine best-fitting Euler vectors, closure-fitting Euler vectors, and a new global model (NUVEL-1) describing the geologically current motion between 12 assumed-rigid plates by inverting plate motion data we have compiled, critically analysed, and tested for self-consistency. We treat Arabia, India and Australia, and North America and South America as distinct plates, but combine Nubia and Somalia into a single African plate because motion between them could not be reliably resolved. The 1122 data from 22 plate boundaries inverted to obtain NUVEL-1 consist of 277 spreading rates, 121 transform fault azimuths, and 724 earthquake slip vectors. We determined all rates over a uniform time interval of 3.0m.y., corresponding to the centre of the anomaly 2A sequence, by comparing synthetic magnetic anomalies with observed profiles. The model fits the data well. Unlike prior global plate motion models, which systematically misfit some spreading rates in the Indian Ocean by 8–12mm yr−1, the systematic misfits by NUVEL-1 nowhere exceed ∼3 mm yr−1. The model differs significantly from prior global plate motion models. For the 30 pairs of plates sharing a common boundary, 29 of 30 P071, and 25 of 30 RM2 Euler vectors lie outside the 99 per cent confidence limits of NUVEL-1. Differences are large in the Indian Ocean where NUVEL-1 plate motion data and plate geometry differ from those used in prior studies and in the Pacific Ocean where NUVEL-1 rates are systematically 5–20 mm yr−1 slower than those of prior models. The strikes of transform faults mapped with GLORIA and Seabeam along the Mid-Atlantic Ridge greatly improve the accuracy of estimates of the direction of plate motion. These data give Euler vectors differing significantly from those of prior studies, show that motion about the Azores triple junction is consistent with plate circuit closure, and better resolve motion between North America and South America. Motion of the Caribbean plate relative to North or South America is about 7 mm yr−1 slower than in prior global models. Trench slip vectors tend to be systematically misfit wherever convergence is oblique, and best-fitting poles determined only from trench slip vectors differ significantly from their corresponding closure-fitting Euler vectors. The direction of slip in trench earthquakes tends to be between the direction of plate motion and the normal to the trench strike. Part of this bias may be due to the neglect of lateral heterogeneities of seismic velocities caused by cold subducting slabs, but the larger part is likely caused by independent motion of fore-arc crust and lithosphere relative to the overriding plate.

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

Contemporary in situ tectonic stress indicators along the San Andreas fault system in central California show northeast-directed horizontal compression that is nearly perpendicular to the strike of the fault. Such compression explains recent uplift of the Coast Ranges and the numerous active reverse faults and folds that trend nearly parallel to the San Andreas and that are otherwise unexplainable in terms of strike-slip deformation. Fault-normal crustal compression in central California is proposed to result from the extremely low shear strength of the San Andreas and the slightly convergent relative motion between the Pacific and North American plates. Preliminary in situ stress data from the Cajon Pass scientific drill hole (located 3.6 kilometers northeast of the San Andreas in southern California near San Bernardino, California) are also consistent with a weak fault, as they show no right-lateral shear stress at approximately 2-kilometer depth on planes parallel to the San Andreas fault.

https://science.sciencemag.org/content/sci/238/4830/1105.full.pdf

New evidence on the state of stress of the San Andreas fault system

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

Measurements of in situ stress orientation and magnitude at the site of the Cajon Pass research borehole have been made from depths of 0.9-3.5 km using the hydraulic fracturing technique and analysis of stress-induced well bore breakouts. The results of these measurements support two important conclusions about the mechanics of crustal faulting. First, the magnitudes of measured in situ stresses indicate ratios of shear to normal stres. s on favorably oriented fault planes that are consistent with predictions based on Mohr-Coulomb theory and laboratory-determined coefficients of friction in the range of 0.6-1.0 assuming hydrostatic pore pressure (this is commonly known as Byerlee's law). Thus the stress measurements indicate that the frictional strength of the crust adjacent to the San Andreas fault is high (i.e., consistent with laboratory-derived friction values) and that the level of shear stress in the crust adjacent to the San Andreas is principally controlled by its frictional strength. However, data on the orientation of maximum horizontal compression in the borehole from 1.75 to 3.5 km (N57oE + 19 o) indicate that the San Andreas must be quite weak as a complete absence of right-lateral shear stress resolved on planes parallel to the --N60oW striking San Andreas fault is observed. The lack of right-lateral shear stress on planes parallel to the San Andreas fault at this site is especially surprising as Cajon Pass is located along a section of the San Andreas which has not had a major earthquake since 1812 and is thus presumably quite "late" in the earthquake cycle. Nevertheless, both the orientation and magnitudes of stresses measured in the well are consistent with the style of active faulting in the area surrounding the drill site, most notably normal faulting and Quaternary age left-lateral slip on the Cleghorn fault that parallels the San Andreas in the vicinity of the drill site (Meisling and Weldon, 1982; Weldon, 1986; R. J. Weldon et al., unpublished report, 1981). We argue that the stress state (and Quaternary fault offsets) observed in the Cajon Pass area could exist only if the San Andreas moved at low shear stresses comparable to seismic stress drops rather than the much higher values predicted by Byerlee's law, a conclusion consistent with the lack of frictionally generated heat flow along the San Andreas system (e.g., Brune et al., 1969; Henyey and Wasserburg, 1971; Lachenbruch and Sass, 1973, 1980). Taken together, the Cajon Pass in situ stress and heat flow measurements (Lachenbruch and Sass, this issue) support a conceptual model of the San Andreas system in which the San Andreas is extremely weak with respect to the surrounding crust.

/pdf/in-situ-stress-measurements-to-3-5-km-depth-in-the-cajon-1rqivcjaws.pdf

John H. Healy

Papers

New evidence on the state of stress of the San Andreas fault system

An Experiment in Earthquake Control at Rangely, Colorado

In situ stress measurements to 3.5 km depth in the Cajon Pass Scientific Research Borehole: Implications for the mechanics of crustal faulting

Friction, faulting and «in situ» stress

Faulting and Crustal Stress at Rangely, Colorado