Search

The India–Asia collision zone is the archetype to calibrate geological responses to continent–continent collision, but hosts a paradox: there is no orogen-wide geological record of oceanic subduction after initial collision around 60–55 Ma, yet thousands of kilometers of post-collisional subduction occurred before the arrival of unsubductable continental lithosphere that currently horizontally underlies Tibet. Kinematically restoring incipient horizontal underthrusting accurately predicts geologically estimated diachronous slab break-off, unlocking the Miocene of Himalaya–Tibet as a natural laboratory for unsubductable lithosphere convergence. Additionally, three endmember paleogeographic scenarios exist with different predictions for the nature of post-collisional subducted lithosphere but each is defended and challenged based on similar data types. This paper attempts at breaking through this impasse by identifying how the three paleogeographic scenarios each challenge paradigms in geodynamics, orogenesis, magmatism or paleogeographic reconstruction and identify opportunities for methodological advances in paleomagnetism, sediment provenance analysis, and seismology to conclusively constrain Greater Indian paleogeography.

The Panthalassa Ocean, which surrounded the late Paleozoic-early Mesozoic Pangea supercontinent, was underlain by multiple tectonic plates that have since been lost to subduction. In this study, we develop an approach to reconstruct plate motions of this subducted lithosphere utilizing paleomagnetic data from accreted Ocean Plate Stratigraphy (OPS). We first establish the boundaries of the Panthalassa domain by using available Indo-Atlantic plate reconstructions and restorations of complex plate boundary deformation at circum-Panthalassa trenches. We reconstruct the Pacific Plate and its conjugates, the Farallon, Phoenix, and Izanagi plates, back to 190 Ma using marine magnetic anomaly records of the modern Pacific. Then, we present new and review published paleomagnetic data from OPS exposed in the accretionary complexes of Cedros Island (Mexico), the Santa Elena Peninsula (Costa Rica), the North Island of New Zealand, and Japan. These data provide paleolatitudinal plate motion components of the Farallon, Phoenix and Izanagi plates, and constrain the trajectories of these plates from their spreading ridges towards the trenches in which they subducted. For 83 to 150 Ma, we use two independent mantle frames to connect the Panthalassa plate system to the Indo-Atlantic plate system and test the feasibility of this approach with the paleomagnetic data. For times prior to 150 Ma, and as far back as Permian time, we reconstruct relative and absolute Panthalassa plate motions such that divergence is maintained between the Izanagi, Farallon and Phoenix plates, convergence is maintained with Pangean continental margins in Japan, Mexico and New Zealand, and paleomagnetic constraints are met. The reconstruction approach developed here enables data-based reconstruction of oceanic plates and plate boundaries in the absence of marine magnetic anomaly data or mantle reference frames, using Ocean Plate Stratigraphy, paleo-magnetism, and constraints on the nature of circum-oceanic plate boundaries. Such an approach is a crucial next step towards quantitative reconstruction of the currently largely unknown tectonic evolution of the Earth's oceanic domains in deep geological time.

International audience; Subduction polarity reversal during arc-continent collision has been proposed as a key mechanism to initiate new subduction zones. Despite often interpreted, well-exposed geological record that document the reversal is sparse. The ophiolitic lithounits of the Andaman and Nicobar Islands have been proposed to have formed during the initiation of a new subduction zone following the collision of the Woyla Arc of Sumatra with Sundaland (Eurasia). We here present new field, petrological and geochronological data to evaluate the timing of the initiation of Andaman subduction. We targeted the previously inferred but unstudied metamorphic sole of the Andaman ophiolites that witnessed juvenile subduction. Thermodynamic modeling reveals that the exposed amphibolites of the sole formed at around 0.9 GPa and 675°C. We dated two samples of the metamorphic sole using the Ar/Ar method on amphibole, giving cooling ages of 106.4 ± 2.1 and 105.3 ± 1.6 Ma. This is similar to published ages from plagioclase xenocrysts in recent Barren Island volcanics and in zircons from a gabbro sample from the Andaman ophiolite, which we interpret as the age of the original ophiolite formation. The Ar/Ar ages are considerably older than arc magmatic gabbros and plagiogranites of the overlying ophiolite previously dated at 99-93 Ma and thought to reflect the ophiolite age but recently reinterpreted as a volcanic arc built on the ophiolite. Combined with the ages of Woyla-Sundaland collision, we argue that subduction polarity reversal occurred in a transient period of perhaps some 10 Myr, similar to recent settings.

The puzzling

The Tibetan Plateau accommodated major upper crustal shortening during Indian Plate oceanic and continental lithosphere subduction. Deciphering whether shortening was continuous or episodic, and how it correlates to major geodynamic changes is challenging. Here we apply anisotropy of magnetic susceptibility (AMS), a sensitive synsedimentary strain indicator, to a ~3 km thick magnetostratigraphically dated sedimentary section (69–41.5 Ma) in eastern Tibet. AMS shows “earliest deformation” fabrics from 69–52 Ma, followed by a sudden change to “pencil structure” fabrics with increasing anisotropy degree at ~52 Ma, dating a sudden increased synsedimentary shortening strain. This change coincides with enhanced sedimentation rates and synsedimentary vertical-axis rotations of the Gonjo Basin, suggesting a causal link to a marked India-Asia convergence rate deceleration. We show that AMS analysis provides a strong tool to distinguish between climatic and tectonic causes of sedimentological change and is an asset in identifying discrete tectonic pulses in intensely deformed terrane.

