Rheology of the upper mantle: a synthesis

Growth of continental crust and lithosphere subduction in the Hadean revealed by geochemistry and geodynamics

The rates of continental crust growth and recycling on early Earth remain unclear due to the lack of information resulting from the extensive alteration of ancient rocks. Melt inclusions trapped and shielded from alteration in Archean high-Mg olivine crystals offer a solution to this problem. We report an unprecedented unradiogenic Sr mantle source component (87Sr/86Sr = 0.69932 ± 0.00024, 95% confidence interval) of melts included in olivine from 3.27 Ga komatiitic lava flows in the Barberton Greenstone Belt, South Africa. This component indicates a model age of 4.31 ± 0.19 Ga and significant chemical fractionation (Nb/U = 36.9 ± 1.5, Ce/Pb=16.7 ± 1.1), suggesting up to 80% ± 16% of the present-day continental crust's mass was extracted by the late Hadean from the whole mantle. Geodynamic models support this finding, explaining geochemical data by producing 40% to 70% of the present-day continental crust mass during the Hadean in a variable tectonic regime with tens of millions of years-long periods of massive impulsive subduction induced by mantle plumes.

Terrane accretion explains thin and hot ocean-continent back-arcs

The origin of hot ocean-continent back-arc regions with very thin mantle lithosphere and very high surface heat flow in both extensional and contractional ocean-continent subduction systems is highly enigmatic and unresolved. These first-order characteristics have often been explained with either convective mantle lithosphere removal or by back-arc extension. However, it is unclear what may cause the proposed convective thinning and/or delamination of eclogitic lower crust over very wide regions, whereas back-arc extension is either not observed or insufficient to explain the observed very thin mantle lithosphere. Notably, many of these ocean-continent systems have a long history of terrane accretion. Here, we show, using thermomechanical model experiments, that terrane accretion provides a consistent explanation for the observed key characteristics and naturally leads to rheologically weak back-arcs with continental crust directly on top of the hot sublithospheric mantle.

Tidal Deformation and Dissipation Processes in Icy Worlds

Tidal interactions play a key role in the dynamics and evolution of icy worlds. The intense tectonic activity of Europa and the eruption activity on Enceladus are clear examples of the manifestation of tidal deformation and associated dissipation. While tidal heating has long been recognized as a major driver in the activity of these icy worlds, the mechanism controlling how tidal forces deform the different internal layers and produce heat by tidal friction still remains poorly constrained. As tidal forcing varies with orbital characteristics (distance to the central planet, eccentricity, obliquity), the contribution of tidal heating to the internal heat budget can strongly change over geological timescales. In some circumstances, the tidally-produced heat can result in internal melting and surface activity taking various forms. Even in the absence of significant heat production, tidal deformation can be used to probe the interior structure, the tidal response of icy moons being strongly sensitive to their hydrosphere structure. In the present paper, we review the methods to compute tidal deformation and dissipation in the different layers composing icy worlds. After summarizing the main principle of tidal deformation and the different rheological models used to model visco-elastic tidal response, we describe the dissipation processes expected in rock-dominated cores, subsurface oceans and icy shells and highlight the potential effects of tidal heating in terms of thermal evolution and activity. We finally anticipate how data collected by future missions to Jupiter's and Saturn's moons could be used to constrain their tidal response and the consequences for past and present activities.

Effect of bridgmanite-ferropericlase grain size evolution on Earth's average mantle viscosity: implications for mantle convection in early and present-day Earth

Recent experimental investigations of grain size evolution in bridgmanite-ferropericlase assemblages have suggested very slow growth for these bimodal phases. Despite numerous speculations on grain size-dependent viscosity, a comprehensive test with realistic grain size evolution parameters compatible with the lower mantle has been lacking. In this study, we develop self-consistent 2-D spherical half-annulus geodynamic models of Earth's evolution using the finite volume code StagYY to assess the role of grain size on lower mantle viscosity. We explore several models with and without grain size evolution to compare their effects on mantle viscosity. In models with grain size evolution, we consider three scenarios: (1) uniform grain growth throughout the entire mantle with a composite rheology, (2) different grain growth in the upper and lower mantle with a composite rheology, and (3) different grain growth in the upper and lower mantle with purely diffusion creep rheology. In the case of different grain size evolution, the upper mantle's grain size evolution law is controlled by forsterite-enstatite grain growth, while the lower mantle's grain size evolution law is controlled by bridgmanite-ferropericlase grain growth. Our results suggest that mantle viscosity is primarily controlled by temperature, whereas grain size has a minor effect compared to the effect of temperature. We attribute two primary reasons for this: First, the bridgmanite-ferropericlase growth is very slow in the lower mantle and the grain size variation is too small to significantly alter the mantle viscosity. Secondly, if grains grow too fast, thus the mantle deforms in the dislocation creep regime, making viscosity grain size-independent. To establish the robustness of this finding we vary several other model parameters, such as surface yield strength, phase transition grain size reset, different transitional stresses for creep mechanisms, pressure dependence on grain growth, and different grain damage parameters. For all our models, we consistently find that grain size has a very limited effect on controlling lower mantle viscosity in the present-day Earth. However, large grain size may have affected the lower mantle viscosity in the early Earth as larger grains of single phase bridgmanite could increase the viscosity of the early mantle delaying the onset of global convection.

Coordination between deformation, precipitation, and erosion during orogenic growth

Crustal thickening associated with orogenic growth elevates topography, causing orographic enhancement of precipitation, which in turn facilitates local erosion and possibly intensifies localization of deformation. How these three processes-deformation, precipitation, and erosion-coordinate during orogenic growth remains unknown. Here, we present a numerical model where tectonics, surface processes, and orographic precipitation are tightly coupled, and explore the impact on low, intermediate, and high erodibility orogens. We show that, for intermediate erodibility models, rock uplift rates and precipitation rates correlate well with erosion rates during the formation of orogenic plateaus with high correlation coefficients of ~0.9 between rock uplift and erosion rates, and ~0.8 between precipitation and erosion rates. We demonstrate a cyclicity of correlation evolution among uplift, precipitation, and erosion rates through the development of new faults propagating outward. These results shed insights into the relative tectonic or climatic control on erosion in active orogens (e.g., Himalayas, Central Andes, and Southern Alps of New Zealand), and provide a plausible explanation for several conflicting data and interpretations in the Himalayas, which depend on the stage of maturity of the newest fault and the relative locations to old faults.

Rapid shear zone weakening during subduction initiation

