Superplume Research Articles

Strong precursor softening in cubic CaSiO3 perovskite

CaSiO[Formula: see text] perovskite (CaPv) is the last major mineral in the Earth's lower mantle whose elasticity remains largely unresolved. Here, we investigate the elasticity of CaPv using ab initio machine-learning force fields (MLFF). At room temperature, the elasticity of tetragonal CaPv determined by MLFF molecular dynamics (MD) agrees well with experimental measurements after considering temperature induced variations in the hydrostatic structure, proving the effectiveness of the method. We use the MLFF MD in the [Formula: see text] ensemble to establish the tetragonal-cubic phase boundary and confirm that in the lower mantle CaPv is in the cubic phase. The elasticity of cubic CaPv shows distinct temperature dependence at different ranges: it is linear at high temperatures, whereas it exhibits anomalous precursor softening near the tetragonal-cubic phase boundary. The temperature interval of precursor softening widens as the pressure increases and overlaps with the temperature profile of subducted cold slabs near the core-mantle boundary. While cubic CaPv is seismically invisible along the average mantle geotherm, it may induce low-velocity zones with negative temperature anomaly, leading to the view that the large low shear velocity provinces (LLSVPs) may be caused by subducted oceanic crust rich in CaPv with temperature lower than ambient mantle. A cool, rigid LLSVP may help explain the preferential formation of mantle plumes at its margins, as well as its weaker seismic anisotropy.

Spherically symmetric Earth models yield no net electron spin

Terrestrial experiments that use electrons in Earth as a spin-polarized source have been demonstrated to provide strong bounds on exotic long-range spin-spin and spin-velocity interactions. These bounds constrain the coupling strength of many proposed ultralight bosonic dark-matter candidates. Recently, it was pointed out that a monopole-dipole coupling between the Sun and the spin-polarized electrons of Earth would result in a modification of the precession of the perihelion of Earth. Using an estimate for the net spin polarization of Earth and experimental bounds on Earth’s perihelion precession, interesting constraints were placed on the magnitude of this monopole-dipole coupling. Here we investigate the spin associated with Earth’s electrons. We find that there are about 6×1041 spin-polarized electrons in the mantle and crust of Earth oriented antiparallel to their local magnetic field. However, when integrated over any spherically symmetric Earth model, we find that the vector sum of these spins is zero. In order to establish a lower bound on the magnitude of the net spin along Earth’s rotation axis we have investigated three of the largest breakdowns of Earth’s spherical symmetry: the large low shear-velocity provinces of the mantle, the crustal composition, and the oblate spheroid of Earth. From these investigations we conclude that there are at least 5×1038 spin-polarized electrons aligned antiparallel to Earth’s rotation axis. This analysis suggests that the bounds on the monopole-dipole coupling that were extracted from Earth’s perihelion precession need to be relaxed by a factor of about 2000. Published by the American Physical Society 2025

The Evolution of Dense ULVZs Originating Outside LLSVPs and Implications for Dynamics at LLSVP Margins

AbstractInteractions between multiple‐scale thermochemical heterogeneities in the lowermost mantle, specifically ultralow velocity zones (ULVZs) and large low shear velocity provinces (LLSVPs), are critical in lower mantle dynamics. However, the evolution of ULVZs formed outside LLSVPs has not been thoroughly explored. Here we perform two‐dimensional numerical experiments to examine the evolution of highly dense ULVZs originating beneath cold downwellings and their interactions with the LLSVP. We find that ULVZs with an intrinsic density anomaly more than 500 kg/m3 compared with the ambient lowermost mantle cannot fully enter the LLSVP and would dwell at LLSVP margins for an indefinitely long time. This suggests that dense ULVZ within LLSVPs might have different sources from those outside LLSVPs. The buoyancy number and compositional viscosity of ULVZs are controlling factors on their dynamics and imprints on the core‐mantle boundary (CMB), such as how much the ULVZ protrudes into the LLSVP and the CMB topography beneath the ULVZ. The excess density of ULVZs dictates their width but not their thickness. The oscillations of ULVZ morphology suggest that various types of plumes occur at the LLSVP margin. The mobility of ULVZ implies that the bottom margin of the LLSVP moves much faster than its center. Hot zones exist within the LLSVP near its margins, which may affect the evolution of ULVZs and subducted material nearby. The CMB topography under dense ULVZs are positive unless the buoyancy number of ULVZs exceeds 6.0. These results have intriguing implications for the distribution of ULVZs as well as thermochemical evolution in the lowermost mantle.