We present a map that correlates tectonic units between Alps and western Turkey accompanied by a text providing access to literature data, explaining the concepts used for defining the mapped tectonic units, and first-order paleogeographic inferences. Along-strike similarities and differences of the Alpine-Eastern Mediterranean orogenic system are discussed. The map allows (1) for superimposing additional information, such as e.g., post-tectonic sedimentary basins, manifestations of magmatic activity, onto a coherent tectonic framework and (2) for outlining the major features of the Alpine-Eastern Mediterranean orogen. Dinarides-Hellenides, Anatolides and Taurides are orogens of opposite subduction polarity and direction of major transport with respect to Alps and Carpathians, and polarity switches across the Mid-Hungarian fault zone. The Dinarides-Hellenides-Taurides (and Apennines) consist of nappes detached from the Greater Adriatic continental margin during Cretaceous and Cenozoic orogeny. Internal units form composite nappes that passively carry ophiolites obducted in the latest Jurassic–earliest Cretaceous or during the Late Cretaceous on top of the Greater Adriatic margin successions. The ophiolites on top of composite nappes do not represent oceanic sutures zones, but root in the suture zones of Neotethys that formed after obduction. Suturing between Greater Adria and the northern and eastern Neotethys margin occupied by the Tisza and Dacia mega-units and the Pontides occurred in the latest Cretaceous along the Sava-Izmir-Ankara-Erzincan suture zones. The Rhodopian orogen is interpreted as a deep-crustal nappe stack formed in tandem with the Carpatho-Balkanides fold-thrust belt, now exposed in a giant core complex exhumed in late Eocene to Miocene times from below the Carpatho-Balkan orogen and the Circum-Rhodope unit. Its tectonic position is similar to that of the Sakarya unit of the Pontides. We infer that the Rhodope nappe stack formed due to north-directed thrusting. Both Rhodopes and Pontides are suspected to preserve the westernmost relics of the suture zone of Paleotethys.

International audience; As subduction zones and their related processes are often studied in 2D, or cylindrical 3D sections, the dynamic effects of trench curvature and its evolution through time remain under-explored. Whereas temporal variations in trench trend may be estimated through restoring upper plate deformation, we investigate the forearc deformation history of the strongly curved northern Lesser Antilles trench, connecting the near-orthogonal Lesser Antilles subduction zone with the Motagua-Cayman transform plate boundary. Our new paleomagnetic dataset consists of 310 cores from Eo-Oligocene magmatic rocks and limestones from St. Barthélemy Island. The limestones yielded a post-folding magnetization containing a similar magnetic direction to those stored in magmatic rocks that intrude the folded carbonates, both indicating a post-Oligocene ~15°, and perhaps up to 25° counterclockwise rotation of the island. Our results highlight that the present-day trench curvature formed progressively during the Cenozoic, allowing us to discuss different tectonic scenarios explaining NE Caribbean plate deformation, and to identify key targets for future research on tectonic architecture and the potential present-day activity of intra-plate deformation that may pose seismic hazards.

Central Anatolia (Turkey) is a small and nascent example of a high orogenic plateau, providing a natural laboratory to study processes driving plateau rise. The 1-km-high plateau interior uplifted since c. 8–5 Ma, with a further phase of kilometre-scale uplift affecting the southern plateau margin since 0.45 Ma. Several causes of plateau rise have been proposed: peeling or dripping delamination of the lithospheric mantle; asthenospheric upwelling through slab gaps created by slab fragmentation or break-off, and; continental underthrusting and crustal shortening below the southern plateau margin. The Neogene history of the plateau has not been diagnostic of the causes of plateau rise. We thus evaluate proposed uplift causes in the context of the Anatolian orogenesis, which formed the plateau lithosphere during subduction since the Cretaceous. We combine this analysis with available constraints on uplift, and geophysical data that illuminate the modern mantle (and crustal) structure. Our analysis suggests that lithospheric dripping, which followed arc magmatism and shortening in the Kirsehir Block (eastern Central Anatolia), is the most likely cause of plateau interior uplift. Lithospheric dripping is, however, an unlikely sole driver of multi-phase uplift along the southern plateau margin. There, underthrusting of the African continental margin, recorded by c. 11–7 Ma thrusting on Cyprus, is a viable cause of uplift since 0.45 Ma, but cannot account for earlier uplift since c. 8–5 Ma. Instead, slab break-off below the southern plateau margin is likely in light of geophysical data. On the SW plateau margin, small-scale peeling delamination of the Central Taurides by the Antalya slab since early Miocene times accounts for >150 km slab retreat with no corresponding upper-plate deformation. A southwest-travelling wave of subsidence and uplift signalled this retreat and may have contributed to coeval oroclinal bending of the western Central Taurides and southeastward thrusting of the Lycian Nappes.

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., Key Points We present new paleomagnetic data from Eocene‐Miocene sedimentary rocks in Kutai Basin, Borneo, and compile all available paleomagnetic data from Mesozoic and Cenozoic rocks in Sundaland and BorneoWe compare reconstructed trenches with seismic tomography to evaluate the regional or local nature of paleomagnetic resultsWe provide a kinematic restoration of Borneo and Sundaland since the Eocene that is consistent with paleomagnetism and tomography

