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

Physical factors likely to affect the genesis of the various fault rocks—frictional properties, temperature, effective stress normal to the fault and differential stress—are examined in relation to the energy budget of fault zones, the main velocity modes of faulting and the type of faulting, whether thrust, wrench, or normal. In a conceptual model of a major fault zone cutting crystalline quartzo-feldspathic crust, a zone of elastico-frictional (EF) behaviour generating random-fabric fault rocks (gouge—breccia—cataclasite series—pseudotachylyte) overlies a region where quasi-plastic (QP) processes of rock deformation operate in ductile shear zones with the production of mylonite series rocks possessing strong tectonite fabrics. In some cases, fault rocks developed by transient seismic faulting can be distinguished from those generated by slow aseismic shear. Random-fabric fault rocks may form as a result of seismic faulting within the ductile shear zones from time to time, but tend to be obliterated by continued shearing. Resistance to shear within the fault zone reaches a peak value (greatest for thrusts and least for normal faults) around the EF/QP transition level, which for normal geothermal gradients and an adequate supply of water, occurs at depths of 10–15 km.

Fault rocks and fault mechanisms

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

Convective removal of lower lithosphere beneath the Tibetan Plateau can account for a rapid increase in the mean elevation of the Tibetan Plateau of 1000 m or more in a few million years. Such uplift seems to be required by abrupt tectonic and environmental changes in Asia and the Indian Ocean in late Cenozoic time. The composition of basaltic volcanism in northern Tibet, which apparently began at about 13 Ma, implies melting of lithosphere, not asthenosphere. The most plausible mechanism for rapid heat transfer to the midlithosphere is by convective removal of deeper lithosphere and its replacement by hotter asthenosphere. The initiation of normal faulting in Tibet at about 8 (± 3) Ma suggests that the plateau underwent an appreciable increase in elevation at that time. An increase due solely to the isostatic response to crustal thickening caused by India's penetration into Eurasia should have been slow and could not have triggered normal faulting. Another process, such as removal of relatively cold, dense lower lithosphere, must have caused a supplemental uplift of the surface. Folding and faulting of the Indo-Australian plate south of India, the most prominent oceanic intraplate deformation on Earth, began between about 7.5 and 8 Ma and indicates an increased north-south compressional stress within the Indo-Australian plate. A Tibetan uplift of only 1000 m, if the result of removal of lower lithosphere, should have increased the compressional stress that the plateau applies to India and that resists India's northward movement, from an amount too small to fold oceanic lithosphere, to one sufficient to do so. The climate of the equatorial Indian Ocean and southern Asia changed at about 6–9 Ma: monsoonal winds apparently strengthened, northern Pakistan became more arid, but weathering of rock in the eastern Himalaya apparently increased. Because of its high altitude and lateral extent, the Tibetan Plateau provides a heat source at midlatitudes that should oppose classical (symmetric) Hadley circulation between the equator and temperate latitudes and that should help to drive an essentially opposite circulation characteristic of summer monsoons. For the simple case of axisymmetric heating (no dependence on longitude) of an atmosphere without dissipation, theoretical analyses by Hou, Lindzen, and Plumb show that an axisymmetric heat source displaced from the equator can drive a much stronger meridianal (monsoonlike) circulation than such a source centered on the equator, but only if heating exceeds a threshold whose level increases with the latitude of the heat source. Because heating of the atmosphere over Tibet should increase monotonically with elevation of the plateau, a modest uplift (1000–2500 m) of Tibet, already of substantial extent and height, might have been sufficient to exceed a threshold necessary for a strong monsoon. The virtual simultaneity of these phenomena suggests that uplift was rapid: approximately 1000 m to 2500 m in a few million years. Moreover, nearly simultaneously with the late Miocene strengthening of the monsoon, the calcite compensation depth in the oceans dropped, plants using the relatively efficient C4 pathway for photosynthesis evolved rapidly, and atmospheric CO2 seems to have decreased, suggesting causal relationships and positive feedbacks among these phenomena. Both a supplemental uplift of the Himalaya, the southern edge of Tibet, and a strengthened monsoon may have accelerated erosion and weathering of silicate rock in the Himalaya that, in turn, enhanced extraction of CO2 from the atmosphere. Thus these correlations offer some support for links between plateau uplift, a downdrawing of CO2 from the atmosphere, and global climate change, as proposed by Raymo, Ruddiman, and Froehlich. Mantle dynamics beneath mountain belts not only may profoundly affect tectonic processes near and far from the belts, but might also play an important role in altering regional and global climates.

Mantle dynamics, uplift of the Tibetan Plateau, and the Indian Monsoon

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

