Lower mantle structure from paleogeographically constrained dynamic Earth models

Seismic tomography reveals two large, low-shear velocity provinces (LLSVPs) beneath Africa and the Pacific Ocean. These structures may have existed for several 100 Myr and are likely compositionally distinct based on observed seismic characteristics interpreted in light of geodynamic models and mineral physics constraints. We investigate the dynamics of the LLSVPs through the use of evolutionary models of thermochemical convection from 250 Ma to present day. We use a spherical convection model in which the anomalous structures have a high bulk modulus, consistent with seismic interpretation. A new progressive assimilation method incorporates constraints from paleogeography using a refined plate history model (with 1 Myr time spacing) to guide the thermal structure of the lithosphere and steer the thermal evolution of slabs in the uppermost mantle. The thermochemical structures deform and migrate along the core-mantle boundary (CMB) through coupling to plate motions and in response to slab stresses. The models produce a ridge-like anomaly beneath Africa and a rounded pile beneath the Pacific Ocean, which at present day agrees with tomography, waveform modeling, and other geodynamic studies. Plumes emanate from the margins of the domes and ridges of thickened boundary layer between the domes. Dense and viscous slabs can undermine the stability of high bulk modulus structures at the CMB. High bulk modulus structures are not necessarily required to satisfy dynamic constraints on the LLSVPs

Insights on the kinematics of the India-Eurasia collision from global geodynamic models

The Eocene India-Eurasia collision is a first order tectonic event whose nature and chronology remains controversial. We test two end-member collision scenarios using coupled global plate motion-subduction models. The first, conventional model, invokes a continental collision soon after ∼60 Ma between a maximum extent Greater India and an Andean-style Eurasian margin. The alternative scenario involves a collision between a minimum extent Greater India and a NeoTethyan back-arc at ∼60 Ma that is subsequently subducted along southern Lhasa at an Andean-style margin, culminating with continent-continent contact at ∼40 Ma. Our numerical models suggest the conventional scenario does not adequately reproduce mantle structure related to Tethyan convergence. The alternative scenario better reproduces the discrete slab volumes and their lateral and vertical distribution in the mantle, and is also supported by the distribution of ophiolites indicative of Tethyan intraoceanic subduction, magmatic gaps along southern Lhasa and a two-stage slowdown of India. Our models show a strong component of southward mantle return flow for the Tethyan region, suggesting that the common assumption of near-vertical slab sinking is an oversimplification with significant consequences for interpretations of seismic tomography in the context of subduction reference frames

A Global Plate Model Including Lithospheric Deformation Along Major Rifts and Orogens Since the Triassic

Global deep‐time plate motion models have traditionally followed a classical rigid plate approach, even though plate deformation is known to be significant. Here we present a global Mesozoic–Cenozoic deforming plate motion model that captures the progressive extension of all continental margins since the initiation of rifting within Pangea at ~240 Ma. The model also includes major failed continental rifts and compressional deformation along collision zones. The outlines and timing of regional deformation episodes are reconstructed from a wealth of published regional tectonic models and associated geological and geophysical data. We reconstruct absolute plate motions in a mantle reference frame with a joint global inversion using hot spot tracks for the last 80 million years and minimizing global trench migration velocities and net lithospheric rotation. In our optimized model, net rotation is consistently below 0.2°/Myr, and trench migration scatter is substantially reduced. Distributed plate deformation reaches a Mesozoic peak of 30 × 106 km2 in the Late Jurassic (~160–155 Ma), driven by a vast network of rift systems. After a mid‐Cretaceous drop in deformation, it reaches a high of 48 x 106 km2 in the Late Eocene (~35 Ma), driven by the progressive growth of plate collisions and the formation of new rift systems. About a third of the continental crustal area has been deformed since 240 Ma, partitioned roughly into 65% extension and 35% compression. This community plate model provides a framework for building detailed regional deforming plate networks and form a constraint for models of basin evolution and the plate‐mantle system

Oceanic plateau subduction beneath North America and its geological and geophysical implications

