The importance of continents, oceans and plate tectonics for the evolution of complex life: implications for finding extraterrestrial civilizations.

Within the uncertainties of involved astronomical and biological parameters, the Drake Equation typically predicts that there should be many exoplanets in our galaxy hosting active, communicative civilizations (ACCs). These optimistic calculations are however not supported by evidence, which is often referred to as the Fermi Paradox. Here, we elaborate on this long-standing enigma by showing the importance of planetary tectonic style for biological evolution. We summarize growing evidence that a prolonged transition from Mesoproterozoic active single lid tectonics (1.6 to 1.0 Ga) to modern plate tectonics occurred in the Neoproterozoic Era (1.0 to 0.541 Ga), which dramatically accelerated emergence and evolution of complex species. We further suggest that both continents and oceans are required for ACCs because early evolution of simple life must happen in water but late evolution of advanced life capable of creating technology must happen on land. We resolve the Fermi Paradox (1) by adding two additional terms to the Drake Equation: foc (the fraction of habitable exoplanets with significant continents and oceans) and fpt (the fraction of habitable exoplanets with significant continents and oceans that have had plate tectonics operating for at least 0.5 Ga); and (2) by demonstrating that the product of foc and fpt is very small (< 0.00003-0.002). We propose that the lack of evidence for ACCs reflects the scarcity of long-lived plate tectonics and/or continents and oceans on exoplanets with primitive life.

Speed of thermal adaptation of terrestrial vegetation alters Earth's long-term climate.

Earth's long-term climate is driven by the cycling of carbon between geologic reservoirs and the atmosphere-ocean system. Our understanding of carbon-climate regulation remains incomplete, with large discrepancies remaining between biogeochemical model predictions and the geologic record. Here, we evaluate the importance of the continuous biological climate adaptation of vegetation as a regulation mechanism in the geologic carbon cycle since the establishment of forest ecosystems. Using a model, we show that the vegetation's speed of adaptation to temperature changes through eco-evolutionary processes can strongly influence global rates of organic carbon burial and silicate weathering. Considering a limited thermal adaptation capacity of the vegetation results in a closer balance of reconstructed carbon fluxes into and out of the atmosphere-ocean system, which is a prerequisite to maintain habitable conditions on Earth's surface on a multimillion-year timescale. We conclude that the long-term carbon-climate system is more sensitive to biological dynamics than previously expected, which may help to explain large shifts in Phanerozoic climate.

Variable plate kinematics promotes changes in back-arc deformation regime along the north-eastern Eurasia plate boundary.

The stretching of the lithosphere leading to back-arc basins formation generally develops behind arc-trench systems and is considered the consequence of slab retreat relative to the upper plate. Here, we examine the deformation regime evolution within the overriding plate due to subduction processes, using thermo-mechanical numerical simulations. We explore the north-eastern Eurasia plate boundary and the mechanisms of subducting Pacific plate since 57 Ma. During this time interval, several extensional basins formed along the Eurasia margin, such as the East China Sea, the Japan Sea, and the Kuril basin. Here, we increased the simulation complexity, with the inclusion of (i) the kinematic variability of the Pacific plate over the geological past with respect to a fixed Eurasia, incorporating time-dependent (i.e., temporally evolving) velocities computed from plate motion reconstructions; (ii) a Low-Velocity Zone within the asthenosphere, and (iii) a horizontal eastward mantle flow. Our results show a crucial role of the mantle flow for the development of lithospheric extension and back-arc basin opening, and a main kinematic control of the subduction trench position, which advances and retreats, into distance intervals in the order of ∼  100 km, and providing stages of compression and extension in a back-arc basin.

Onset of double subduction controls plate motion reorganisation.

Face-to-face double subduction systems, in which two oceanic plates subduct toward each other, are essential elements of plate tectonics. Two subduction zones in such systems are typically uneven in age and their spatially and temporally variable dynamics remain enigmatic. Here, with 2D numerical modelling, we demonstrate that the onset of the younger subduction zone strongly changes the dynamics of the older subduction zone. The waxing younger subduction may gradually absorb plate convergence from the older one, resulting in older subduction waning featured by the dramatic decrease in subduction rate and trench retreat. The dynamical transformation of subduction predominance alters the intraplate stress and mantle flow, regulating the relative motion among the three different plates. The process of waxing and waning of subduction zones controls plate motion reorganisation, providing a reference to interpret the past, present, and future evolution of several key double subduction regions found on the modern Earth.

Vertical tearing of subducting plates controlled by geometry and rheology of oceanic plates.