Subduction zones play a pivotal role in the mechanics of plate tectonics by providing the driving force through slab pull and weak megathrusts that facilitate the relative motion between tectonic plates. The initiation of subduction zones is intricately linked to the accumulation of slab pull and development of weakness at plate boundaries and, by consequence, the largest changes in the energetics of mantle convection. However, the transient nature of subduction initiation accompanied by intense subsequent tectonic activity, leaves critical evidence poorly preserved and making subduction initiation difficult to constrain. We overcome these limitations through a comprehensive analysis focused on Puysegur, a well-constrained extant example of subduction initiation offshore South Island, New Zealand. Through time-dependent, three-dimensional thermo-mechanical computations and quantitative comparison to new geophysical and geological observations, including topography, stratigraphy, and seismicity, we demonstrate that subduction initiation develops with a fast strain weakening described with a small characteristic displacement ([Formula: see text] 4 to 8 km). Potential physical mechanisms contributing to the strain weakening are explored and we find that the observed fast weakening may arise through a combination of grain-size reduction within the lower lithosphere and fluid pressurization at shallower depths. With the shared commonality in the underlying physics of tectonic processes, the rapid strain weakening constrained at Puysegur offers insights into the formation of the first subduction during early Earth and the onset of plate tectonics.

Coevolution of craton margins and interiors during continental break-up

Many cratonic continental fragments dispersed during the rifting and break-up of Gondwana are bound by steep topographic landforms known as 'great escarpments'1-4, which rim elevated plateaus in the craton interior5,6. In terms of formation, escarpments and plateaus are traditionally considered distinct owing to their spatial separation, occasionally spanning more than a thousand kilometres. Here we integrate geological observations, statistical analysis, geodynamic simulations and landscape-evolution models to develop a physical model that mechanistically links both phenomena to continental rifting. Escarpments primarily initiate at rift-border faults and slowly retreat at about 1 km Myr-1 through headward erosion. Simultaneously, rifting generates convective instabilities in the mantle7-10 that migrate cratonward at a faster rate of about 15-20 km Myr-1 along the lithospheric root, progressively removing cratonic keels11, driving isostatic uplift of craton interiors and forming a stable, elevated plateau. This process forces a synchronized wave of denudation, documented in thermochronology studies, which persists for tens of millions of years and migrates across the craton at a comparable or slower pace. We interpret the observed sequence of rifting, escarpment formation and exhumation of craton interiors as an evolving record of geodynamic mantle processes tied to continental break-up, upending the prevailing notion of cratons as geologically stable terrains.

Geologic Constraints on the Formation and Evolution of Saturn's Mid-Sized Moons

Saturn's mid-sized icy moons have complex relationships with Saturn's interior, the rings, and with each other, which can be expressed in their shapes, interiors, and geology. Observations of their physical states can, thus, provide important constraints on the ages and formation mechanism(s) of the moons, which in turn informs our understanding of the formation and evolution of Saturn and its rings. Here, we describe the cratering records of the mid-sized moons and the value and limitations of their use for constraining the histories of the moons. We also discuss observational constraints on the interior structures of the moons and geologically-derived inferences on their thermal budgets through time. Overall, the geologic records of the moons (with the exception of Mimas) include evidence of epochs of high heat flows, short- and long-lived subsurface oceans, extensional tectonics, and considerable cratering. Curiously, Mimas presents no clear evidence of an ocean within its surface geology, but its rotation and orbit indicate a present-day ocean. While the moons need not be primordial to produce the observed levels of interior evolution and geologic activity, there is likely a minimum age associated with their development that has yet to be determined. Uncertainties in the populations impacting the moons makes it challenging to further constrain their formation timeframes using craters, whereas the characteristics of their cores and other geologic inferences of their thermal evolutions may help narrow down their potential histories. Disruptive collisions may have also played an important role in the formation and evolution of Saturn's mid-sized moons, and even the rings of Saturn, although more sophisticated modeling is needed to determine the collision conditions that produce rings and moons that fit the observational constraints. Overall, the existence and physical characteristics of Saturn's mid-sized moons provide critical benchmarks for the development of formation theories.

Constraining Earth's nonlinear mantle viscosity using plate-boundary resolving global inversions

Variable viscosity in Earth's mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 [Formula: see text] 0.25 (mean [Formula: see text] SD). A relatively low yield stress of 151 [Formula: see text] 19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes.

Understanding subduction infancy to mature subduction in Southwest Japan via the self-consistent formation of a weak slab interface

The weak slab interface controls long-term subduction dynamics. A weak hydrous layer at the slab interface promotes mechanical decoupling between the forearc mantle and the subducting slab and converts a hot forearc mantle to a cold mantle. Often referred to as a cold nose, the cold forearc mantle, plays a key role in the transition from subduction infancy to mature subduction. This study was the first to numerically demonstrate the self-consistent formation of a weak hydrous layer with permeability anisotropy based on the Southwest Japan subduction zone case, where transition-related geological features were present. Our models showed that mechanical decoupling by spontaneous downdip growth of the weak hydrous layer created a cold nose by converting a hot forearc mantle to a cold mantle. The emergence of a cold nose explained the migration of the forearc-to-arc volcanic front, expressed as the formation of mid-Miocene forearc high-magnesium andesite and Quaternary arc adakite. Furthermore, the weak hydrous layer providing a pathway for free-water transport toward the mantle wedge tip elucidates slab/mantle-derived geochemical components in deep groundwater as well as large S-wave delay times and non-volcanic seismic tremors in the forearc.

Paleogene India-Eurasia collision constrained by observed plate rotation

The Cenozoic India-Eurasia collision has had profound impacts on shaping the Tibetan plateau, but its early history remains controversial due to uneven availability of constraints. Recent plate reconstructions reveal two prominent counterclockwise rotation (azimuthal change) rate peaks of the Indian plate at 52-44 and 33-20 Ma, respectively, which could bear key information about this collision history. Using fully dynamic three-dimensional numerical modeling, we show that the first rotation rate peak reflected the initial diachronous collision from the western-central to eastern Indian front, and the second peak reflected the full collision leading to strong coupling between India and Eurasia. Further comparison with observation suggests that the initial and complete India-Eurasia collision likely occurred at 55 ± 5 and 40 ± 5 Ma, respectively, an inference consistent with key geological observations. We suggest that this collision history is instructive for studying the tectonic history of the Tibetan plateau and its surrounding areas.

Geophysical evidence for an enriched molten silicate layer above Mars's core

The detection of deep reflected S waves on Mars inferred a core size of 1,830 ± 40 km (ref. 1), requiring light-element contents that are incompatible with experimental petrological constraints. This estimate assumes a compositionally homogeneous Martian mantle, at odds with recent measurements of anomalously slow propagating P waves diffracted along the core-mantle boundary2. An alternative hypothesis is that Mars's mantle is heterogeneous as a consequence of an early magma ocean that solidified to form a basal layer enriched in iron and heat-producing elements. Such enrichment results in the formation of a molten silicate layer above the core, overlain by a partially molten layer3. Here we show that this structure is compatible with all geophysical data, notably (1) deep reflected and diffracted mantle seismic phases, (2) weak shear attenuation at seismic frequency and (3) Mars's dissipative nature at Phobos tides. The core size in this scenario is 1,650 ± 20 km, implying a density of 6.5 g cm-3, 5-8% larger than previous seismic estimates, and can be explained by fewer, and less abundant, alloying light elements than previously required, in amounts compatible with experimental and cosmochemical constraints. Finally, the layered mantle structure requires external sources to generate the magnetic signatures recorded in Mars's crust.