We use two independent approaches, inverse models of mantle convection and plate reconstructions, to predict the temporal and spatial association of the Laramide events to subduction of oceanic plateaus. Inverse convection models, consistent with vertical motions in western US, recover two prominent anomalies on the Farallon plate during the Late Cretaceous that coincide with paleogeographically restored Shatsky and Hess conjugate plateaus when they collided with North America. The distributed deformation of the Laramide orogeny closely tracked the passage of the Shatsky conjugate massif, suggesting that subduction of this plateau dominated the distinctive geology of the western United States. Subduction of the Hess conjugate corresponds to termination of a Latest Cretaceous arc magmatism and intense crustal shortening in Early Paleogene in northwest Mexico. At present, conjugates of the Shatsky and Hess plateaus are located beneath the east coast of North America, and we predict that +4% seismic anomalies in P and S velocities are associated with the remnant plateaus with sharp lateral boundaries detectable by the USArray seismic experiment. Flat subduction of the Shatsky conjugate caused drastic subsidence/uplift and tilt of the Colorado Plateau (CP). From the inverse convection calculations, we find that with the arrival of the flat slab, dynamic subsidence starts at the southwestern CP and reaches a maximum at ~86 Ma. Two stages of uplift follow the removal of the Farallon slab: one in Latest Cretaceous and the other in Eocene with a cumulative uplift of ~1.2 km. The southwestern plateau reaches a high dynamic topography in the Eocene which is sustained to the present. Both the descent of the slab and buoyant upwelling may have contributed to late Cenozoic plateau uplift. The CP tilts downward to the NE before the Oligocene, caused by NE trending subduction of the Farallon slab. The NE tilt diminishes and switches to a SW tilt during the Miocene when buoyant mantle upwellings occur

Ridge subduction sparked reorganization of the Pacific plate-mantle system 60-50 million years ago

A reorganization centered on the Pacific plate occurred ~53–47 million years ago. A “top-down” plate tectonic mechanism, complete subduction of the Izanagi plate, as opposed to a “bottom-up” mantle flow mechanism, has been proposed as the main driver. Verification based on marine geophysical observations is impossible as most ocean crust recording this event has been subducted. Using a forward modeling approach, which assimilates surface plate velocities and shallow thermal structure of slabs into mantle flow models, we show that complete Izanagi plate subduction and margin-wide slab detachment induced a major change in sub-Pacific mantle flow, from dominantly southward before 60 Ma to north-northeastward after 50 Ma. Our results agree with onshore geology, mantle tomography, and the inferred motion of the Hawaiian hot spot and are consistent with a plate tectonic process driving the rapid plate-mantle reorganization in the Pacific hemisphere between 60 and 50 Ma. This reorganization is reflected in tectonic changes in the Pacific and surrounding ocean basins

GPlates – Building a Virtual Earth Through Deep Time

GPlates is an open‐source, cross‐platform plate tectonic geographic information system, enabling the interactive manipulation of plate‐tectonic reconstructions and the visualization of geodata through geological time. GPlates allows the building of topological plate models representing the mosaic of evolving plate boundary networks through time, useful for computing plate velocity fields as surface boundary conditions for mantle convection models and for investigating physical and chemical exchanges of material between the surface and the deep Earth along tectonic plate boundaries. The ability of GPlates to visualize subsurface 3‐D scalar fields together with traditional geological surface data enables researchers to analyze their relationships through geological time in a common plate tectonic reference frame. To achieve this, a hierarchical cube map framework is used for rendering reconstructed surface raster data to support the rendering of subsurface 3‐D scalar fields using graphics‐hardware‐accelerated ray‐tracing techniques. GPlates enables the construction of plate deformation zones—regions combining extension, compression, and shearing that accommodate the relative motion between rigid blocks. Users can explore how strain rates, stretching/shortening factors, and crustal thickness evolve through space and time and interactively update the kinematics associated with deformation. Where data sets described by geometries (points, lines, or polygons) fall within deformation regions, the deformation can be applied to these geometries. Together, these tools allow users to build virtual Earth models that quantitatively describe continental assembly, fragmentation and dispersal and are interoperable with many other mapping and modeling tools, enabling applications in tectonics, geodynamics, basin evolution, orogenesis, deep Earth resource exploration, paleobiology, paleoceanography, and paleoclimate

Community infrastructure and repository for marine magnetic identifications

Magnetic anomaly identifications underpin plate tectonic reconstructions and form the primary data set from which the age of the oceanic lithosphere and seafloor spreading regimes in the ocean basins can be determined. Although these identifications are an invaluable resource, their usefulness to the wider scientific community has been limited due to the lack of a central community infrastructure to organize, host, and update these interpretations. We have developed an open-source, community-driven online infrastructure as a repository for quality-checked magnetic anomaly identifications from all ocean basins. We provide a global sample data set that comprises 96,733 individually picked magnetic anomaly identifications organized by ocean basin and publication reference, and provide accompanying Hellingerformat files, where available. Our infrastructure is designed to facilitate research in plate tectonic reconstructions or research that relies on an assessment of plate reconstructions, for both experts and nonexperts alike. To further enhance the existing repository and strengthen its value, we encourage others in the community to contribute to this effort

A global plate model including lithospheric deformation along major rifts and orogens since the Triassic

