Evidence for deep mantle circulation from global tomography

Gondwanan flood basalts linked seismically to plume-induced lithosphere instability

Delamination of the continental lithospheric mantle is well recorded beneath several continents. However, the fate of the removed continental lithosphere has been rarely noted, unlike subducted slabs reasonably well imaged in the upper and mid mantle. Beneath former Gondwana, recent seismic tomographic models indicate the presence of at least 5 horizontal fast-wavespeed anomalies at ~600 km depths that do not appear to be related to slab subduction, including fast structures in locations consistent with delamination associated with the Paraná Flood Basalt event at ~134 Ma and the Deccan Traps event at ~66 Ma. These fast-wavespeed anomalies often lie above broad slow seismic wavespeed trunks at 500 to 700 km depths beneath former Gondwana, with slow wavespeed anomalies branching around them. Numerical experiments indicate that delaminated lithosphere tends to stagnate in the transition zone and mid-mantle above a mantle plume where it shapes subsequent plume upwelling. For hot plumes, the melt volume generated during plume-influenced delamination can easily reach ~2 to 4 × 106 km3, consistent with the basalt eruption volume at the Deccan Traps. This seismic and numerical evidence suggests that observed high-wavespeed mid-mantle anomalies beneath the locations of former flood basalts are delaminated fragments of former continental lithosphere, and that lithospheric delamination events in the presence of subcontinental plumes induced several of the continental flood basalts associated with the multiple breakup stages of Gondwanaland. Continued upwelling in these plumes can also have entrained subcontinental lithosphere in the mid-mantle to bring its distinctive geochemical signal to the modern mid-ocean spreading centers that surround southern and western Africa.

Fragmentation of continental subduction is ending the Himalayan orogeny

After two continents collide, plate convergence and orogenesis are sustained because subducted continental lithosphere continues pulling the surface plate. It remains controversial how, why, and when continental plate convergence and collision slow down and eventually cease. We use an unprecedented data coverage and present a regional-scale seismic tomographic image of the mantle structure beneath the Tibetan Plateau. In the mantle transition zone, we identify multiple high-velocity anomalies and interpret them as detached pieces of the Indian continental slab. Facilitated by internal heterogeneity of the continental lithosphere, piecewise slab detachments could reduce the slab pull force, resulting in the Miocene slowdown of the India-Eurasia convergence and coeval diachronous potassic volcanism in southern Tibet. We propose that slab detachment is a mechanism that eventually will lead to the end of the Indo-Eurasian continental collision and the Himalayan orogeny.

Earth's evolving geodynamic regime recorded by titanium isotopes

Earth's mantle has a two-layered structure, with the upper and lower mantle domains separated by a seismic discontinuity at about 660 km (refs. 1,2). The extent of mass transfer between these mantle domains throughout Earth's history is, however, poorly understood. Continental crust extraction results in Ti-stable isotopic fractionation, producing isotopically light melting residues3-7. Mantle recycling of these components can impart Ti isotope variability that is trackable in deep time. We report ultrahigh-precision 49Ti/47Ti ratios for chondrites, ancient terrestrial mantle-derived lavas ranging from 3.8 to 2.0 billion years ago (Ga) and modern ocean island basalts (OIBs). Our new Ti bulk silicate Earth (BSE) estimate based on chondrites is 0.052 ± 0.006‰ heavier than the modern upper mantle sampled by normal mid-ocean ridge basalts (N-MORBs). The 49Ti/47Ti ratio of Earth's upper mantle was chondritic before 3.5 Ga and evolved to a N-MORB-like composition between approximately 3.5 and 2.7 Ga, establishing that more continental crust was extracted during this epoch. The +0.052 ± 0.006‰ offset between BSE and N-MORBs requires that <30% of Earth's mantle equilibrated with recycled crustal material, implying limited mass exchange between the upper and lower mantle and, therefore, preservation of a primordial lower-mantle reservoir for most of Earth's geologic history. Modern OIBs record variable 49Ti/47Ti ratios ranging from chondritic to N-MORBs compositions, indicating continuing disruption of Earth's primordial mantle. Thus, modern-style plate tectonics with high mass transfer between the upper and lower mantle only represents a recent feature of Earth's history.