Lateral non-uniform subduction is impacted by continuous plate segmentation owing to vertical tearing of the subducting plate. However, the dynamics and physical controls of vertical tearing remain controversial. Here, we employed 3D numerical models to investigate the effects of trench geometry (offset by a transform boundary) and plate rheology (plate age and the magnitude of brittle/plastic strain weakening) on the evolution of shear stress-controlled vertical tearing within a homogenous subducting oceanic plate. Numerical results suggest that the trench offset geometry could result in self-sustained vertical tearing as a narrow shear zone within the intact subducting oceanic plate, and that this process of tearing could operate throughout the entire subduction process. Further, the critical trench offset length for the maturation of vertical tearing is impacted by plate rheology. Comparison between numerical modelling results and natural observations suggests that vertical tearing attributed to trench offset geometry is broadly developed in modern subduction and collision systems worldwide.

Lithospheric delamination as the driving mechanism of intermediate-depth seismicity in the Bucaramanga Nest, Colombia.

The Bucaramanga nest (BN) is an area of exceptionally strong intermediate-depth seismicity localized in a narrow zone at 150-170 km depth beneath the continental plate in Colombia. To explain the very unusual mantle seismicity cluster in this area, we built a seismic velocity model in the vicinity of BN with the use of local earthquake tomography and developed a numerical hydromechanical model. Our seismic model shows a strong high-velocity anomaly at 130-167 km coinciding with the BN seismicity. The relocated seismicity can be separated in two clusters. We propose that the upper BN cluster at ~ 130 km depth is attributed to dehydration embrittlement, whereas the lower BN cluster at ~ 150 km depth coinciding with the high-velocity body is caused by lithospheric delamination, creating a "drip" that falls down over the subducting oceanic plate, enhancing fluid release from the slab, potentially increasing seismicity.

Low-degree mantle melting controls the deep seismicity and explosive volcanism of the Gakkel Ridge.

The world's strongest known spreading-related seismicity swarm occurred in 1999 in a segment of the Gakkel Ridge located at 85°E as a consequence of an effusive-explosive submarine volcanic eruption. The data of a seismic network deployed on ice floes were used to locate hundreds of local earthquakes down to ∼25 km depth and to build a seismic tomography model under the volcanic area. Here we show the seismicity and the distribution of seismic velocities together with the 3D magmatic-thermomechanical numerical model, which demonstrate how a magma reservoir under the Gakkel Ridge may form, rise and trigger volcanic eruptions in the rift valley. The ultraslow spreading rates with low mantle potential temperatures appear to be a critical factor in the production of volatile-rich, low-degree mantle melts that are focused toward the magma reservoirs within narrow magmatic sections. The degassing of these melts is the main cause of the explosive submarine eruptions.

Inversion in the permeability evolution of deforming Westerly granite near the brittle-ductile transition.

Fluid flow through crustal rocks is controlled by permeability. Underground fluid flow is crucial in many geotechnical endeavors, such as CO2 sequestration, geothermal energy, and oil and gas recovery. Pervasive fluid flow and pore fluid pressure control the strength of a rock and affect seismicity in tectonic and geotechnical settings. Despite its relevance, the evolution of permeability with changing temperature and during deformation remains elusive. In this study, the permeability of Westerly granite at an effective pressure of 100 MPa was measured under conditions near its brittle-ductile transition, between 650 °C and 850 °C, with a strain rate on the order of 2·10-6 s-1. To capture the evolution of permeability with increasing axial strain, the samples were continuously deformed in a Paterson gas-medium triaxial apparatus. The microstructures of the rock were studied after testing. The experiments reveal an inversion in the permeability evolution: an initial decrease in permeability due to compaction and then an increase in permeability shortly before and immediately after failure. The increase in permeability after failure, also present at high temperatures, is attributed to the creation of interconnected fluid pathways along the induced fractures. This systematic increase demonstrates the subordinate role that temperature dilatancy plays in permeability control compared to stress and its related deformation. These new experimental results thus demonstrate that permeability enhancement under brittle-ductile conditions unveils the potential for EGS exploitation in high-temperature rocks.

Transition from continental rifting to oceanic spreading in the northern Red Sea area.

Lithosphere extension, which plays an essential role in plate tectonics, occurs both in continents (as rift systems) and oceans (spreading along mid-oceanic ridges). The northern Red Sea area is a unique natural geodynamic laboratory, where the ongoing transition from continental rifting to oceanic spreading can be observed. Here, we analyze travel time data from a merged catalogue provided by the Egyptian and Saudi Arabian seismic networks to build a three-dimensional model of seismic velocities in the crust and uppermost mantle beneath the northern Red Sea and surroundings. The derived structures clearly reveal a high-velocity anomaly coinciding with the Red Sea basin and a narrow low-velocity anomaly centered along the rift axis. We interpret these structures as a transition of lithospheric extension from continental rifting to oceanic spreading. The transitional lithosphere is manifested by a dominantly positive seismic anomaly indicating the presence of a 50-70-km-thick and 200-300-km-wide cold lithosphere. Along the forming oceanic ridge axis, an elongated low-velocity anomaly marks a narrow localized nascent spreading zone that disrupts the transitional lithosphere. Along the eastern margins of the Red Sea, several low-velocity anomalies may represent crustal zone of massive Cenozoic basaltic magmatism.