International audience; Abstract. The geotherm in subduction zones is thought to vary as a function of the subduction rate and the age of the subducting lithosphere. Along a single subduction zone the rate of subduction may strongly vary due to changes in the angle between the trench and the plate convergence vector, i.e., the subduction obliquity, due to trench curvature. We currently observe such curvature in, e.g., the Marianas, Chile and Aleutian trenches. Recently, strong along-strike variations in subduction obliquity were proposed to have caused a major temperature contrast between Cretaceous geological records of western and central Turkey. We test here whether first-order temperature variation in a subduction zone may be caused by variation in the trench geometry using simple thermo-kinematic finite-element 3-D numerical models. We prescribe the trench geometry by means of a simple mathematical function and compute the mantle flow in the mantle wedge by solving the equation of mass and momentum conservation. We then solve the energy conservation equation until steady state is reached. We analyze the results (i) in terms of mantle wedge flow with emphasis on the trench-parallel component and (ii) in terms of temperature along the plate interface by means of maps and the depth–temperature path at the interface. In our experiments, the effect of the trench curvature on the geotherm is substantial. A small obliquity yields a small but not negligible trench-parallel mantle flow, leading to differences of 30 °C along-strike of the model. Advected heat causes such temperature variations (linked to the magnitude of the trench-parallel component of velocity). With increasing obliquity, the trench-parallel component of the velocity consequently increases and the temperature variation reaches 200 °C along-strike. Finally, we discuss the implication of our simulations for the ubiquitous oblique systems that are observed on Earth and the limitations of our modeling approach. Lateral variations in plate sinking rate associated with curvature will further enhance this temperature contrast. We conclude that the synchronous metamorphic temperature contrast between central and western Turkey may well have resulted from reconstructed major variations in subduction obliquity.

The recent tectonic evolution of the western Mediterranean region is enigmatic. The causes for the closure of the Moroccan marine gateway prior to the Messinian salinity crisis, for the ongoing shortening of the Moroccan Rif and for the origin of the seismogenic Trans-Alboran shear zone and eastern Betics extension are unclear. These puzzling tectonic features cannot be fully explained by subduction of the east-dipping Gibraltar slab in the context of the regional relative plate motion frame. Here we use a combination of geological and geodetic data, as well as three-dimensional numerical modelling of subduction, to show that these unusual tectonic features could be the consequence of slab dragging—the north to north-eastward dragging of the Gibraltar slab by the absolute motion of the African Plate. Comparison of our model results to patterns of deformation in the western Mediterranean constrained by geological and geodetic data confirm that slab dragging provides a plausible mechanism for the observed deformation. Our results imply that the impact of absolute plate motion on subduction is identifiable from crustal observations. Identifying such signatures elsewhere may improve the mantle reference frame and provide insights on subduction evolution and associated crustal deformation.

Most ophiolites have geochemical signatures that indicate formation by suprasubduction seafloor spreading above newly initiated subduction zones, and hence they record fore-arc processes operating following subduction initiation. They are frequently underlain by a metamorphic sole formed at the top of the downgoing plate and accreted below the overlying suprasubduction zone lithosphere immediately following ophiolite formation. Paleomagnetic analyses of ophiolites can provide important insights into the enigmatic geodynamic processes operating in this setting via identification of tectonic rotations related to upper plate extension. Here we present net tectonic rotation results from the Late Cretaceous Mersin ophiolite of southern Turkey that document rapid and progressive rotation of ophiolitic rocks and their associated metamorphic sole. Specifically, we demonstrate that lower crustal cumulate rocks and early dykes intruded into the underlying mantle section have undergone extreme rotation around ridge-parallel, shallowly-plunging axes, consistent with oceanic detachment faulting during spreading. Importantly, later dykes cutting the metamorphic sole experienced rotation around the same axis but with a lower magnitude. We show that these rotations occurred via a common mechanism in a pre-obduction, fore-arc setting, and are best explained by combining (hyper)extension resulting from detachment-mode, amagmatic suprasubduction zone spreading in a fore-arc environment with a recently proposed mechanism for exhumation of metamorphic soles driven by upper plate extension. Available age constraints demonstrate that extreme rotation of these units was accommodated rapidly by these processes over a time period of