Seismic evidence for global basalt accumulation in the mantle transition zone

The mantle's compositional structure reflects the thermochemical evolution of Earth. Yet, even the radial average composition of the mantle remains debated. Here, we analyze a global dataset of shear and compressional waves reflecting off the 410- and 660-km discontinuities that is 10 times larger than any previous studies. Our array analysis retrieves globally averaged amplitude-distance trends in SS and PP precursor reflectivity from which we infer relative wavespeed and density contrasts and associated mantle composition. Our results are best matched by a basalt-enriched mantle transition zone, with higher basalt fractions near 660 (~40%) than 410 (~18-31%). These are consistent with mantle-convection/plate-recycling simulations, which predict that basaltic crust accumulates in the mantle transition zone, with basalt fractions peaking near the 660. Basalt segregation in the mantle transition zone also implies that the overall mantle is more silica enriched than the often-assumed pyrolitic mantle reference composition.

A neural network based global traveltime function (GlobeNN)

Global traveltime modeling is an essential component of modern seismological studies with a whole gamut of applications ranging from earthquake source localization to seismic velocity inversion. Emerging acquisition technologies like distributed acoustic sensing (DAS) promise a new era of seismological discovery by allowing a high-density of seismic observations. Conventional traveltime computation algorithms are unable to handle virtually millions of receivers made available by DAS arrays. Therefore, we develop GlobeNN-a neural network based traveltime function that can provide seismic traveltimes obtained from the cached realistic 3-D Earth model. We train a neural network to estimate the traveltime between any two points in the global mantle Earth model by imposing the validity of the eikonal equation through the loss function. The traveltime gradients in the loss function are computed efficiently using automatic differentiation, while the P-wave velocity is obtained from the vertically polarized P-wave velocity of the GLAD-M25 model. The network is trained using a random selection of source and receiver pairs from within the computational domain. Once trained, the neural network produces traveltimes rapidly at the global scale through a single evaluation of the network. As a byproduct of the training process, we obtain a neural network that learns the underlying velocity model and, therefore, can be used as an efficient storage mechanism for the huge 3-D Earth velocity model. These exciting features make our proposed neural network based global traveltime computation method an indispensable tool for the next generation of seismological advances.

The combined Hf and Nd isotope evolution of the depleted mantle requires Hadean continental formation

The onset and rates of continental growth are first-order indicators of early Earth dynamics, and whether substantial crust existed in the Hadean or much later has long been debated. Here, we present a theoretical analysis of published Hf and Nd isotopic data representing the depleted mantle and demonstrate that continental growth must have started in the early Hadean. Whereas the traditional interpretation of depleted mantle signatures in crustal rocks assumes unrealistic instantaneous mantle mixing, our modeling incorporates the effect of a finite mixing time over which these signatures are recorded in rocks produced through mantle melting. This effect is shown to delay, by as much as 0.65 to 0.75 billion years, the appearance of the earliest depleted mantle signatures in continental crust. Our results suggest that published observations of εHf, ε¹⁴³Nd, and μ¹⁴²Nd require Hadean growth of continental crust, with a minimum of 50% of today's continental volume already existing by the end of Hadean.

Primitive noble gases sampled from ocean island basalts cannot be from the Earth's core

Noble gas isotopes in plumes require a source of primitive volatiles largely isolated in the Earth for 4.5 Gyrs. Among the proposed reservoirs, the core is gaining interest in the absence of robust geochemical and geophysical evidence for a mantle source. This is supported by partitioning data showing that sufficient He and Ne could have been incorporated into the core to source plumes today. Here we perform ab initio calculations on the partitioning of He, Ne, Ar, Kr and Xe between liquid iron and silicate melt under core forming conditions. For He our results are consistent with previous studies allowing for substantial amounts of He in the core. In contrast, the partition coefficient for Ne is three orders of magnitude lower than He. This very low partition coefficient would result in a ³He/²²Ne ratio of ~10³ in the core, far higher than observed in ocean island basalts (OIBs). We conclude that the core is not the source of noble gases in OIBs.