Influence of the asthenosphere on earth dynamics and evolution

The existence of a thin, weak asthenospheric layer beneath Earth's lithospheric plates is consistent with existing geological and geophysical constraints, including Pleistocene glacio-isostatic adjustment, modeling of gravity anomalies, studies of seismic anisotropy, and post-seismic rebound. Mantle convection models suggest that a pronounced weak zone beneath the upper thermal boundary layer (lithosphere) may be essential to the plate tectonic style of convection found on Earth. The asthenosphere is likely related to partial melting and the presence of water in the sub-lithospheric mantle, further implying that the long-term evolution of the Earth may be controlled by thermal regulation and volatile recycling that maintain a geotherm that approaches the wet mantle solidus at asthenospheric depths.

The lunar solid inner core and the mantle overturn

Seismological models from Apollo missions provided the first records of the Moon inner structure with a decrease in seismic wave velocities at the core-mantle boundary^1-3. The resolution of these records prevents a strict detection of a putative lunar solid inner core and the impact of the lunar mantle overturn in the lowest part of the Moon is still discussed^4-7. Here we combine geophysical and geodesic constraints from Monte Carlo exploration and thermodynamical simulations for different Moon internal structures to show that only models with a low viscosity zone enriched in ilmenite and an inner core present densities deduced from thermodynamic constraints compatible with densities deduced from tidal deformations. We thus obtain strong indications in favour of the lunar mantle overturn scenario and, in this context, demonstrate the existence of the lunar inner core with a radius of 258 ± 40 km and density 7,822 ± 1,615 kg m^-3. Our results question the evolution of the Moon magnetic field thanks to its demonstration of the existence of the inner core and support a global mantle overturn scenario that brings substantial insights on the timeline of the lunar bombardment in the first billion years of the Solar System⁸.

Ghost-arc geochemical anomaly at a spreading ridge caused by supersized flat subduction

The Southern Atlantic-Southwest Indian ridges (SASWIR) host mid-ocean ridge basalts with a residual subduction-related geochemical fingerprint (i.e., a ghost-arc signature) of unclear origin. Here, we show through an analysis of plate kinematic reconstructions and seismic tomography models that the SASWIR subduction-modified mantle source formed in the Jurassic close to the Georgia Islands slab (GI) and remained near-stationary in the mantle reference frame. In this analysis, the GI lies far inboard the Jurassic Patagonian-Antarctic Peninsula active margin. This was formerly attributed to a large-scale flat subduction event in the Late Triassic-Early Jurassic. We propose that during this flat slab stage, the subduction-modified mantle areas beneath the Mesozoic active margin and surrounding sutures zones may have been bulldozed inland by >2280 km. After the demise of the flat slab, this mantle anomaly remained near-stationary and was sampled by the Karoo mantle plume 183 Million years (Myr) ago and again since 55 Myr ago by the SASWIR. We refer to this process as asthenospheric anomaly telescoping. This study provides a hitherto unrecognized geodynamic effect of flat subduction, the viability of which we support through numerical modeling.

A Positive Feedback Between Crustal Thickness and Melt Extraction for the Origin of the Martian Dichotomy

A North/South difference in crustal thickness is likely at the origin of the Martian dichotomy in topography. Recent crustal thickness maps were obtained by inversion of topography and gravity data seismically anchored at the InSight station. On average, the Martian crust is 51-71 km thick with a southern crust thicker by 18-28 km than the northern one. The origin of this crustal dichotomy is still debated although the hypothesis of a large impact is at present very popular. Here, we propose a new mechanism for the formation of this dichotomy that involves a positive feedback between crustal growth and mantle melting. As the crust is enriched in heat-producing elements, the lid of a one-plate planet is hotter and thinner where the crust is thicker, inducing a larger amount of partial melt below the lid and hence a larger rate of melt extraction and crustal growth. We first demonstrate analytically that larger wavelength perturbations, that is, hemispherical perturbations, grow faster because smaller wavelengths are more attenuated by thermal diffusion. We then use a parameterized thermal evolution model with a well-mixed mantle topped by two different lids characterized by their thermal structures and thicknesses to study the growth of the crust in the two hemispheres. Our results demonstrate that this positive feedback can generate a significant crustal dichotomy.

The role of gravitational body forces in the development of metamorphic core complexes

Within extreme continental extension areas, ductile middle crust is exhumed at the surface as metamorphic core complexes. Sophisticated quantitative models of extreme extension predicted upward transport of ductile middle-lower crust through time. Here we develop a general model for metamorphic core complexes formation and demonstrate that they result from the collapse of a mountain belt supported by a thickened crustal root. We show that gravitational body forces generated by topography and crustal root cause an upward flow pattern of the ductile lower-middle crust, facilitated by a detachment surface evolving into low-angle normal fault. This detachment surface acquires large amounts of finite strain, consistent with thick mylonite zones found in metamorphic core complexes. Isostatic rebound exposes the detachment in a domed upwarp, while the final Moho discontinuity across the extended region relaxes to a flat geometry. This work suggests that belts of metamorphic core complexes are a fossil signature of collapsed highlands.

Glacial isostatic adjustment: physical models and observational constraints

By far the most prescient insights into the interior structure of the planet have been provided on the basis of elastic wave seismology. Analysis of the travel times of shear or compression wave phases excited by individual earthquakes, or through analysis of the elastic gravitational free oscillations that individual earthquakes of sufficiently large magnitude may excite, has been the central focus of Earth physics research for more than a century. Unfortunately, data provide no information that is directly relevant to understanding the solid state 'flow' of the polycrystalline outer 'mantle' shell of the planet that is involved in the thermally driven convective circulation that is responsible for powering the 'drift' of the continents and which controls the rate of planetary cooling on long timescales. For this reason, there has been an increasing focus on the understanding of physical phenomenology that is unambiguously associated with mantle flow processes that are distinct from those directly associated with the convective circulation itself. This paper reviews the past many decades of work that has been invested in understanding the most important of such processes, namely that which has come to be referred to as 'glacial isostatic adjustment' (GIA). This process concerns the response of the planet to the loading and unloading of the high latitude continents by the massive accumulations of glacial ice that have occurred with almost metronomic regularity over the most recent million years of Earth history. Forced by the impact of gravitationaln-body effects on the geometry of Earth's orbit around the Sun through the impact upon the terrestrial regime of received solar insolation, these surface mass loads on the continents have left indelible records of their occurrence in the 'Earth system' consisting of the oceans, continents, and the great polar ice sheets on Greenland and Antarctica themselves. Although this ice-age phenomenology has been clearly recognized since early in the last century, it was for over 50 years considered to be no more than an interesting curiosity, the understanding of which remained on the periphery of the theoretical physics of the Earth. This was the case in part because no globally applicable theory was available that could be applied to rigorously interpret the observations. Equally important to understanding the scientific lethargy that held back the understanding of this phenomenon involving mantle flow processes was the lack of appreciation of the wide range of observations that were in fact related to GIA physics. This paper is devoted to a review of the global theories of the GIA process that have since been developed as a means of interpreting the extensive variety of observations that are now recognized as being involved in the response of the planet to the loading and unloading of its surface by glacial ice. The paper will also provide examples of the further analyses of Earth physics and climate related processes that applications of the modern theoretical structures have enabled.

