Assessing plate reconstruction models using plate driving force consistency tests

Plate reconstruction models are constructed to fit constraints such as magnetic anomalies, fracture zones, paleomagnetic poles, geological observations and seismic tomography. However, these models do not consider the physical equations of plate driving forces when reconstructing plate motion. This can potentially result in geodynamically-implausible plate motions, which has implications for a range of work based on plate reconstruction models. We present a new algorithm that calculates time-dependent slab pull, ridge push (GPE force) and mantle drag resistance for any topologically closed reconstruction, and evaluates the residuals-or missing components-required for torques to balance given our assumed plate driving force relationships. In all analyzed models, residual torques for the present-day are three orders of magnitude smaller than the typical driving torques for oceanic plates, but can be of the same order of magnitude back in time-particularly from 90 to 50 Ma. Using the Pacific plate as an example, we show how our algorithm can be used to identify areas and times with high residual torques, where either plate reconstructions have a high degree of geodynamic implausibility or our understanding of the underlying geodynamic forces is incomplete. We suggest strategies for plate model improvements and also identify times when other forces such as active mantle flow were likely important contributors. Our algorithm is intended as a tool to help assess and improve plate reconstruction models based on a transparent and expandable set of a priori dynamic constraints.

Southern Ocean biogenic blooms freezing-in Oligocene colder climates

Crossing a key atmospheric CO₂ threshold triggered a fundamental global climate reorganisation ~34 million years ago (Ma) establishing permanent Antarctic ice sheets. Curiously, a more dramatic CO₂ decline (~800-400 ppm by the Early Oligocene(~27 Ma)), postdates initial ice sheet expansion but the mechanisms driving this later, rapid drop in atmospheric carbon during the early Oligocene remains elusive and controversial. Here we use marine seismic reflection and borehole data to reveal an unprecedented accumulation of early Oligocene strata (up to 2.2 km thick over 1500 × 500 km) with a major biogenic component in the Australian Southern Ocean. High-resolution ocean simulations demonstrate that a tectonically-driven, one-off reorganisation of ocean currents, caused a unique period where current instability coincided with high nutrient input from the Antarctic continent. This unrepeated and short-lived environment favoured extreme bioproductivity and enhanced sediment burial. The size and rapid accumulation of this sediment package potentially holds ~1.067 × 10¹⁵ kg of the 'missing carbon' sequestered during the decline from an Eocene high CO₂-world to a mid-Oligocene medium CO₂-world, highlighting the exceptional role of the Southern Ocean in modulating long-term climate.

Rheology of Naturally Deformed Antigorite Serpentinite: Strain and Strain-Rate Dependence at Mantle-Wedge Conditions

Antigorite serpentinite is expected to occur in parts of subduction plate boundaries, and may suppress earthquake slip, but the dominant deformation mechanisms and resultant rheology of antigorite are unclear. An exhumed plate boundary shear zone exposed near Nagasaki, Japan, contains antigorite deformed at 474°C ± 30°C. Observations indicate that a foliation defined by (001) crystal facets developed during plate-boundary shear. Microstructures indicating grain-scale dissolution at high-stress interfaces and precipitation in low-stress regions suggest that dissolution-precipitation creep contributed to foliation development. Analysis of crystal orientations indicate a small contribution from dislocation activity. We suggest a frictional-viscous rheology for antigorite, where dissolution-precipitation produces a foliation defined by (001) crystal facets and acts to resolve strain incompatibilities, allowing for efficient face-to-face sliding between facets. This rheology can not only explain aseismic behavior at ambient plate boundary conditions, but also some of the contrasting behaviors shown by previous field and laboratory studies.

Hydrous oceanic crust hosts megathrust creep at low shear stresses

The rheology of the metamorphosed oceanic crust may be a critical control on megathrust strength and deformation style. However, little is known about the strength and deformation style of metamorphosed basalt. Exhumed megathrust shear zones exposed on Kyushu, SW Japan, contain hydrous metabasalts deformed at temperatures between ~300° and ~500°C, spanning the inferred temperature-controlled seismic-aseismic transition. Field and microstructural observations of these shear zones, combined with quartz grain-size piezometry, indicate that metabasalts creep at shear stresses <100 MPa at ~370°C and at shear stresses <30 MPa at ~500°C. These values are much lower than those suggested by viscous flow laws for basalt. The implication is that relatively weak, hydrous, metamorphosed oceanic crust can creep at low viscosities over a wide shear zone and have a critical influence on plate interface strength and deformation style around the seismic-aseismic transition.