Depressed 660-km discontinuity caused by akimotoite-bridgmanite transition

The 660-kilometre seismic discontinuity is the boundary between the Earth’s lower mantle and transition zone and is commonly interpreted as being due to the dissociation of ringwoodite to bridgmanite plus ferropericlase (post-spinel transition)^1–3. A distinct feature of the 660-kilometre discontinuity is its depression to 750 kilometres beneath subduction zones^4–10. However, in situ X-ray diffraction studies using multi-anvil techniques have demonstrated negative but gentle Clapeyron slopes (that is,  the ratio between pressure and temperature changes) of the post-spinel transition that do not allow a significant depression^11–13. On the other hand, conventional high-pressure experiments face difficulties in accurate phase identification due to inevitable pressure changes during heating and the persistent presence of metastable phases^1,3. Here we determine the post-spinel and akimotoite–bridgmanite transition boundaries by multi-anvil experiments using in situ X-ray diffraction, with the boundaries strictly based on the definition of phase equilibrium. The post-spinel boundary has almost no temperature dependence, whereas the akimotoite–bridgmanite transition has a very steep negative boundary slope at temperatures lower than ambient mantle geotherms. The large depressions of the 660-kilometre discontinuity in cold subduction zones are thus interpreted as the akimotoite–bridgmanite transition. The steep negative boundary of the akimotoite–bridgmanite transition will cause slab stagnation (a stalling of the slab’s descent) due to significant upward buoyancy^14,15.

X-ray diffraction experiments indicate that the depression of the Earth’s 660-kilometre seismic discontinuity beneath cold subduction zones is caused by a phase transition from akimotoite to bridgmanite, leading to slab stagnation.

Structural Diversity of Magnetite and Products of Its Decomposition at Extreme Conditions

Magnetite, Fe₃O₄, is the oldest known magnetic mineral and archetypal mixed-valence oxide. Despite its recognized role in deep Earth processes, the behavior of magnetite at extreme high-pressure high-temperature (HPHT) conditions remains insufficiently studied. Here, we report on single-crystal synchrotron X-ray diffraction experiments up to ∼80 GPa and 5000 K in diamond anvil cells, which reveal two previously unknown Fe₃O₄ polymorphs, γ-Fe₃O₄ with the orthorhombic Yb₃S₄-type structure and δ-Fe₃O₄ with the modified Th₃P₄-type structure. The latter has never been predicted for iron compounds. The decomposition of Fe₃O₄ at HPHT conditions was found to result in the formation of exotic phases, Fe₅O₇ and Fe₂₅O₃₂, with complex structures. Crystal-chemical analysis of iron oxides suggests the high-spin to low-spin crossover in octahedrally coordinated Fe³⁺ in the pressure interval between 43 and 51 GPa. Our experiments demonstrate that HPHT conditions promote the formation of ferric-rich Fe-O compounds, thus arguing for the possible involvement of magnetite in the deep oxygen cycle.

The evolution of basal mantle structure in response to supercontinent aggregation and dispersal

Seismic studies have revealed two Large Low-Shear Velocity Provinces (LLSVPs) in the lowermost mantle. Whether these structures remain stable over time or evolve through supercontinent cycles is debated. Here we analyze a recently published mantle flow model constrained by a synthetic plate motion model extending back to one billion years ago, to investigate how the mantle evolves in response to changing plate configurations. Our model predicts that sinking slabs segment the basal thermochemical structure below an assembling supercontinent, and that this structure eventually becomes unified due to slab push from circum-supercontinental subduction. In contrast, the basal thermochemical structure below the superocean is generally coherent due to the persistence of a superocean in our imposed plate reconstruction. The two antipodal basal thermochemical structures exchange material several times when part of one of the structures is carved out and merged with the other one, similarly to “exotic” tectonic terranes. Plumes mostly rise from thick basal thermochemical structures and in some instances migrate from the edges towards the interior of basal thermochemical structures due to slab push. Our results suggest that the topography of basal structures and distribution of plumes change over time due to the changing subduction network over supercontinent cycles.

