A digital bathymetric map of the oceans with a horizontal resolution of 1 to 12 kilometers was derived by combining available depth soundings with high-resolution marine gravity information from the Geosat and ERS-1 spacecraft. Previous global bathymetric maps lacked features such as the 1600-kilometer-long Foundation Seamounts chain in the South Pacific. This map shows relations among the distributions of depth, sea floor area, and sea floor age that do not fit the predictions of deterministic models of subsidence due to lithosphere cooling but may be explained by a stochastic model in which randomly distributed reheating events warm the lithosphere and raise the ocean floor.

Global Sea Floor Topography from Satellite Altimetry and Ship Depth Soundings

http://geo.mff.cuni.cz/seismosoft/clanky/2002_EPSL_Stampfli_Borel.pdf

A plate tectonic model for the Paleozoic and Mesozoic constrained by dynamic plate boundaries and restored synthetic oceanic isochrons

Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations

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

[1] A global set of present plate boundaries on the Earth is presented in digital form. Most come from sources in the literature. A few boundaries are newly interpreted from topography, volcanism, and/or seismicity, taking into account relative plate velocities from magnetic anomalies, moment tensor solutions, and/or geodesy. In addition to the 14 large plates whose motion was described by the NUVEL-1A poles (Africa, Antarctica, Arabia, Australia, Caribbean, Cocos, Eurasia, India, Juan de Fuca, Nazca, North America, Pacific, Philippine Sea, South America), model PB2002 includes 38 small plates (Okhotsk, Amur, Yangtze, Okinawa, Sunda, Burma, Molucca Sea, Banda Sea, Timor, Birds Head, Maoke, Caroline, Mariana, North Bismarck, Manus, South Bismarck, Solomon Sea, Woodlark, New Hebrides, Conway Reef, Balmoral Reef, Futuna, Niuafo'ou, Tonga, Kermadec, Rivera, Galapagos, Easter, Juan Fernandez, Panama, North Andes, Altiplano, Shetland, Scotia, Sandwich, Aegean Sea, Anatolia, Somalia), for a total of 52 plates. No attempt is made to divide the Alps-Persia-Tibet mountain belt, the Philippine Islands, the Peruvian Andes, the Sierras Pampeanas, or the California-Nevada zone of dextral transtension into plates; instead, they are designated as “orogens” in which this plate model is not expected to be accurate. The cumulative-number/area distribution for this model follows a power law for plates with areas between 0.002 and 1 steradian. Departure from this scaling at the small-plate end suggests that future work is very likely to define more very small plates within the orogens. The model is presented in four digital files: a set of plate boundary segments; a set of plate outlines; a set of outlines of the orogens; and a table of characteristics of each digitization step along plate boundaries, including estimated relative velocity vector and classification into one of 7 types (continental convergence zone, continental transform fault, continental rift, oceanic spreading ridge, oceanic transform fault, oceanic convergent boundary, subduction zone). Total length, mean velocity, and total rate of area production/destruction are computed for each class; the global rate of area production and destruction is 0.108 m2/s, which is higher than in previous models because of the incorporation of back-arc spreading.

An updated digital model of plate boundaries

We have created a digital age grid of the ocean floor with a grid node interval of 6 arc min using a self-consistent set of global isochrons and associated plate reconstruction poles. The age at each grid node was determined by linear interpolation between adjacent isochrons in the direction of spreading. Ages for ocean floor between the oldest identified magnetic anomalies and continental crust were interpolated by estimating the ages of passive continental margin segments from geological data and published plate models. We have constructed an age grid with error estimates for each grid cell as a function of (1) the error of ocean floor ages identified from magnetic anomalies along ship tracks and the age of the corresponding grid cells in our age grid, (2) the distance of a given grid cell to the nearest magnetic anomaly identification, and (3) the gradient of the age grid: i.e., larger errors are associated with high age gradients at fracture zones or other age discontinuities. Future applications of this digital grid include studies of the thermal and elastic structure of the lithosphere, the heat loss of the Earth, ridge-push forces through time, asymmetry of spreading, and providing constraints for seismic tomography and mantle convection models.

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

Digital isochrons of the world's ocean floor

We use an updated model for global relative plate motions during the past 130 m.y. together with a compilation of bathymetry and recently published radiometric dates of major hotspot tracks to derive a plate-motion model relative to major hotspots in the Atlantic and Indian oceans. Interactive computer graphics were used to find the best fit of dated hotspot tracks on the Australian, Indian, African, and North and South American plates relative to present-day hotspots assumed fixed in the mantle. One set of rotation parameters can be found that satisfies all data constraints back to chron 34 (84 Ma) and supports little motion between the major hotspots in this hemisphere. For times between 130 and 84 Ma, the plate model is based solely on the trails of the Tristan da Cunha and Great Meteor hotspots. This approach results in a location of the Kerguelen hotspot distinct from and south of the Rajmahal Traps for this time interval. Between 115 and 105 Ma, our model locates the hotspot underneath the southern Kerguelen Plateau, which is compatible with an age estimate of this part of the plateau of 115-95 Ma. Our model suggests that the 85°E ridge between lat 10°N and the Afanasiy Nikitin seamounts may have been formed by a hotspot now located underneath the eastern Conrad rise.

Revised plate motions relative to the hotspots from combined Atlantic and Indian Ocean hotspot tracks

A model of the tectonic evolution of the eastern Indian Ocean from the Late Cretaceous to the present is constructed. Reinterpretations of seafloor magnetic anomalies are used to place constraints on spreading rates and the timing of major reorganizations. Vertical deflection profiles from 22 repeat cycles of Geosat altimeter data are used to determine the tectonic fabric associated with fracture zones. The Geosat data place constraints on paleospreading directions. Consideration is given to the implications of the model for the relative motions between India, Australia, and Antarctica.

Evolution of the eastern Indian Ocean since the Late Cretaceous: Constraints from Geosat altimetry

The motions between the Somalian, Antarctic, and Australian plates—the three plates believed to meet at the Rodrigues triple junction in the Indian Ocean—are inconsistent with the assumption that all three plates are rigid. The discrepancy is best explained if the Australian plate contains two component plates. Thus, the traditionally defined Indo-Australian plate consists of three component plates and multiple diffuse plate boundaries. The pattern of present deformation indicates that the boundaries between the three component plates are two unconnected zones accommodating divergence and a larger zone, which we interpret as three diffuse convergent plate boundaries and a diffuse triple junction.

https://science.sciencemag.org/content/sci/277/5330/1268.full.pdf

J. Y. Royer

Papers

Digital isochrons of the world's ocean floor

Digital isochrons of the world's ocean floor

Revised plate motions relative to the hotspots from combined Atlantic and Indian Ocean hotspot tracks

Evolution of the eastern Indian Ocean since the Late Cretaceous: Constraints from Geosat altimetry

The Motion and Boundary Between the Capricorn and Australian Plates