Global deep‐time plate motion models have traditionally followed a classical rigid plate approach, even though plate deformation is known to be significant. Here we present a global Mesozoic–Cenozoic deforming plate motion model that captures the progressive extension of all continental margins since the initiation of rifting within Pangea at ~240 Ma. The model also includes major failed continental rifts and compressional deformation along collision zones. The outlines and timing of regional deformation episodes are reconstructed from a wealth of published regional tectonic models and associated geological and geophysical data. We reconstruct absolute plate motions in a mantle reference frame with a joint global inversion using hot spot tracks for the last 80 million years and minimizing global trench migration velocities and net lithospheric rotation. In our optimized model, net rotation is consistently below 0.2°/Myr, and trench migration scatter is substantially reduced. Distributed plate deformation reaches a Mesozoic peak of 30 × 10^6 km^2 in the Late Jurassic (~160–155 Ma), driven by a vast network of rift systems. After a mid‐Cretaceous drop in deformation, it reaches a high of 48 x 10^6 km^2 in the Late Eocene (~35 Ma), driven by the progressive growth of plate collisions and the formation of new rift systems. About a third of the continental crustal area has been deformed since 240 Ma, partitioned roughly into 65% extension and 35% compression. This community plate model provides a framework for building detailed regional deforming plate networks and form a constraint for models of basin evolution and the plate‐mantle system

Subduction controls the distribution and fragmentation of Earth’s tectonic plates

International audienceThe theory of plate tectonics describes how the surface of the Earth is split into an organized jigsaw of seven large plates 1 of similar sizes and a population of smaller plates, whose areas follow a fractal distribution 2,3. The reconstruction of global tectonics during the past 200 My 4 suggests that this layout is probably a long-term feature of our planet, but the forces governing it are unknown. Previous studies 3,5,6 , primarily based on statistical properties of plate distributions, were unable to resolve how the size of plates is determined by lithosphere properties and/or underlying mantle convection. Here, we demonstrate that the plate layout of the Earth is produced by a dynamic feedback between mantle convection and the strength of the lithosphere. Using 3D spherical models of mantle convection with plate-like behaviour that match the plate size-frequency distribution observed for Earth, we show that subduction geometry drives the tectonic fragmentation that generates plates. The spacing between slabs controls the layout of large plates, and the stresses caused by the bending of trenches, break plates into smaller fragments. Our results explain why the fast evolution in small back-arc plates 7,8 reflects the dramatic changes in plate motions during times of major reorganizations. Our study opens the way to use convection simulations with plate-like behaviour to unravel how global tectonics and mantle convection are dynamically connected

Kinematics and extent of the Piemont–Liguria Basin – implications for subduction processes in the Alps

Assessing the size of a former ocean of which only remnants are found in mountain belts is challenging but crucial to understanding subduction and exhumation processes. Here we present new constraints on the opening and width of the Piemont–Liguria (PL) Ocean, known as the Alpine Tethys together with the Valais Basin. We use a regional tectonic reconstruction of the Western Mediterranean–Alpine area, implemented into a global plate motion model with lithospheric deformation, and 2D thermo-mechanical modeling of the rifting phase to test our kinematic reconstructions for geodynamic consistency. Our model fits well with independent datasets (i.e., ages of syn-rift sediments, rift-related fault activity, and mafic rocks) and shows that, between Europe and northern Adria, the PL Basin opened in four stages: (1) rifting of the proximal continental margin in the Early Jurassic (200–180 Ma), (2) hyper-extension of the distal margin in the Early to Middle Jurassic (180–165 Ma), (3) ocean–continent transition (OCT) formation with mantle exhumation and MORB-type magmatism in the Middle–Late Jurassic (165–154 Ma), and (4) breakup and mature oceanic spreading mostly in the Late Jurassic (154–145 Ma). Spreading was slow to ultra-slow (max. 22 mm yr−1, full rate) and decreased to ∼5 mm yr−1 after 145 Ma while completely ceasing at about 130 Ma due to the motion of Iberia relative to Europe during the opening of the North Atlantic. The final width of the PL mature (“true”) oceanic crust reached a maximum of 250 km along a NW–SE transect between Europe and northwestern Adria. Plate convergence along that same transect has reached 680 km since 84 Ma (420 km between 84–35 Ma, 260 km between 35–0 Ma), which greatly exceeds the width of the ocean. We suggest that at least 63 % of the subducted and accreted material was highly thinned continental lithosphere and most of the Alpine Tethys units exhumed today derived from OCT zones. Our work highlights the significant proportion of distal rifted continental margins involved in subduction and exhumation processes and provides quantitative estimates for future geodynamic modeling and a better understanding of the Alpine Orogeny