Seismological expression of the iron spin crossover in ferropericlase in the Earth's lower mantle

The two most abundant minerals in the Earth’s lower mantle are bridgmanite and ferropericlase. The bulk modulus of ferropericlase (Fp) softens as iron d-electrons transition from a high-spin to low-spin state, affecting the seismic compressional velocity but not the shear velocity. Here, we identify a seismological expression of the iron spin crossover in fast regions associated with cold Fp-rich subducted oceanic lithosphere: the relative abundance of fast velocities in P- and S-wave tomography models diverges in the ~1,400-2,000 km depth range. This is consistent with a reduced temperature sensitivity of P-waves throughout the iron spin crossover. A similar signal is also found in seismically slow regions below ~1,800 km, consistent with broadening and deepening of the crossover at higher temperatures. The corresponding inflection in P-wave velocity is not yet observed in 1-D seismic profiles, suggesting that the lower mantle is composed of non-uniformly distributed thermochemical heterogeneities which dampen the global signature of the Fp spin crossover.

This study identifies the predicted seismic expression of the high-to-low iron spin crossover in the deep Earth mineral ferropericlase. A depth-dependent signal is detected in the fastest and slowest regions, related to lateral temperature variations, of several global seismic tomography models.

Segregated oceanic crust trapped at the bottom mantle transition zone revealed from ambient noise interferometry

The recycling of oceanic crust, with distinct isotopic and chemical signature from the pyrolite mantle, plays a critical role in the chemical evolution of the Earth with insights into mantle circulation. However, the role of the mantle transition zone during this recycling remains ambiguous. We here combine the unique resolution reflected body waves (P410P and P660P) retrieved from ambient noise interferometry with mineral physics modeling, to shed new light on transition zone physics. Our joint analysis reveals a generally sharp 660-km discontinuity and the existence of a localized accumulation of oceanic crust at the bottom mantle transition zone just ahead of the stagnant Pacific slab. The basalt accumulation is plausibly derived from the segregation of oceanic crust and depleted mantle of the adjacent stagnant slab. Our findings provide direct evidence of segregated oceanic crust trapped within the mantle transition zone and new insights into imperfect whole mantle circulation.

Seismic evidence for subduction-induced mantle flows underneath Middle America

Laboratory experiments and geodynamic simulations demonstrate that poloidal- and toroidal-mode mantle flows develop around subduction zones. Here, we use a new 3-D azimuthal anisotropy model constructed by full waveform inversion, to infer deep subduction-induced mantle flows underneath Middle America. At depths shallower than 150 km, poloidal-mode flow is perpendicular to the trajectory of the Middle American Trench. From 300 to 450 km depth, return flows surround the edges of the Rivera and Atlantic slabs, while escape flows are inferred through slab windows beneath Panama and central Mexico. Furthermore, at 700 km depth, the study region is dominated by the Farallon anomaly, with fast axes perpendicular to its strike, suggesting the development of lattice-preferred orientations by substantial stress. These observations provide depth-dependent seismic anisotropy for future mantle flow simulations, and call for further investigations about the deformation mechanisms and elasticity of minerals in the transition zone and uppermost lower mantle.

The motions of subducted slabs are expected to drive mantle flow around slab edges, however, evidence of deep mantle flow has so far remained elusive. Here, the authors present a Full Waveform Inversion 3-D anisotropy model which allows them to infer deep subduction-induced mantle flows underneath the Mid-Americas and the Caribbean.

Kimberlites reveal 2.5-billion-year evolution of a deep, isolated mantle reservoir