The Iron Spin Transitions in Hydrous Fe3+‐Bearing Bridgmanite and Its Geophysical Properties in the Lower Mantle

AbstractHydrous Fe3+‐bearing bridgmanite (Bdg) is potentially a critical water host in the lowermost mantle. The spin transition behaviors of such materials are pivotal for understanding geophysical heterogeneity in the deep Earth but are poorly understood. Here, we investigated the spin transition and related geophysical properties of Fe3+ with associated H defects [Fe3+‐H] at high P‐T conditions using first‐principles simulations. Our calculations predict that the presence of hydrogen reduces the onset pressure of the spin transition of Fe3+ in Bdg, leading to higher fractions of low spin Fe3+. Along standard geotherms, spin transition is predicted to remain incomplete even at the core‐mantle boundary (CMB), and lateral temperature variations would significantly affect the proportions of high and low spin Fe and related properties like elasticity. The thermoelastic property of hydrous Fe3+‐bearing bridgmanite exhibit stronger softening anomalies at the lower mantle conditions compared to dry system, which potentially enhancing the seismic detectability of the hydrous Bdg in the deep earth. Density profile of hydrous Fe3+‐bearing bridgmanite indicates that the [Fe3+‐H] defect modestly increases the system's density, but much less than that caused by incorporating an equivalent amount of iron alone. This is crucial for understanding regions like Large Low Shear Velocity Provinces (LLSVPs), which exhibits large velocity drops but only minor density changes. The co‐adsorption of Fe and H allows for the introduction of Fe to induce velocity drops without the concomitant sharp increase in density, as pure iron would, thus enabling Fe‐H enrichment as a potential source of LLSVPs.

Synchronous motion of the Easter mantle plume and the East Pacific Rise

The Easter mantle plume has produced one of the longest hotspot tracks in the Pacific Ocean. While previous studies have focused on the eastern side extending across the Nazca Plate, we use 40Ar/39Ar isotopic and geochemical data to investigate the less explored western side around the Easter Microplate. We propose a dynamic model in which a deeper (600 km-depth), less buoyant mantle exerts a westward force on the East Pacific Rise (EPR), while a more buoyant plume region drives Easter hotspot volcanism and a localised acceleration in seafloor spreading. Our findings suggest that the Easter hotspot is the more focused surface expression of the most buoyant region of a vast, deep-seated mantle plume extending from the Pacific Large Low Shear Velocity Province (LLSVP). This challenges the traditional view of hotspots as isolated phenomena and suggests they are part of broader LLSVP-related mantle structures. Our results imply a more intricate, large-scale relationship between hotspots, mantle plumes, spreading ridges, and mantle dynamics.

Refining Heterogeneities near the Core–Mantle Boundary Beneath East Pacific Regions: Enhanced Differential Travel-Time Analysis Using USArray

