Mantle convection and plate tectonics: toward an integrated physical and chemical theory

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.

Thermodynamics of continental deformation

Continental deformation is known to be controlled by the interplay between tectonic and gravitational forces modulated by thermal relaxation-controlled lithospheric strength leading to oscillations around an equilibrium state, or to runaway extension. Using data-driven thermomechanical modelling of the Alpine Himalayan Collision Zone, we demonstrate how deviations from an equilibrium between mantle dynamics, plate-boundary forces, and the thermochemical configuration of the lithosphere control continental deformation. We quantify such balance between the internal energy of the plate and tectonic forces in terms of a critical crustal thickness, that match the global average of present-day continental crust. It follows that thicker intraplate domains than the critical crust (orogens) must undergo weakening due to their increased internal energy, and, in doing so, they dissipate the acquired energy within a diffused zone of deformation, unlike the localized deformation seen along plate boundaries. This evolution is controlled by a dissipative thermodynamic feedback loop between thermal and mechanical relaxation of the driving energy in the orogenic lithosphere. Exponentially growing energy states, leading to runaway extension are efficiently dampened by enhanced dissipation from radioactive heat sources. This ultimately drives orogens with their thickened radiogenic crust towards a final equilibrium state. Our results suggest a genetic link between the thermochemical state of the crust and the tectonic evolution of silicate Earth-like planets.

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.

The Role of N2 as a Geo-Biosignature for the Detection and Characterization of Earth-like Habitats

Since the Archean, N₂ has been a major atmospheric constituent in Earth's atmosphere. Nitrogen is an essential element in the building blocks of life; therefore, the geobiological nitrogen cycle is a fundamental factor in the long-term evolution of both Earth and Earth-like exoplanets. We discuss the development of Earth's N₂ atmosphere since the planet's formation and its relation with the geobiological cycle. Then we suggest atmospheric evolution scenarios and their possible interaction with life-forms: first for a stagnant-lid anoxic world, second for a tectonically active anoxic world, and third for an oxidized tectonically active world. Furthermore, we discuss a possible demise of present Earth's biosphere and its effects on the atmosphere. Since life-forms are the most efficient means for recycling deposited nitrogen back into the atmosphere at present, they sustain its surface partial pressure at high levels. Also, the simultaneous presence of significant N₂ and O₂ is chemically incompatible in an atmosphere over geological timescales. Thus, we argue that an N₂-dominated atmosphere in combination with O₂ on Earth-like planets within circumstellar habitable zones can be considered as a geo-biosignature. Terrestrial planets with such atmospheres will have an operating tectonic regime connected with an aerobic biosphere, whereas other scenarios in most cases end up with a CO₂-dominated atmosphere. We conclude with implications for the search for life on Earth-like exoplanets inside the habitable zones of M to K stars.

High-pressure, High-temperature Deformation Experiment Using the New Generation Griggs-type Apparatus

In order to address geological processes at great depths, rock deformation should ideally be tested at high pressure (> 0.5 GPa) and high temperature (> 300 °C). However, because of the low stress resolution of current solid-pressure-medium apparatuses, high-resolution measurements are today restricted to low-pressure deformation experiments in the gas-pressure-medium apparatus. A new generation of solid-medium piston-cylinder ("Griggs-type") apparatus is here described. Able to perform high-pressure deformation experiments up to 5 GPa and designed to adapt an internal load cell, such a new apparatus offers the potential to establish a technological basis for high-pressure rheology. This paper provides video-based detailed documentation of the procedure (using the "conventional" solid-salt assembly) to perform high-pressure, high-temperature experiments with the newly designed Griggs-type apparatus. A representative result of a Carrara marble sample deformed at 700 °C, 1.5 GPa and 10^-5 s^-1 with the new press is also given. The related stress-time curve illustrates all steps of a Griggs-type experiment, from increasing pressure and temperature to sample quenching when deformation is stopped. Together with future developments, the critical steps and limitations of the Griggs apparatus are then discussed.

Size effects resolve discrepancies in 40 years of work on low-temperature plasticity in olivine