Rapid subduction initiation and magmatism in the Western Pacific driven by internal vertical forces

Plate tectonics requires the formation of plate boundaries. Particularly important is the enigmatic initiation of subduction: the sliding of one plate below the other, and the primary driver of plate tectonics. A continuous, in situ record of subduction initiation was recovered by the International Ocean Discovery Program Expedition 352, which drilled a segment of the fore-arc of the Izu-Bonin-Mariana subduction system, revealing a distinct magmatic progression with a rapid timescale (approximately 1 million years). Here, using numerical models, we demonstrate that these observations cannot be produced by previously proposed horizontal external forcing. Instead a geodynamic evolution that is dominated by internal, vertical forces produces both the temporal and spatial distribution of magmatic products, and progresses to self-sustained subduction. Such a primarily internally driven initiation event is necessarily whole-plate scale and the rock sequence generated (also found along the Tethyan margin) may be considered as a smoking gun for this type of event.

The magmatic progression produced during the initiation of the Izu-Bonin-Marianas subduction zone took place rapidly over 1 million years, but it has been unclear why. Here, using numerical models, the authors show that subduction initiation was dominated by vertical forces, internal to the system itself, progressing to self-sustained subduction.

Postmagmatic Tectonic Evolution of the Outer Izu-Bonin Forearc Revealed by Sediment Basin Structure and Vein Microstructure Analysis: Implications for a 15 Ma Hiatus Between Pacific Plate Subduction Initiation and Forearc Extension

International Ocean Discovery Program Expedition 352 recovered sedimentary-volcaniclastic successions and extensional structures (faults and extensional veins) that allow the reconstruction of the Izu-Bonin forearc tectonic evolution using a combination of shipboard core data, seismic reflection images, and calcite vein microstructure analysis. The oldest recorded biostratigraphic ages within fault-bounded sedimentary basins (Late Eocene to Early Oligocene) imply a ~15 Ma hiatus between the formation of the igneous basement (52 to 50 Ma) and the onset of sedimentation. At the upslope sites (U1439 and U1442) extension led to the formation of asymmetric basins reflecting regional stretch of ~16-19% at strain rates of ~1.58 × 10-16 to 4.62 × 10-16 s-1. Downslope Site U1440 (closer to the trench) is characterized by a symmetric graben bounded by conjugate normal faults reflecting regional stretch of ~55% at strain rates of 4.40 × 10-16 to 1.43 × 10-15 s-1. Mean differential stresses are in the range of ~70-90 MPa. We infer that upper plate extension was triggered by incipient Pacific Plate rollback ~15 Ma after subduction initiation. Extension was accommodated by normal faulting with syntectonic sedimentation during Late Eocene to Early Oligocene times. Backarc extension was assisted by magmatism with related Shikoku and Parece-Vela Basin spreading at ~25 Ma, so that parts of the arc and rear arc, and the West Philippine backarc Basin were dismembered from the forearc. This was followed by slow-rift to postrift sedimentation during the transition from forearc to arc rifting to spreading within the Shikoku-Parece-Vela Basin system.

Magmatic Response to Subduction Initiation: Part 1. Fore-arc Basalts of the Izu-Bonin Arc From IODP Expedition 352

The Izu-Bonin-Mariana (IBM) fore arc preserves igneous rock assemblages that formed during subduction initiation circa 52 Ma. International Ocean Discovery Program (IODP) Expedition 352 cored four sites in the fore arc near the Ogasawara Plateau in order to document the magmatic response to subduction initiation and the physical, petrologic, and chemical stratigraphy of a nascent subduction zone. Two of these sites (U1440 and U1441) are underlain by fore-arc basalt (FAB). FABs have mid-ocean ridge basalt (MORB)-like compositions, however, FAB are consistently lower in the high-field strength elements (TiO2, P2O5, Zr) and Ni compared to MORB, with Na2O at the low end of the MORB field and FeO* at the high end. Almost all FABs are light rare earth element depleted, with low total REE, and have low ratios of highly incompatible to less incompatible elements (Ti/V, Zr/Y, Ce/Yb, and Zr/Sm) relative to MORB. Chemostratigraphic trends in Hole U1440B are consistent with the uppermost lavas forming off axis, whereas the lower lavas formed beneath a spreading center axis. Axial magma of U1440B becomes more fractionated upsection; overlying off-axis magmas return to more primitive compositions. Melt models require a two-stage process, with early garnet field melts extracted prior to later spinel field melts, with up to 23% melting to form the most depleted compositions. Mantle equilibration temperatures are higher than normal MORB (1,400 °C-1,480 °C) at relatively low pressures (1-2 GPa), which may reflect an influence of the Manus plume during subduction initiation. Our data support previous models of FAB origin by decompression melting but imply a source more depleted than normal MORB source mantle.

