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.

[1] We present a new global model for Recent plate velocities, REVEL, describing the relative velocities of 19 plates and continental blocks. The model is derived from publicly available space geodetic (primarily GPS) data for the period 1993–2000. We include an independent and rigorous estimate for GPS velocity uncertainties to assess plate rigidity and propagate these uncertainties to the velocity estimates. The velocity fields for North America, Eurasia, and Antarctica clearly show the effects of glacial isostatic adjustment, and Australia appears to depart from rigid plate behavior in a manner consistent with the mapped intraplate stress field. Two thirds of tested plate pairs agree with the NUVEL-1A geologic (3 Myr average) velocities within uncertainties. Three plate pairs (Caribbean–North America, Caribbean–South America, and North America–Pacific) exhibit significant differences between the geodetic and geologic model that may reflect systematic errors in NUVEL-1A due to the use of seafloor magnetic rate data that do not reflect the full plate rate because of tectonic complexities. Most other differences probably reflect real velocity changes over the last few million years. Several plate pairs (Arabia–Eurasia, Arabia–Nubia, Eurasia–India) move more slowly than the 3 Myr NUVEL-1A average, perhaps reflecting long-term deceleration associated with continental collision. Several other plate pairs, including Nazca–Pacific, Nazca–South America and Nubia–South America, are experiencing slowing that began ∼25 Ma, the beginning of the current phase of Andean crustal shortening.

/pdf/revel-a-model-for-recent-plate-velocities-from-space-geodesy-18oq1lu9ey.pdf

REVEL: A model for Recent plate velocities from space geodesy

SUMMARY In this paper we present a global model (GSRM-1) of both horizontal velocities on the Earth’s surface and horizontal strain rates for almost all deforming plate boundary zones. A model strain rate field is obtained jointly with a global velocity field in the process of solving for a global velocity gradient tensor field. In our model we perform a least-squares fit between model velocities and observed geodetic velocities, as well as between model strain rates and observed geological strain rates. Model velocities and strain rates are interpolated over a spherical Earth using bi-cubic Bessel splines. We include 3000 geodetic velocities from 50 different, mostly published, studies. Geological strain rates are obtained for central Asia only and they are inferred from Quaternary fault slip rates. For all areas where no geological information is included a priori constraints are placed on the style and direction (but not magnitude) of the model strain rate field. These constraints are taken from a seismic strain rate field inferred from (normalized) focal mechanisms of shallow earthquakes. We present a global solution of the second invariant of the model strain rate field as well as strain rate solutions for a few selected plate boundary zones. Generally, the strain rate tensor field is consistent with geological and seismological data. With the assumption of plate rigidity for all areas other than the plate boundary zones we also present relative angular velocities for most plate pairs. We find that in general there is a good agreement between the present-day plate motions we obtain and longterm plate motions, but a few significant differences exist. The rotation rates for the Indian, Arabian and Nubian plates relative to Eurasia are 30, 13 and 50 per cent slower than the NUVEL1A estimate, respectively, and the rotation rate for the Nazca Plate relative to South America is 17 per cent slower. On the other hand, Caribbean‐North America motion is 76 per cent faster than the NUVEL-1A estimate. While crustal blocks in the India‐Eurasia collision zone move significantly and self-consistently with respect to bounding plates, only a very small motion is predicted between the Nubian and Somalian plates. By integrating plate boundary zone deformation with the traditional modelling of angular velocities of rigid plates we have obtained a model that has already been proven valuable in, for instance, redefining a no-net-rotation model of surface motions and by confirming a global correlation between seismicity rates and tectonic moment rates along subduction zones and within zones of continental deformation.

/pdf/an-integrated-global-model-of-present-day-plate-motions-and-57ic0m6oev.pdf

An integrated global model of present-day plate motions and plate boundary deformation

A significant discrepancy exists between the heat flow measured at the seafloor and the higher values predicted by thermal models of the cooling lithosphere. This discrepancy is generally interpreted as indicating that the upper oceanic crust is cooled significantly by hydrothermal circulation. The magnitude of this heat flow discrepancy is the primary datum used to estimate the volume of hydrothermal flow, and the variation in the discrepancy with lithospheric age is the primary constraint on how the hydrothermal flux is divided between near-ridge and off-ridge environments. The resulting estimates are important for investigation of both the thermal structure of the lithosphere and the chemistry of the oceans. We reevaluate the magnitude and age variation of the discrepancy using a global heat flow data set substantially larger than in earlier studies, and the GDHI (Global Depth and Heat Flow) model that better predicts the heat flow. We estimate that of the predicted global oceanic heat flux of 32 x 10(exp 12) W, 34% (11 x 10(exp 12) W) occurs by hydrothermal flow. Approximately 30% of the hydrothermal heat flux occurs in crust younger than 1 Ma, so the majority of this flux is off-ridge. These hydrothermal heat flux estimates are upper bounds, because heat flow measurements require sediment at the site and so are made preferentially at topographic lows, where heat flow may be depressed. Because the water temperature for the near-ridge flow exceeds that for the off-ridge flow, the near-ridge water flow will be even a smaller fraction of the total water flow. As a result, in estimating fluxes from geochemical data, use of the high water temperatures appropriate for the ridge axis may significantly overestimate the heat flux for an assumed water flux or underestimate the water flux for an assumed heat flux. Our data also permit improved estimates of the 'sealing' age, defined as the age where the observed heat flow approximately equals that predicted, suggesting that hydrothermal heat transfer has largely ceased. Although earlier studies suggested major differences in sealing ages for different ocean basins, we find that the sealing ages for the Atlantic, Pacific, and Indian oceans are similar and consistent with the sealing age for the entire data set, 65 +/- 10 Ma. The previous inference of a young (approximately 20 Ma) sealing age for the Pacific appears to have biased downward several previous estimates of the global hydrothermal flux. The heat flow data also provide indirect evidence for the mechanism by which the hydrothermal heat flux becomes small, which has often been ascribed to isolation of the igneous crust from seawater due to the hydraulic conductivity of the intervening sediment. We find, however, that even the least sedimented sites show the systematic increase of the ratio of observed to predicted heat flow with age, although the more sedimented sites have a younger sealing age. Moreover, the heat flow discrepancy persists at heavily sedimented sites until approximately 50 Ma. It thus appears that approximately 100-200 m of sediment is neither necessary nor sufficient to stop hydrothermal heat transfer. We therefore conclude that the age of the crust is the primary control on the fraction of heat transported by hydrothermal flow and that sediment thickness has a lesser effect. This inference is consistent with models in which hydrothermal flow decreases with age due to reduced crustal porosity and hence permeability.

Constraints on hydrothermal heat flux through the oceanic lithosphere from global heat flow

It is suggested that motion along the virtually aseismic Owen fracture zone is negligible, so that Arabia and India are contained within a single Indo-Arabian plate divided from the Australian plate by a diffuse boundary. The boundary is a zone of concentrated seismicity and deformation commonly characterized as 'intraplate'. The rotation vector of Australia relative to Indo-Arabia is consistent with the seismologically observed 2 cm/yr of left-lateral strike-slip along the Ninetyeast Ridge, north-south compression in the Central Indian Ocean, and the north-south extension near Chagos.

Stuart Weinstein

Papers

A diffuse plate boundary model for Indian Ocean tectonics