Recent advancements in seismic data analysis have enhanced our grasp of the seismic heterogeneities near the core–mantle boundary (CMB). Through seismic tomography, persistent lower-mantle structures like the large low shear velocity provinces (LLSVPs) beneath the Pacific and South Africa have been identified. However, variations in the finer-scale features across different models raise questions about their origins. This study utilizes differential travel-time measurements from the USArray, operational across the contiguous United States from 2007 to 2014, to examine the impact of upper-mantle heterogeneities on tomographic models. By averaging the P-wave travel times and calibrating them with diffracted P-waves at the same stations, we mitigate the effects of shallow heterogeneities. The findings confirm that this method accurately maps the seismic anomalies beneath the USArray, consistent with other North American studies. Calibrated Pdiff travel-time data indicate significant anomalies in the mid-Pacific and Bering Sea and lesser anomalies in the northern Pacific, aligning with the global tomographic images. Moreover, the study highlights sharp travel-time variations over short distances, such as those across the northern boundary of the mid-Pacific anomaly, suggesting a chemically heterogeneous Pacific LLSVP.

Deep mantle plumes feeding periodic alignments of asthenospheric fingers beneath the central and southern Atlantic Ocean

High-resolution full waveform seismic tomography of the Earth's mantle beneath the south and central Atlantic Ocean brings into focus a series of asthenospheric low shear velocity channels, or "fingers" on both sides of the southern and central mid-Atlantic ridge (MAR), elongated in the direction of absolute plate motion with a spacing of [Formula: see text]1,800 to 2,000 km, and associated with bands of shallower residual seafloor depth anomalies that suggest channeled flow over thousands of kilometers. Each of the three most clearly resolved fingers on the African side of the MAR corresponds to a separate group of whole mantle plumes rooted in distinct patches at the core-mantle boundary, feeding hotspots, and volcanic lines with distinct isotopic signatures. Plumes of a given group appear to merge at the top of the lower mantle before separating again, suggesting interaction of deep mantle flow with a more vigorous mesoscale circulation in the upper mantle. The corresponding hotspots are generally offset from the location of the deep mantle plume roots. The distinct isotopic signatures of these hotspot groups are also detected in the mid-ocean ridge basalts at the location where the fingers meet the ridge. Meanwhile, at least some of the variability within each plume group could originate in the upper mantle and extended transition zone where plumes in a given group appear to merge and pond. This study also adds to mounting evidence that the African large low shear velocity province is not a uniform, unbroken pile of dense material rising high above the core-mantle boundary, but rather a collection of mantle plumes rooted in patches of distinct composition.

Influence of Grain Size Evolution on Mantle Plume and LLSVP Dynamics

AbstractRecent seismic tomography models suggest large‐radius primary plumes originating from the core‐mantle boundary, with grain size variations potentially explaining these observations. Additionally, grain size variations are thought to enhance the long‐term stability of Large Low Shear Velocity Provinces (LLSVPs), identified as thermochemical piles near the core‐mantle boundary. Nevertheless, geodynamic models investigating these hypotheses remain limited. To address this gap, we constructed a series of geodynamic numerical models incorporating grain size evolution, plate tectonics, and the spontaneous generation of deep mantle plumes above LLSVPs. Our results reveal that grain size evolution does not significantly affect the plume width, primarily because the increased strain rate in the mantle plume suppresses both its grain size and viscosity. The region adjacent to the plumes, characterized by the accumulation of mantle materials with larger grain size and low‐temperature remnants of subducted slabs, displays a higher viscosity compared to the area near the subducted slabs. Furthermore, grain size evolution plays a crucial role in enhancing the stability of LLSVPs by increasing the viscosity ratio between LLSVPs and the ambient mantle. These findings underscore the need for incorporating grain size evolution in geodynamic models to gain a better understanding of the dynamics of plumes and lower mantle.

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.

Ultra‐Low Velocity Zone Beneath the Atlantic Near St. Helena