Cenozoic Rotation History of Borneo and Sundaland, SE Asia Revealed by Paleomagnetism, Seismic Tomography, and Kinematic Reconstruction

SE Asia comprises a heterogeneous assemblage of fragments derived from Cathaysia (Eurasia) in the north and Gondwana in the south, separated by suture zones representing closed former ocean basins. The western part of the region comprises Sundaland, which was formed by Late Permian-Triassic amalgamation of continental and arc fragments now found in Indochina, the Thai Penisula, Peninsular Malaysia, and Sumatra. On Borneo, the Kuching Zone formed the eastern margin of Sundaland since the Triassic. To the SE of the Kuching Zone, the Gondwana-derived continental fragments of SW Borneo and East Kalimantan accreted in the Cretaceous. South China-derived fragments accreted to north of the Kuching Zone in the Miocene. Deciphering this complex geodynamic history of SE Asia requires restoration of its deformation history, but quantitative constraints are often sparse. Paleomagnetism may provide such constraints. Previous paleomagnetic studies demonstrated that Sundaland and fragments in Borneo underwent vertical axis rotations since the Cretaceous. We provide new paleomagnetic data from Eocene-Miocene sedimentary rocks in the Kutai Basin, east Borneo, and critically reevaluate the published database, omitting sites that do not pass widely used, up-to-date reliability criteria. We use the resulting database to develop an updated kinematic restoration. We test the regional or local nature of paleomagnetic rotations against fits between the restored orientation of the Sunda Trench and seismic tomography images of the associated slabs. Paleomagnetic data and mantle tomography of the Sunda slab indicate that Sundaland did not experience significant vertical axis rotations since the Late Jurassic. Paleomagnetic data show that Borneo underwent a ~35° counterclockwise rotation constrained to the Late Eocene and an additional ~10° counterclockwise rotation since the Early Miocene. How this rotation was accommodated relative to Sundaland is enigmatic but likely involved distributed extension in the West Java Sea between Borneo and Sumatra. This Late Eocene-Early Oligocene rotation is contemporaneous with and may have been driven by a marked change in motion of Australia relative to Eurasia, from eastward to northward, which also has led to the initiation of subduction along the eastern Sunda trench and the proto-South China Sea to the south and north of Borneo, respectively.

Southwest Pacific Absolute Plate Kinematic Reconstruction Reveals Major Cenozoic Tonga-Kermadec Slab Dragging

Tectonic plates subducting at trenches having strikes oblique to the absolute subducting plate motion undergo trench-parallel slab motion through the mantle, recently defined as a form of "slab dragging." We investigate here long-term slab-dragging components of the Tonga-Kermadec subduction system driven by absolute Pacific plate motion. To this end we develop a kinematic restoration of Tonga-Kermadec Trench motion placed in a mantle reference frame and compare it to tomographically imaged slabs in the mantle. Estimating Tonga-Kermadec subduction initiation is challenging because another (New Caledonia) subduction zone existed during the Paleogene between the Australia and Pacific plates. We test partitioning of plate convergence across the Paleogene New Caledonia and Tonga-Kermadec subduction zones against resulting mantle structure and show that most, if not all, Tonga-Kermadec subduction occurred after ca. 30 Ma. Since then, Tonga-Kermadec subduction has accommodated 1,700 to 3,500 km of subduction along the southern and northern ends of the trench, respectively. When placed in a mantle reference frame, the predominantly westward directed subduction evolved while the Tonga-Kermadec Trench underwent ~1,200 km of northward absolute motion. We infer that the entire Tonga-Kermadec slab was laterally transported through the mantle over 1,200 km. Such slab dragging by the Pacific plate may explain observed deep-slab deformation and may also have significant effects on surface tectonics, both resulting from the resistance to slab dragging by the viscous mantle.