The widely accepted paradigm of Earth's geochemical evolution states that the successive extraction of melts from the mantle over the past 4.5 billion years formed the continental crust, and produced at least one complementary melt-depleted reservoir that is now recognized as the upper-mantle source of mid-ocean-ridge basalts¹. However, geochemical modelling and the occurrence of high ³He/⁴He (that is, primordial) signatures in some volcanic rocks suggest that volumes of relatively undifferentiated mantle may reside in deeper, isolated regions². Some basalts from large igneous provinces may provide temporally restricted glimpses of the most primitive parts of the mantle^3,4, but key questions regarding the longevity of such sources on planetary timescales-and whether any survive today-remain unresolved. Kimberlites, small-volume volcanic rocks that are the source of most diamonds, offer rare insights into aspects of the composition of the Earth's deep mantle. The radiogenic isotope ratios of kimberlites of different ages enable us to map the evolution of this domain through time. Here we show that globally distributed kimberlites originate from a single homogeneous reservoir with an isotopic composition that is indicative of a uniform and pristine mantle source, which evolved in isolation over at least 2.5 billion years of Earth history-to our knowledge, the only such reservoir that has been identified to date. Around 200 million years ago, extensive volumes of the same source were perturbed, probably as a result of contamination by exogenic material. The distribution of affected kimberlites suggests that this event may be related to subduction along the margin of the Pangaea supercontinent. These results reveal a long-lived and globally extensive mantle reservoir that underwent subsequent disruption, possibly heralding a marked change to large-scale mantle-mixing regimes. These processes may explain why uncontaminated primordial mantle is so difficult to identify in recent mantle-derived melts.

Primordial and recycled helium isotope signatures in the mantle transition zone

Isotope compositions of basalts provide information about the chemical reservoirs in Earth's interior and play a critical role in defining models of Earth's structure. However, the helium isotope signature of the mantle below depths of a few hundred kilometers has been difficult to measure directly. This information is a vital baseline for understanding helium isotopes in erupted basalts. We measured He-Sr-Pb isotope ratios in superdeep diamond fluid inclusions from the transition zone (depth of 410 to 660 kilometers) unaffected by degassing and shallow crustal contamination. We found extreme He-C-Pb-Sr isotope variability, with high ³He/⁴He ratios related to higher helium concentrations. This indicates that a less degassed, high-³He/⁴He deep mantle source infiltrates the transition zone, where it interacts with recycled material, creating the diverse compositions recorded in ocean island basalts.

Slab morphology and deformation beneath Izu-Bonin

Seismic tomography provides unique constraints on the morphology, the deformation, and (indirectly) the rheology of subducting slabs. We use teleseismic double-difference P-wave tomography to image with unprecedented clarity the structural complexity of the Izu-Bonin slab. We resolve a tear in the slab in the mantle transition zone (MTZ) between 26.5° N and 28° N. North of the tear, the slab is folded in the MTZ. Immediately above the fold hinge, a zone of reduced P-wavespeed may result from viscous dissipation within an incipient shear zone. To the south of the tear, the slab overturns and lies flat at the base of the MTZ. The ~680 km deep 2015 Bonin earthquake (Mw~7.9) is located at the northernmost edge of the overturning part of the slab. The localised tearing, shearing and buckling of the Izu-Bonin slab indicates that it remains highly viscous throughout the upper mantle and transition zone.

In the 1000 km long Izu-Bonin subduction zone to the south of Tokyo, the Pacific Plate descends beneath the Philippine Sea Plate. Here the authors use teleseismic double-difference tomography to image the complex morphology of the Izu-Bonin slab, especially in the mantle transition zone.

Inferring Earth's discontinuous chemical layering from the 660-kilometer boundary topography

Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth's upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle.

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

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

Insights on Upper Mantle Melting, Rheology, and Anelastic Behavior From Seismic Shear Wave Tomography