Ongoing formation of felsic lower crustal channel by relamination in Zagros collision zone revealed from regional tomography.

Complex interaction of rheologically contrasting layers within the lithosphere during the collision of continental plates leads to active faulting, which represents a serious hazard to the population and infrastructure. One of the collision scenarios presumes the existence of a middle-lower crustal channel composed of subducted silicic upper crustal rocks, which is thought to exist in the Tibetan-Himalayan system. Based on the results of seismic tomography, we argue that a similar mechanism of crustal channeling takes place beneath the Zagros mountain system in southwestern Iran. The 3D seismic velocity model reveals an inverted crustal architecture of the collision zone, in which the low-velocity felsic (granitic and sedimentary) upper crustal rocks of the Arabian plate form a seismically inactive lower crustal channel below the higher-velocity mafic (basaltic) middle-upper crustal layer of the Iranian crust. Based on existing numerical models, we suggest that the formation of the felsic channel is likely governed by separation (delamination) of the weak felsic upper crust of the subducting Arabian lithosphere and its ductile underplating under rheologically stronger upper-middle crust of the Iranian plate.

Oceanic crust recycling controlled by weakening at slab edges.

Retreating subduction zones such as the Lesser Antilles, Gibraltar and Scotia have been migrating towards the Atlantic Ocean by cutting their way through the oceanic crust. This spontaneously retreating subduction is enabled by the development of faults at the edges of the slab, but the physical mechanisms controlling fault propagation and direction remain unknown. Here, using 3D numerical subduction models we show that oceanic lithosphere recycling is mainly controlled by the intensity of strain-induced weakening of fractures forming at the edges of the slab. Intense strain-induced weakening causes predominantly brittle fault propagation and slab narrowing until detachment. Without weakening, preponderantly ductile slab edge propagation occurs, which causes slab widening. This rheological control is not affected by the proximity of non-weakened passive continental margins. Natural examples suggest that slab edges follow convergent paths that could be controlled by fractures weakening due to deep water penetration into the oceanic lithosphere.

Transient stripping of subducting slabs controls periodic forearc uplift.

Topography in forearc regions reflects tectonic processes along the subduction interface, from seismic cycle-related transients to long-term competition between accretion and erosion. Yet, no consensus exists about the topography drivers, especially as the contribution of deep accretion remains poorly constrained. Here, we use thermo-mechanical simulations to show that transient slab-top stripping events at the base of the forearc crust control uplift-then-subsidence sequences. This 100s-m-high topographic signal with a Myr-long periodicity, mostly inaccessible to geodetic and geomorphological records, reflects the nature and influx rate of material involved in the accretion process. The protracted succession of stripping events eventually results in the pulsing rise of a large, positive coastal topography. Trench-parallel alternation of forearc highs and depressions along active margins worldwide may reflect temporal snapshots of different stages of these surface oscillations, implying that the 3D shape of topography enables tracking deep accretion and associated plate-interface frictional properties in space and time.

Lateral propagation-induced subduction initiation at passive continental margins controlled by preexisting lithospheric weakness.

Understanding the conditions for forming new subduction zones at passive continental margins is important for understanding plate tectonics and the Wilson cycle. Previous models of subduction initiation (SI) at passive margins generally ignore effects due to the lateral transition from oceanic to continental lithosphere. Here, we use three-dimensional numerical models to study the possibility of propagating convergent plate margins from preexisting intraoceanic subduction zones along passive margins [subduction propagation (SP)]. Three possible regimes are achieved: (i) subducting slab tearing along a STEP fault, (ii) lateral propagation-induced SI at passive margin, and (iii) aborted SI with slab break-off. Passive margin SP requires a significant preexisting lithospheric weakness and a strong slab pull from neighboring subduction zones. The Atlantic passive margin to the north of Lesser Antilles could experience SP if it has a notable lithospheric weakness. In contrast, the Scotia subduction zone in the Southern Atlantic will most likely not propagate laterally.

Stress-driven fluid flow controls long-term megathrust strength and deep accretionary dynamics.

The heterogeneity of frictional strength along the megathrust earthquake zone critically controls plate coupling and long-term subduction dynamics. However, the persistence and distribution of high-friction segments through space and time remain poorly constrained. Here, we show that accretion processes, such as tectonic underplating (i.e., basal accretion of material below the fore-arc region), can be used as a proxy to characterize the long-term frictional zonation of the subduction interface. We carry out numerical thermo-mechanical experiments, which predict a first-order control of tectonic-stress variations on fluid transport in deep fore-arc regions. Accordingly, positive feedback between fluid distribution and effective stress favours the stability of the interface frictional properties at Myr-scale which, in turn, controls the deep accretionary dynamics. We propose that the recognition of thick duplex structures resulting from successive underplating events over tens of Myr, allows for tracking subduction segments exhibiting an increasing frictional behaviour. Our numerical results help ascertain the long-term hydro-mechanical properties and distribution of coupling/decoupling segments of megathrust earthquake zones worldwide where active tectonic underplating is recognized.

