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

To understand whether the 1992 M = 7.4 Landers earthquake changed the proximity to failure on the San Andreas fault system, we examine the general problem of how one earthquake might trigger another. The tendency of rocks to fail in a brittle manner is thought to be a function of both shear and confining stresses, commonly formulated as the Coulomb failure criterion. Here we explore how changes in Coulomb conditions associated with one or more earthquakes may trigger subsequent events. We first consider a Coulomb criterion appropriate for the production of aftershocks, where faults most likely to slip are those optimally orientated for failure as a result of the prevailing regional stress field and the stress change caused by the mainshock. We find that the distribution of aftershocks for the Landers earthquake, as well as for several other moderate events in its vicinity, can be explained by the Coulomb criterion as follows: aftershocks are abundant where the Coulomb stress on optimally orientated faults rose by more than one-half bar, and aftershocks are sparse where the Coulomb stress dropped by a similar amount. Further, we find that several moderate shocks raised the stress at the future Landers epicenter and along much of the Landers rupture zone by about a bar, advancing the Landers shock by 1 to 3 centuries. The Landers rupture, in turn, raised the stress at site of the future M = 6.5 Big Bear aftershock site by 3 bars. The Coulomb stress change on a specified fault is independent of regional stress but depends on the fault geometry, sense of slip, and the coefficient of friction. We use this method to resolve stress changes on the San Andreas and San Jacinto faults imposed by the Landers sequence. Together the Landers and Big Bear earthquakes raised the stress along the San Bernardino segment of the southern San Andreas fault by 2 to 6 bars, hastening the next great earthquake there by about a decade.

/pdf/static-stress-changes-and-the-triggering-of-earthquakes-4r29mmr407.pdf

Static stress changes and the triggering of earthquakes

The traditional view of tectonics is that the lithosphere comprises a strong brittle layer overlying a weak ductile layer, which gives rise to two forms of deformation: brittle fracture, accompanied by earth- quakes, in the upper layer, and aseismic ductile flow in the layer beneath Although this view is not incorrect, it is imprecise, and in ways that can lead to serious misunderstandings The term ductility, for example, can apply equally to two common rock deformation mechanisms: crystal plasticity, which occurs in rock above a critical temperature, and cataclastic flow, a type of granular deformation which can occur in poorly consolidated sediments Although both exhibit ductility, these two deformation mechanisms have very different rheologies Earthquakes, in turn, are associated with strength and brittleness—associations that are likewise sufficiently imprecise that, if taken much beyond the generality implied in the opening sentence, they can lead to serious misinterpretations about earthquake mechanics Lately, a newer, much more precise and predictive model for the earthquake mechanism has emerged, which has its roots in the observation that tectonic earthquakes seldom if ever occur by the sudden appearance and propagation of a new shear crack (or 'fault') Instead, they occur by sudden slippage along a pre-existing fault or plate interface They are therefore a frictional, rather than fracture, phenomenon, with brittle fracture playing a secondary role in the lengthening of faults 1 and frictional wear 2  This distinction was noted by several early workers 3 , but it was not until 1966 that Brace and Byerlee 4 pointed out that earthquakes must be the result of a stick-slip frictional instability Thus, the earthquake is the 'slip', and the 'stick' is the interseismic period of elastic strain accumula- tion Subsequently, a complete constitutive law for rock friction has been developed based on laboratory studies A surprising result is that a great many other aspects of earthquake phenomena also now seem to result from the nature of the friction on faults The properties traditionally thought to control these processes— strength, brittleness and ductility—are subsumed within the over- arching concept of frictional stability regimes Constitutive law of rock friction In the standard model of stick-slip friction it is assumed that sliding begins when the ratio of shear to normal stress on the surface reaches a value ms, the static friction coefficient Once sliding initiates, frictional resistance falls to a lower dynamic friction coefficient, md, and this weakening of sliding resistance may,

/pdf/earthquakes-and-friction-laws-22mbgrykdx.pdf

Earthquakes and friction laws

The concept of self-organized criticality was introduced to explain the behaviour of the sandpile model. In this model, particles are randomly dropped onto a square grid of boxes. When a box accumulates four particles they are redistributed to the four adjacent boxes or lost off the edge of the grid. Redistributions can lead to further instabilities with the possibility of more particles being lost from the grid, contributing to the size of each ‘avalanche’. These model ‘avalanches’ satisfied a power-law frequency‐area distribution with a slope near unity. Other cellular-automata models, including the slider-block and forest-fire models, are also said to exhibit self-organized critical behaviour. It has been argued that earthquakes, landslides, forest fires, and species extinctions are examples of self-organized criticality in nature. In addition, wars and stock market crashes have been associated with this behaviour. The forest-fire model is particularly interesting in terms of its relation to the critical-point behaviour of the sitepercolation model. In the basic forest-fire model, trees are randomly planted on a grid of points. Periodically in time, sparks are randomly dropped on the grid. If a spark drops on a tree, that tree and adjacent trees burn in a model fire. The fires are the ‘avalanches’ and they are found to satisfy power-law frequency‐area distributions with slopes near unity. This forest-fire model is closely related to the site-percolation model, that exhibits critical behaviour. In the forest-fire model there is an inverse cascade of trees from small clusters to large clusters, trees are lost primarily from model fires that destroy the largest clusters. This quasi steady-state cascade gives a power-law frequency‐area distribution for both clusters of trees and smaller fires. The site-percolation model is equivalent to the forest-fire model without fires. In this case there is a transient cascade of trees from small to large clusters and a power-law distribution is found only at a critical density of trees.