The strength of olivine at low temperatures and high stresses in Earth's lithospheric mantle exerts a critical control on many geodynamic processes, including lithospheric flexure and the formation of plate boundaries. Unfortunately, laboratory-derived values of the strength of olivine at lithospheric conditions are highly variable and significantly disagree with those inferred from geophysical observations. We demonstrate via nanoindentation that the strength of olivine depends on the length scale of deformation, with experiments on smaller volumes of material exhibiting larger yield stresses. This "size effect" resolves discrepancies among previous measurements of olivine strength using other techniques. It also corroborates the most recent flow law for olivine, which proposes a much weaker lithospheric mantle than previously estimated, thus bringing experimental measurements into closer alignment with geophysical constraints. Further implications include an increased difficulty of activating plasticity in cold, fine-grained shear zones and an impact on the evolution of fault surface roughness due to the size-dependent deformation of nanometer- to micrometer-sized asperities.

The fate of water within Earth and super-Earths and implications for plate tectonics

The Earth is likely to have acquired most of its water during accretion. Internal heat of planetesimals by short-lived radioisotopes would have caused some water loss, but impacts into planetesimals were insufficiently energetic to produce further drying. Water is thought to be critical for the development of plate tectonics, because it lowers viscosities in the asthenosphere, enabling subduction. The following issue persists: if water is necessary for plate tectonics, but subduction itself hydrates the upper mantle, how is the upper mantle initially hydrated? The giant impacts of late accretion created magma lakes and oceans, which degassed during solidification to produce a heavy atmosphere. However, some water would have remained in the mantle, trapped within crystallographic defects in nominally anhydrous minerals. In this paper, we present models demonstrating that processes associated with magma ocean solidification and overturn may segregate sufficient quantities of water within the upper mantle to induce partial melting and produce a damp asthenosphere, thereby facilitating plate tectonics and, in turn, the habitability of Earth-like extrasolar planets.This article is part of the themed issue 'The origin, history and role of water in the evolution of the inner Solar System'.

Episodic entrainment of deep primordial mantle material into ocean island basalts

Chemical differences between mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs) provide critical evidence that the Earth's mantle is compositionally heterogeneous. MORBs generally exhibit a relatively low and narrow range of (3)He/(4)He ratios on a global scale, whereas OIBs display larger variability in both time and space. The primordial origin of (3)He in OIBs has motivated hypotheses that high (3)He/(4)He ratios are the product of mantle plumes sampling chemically distinct material, but do not account for lower MORB-like (3)He/(4)He ratios in OIBs, nor their observed spatial and temporal variability. Here we perform thermochemical convection calculations which show the variable (3)He/(4)He signature of OIBs can be reproduced by deep isolated mantle reservoirs of primordial material that are viscously entrained by thermal plumes. Entrainment is highly time-dependent, producing a wide range of (3)He/(4)He ratios similar to that observed in OIBs worldwide and indicate MORB-like (3)He/(4)He ratios in OIBs cannot be used to preclude deep mantle-sourced hotspots.

Spin crossover in ferropericlase and velocity heterogeneities in the lower mantle

Deciphering the origin of seismic velocity heterogeneities in the mantle is crucial to understanding internal structures and processes at work in the Earth. The spin crossover in iron in ferropericlase (Fp), the second most abundant phase in the lower mantle, introduces unfamiliar effects on seismic velocities. First-principles calculations indicate that anticorrelation between shear velocity (VS) and bulk sound velocity (Vφ) in the mantle, usually interpreted as compositional heterogeneity, can also be produced in homogeneous aggregates containing Fp. The spin crossover also suppresses thermally induced heterogeneity in longitudinal velocity (VP) at certain depths but not in VS. This effect is observed in tomography models at conditions where the spin crossover in Fp is expected in the lower mantle. In addition, the one-of-a-kind signature of this spin crossover in the RS/P (∂ ln VS/∂ ln VP) heterogeneity ratio might be a useful fingerprint to detect the presence of Fp in the lower mantle.

Pattern formation and coarsening dynamics in three-dimensional convective mixing in porous media

Geological carbon dioxide (CO2) sequestration entails capturing and injecting CO2 into deep saline aquifers for long-term storage. The injected CO2 partially dissolves in groundwater to form a mixture that is denser than the initial groundwater. The local increase in density triggers a gravitational instability at the boundary layer that further develops into columnar plumes of CO2-rich brine, a process that greatly accelerates solubility trapping of the CO2. Here, we investigate the pattern-formation aspects of convective mixing during geological CO2 sequestration by means of high-resolution three-dimensional simulation. We find that the CO2 concentration field self-organizes as a cellular network structure in the diffusive boundary layer at the top boundary. By studying the statistics of the cellular network, we identify various regimes of finger coarsening over time, the existence of a non-equilibrium stationary state, and a universal scaling of three-dimensional convective mixing.

