Thermochemical structures beneath Africa and the Pacific Ocean

Mesozoic intraoceanic subduction shaped the lower mantle beneath the East Pacific Rise

The Pacific large low-shear-velocity province (LLSVP), as revealed by cluster analysis of global tomographic models, hosts multiple internal anomalies, including a notable gap (~20° wide) between the central and eastern Pacific. The cause of the structural gap remains unconstrained. Directly above this structural gap, we identify an anomalously thick mantle transition zone east of the East Pacific Rise, the fastest-spreading ocean ridge in the world, using a dense set of SS precursors. The area of the thickened transition zone exhibits faster-than-average velocities according to recent tomographic images, suggesting perturbed postolivine phase boundaries shifting in response to lowered temperatures. We attribute this observation to episodes of Mesozoic-aged (250 to 120 million years ago) intraoceanic subduction beneath the present-day Nazca Plate. The eastern portion of the Pacific LLSVP was separated by downwelling because of this ancient oceanic slab. Our discovery provides a unique perspective on linking deep Earth structures with surface subduction.

Moon-forming impactor as a source of Earth's basal mantle anomalies

Seismic images of Earth's interior have revealed two continent-sized anomalies with low seismic velocities, known as the large low-velocity provinces (LLVPs), in the lowermost mantle1. The LLVPs are often interpreted as intrinsically dense heterogeneities that are compositionally distinct from the surrounding mantle2. Here we show that LLVPs may represent buried relics of Theia mantle material (TMM) that was preserved in proto-Earth's mantle after the Moon-forming giant impact3. Our canonical giant-impact simulations show that a fraction of Theia's mantle could have been delivered to proto-Earth's solid lower mantle. We find that TMM is intrinsically 2.0-3.5% denser than proto-Earth's mantle based on models of Theia's mantle and the observed higher FeO content of the Moon. Our mantle convection models show that dense TMM blobs with a size of tens of kilometres after the impact can later sink and accumulate into LLVP-like thermochemical piles atop Earth's core and survive to the present day. The LLVPs may, thus, be a natural consequence of the Moon-forming giant impact. Because giant impacts are common at the end stages of planet accretion, similar mantle heterogeneities caused by impacts may also exist in the interiors of other planetary bodies.

Mapping global kimberlite potential from reconstructions of mantle flow over the past billion years

Kimberlites are the primary source of economic grade diamonds. Their geologically rapid eruptions preferentially occur near or through thick and ancient continental lithosphere. Studies combining tomographic models with tectonic reconstructions and kimberlite emplacement ages and locations have revealed spatial correlations between large low shear velocity provinces in the lowermost mantle and reconstructed global kimberlite eruption locations over the last 320 Myr. These spatial correlations assume that the lowermost mantle structure has not changed over time, which is at odds with mantle flow models that show basal thermochemical structures to be mobile features shaped by cold sinking oceanic lithosphere. Here we investigate the match to the global kimberlite record of stationary seismically slow basal mantle structures (as imaged through tomographic modelling) and mobile hot basal structures (as predicted by reconstructions of mantle flow over the past billion years). We refer to these structures as “basal mantle structures” and consider their intersection with reconstructed thick or ancient lithosphere to represent areas with a high potential for past eruptions of kimberlites, and therefore areas of potential interest for diamond exploration. We use the distance between reconstructed kimberlite eruption locations and kimberlite potential maps as an indicator of model success, and we find that mobile lowermost mantle structures are as close to reconstructed kimberlites as stationary ones. Additionally, we find that mobile lowermost mantle structures better fit major kimberlitic events, such as the South African kimberlite bloom around 100 Ma. Mobile basal structures are therefore consistent with both solid Earth dynamics and with the kimberlite record.

Assembly of the basal mantle structure beneath Africa

Plate tectonics shapes Earth's surface, and is linked to motions within its deep interior^1,2. Cold oceanic lithosphere sinks into the mantle, and hot mantle plumes rise from the deep Earth, leading to volcanism^3,4. Volcanic eruptions over the past 320 million years have been linked to two large structures at the base of the mantle presently under Africa and the Pacific Ocean^5,6. This has led to the hypothesis that these basal mantle structures have been stationary over geological time^7,8, in contrast to observations and models suggesting that tectonic plates^9,10, subduction zones^11-14 and mantle plumes^15,16 have been mobile, and that basal mantle structures are presently deforming^17,18. Here we reconstruct mantle flow from one billion years ago to the present day to show that the history of volcanism is statistically as consistent with mobile basal mantle structures as with fixed ones. In our reconstructions, cold lithosphere sank deep into the African hemisphere between 740 and 500 million years ago, and from 400 million years ago the structure beneath Africa progressively assembled, pushed by peri-Gondwana slabs, to become a coherent structure as recently as 60 million years ago. Our mantle flow models suggest that basal mantle structures are mobile, and aggregate and disperse over time, similarly to continents at Earth's surface⁹. Our models also predict the presence of continental material in the mantle beneath Africa, consistent with geochemical data^19,20.

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.