The modern physiography of central Turkey is dominated by the 1-km-high Central Anatolian Plateau and the Central Tauride mountains that form the southern plateau margin. These correspond to a Cretaceous–Eocene backarc extensional province and forearc fold-thrust belt, respectively. The extent to which the morphology of the Miocene plateau was inherited from the physiography of the Cretaceous–Eocene subduction zone that assembled the Anatolian crust has not been tested but is important if we are to isolate the signal of Miocene and younger subduction dynamics in the formation of the modern plateau margin. There is no known stratigraphic record of the post-Eocene pre-Miocene evolution of the Taurides. We therefore collected rock samples across the Taurides and used zircon (U-Th)/He (ZHe), apatite (U-Th)/He (AHe), and apatite fission-track (AFT) low-temperature thermochronometers to constrain cooling; we interpret these thermochronometers to signal erosional exhumation. We use inverse thermal modeling to aid interpretation of our results and find that: (1) thermochronometers across the Taurides were reset as a result of heating by the emplacement of the Antalya and Bozkır nappes; (2) AFT and ZHe Eocene cooling ages are related to structurally driven uplift and erosional exhumation on major thrust culminations; (3) dispersed AHe ages record low rates of Oligocene–early Miocene cooling and hence low rates of erosional exhumation; and (4) fast rates of cooling were determined for samples along the margin of the Köprüçay Basin. We interpret that early Miocene cooling is a signal of active erosion of the western Central Taurides at a time of marine sedimentation in the Mut Basin on the southern Central Taurides, and these differing histories may reflect evolution above the Antalya and Cyprus slabs. Our thermochronological data, the enigmatic development of the Antalya Basin, and thrusting within the basin may be explained as the surface expression of stepwise delamination of the Antalya slab from the Tauride hinterland to its current position below the Gulf of Antalya since early Miocene time over a distance of ∼150 km.

The archetypal Semail ophiolite of Oman has inspired much thought on the dynamics of initiation of intra-oceanic subduction zones. Current models invoke subduction initiation at a mid-oceanic ridge located sufficiently close to the Arabian passive margin to allow initiation of continental subduction below the ophiolite within ∼10–15 Myr after the 96–95 Ma age of formation of supra-subduction zone crust. Here, we perform an extensive paleomagnetic analysis of sheeted dyke sections across the Semail ophiolite to restore the orientation of the supra-subduction zone ridge during spreading. Our results consistently indicate that the ridge was oriented NNE–SSW, and we infer that the associated trench, close to the modern obduction front, had the same orientation. Our data are consistent with a previously documented ∼150° clockwise rotation of the ophiolite, and we reconstruct that the original subduction zone was WNW-ward dipping and NNE–SSW striking. Initial subduction likely occurred in the ocean adjacent and parallel to a transform margin of the part of the Arabian continent now underthrust below Iran that originally underpinned the nappes of the Zagros fold-thrust belt. Subduction thus likely initiated along an ancient, continental margin-parallel fracture zone, as also recently inferred from near-coeval ophiolites from the eastern Mediterranean and NW Arabian regions. Subduction initiation was therefore likely induced by (WN)W–(ES)E contraction and this constraint may help the future identification of the dynamic triggers of Neotethyan subduction initiation in the Late Cretaceous.

Key in understanding the geodynamics governing subduction and orogeny is reconstructing the paleogeography of ‘Greater India’ the Indian plate lithosphere that subducted since Tethyan Himalayan continental collision with Asia. Here, we discuss this reconstruction from paleogeographic, kinematic, and geodynamic perspectives and isolate the evolution scenario that is consistent with all three. We follow recent constraints advocating a ~58 Ma initial collision and update a previous kinematic restoration of intra-Asian shortening with a recently proposed model that reconciles long-debated large and small estimates of Indochina extrusion. Our new reconstruction is tested against paleomagnetic data, and against seismic tomographic constraints on paleo-subduction zone locations. The resulting restoration shows ~1000–1200 km of post-collisional intra-Asian shortening, leaving a 2600–3400 km wide Greater India. From a paleogeographic, sediment provenance perspective, Eocene sediments in the Lesser Himalaya and on undeformed India may be derived from Tibet, suggesting that all Greater Indian lithosphere was continental, but may alternatively be sourced from the contemporaneous western Indian orogen unrelated to India-Asia collision. A quantitative kinematic, paleomagnetic perspective prefers major Cretaceous extension and a ‘Greater India Basin’ opening within Greater India, but data uncertainty may speculatively allow for minimal extension. Finally, from a geodynamic perspective, assuming a fully continental Greater India would require that subduction rates close to 20 cm/yr was driven by a down-going lithosphere-crust assemblage more buoyant than the mantle, which seems physically improbable. We conclude that the Greater India Basin scenario is the only sustainable one from all three perspectives. We infer that old pre-collisional lithosphere rapidly entered the lower mantle sustaining high subduction rates, whilst post-collisional continental and young Greater India basin lithosphere did not, inciting the rapid India-Asia convergence deceleration ~8 Myr after collision. Subsequent absolute northward slab migration and overturning caused flat slab subduction, Tibetan shortening, arc migration and arc volume decrease.