Search for magnetic monopoles in polar volcanic rocks

For a broad range of values of magnetic monopole mass and charge, the abundance of monopoles trapped inside Earth would be expected to be enhanced in the mantle beneath the geomagnetic poles. A search for magnetic monopoles was conducted using the signature of an induced persistent current following the passage of igneous rock samples through a SQUID-based magnetometer. A total of 24.6 kg of rocks from various selected sites, among which 23.4 kg are mantle-derived rocks from the Arctic and Antarctic areas, was analyzed. No monopoles were found, and a 90% confidence level upper limit of 9.8 × 10(-5)/g is set on the monopole density in the search samples.

Basics of lava-lamp convection

Laboratory experiments are reported in an immiscible two-fluid system, where thermal convection is initiated by heating at the bottom and cooling at the top. The lava-lamp regime is characterized by a robust periodic exchange process where warm blobs rise from the bottom, attach to the top surface for a while, then cold blobs sink down again. Immiscibility allows to reach real steady (dynamical equilibrium) states which can be sustained for several days. Two modes of lava-lamp convection could be identified by recording and evaluating temperature time series at the bottom and at the top of the container: a "slow" mode is determined by an effective heat transport speed at a given temperature gradient, while a second mode of constant periodicity is viscosity limited. Contrasting of laboratory and geophysical observations yields the conclusion that the frequently suggested lava-lamp analogy fails for the accepted models of mantle convection.

Palaeoproterozoic supercontinents and global evolution: correlations from core to atmosphere

The Palaeoproterozoic era was a time of profound change in Earth evolution and represented perhaps the first supercontinent cycle, from the amalgamation and dispersal of a possible Neoarchaean supercontinent to the formation of the 1.9–1.8 Ga supercontinent Nuna. This supercontinent cycle, although currently lacking in palaeogeographic detail, can in principle provide a contextual framework to investigate the relationships between deep-Earth and surface processes. In this article, we graphically summarize secular evolution from the Earth's core to its atmosphere, from the Neoarchaean to the Mesoproterozoic eras (specifically 3.0–1.2 Ga), to reveal intriguing temporal relationships across the various ‘spheres’ of the Earth system. At the broadest level our compilation confirms an important deep-Earth event at c. 2.7 Ga that is manifested in an abrupt increase in geodynamo palaeointensity, a peak in the global record of large igneous provinces, and a broad maximum in several mantle-depletion proxies. Temporal coincidence with juvenile continental crust production and orogenic gold, massive-sulphide and porphyry copper deposits, indicate enhanced mantle convection linked to a series of mantle plumes and/or slab avalanches. The subsequent stabilization of cratonic lithosphere, the possible development of Earth's first supercontinent and the emergence of the continents led to a changing surface environment in which voluminous banded iron-formations could accumulate on the continental margins and photosynthetic life could flourish. This in turn led to irreversible atmospheric oxidation at 2.4–2.3 Ga, extreme events in global carbon cycling, and the possible dissipation of a former methane greenhouse atmosphere that resulted in extensive Palaeoproterozoic ice ages. Following the great oxidation event, shallow marine sulphate levels rose, sediment-hosted and iron-oxide-rich metal deposits became abundant, and the transition to sulphide-stratified oceans provided the environment for early eukaryotic evolution. Recent advances in the geochronology of the global stratigraphic record have made these inferences possible. Frontiers for future research include more refined modelling of Earth's thermal and geodynamic evolution, palaeomagnetic studies of geodynamo intensity and continental motions, further geochronology and tectonic syntheses at regional levels, development of new isotopic systems to constrain geochemical cycles, and continued innovation in the search for records of early life in relation to changing palaeoenvironments.

What CO2 well gases tell us about the origin of noble gases in the mantle and their relationship to the atmosphere