In seismic tomography we map the wave speed structure inside the Earth, but we ultimately seek to interpret those images in terms of physical parameters. This is challenging because many parameters can trade-off with each other to produce a given wave speed. The problem is compounded by the convention of mapping seismic structures as perturbations relative to a 1-D reference model, rather than absolute wave speeds. Using a full waveform tomography model of Europe as a case study, we quantify the extent to which thermochemical and dynamic properties can be constrained using only S wave speed, expressed in absolute values. The wave speed distributions of this tomography model are compared with 4 million thermochemical models, whose seismic properties are computed via thermodynamic modeling. These models sample the full range of realistic mantle compositions, including variable water and melt contents, and mineral intrinsic anelasticity is taken into account. Intrinsic anelasticity causes waves to travel more slowly at higher temperatures, leading to seismic attenuation, but the sensitivity of the wave speed reduction to temperature is, in turn, controlled by the wave frequency. Global studies of surface waves indicate an anticorrelation between S wave speed and attenuation. We therefore only retain thermochemical models satisfying this anticorrelation. Our study indicates that the frequency dependence of anelasticity, α, depends on temperature or rheology, with α ≈ 0.1 being most appropriate in cold or lithospheric mantle and α ≈ 0.3 in warmer regions (i.e., the asthenosphere). Additionally, the slowest regions require specific compositions and/or a velocity-weakening mechanism, such as partial melting, elastically accommodated grain boundary sliding, or water.

On the consistency of seismically imaged lower mantle slabs

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000-1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.

Stability of ferrous-iron-rich bridgmanite under reducing midmantle conditions

Our current understanding of the electronic state of iron in lower-mantle minerals leads to a considerable disagreement in bulk sound speed with seismic measurements if the lower mantle has the same composition as the upper mantle (pyrolite). In the modeling studies, the content and oxidation state of Fe in the minerals have been assumed to be constant throughout the lower mantle. Here, we report high-pressure experimental results in which Fe becomes dominantly Fe²⁺ in bridgmanite synthesized at 40-70 GPa and 2,000 K, while it is in mixed oxidation state (Fe³⁺/∑Fe = 60%) in the samples synthesized below and above the pressure range. Little Fe³⁺ in bridgmanite combined with the strong partitioning of Fe²⁺ into ferropericlase will alter the Fe content for these minerals at 1,100- to 1,700-km depths. Our calculations show that the change in iron content harmonizes the bulk sound speed of pyrolite with the seismic values in this region. Our experiments support no significant changes in bulk composition for most of the mantle, but possible changes in physical properties and processes (such as viscosity and mantle flow patterns) in the midmantle.

Structural complexity of simple Fe2O3 at high pressures and temperatures

Although chemically very simple, Fe2O3 is known to undergo a series of enigmatic structural, electronic and magnetic transformations at high pressures and high temperatures. So far, these transformations have neither been correctly described nor understood because of the lack of structural data. Here we report a systematic investigation of the behaviour of Fe2O3 at pressures over 100 GPa and temperatures above 2,500 K employing single crystal X-ray diffraction and synchrotron Mössbauer source spectroscopy. Crystal chemical analysis of structures presented here and known Fe(II, III) oxides shows their fundamental relationships and that they can be described by the homologous series nFeO·mFe2O3. Decomposition of Fe2O3 and Fe3O4 observed at pressures above 60 GPa and temperatures of 2,000 K leads to crystallization of unusual Fe5O7 and Fe25O32 phases with release of oxygen. Our findings suggest that mixed-valence iron oxides may play a significant role in oxygen cycling between earth reservoirs.

Compositional mantle layering revealed by slab stagnation at ~1000-km depth

Improved constraints on lower-mantle composition are fundamental to understand the accretion, differentiation, and thermochemical evolution of our planet. Cosmochemical arguments indicate that lower-mantle rocks may be enriched in Si relative to upper-mantle pyrolite, whereas seismic tomography images suggest whole-mantle convection and hence appear to imply efficient mantle mixing. This study reconciles cosmochemical and geophysical constraints using the stagnation of some slab segments at ~1000-km depth as the key observation. Through numerical modeling of subduction, we show that lower-mantle enrichment in intrinsically dense basaltic lithologies can render slabs neutrally buoyant in the uppermost lower mantle. Slab stagnation (at depths of ~660 and ~1000 km) and unimpeded slab sinking to great depths can coexist if the basalt fraction is ~8% higher in the lower mantle than in the upper mantle, equivalent to a lower-mantle Mg/Si of ~1.18. Global-scale geodynamic models demonstrate that such a moderate compositional gradient across the mantle can persist can in the presence of whole-mantle convection.