Continent-continent collision drives crustal deformation, topographic rise and geodynamic change. Africa-Eurasia convergence accommodated in the Eastern Mediterranean involved subduction of the Neotethyan oceanic lithosphere in Anatolia. Subduction was followed by collision of continental crust of Greater Adria with Eurasia to form the Izmir-Ankara-Erzincan suture zone. Discerning the effects of this collision from pre-collisional ophiolite obduction-related orogeny of Greater Adria is notoriously difficult. Estimates on the timing of collision in Central Anatolia are based on a forearc-to-foreland basin transition along the Eurasian margin and suggest a ~60 Ma age of initial collision. Here, we assess whether this age is also representative for collision in Eastern Anatolia and across the Cenozoic Sivas Basin that straddles the Greater Adria-Europe suture. To this end we retro-deform regional block rotations in the Pontides, the Kirsehir Block and the Taurides, building a first-order regional ‘block circuit’ around the Sivas Basin. We show that up to ~700 km of convergence must have been accommodated across the Sivas Basin after Central Anatolian Kirsehir-Pontide collision at ~60 Ma – an order of magnitude more than estimated crustal shortening, and that wholesale lithospheric subduction must have occurred throughout much of the Cenozoic. Paleocene collision would require that this subduction consumed continental lithosphere, which is unlikely. We consequently infer that oceanic subduction continued much longer in Eastern Anatolia, perhaps well into the Miocene. We postulate that prolonged oceanic subduction and slab pull drew the Eastern Taurides north relative to the Central Taurides, leading to shortening and oroclinal bending in Central Anatolia. The diachronous demise of the Neotethys Ocean in Anatolia, as a function of its paleogeography, is thus a likely driver for the strong non-cylindricity of the Cenozoic Anatolian collisional orogen.

During evolution of the South Sandwich subduction zone, which has consumed South American plate oceanic lithosphere, somehow continental crust of both the South American and Antarctic plates have become incorporated into its upper plate. Continental fragments of both plates are currently separated by small oceanic basins in the upper plate above the South Sandwich subduction zone, in the Scotia Sea region, but how fragments of both continents became incorporated in the same upper plate remains enigmatic. Here we present an updated kinematic reconstruction of the Scotia Sea region using the latest published marine magnetic anomaly constraints, and place this in a South America-Africa-Antarctica plate circuit in which we take intracontinental deformation into account. We show that a change in fracture zone orientation in the Weddell Sea requires that previously inferred initiation of subduction of South American oceanic crust of the northern Weddell below the eastern margin of South Orkney Islands continental crust, then still attached to the Antarctic Peninsula, already occurred around 80 Ma. We propose that subsequently, between ~71-50 Ma, the trench propagated northwards into South America by delamination of South American lithosphere: this resulted in the transfer of delaminated South American continental crust to the overriding plate of the South Sandwich subduction zone. We show continental delamination may have been facilitated by absolute southward motion of South America that was resisted by South Sandwich slab dragging. Pre-drift extension preceding the oceanic Scotia Sea basins led around 50 Ma to opening of the Drake Passage, preconditioning the southern ocean for the Antarctic Circumpolar Current. This 50 Ma extension was concurrent with a strong change in absolute plate motion of the South American Plate that changed from S to WNW, leading to upper plate retreat relative to the more or less mantle stationary South Sandwich Trench that did not partake in the absolute plate motion change. While subduction continued, this mantle-stationary trench setting lasted until ~30 Ma, after which rollback started to contribute to back-arc extension. We find that roll-back and upper plate retreat have contributed more or less equally to the total amount of ~2000 km of extension accommodated in the Scotia Sea basins. We highlight that viewing tectonic motions in a context of absolute plate motion is key for identifying slab motion (e.g. rollback, trench-parallel slab dragging) and consequently mantle-forcing of geological processes.

Subduction zones are unique to Earth and fundamental in its evolution, yet we still know little about the causes and mechanisms of their initiation. Numerical models show that far-field forcing may cause subduction initiation at weak pre-existing structures, while inferences from modern subduction zones suggest initiation through spontaneous lithospheric gravitational collapse. For both endmembers, the timing of subduction inception corresponds with initial lower plate burial, whereas coeval or delayed extension in the upper plate are diagnostic of spontaneous or forced subduction initiation, respectively. In modern systems, the earliest extension-related upper plate rocks are found in forearcs, but lower plate rocks that recorded initial burial have been subducted and are inaccessible. Here, we investigate a fossil system, the archetypal Semail Ophiolite of Oman, which exposes both lower and upper plate relics of incipient subduction stages. We show with Lu–Hf and U–Pb geochronology of the lower and upper plate material that initial burial of the lower plate occurred before 104 million years ago, predating upper plate extension and the formation of Semail oceanic crust by at least 8 Myr. Such a time lag reveals far-field forced subduction initiation and provides unequivocal, direct evidence for a subduction initiation mechanism in the geological record.

The North American Cordillera has been shaped by a long history of accretion of arcs and other buoyant crustal fragments to the western margin of the North American plate since early Mesozoic time. The southernmost accreted terrane is the Guerrero terrane of southwestern Mexico, a latest Jurassic-Cretaceous volcanic arc built on a Triassic accretionary prism. Interpretations of the origin of the Guerrero terrane vary: Some authors consider it a far-traveled, exotic intra-oceanic island arc, while others view it as the (par)autochthonous, extended North American continental margin. We present new paleomagnetic and U-Pb zircon data from Lower Cretaceous sedimentary rocks of the Guerrero terrane. These data show that the Guerrero terrane has a latitudinal plate-motion history equal to that of the North America plate, both before and after accretion. This confirms paleogeographic models in which the Guerrero arc successions formed on North American crust that rifted away from the Mexican mainland by approximately east-west opening of a back-arc basin above an eastward-dipping subduction zone. Additionally, it renders alternative paleogeographical models in which the Guerrero terrane is considered to be exotic to the North American continent unlikely. The phase of back-arc spreading resulted in the short-lived existence of an additional "Guerrero" tectonic plate between the North American and Farallon plates and, upon closure of the back-arc basin, the growth of the North American continent.

