Dynamic slab segmentation due to brittle-ductile damage in the outer rise.

Subduction is the major plate driving force, and the strength of the subducting plate controls many aspects of the thermochemical evolution of Earth. Each subducting plate experiences intense normal faulting1-9 during bending that accommodates the transition from horizontal to downwards motion at the outer rise at trenches. Here we investigate the consequences of this bending-induced plate damage using numerical subduction models in which both brittle and ductile deformation, including grain damage, are tracked and coupled self-consistently. Pervasive slab weakening and pronounced segmentation can occur at the outer-rise region owing to the strong feedback between brittle and ductile damage localization. This slab-damage phenomenon explains the subduction dichotomy of strong plates and weak slabs10, the development of large-offset normal faults6,7 near trenches, the occurrence of segmented seismic velocity anomalies11 and distinct interfaces imaged within subducted slabs12,13, and the appearance of deep, localized intraplate areas of reduced effective viscosity14 observed at trenches. Furthermore, brittle-viscously damaged slabs show a tendency for detachment at elevated mantle temperatures. Given Earth's planetary cooling history15, this implies that intermittent subduction with frequent slab break-off episodes16 may have been characteristic for Earth until more recent times than previously suggested17.

Depletion of the upper mantle by convergent tectonics in the Early Earth.

Partial melting of mantle peridotites at spreading ridges is a continuous global process that forms the oceanic crust and refractory, positively buoyant residues (melt-depleted mantle peridotites). In the modern Earth, these rocks enter subduction zones as part of the oceanic lithosphere. However, in the early Earth, the melt-depleted peridotites were 2-3 times more voluminous and their role in controlling subduction regimes and the composition of the upper mantle remains poorly constrained. Here, we investigate styles of lithospheric tectonics, and related dynamics of the depleted mantle, using 2-D geodynamic models of converging oceanic plates over the range of mantle potential temperatures (Tp = 1300-1550 °C, ∆T = T - Tmodern = 0-250 °C) from the Archean to the present. Numerical modeling using prescribed plate convergence rates reveals that oceanic subduction can operate over this whole range of temperatures but changes from a two-sided regime at ∆T = 250 °C to one-sided at lower mantle temperatures. Two-sided subduction creates V-shaped accretionary terrains up to 180 km thick, composed mainly of highly hydrated metabasic rocks of the subducted oceanic crust, decoupled from the mantle. Partial melting of the metabasic rocks and related formation of sodic granitoids (Tonalite-Trondhjemite-Granodiorite suites, TTGs) does not occur until subduction ceases. In contrast, one sided-subduction leads to volcanic arcs with or without back-arc basins. Both subduction regimes produce over-thickened depleted upper mantle that cannot subduct and thus delaminates from the slab and accumulates under the oceanic lithosphere. The higher the mantle temperature, the larger the volume of depleted peridotites stored in the upper mantle. Extrapolation of the modeling results reveals that oceanic plate convergence at ∆T = 200-250 °C might create depleted peridotites (melt extraction of > 20%) constituting more than half of the upper mantle over relatively short geological times (~ 100-200 million years). This contrasts with the modeling results at modern mantle temperatures, where the amount of depleted peridotites in the upper mantle does not increase significantly with time. We therefore suggest that the bulk chemical composition of upper mantle in the Archean was much more depleted than the present mantle, which is consistent with the composition of the most ancient lithospheric mantle preserved in cratonic keels.

Building cratonic keels in Precambrian plate tectonics.

The ancient cores of continents (cratons) are underlain by mantle keels-volumes of melt-depleted, mechanically resistant, buoyant and diamondiferous mantle up to 350 kilometres thick, which have remained isolated from the hotter and denser convecting mantle for more than two billion years. Mantle keels formed only in the Early Earth (approximately 1.5 to 3.5 billion years ago in the Precambrian eon); they have no modern analogues1-4. Many keels show layering in terms of degree of melt depletion5-7. The origin of such layered lithosphere remains unknown and may be indicative of a global tectonics mode (plate rather than plume tectonics) operating in the Early Earth. Here we investigate the possible origin of mantle keels using models of oceanic subduction followed by arc-continent collision at increased mantle temperatures (150-250 degrees Celsius higher than the present-day values). We demonstrate that after Archaean plate tectonics began, the hot, ductile, positively buoyant, melt-depleted sublithospheric mantle layer located under subducting oceanic plates was unable to subduct together with the slab. The moving slab left behind craton-scale emplacements of viscous protokeel beneath adjacent continental domains. Estimates of the thickness of this sublithospheric depleted mantle show that this mechanism was efficient at the time of the major statistical maxima of cratonic lithosphere ages. Subsequent conductive cooling of these protokeels would produce mantle keels with their low modern temperatures, which are suitable for diamond formation. Precambrian subduction of oceanic plates with highly depleted mantle is thus a prerequisite for the formation of thick layered lithosphere under the continents, which permitted their longevity and survival in subsequent plate tectonic processes.

Horizontal mantle flow controls subduction dynamics.

It is generally accepted that subduction is driven by downgoing-plate negative buoyancy. Yet plate age -the main control on buoyancy- exhibits little correlation with most of the present-day subduction velocities and slab dips. "West"-directed subduction zones are on average steeper (~65°) than "East"-directed (~27°). Also, a "westerly"-directed net rotation of the lithosphere relative to the mantle has been detected in the hotspot reference frame. Thus, the existence of an "easterly"-directed horizontal mantle wind could explain this subduction asymmetry, favouring steepening or lifting of slab dip angles. Here we test this hypothesis using high-resolution two-dimensional numerical thermomechanical models of oceanic plate subduction interacting with a mantle flow. Results show that when subduction polarity is opposite to that of the mantle flow, the descending slab dips subvertically and the hinge retreats, thus leading to the development of a back-arc basin. In contrast, concordance between mantle flow and subduction polarity results in shallow dipping subduction, hinge advance and pronounced topography of the overriding plate, regardless of their age-dependent negative buoyancy. Our results are consistent with seismicity data and tomographic images of subduction zones. Thus, our models may explain why subduction asymmetry is a common feature of convergent margins on Earth.

Plate tectonics on the Earth triggered by plume-induced subduction initiation.

Scientific theories of how subduction and plate tectonics began on Earth--and what the tectonic structure of Earth was before this--remain enigmatic and contentious. Understanding viable scenarios for the onset of subduction and plate tectonics is hampered by the fact that subduction initiation processes must have been markedly different before the onset of global plate tectonics because most present-day subduction initiation mechanisms require acting plate forces and existing zones of lithospheric weakness, which are both consequences of plate tectonics. However, plume-induced subduction initiation could have started the first subduction zone without the help of plate tectonics. Here, we test this mechanism using high-resolution three-dimensional numerical thermomechanical modelling. We demonstrate that three key physical factors combine to trigger self-sustained subduction: (1) a strong, negatively buoyant oceanic lithosphere; (2) focused magmatic weakening and thinning of lithosphere above the plume; and (3) lubrication of the slab interface by hydrated crust. We also show that plume-induced subduction could only have been feasible in the hotter early Earth for old oceanic plates. In contrast, younger plates favoured episodic lithospheric drips rather than self-sustained subduction and global plate tectonics.