AbstractThere are various hotspots in the Atlantic Ocean, which are underlain by mantle plumes that likely cross the mantle and originate at the core‐mantle boundary. We use teleseismic core‐diffracted shear waves to look for an Ultra‐Low Velocity Zone (ULVZ) at the potential base of central Atlantic mantle plumes. Our data set shows delayed postcursory phases after the core‐diffracted shear waves. The observed patterns are consistent in frequency dependence, delay time, and scatter pattern with those caused by mega‐ULVZs previously modeled elsewhere. Synthetic modeling of a cylindrical structure on the core‐mantle boundary below St. Helena provides a good fit to the data. The preferred model is 600 km across and 20 km high, centered at approximately 15° South, 15° West, and with a 30% S‐wave velocity reduction. Significant uncertainties and trade‐offs do remain to these parameters, but a large ULVZ is needed to explain the data. The location is west of St. Helena and south of Ascension. Helium and neon isotopic systematics observed in samples from this region could point to a less‐outgassed mantle component mixed in with the dominant signature of recycled material. These observations could be explained by a contribution from the Large Low Shear Velocity Province (LLSVP). Tungsten isotopic measurements would be needed to understand whether a contribution from the mega‐ULVZ is also required at St. Helena or Ascension.

Present day mantle structure from global mantle convection models since the Cretaceous

SUMMARY Using forward mantle convection models starting at 140 Ma, and assimilating plate reconstructions as surface velocity boundary condition, we predict present-day mantle structure and compare them with tomography models, using geoid as an additional constraint. We explore a wide model parameter space, such as different values of Clapeyron slope and density change across 660 km, density and viscosity of the thermochemical piles at the core–mantle boundary (CMB), internal heat generation rate, and model initiation age. We also investigate the effects of different strengths of a weak layer below 660 km and weaker asthenosphere and slabs. Our results suggest that slab structures at different subduction zones are sensitive to the viscosity of the asthenosphere, strength of slabs, values of Clapeyron slope and the density and viscosity of the thermochemical piles, while different internal heat generation rates do not affect the slab structures. We find that with a moderately weak asthenosphere ($10^{20}$ Pa·s) and strong slabs, the predicted slab structures are consistent with the tomography models, and the observed geoid is also matched well. Moreover, our models successfully reproduce the degree-2 structure of the lower mantle beneath Africa and the Pacific, also known as Large Low Shear Velocity provinces (LLSVPs). A moderate Clapeyron slope of −2.5 MPa K−1 at 660 km aids in slab stagnation while higher values result in massive slab accumulation at that depth, ultimately leading to slab avalanches. We also find that the convective patterns in the thermal and thermochemical cases with slightly denser LLSVPs are similar, although the geoid amplitudes are lower for the latter. However, with more dense LLSVPs, the slabs cannot perturb them and no plumes are generated. Plumes arise as thermal instabilities from the edges of the LLSVPs, when cold and viscous slabs perturb them. While our predicted plume locations are consistent with the observed hotspot locations, matching the plume structures in tomography models is difficult. These plumes are essential in fitting the finer features of the observed geoid. In longer-duration models, more voluminous subducted material reaches the CMB, which tends to erode the LLSVPs significantly, and yields a poor fit to the observed geoid. Our results suggest that with the presence of a thin, moderately weak layer below 660 km, a slightly dense LLSVP, and Clapeyron slope of −2.5 MPa K−1, the velocity anomalies in seismic tomography and the long-wavelength geoid can be matched well. One of the limitations of our models is that the assimilated plate motion history may be too short to overcome arbitrary initial conditions effects. Also, assimilated true plate velocities in our models may not represent the true convective vigour of the Earth.

A Giant Impact Origin for the First Subduction on Earth

AbstractHadean zircons provide a potential record of Earth's earliest subduction 4.3 billion years ago. It remains enigmatic how subduction could be initiated so soon after the presumably Moon‐forming giant impact (MGI). Earlier studies found an increase in Earth's core‐mantle boundary (CMB) temperature due to the accumulation of the impactor's core, and our recent work shows Earth's lower mantle remains largely solid, with some of the impactor's mantle potentially surviving as the large low‐shear velocity provinces (LLSVPs). Here, we show that a hot post‐impact CMB drives the initiation of strong mantle plumes that can induce subduction initiation ∼200 Myr after the MGI. 2D and 3D thermomechanical computations show that a high CMB temperature is the primary factor triggering early subduction, with enrichment of heat‐producing elements in LLSVPs as another potential factor. The models link the earliest subduction to the MGI with implications for understanding the diverse tectonic regimes of rocky planets.