To quantitatively reconstruct the kinematic evolution of Central and Eastern Anatolia within the framework of Neotethyan subduction accommodating Africa–Eurasia convergence, we paleomagnetically assess the timing and amount of vertical axis rotations across the Ulukışla and Sivas regions. We show paleomagnetic results from ∼ 30 localities identifying a coherent rotation of a SE Anatolian rotating block comprised of the southern Kırşehir Block, the Ulukışla Basin, the Central and Eastern Taurides, and the southern part of the Sivas Basin. Using our new and published results, we compute an apparent polar wander path (APWP) for this block since the Late Cretaceous, showing that it experienced a ∼ 30–35° counterclockwise vertical axis rotation since the Oligocene time relative to Eurasia. Sediments in the northern Sivas region show clockwise rotations. We use the rotation patterns together with known fault zones to argue that the counterclockwise-rotating domain of south-central Anatolia was bounded by the Savcılı Thrust Zone and Deliler–Tecer Fault Zone in the north and by the African–Arabian trench in the south, the western boundary of which is poorly constrained and requires future study. Our new paleomagnetic constraints provide a key ingredient for future kinematic restorations of the Anatolian tectonic collage.

Across the entire mantle we interpret 94 positive seismic wave-speed anomalies as subducted lithosphere and associate these slabs with their geological record. We document this as the Atlas of the Underworld, also accessible online at www.atlas-of-the-underworld.org, a compilation comprising subduction systems active in the past ~ 300 Myr. Deeper slabs are correlated to older geological records, assuming no relative horizontal motions between adjacent slabs following break-off, using knowledge of global plate circuits, but without assuming a mantle reference frame. The longest actively subducting slabs identified reach the depth of ~ 2500 km and some slabs have impinged on Large Low Shear Velocity Provinces in the deepest mantle. Anomously fast sinking of some slabs occurs in regions affected by long-term plume rising. We conclude that slab remnants eventually sink from the upper mantle to the core-mantle boundary. The range in subduction-age versus – depth in the lower mantle is largely inherited from the upper mantle history of subduction. We find a significant depth variation in average sinking speed of slabs. At the top of the lower mantle average slab sinking speeds are between 10 and 40 mm/yr, followed by a deceleration to 10–15 mm/yr down to depths around 1600–1700 km. In this interval, in situ time-stationary sinking rates suggest deceleration from 20 to 30 mm/yr to 4–8 mm/yr, increasing to 12–15 mm/yr below 2000 km. This corroborates the existence of a slab deceleration zone but we do not observe long-term (> 60 My) slab stagnation, excluding long-term stagnation due to compositional effects. Conversion of slab sinking profiles to viscosity profiles shows the general trend that mantle viscosity increases in the slab deceleration zone below which viscosity slowly decreases in the deep mantle. This is at variance with most published viscosity profiles that are derived from different observations, but agrees qualitatively with recent viscosity profiles suggested from material experiments.

Central Anatolia exposes previously buried and metamorphosed, continent-derived rocks – the Kirsehir and Afyon zones – now covering an area of ∼300 × 400 km. So far, the exhumation history of these rocks has been poorly constrained. We show for the first time that the major, >120 km long, top-NE ‘Ivriz’ Detachment controlled the exhumation of the HP/LT metamorphic Afyon Zone in southern Central Anatolia. We date its activity at between the latest Cretaceous and early Eocene times. Combined with previously documented isolated extensional detachments found in the Kirsehir Block, our results suggest that a major province governed by extensional exhumation was active throughout Central Anatolia between ∼80 and ∼48 Ma. Although similar in dimension to the Aegean extensional province to the east, the Central Anatolian extensional province is considerably older and was controlled by a different extension direction. From this, we infer that the African slab(s) that subducted below Anatolia must have rolled back relative to the Aegean slab since at least the latest Cretaceous, suggesting that these regions were underlain by a segmented slab. Whether or not these early segments already corresponded to the modern Aegean, Antalya, and Cyprus slab segments remains open for debate, but slab segmentation must have occurred much earlier than previously thought.

