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

Recent revisions to the geomagnetic time scale indicate that global plate motion model NUVEL-1 should be modified for comparison with other rates of motion including those estimated from space geodetic measurements. The optimal recalibration, which is a compromise among slightly different calibrations appropriate for slow, medium, and fast rates of seafloor spreading, is to multiply NUVEL-1 angular velocities by a constant, α, of 0.9562. We refer to this simply recalibrated plate motion model as NUVEL-1A, and give correspondingly revised tables of angular velocities and uncertainties. Published work indicates that space geodetic rates are slower on average than those calculated from NUVEL-1 by 6±1%. This average discrepancy is reduced to less than 2% when space geodetic rates are instead compared with NUVEL-1A.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/94GL02118%4010.1002/%28ISSN%291944-8007.GRL40

Effect of recent revisions to the geomagnetic reversal time scale on estimates of 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.

Networks of dozens to hundreds of permanently operating precision Global Positioning System (GPS) receivers are emerging at spatial scales that range from 10(exp 0) to 10(exp 3) km. To keep the computational burden associated with the analysis of such data economically feasible, one approach is to first determine precise GPS satellite positions and clock corrections from a globally distributed network of GPS receivers. Their, data from the local network are analyzed by estimating receiver- specific parameters with receiver-specific data satellite parameters are held fixed at their values determined in the global solution. This "precise point positioning" allows analysis of data from hundreds to thousands of sites every (lay with 40-Mflop computers, with results comparable in quality to the simultaneous analysis of all data. The reference frames for the global and network solutions can be free of distortion imposed by erroneous fiducial constraints on any sites.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/96JB03860

Precise point positioning for the efficient and robust analysis of GPS data from large networks

SUMMARY
We describe best-fitting angular velocities and MORVEL, a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface. Seafloor spreading rates and fault azimuths are used to determine the motions of 19 plates bordered by mid-ocean ridges, including all the major plates. Six smaller plates with little or no connection to the mid-ocean ridges are linked to MORVEL with GPS station velocities and azimuthal data. By design, almost no kinematic information is exchanged between the geologically determined and geodetically constrained subsets of the global circuit—MORVEL thus averages motion over geological intervals for all the major plates. Plate geometry changes relative to NUVEL-1A include the incorporation of Nubia, Lwandle and Somalia plates for the former Africa plate, Capricorn, Australia and Macquarie plates for the former Australia plate, and Sur and South America plates for the former South America plate. MORVEL also includes Amur, Philippine Sea, Sundaland and Yangtze plates, making it more useful than NUVEL-1A for studies of deformation in Asia and the western Pacific. Seafloor spreading rates are estimated over the past 0.78 Myr for intermediate and fast spreading centres and since 3.16 Ma for slow and ultraslow spreading centres. Rates are adjusted downward by 0.6–2.6 mm yr−1 to compensate for the several kilometre width of magnetic reversal zones. Nearly all the NUVEL-1A angular velocities differ significantly from the MORVEL angular velocities. The many new data, revised plate geometries, and correction for outward displacement thus significantly modify our knowledge of geologically current plate motions. MORVEL indicates significantly slower 0.78-Myr-average motion across the Nazca–Antarctic and Nazca–Pacific boundaries than does NUVEL-1A, consistent with a progressive slowdown in the eastward component of Nazca plate motion since 3.16 Ma. It also indicates that motions across the Caribbean–North America and Caribbean–South America plate boundaries are twice as fast as given by NUVEL-1A. Summed, least-squares differences between angular velocities estimated from GPS and those for MORVEL, NUVEL-1 and NUVEL-1A are, respectively, 260 per cent larger for NUVEL-1 and 50 per cent larger for NUVEL-1A than for MORVEL, suggesting that MORVEL more accurately describes historically current plate motions. Significant differences between geological and GPS estimates of Nazca plate motion and Arabia–Eurasia and India–Eurasia motion are reduced but not eliminated when using MORVEL instead of NUVEL-1A, possibly indicating that changes have occurred in those plate motions since 3.16 Ma. The MORVEL and GPS estimates of Pacific–North America plate motion in western North America differ by only 2.6 ± 1.7 mm yr−1, ≈25 per cent smaller than for NUVEL-1A. The remaining difference for this plate pair, assuming there are no unrecognized systematic errors and no measurable change in Pacific–North America motion over the past 1–3 Myr, indicates deformation of one or more plates in the global circuit. Tests for closure of six three-plate circuits indicate that two, Pacific–Cocos–Nazca and Sur–Nubia–Antarctic, fail closure, with respective linear velocities of non-closure of 14 ± 5 and 3 ± 1 mm yr−1 (95 per cent confidence limits) at their triple junctions. We conclude that the rigid plate approximation continues to be tremendously useful, but—absent any unrecognized systematic errors—the plates deform measurably, possibly by thermal contraction and wide plate boundaries with deformation rates near or beneath the level of noise in plate kinematic data.

Geologically current plate motions

A new model of the last deglaciation event of the Late Quaternary ice age is here described and denoted as ICE-6G_C (VM5a). It differs from previously published models in this sequence in that it has been explicitly refined by applying all of the available Global Positioning System (GPS) measurements of vertical motion of the crust that may be brought to bear to constrain the thickness of local ice cover as well as the timing of its removal. Additional space geodetic constraints have also been applied to specify the reference frame within which the GPS data are described. The focus of the paper is upon the three main regions of Last Glacial Maximum ice cover, namely, North America, Northwestern Europe/Eurasia, and Antarctica, although Greenland and the British Isles will also be included, if peripherally, in the discussion. In each of the three major regions, the model predictions of the time rate of change of the gravitational field are also compared to that being measured by the Gravity Recovery and Climate Experiment satellites as an independent means of verifying the improvement of the model achieved by applying the GPS constraints. Several aspects of the global characteristics of this new model are also discussed, including the nature of relative sea level history predictions at far-field locations, in particular the Caribbean island of Barbados, from which especially high-quality records of postglacial sea level change are available but which records were not employed in the development of the model. Although ICE-6G_C (VM5a) is a significant improvement insofar as the most recently available GPS observations are concerned, comparison of model predictions with such far-field relative sea level histories enables us to identify a series of additional improvements that should follow from a further stage of model iteration.

Donald F. Argus

Papers

Current plate motions

Effect of recent revisions to the geomagnetic reversal time scale on estimates of current plate motions

Current plate motions

Geologically current plate motions

Space geodesy constrains ice age terminal deglaciation: The global ICE‐6G_C (VM5a) model