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.

Massive lithospheric delamination in southeastern Tibet facilitating continental extrusion

Significant left-lateral movement along the Ailao Shan-Red River fault accommodated a substantial amount of the late Eocene to early Miocene India-Asia convergence. However, the activation of this critical strike-slip fault remains poorly understood. Here, we show key seismic evidence for the occurrence of massive lithospheric delamination in southeastern Tibet. Our novel observation of reflected body waves (e.g. P₄₁₀P and P₆₆₀P) retrieved from ambient noise interferometry sheds new light on the massive foundered lithosphere currently near the bottom of the mantle transition zone beneath southeastern Tibet. By integrating the novel seismic and pre-existing geochemical observations, we highlight a linkage between massive lithospheric delamination shortly after the onset of hard collision and activation of continental extrusion along the Ailao Shan-Red River fault. This information provides critical insight into the early-stage evolution of the India-Asia collision in southeastern Tibet, which has significant implications for continental collision and its intracontinental response.

Mantle-flow diversion beneath the Iranian plateau induced by Zagros' lithospheric keel

Previous investigation of seismic anisotropy indicates the presence of a simple mantle flow regime beneath the Turkish-Anatolian Plateau and Arabian Plate. Numerical modeling suggests that this simple flow is a component of a large-scale global mantle flow associated with the African superplume, which plays a key role in the geodynamic framework of the Arabia-Eurasia continental collision zone. However, the extent and impact of the flow pattern farther east beneath the Iranian Plateau and Zagros remains unclear. While the relatively smoothly varying lithospheric thickness beneath the Anatolian Plateau and Arabian Plate allows progress of the simple mantle flow, the variable lithospheric thickness across the Iranian Plateau is expected to impose additional boundary conditions on the mantle flow field. In this study, for the first time, we use an unprecedented data set of seismic waveforms from a network of 245 seismic stations to examine the mantle flow pattern and lithospheric deformation over the entire region of the Iranian Plateau and Zagros by investigation of seismic anisotropy. We also examine the correlation between the pattern of seismic anisotropy, plate motion using GPS velocities and surface strain fields. Our study reveals a complex pattern of seismic anisotropy that implies a similarly complex mantle flow field. The pattern of seismic anisotropy suggests that the regional simple mantle flow beneath the Arabian Platform and eastern Turkey deflects as a circular flow around the thick Zagros lithosphere. This circular flow merges into a toroidal component beneath the NW Zagros that is likely an indicator of a lateral discontinuity in the lithosphere. Our examination also suggests that the main lithospheric deformation in the Zagros occurs as an axial shortening across the belt, whereas in the eastern Alborz and Kopeh-Dagh a belt-parallel horizontal lithospheric deformation plays a major role.

Ongoing formation of felsic lower crustal channel by relamination in Zagros collision zone revealed from regional tomography

Complex interaction of rheologically contrasting layers within the lithosphere during the collision of continental plates leads to active faulting, which represents a serious hazard to the population and infrastructure. One of the collision scenarios presumes the existence of a middle-lower crustal channel composed of subducted silicic upper crustal rocks, which is thought to exist in the Tibetan-Himalayan system. Based on the results of seismic tomography, we argue that a similar mechanism of crustal channeling takes place beneath the Zagros mountain system in southwestern Iran. The 3D seismic velocity model reveals an inverted crustal architecture of the collision zone, in which the low-velocity felsic (granitic and sedimentary) upper crustal rocks of the Arabian plate form a seismically inactive lower crustal channel below the higher-velocity mafic (basaltic) middle-upper crustal layer of the Iranian crust. Based on existing numerical models, we suggest that the formation of the felsic channel is likely governed by separation (delamination) of the weak felsic upper crust of the subducting Arabian lithosphere and its ductile underplating under rheologically stronger upper-middle crust of the Iranian plate.

The convergence history of India-Eurasia records multiple subduction dynamics processes

During the Cretaceous, the Indian plate moved towards Eurasia at the fastest rates ever recorded. The details of this journey are preserved in the Indian Ocean seafloor, which document two distinct pulses of fast motion, separated by a noticeable slowdown. The nature of this rapid acceleration, followed by a rapid slowdown and then succeeded by a second speedup, is puzzling to explain. Using an extensive observation dataset and numerical models of subduction, we show that the arrival of the Reunion mantle plume started a sequence of events that can explain this history of plate motion. The forces applied by the plume initiate an intra-oceanic subduction zone, which eventually adds enough additional force to drive the plates at the anomalously fast speeds. The two-stage closure of a double subduction system, including accretion of an island arc at 50 million years ago, may help reconcile geological evidence for a protracted India-Eurasia collision.

Southward propagation of Nazca subduction along the Andes