Formation of large low shear velocity provinces through the decomposition of oxidized mantle

Large Low Shear Velocity Provinces (LLSVPs) in the lowermost mantle are key to understanding the chemical composition and thermal structure of the deep Earth, but their origins have long been debated. Bridgmanite, the most abundant lower-mantle mineral, can incorporate extensive amounts of iron (Fe) with effects on various geophysical properties. Here our high-pressure experiments and ab initio calculations reveal that a ferric-iron-rich bridgmanite coexists with an Fe-poor bridgmanite in the 90 mol% MgSiO₃-10 mol% Fe₂O₃ system, rather than forming a homogeneous single phase. The Fe³⁺-rich bridgmanite has substantially lower velocities and a higher V_P/V_S ratio than MgSiO₃ bridgmanite under lowermost-mantle conditions. Our modeling shows that the enrichment of Fe³⁺-rich bridgmanite in a pyrolitic composition can explain the observed features of the LLSVPs. The presence of Fe³⁺-rich materials within LLSVPs may have profound effects on the deep reservoirs of redox-sensitive elements and their isotopes.

A thin mantle transition zone beneath the equatorial Mid-Atlantic Ridge

The location and degree of material transfer between the upper and lower mantle are key to the Earth's thermal and chemical evolution. Sinking slabs and rising plumes are generally accepted as locations of transfer^1,2, whereas mid-ocean ridges are not typically assumed to have a role³. However, tight constraints from in situ measurements at ridges have proved to be challenging. Here we use receiver functions that reveal the conversion of primary to secondary seismic waves to image the discontinuities that bound the mantle transition zone, using ocean bottom seismic data from the equatorial Mid-Atlantic Ridge. Our images show that the seismic discontinuity at depths of about 660 kilometres is broadly uplifted by 10 ± 4 kilometres over a swath about 600 kilometres wide and that the 410-kilometre discontinuity is depressed by 5 ± 4 kilometres. This thinning of the mantle transition zone is coincident with slow shear-wave velocities in the mantle, from global seismic tomography^4-7. In addition, seismic velocities in the mantle transition zone beneath the Mid-Atlantic Ridge are on average slower than those beneath older Atlantic Ocean seafloor. The observations imply material transfer from the lower to the upper mantle-either continuous or punctuated-that is linked to the Mid-Atlantic Ridge. Given the length and longevity of the mid-ocean ridge system, this implies that whole-mantle convection may be more prevalent than previously thought, with ridge upwellings having a role in counterbalancing slab downwellings.

Velocity and density characteristics of subducted oceanic crust and the origin of lower-mantle heterogeneities

Seismic heterogeneities detected in the lower mantle were proposed to be related to subducted oceanic crust. However, the velocity and density of subducted oceanic crust at lower-mantle conditions remain unknown. Here, we report ab initio results for the elastic properties of calcium ferrite-type phases and determine the velocities and density of oceanic crust along different mantle geotherms. We find that the subducted oceanic crust shows a large negative shear velocity anomaly at the phase boundary between stishovite and CaCl₂-type silica, which is highly consistent with the feature of mid-mantle scatterers. After this phase transition in silica, subducted oceanic crust will be visible as high-velocity heterogeneities as imaged by seismic tomography. This study suggests that the presence of subducted oceanic crust could provide good explanations for some lower-mantle seismic heterogeneities with different length scales except large low shear velocity provinces (LLSVPs).

Magma oceans as a critical stage in the tectonic development of rocky planets