Study of commercially produced volcanic CO2 gas associated with the Colorado Plateau, USA, has revealed substantial new information about the noble gas isotopic composition and elemental abundance pattern of the mantle. Combined with published data from mid-ocean ridge basalts, it is now clear that the convecting mantle has a maximum (20)Ne/(22)Ne isotopic composition, indistinguishable from that attributed to solar wind-implanted (SWI) neon in meteorites. This is distinct from the higher (20)Ne/(22)Ne isotopic value expected for solar nebula gases. The non-radiogenic xenon isotopic composition of the well gases shows that 20 per cent of the mantle Xe is 'solar-like' in origin, but cannot resolve the small isotopic difference between the trapped meteorite 'Q'-component and solar Xe. The mantle primordial (20)Ne/(132)Xe is approximately 1400 and is comparable with the upper end of that observed in meteorites. Previous work using the terrestrial (129)I - (129)Xe mass balance demands that almost 99 per cent of the Xe (and therefore other noble gases) has been lost from the accreting solids and that Pu-I closure age models have shown this to have occurred in the first ca 100Ma of the Earth's history. The highest concentrations of Q-Xe and solar wind-implanted (SWI)-Ne measured in meteorites allow for this loss and these high-abundance samples have a Ne/Xe ratio range compatible with the 'recycled-air-corrected' terrestrial mantle. These observations do not support models in which the terrestrial mantle acquired its volatiles from the primary capture of solar nebula gases and, in turn, strongly suggest that the primary terrestrial atmosphere, before isotopic fractionation, is most probably derived from degassed trapped volatiles in accreting material.By contrast, the non-radiogenic argon, krypton and 80 per cent of the xenon in the convecting mantle have the same isotopic composition and elemental abundance pattern as that found in seawater with a small sedimentary Kr and Xe admix. These mantle heavy noble gases are dominated by recycling of air dissolved in seawater back into the mantle. Numerical simulations suggest that plumes sampling the core-mantle boundary would be enriched in seawater-derived noble gases compared with the convecting mantle, and therefore have substantially lower (40)Ar/(36)Ar. This is compatible with observation. The subduction process is not a complete barrier to volatile return to the mantle.

Reduced Radiative Conductivity of Low-Spin (Mg,Fe)O in the Lower Mantle

Optical absorption spectra have been measured at pressures up to 80 gigapascals (GPa) for the lower-mantle oxide magnesiowüstite (Mg,Fe)O. Upon reaching the high-spin to low-spin transition of Fe2+ at about 60 GPa, we observed enhanced absorption in the mid- and near-infrared spectral range, whereas absorption in the visible-ultraviolet was reduced. The observed changes in absorption are in contrast to prediction and are attributed to d-d orbital charge transfer in the Fe2+ ion. The results indicate that low-spin (Mg,Fe)O will exhibit lower radiative thermal conductivity than high-spin (Mg,Fe)O, which needs to be considered in future geodynamic models of convection and plume stabilization in the lower mantle.

Thermochemical structures beneath Africa and the Pacific Ocean

Large low-velocity seismic anomalies have been detected in the Earth's lower mantle beneath Africa and the Pacific Ocean that are not easily explained by temperature variations alone. The African anomaly has been interpreted to be a northwest-southeast-trending structure with a sharp-edged linear, ridge-like morphology. The Pacific anomaly, on the other hand, appears to be more rounded in shape. Mantle models with heterogeneous composition have related these structures to dense thermochemical piles or superplumes. It has not been shown, however, that such models can lead to thermochemical structures that satisfy the geometrical constraints, as inferred from seismological observations. Here we present numerical models of thermochemical convection in a three-dimensional spherical geometry using plate velocities inferred for the past 119 million years. We show that Earth's subduction history can lead to thermochemical structures similar in shape to the observed large, lower-mantle velocity anomalies. We find that subduction history tends to focus dense material into a ridge-like pile beneath Africa and a relatively more-rounded pile under the Pacific Ocean, consistent with seismic observations.

Forced chemical mixing in alloys driven by plastic deformation

Molecular dynamics simulations of forced atomic mixing in crystalline binary alloys during plastic deformation at 100 K are performed. Nearly complete atomic mixing is observed in systems that have a large positive heat mixing and in systems with a large lattice mismatch. Only systems that contained a hard precipitate in a soft matrix do not mix. The amount of mixing is quantified by defining a mean square relative displacement of pairs of atoms, sigma(2)(R,t), that were initially separated by a distance R. Analysis of sigma(2)(R,t) and visual inspection of the displacement fields reveal that forced mixing results from dislocation glide, and that it resembles the forced mixing of a substance advected by a turbulent flow. Consideration of sigma(2)(R,t) also provides a rationalization of compositional self-organization during plastic deformation at higher temperatures.