Multiple seismic reflectors in Earth's lowermost mantle

The modern view of Earth's lowermost mantle considers a D″ region of enhanced (seismologically inferred) heterogeneity bounded by the core-mantle boundary and an interface some 150-300 km above it, with the latter often attributed to the postperovskite phase transition (in MgSiO3). Seismic exploration of Earth's deep interior suggests, however, that this view needs modification. So-called ScS and SKKS waves, which probe the lowermost mantle from above and below, respectively, reveal multiple reflectors beneath Central America and East Asia, two areas known for subduction of oceanic plates deep into Earth's mantle. This observation is inconsistent with expectations from a thermal response of a single isochemical postperovskite transition, but some of the newly observed structures can be explained with postperovskite transitions in differentiated slab materials. Our results imply that the lowermost mantle is more complex than hitherto thought and that interfaces and compositional heterogeneity occur beyond the D″ region sensu stricto.

Deep mantle cycling of oceanic crust: evidence from diamonds and their mineral inclusions

A primary consequence of plate tectonics is that basaltic oceanic crust subducts with lithospheric slabs into the mantle. Seismological studies extend this process to the lower mantle, and geochemical observations indicate return of oceanic crust to the upper mantle in plumes. There has been no direct petrologic evidence, however, of the return of subducted oceanic crustal components from the lower mantle. We analyzed superdeep diamonds from Juina-5 kimberlite, Brazil, which host inclusions with compositions comprising the entire phase assemblage expected to crystallize from basalt under lower-mantle conditions. The inclusion mineralogies require exhumation from the lower to upper mantle. Because the diamond hosts have carbon isotope signatures consistent with surface-derived carbon, we conclude that the deep carbon cycle extends into the lower mantle.

Upside-down differentiation and generation of a 'primordial' lower mantle

Except for the first 50-100 million years or so of the Earth's history, when most of the mantle may have been subjected to melting, the differentiation of Earth's silicate mantle has been controlled by solid-state convection. As the mantle upwells and decompresses across its solidus, it partially melts. These low-density melts rise to the surface and form the continental and oceanic crusts, driving the differentiation of the silicate part of the Earth. Because many trace elements, such as heat-producing U, Th and K, as well as the noble gases, preferentially partition into melts (here referred to as incompatible elements), melt extraction concentrates these elements into the crust (or atmosphere in the case of noble gases), where nearly half of the Earth's budget of these elements now resides. In contrast, the upper mantle, as sampled by mid-ocean ridge basalts, is highly depleted in incompatible elements, suggesting a complementary relationship with the crust. Mass balance arguments require that the other half of these incompatible elements be hidden in the Earth's interior. Hypotheses abound for the origin of this hidden reservoir. The most widely held view has been that this hidden reservoir represents primordial material never processed by melting or degassing. Here, we suggest that a necessary by-product of whole-mantle convection during the Earth's first billion years is deep and hot melting, resulting in the generation of dense liquids that crystallized and sank into the lower mantle. These sunken lithologies would have 'primordial' chemical signatures despite a non-primordial origin.

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.

Preserving noble gases in a convecting mantle

High (3)He/(4)He ratios sampled at many ocean islands are usually attributed to an essentially undegassed lower-mantle reservoir with high (3)He concentrations. A large and mostly undegassed mantle reservoir is also required to balance the Earth's (40)Ar budget, because only half of the (40)Ar produced from the radioactive decay of (40)K is accounted for by the atmosphere and upper mantle. However, geophysical and geochemical observations suggest slab subduction into the lower mantle, implying that most or all of Earth's mantle should have been processed by partial melting beneath mid-ocean ridges and hotspot volcanoes. This should have left noble gases in both the upper and the lower mantle extensively outgassed, contrary to expectations from (3)He/(4)He ratios and the Earth's (40)Ar budget. Here we suggest a simple solution: recycling and mixing of noble-gas-depleted slabs dilutes the concentrations of noble gases in the mantle, thereby decreasing the rate of mantle degassing and leaving significant amounts of noble gases in the processed mantle. As a result, even when the mass flux across the 660-km seismic discontinuity is equivalent to approximately one lower-mantle mass over the Earth's history, high (3)He contents, high (3)He/(4)He ratios and (40)Ar concentrations high enough to satisfy the (40)Ar mass balance of the Earth can be preserved in the lower mantle. The differences in (3)He/(4)He ratios between mid-ocean-ridge basalts and ocean island basalts, as well as high concentrations of (3)He and (40)Ar in the mantle source of ocean island basalts, can be explained within the framework of different processing rates for the upper and the lower mantle. Hence, to preserve primitive noble gas signatures, we find no need for hidden reservoirs or convective isolation of the lower mantle for any length of time.