Magma oceans are a common result of the high degree of heating that occurs during planet formation. It is thought that almost all of the large rocky bodies in the Solar System went through at least one magma ocean phase. In this paper, we review some of the ways in which magma ocean models for the Earth, Moon and Mars match present-day observations of mantle reservoirs, internal structure and primordial crusts, and then we present new calculations for the oxidation state of the mantle produced during the magma ocean phase. The crystallization of magma oceans probably leads to a massive mantle overturn that may set up a stably stratified mantle. This may lead to significant delays or total prevention of plate tectonics on some planets. We review recent models that may help alleviate the mantle stability issue and lead to earlier onset of plate tectonics.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

Tidal tomography constrains Earth's deep-mantle buoyancy

Earth's body tide-also known as the solid Earth tide, the displacement of the solid Earth's surface caused by gravitational forces from the Moon and the Sun-is sensitive to the density of the two Large Low Shear Velocity Provinces (LLSVPs) beneath Africa and the Pacific. These massive regions extend approximately 1,000 kilometres upward from the base of the mantle and their buoyancy remains actively debated within the geophysical community. Here we use tidal tomography to constrain Earth's deep-mantle buoyancy derived from Global Positioning System (GPS)-based measurements of semi-diurnal body tide deformation. Using a probabilistic approach, we show that across the bottom two-thirds of the two LLSVPs the mean density is about 0.5 per cent higher than the average mantle density across this depth range (that is, its mean buoyancy is minus 0.5 per cent), although this anomaly may be concentrated towards the very base of the mantle. We conclude that the buoyancy of these structures is dominated by the enrichment of high-density chemical components, probably related to subducted oceanic plates or primordial material associated with Earth's formation. Because the dynamics of the mantle is driven by density variations, our result has important dynamical implications for the stability of the LLSVPs and the long-term evolution of the Earth system.

Compositionally-distinct ultra-low velocity zones on Earth's core-mantle boundary

The Earth’s lowermost mantle large low velocity provinces are accompanied by small-scale ultralow velocity zones in localized regions on the core-mantle boundary. Large low velocity provinces are hypothesized to be caused by large-scale compositional heterogeneity (i.e., thermochemical piles). The origin of ultralow velocity zones, however, remains elusive. Here we perform three-dimensional geodynamical calculations to show that the current locations and shapes of ultralow velocity zones are related to their cause. We find that the hottest lowermost mantle regions are commonly located well within the interiors of thermochemical piles. In contrast, accumulations of ultradense compositionally distinct material occur as discontinuous patches along the margins of thermochemical piles and have asymmetrical cross-sectional shape. Furthermore, the lateral morphology of these patches provides insight into mantle flow directions and long-term stability. The global distribution and large variations of morphology of ultralow velocity zones validate a compositionally distinct origin for most ultralow velocity zones.

Ultralow velocity zones are detected on the core-mantle boundary, but their origin is enigmatic. Here, the authors find that the global distribution and large variations of morphology of ultralow velocity zones are consistent with most having a compositionally-distinct origin.

Density structure of Earth's lowermost mantle from Stoneley mode splitting observations

Advances in our understanding of Earth's thermal evolution and the style of mantle convection rely on robust seismological constraints on lateral variations of density. The large-low-shear-wave velocity provinces (LLSVPs) atop the core-mantle boundary beneath Africa and the Pacific are the largest structures in the lower mantle, and hence severely affect the convective flow. Here, we show that anomalous splitting of Stoneley modes, a unique class of free oscillations that are perturbed primarily by velocity and density variations at the core-mantle boundary, is explained best when the overall density of the LLSVPs is lower than the surrounding mantle. The resolved density variations can be explained by the presence of post-perovskite, chemical heterogeneity or a combination of the two. Although we cannot rule out the presence of a ∼100-km-thick denser-than-average basal structure, our results support the hypothesis that LLSVPs signify large-scale mantle upwelling in two antipodal regions of the mantle.

Origin and evolution of the deep thermochemical structure beneath Eurasia

A unique structure in the Earth's lowermost mantle, the Perm Anomaly, was recently identified beneath Eurasia. It seismologically resembles the large low-shear velocity provinces (LLSVPs) under Africa and the Pacific, but is much smaller. This challenges the current understanding of the evolution of the plate–mantle system in which plumes rise from the edges of the two LLSVPs, spatially fixed in time. New models of mantle flow over the last 230 million years reproduce the present-day structure of the lower mantle, and show a Perm-like anomaly. The anomaly formed in isolation within a closed subduction network ∼22,000 km in circumference prior to 150 million years ago before migrating ∼1,500 km westward at an average rate of 1 cm year⁻¹, indicating a greater mobility of deep mantle structures than previously recognized. We hypothesize that the mobile Perm Anomaly could be linked to the Emeishan volcanics, in contrast to the previously proposed Siberian Traps.