Many mesothermal gold-quartz deposits are localized along high-angle reverse or reverse-oblique shear zones within greenstone belt terrains. Characteristically, these fault-hosted vein deposits exhibit a mixed "brittle-ductile" style of deformation (discrete shears and vein fractures as well as a schistose shear-zone fabric) developed under greenschist facies metamorphic conditions. Many of the vein systems are of considerable vertical extent (>2 km); they include steeply dipping fault veins (lenticular veins subparallel to the shear-zone schistosity) and, in some cases, associated flats (subhorizontal extensional veins). Textures of both vein sets record histories of incremental deposition. We infer that the vein sets developed near the roofs of active metamorphic/magmatic systems and represent the roots of brittle, high-angle reverse fault systems extending upward through the seismogenic regime. Friction theory and field relations suggest that the high-angle reverse faults acted as valves , promoting cyclic fluctuations in fluid pressure from supralithostatic to hydrostatic values. Because of their unfavorable orientation in the prevailing stress field, reactivation of the faults could only occur when fluid pressure exceeded the lithostatic load. Seismogenic fault failure then created fracture permeability within the rupture zone, allowing sudden draining of the geopressured reservoir at depth. Incremental opening of flats is attributed to the prefailure stage of supralithostatic fluid pressures; deposition within fault veins is attributed to the immediate postfailure discharge phase. Hydrothermal self-sealing leads to reaccumulation of fluid pressure and a repetition of the cycle. Mutual crosscutting relations between the two vein sets are a natural consequence of the cyclicity of the process. Abrupt fluid-pressure fluctuations from this fault-valve behavior of reverse faults seem likely to be integral to the mineralizing process at this structural level.

High-angle reverse faults, fluid-pressure cycling, and mesothermal gold-quartz deposits

Structural permeability of fluid-driven fault-fracture meshes

Models of fault zones in continental crust, based on the analysis of rock deformation textures, suggest that the depth of seismic activity is controlled by the passage from a pressure-sensitive, dominantly frictional regime to strongly temperature-dependent, quasi-plastic mylonitization at greenschist and higher grades of metamorphism. Sufficient knowledge now exists concerning the frictional and rheological properties of quartz-bearing rocks to construct crude strength-depth curves for different geotherms. In such models, shear resistance peaks sharply at the inferred seismic-aseismic transition. The maximum depth of microseismic activity in various heat flow provinces of the conterminous United States generally correlates well with the frictional to quasi-plastic transition modeled for the different geotherms. Larger earthquakes ( M L > 5.5) also tend to nucleate near the base of the seismogenic zone. This region is postulated to have the highest concentration of distortional strain energy for stress levels at failure, and can be regarded as the prime asperity in crustal fault zones.

Fault zone models, heat flow, and the depth distribution of earthquakes in the continental crust of the United States

Summary
Pseudotachylyte occurs as vein material infilling highly brittle shear and extensional fractures developed along the western margin of the late Caledonian, Outer Hebrides Thrust zone in NW Scotland. Vein geometries and textures show clearly that the pseudotachylyte has been through a melt phase. From the composition of the pseudotachylyte matrix which is close to that of a basaltic andesite, probable melt temperatures of around 1100°C are inferred. Field and theoretical studies demonstrate that the pseudotachylyte was generated by relatively high stress seismic faulting in crystalline sialic crust devoid of an intergranular fluid, most probably at the time of thrust inception and at a depth of around 4–5 km. A study of pseudotachylyte-bearing ‘single-jerk’ microfaults shows that the slip (d) is related to the thickness of the pseudotachylyte layer (a) by the equation, 



where d and a are measured in centimetres. Work-energy calculations based on this empirical relationship suggest that the pre-failure shear stress on the microfaults must have been as high as 1.6 kbar to overcome the initial frictional resistance (τf), which decreases with increasing slip during a single movement according to the relationship, 



which may arise solely from the viscous shear resistance of the melt layer.

Delineation of palaeoseismic zones by the recognition of those cataclastic rocks which are necessarily the products of earthquake faulting, may assist in the determination of ancient plate boundaries.

/pdf/generation-of-pseudotachylyte-by-ancient-seismic-faulting-26w4kjehmz.pdf

Richard H. Sibson

Papers

Fault rocks and fault mechanisms

High-angle reverse faults, fluid-pressure cycling, and mesothermal gold-quartz deposits

Structural permeability of fluid-driven fault-fracture meshes

Fault zone models, heat flow, and the depth distribution of earthquakes in the continental crust of the United States

Generation of Pseudotachylyte by Ancient Seismic Faulting