On the Relationship Between Oceanic Plate Speed, Tectonic Stress, and Seismic Anisotropy

Seismic radial anisotropy (the squared ratio between the speeds of horizontally and vertically polarized shear waves,ξ = V S H 2 V S V 2 ) is a powerful tool to probe the direction of mantle flow and accumulated strain. While previous studies have confirmed the dependence of azimuthal anisotropy on plate speed, the first order control on radial anisotropy is unclear. In this study, we develop 2D ridge flow models combined with mantle fabric calculations to report that faster plates generate higher tectonics stresses and strain rates which lower the dislocation creep viscosity and lead to deeper anisotropy than beneath slower plates. We apply the SGLOBE-rani tomographic filter, resulting in a flat depth-age trend and stronger anisotropy beneath faster plates, which correlates well with 3D global anisotropic mantle models. Our predictions and observations suggest that as plate speed increases from 2 to 8 cm/yr, radial anisotropy increases by ∼0.01-0.025 in the upper 100-200 km of the mantle between 10 and 60 Ma.

Coupled influence of tectonics, climate, and surface processes on landscape evolution in southwestern North America

The Cenozoic landscape evolution in southwestern North America is ascribed to crustal isostasy, dynamic topography, or lithosphere tectonics, but their relative contributions remain controversial. Here we reconstruct landscape history since the late Eocene by investigating the interplay between mantle convection, lithosphere dynamics, climate, and surface processes using fully coupled four-dimensional numerical models. Our quantified depth-dependent strain rate and stress history within the lithosphere, under the influence of gravitational collapse and sub-lithospheric mantle flow, show that high gravitational potential energy of a mountain chain relative to a lower Colorado Plateau can explain extension directions and stress magnitudes in the belt of metamorphic core complexes during topographic collapse. Profound lithospheric weakening through heating and partial melting, following slab rollback, promoted this extensional collapse. Landscape evolution guided northeast drainage onto the Colorado Plateau during the late Eocene-late Oligocene, south-southwest drainage reversal during the late Oligocene-middle Miocene, and southwest drainage following the late Miocene.

Evidence for active upper mantle flow in the Atlantic and Indo-Australian realms since the Upper Jurassic from hiatus maps and spreading rate changes

Histories of large-scale horizontal and vertical lithosphere motion hold important information on mantle convection. Here, we compare continent-scale hiatus maps as a proxy for mantle flow induced dynamic topography and plate motion variations in the Atlantic and Indo-Australian realms since the Upper Jurassic, finding they frequently correlate, except when plate boundary forces may play a significant role. This correlation agrees with descriptions of asthenosphere flow beneath tectonic plates in terms of Poiseuille/Couette flow, as it explicitly relates plate motion changes, induced by evolving basal shear forces, to non-isostatic vertical motion of the lithosphere. Our analysis reveals a timescale, on the order of a geological series, between the occurrence of continent-scale hiatus and plate motion changes. This is consistent with the presence of a weak upper mantle. It also shows a spatial scale for interregional hiatus, on the order of 2000-3000 km in diameter, which can be linked by fluid dynamic analysis to active upper mantle flow regions. Our results suggest future studies should pursue large-scale horizontal and vertical lithosphere motion in combination, to track the expressions of past mantle flow. Such studies would provide powerful constraints for adjoint-based geodynamic inverse models of past mantle convection.

Archean cratonic mantle recycled at a mid-ocean ridge

Basalts and mantle peridotites of mid-ocean ridges are thought to sample Earth's upper mantle. Osmium isotopes of abyssal peridotites uniquely preserve melt extraction events throughout Earth history, but existing records only indicate ages up to ~2 billion years (Ga) ago. Thus, the memory of the suspected large volumes of mantle lithosphere that existed in Archean time (>2.5 Ga) has apparently been lost somehow. We report abyssal peridotites with melt-depletion ages up to 2.8 Ga, documented by extremely unradiogenic ¹⁸⁷Os/¹⁸⁸Os ratios (to as low as 0.1095) and refractory major elements that compositionally resemble the deep keels of Archean cratons. These oceanic rocks were thus derived from the once-extensive Archean continental keels that have been dislodged and recycled back into the mantle, the feasibility of which we confirm with numerical modeling. This unexpected connection between young oceanic and ancient continental lithosphere indicates an underappreciated degree of compositional recycling over time.

Topography of mountain belts controlled by rheology and surface processes