The Andean margin is the plate-tectonic paradigm for long-lived, continuous subduction, yet its geology since the late Mesozoic era (the past 100 million years or so) has been far from steady state. The episodic deformation and magmatism have been attributed to cyclic changes in the dip angle of the subducting slab, slab break-off and the penetration of the slab into the lower mantle; the role of plate tectonics remains unclear, owing to the extensive subduction of the Nazca-Farallon plate (which has resulted in more than 5,500 kilometres of lithosphere being lost to the mantle). Here, using tomographic data, we recreate the plate-tectonic geometry of the subducted Nazca slab, which enables us to reconstruct Andean plate tectonics since the late Mesozoic. Our model suggests that the current phase of Nazca subduction began at the northern Andes (5° S) during the late Cretaceous period (around 80 million years ago) and propagated southwards, reaching the southern Andes (40° S) by the early Cenozoic era (around 55 million year ago). Thus, contrary to the current paradigm, Nazca subduction has not been fully continuous since the Mesozoic but instead included episodic divergent phases. In addition, we find that foredeep sedimentation and the initiation of Andean compression are both linked to interactions between the Nazca slab and the lower mantle, consistent with previous modelling.

Stratigraphy, Tectonics and Hydrocarbon Habitat of the Abadan Plain Basin: A Geological Review of a Prolific Middle Eastern Hydrocarbon Province

The Abadan Plain Basin is located in the Middle East region which is host to some of the world’s largest oil and gas fields around the Persian Gulf. This basin is a foredeep basin to the southwest of the Zagros Fold-Thrust-Belt, bounded along its northern and eastern margins by the Dezful Embayment. Most of the rocks in this basin have been deposited in a carbonate environment, and existing fractures have made the formations a favourable place for hydrocarbon accumulations. The basin is enriched by oil and, therefore, gas reservoirs are few, and some of the explored reservoirs exhibit significant degrees of overpressure. This paper has compiled several aspects of the Abadan Plain Basin tectonics, structural geology and petroleum systems to provide a better understanding of the opportunities and risks of development activities in this region. In addition to the existing knowledge, this paper provides a basin-wide examination of pore pressure, vertical stress, temperature gradient, and wellbore stability issues.

Mantle Sources of Recent Anatolian Intraplate Magmatism: A Regional Plume or Local Tectonic Origin?

We present an extensive study of rehomogenized olivine-hosted melt inclusions, olivine phenocrysts, and chromian spinel inclusions to explore the link between geodynamic conditions and the origin and composition of Pliocene-Quaternary intraplate magmatism in Anatolia at Kula, Ceyhan-Osmaniye, and Karacadağ. Exceptional compositional variability of these products reveals early and incomplete mixing of distinct parental melts in each volcanic center, reflecting asthenospheric and lithospheric mantle sources. The studied primitive magmas consist of (1) two variably enriched ocean island basalt (OIB)-type melts in Kula; (2) both OIB-type and plume mid-ocean ridge basalt (P-MORB)-like melts beneath Toprakkale and Üçtepeler (Ceyhan-Osmaniye); and (3) two variably enriched OIB-type melts beneath Karacadağ. Estimated conditions of primary melt generation are 23-9 kbar, 75-30 km, and 1415-1215 °C for Kula; 28-19 kbar, 90-65 km, and 1430-1350 °C for Toprakkale; 23-18 kbar, 75-60 km, and 1400-1355 °C for Üçtepeler; and 35-27 kbar, 115-90 km, and 1530-1455 °C for Karacadağ, the deepest levels of which correspond to the depth of the lithosphere-asthenosphere boundary in all regions. Although magma ascent was likely facilitated by local deformation structures, recent Anatolian intraplate magmatism seems to be triggered by large-scale mantle flow that also affects the wider Arabian and North African regions. We infer that these volcanics form part of a much wider Arabian-North African intraplate volcanic province, which was able to invade the Anatolian upper plate through slab gaps.

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.

On the consistency of seismically imaged lower mantle slabs

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000-1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.

Pangea breakup and northward drift of the Indian subcontinent reproduced by a numerical model of mantle convection

Since around 200 Ma, the most notable event in the process of the breakup of Pangea has been the high speed (up to 20 cm yr(-1)) of the northward drift of the Indian subcontinent. Our numerical simulations of 3-D spherical mantle convection approximately reproduced the process of continental drift from the breakup of Pangea at 200 Ma to the present-day continental distribution. These simulations revealed that a major factor in the northward drift of the Indian subcontinent was the large-scale cold mantle downwelling that developed spontaneously in the North Tethys Ocean, attributed to the overall shape of Pangea. The strong lateral mantle flow caused by the high-temperature anomaly beneath Pangea, due to the thermal insulation effect, enhanced the acceleration of the Indian subcontinent during the early stage of the Pangea breakup. The large-scale hot upwelling plumes from the lower mantle, initially located under Africa, might have contributed to the formation of the large-scale cold mantle downwelling in the North Tethys Ocean.

Plate tectonic controls on atmospheric CO2 levels since the Triassic