Electronic and magnetic structures of the postperovskite-type Fe2O3 and implications for planetary magnetic records and deep interiors

Recent studies have shown that high pressure (P) induces the metallization of the Fe(2+)-O bonding, the destruction of magnetic ordering in Fe, and the high-spin (HS) to low-spin (LS) transition of Fe in silicate and oxide phases at the deep planetary interiors. Hematite (Fe(2)O(3)) is an important magnetic carrier mineral for deciphering planetary magnetism and a proxy for Fe in the planetary interiors. Here, we present synchrotron Mössbauer spectroscopy and X-ray diffraction combined with ab initio calculations for Fe(2)O(3) revealing the destruction of magnetic ordering at the hematite --> Rh(2)O(3)-II type (RhII) transition at 70 GPa and 300 K, and then the revival of magnetic ordering at the RhII --> postperovskite (PPv) transition after laser heating at 73 GPa. At the latter transition, at least half of Fe(3+) ions transform from LS to HS and Fe(2)O(3) changes from a semiconductor to a metal. This result demonstrates that some magnetic carrier minerals may experience a complex sequence of magnetic ordering changes during impact rather than a monotonic demagnetization. Also local Fe enrichment at Earth's core-mantle boundary will lead to changes in the electronic structure and spin state of Fe in silicate PPv. If the ultra-low-velocity zones are composed of Fe-enriched silicate PPv and/or the basaltic materials are accumulated at the lowermost mantle, high electrical conductivity of these regions will play an important role for the electromagnetic coupling between the mantle and the core.

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.

Some recent advances in understanding the mineralogy of Earth's deep mantle

Understanding planetary structure and evolution requires a detailed knowledge of the properties of geological materials under the conditions of deep planetary interiors. Experiments under the extreme pressure-temperature conditions of the deep mantle are challenging, and many fundamental properties remain poorly constrained or are inferred only through uncertain extrapolations from lower pressure-temperature states. Nevertheless, the last several years have witnessed a number of new developments in this area, and a broad overview of the current understanding of the Earth's lower mantle is presented here. Some recent experimental and theoretical advances related to the lowermost mantle are highlighted. Measurements of the equation of state and deformation behaviour of (Mg,Fe)SiO3 in the CaIrO3-type (post-perovskite) structure yield insights into the nature of the core-mantle boundary region. Theoretical studies of the behaviour of MgSiO3 liquids under high pressure-temperature conditions provide constraints on melt volumes, diffusivities and viscosities that are relevant to understanding both the early Earth (e.g. deep magma oceans) and seismic structure observed in the present Earth (e.g. ultra-low-velocity zones).

Structure and dynamics of Earth's lower mantle

Processes within the lowest several hundred kilometers of Earth's rocky mantle play a critical role in the evolution of the planet. Understanding Earth's lower mantle requires putting recent seismic and mineral physics discoveries into a self-consistent, geodynamically feasible context. Two nearly antipodal large low-shear-velocity provinces in the deep mantle likely represent chemically distinct and denser material. High-resolution seismological studies have revealed laterally varying seismic velocity discontinuities in the deepest few hundred kilometers, consistent with a phase transition from perovskite to post-perovskite. In the deepest tens of kilometers of the mantle, isolated pockets of ultralow seismic velocities may denote Earth's deepest magma chamber.