Kinematic Boundary Conditions Favouring Subduction Initiation at Passive Margins Over Subduction at Mid-oceanic Ridges

We perform numerical modelling to simulate the shortening of an oceanic basin and the adjacent continental margins in order to discuss the relationship between compressional stresses acting on the lithosphere and the time dependent strength of the mid-oceanic ridges within the frame of subduction initiation. We focus on the role of stress regulating mechanisms by testing the stress–strain-rate response to convergence rate, and the thermo-tectonic age of oceanic and continental lithospheres. We find that, upon compression, subduction initiation at passive margin is favoured for thermally thin (Palaeozoic or younger) continental lithospheres (<160 km) over cratons (>180 km), and for oceanic basins younger than 60 Myr (after rifting). The results also highlight the importance of convergence rate that controls stress distribution and magnitudes in the oceanic lithosphere. Slow convergence (<0.9 cm/yr) favours strengthening of the ridge and build-up of stress at the ocean-continent transition allowing for subduction initiation at passive margins over subduction at mid-oceanic ridges. The results allow for identifying geodynamic processes that fit conditions for subduction nucleation at passive margins, which is relevant for the unique case of the Alps. We speculate that the slow Africa–Europe convergence between 130 and 85 Ma contributes to the strengthening of the mid-oceanic ridge, leading to subduction initiation at passive margin 60–70 Myr after rifting and passive margin formation.

Hemispheric geochemical dichotomy of the mantle is a legacy of austral supercontinent assembly and deep subduction of continental crust

Oceanic hotspots with extreme enriched mantle radiogenic isotopic signatures—including high 87Sr/86Sr and low 143Nd/144Nd indicative of ancient subduction of continental crust—are restricted to the southern hemispheric mantle. However, the mechanisms responsible for concentrating subducted continental crust in the austral mantle are unknown. We show subduction of sediments and subduction eroded material, and lower continental crust delamination, cannot generate this spatially coherent austral domain. However, late Neoproterozoic to Paleozoic continental collisions—associated with the assembly of Gondwana and Pangea—were positioned predominantly in the southern hemisphere during the late Neoproterozoic appearance of widespread continental ultra-high-pressure (UHP, >2.7 gigapascals) metamorphic terranes, which marked the onset of deep subduction of upper continental crust. We propose that deep subduction of upper continental crust at ancient rifted-passive margins during austral supercontinent assembly, from 650-300 Ma, resulted in enhanced upper continental crust delivery into the southern hemisphere mantle. In contrast, EM domains are absent in boreal hotspots, for two reasons. First, continental crust subducted after 300—when the continents drifted into the northern hemisphere—has had insufficient time to return to the surface in plumes feeding northern hemisphere hotspots. Second, before the appearance of continental UHP rocks at 650 Ma, upper continental crust was not subducted to great depths, thus precluding its subduction into the northern hemisphere mantle during the Precambrian when continents may have been located in the northern hemisphere. Our model implies a recent formation of the austral EM domain, explains the geochemical dichotomy between austral and boreal hotspots, and may explain why austral hotspots outnumber boreal hotspots.

Evidence of oceanic plate delamination in the Northern Atlantic

The Earth’s surface is constantly being recycled by plate tectonics. Subduction of oceanic lithosphere and delamination of continental lithosphere constitute the two most important mechanisms by which the Earth’s lithosphere is recycled into the mantle. Delamination or detachment in continental regions typically occurs below mountain belts due to a weight excess of overthickened lithospheric mantle, which detaches from overlying lighter crust, aided by the existence of weak layers within the continental lithosphere. Oceanic lithosphere is classically pictured as a rigid plate with a strong core that does not allow for delamination to occur. Here, we propose that active delamination of oceanic lithosphere occurs offshore Southwest Iberia. The process is assisted by the existence of a lithospheric serpentinized layer that allows the lower part of the lithosphere to decouple from the overlying crust. Tomography images reveal a sub-lithospheric high-velocity anomaly below this region, which we interpret as a delaminating block of old oceanic lithosphere. We present numerical models showing that for a geological setting mimicking offshore Southwest Iberia delamination of oceanic lithosphere is possible and may herald subduction initiation, which is a long-unsolved problem in the theory of plate tectonics. We further propose that such oceanic delamination is responsible for the highest-magnitude earthquakes in Europe, including the M8.5-8.7 Great Lisbon Earthquake of 1755 and the M7.9 San Vincente earthquake of 1969. In particular, our numerical models, in combination with calculations on seismic potential, provide a solution for the instrumentally recorded 1969 event below the flat Horseshoe abyssal plain, away from mapped tectonics faults. Delamination of old oceanic lithosphere near passive margins constitutes a new class of subduction initiation mechanisms, with fundamental implications for the dynamics of the Wilson cycle.

Escarpment retreat rates derived from detrital cosmogenic nuclide concentrations

High-relief great escarpments at passive margins present a paradoxical combination of high-relief topography but low erosion rates suggesting low rates of landscape change. However, vertical erosion rates do not offer a straightforward metric of horizontal escarpment retreat rates, so we attempt to address this problem in this paper. We show that detrital cosmogenic nuclide concentrations can be interpreted as a directionally dependent mass flux to characterize patterns of non-vertical landscape evolution, e.g., an escarpment characterized by horizontal retreat. We present two methods for converting cosmogenic nuclide concentrations into escarpment retreat rates and calculate the retreat rates of escarpments with published cosmogenic 10Be concentrations from the Western Ghats of India. Escarpment retreat rates of the Western Ghats inferred from this study vary within a range of hundreds to thousands of meters per Myr. We show that the current position and morphology of the Western Ghats are consistent with an escarpment retreating at a near-constant rate from the coastline since rifting.

The closure of the Vardar ocean (the western domain of the northern Neotethys) from early Middle Jurassic to Paleocene time, based on surface geology of eastern Pelagonia and the Vardar zone, biostratigraphy, and seismic-tomographic images of the mantle below the Central Hellenides

Seismic tomographic images of the mantle below the Hellenides indicate that the Vardar ocean probably had a composite width of over 3000 kilometres. From surface geology we know that this ocean was initially located between two passive margins: Pelagonian Adria in the west and Serbo-Macedonian-Eurasia in the east. Pelagonia was covered by a carbonate platform that accumulated, during Late Triassic to Early Cretaceous time, where highly diversified carbonate sedimentary environments evolved and reacted to the adjacent, converging Vardar ocean plate. We conceive that on the east side of the Vardar ocean, a Cretaceous carbonate platform evolved from Aptian to Maastrichtian time in the forearc basin of the Vardar supra-subduction volcanic arc complex. The closure of the Vardar ocean occurred in one episode of ophiolite obduction and in two episodes of intra-oceanic subduction.