Tectonic evolution and deep mantle structure of the eastern Tethys since the latest Jurassic

Abstract

The breakup of Pangea in the Jurassic saw the opening of major ocean basins at the expense of older Tethyan and Pacific oceanic plates. Although the Tethyan seafloor spreading history has been lost to subduction, proxy indicators from multiple generations of Tethyan ribbon terranes, as well as the active margin geological histories of volcanism and ophiolite obduction events can be used to reconstruct these ancient oceanic plates. The plate reconstructions presented in this study reconcile observations from ocean basins and the onshore geological record to provide a regional synthesis, embedded in a global plate motion model, of the IndiaEurasia convergence history, the accretionary growth of Southeast Asia and the Tethyan-Pacific tectonic link through the New Guinea margin. The global plate motion model presented in this study captures the timedependent evolution of plates and their tectonic boundaries since 160 Ma, which are assimilated as surface boundary conditions for numerical experiments of mantle convection. We evaluate subducted slab locations and geometries predicted by forward mantle flow models against P-and S-wave seismic tomography models. This approach harnesses modern plate reconstruction techniques, mantle convection models with imposed one-sided subduction, and constraints from the surface geology to address a number of unresolved Tethyan geodynamic controversies. Our synthesis reveals that north-dipping subduction beneath Eurasia in the latest Jurassic consumed the Meso-Tethys, and suggests that northward slab pull opened the younger Neo-Tethyan ocean basin from ~ 155 Ma. We model the rifting of ' Argoland', representing the East Java and West Sulawesi continental fragments, as a northward transfer of continental terranes in the latest Jurassic from the northwest Australian shelf -likely colliding first with parts of the Woyla intra-oceanic arc in the mid-Cretaceous, and accreting to the Borneo (Sundaland) core by ~ 80 Ma. The Neo-Tethyan ridge was likely consumed along an intra-oceanic subduction zone south of Eurasia from ~ 105 Ma, leading to a major change in the motion of the Indian Plate by ~ 100 Ma, as observed in the Wharton Basin fracture zone bends. We investigate the geodynamic consequences of long-lived intra-oceanic subduction within the Neo-Tethys, requiring a twostage India-Eurasia collision involving first contact between Greater India and the Kohistan-Ladakh Arc sometime between ~ 60 and 50 Ma, followed by continent-continent collision from ~ 47 Ma. Our models suggest that the Sunda slab kink beneath northwest Sumatra in the mantle transition zone results from the rotation and extrusion of Indochina from ~ 30 Ma. Our results are also the first to reproduce the enigmatic Proto South China Sea slab beneath northern Borneo, as well as the Tethyan/Woyla slab that is predicted at mid-mantle depths south of Sumatra. Further east, our revised reconstructions of the New Guinea margin, notably the evolution of the Sepik composite terrane and the Maramuni subduction zone, produce...

Slab Horizontal Subduction and Slab Tearing Beneath East Asia

“…We modified a recently published global reconstruction (Zahirovic et al, ) with GPlates (http://www.gplates.org), which incorporated compiled global seafloor‐spreading histories with a hybrid reference framework to tie the global tectonic system with an underlying mantle. Our reconstruction combined geological and geophysical observations from East Asia margin (supporting information Table S1) to fully depict the evolving plate boundaries and plate kinematics.…”

“…Similar to the boundary conditions, the initial temperature field was reconstruction‐based (see Supporting Information S1 for more details). All the implementations ensured a general agreement between the predicted buoyancy distribution in the mantle and geophysical observations and helped build a fundamental platform to assess alternative plate models and associated geodynamic scenarios (Zahirovic et al, ).…”

“…Structural map of East Asia. The seafloor crust ages of oceanic basins are rendered according to a global age model (Zahirovic et al, ). The details of the Philippine Sea slab that is subducting beneath SW Japan were adopted from Asamori and Zhao () and are colored by depth.…”

SubMachine: Web‐Based Tools for Exploring Seismic Tomography and Other Models of Earth's Deep Interior

“…This currently include models of the geoid (Pavlis et al, ), free‐air gravity (Sandwell et al, ), vertical gravity gradient (Sandwell et al, ), and topography and bathymetry (Olson et al, ), all in adjustable spherical harmonic degrees. Also included are global crustal model CRUST1.0 (Laske et al, ), normal mode observations and their sensitivities for a selection of upper mantle and lower mantle modes (Koelemeijer, ; Koelemeijer et al, ), and a suite of time‐dependent plate reconstruction models (Matthews et al, ; Seton et al, ; Zahirovic et al, ). Section 4 raises discussion points about SubMachine's current performance and future directions.…”

Arc volcanism, carbonate platform evolution and palaeo-atmospheric CO<sub>2</sub>: Components and interactions in the deep carbon cycle

“…According to the Accumulation model, the dramatic change that characterised the reduction in CIA 15 lengths ~75 Ma included the termination of a ~3 000-km-long section of a north-dipping intra-oceanic subduction zone in the Northern Tethys ocean, corresponding to the Trans-Tethyan Subduction System (Jagoutz et al, 2016), and the cessation of subduction along the Sundaland margin (McCourt et al, 1996;Zahirovic et al, 2016) (Figs. 5c, 5d).…”