Bimodal seismicity in the Himalaya controlled by fault friction and geometry.

There is increasing evidence that the Himalayan seismicity can be bimodal: blind earthquakes (up to Mw ~ 7.8) tend to cluster in the downdip part of the seismogenic zone, whereas infrequent great earthquakes (Mw 8+) propagate up to the Himalayan frontal thrust. To explore the causes of this bimodal seismicity, we developed a two-dimensional, seismic cycle model of the Nepal Himalaya. Our visco-elasto-plastic simulations reproduce important features of the earthquake cycle, including interseismic strain and a bimodal seismicity pattern. Bimodal seismicity emerges as a result of relatively higher friction and a non-planar geometry of the Main Himalayan Thrust fault. This introduces a region of large strength excess that can only be activated once enough stress is transferred upwards by blind earthquakes. This supports the view that most segments of the Himalaya might produce complete ruptures significantly larger than the 2015 Mw 7.8 Gorkha earthquake, which should be accounted for in future seismic hazard assessments.

Afar triple junction triggered by plume-assisted bi-directional continental break-up.

Divergent ridge-ridge-ridge (R-R-R) triple junctions are one of the most remarkable, yet largely enigmatic, features of plate tectonics. The juncture of the Arabian, Nubian, and Somalian plates is a type-example of the early development stage of a triple junction where three active rifts meet at a 'triple point' in Central Afar. This structure may result from the impingement of the Afar plume into a non-uniformly stressed continental lithosphere, but this process has never been reproduced by self-consistent plume-lithosphere interaction experiments. Here we use 3D thermo-mechanical numerical models to examine the initiation of plume-induced rift systems under variable far-field stress conditions. Whereas simple linear rift structures are preferred under uni-directional extension, we find that more complex patterns form in response to bi-directional extension, combining one or several R-R-R triple junctions. These triple junctions optimize the geometry of continental break-up by minimizing the amount of dissipative mechanical work required to accommodate multi-directional extension. Our models suggest that Afar-like triple junctions are an end-member mode of plume-induced bi-directional rifting that combines asymmetrical northward pull and symmetrical EW extension at similar rates.

Multi-terrane structure controls the contrasting lithospheric evolution beneath the western and central-eastern Tibetan plateau.

The Tibetan plateau is manifested by contrasting along-strike lithospheric structures, but its formation mechanism and the relationship with the heterogeneous multi-terrane configuration is a challenging problem. Here we conduct systematic numerical modeling to explore the roles of width, density, and rheological properties of the multiple terranes in the lithospheric evolution of the Tibetan plateau, which reveals two distinct collision modes. In Mode-I, the lithospheric mantles of both the strong and weak terranes in the Tibetan plate are completely detached, followed by the underthrusting of Indian lithosphere beneath the whole plateau. Alternatively, Mode-II is characterized by full detachment of the weak terranes, but (partial) residue of the strong terranes during collision. These two contrasting modes, broadly consistent with the lithospheric structures of western and central-eastern Tibetan plateau, respectively, are strongly dependent on the along-strike variation of the width of the strong Lhasa-Qiangtang terranes.

Long-distance impact of Iceland plume on Norway's rifted margin.

Results of a 3D modeling study inspired by recent seismic tomography of the Northern Atlantic mantle suggest that a complex pattern of hot mantle distribution with long horizontal flows originating from the Iceland mantle plume has been the norm in the geological past. In the Northern Atlantic the Iceland plume has a strong long-distance impact on intraplate deformation affecting both onshore and offshore parts of Norway's rifted margin. As a result, this margin is characterized by large magnitude differential topography sustained over at least several tens of Myr. Here we use high-resolution 3D thermo-mechanical modeling to demonstrate that the long-distance plume impact can be explained by its fast lateral propagation controlled by pre-existing lithosphere structures. Numerical models show that these structures strongly affect the style of horizontal flow of plume head material. This results in long-distance propagation of hot material emplaced at the lithosphere-asthenosphere boundary causing long-wavelength anomalies in onshore topography of Norway's rifted margin. Short-wavelength offshore topographic domes are likely caused by joint occurrence of plume-related thermal perturbations and gravitational forces related to plate thickening (ridge push). Our 3D modeling brings together plume impingement, spreading ridge dynamics, and the formation of anomalous intraplate structures offshore Norway in one scenario.