Coherent Seismic Anisotropy Pattern Across Southern Africa Revealed by Shear Wave Splitting Measurements

AbstractWe report new PKS, SKS, and SKKS splitting measurements for 88 seismic stations in Namibia, Botswana, South Africa, and Mozambique. When combined with measurements from previous studies, the ensemble of measurements shows a fairly uniform NNE to NE (∼41° on average) fast‐polarization direction (ϕ) and delay time (δt) (∼0.7 s on average) across the entire southern African subcontinent. It is difficult to attribute the NNE‐NE ϕ direction to just one source of anisotropy either within the lithospheric or sublithospheric mantle. We instead propose the observed anisotropy pattern could result from a combination of several sources that together give rise to a pervasive NNE‐NE ϕ direction; (a) fossil anisotropy in the lithospheric mantle resulting from the Neoproterozoic collision of the Congo and Kalahari cratons to form the Damara Belt, (b) movement of the African plate over the asthenosphere, and (c) flow in the upper mantle induced by the African Superplume. In addition, a contribution from anisotropy in the lowermost mantle in the vicinity of the African large low shear velocity province cannot be ruled out.

Structure and equation of state of Ti-bearing davemaoite: new insights into the chemical heterogeneity in the lower mantle

Abstract Davemaoite (CaSiO3 perovskite) is considered the third most abundant phase in the pyrolytic lower mantle and the second most abundant phase in the subducted mid-ocean ridge basalt (MORB). During the partial melting of the pyrolytic upper mantle, incompatible titanium (Ti) becomes enriched in the basaltic magma, forming Ti-rich MORB. Davemaoite is considered an important Ti-bearing mineral in subducted slabs by forming a Ca(Si,Ti)O3 solid solution. However, the crystal structure and compressibility of Ca(Si,Ti)O3 perovskite solid solution at relevant pressure and temperature conditions were not systematically investigated. In this study, we investigated the structure and equations of state of Ca(Si0.83Ti0.17)O3 and Ca(Si0.75Ti0.25)O3 perovskites at room temperature up to 82 GPa and 64 GPa, respectively, by synchrotron X-ray diffraction (XRD). We found that both Ca(Si0.83Ti0.17)O3 and Ca(Si0.75Ti0.25)O3 perovskites have a tetragonal structure up to the maximum pressures investigated. Based on the observed data and compared to pure CaSiO3 davemaoite, both Ca(Si0.83Ti0.17)O3 and Ca(Si0.75Ti0.25)O3 perovskites are expected to be less dense up to the core-mantle boundary (CMB), and specifically ~1-2 % less dense than CaSiO3 davemaoite in the pressure range of the transition zone (15-25 GPa). Our results suggest that the presence of Ti-bearing davemaoite phases may result in a reduction in the average density of the subducting slabs, which in turn promotes their stagnation in the lower mantle. The presence of low-density Ti-bearing davemaoite phases and subduction of MORB in the lower mantle may also explain the seismic heterogeneity in the lower mantle, such as large low shear velocity provinces (LLSVPs).

Structure and elasticity of CaC2O5 suggests carbonate contribution to the seismic anomalies of Earth’s mantle

Knowledge of carbonate compounds under high pressure inside Earth is key to understanding the internal structure of the Earth, the deep carbon cycle and major geological events. Here we use first-principles simulations to calculate the structure and elasticity of CaC2O5-minerals with different symmetries under high pressure. Our calculations show that CaC2O5-minerals represent a group of low-density low-seismic-wave velocity mantle minerals. Changes in seismic wave velocity caused by the phase transformation of CaC2O5-Cc to CaC2O5-I4¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\bar{4}$$\\end{document}2d (CaC2O5-C2-l) agree with wave velocity discontinuity at a depth of 660 km in the mantle transition zone. Moreover, when CaC2O5-Fdd2 transforms into CaC2O5-C2 under 70 GPa, its shear wave velocity decreases by 7.4%, and its density increases by 5.8%, which is consistent with the characteristics of large low-shear-velocity provinces (LLSVPs). Furthermore, the shear wave velocity of CaC2O5-I4¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\bar{4}$$\\end{document}2d is very similar to that of cubic Ca-perovskite, which is one of the main constituents of the previously detected LLSVPs. Therefore, we propose that CaC2O5 and its high-pressure polymorphs may be a main component of LLSVPs.

