Indian plate paleogeography, subduction and horizontal underthrusting below Tibet: paradoxes, controversies and opportunities.

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.

Correction: Eocene intra-plate shortening responsible for the rise of a faunal pathway in the northeastern Caribbean realm.

[This corrects the article DOI: 10.1371/journal.pone.0241000.].

Eocene intra-plate shortening responsible for the rise of a faunal pathway in the northeastern Caribbean realm.

Intriguing latest Eocene land-faunal dispersals between South America and the Greater Antilles (northern Caribbean) has inspired the hypothesis of the GAARlandia (Greater Antilles Aves Ridge) land bridge. This landbridge, however, should have crossed the Caribbean oceanic plate, and the geological evolution of its rise and demise, or its geodynamic forcing, remain unknown. Here we present the results of a land-sea survey from the northeast Caribbean plate, combined with chronostratigraphic data, revealing a regional episode of mid to late Eocene, trench-normal, E-W shortening and crustal thickening by ∼25%. This shortening led to a regional late Eocene-early Oligocene hiatus in the sedimentary record revealing the location of an emerged land (the Greater Antilles-Northern Lesser Antilles, or GrANoLA, landmass), consistent with the GAARlandia hypothesis. Subsequent submergence is explained by combined trench-parallel extension and thermal relaxation following a shift of arc magmatism, expressed by a regional early Miocene transgression. We tentatively link the NE Caribbean intra-plate shortening to a well-known absolute and relative North American and Caribbean plate motion change, which may provide focus for the search of the remaining connection between 'GrANoLA' land and South America, through the Aves Ridge or Lesser Antilles island arc. Our study highlights the how regional geodynamic evolution may have driven paleogeographic change that is still reflected in current biology.

First Balanced Cross Section Across the Taurides Fold-Thrust Belt: Geological Constraints on the Subduction History of the Antalya Slab in Southern Anatolia.

Eastern Mediterranean subduction accommodated Africa-Eurasia convergence since Mesozoic time and produced multiple subducted slab fragments in the mantle below Anatolia. These included the north dipping Cyprus and ENE-dipping Antalya slabs, which are currently separated by an upper mantle slab gap. Segmentation of these slabs, and associated mantle flow, may have contributed to <8 Ma uplift of the Central Anatolian Plateau. The western Central Taurides fold-thrust belt in southern Turkey is in the upper plate above the Antalya slab and contains a geological record of its subduction. We present the first orogen-scale balanced cross section of the Taurides and find that it formed in two stages: (1) Cretaceous to middle Eocene thrusting resulted in a minimum of 73-km shortening, and (2) Mio-Pliocene thrusting resulted in a minimum of 17.5-km shortening. Eocene shortening accounts for only ~5 Myr of Africa-Eurasia plate convergence. It is unlikely that >400 km of post to middle Eocene plate convergence was accommodated between the Taurides and its Beydaglari platform foreland and instead must have been accommodated south of Beydaglari. The associated southward plate boundary jump separated the Antalya slab from the African plate and the Cyprus slab. The isolated Antalya slab was left in an intraplate setting and is probably still attached to Beydaglari today. We suggest the continental composition of the Antalya slab may have prevented its detachment. Finally, the gap between the Antalya and Cyprus slabs existed since at least Eocene time; their decoupling likely did not contribute to late Neogene Central Anatolian Plateau uplift.

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.

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

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

South-American plate advance and forced Andean trench retreat as drivers for transient flat subduction episodes.

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.

On the enigmatic birth of the Pacific Plate within the Panthalassa Ocean.

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.

A Paleolatitude Calculator for Paleoclimate Studies.

Realistic appraisal of paleoclimatic information obtained from a particular location requires accurate knowledge of its paleolatitude defined relative to the Earth's spin-axis. This is crucial to, among others, correctly assess the amount of solar energy received at a location at the moment of sediment deposition. The paleolatitude of an arbitrary location can in principle be reconstructed from tectonic plate reconstructions that (1) restore the relative motions between plates based on (marine) magnetic anomalies, and (2) reconstruct all plates relative to the spin axis using a paleomagnetic reference frame based on a global apparent polar wander path. Whereas many studies do employ high-quality relative plate reconstructions, the necessity of using a paleomagnetic reference frame for climate studies rather than a mantle reference frame appears under-appreciated. In this paper, we briefly summarize the theory of plate tectonic reconstructions and their reference frames tailored towards applications of paleoclimate reconstruction, and show that using a mantle reference frame, which defines plate positions relative to the mantle, instead of a paleomagnetic reference frame may introduce errors in paleolatitude of more than 15° (>1500 km). This is because mantle reference frames cannot constrain, or are specifically corrected for the effects of true polar wander. We used the latest, state-of-the-art plate reconstructions to build a global plate circuit, and developed an online, user-friendly paleolatitude calculator for the last 200 million years by placing this plate circuit in three widely used global apparent polar wander paths. As a novelty, this calculator adds error bars to paleolatitude estimates that can be incorporated in climate modeling. The calculator is available at www.paleolatitude.org. We illustrate the use of the paleolatitude calculator by showing how an apparent wide spread in Eocene sea surface temperatures of southern high latitudes may be in part explained by a much wider paleolatitudinal distribution of sites than previously assumed.

Sedimentary geology of the middle Carboniferous of the Donbas region (Dniepr-Donets Basin, Ukraine).

The Paleozoic Dniepr-Donets Basin in Belarus, Ukraine, and Russia forms a major hydrocarbon province. Although well- and seismic data have established a 20 km thick stratigraphy, field-studies of its sediments are scarce. The inverted Donbas segment (Ukraine) exposes the middle Carboniferous part of the basin's stratigraphy. Here, we provide detailed sedimentological data from 13 sections that cover 1.5 of the total of 5 km of the Bashkirian and Moscovian stages and assess the paleoenvironment and paleo-current directions. Middle Carboniferous deposition occurred in a shelf environment, with coal deposition, subordinate fluvial facies, and abundant lower and middle shoreface facies, comprising an intercalated package of potential source and reservoir rocks. Sedimentary facies indicate a paleodepth range from below storm wave base to near-coastal swamp environments. Sedimentation and subsidence were hence in pace, with subtle facies changes likely representing relative sea-level changes. Paleocurrent directions are remarkably consistently southeastward in time and space in the different sedimentary facies across the Donbas Fold Belt, illustrating a dominant sedimentary infill along the basin axis, with little basin margin influence. This suggests that the middle Carboniferous stratigraphy of the Dniepr-Donets basin to the northwest probably contains significant amounts of fluvial sandstones, important for assessing hydrocarbon reservoir potential.

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.

Direct terrestrial-marine correlation demonstrates surprisingly late onset of the last interglacial in central Europe.

An interdisciplinary study of a small sedimentary basin at Neumark Nord 2 (NN2), Germany, has yielded a high-resolution record of the palaeomagnetic Blake Event, which we are able to place at the early part of the last interglacial pollen sequence documented from the same section. We use this data to calculate the duration of this stratigraphically important event at 3400 ± 350 yr. More importantly, the Neumark Nord 2 data enables precise terrestrial-marine correlation for the Eemian stage in central Europe. This shows a remarkably large time lag of ca. 5000 yr between the MIS 5e 'peak' in the marine record and the start of the last interglacial in this region.