Climate trends on timescales of 10s to 100s of millions of years are controlled by changes in solar luminosity, continent distribution, and atmosphere composition. Plate tectonics affect geography, but also atmosphere composition through volcanic degassing of CO2 at subduction zones and midocean ridges. So far, such degassing estimates were based on reconstructions of ocean floor production for the last 150 My and indirectly, through sea level inversion before 150 My. Here we quantitatively estimate CO2 degassing by reconstructing lithosphere subduction evolution, using recent advances in combining global plate reconstructions and present-day structure of the mantle. First, we estimate that since the Triassic (250-200 My) until the present, the total paleosubduction-zone length reached up to ∼200% of the present-day value. Comparing our subduction-zone lengths with previously reconstructed ocean-crust production rates over the past 140 My suggests average global subduction rates have been constant, ∼6 cm/y: Higher ocean-crust production is associated with longer total subduction length. We compute a strontium isotope record based on subduction-zone length, which agrees well with geological records supporting the validity of our approach: The total subduction-zone length is proportional to the summed arc and ridge volcanic CO2 production and thereby to global volcanic degassing at plate boundaries. We therefore use our degassing curve as input for the GEOCARBSULF model to estimate atmospheric CO2 levels since the Triassic. Our calculated CO2 levels for the mid Mesozoic differ from previous modeling results and are more consistent with available proxy data.

Greater India Basin hypothesis and a two-stage Cenozoic collision between India and Asia

Cenozoic convergence between the Indian and Asian plates produced the archetypical continental collision zone comprising the Himalaya mountain belt and the Tibetan Plateau. How and where India-Asia convergence was accommodated after collision at or before 52 Ma remains a long-standing controversy. Since 52 Ma, the two plates have converged up to 3,600 ± 35 km, yet the upper crustal shortening documented from the geological record of Asia and the Himalaya is up to approximately 2,350-km less. Here we show that the discrepancy between the convergence and the shortening can be explained by subduction of highly extended continental and oceanic Indian lithosphere within the Himalaya between approximately 50 and 25 Ma. Paleomagnetic data show that this extended continental and oceanic "Greater India" promontory resulted from 2,675 ± 700 km of North-South extension between 120 and 70 Ma, accommodated between the Tibetan Himalaya and cratonic India. We suggest that the approximately 50 Ma "India"-Asia collision was a collision of a Tibetan-Himalayan microcontinent with Asia, followed by subduction of the largely oceanic Greater India Basin along a subduction zone at the location of the Greater Himalaya. The "hard" India-Asia collision with thicker and contiguous Indian continental lithosphere occurred around 25-20 Ma. This hard collision is coincident with far-field deformation in central Asia and rapid exhumation of Greater Himalaya crystalline rocks, and may be linked to intensification of the Asian monsoon system. This two-stage collision between India and Asia is also reflected in the deep mantle remnants of subduction imaged with seismic tomography.

Oligocene–Miocene basin evolution in SE Anatolia, Turkey: constraints on the closure of the eastern Tethys gateway

The Oligocene–Miocene was a time characterized by major climate changes as well as changing plate configurations. The Middle Miocene Climate Transition (17 to 11 Ma) may even have been triggered by a plate tectonic event: the closure of the eastern Tethys gateway, the marine connection between the Mediterranean and Indian Ocean. To address this idea, we focus on the evolution of Oligocene and Miocene foreland basins in the southernmost part of Turkey, the most likely candidates to have formed this gateway. In addition, we take the geodynamic evolution of the Arabian–Eurasian collision into account.

The Muş and Elazığ basins, located to the north of the Bitlis–Zagros suture zone, were most likely connected during the Oligocene. The deepening of both basins is biostratigraphically dated by us to occur during the Rupelian (Early Oligocene). Deep marine conditions (between 350 and 750 m) prevailed until the Chattian (Late Oligocene), when the basins shoaled rapidly to subtidal/intertidal environment in tropical to subtropical conditions, as indicated by the macrofossil assemblages. We conclude that the emergence of this basin during the Chattian severely restricted the marine connection between an eastern (Indian Ocean) and western (Mediterranean) marine domain. If a connection persisted it was likely located south of the Bitlis–Zagros suture zone. The Kahramanmaraş basin, located on the northern Arabian promontory south of the Bitlis–Zagros suture zone, was a foreland basin during the Middle and Late Miocene, possibly linked to the Hatay basin to the west and the Lice basin to the east. Our data indicates that this foreland basin experienced shallow marine conditions during the Langhian, followed by a rapid deepening during Langhian/Serravallian and prevailing deep marine conditions (between 350 and 750 m) until the early Tortonian. We have dated the youngest sediments underneath a subduction-related thrust at c. 11 Ma and suggest that this corresponds to the end of underthrusting in the Kahramanmaraş region, i.e. the end of subduction of Arabia. This age coincides in time with the onset of eastern Anatolian volcanism, uplift of the East Anatolian Accretionary Complex, and the onset of the North and East Anatolian Fault Zones accommodating westward escape tectonics of Anatolia. After c. 11 Ma, the foreland basin south of the Bitlis formed not (or no longer) a deep marine connection along the northern margin of Arabia between the Mediterranean Sea and the Indian Ocean. We finally conclude that a causal link between gateway closure and global climate change to a cooler mode, recorded in the Mi3b event (δ18O increase) dated at 13.82 Ma, cannot be supported.