Superplume Research Articles

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

AbstractThe 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.

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.

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity

The compositional and thermal state of Earth's mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earth's early history.

How the Indian Ocean Geoid Low Was Formed

AbstractThe origin of the Earth's lowest geoid, the Indian Ocean geoid low (IOGL) has been controversial. The geoid predicted from present‐day tomography models has shown that mid to upper mantle hot anomalies are integral in generating the IOGL. Here we assimilate plate reconstruction in global mantle convection models starting from 140 Ma and show that sinking Tethyan slabs perturbed the African Large Low Shear Velocity province and generated plumes beneath the Indian Ocean, which led to the formation of this negative geoid anomaly. We also show that this low can be reproduced by surrounding mantle density anomalies, without having them present directly beneath the geoid low. We tune the density and viscosity of thermochemical piles at core‐mantle boundary, Clapeyron slope and density jump at 660 km discontinuity, and the strength of slabs, to control the rise of plumes, which in turn determine the shape and amplitude of the geoid low.

Tomographic filtering of shear and compressional wave models reveals uncorrelated variations in the lowermost mantle

SUMMARY Seismic tomography models reveal differences in the geographic distribution and magnitude of P- and S-wave velocity variations (VP and VS, respectively) below ∼2200 km depth in the Earth’s mantle. In particular, large low shear velocity provinces (LLSVPs) beneath the Pacific and Africa exhibit a distinct low velocity population in the distribution of VS that does not stand out in VP models, carrying important implications for the origin of these features. However, it is possible that the absence of a distinct low velocity feature in VP models is an artefact of VP models having lower resolution compared to VS models owing to differences in coverage. Here, we use ‘tomographic filters’ computed from the singular value decomposition of the sensitivity matrices for a pair of VP and VS models in order to test whether such low velocity features are suppressed in VP models. Our ‘cross-filtered’ results show that resolution alone cannot explain the absence of a corresponding low VP population. We additionally apply the joint VP and VS tomographic filter technique to thermochemical mantle convection models to show that cases with distinct phase and/or composition may be differentiated from cases where only temperature varies. We then develop a new proxy for exploring uncorrelated VP and VS more broadly using the difference between the observed VP model and the filtered VS model input. Our results show that ‘large uncorrelated modulus provinces’ (LUMPs) extend beyond the boundaries of LLSVPs, and exhibit anomalies in both fast and slow regions.

A peridotite source for strongly alkalic ultrabasic HIMU lavas of the Oslo Rift, Norway