Basal Mantle Flow Over LLSVPs Explains Differences in Pacific and Indo‐Atlantic Hotspot Motions

AbstractSurface hotspot motions are approximately a factor of two faster in the Pacific than the Indo‐Atlantic, and the Indo‐Atlantic large low shear velocity province (LLSVP) appears to be significantly taller than the Pacific LLSVP. Hypothesizing that surface hotspot motions are correlated with the motion of plume sources on the upper surface of chemically distinct, intrinsically dense LLSVPs, we use 3D spherical mantle convection models to compute the velocity of plume sources and compare with observed surface hotspot motions. No contrast in the mean speed of Pacific and Indo‐Atlantic hotspots is predicted if the LLSVPs are treated as purely thermal anomalies and plume sources move laterally across the core‐mantle boundary. However, when LLSVP topography is included in the model, the predicted hotspot speeds are, on average, faster in the Pacific than the Indo‐Atlantic, even when modest topography is assigned to both LLSVPs (e.g., 100–300 km). The difference in mean hotspot speed increases to a factor of two for larger and laterally variable LLSVP topography estimated from seismic tomographic model S40RTS (up to 1,100–1,500 km for the Indo‐Atlantic region vs. 700–1,400 km for the Pacific region) and our results also broadly reproduce the convergence of Pacific hotspots toward the center of the Pacific LLSVP. These largescale features of global hotspot motions are only reproduced when ambient mantle material flows over large, relatively stable topographical features, suggesting that LLSVPs are chemically distinct and intrinsically dense relative to ambient mantle material.

Perovskite-bearing crystal-controlled oxide-silicate mantle xenoliths: Resolution to controversial origins?

Classic lamellar clinopyroxene-ilmenite intergrowths (type 1) are extended to include discovery of olivine-ilmenite-perovskite-wüstite (type 2) and olivine-spinel-perovskite (type 3) xenoliths in kimberlites from Liberia. Low titanium solubilities in olivine, garnet, and pyroxene cannot account for exsolution-like relations. Because the oxides coexist with high-pressure perovskite-structured silicate minerals in diamond, a permissive conclusion is that type 1 to type 3 xenoliths are of super-deep origin. Phase equilibria and thermodynamic studies show that type 1 xenoliths are stable at P > 80 GPa, with type 2 and type 3 at 35 to 50 GPa consistent with an origin in anomalous large low shear velocity province bodies anchored at the core-mantle boundary. Dissociated precursor perovskite-structured Ca-Fe-Ti bridgmanite is proposed and is indirectly supported by the copresence of type II diamonds with a sublithospheric lower mantle origin.

Decoupling of short-lived radiogenic and helium isotopes in the Marquesas hotspot