The plate kinematic history of the South China Sea opening is key to reconstructing how the Mesozoic configuration of Panthalassa and Tethyan subduction systems evolved into today's complex Southeast Asian tectonic collage. The South China Sea is currently flanked by the Palawan Continental Terrane in the south and South China in the north and the two blocks have long been assumed to be conjugate margins. However, the paleogeographic history of the Palawan Continental Terrane remains an issue of uncertainty and controversy, especially regarding the questions of where and when it was separated from South China. Here we employ detrital zircon U–Pb geochronology and heavy mineral analysis on Cretaceous and Eocene strata from the northern South China Sea and Palawan to constrain the Late Mesozoic–Early Cenozoic provenance and paleogeographic evolution of the region testing possible connection between the Palawan Continental Terrane and the northern South China Sea margin. In addition to a revision of the regional stratigraphic framework using the youngest zircon U–Pb ages, these analyses show that while the Upper Cretaceous strata from the Palawan Continental Terrane are characterized by a dominance of zircon with crystallization ages clustering around the Cretaceous, the Eocene strata feature a large range of zircon ages and a new mineral group of rutile, anatase, and monazite. On the one hand, this change of sediment compositions seems to exclude the possibility of a latest Cretaceous drift of the Palawan Continental Terrane in response to the Proto-South China Sea opening as previously inferred. On the other hand, the zircon age signatures of the Cretaceous–Eocene strata from the Palawan Continental Terrane are largely comparable to those of contemporary samples from the northeastern South China Sea region, suggesting a possible conjugate relationship between the Palawan Continental Terrane and the eastern Pearl River Mouth Basin. Thus, the Palawan Continental Terrane is interpreted to have been attached to the South China margin from the Cretaceous until the Oligocene oceanization of the South China Sea. In our preferred paleogeographic scenario, the sediment provenance in the northeastern South China Sea region changed from dominantly nearby Cretaceous continental arcs of the South China margin to more distal southeastern South China in the Eocene.

Formation of new subduction zones represents one of the cornerstones of plate tectonics, yet both the kinematics and geodynamics governing this process remain enigmatic. A major subduction initiation event occurred in the Late Cretaceous, within the Neo-Tethys Ocean between Gondwana and Eurasia. Suprasubduction zone ophiolites (i.e., emerged fragments of ancient oceanic lithosphere formed at suprasubduction spreading centers) were generated during this subduction event and are today distributed in the eastern Mediterranean region along three ~E-W trending ophiolitic belts. Several models have been proposed to explain the formation of these ophiolites and the evolution of the associated intra-Neo-Tethyan subduction zone. Here we present new paleospreading directions from six Upper Cretaceous ophiolites of Turkey, Cyprus, and Syria, calculated by using new and published paleomagnetic data from sheeted dyke complexes. Our results show that ~NNE-SSW subduction zones were formed within the Neo-Tethys during the Late Cretaceous, which we propose were part of a major step-shaped subduction system composed of ~NNE-SSW and ~WNW-ESE segments. We infer that this subduction system developed within old (Triassic?) lithosphere, along fracture zones and perpendicular weakness zones, since the Neo-Tethyan spreading ridge formed during Gondwana fragmentation would have already been subducted at the Pontides subduction zone by the Late Cretaceous. Our new results provide an alternative kinematic model of Cretaceous Neo-Tethyan subduction initiation and call for future research on the mechanisms of subduction inception within old (and cold) lithosphere and the formation of metamorphic soles below suprasubduction zone ophiolites in the absence of nearby spreading ridges.

At two trench segments below the Andes, the Nazca Plate is subducting sub-horizontally over ∼200–300 km, thought to result from a combination of buoyant oceanic-plateau subduction and hydrodynamic mantle-wedge suction. Whether the actual conditions for both processes to work in concert existed is uncertain. Here we infer from a tectonic reconstruction of the Andes constructed in a mantle reference frame that the Nazca slab has retreated at ∼2 cm per year since ∼50 Ma. In the flat slab portions, no rollback has occurred since their formation at ∼12 Ma, generating ‘horse-shoe' slab geometries. We propose that, in concert with other drivers, an overpressured sub-slab mantle supporting the weight of the slab in an advancing upper plate-motion setting can locally impede rollback and maintain flat slabs until slab tearing releases the overpressure. Tear subduction re-establishes a continuous slab and allows the process to recur, providing a mechanism for the transient character of flat slabs., How flat slabs at subduction zones are created remains unclear. Here, the authors show that the Nazca slab has retreated at ∼2 cm per year since 50 Ma but no rollback has occurred in the last ∼12 Myr in the flat slab, implying that an overpressured sub-slab mantle can impede rollback and maintain a flat slab.

The extended northern continental margin of Africa is currently colliding with and under-thrusting northward below Cyprus. The age of onset of this collision is poorly constrained, but is critical if we are to quantify continental subduction and evaluate its role during spectacular recent uplift of southern Turkey from ~7 Ma to the present. Here, we reevaluate the evolution of northern Cyprus and document the modern structure of the Kyrenia fold-thrust belt in a balanced cross section for the first time, to determine the timing and amount of shortening. The belt deformed an Upper Cretaceous to Miocene stratigraphy, which was deposited onto a metamorphic basement. The fold-thrust belt was previously proposed to have formed during two stages: in Eocene and late Miocene time, based on stratigraphic and limited structural evidence. We revaluate evidence for an Eocene phase of thrusting in Kyrenia and find that 1) repetitions of Eocene and older rocks are explained by Miocene thrusting, and newly documented olistoliths were misinterpreted as pre-Miocene thrust slices 2) uplift in Eocene time was not isolated to Kyrenia but also affected northern Arabia and results from regional dynamic topography or forebulge formation. Based on kinematic constraints from a plate reconstruction, we interpret that metamorphic rocks in Kyrenia were metamorphosed by Cretaceous burial below the Troodos ophiolite and were subsequently exhumed by upper plate extension during latest Cretaceous obduction onto the African margin. The latest Cretaceous-Miocene plate boundary between Africa and the Taurides was located to the north of Kyrenia. The Kyrenia fold-thrust belt accommodated a minimum of 17.5 km of shortening as a result of ~9–6 Ma thrusting as a result of initial Africa-Taurides collision. The northern edge of African continental crust arrived below the Taurides, where it has been seismologically imaged, in late Pleistocene time and may have caused or contributed to recently recognised major Tauride uplift since ~0.5 Ma.