The Perm anomaly is found in the lower mantle beneath Eurasia, but how this structure formed has remained unclear. Here, the authors show that the anomaly has been mobile since it formed in isolation within a closed subduction network and propose that the anomaly is linked to the Emeishan volcanics.

Kinematics and dynamics of the East Pacific Rise linked to a stable, deep-mantle upwelling

Earth's tectonic plates are generally considered to be driven largely by negative buoyancy associated with subduction of oceanic lithosphere. In this context, mid-ocean ridges (MORs) are passive plate boundaries whose divergence accommodates flow driven by subduction of oceanic slabs at trenches. We show that over the past 80 million years (My), the East Pacific Rise (EPR), Earth's dominant MOR, has been characterized by limited ridge-perpendicular migration and persistent, asymmetric ridge accretion that are anomalous relative to other MORs. We reconstruct the subduction-related buoyancy fluxes of plates on either side of the EPR. The general expectation is that greater slab pull should correlate with faster plate motion and faster spreading at the EPR. Moreover, asymmetry in slab pull on either side of the EPR should correlate with either ridge migration or enhanced plate velocity in the direction of greater slab pull. Based on our analysis, none of the expected correlations are evident. This implies that other forces significantly contribute to EPR behavior. We explain these observations using mantle flow calculations based on globally integrated buoyancy distributions that require core-mantle boundary heat flux of up to 20 TW. The time-dependent mantle flow predictions yield a long-lived deep-seated upwelling that has its highest radial velocity under the EPR and is inferred to control its observed kinematics. The mantle-wide upwelling beneath the EPR drives horizontal components of asthenospheric flows beneath the plates that are similarly asymmetric but faster than the overlying surface plates, thereby contributing to plate motions through viscous tractions in the Pacific region.

A rapid burst in hotspot motion through the interaction of tectonics and deep mantle flow

Volcanic hotspot tracks featuring linear progressions in the age of volcanism are typical surface expressions of plate tectonic movement on top of narrow plumes of hot material within Earth's mantle. Seismic imaging reveals that these plumes can be of deep origin--probably rooted on thermochemical structures in the lower mantle. Although palaeomagnetic and radiometric age data suggest that mantle flow can advect plume conduits laterally, the flow dynamics underlying the formation of the sharp bend occurring only in the Hawaiian-Emperor hotspot track in the Pacific Ocean remains enigmatic. Here we present palaeogeographically constrained numerical models of thermochemical convection and demonstrate that flow in the deep lower mantle under the north Pacific was anomalously vigorous between 100 million years ago and 50 million years ago as a consequence of long-lasting subduction systems, unlike those in the south Pacific. These models show a sharp bend in the Hawaiian-Emperor hotspot track arising from the interplay of plume tilt and the lateral advection of plume sources. The different trajectories of the Hawaiian and Louisville hotspot tracks arise from asymmetric deformation of thermochemical structures under the Pacific between 100 million years ago and 50 million years ago. This asymmetric deformation waned just before the Hawaiian-Emperor bend developed, owing to flow in the deepest lower mantle associated with slab descent in the north and south Pacific.

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.

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.

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.

A Post-Perovskite Lens and D'' Heat Flux Beneath the Central Pacific

Temperature gradients in a low-shear-velocity province in the lowermost mantle (D'' region) beneath the central Pacific Ocean were inferred from the observation of a rapid S-wave velocity increase overlying a rapid decrease. These paired seismic discontinuities are attributed to a phase change from perovskite to post-perovskite and then back to perovskite as the temperature increases with depth. Iron enrichment could explain the occurrence of post-perovskite several hundred kilometers above the core-mantle boundary in this warm, chemically distinct province. The double phase-boundary crossing directly constrains the lowermost mantle temperature gradients. Assuming a standard but unconstrained choice of thermal conductivity, the regional core-mantle boundary heat flux (approximately 85 +/- 25 milliwatts per square meter), comparable to the average at Earth's surface, was estimated, along with a lower bound on global core-mantle boundary heat flow in the range of 13 +/- 4 terawatts. Mapped velocity-contrast variations indicate that the lens of post-perovskite minerals thins and vanishes over 1000 kilometers laterally toward the margin of the chemical distinct region as a result of a approximately 500-kelvin temperature increase.