The short-lived radiogenic 182Hf/182W (t½ = 9 Ma) and 146Sm/142Nd (t½ = 103 Ma) systems are useful tools to investigate differentiation processes in Earth's early history. Mantle reservoirs that have been preserved over long geological timescales could carry signatures formed early on in Earth's history and later could act as sources to modern mantle plumes. Negative μ182W anomalies in combination with high 3He/4He ratios have been detected globally as differently sloping negative trends in plume-derived ocean island basalts (OIB) from individual hotspots. Large low shear velocity provinces (LLSVP) and ultra-low velocity zones (ULVZ) are seismically anomalous regions in the lower mantle that have been suggested as sources of primordial mantle components near the core-mantle boundary (CMB) and accessed by deep-rooted mantle plumes.In this work, seven samples from the islands Hiva Oa and Fatu Hiva of the Marquesas Archipelago, showing a range in 3He/4He ratios, have been investigated motivated by its location above the Pacific LLSVP and a recently discovered “mega-ULVZ”. Despite elevated 3He/4He ratios of up to 14.4 R/RA, new high precision measurements using thermal ionization mass spectrometry revealed non-anomalous μ182W (−0.9 ± 4.9 ppm) and μ142Nd (1.3 ± 2.4 ppm, 95% confidence interval (c.i.)). These mark the Marquesas Archipelago as the first hotspot to show a decoupling of 3He/4He ratios and μ182W in OIB. A mixing model between an elevated 3He/4He-negative-μ182W and an EM2 or a HIMU endmember suggest that the addition of recycled material is not responsible for this decoupling. Instead, like many previously studied OIB, the Marquesas plume is most likely dominated by an LLSVP component proposed to contain high 3He/4He ratios and μ182W of 0. The unique observed absence of negative μ182W anomalies, on the other hand, implies a lack of involvement of the ULVZ source in the plume's melting region, despite the hotspot being located above a recently discovered mega-ULVZ. This may imply that the expression of ancient mantle domains in hotspot geochemistry is temporally variable, thus, delaying the arrival or involvement of ULVZ material in the origin of the Marquesas OIB.

Mapping structures on the core–mantle boundary using Sdiff postcursors: Part II. Application to the Hawaiian ULVZ

SUMMARY We present a new data set of nearly 100 earthquakes which show clear evidence of the Hawaiian ultra-low velocity zone (ULVZ) in the S core-diffracted phase (Sdiff), representing the most comprehensive Sdiff data set of a ULVZ to date. Using a Bayesian inversion approach, as outlined in Martin et al., and a subset of the data set, we characterise the 2-D morphology and velocity of the Hawaiian ULVZ. The results suggest that the ULVZ is smaller than previously estimated, with an elliptical shape, and oriented along the direction of the large low-shear velocity province boundary. Using forward modelling, we infer that the ULVZ has a height of 25 km and shear velocity reduction of 25 %.

Inferring the relationship between core-mantle heat flux and seismic tomography from mantle convection simulations

The heat flux pattern at Earth’s core-mantle boundary (CMB) imposes a heterogeneous boundary condition on core dynamics that may profoundly affect the geodynamo. Owing to the expected temperature dependence of seismic velocities, this pattern is classically approximated as proportional to the lowermost layer of seismic tomography models for the global mantle. Two biases however undermine such a simple linear relationship: 1) other contributions than thermal (compositional and mineralogical) influence seismic velocities and 2) the radial average is inherent to tomographic models whereas the local thermal state at the CMB is relevant for the heat flux. We analyze here simulations of thermochemical mantle convection where, owing to their spatial characteristics, specific mantle components are readily identified: hot thermochemical piles (TCPs), “normal” mantle (NM) and, when post-peroskite (pPv) is included, a cold region where this phase is present. Synthetic seismic velocities (i.e. from the mantle simulations) are then computed based on thermal, compositional and mineralogical sensitivities. A formalism to infer the CMB heat flux from these seismic shear velocity anomalies is derived. In this formalism, within each mantle population (i.e. TCPs, NM or pPv) the CMB heat flux vs. seismic anomalies follows a unique fitting function. The transition from one mantle population to another is marked by a jump in the seismic anomaly, i.e. a range of seismic anomalies in between two mantle populations corresponds to a similar CMB heat flux. Applying our formalism to the seismic anomalies from the mantle convection simulations provides far superior fits than the commonly used linear fits. The results highlight reduced negative heat flux anomalies beneath large low shear velocity provinces (LLSVPs), while positive heat flux anomalies are enhanced, both with respect to the classical linear interpretation.