Reconstruction of how the Pacific Plate originated at a point improves plate reconstructions of the vast ocean around Pangea., The oceanic Pacific Plate started forming in Early Jurassic time within the vast Panthalassa Ocean that surrounded the supercontinent Pangea, and contains the oldest lithosphere that can directly constrain the geodynamic history of the circum-Pangean Earth. We show that the geometry of the oldest marine magnetic anomalies of the Pacific Plate attests to a unique plate kinematic event that sparked the plate’s birth at virtually a point location, surrounded by the Izanagi, Farallon, and Phoenix Plates. We reconstruct the unstable triple junction that caused the plate reorganization, which led to the birth of the Pacific Plate, and present a model of the plate tectonic configuration that preconditioned this event. We show that a stable but migrating triple junction involving the gradual cessation of intraoceanic Panthalassa subduction culminated in the formation of an unstable transform-transform-transform triple junction. The consequent plate boundary reorganization resulted in the formation of a stable triangular three-ridge system from which the nascent Pacific Plate expanded. We link the birth of the Pacific Plate to the regional termination of intra-Panthalassa subduction. Remnants thereof have been identified in the deep lower mantle of which the locations may provide paleolongitudinal control on the absolute location of the early Pacific Plate. Our results constitute an essential step in unraveling the plate tectonic evolution of “Thalassa Incognita” that comprises the comprehensive Panthalassa Ocean surrounding Pangea.

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.Tonga‐Kermadec subduction zone is likely only ~30 Ma oldNorthward Australian absolute plate motion is illustrated by 1,500 km offset between South Loyalty Basin slab and New Caledonian sutureTonga‐Kermadec slab was dragged 1,200 km northward through the mantle, alongside eastward rollback.

Jurassic subduction initiation in the Neo-Tethys Ocean eventually led to the collision of the Adria-Africa and Eurasia continents and the formation of an ~6,000 km long Alpine orogen spanning from Iberia to Iran. Reconstructing the location and geometry of the plate boundaries of the now disappeared Neo-Tethys during the initial moments of its closure is instrumental to perform more realistic plate reconstructions of this region, of ancient ocean basins in general, and on the process of subduction initiation. Neo-Tethyan relics are preserved in an ophiolite belt distributed above the Dinaric-Hellenic fold-thrust belt. Here we provide the first quantitative constraints on the geometry of the spreading ridges and trenches active in the Jurassic Neo-Tethys using a paleomagnetically based net tectonic rotation analysis of sheeted dykes and dykes from the West and East Vardar Ophiolites of Serbia (Maljen and Ibar) and Greece (Othris, Pindos, Vourinos, and Guevgueli). Based on our results and existing geological evidence, we show that initial Middle Jurassic (~175 Ma) closure of the western Neo-Tethys was accommodated at a N-S trending, west dipping subduction zone initiated near and parallel to the spreading ridge. The West Vardar Ophiolites formed in the forearc parallel to this new trench. Simultaneously, the East Vardar Ophiolites formed above a second N-S to NW-SE trending subduction zone located close to the European passive margin. We tentatively propose that this second subduction zone had been active since at least the Middle Triassic, simultaneously accommodating the closure of the Paleo-Tethys and the back-arc opening of Neo-Tethys.

Models of arc-continent accretion often assume that the period of subduction of continental lithosphere before plate boundary reorganization is fairly short lived, yet the timescale of this period is poorly constrained by observations in the geologic record. The island of Timor is the uplifted accretionary complex resulting from the active collision of the Banda volcanic arc with the Australian continental margin. The exposure of underplated and exhumed Australian strata on Timor allows for the characterization of the structural history of accretion of uppermost Australian crust and the quantification of subduction of its original continental lithospheric underpinnings. New structural mapping in East Timor (Timor-Leste) reveals that duplexing of a 2-km-thick package of Australian continental strata has built the majority of the structural elevation of the Timor orogen. Coupling new structural observations with previous thermochronology results reveals the sequence of deformation within the orogen, the presence of subsurface duplexing below the hinterland slate belt, and motion along a foreland subsurface thrust ramp. Construction of balanced cross sections allows for the quantification of the amount of shortening in the orogen, and from that, the length of the subducted Australian continental lithosphere. Two balanced cross sections in East Timor reveal 326–362 km of shortening and that 215–229 km of Australian continental lithosphere have been subducted below the Banda forearc. These results highlight the fact that considerable amounts of continental lithosphere can be subducted while accreting only a thin section of uppermost crust. Continental subduction may have been favorable at Timor because of fast subduction rates, old oceanic crust at the consumed Australian margin, and subduction of some length of transitional crust. These results provide quantitative constraints for future numerical modeling of the geodynamics of continental subduction and arc-continent collision.