Intraoceanic subduction spanned the Pacific in the Late Cretaceous-Paleocene

The notorious ~60° bend separating the Hawaiian and Emperor chains marked a prominent change in the motion of the Pacific plate at ~47 Ma (million years ago), but the origin of that change remains an outstanding controversy that bears on the nature of major plate reorganizations. Lesser known but equally significant is a conundrum posed by the pre-bend (~80 to 47 Ma) motion of the Pacific plate, which, according to conventional plate models, was directed toward a fast-spreading ridge, in contradiction to tectonic forcing expectations. Using constraints provided by seismic tomography, paleomagnetism, and continental margin geology, we demonstrate that two intraoceanic subduction zones spanned the width of the North Pacific Ocean in Late Cretaceous through Paleocene time, and we present a simple plate tectonic model that explains how those intraoceanic subduction zones shaped the ~80 to 47 Ma kinematic history of the Pacific realm and drove a major plate reorganization.

Kinematics and dynamics of the East Pacific Rise linked to a stable, deep-mantle upwelling

Earth's tectonic plates are generally considered to be driven largely by negative buoyancy associated with subduction of oceanic lithosphere. In this context, mid-ocean ridges (MORs) are passive plate boundaries whose divergence accommodates flow driven by subduction of oceanic slabs at trenches. We show that over the past 80 million years (My), the East Pacific Rise (EPR), Earth's dominant MOR, has been characterized by limited ridge-perpendicular migration and persistent, asymmetric ridge accretion that are anomalous relative to other MORs. We reconstruct the subduction-related buoyancy fluxes of plates on either side of the EPR. The general expectation is that greater slab pull should correlate with faster plate motion and faster spreading at the EPR. Moreover, asymmetry in slab pull on either side of the EPR should correlate with either ridge migration or enhanced plate velocity in the direction of greater slab pull. Based on our analysis, none of the expected correlations are evident. This implies that other forces significantly contribute to EPR behavior. We explain these observations using mantle flow calculations based on globally integrated buoyancy distributions that require core-mantle boundary heat flux of up to 20 TW. The time-dependent mantle flow predictions yield a long-lived deep-seated upwelling that has its highest radial velocity under the EPR and is inferred to control its observed kinematics. The mantle-wide upwelling beneath the EPR drives horizontal components of asthenospheric flows beneath the plates that are similarly asymmetric but faster than the overlying surface plates, thereby contributing to plate motions through viscous tractions in the Pacific region.

Indian and African plate motions driven by the push force of the Réunion plume head

Mantle plumes are thought to play an important part in the Earth's tectonics, yet it has been difficult to isolate the effect that plumes have on plate motions. Here we analyse the plate motions involved in two apparently disparate events--the unusually rapid motion of India between 67 and 52 million years ago and a contemporaneous, transitory slowing of Africa's motion--and show that the events are coupled, with the common element being the position of the Indian and African plates relative to the location of the Réunion plume head. The synchroneity of these events suggests that they were both driven by the force of the Réunion plume head. The recognition of this plume force has substantial tectonic implications: the speed-up and slowdown of India, the possible cessation of convergence between Africa and Eurasia in the Palaeocene epoch and the enigmatic bends of the fracture zones on the Southwest Indian Ridge can all be attributed to the Réunion plume.

Response to Comment on "Major Australian-Antarctic Plate Reorganization at Hawaiian-Emperor Bend Time"

Accurately locating boundaries between continental and oceanic crust is topical in view of locating offshore boundaries relevant to margin formation models, plate kinematics, and frontier resource exploration. Although we disagree with Tikku and Direen's interpretations, the associated controversies reflect an absence of agreed-upon geophysical criteria for distinguishing stretched continental from oceanic crust, and a lack of samples from nonvolcanic margins.

Comment on "Major Australian-Antarctic plate reorganization at Hawaiian-Emperor bend time"

Whittaker et al. (Reports, 5 October 2007, p. 83) presented reconstructions for Australia and Antarctica showing a change in relative plate motion approximately 53 million years ago, coincident with an inferred major global plate reorganization. This comment addresses problematic areas in their assumptions and the geological consequences of their reconstructions.