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

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

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

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

Topographic asymmetry of the South Atlantic from global models of mantle flow and lithospheric stretching

The relief of the South Atlantic is characterized by elevated passive continental margins along southern Africa and eastern Brazil, and by the bathymetric asymmetry of the southern oceanic basin where the western flank is much deeper than the eastern flank. We investigate the origin of these topographic features in the present and over time since the Jurassic with a model of global mantle flow and lithospheric deformation. The model progressively assimilates plate kinematics, plate boundaries and lithospheric age derived from global tectonic reconstructions with deforming plates, and predicts the evolution of mantle temperature, continental crustal thickness, long-wavelength dynamic topography, and isostatic topography. Mantle viscosity and the kinematics of the opening of the South Atlantic are adjustable parameters in thirteen model cases. Model predictions are compared to observables both for the present-day and in the past. Present-day predictions are compared to topography, mantle tomography, and an estimate of residual topography. Predictions for the past are compared to tectonic subsidence from backstripped borehole data along the South American passive margin, and to dynamic uplift as constrained by thermochronology in southern Africa. Comparison between model predictions and observations suggests that the first-order features of the topography of the South Atlantic are due to long-wavelength dynamic topography, rather than to asthenospheric processes. The uplift of southern Africa is best reproduced with a lower mantle that is at least 40 times more viscous than the upper mantle

Geodynamic reconstruction of an accreted Cretaceous back-arc basin in the Northern Andes

A complex history of subduction, back-arc basin formation, terrane accretion and transpressional shearing characterizes the evolution of the Caribbean and northern South American margin since Jurassic times. Quantitative plate tectonic reconstructions of the area do not include Jurassic-Cretaceous back-arc terranes of which there are both geological and geophysical observations. We developed a revised plate tectonic reconstruction based on geological observations and seismic tomography models to constrain the Jurassic-Cretaceous subduction history of eastern Panthalassa, along the western margin of the Caribbean region. This reconstruction considers the opening of a Northern Andean back-arc basin at 145 Ma, the Quebradagrande back-arc, closing at 120 Ma and followed by terrane accretion and northward translation along the South American margin starting at 100 Ma. This kinematic reconstruction is tested against two previously published tectonic reconstructions via coupling with global numerical mantle convection models using CitcomS. A comparison of modelled versus tomographically imaged mantle structure reveals that subduction outboard of the South American margin, lacking in previous tectonic models, is required to reproduce mid-mantle positive seismic anomalies imaged in P- and S-wave seismic tomography beneath South America, 500–2000 km in depth. Furthermore, we show that this subduction zone is likely produced by a back-arc basin that developed along the northern Andes during the Cretaceous via trench roll-back from 145 Ma and was closed at 100 Ma. The contemporaneous opening of the Quebradagrande back-arc basin with the Rocas Verdes back-arc basin in the southern Andes is consistent with a model that invokes return flow of mantle material behind a retreating slab and may explain why extension along the Peruvian and Chilean sections of the Andean margin did not experience full crustal break-up and back-arc opening during the late Jurassic-early Cretaceous Period.This research was undertaken with the assistance of the Sydney Informatics Hub in accessing resources from the National Computational Infrastructure (NCI), which is supported by the Australian Government. The authors were supported by Australian Research Council grants IH130200012 (RDM, NF), FT130101564 (MS) and DE160101020 (NF). We would like to thank Wouter Schellart and an anonymous reviewer for their helpful comments, which improved the manuscript

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

Tectonic evolution of Western Tethys from Jurassic to present day: coupling geological and geophysical data with seismic tomography models

The geodynamic evolution of the Western Tethys is characterized by multiple phases of rifting, seafloor spreading, subduction, and collisional events. Regional reconstructions are highly dependent on the kinematic history of the major plates bounding the Atlantic and Tethyan tectonic domains, as well as small micro-plates resulted from the fragmentation of northern Gondwanaland. The complexity of tectonic events in this area leads to major discrepancies between competing models about the timing, location, and polarity of subduction zones, for both the Cenozoic evolution and earlier phases. We focus on unravelling the Mesozoic evolution of the Western Tethys. We first reassessed kinematic models for the Early Jurassic–Late Cretaceous opening of the central, north central, and north Atlantic and used these as boundary conditions on the kinematic reconstructions of the Tethyan realm. We combined reconstructions of rifting and early seafloor spreading in northern Pangea that incorporate quantitative estimates of continental extension, and suggest a transtensional motion of Iberia relative to Europe in Early Cretaceous time to fit within the refined plate configuration of Central North Atlantic. We combined this regional framework with a recently published model for the motion of smaller blocks within the Western Tethys; from this model, we created synthetic isochrons for extinct oceanic basins and built evolving topological plate boundaries based on the new rigid plate model to derive a self-consistent and time-dependant model for the last 200 million years. We then examined the consistency of subduction history implied by the kinematic reconstructions, by comparing reconstructed plate boundary configurations to mantle velocity structure imaged by a range of seismic tomography models. Our results show that a satisfactory match can be made between Cenozoic subduction events in the Western Tethys region and observed shallow tomographic high-velocity material. However, the match is less clear for older subducted material. Correlations between surface reconstructions and deep Earth structure suggest that mid-deep mantle seismic features under present day Northeast-Central and Northwest Africa-Arabia may correspond with the Mesozoic subduction systems in the Vardar Ocean, Alpine Tethys, and Western Neotethys, respectively. These correlations support a model with intra-oceanic subduction of the Vardar Ocean from Middle Jurassic to Early Cretaceous, and mid-Early Cretaceous initiation of oceanic subduction in the Ligurian-Piemont Ocean. The results from the plate-tomography comparison suggest the existence of oceanic subduction in the Alpine Tethys Oceans in late early Cretaceous time. We investigated the uncertainties in the tectonic model in terms of absolute and relative plate motion and surface velocities and showed that choice of absolute reference frame can partially account for the lateral offset between the Vardar subduction zone and associated slab material in deep mantle; additional mismatches may be attributable to the limitations of our methodology, such as the assumption that slabs sink vertically. © 2016 Informa UK Limited, trading as Taylor & Francis Group.Australian Research Council via [Grant number FL0992245]; (MH, SW, RDM) and [grant number FT130101564] (MS)

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