It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics^1-4. However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process^5,6 and mantle-scale tectonic model⁷. End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state⁸, strength controlled (Bm ≈ 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.

Mantle exhumation at magma-poor rifted margins controlled by frictional shear zones

The transition zone from continental crust to the mature mid-ocean ridge spreading center of the Iberia-Newfoundland magma-poor rifted margins is mostly composed of exhumed mantle characterized by highs and domes with varying elevation, spacing and shape. The mechanism controlling strain localization and fault migration explaining the geometry of these peridotite ridges is poorly understood. Here we show using forward geodynamic models that multiple out-of-sequence detachments with recurring dip reversal form during magma-poor rifting and mantle exhumation as a consequence of the strength competition between weak frictional-plastic shear zones and the thermally weakened necking domain beneath the exhuming footwall explaining geometry of these peridotite ridges. Model behaviour also shows that fault types and detachment styles vary with spreading rate and fault strength and confirm that these results can be compared to other magma poor passive margins such as along Antarctica-Australia and to ultra-slow mid-ocean spreading systems as the South-West Indian Ridge.

Formation of the Lunar Primary Crust From a Long-Lived Slushy Magma Ocean

Classical fractional crystallization scenarios of early lunar evolution suggest crustal formation by the flotation of light anorthite minerals from a liquid magma ocean. However, this model is challenged by the> 200 Myr age range of primitive ferroan anorthosites, their concordance with Mg-suite magmatism and by the compositional diversity observed in lunar anorthosites. Here, we propose a new model of slushy magma ocean crystallization in which crystals remain suspended in the lunar interior and crust formation only begins once a critical crystal content is reached. Thereafter crustal formation occurs by buoyant melt extraction and magmatism. The mixture viscosity strongly depends on temperature and solid fraction driving the development of a surface stagnant lid where enhanced solidification and buoyant ascent of melt lead to an anorthite-enriched crust. This model explains lunar anorthosites heterogeneity and suggests a crustal formation timescale of 100s Ma, reconciling anorthosite ages with an early age of the Moon.

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 faulting^1-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 slabs¹⁰, the development of large-offset normal faults^6,7 near trenches, the occurrence of segmented seismic velocity anomalies¹¹ and distinct interfaces imaged within subducted slabs^12,13, and the appearance of deep, localized intraplate areas of reduced effective viscosity¹⁴ observed at trenches. Furthermore, brittle-viscously damaged slabs show a tendency for detachment at elevated mantle temperatures. Given Earth's planetary cooling history¹⁵, this implies that intermittent subduction with frequent slab break-off episodes¹⁶ may have been characteristic for Earth until more recent times than previously suggested¹⁷.

Was There Land on the Early Earth?

The presence of exposed land on the early Earth is a prerequisite for a certain type of prebiotic chemical evolution in which the oscillating activity of water, driven by short-term, day–night, and seasonal cycles, facilitates the synthesis of proto-biopolymers. Exposed land is, however, not guaranteed to exist on the early Earth, which is likely to have been drastically different from the modern Earth. This mini-review attempts to provide an up-to-date account on the possibility of exposed land on the early Earth by integrating recent geological and geophysical findings. Owing to the competing effects of the growing ocean and continents in the Hadean, a substantial expanse of the Earth’s surface (∼20% or more) could have been covered by exposed continents in the mid-Hadean. In contrast, exposed land may have been limited to isolated ocean islands in the late Hadean and early Archean. The importance of exposed land during the origins of life remains an open question.

The Role of Adria Plate Lithospheric Structures on the Recent Dynamics of the Central Mediterranean Region

The Tertiary tectonic evolution of the Central Mediterranean has been relatively well constrained by abundant geological data. Yet, several uncertainties persist about the mechanisms that led to the present-day surface morphology and deep slab geometry. Here, we combine geodynamic and seismological numerical modeling techniques to reproduce the recent large-scale evolution of the Central Mediterranean and the associated strain-induced upper mantle fabrics and seismic anisotropy. 3D thermo-mechanical subduction models were designed and calibrated according to paleogeographic-tectonic reconstructions and seismological observations available in the literature. It is found that, although the opening of back-arc extensional basins in response to the retreat of the Ionian slab is a common feature in all models, structural heterogeneities within the Adria plate and/or the geometry of its Tyrrhenian passive margin profoundly impact on the segmentation of the subducting slab and the amount of Ionian trench retreat. More, in general, this study highlights the importance of coupling geodynamic and seismological modeling to better constrain the tectonic evolution of complex convergent margins such as the Central Mediterranean.

Differentiating induced versus spontaneous subduction initiation using thermomechanical models and metamorphic soles

Despite the critical role of subduction in plate tectonics, the dynamics of its initiation remains unclear. High-temperature low-pressure metamorphic soles are vestiges of subduction initiation, providing records of the pressure and temperature conditions along the subducting slab surface during subduction initiation that can possibly differentiate the two end-member subduction initiation modes: spontaneous and induced. Here, using numerical models, we show that the slab surface temperature reaches 800-900 °C at ~1 GPa over a wide range of parameter values for spontaneous subduction initiation whereas for induced subduction initiation, such conditions can be reached only if the age of the overriding plate is <5 Ma. These modeling results indicate that spontaneous subduction initiation would be more favorable for creating high-temperature conditions. However, the synthesis of our modeling results and geological observations indicate that the majority of the metamorphic soles likely formed during induced subduction initiation that involved a young overriding plate.

Molecular Insight into the Possible Mechanism of Drag Reduction of Surfactant Aqueous Solution in Pipe Flow

The phenomenon of drag reduction (known as the “Toms effect”) has many industrial and engineering applications, but a definitive molecular-level theory has not yet been constructed. This is due both to the multiscale nature of complex fluids and to the difficulty of directly observing self-assembled structures in nonequilibrium states. On the basis of a large-scale coarse-grained molecular simulation that we conducted, we propose a possible mechanism of turbulence suppression in surfactant aqueous solution. We demonstrate that maintaining sufficiently large micellar structures and a homogeneous radial distribution of surfactant molecules is necessary to obtain the drag-reduction effect. This is the first molecular-simulation evidence that a micellar structure is responsible for drag reduction in pipe flow, and should help in understanding the mechanisms underlying drag reduction by surfactant molecules under nonequilibrium conditions.

Role of warm subduction in the seismological properties of the forearc mantle: An example from southwest Japan

A warm slab thermal structure plays an important role in controlling seismic properties of the slab and mantle wedge. Among warm subduction zones, most notably in southwest Japan, the spatial distribution of large S-wave delay times and deep nonvolcanic tremors in the forearc mantle indicate the presence of a serpentinite layer along the slab interface. However, the conditions under which such a layer is generated remains unclear. Using numerical models, we here show that a serpentinite layer begins to develop by the slab-derived fluids below the deeper end of the slab-mantle decoupling interface and grows toward the corner of the mantle wedge along the interface under warm subduction conditions only, explaining the large S-wave delay times in the forearc mantle. The serpentinite layer then allows continuous free-fluid flow toward the corner of the mantle wedge, presenting possible mechanisms for the deep nonvolcanic tremors in the forearc mantle.

Melt volume at Atlantic volcanic rifted margins controlled by depth-dependent extension and mantle temperature

Breakup volcanism along rifted passive margins is highly variable in time and space. The factors controlling magmatic activity during continental rifting and breakup are not resolved and controversial. Here we use numerical models to investigate melt generation at rifted margins with contrasting rifting styles corresponding to those observed in natural systems. Our results demonstrate a surprising correlation of enhanced magmatism with margin width. This relationship is explained by depth-dependent extension, during which the lithospheric mantle ruptures earlier than the crust, and is confirmed by a semi-analytical prediction of melt volume over margin width. The results presented here show that the effect of increased mantle temperature at wide volcanic margins is likely over-estimated, and demonstrate that the large volumes of magmatism at volcanic rifted margin can be explained by depth-dependent extension and very moderate excess mantle potential temperature in the order of 50-80 °C, significantly smaller than previously suggested.

Amagmatic hydrothermal systems on Mars from radiogenic heat

Long-lived hydrothermal systems are prime targets for astrobiological exploration on Mars. Unlike magmatic or impact settings, radiogenic hydrothermal systems can survive for >100 million years because of the Ga half-lives of key radioactive elements (e.g., U, Th, and K), but remain unknown on Mars. Here, we use geochemistry, gravity, topography data, and numerical models to find potential radiogenic hydrothermal systems on Mars. We show that the Eridania region, which once contained a vast inland sea, possibly exceeding the combined volume of all other Martian surface water, could have readily hosted a radiogenic hydrothermal system. Thus, radiogenic hydrothermalism in Eridania could have sustained clement conditions for life far longer than most other habitable sites on Mars. Water radiolysis by radiogenic heat could have produced H2, a key electron donor for microbial life. Furthermore, hydrothermal circulation may help explain the region's high crustal magnetic field and gravity anomaly.

Thermochemical lithosphere differentiation and the origin of cratonic mantle

Cratons record the early history of continental lithosphere formation, yet how they became the most enduring part of the lithosphere on Earth remains unknown¹. Here we propose a mechanism for the formation of large volumes of melt-depleted cratonic lithospheric mantle (CLM) and its evolution to stable cratons. Numerical models show large decompression melting of a hot, early Earth mantle beneath a stretching lithosphere, where melt extraction leaves large volumes of depleted mantle at depth. The dehydrated, stiffer mantle resists further deformation, forcing strain migration and cooling, thereby assimilating depleted mantle into the lithosphere. The negative feedback between strain localization and stiffening sustains long-term diffused extension and emplacement of large amounts of depleted CLM. The formation of CLM at low pressure and its deeper re-equilibration reproduces the evolution of Archaean lithosphere constrained by depth-temperature conditions^1,2, whereas large degrees of depletion^3,4 and melt volumes⁵ in Archaean cratons are best matched by models with lower lithospheric strength. Under these conditions, which are otherwise viable for plate tectonics^6,7, thermochemical differentiation effectively prevents yielding and formation of margins: rifting and lithosphere subduction are short lived and embedded in the cooling CLM as relict structures, reproducing the recycling and reworking environments that are found in Archaean cratons^8,9. Although they undergo major melting and extensive recycling during an early stage lasting approximately 500 million years, the modelled lithospheres progressively differentiate and stabilize, and then recycling and reworking become episodic. Early major melting and recycling events explain the production and loss of primordial Hadean lithosphere and crust¹⁰, whereas later stabilization and episodic reworking provides a context for the creation of continental cratons in the Archaean era^4,8.

Victoria continental microplate dynamics controlled by the lithospheric strength distribution of the East African Rift

The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria’s rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria’s GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress.

One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.

Rapid subduction initiation and magmatism in the Western Pacific driven by internal vertical forces

Plate tectonics requires the formation of plate boundaries. Particularly important is the enigmatic initiation of subduction: the sliding of one plate below the other, and the primary driver of plate tectonics. A continuous, in situ record of subduction initiation was recovered by the International Ocean Discovery Program Expedition 352, which drilled a segment of the fore-arc of the Izu-Bonin-Mariana subduction system, revealing a distinct magmatic progression with a rapid timescale (approximately 1 million years). Here, using numerical models, we demonstrate that these observations cannot be produced by previously proposed horizontal external forcing. Instead a geodynamic evolution that is dominated by internal, vertical forces produces both the temporal and spatial distribution of magmatic products, and progresses to self-sustained subduction. Such a primarily internally driven initiation event is necessarily whole-plate scale and the rock sequence generated (also found along the Tethyan margin) may be considered as a smoking gun for this type of event.

The magmatic progression produced during the initiation of the Izu-Bonin-Marianas subduction zone took place rapidly over 1 million years, but it has been unclear why. Here, using numerical models, the authors show that subduction initiation was dominated by vertical forces, internal to the system itself, progressing to self-sustained subduction.

An improved setup for radial diffraction experiments at high pressures and high temperatures in a resistive graphite-heated diamond anvil cell

We present an improved setup for the experimental study of deformation of solids at simultaneous high pressures and temperatures by radial x-ray diffraction. This technique employs a graphite resistive heated Mao-Bell type diamond anvil cell for radial x-ray diffraction in combination with a water-cooled vacuum chamber. The new chamber has been developed by the sample environment group at PETRA III and implemented at the Extreme Conditions Beamline P02.2 at PETRA III, DESY (Hamburg, Germany). We discuss applications of the new setup to study deformation of a variety of materials, including ferropericlase, calcium perovskite, bridgmanite, and tantalum carbide, at high-pressure/temperature.

Intraplate volcanism originating from upwelling hydrous mantle transition zone

Most magmatism occurring on Earth is conventionally attributed to passive mantle upwelling at mid-ocean ridges, to slab devolatilization at subduction zones, or to mantle plumes. However, the widespread Cenozoic intraplate volcanism in northeast China^1-3 and the young petit-spot volcanoes^4-7 offshore of the Japan Trench cannot readily be associated with any of these mechanisms. In addition, the mantle beneath these types of volcanism is characterized by zones of anomalously low seismic velocity above and below the transition zone^8-12 (a mantle level located at depths between 410 and 660 kilometres). A comprehensive interpretation of these phenomena is lacking. Here we show that most (or possibly all) of the intraplate and petit-spot volcanism and low-velocity zones around the Japanese subduction zone can be explained by the Cenozoic interaction of the subducting Pacific slab with a hydrous mantle transition zone. Numerical modelling indicates that 0.2 to 0.3 weight per cent of water dissolved in mantle minerals that are driven out from the transition zone in response to subduction and retreat of a tectonic plate is sufficient to reproduce the observations. This suggests that a critical amount of water may have accumulated in the transition zone around this subduction zone, as well as in others of the Tethyan tectonic belt¹³ that are characterized by intraplate or petit-spot volcanism and low-velocity zones in the underlying mantle.

Influence of the Thickness of the Overriding Plate on Convergence Zone Dynamics

The important role played by the upper plate in convergence zones dynamics has long been underestimated but is now more and more emphasized. However, the influence of its thickness and/or strength on orogenic systems evolution remains largely unknown. Here we present results from 3D thermo-mechanical numerical simulations of convergence zones (including oceanic subduction followed by continental subduction/collision), in which we vary the rheological profile of the overriding plate (OP). For this, we systematically modify the crustal thickness of the overriding lithosphere and the temperature at the Moho to obtain a thermal thickness of the overriding lithosphere ranging from 80 to 180 km. While all models share a common global evolution (i.e., slab sinking, interaction between slab and the 660 km discontinuity, continental subduction/collision, and slab breakoff), they also highlight first-order differences arising from the variations in the OP strength (thermal thickness). With a thin/weak OP, slab rollback is favored, the slab dip is low, the mantle flow above the slab is vigorous, and the trench migrates at a high rate compared to a thick/strong OP. In addition, slab breakoff and back-arc basin formation events occur significantly earlier than in models involving a thick OP. Our models therefore highlight the major role played by the thickness/strength of the OP on convergence zone dynamics and illustrate its influence in a quantitative way.

Discharge of deeply rooted fluids from submarine mud volcanism in the Taiwan accretionary prism

Qualitative and quantitative assessments of fluid cycling are essential to address the role and transport of deeply sourced fluids in subduction systems. In this study, sediment cores distributed across a submarine mud volcano (SMV) offshore southwestern Taiwan were investigated to determine the characteristics of fluids generated through the convergence between the Eurasian and Phillippine Sea Plates. The low dissolved chloride concentration combined with the enrichment of ¹⁸O, and depletion of ²H of pore fluids suggest the discharge of deep freshwater formed by smectite dehydration at an equilibrium temperature of 100 to 150 °C. The upward fluid velocities, decreasing from 2.0 to 5.0 cm yr⁻¹ at the center to a negligible value at margin sites, varied with the rate and efficiency of anaerobic methanotrophy, demonstrating the impact of fluid migration on biogeochemical processes and carbon cycling. By extrapolating the velocity pattern, the flux of fluids exported from 13 SMVs into seawater amounted up to 1.3–2.5 × 10⁷ kg yr⁻¹, a quantity accounting for 1.1–28.6% of the smectite-bound water originally stored in the incoming sediments. Our results imply that SMVs could act as a conduit to channel the fluids produced from great depth/temperature into seafloor environments in a subduction system of the western Pacific Ocean.

Solid Earth change and the evolution of the Antarctic Ice Sheet

Recent studies suggest that Antarctica has the potential to contribute up to ~15 m of sea-level rise over the next few centuries. The evolution of the Antarctic Ice Sheet is driven by a combination of climate forcing and non-climatic feedbacks. In this review we focus on feedbacks between the Antarctic Ice Sheet and the solid Earth, and the role of these feedbacks in shaping the response of the ice sheet to past and future climate changes. The growth and decay of the Antarctic Ice Sheet reshapes the solid Earth via isostasy and erosion. In turn, the shape of the bed exerts a fundamental control on ice dynamics as well as the position of the grounding line—the location where ice starts to float. A complicating issue is the fact that Antarctica is situated on a region of the Earth that displays large spatial variations in rheological properties. These properties affect the timescale and strength of feedbacks between ice-sheet change and solid Earth deformation, and hence must be accounted for when considering the future evolution of the ice sheet.

The evolution of the Antarctic Ice Sheet is driven by a combination of climate forcing and non-climatic feedbacks. In this review, the authors focus on feedbacks between the Antarctic Ice Sheet and the solid Earth, and the role of these feedbacks in shaping the response of the ice sheet to past and future climate changes.

LIP formation and protracted lower mantle upwelling induced by rifting and delamination

Large Igneous Provinces (LIPs) are commonly attributed to mantle plumes, hot upwellings from the deep lower mantle, apparently unrelated to plate motions. However, LIPs often form in association with rifting and breakup. Using numerical modelling, we introduce a novel idea that explains plume-like mantle upwelling by plate tectonic processes. Our model indicates that rifting-induced delamination of orogenic lithosphere can perturb the thermochemical mantle stratification and induce lower mantle upwelling which causes syn-rift LIP formation followed by protracted and enhanced mid ocean ridge basalt (MORB) generation. Our model provides an explanation for the geographical correlation between the Caledonian suture, the North Atlantic Igneous Province (NAIP) and present-day Icelandic magmatism.

Crustal evolution and mantle dynamics through Earth history

Resolving the modes of mantle convection through Earth history, i.e. when plate tectonics started and what kind of mantle dynamics reigned before, is essential to the understanding of the evolution of the whole Earth system, because plate tectonics influences almost all aspects of modern geological processes. This is a challenging problem because plate tectonics continuously rejuvenates Earth's surface on a time scale of about 100 Myr, destroying evidence for its past operation. It thus becomes essential to exploit indirect evidence preserved in the buoyant continental crust, part of which has survived over billions of years. This contribution starts with an in-depth review of existing models for continental growth. Growth models proposed so far can be categorized into three types: crust-based, mantle-based and other less direct inferences, and the first two types are particularly important as their difference reflects the extent of crustal recycling, which can be related to subduction. Then, a theoretical basis for a change in the mode of mantle convection in the Precambrian is reviewed, along with a critical appraisal of some popular notions for early Earth dynamics. By combining available geological and geochemical observations with geodynamical considerations, a tentative hypothesis is presented for the evolution of mantle dynamics and its relation to surface environment; the early onset of plate tectonics and gradual mantle hydration are responsible not only for the formation of continental crust but also for its preservation as well as its emergence above sea level. Our current understanding of various material properties and elementary processes is still too premature to build a testable, quantitative model for this hypothesis, but such modelling efforts could potentially transform the nature of the data-starved early Earth research by quantifying the extent of preservation bias.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

When crust comes of age: on the chemical evolution of Archaean, felsic continental crust by crustal drip tectonics

The secular evolution of the Earth's crust is marked by a profound change in average crustal chemistry between 3.2 and 2.5 Ga. A key marker for this change is the transition from Archaean sodic granitoid intrusions of the tonalite-trondhjemite-granodiorite (TTG) series to potassic (K) granitic suites, akin (but not identical) to I-type granites that today are associated with subduction zones. It remains poorly constrained as to how and why this change was initiated and if it holds clues about the geodynamic transition from a pre-plate tectonic mode, often referred to as stagnant lid, to mobile plate tectonics. Here, we combine a series of proposed mechanisms for Archaean crustal geodynamics in a single model to explain the observed change in granitoid chemistry. Numeric modelling indicates that upper mantle convection drives crustal flow and subsidence, leading to profound diversity in lithospheric thickness with thin versus thick proto-plates. When convecting asthenospheric mantle interacts with lower lithosphere, scattered crustal drips are created. Under increasing P-T conditions, partial melting of hydrated meta-basalt within these drips produces felsic melts that intrude the overlying crust to form TTG. Dome structures, in which these melts can be preserved, are a positive diapiric expression of these negative drips. Transitional TTG with elevated K mark a second evolutionary stage, and are blends of subsided and remelted older TTG forming K-rich melts and new TTG melts. Ascending TTG-derived melts from asymmetric drips interact with the asthenospheric mantle to form hot, high-Mg sanukitoid. These melts are small in volume, predominantly underplated, and their heat triggered melting of lower crustal successions to form higher-K granites. Importantly, this evolution operates as a disseminated process in space and time over hundreds of millions of years (greater than 200 Ma) in all cratons. This focused ageing of the crust implies that compiled geochemical data can only broadly reflect geodynamic changes on a global or even craton-wide scale. The observed change in crustal chemistry does mark the lead up to but not the initiation of modern-style subduction.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

The dependence of planetary tectonics on mantle thermal state: applications to early Earth evolution

For plate tectonics to operate on a planet, mantle convective forces must be capable of forming weak, localized shear zones in the lithosphere that act as plate boundaries. Otherwise, a planet's mantle will convect in a stagnant lid regime, where subduction and plate motions are absent. Thus, when and how plate tectonics initiated on the Earth is intrinsically tied to the ability of mantle convection to form plate boundaries; however, the physics behind this process are still uncertain. Most mantle convection models have employed a simple pseudoplastic model of the lithosphere, where the lithosphere 'fails' and develops a mobile lid when stresses in the lithosphere reach the prescribed yield stress. With pseudoplasticity high mantle temperatures and high rates of internal heating, conditions relevant for the early Earth, impede plate boundary formation by decreasing lithospheric stresses, and hence favour a stagnant lid for the early Earth. However, when a model for shear zone formation based on grain size reduction is used, early Earth thermal conditions do not favour a stagnant lid. While lithosphere stress drops with increasing mantle temperature or heat production rate, the deformational work, which drives grain size reduction, increases. Thus, the ability of convection to form weak plate boundaries is not impeded by early Earth thermal conditions. However, mantle thermal state does change the style of subduction and lithosphere mobility; high mantle temperatures lead to a more sluggish, drip-like style of subduction. This 'sluggish lid' convection may be able to explain many of the key observations of early Earth crust formation processes preserved in the geologic record. Moreover, this work highlights the importance of understanding the microphysics of plate boundary formation for assessing early Earth tectonics, as different plate boundary formation mechanisms are influenced by mantle thermal state in fundamentally different ways.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

Strong plates enhance mantle mixing in early Earth

In the present-day Earth, some subducting plates (slabs) are flattening above the upper–lower mantle boundary at ~670 km depth, whereas others go through, indicating a mode between layered and whole-mantle convection. Previous models predicted that in a few hundred degree hotter early Earth, convection was likely more layered due to dominant slab stagnation. In self-consistent numerical models where slabs have a plate-like rheology, strong slabs and mobile plate boundaries favour stagnation for old and penetration for young slabs, as observed today. Here we show that such models predict slabs would have penetrated into the lower mantle more easily in a hotter Earth, when a weaker asthenosphere and decreased plate density and strength resulted in subduction almost without trench retreat. Thus, heat and material transport in the Earth’s mantle was more (rather than less) efficient in the past, which better matches the thermal evolution of the Earth.

The subducting plates can either penetrate straight into the lower mantle or flatten in the mantle transition zone, yet slab dynamics in the past remains unclear. Here, using subduction models, the authors predict that a hotter early Earth was probably more favourable to lower mantle slab penetration.

Carbon Cycling and Habitability of Earth-Sized Stagnant Lid Planets

Models of thermal evolution, crustal production, and CO₂ cycling are used to constrain the prospects for habitability of rocky planets, with Earth-like size and composition, in the stagnant lid regime. Specifically, we determine the conditions under which such planets can maintain rates of CO₂ degassing large enough to prevent global surface glaciation but small enough so as not to exceed the upper limit on weathering rates provided by the supply of fresh rock, a situation which would lead to runaway atmospheric CO₂ accumulation and an inhospitably hot climate. The models show that stagnant lid planets with initial radiogenic heating rates of 100-250 TW, and with total CO₂ budgets ranging from ∼10^-2 to 1 times Earth's estimated CO₂ budget, can maintain volcanic outgassing rates suitable for habitability for ≈1-5 Gyr; larger CO₂ budgets result in uninhabitably hot climates, while smaller budgets result in global glaciation. High radiogenic heat production rates favor habitability by sustaining volcanism and CO₂ outgassing longer. Thus, the results suggest that plate tectonics may not be required for establishing a long-term carbon cycle and maintaining a stable, habitable climate. The model is necessarily highly simplified, as the uncertainties with exoplanet thermal evolution and outgassing are large. Nevertheless, the results provide some first-order guidance for future exoplanet missions, by predicting the age at which habitability becomes unlikely for a stagnant lid planet as a function of initial radiogenic heat budget. This prediction is powerful because both planet heat budget and age can potentially be constrained from stellar observations. Key Words: Exoplanets-Habitability-Stagnant lid tectonics-Carbon cycle-Volcanism. Astrobiology 18, 873-896.

Transient rheology of the Sumatran mantle wedge revealed by a decade of great earthquakes

Understanding the rheological properties of the upper mantle is essential to develop a consistent model of mantle dynamics and plate tectonics. However, the spatial distribution and temporal evolution of these properties remain unclear. Here, we infer the rheological properties of the asthenosphere across multiple great megathrust earthquakes between 2004 and 2014 along the Sumatran subduction zone, taking advantage of decade-long continuous GPS and tide-gauge measurements. We observe transient mantle wedge flow following these earthquakes, and infer the temporal evolution of the effective viscosity. We show that the evolution of stress and strain rate following these earthquakes is better matched by a bi-viscous than by a power-law rheology model, and we estimate laterally heterogeneous transient and background viscosities on the order of ~10¹⁷ and ~10¹⁹ Pa s, respectively. Our results constitute a preliminary rheological model to explain stress evolution within earthquake cycles and the development of seismic hazard in the region.

The rheology of the upper mantle is key to understanding how plate tectonics may evolve. Here, using GPS and tide-gauge measurements along the Sumatran subduction zone, the authors’ show that a bi-viscous rheology model is needed to explain the stress and strain evolution of the upper mantle following earthquakes.

Characterization by Scanning Precession Electron Diffraction of an Aggregate of Bridgmanite and Ferropericlase Deformed at HP-HT

Scanning precession electron diffraction is an emerging promising technique for mapping phases and crystal orientations with short acquisition times (10-20 ms/pixel) in a transmission electron microscope similarly to the Electron Backscattered Diffraction (EBSD) or Transmission Kikuchi Diffraction (TKD) techniques in a scanning electron microscope. In this study, we apply this technique to the characterization of deformation microstructures in an aggregate of bridgmanite and ferropericlase deformed at 27 GPa and 2,130 K. Such a sample is challenging for microstructural characterization for two reasons: (i) the bridgmanite is very unstable under electron irradiation, (ii) under high stress conditions, the dislocation density is so large that standard characterization by diffraction contrast are limited, or impossible. Here we show that detailed analysis of intracrystalline misorientations sheds some light on the deformation mechanisms of both phases. In bridgmanite, deformation is accommodated by localized, amorphous, shear deformation lamellae whereas ferropericlase undergoes large strains leading to grain elongation in response to intense dislocation activity with no evidence for recrystallization. Plastic strain in ferropericlase can be semiquantitatively assessed by following kernel average misorientation distributions.