/pdf/self-organized-criticality-4v87dkoa15.pdf

Self-organized criticality

Characterizing flow and transport in fractured geological media: A review

Scaling in fracture systems has become an active field of research in the last 25 years motivated by practical applications in hazardous waste disposal, hy- drocarbon reservoir management, and earthquake haz- ard assessment. Relevant publications are therefore spread widely through the literature. Although it is rec- ognized that some fracture systems are best described by scale-limited laws (lognormal, exponential), it is now recognized that power laws and fractal geometry provide widely applicable descriptive tools for fracture system characterization. A key argument for power law and fractal scaling is the absence of characteristic length scales in the fracture growth process. All power law and fractal characteristics in nature must have upper and lower bounds. This topic has been largely neglected, but recent studies emphasize the importance of layering on all scales in limiting the scaling characteristics of natural fracture systems. The determination of power law expo- nents and fractal dimensions from observations, al- though outwardly simple, is problematic, and uncritical use of analysis techniques has resulted in inaccurate and even meaningless exponents. We review these tech- niques and suggest guidelines for the accurate and ob- jective estimation of exponents and fractal dimensions. Syntheses of length, displacement, aperture power law exponents, and fractal dimensions are found, after crit- ical appraisal of published studies, to show a wide vari- ation, frequently spanning the theoretically possible range. Extrapolations from one dimension to two and from two dimensions to three are found to be nontrivial, and simple laws must be used with caution. Directions for future research include improved techniques for gathering data sets over great scale ranges and more rigorous application of existing analysis methods. More data are needed on joints and veins to illuminate the differences between different fracture modes. The phys- ical causes of power law scaling and variation in expo- nents and fractal dimensions are still poorly understood.

/pdf/scaling-of-fracture-systems-in-geological-media-1uq6op1ap8.pdf

Scaling of fracture systems in geological media

Physical explanation for the displacement-length relationship of faults using a post-yield fracture mechanics model

Displacement-length scaling relationship for faults: data synthesis and discussion

We report progress made in the last few years on the general problem of the mechanism of fault growth and the scaling laws that result Results are now conclusive that fault growth is a self-similar process in which fault displacement d scales linearly with fault length L Both this result and the overall nature of along-strike fault displacement profiles are consistent with the Dugdale-Barenblatt elastic-plastic fracture mechanics model In this model there is a region of inelastic deformation near the crack tip in which there is a breakdown from the yield strength of the unfractured rock to the residual frictional strength of the fault over a breakdown length S and displacement d0 Limited data also indicate that S and d0 also scale linearly with L, which implies that fracture energy G increases linearly with L The scaling parameters in these relationships depend on rock properties and are therefore not universal In our prime field locality, the Volcanic Tableland of eastern California, we have collected data over 2 orders of magnitude in scale range that show that faults obey a power law size distribution in which the exponent C in the cumulative distribution is ∼13 If the fault is growing within the brittle field, the zone of inelastic deformation consists of a brittle process zone which leaves a wake of fractured rock adjacent to the fault Preliminary results of modeling the process zone are consistent with observations now in hand both in predicting the preferred orientation of cracks in the process zone wake and the rate of falloff of crack density as a function of distance from the fault The preferred orientation of these cracks may be used to infer the mode and direction of propagation of the fault tip past the point in question According to the model, the width of the process zone wake may be used to infer the length of the fault at the time its tip passed the measurement point, but data have not yet been collected to verify this prediction If the fault displacement has been accumulated by repeated seismic slips, each of these will sweep the fault with a crack tip stress field of a smaller spatial extent than that of the fault tip stress field, producing an inner, more intensely fractured, process zone wake This may be the mechanism that creates the cataclasite zone, rather than simple frictional wear, as has been previously supposed

Patience A. Cowie

Papers

Scaling of fracture systems in geological media

Physical explanation for the displacement-length relationship of faults using a post-yield fracture mechanics model

Displacement-length scaling relationship for faults: data synthesis and discussion

Fault growth and fault scaling laws: Preliminary results

Damage zone and slip-surface evolution over μm to km scales in high-porosity Navajo sandstone, Utah