Electronic Transitions in Perovskite: Possible Nonconvecting Layers in the Lower Mantle

We measured the spin state of iron in magnesium silicate perovskite (Mg(0.9),Fe(0.1))SiO(3) at high pressure and found two electronic transitions occurring at 70 gigapascals and at 120 gigapascals, corresponding to partial and full electron pairing in iron, respectively. The proportion of iron in the low spin state thus grows with depth, increasing the transparency of the mantle in the infrared region, with a maximum at pressures consistent with the D" layer above the core-mantle boundary. The resulting increase in radiative thermal conductivity suggests the existence of nonconvecting layers in the lowermost mantle.

Finite-frequency tomography reveals a variety of plumes in the mantle

We present tomographic evidence for the existence of deep-mantle thermal convection plumes. P-wave velocity images show at least six well-resolved plumes that extend into the lowermost mantle: Ascension, Azores, Canary, Easter, Samoa, and Tahiti. Other less well-resolved plumes, including Hawaii, may also reach the lowermost mantle. We also see several plumes that are mostly confined to the upper mantle, suggesting that convection may be partially separated into two depth regimes. All of the observed plumes have diameters of several hundred kilometers, indicating that plumes convey a substantial fraction of the internal heat escaping from Earth.

Geodynamic and seismic constraints on the thermochemical structure and dynamics of convection in the deep mantle

We revisit a recent study by Forte & Mitrovica in which global geophysical observables associated with mantle convection were inverted and the existence of a strong increase in viscosity near a depth of 2000 km was inferred. Employing mineral-physics data and theory we also showed that, although there are chemical anomalies in the lowermost mantle, they are unable to inhibit the dominant thermal buoyancy of the deep-mantle mega-plumes below Africa and the Pacific Ocean. New Monte Carlo simulations are employed to explore the impact of uncertainties in current mineral-physics constraints on inferences of deep-mantle thermochemical structure. To explore the impact of the high-viscosity peak at a depth of 2000 km on the evolution of lower-mantle structure, we carried out time-dependent convection simulations. The latter show that the stability and longevity of the dominant long-wavelength heterogeneity in the lowermost mantle are controlled by this viscosity peak.

Zoned mantle convection

We review the present state of our understanding of mantle convection with respect to geochemical and geophysical evidence and we suggest a model for mantle convection and its evolution over the Earth's history that can reconcile this evidence. Whole-mantle convection, even with material segregated within the D" region just above the core-mantle boundary, is incompatible with the budget of argon and helium and with the inventory of heat sources required by the thermal evolution of the Earth. We show that the deep-mantle composition in lithophilic incompatible elements is inconsistent with the storage of old plates of ordinary oceanic lithosphere, i.e. with the concept of a plate graveyard. Isotopic inventories indicate that the deep-mantle composition is not correctly accounted for by continental debris, primitive material or subducted slabs containing normal oceanic crust. Seismological observations have begun to hint at compositional heterogeneity in the bottom 1000 km or so of the mantle, but there is no compelling evidence in support of an interface between deep and shallow mantle at mid-depth. We suggest that in a system of thermochemical convection, lithospheric plates subduct to a depth that depends - in a complicated fashion - on their composition and thermal structure. The thermal structure of the sinking plates is primarily determined by the direction and rate of convergence, the age of the lithosphere at the trench, the sinking rate and the variation of these parameters over time (i.e. plate-tectonic history) and is not the same for all subduction systems. The sinking rate in the mantle is determined by a combination of thermal (negative) and compositional buoyancy and as regards the latter we consider in particular the effect of the loading of plates with basaltic plateaux produced by plume heads. Barren oceanic plates are relatively buoyant and may be recycled preferentially in the shallow mantle. Oceanic plateau-laden plates have a more pronounced negative buoyancy and can more easily founder to the very base of the mantle. Plateau segregation remains statistical and no sharp compositional interface is expected from the multiple fate of the plates. We show that the variable depth subduction of heavily laden plates can prevent full vertical mixing and preserve a vertical concentration gradient in the mantle. In addition, it can account for the preservation of scattered remnants of primitive material in the deep mantle and therefore for the Ar and (3)He observations in ocean-island basalts.