High alkali (Na2O + K2O) low SiO2 lavas generally have highly fractionated trace element compositions and their petrogenesis is debated. Both pyroxene-rich (pyroxenite) and olivine-rich (peridotite) mantle sources have been proposed for such lavas. HIMU (high 238U/204Pb) lavas usually have such compositions and are often interpreted as being derived from garnet-pyroxenite (eclogite) residue of subducted recycled oceanic crust, possibly including a carbonate component. In some cases, an origin from a garnet peridotite source has also been considered.We present new Mg, Pb, Nd, and Sr isotopic data and elemental data for alkalic ultrabasic lavas (melilitites, nephelinites) from the Oslo Rift, Norway. The magmatic province of the Permo-Carboniferous Oslo Rift is interpreted as part of the much larger Skagerrak-Centered Large Igneous Province (SCLIP) which may be related to a mantle plume source originating from the African Large Low Shear Velocity Province (LLSVP) at the core-mantle boundary. The range of δ26Mg is narrow (−0.32 to −0.28) and is similar to typical mantle values, and therefore is inconsistent with a carbonated source as has been suggested for other similar lavas. The 206Pb/204Pb values of all the lavas are high, which identifies them as HIMU lavas. The εNd of +1 to +2 are substantially lower than for typical HIMU lavas, and in the εNd-εSr diagram they plot on the LoNd array between HIMU and EMI. The lavas are more like archetypal kimberlites (formerly, Group I kimberlites), as these also have εNd slightly higher than the chondritic value and Pb isotopes consistent with a HIMU source. The REE patterns of the lavas are also almost as fractionated as in Group I kimberlites.We used major and trace elements to infer the primary magma composition of these lavas. This primary magma composition is compared with experimental melts of both carbonated peridotite and carbonated eclogite source compositions and shows that their major element compositions are only consistent with either a very high degree of melting of eclogite or a very low degree melting of garnet peridotite. The strongly fractionated REE patterns of the lavas can only be made by low-degree melting of either eclogite or garnet peridotite, which therefore effectively rules out eclogite as a major component of the mantle source of these lavas. Further trace element modeling points to a mantle source composition with >97% garnet peridotite that was lightly carbonated (<0.5 wt% CO2) and melted to a very low degree (∼0.4%) to produce these lavas. The modeled source trace element pattern was constrained by the Nd, Sr, and Pb isotopic compositions of the lavas. The measured high Ti contents of these lavas (3.8–6 wt%) are quantitatively shown to be the result of enrichment by both partial melting and fractional crystallization, using experimentally determined DTi for melting of this type of source. We also show that a low degree of melting of a slightly carbonated source is sufficient to result in the likely very high CO2 content (10–15 wt%) of the primary magma.The rather uniform Pb isotope composition of HIMU magmas suggests that the HIMU source is rather uniform in composition, but the Pb isotope composition requires this source to form late in Earth's history. The similarity of these lavas and kimberlites suggests a very deep source, possibly in the LLSVP that gave origin to the SCLIP. This means that LLSVPs do not represent some primordial layer around the core that escaped homogenization of the mantle by the Moon-forming giant impact. Instead, this layer must have formed late in Earth's history. The constraint that this source contains <3% of garnet-pyroxenite (eclogite) residue of recycled basaltic oceanic crust shows that this layer is primarily not a subducted slab graveyard. The process by which the LLSVPs formed, however, remains obscure.

Antipodal waveform observations of seismic waves diffracting and refracting at the base of Earth's outer core

We approach this analysis of the inner-outer core boundary (IOCB) with antipodal waveform data in the distance range 179.0°–180° to test the hypothesis whether propagation at the base of the outer core is commensurate with diffraction and/or refraction. The propagation paths observed cover about two-thirds of the IOCB surface. Seismic data from seven diameters are examined—Tonga to Algeria (station code TAM), Sulawesi to Amazon (PTGA), northern Chile to Hainan Island (QIZ), two between central Chile and the mainland China (XAN, ENH), and two diameters from New Zealand to both Portugal (PTO) and Spain (ECAL). The adequacy of global Earth models—both 1-D and 1-D core +3D mantle—in fitting the antipodal observations is found to be deficient. Lateral heterogeneity among the propagation paths are mapped and projected to the Earth's surface for context. We stack data to increase signal-to-noise, and model waveform data via the 3D spectral element method on the EarthSimulator4. With the exception of one path to ENH China, none of the global Earth models match amplitudes of diffracted waves or the stacked amplitudes relative to PKIKP, which travels along the antipodal diameter. Energetic arrivals observed in the waveform data set are modeled as a combination of refraction within a low velocity zone (LVZ) at the base of the outer core, and diffraction around this structure. The observations at the base of the outer core may be further subdivided between paths with a thin low velocity zone (QIZ, XAN, PTGA), and those with a thicker, slower LVZ (TAM). The former paths are characterized by a thin LVZ (20–50 km thick, layer–gradient, respectively) with a velocity of approximately 10.0 km/s. The second group of paths shows large amplitude, complicated waveforms for TAM, corroborated by data from PTO and ECAL. TAM waveforms are modeled with a negative gradient from 100 to 50–75 km thickness above the IOCB, to a basal velocity of 8.8 km/s at IOCB. We also modeled the prospect of non-zero rigidity at the bottom of the outer core, but found no evidence of sensitivity in the core data. The structures fitting the data correlate across the inner–outer core boundary, with a velocity anomaly in PKPbc data juxtaposed immediately above the TAM path beneath SE Asia. More generally, the TAM IOCB path when projected up to the core-mantle boundary (CMB) apparently aligns with regions between the large low shear velocity provinces (LLSVP) at the CMB. Other paths with thinner or negligible LVZ at the IOCB partially underlie the LLSVP when projected to the CMB.