Lowermost Mantle Research Articles

Partial Deoxygenation and Dehydration of Ferric Oxyhydroxide in Earth's Subducting Slabs

AbstractThe thermal stability of hydrous minerals in Earth's deep interior is key to understanding the evolution and physicochemical states of the planet. The recently discovered pyrite‐type (Py) FeO2Hx (x ≤ 1) phase, which can be transformed from α/ε‐FeOOH at ∼80 GPa, is believed to be a crucial candidate in transporting water and hydrogen to the lowermost mantle through subducting slabs. Here, we examined the stability and decomposition behavior of FeOOH through a set of shock‐recovery experiments up to ∼70 GPa and ∼2,750 K. Our results show that FeOOH partially decomposes to iron oxides Fe2O3 and Fe3O4 at 35–70 GPa and 1,150–2,750 K, which indicates that H2O and O2 are released during the decomposition of FeOOH in subducting slabs. The released H2O and O2 may have altered the physical and chemical properties of the surrounding mantle and contributed to the oxidation of surface Earth.

Vibrational and Thermodynamic Properties of Hydrous Iron-Bearing Lowermost Mantle Minerals

The vibrational and thermodynamic properties of minerals are key to understanding the phase stability and the thermal structure of the Earth’s mantle. In this study, we modeled hydrous iron-bearing bridgmanite (Brg) and post-perovskite (PPv) with different [Fe3+-H] defect configurations using first-principles calculations combined with quasi-harmonic approximations (QHA). Fe3+-H configurations can be vibrationally stable in Brg and PPv; the site occupancy of this defect will strongly affect its thermodynamic properties and particularly its response to pressure. The presence of Fe3+-H introduces distinctive high-frequency vibrations to the crystal. The frequency of these peaks is configuration dependence. Of the two defect configurations, [FeSi′+OH·] makes large effects on the thermodynamic properties of Brg and PPv, whereas [VMg″+FeMg·+OH·] has negligible effects. With an expected lower mantle water concentrations of <1000 wt. ppm the effect of Fe3+-H clusters on properties such as heat capacity and thermal expansion is negligible, but the effect on the Grüneisen parameter γ can be significant (~1.2%). This may imply that even a small amount of water may affect the anharmonicity of Fe3+-bearing MgSiO3 in lower mantle conditions and that when calculating the adiabaticity of the mantle, water concentrations need to be considered.

Sensitivity of SK(K)S and ScS phases to heterogeneous anisotropy in the lowermost mantle from global wavefield simulations

SUMMARY Observations of seismic anisotropy at the base of the mantle are abundant. Given recent progress in understanding how deformation relates to anisotropy in lowermost mantle minerals at the relevant pressure and temperature conditions, these observations can be used to test specific geodynamic scenarios, and have the potential to reveal patterns of flow at the base of the mantle. For example, several recent studies have sought to reproduce measurements of shear wave splitting due to D″ anisotropy using models that invoke specific flow and texture development geometries. A major limitation in such studies, however, is that the forward modelling is nearly always carried out using a ray theoretical framework, and finite-frequency wave propagation effects are not considered. Here we present a series of numerical wave propagation simulation experiments that explore the finite-frequency sensitivity of SKS, SKKS and ScS phases to laterally varying anisotropy at the base of the mantle. We build on previous work that developed forward modelling capabilities for anisotropic lowermost mantle models using the AxiSEM3D spectral element solver, which can handle arbitrary anisotropic geometries. This approach enables us to compute seismograms for relatively short periods (∼4 s) for models that include fully 3-D anisotropy at moderate computational cost. We generate synthetic waveforms for a suite of anisotropic models with increasing complexity. We first test a variety of candidate elastic tensors in laterally homogeneous models to understand how different lowermost mantle elasticity scenarios express themselves in shear wave splitting measurements. We then consider a series of laterally heterogeneous models of increasing complexity, exploring how splitting behaviour varies across the edges of anisotropic blocks and investigating the minimum sizes of anisotropic heterogeneities that can be reliably detected using SKS, SKKS and ScS splitting. Finally, we apply our modelling strategy to a previously published observational study of anisotropy at the base of the mantle beneath Iceland. Our results show that while ray theory is often a suitable approximation for predicting splitting, particularly for SK(K)S phases, full-wave effects on splitting due to lowermost mantle anisotropy can be considerable in some circumstances. Our simulations illuminate some of the challenges inherent in reliably detecting deep mantle anisotropy using body wave phases, and point to new strategies for interpreting SKS, SKKS and ScS waveforms that take full advantage of newly available computational techniques in seismology.

Sound Velocities in FeSi at Lower Mantle Conditions and the Origin of Ultralow‐Velocity Zones

AbstractThe origin of ultralow‐velocity zones (ULVZs) remains an open question despite recent advances in mineral physics and seismology. Here, we examine the hypothesis that FeSi formed from core‐mantle chemical reactions is a plausible source of ULVZs at the core‐mantle boundary (CMB). The sound velocities of B2‐structured FeSi were measured up to 115(2) GPa and 1600(200) K by nuclear inelastic scattering (NIS) in laser‐heated diamond anvil cells (LH‐DACs). Within uncertainties, the sound velocities of B2‐FeSi display negligible anharmonicity, hence validating the extrapolation of velocity‐density relations (Birch's law) to P‐T conditions of the CMB. The sound velocities of B2‐FeSi are significantly lower compared to other candidate phases in a lowermost mantle assemblage, and the Preliminary Reference Earth Model at CMB conditions. Less than 8.4 vol% of FeSi in the aggregate is thus sufficient to explain both the velocity decrements and the high density anomaly observed in a wide range of ULVZs.

Seismic Evidence for Thermal and Chemical Heterogeneities in D″ Region Beneath Central America From Grid Search Modeling

AbstractLarge lateral variations in the depth, sharpness, and magnitude of the D″ discontinuity are observed beneath Central America and the Caribbean using a new innovative grid search modeling method. The strong correlation between the D″ topography and the underlying seismic velocities suggests a large lateral gradient in the thermal condition in the lowermost mantle. Low‐velocity patches were identified among the high‐velocity areas, coinciding with regions showing a disruption in the D″ discontinuity. The low‐velocity zones may be attributed to the heat trapped in the middle of the Farallon slab remnant while those farther to the west are likely to be associated with the upwelling generated at the edge of the ancient slab. A broad range in the sharpness of the D″ discontinuity was also observed, suggesting the existence of different chemical components in the region, as a result of a disparate mixture of basalt, harzburgite, and pyrolite.

Comparing global seismic tomography models using varimax principal component analysis

Abstract. Global seismic tomography has greatly progressed in the past decades, with many global Earth models being produced by different research groups. Objective, statistical methods are crucial for the quantitative interpretation of the large amount of information encapsulated by the models and for unbiased model comparisons. Here we propose using a rotated version of principal component analysis (PCA) to compress the information in order to ease the geological interpretation and model comparison. The method generates between 7 and 15 principal components (PCs) for each of the seven tested global tomography models, capturing more than 97 % of the total variance of the model. Each PC consists of a vertical profile, with which a horizontal pattern is associated by projection. The depth profiles and the horizontal patterns enable examining the key characteristics of the main components of the models. Most of the information in the models is associated with a few features: large low-shear-velocity provinces (LLSVPs) in the lowermost mantle, subduction signals and low-velocity anomalies likely associated with mantle plumes in the upper and lower mantle, and ridges and cratons in the uppermost mantle. Importantly, all models highlight several independent components in the lower mantle that make between 36 % and 69 % of the total variance, depending on the model, which suggests that the lower mantle is more complex than traditionally assumed. Overall, we find that varimax PCA is a useful additional tool for the quantitative comparison and interpretation of tomography models.

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

Modelling Mie scattering in pyrolite in the laser-heated diamond anvil cell: Implications for the core-mantle boundary temperature determination

The core-mantle boundary (CMB) at ~2900 km depth and ~135 GPa is the interface between the solid silicate mantle and the liquid metallic outer core. The temperature at CMB (TCMB) is a key parameter in assessing the heat flow out of the core and is thus important for our understanding of the geodynamo and mantle convection. The upper bound on TCMB is constrained by the solidus temperature (onset of melting) of pyrolite, a chemical proxy of the mantle rock, at ~135 GPa because the lowermost mantle is predominantly solid. Previous laser-heated diamond anvil cell (LH DAC) studies reported pyrolite solidus temperatures of 3430–4180 K (± <200K) at P ~ 135 GPa (Andrault et al., 2011; Fiquet et al., 2010; Kim et al., 2020; Nomura et al., 2014). These temperatures were determined by pyrometry assuming that the optical properties of the sample are wavelength-independent in the visible to near-IR spectral range. Here we use Mie theory to show that this assumption is grossly inadequate for a polyphase sample with the grain sizes and refractive indices representative of those in a pyrolite compressed in the diamond anvil cell. We compute the wavelength-dependent light extinction coefficient of pyrolite samples and show that the previous melting studies may have underestimated pyrolite solidus temperature by up to 20–30%. This implies that the upper bound on the TCMB is of the order of ~5000 K.

A Compositional Component to the Samoa Ultralow‐Velocity Zone Revealed Through 2‐ and 3‐D Waveform Modeling of SKS and SKKS Differential Travel‐Times and Amplitudes

AbstractWe analyzed 4,754 broadband seismic recordings of the SKS, SKKS, and SPdKS wavefield from 13 high quality events sampling the Samoa ultralow‐velocity zone (ULVZ). We measured differential travel‐times and amplitudes between the SKKS and SKS arrivals, which are highly sensitive to the emergence of the SPdKS seismic phase, which is in turn highly sensitive to lowermost mantle velocity perturbations such as generated by ULVZs. We modeled these data using a 2‐D axi‐symmetric waveform modeling approach and are able to explain these data with a single ULVZ. In order to predict both travel‐time and amplitude perturbations we found that a large ULVZ length in the great circle arc direction on the order of 10° or larger is required. The large ULVZ length limits acceptable ULVZ elastic parameters. Here we find that δVS and δVP reductions from 20% to 22% and 15% to 17% respectively gives us the best fit, with a thickness of 26 km. Initial 3‐D modeling efforts do not recover the extremes in the differential measurements, demonstrating that 3‐D effects are important and must be considered in the future. However, the 3‐D modeling is generally consistent with the velocity reductions recovered with the 2‐D modeling. These velocity reductions are compatible with a compositional component to the ULVZ. Furthermore, geodynamic predictions for a compositional ULVZ that is moving predict a long linear shape similar to the shape of the Samoa ULVZ we confirm in this study.

The morphology, evolution and seismic visibility of partial melt at the core–mantle boundary: implications for ULVZs

SUMMARY Seismic observations indicate that the lowermost mantle above the core–mantle boundary (CMB) is strongly heterogeneous. Body waves reveal a variety of ultra-low velocity zones (ULVZs), which extend not more than 100 km above the CMB and have shear velocity reductions of up to 30 per cent. While the nature and origin of these ULVZs remain uncertain, some have suggested they are evidence of partial melting at the base of mantle plumes. Here we use coupled geodynamic/thermodynamic modelling to explore the hypothesis that present-day deep mantle melting creates ULVZs and introduces compositional heterogeneity in the mantle. Our models explore the generation and migration of melt in a deforming and compacting host rock at the base of a plume in the lowermost mantle. We test whether the balance of gravitational and viscous forces can generate partially molten zones that are consistent with the seismic observations. We find that for a wide range of plausible melt densities, permeabilities and viscosities, lower mantle melt is too dense to be stirred into convective flow and instead sinks down to form a completely molten layer, which is inconsistent with observations of ULVZs. Only if melt is less dense or at most ca. 1 per cent more dense than the solid, or if melt pockets are trapped within the solid, can melt remain suspended in the partial melt zone. In these cases, seismic velocities would be reduced in a cone at the base of the plume. Generally, we find partial melt alone does not explain the observed ULVZ morphologies and solid-state compositional variation is required to explain the anomalies. Our findings provide a framework for testing whether seismically observed ULVZ shapes are consistent with a partial melt origin, which is an important step towards constraining the nature of the heterogeneities in the lowermost mantle and their influence on the thermal, compositional and dynamic evolution of the Earth.

A tree of Indo-African mantle plumes imaged by seismic tomography

Mantle plumes were conceived as thin, vertical conduits in which buoyant, hot rock from the lowermost mantle rises to Earth’s surface, manifesting as hotspot-type volcanism far from plate boundaries. Spatially correlated with hotspots are two vast provinces of slow seismic wave propagation in the lowermost mantle, probably representing the heat reservoirs that feed plumes. Imaging plume conduits has proved difficult because most are located beneath the non-instrumented oceans, and they may be thin. Here we combine new seismological datasets to resolve mantle upwelling across all depths and length scales, centred on Africa and the Indian and Southern oceans. Using seismic waves that sample the deepest mantle extensively, we show that mantle upwellings are arranged in a tree-like structure. From a central, compact trunk below ~1,500 km depth, three branches tilt outwards and up towards various Indo-Austral hotspots. We propose that each tilting branch represents an alignment of vertically rising blobs or proto-plumes, which detached in a linear staggered sequence from their underlying low-velocity corridor at the core–mantle boundary. Once a blob reaches the viscosity discontinuity between lower and upper mantle, it spawns a ‘classical’ plume-head/plume-tail sequence.

Imaging paleoslabs and inferring the Clapeyron slope in D″ beneath the northern Pacific based on high-resolution inversion of seismic waveforms for 3-D transversely isotropic structure

We use waveform inversion to infer high-resolution 3-D models of the S-velocity perturbation δVS and the anisotropy parameter δξ in the lowermost 400 km of the mantle beneath the northern Pacific. Our inferred models show three prominent features: (i) a widespread positive δVS layer with ~100 km thickness that exists ~250 km above the core-mantle boundary (CMB), which could be due to a bridgmanite to post-perovskite phase transition related to the D″ discontinuity; (ii) distinct positive -δVS anomalies with positive δξ resolved from 100 km to more than 400 km above the CMB, which could be the subducted Izanagi, Farallon, and Telkhinia slabs; (iii) a vertically continuous low-δVS anomaly with negative δξ at the edge of the subducted slab, which we interpret as an upwelling plume induced by slab sinking. We infer the Clapeyron slope of the post-perovskite phase transition in the lowermost mantle beneath the northern Pacific to be 10.2 ± 1.3 MPa/K, based on the δVS gradient as a function of depth for the inferred 3-D VS structure. Our estimated large positive Clapeyron slope at the CMB suggests vigorous convection in the lowermost mantle.

Mantle plumes and their role in Earth processes

The existence of mantle plumes was first proposed in the 1970s to explain intra-plate, hotspot volcanism, yet owing to difficulties in resolving mantle upwellings with geophysical images and discrepancies in interpretations of geochemical and geochronological data, the origin, dynamics and composition of plumes and their links to plate tectonics are still contested. In this Review, we discuss progress in seismic imaging, mantle flow modelling, plate tectonic reconstructions and geochemical analyses that have led to a more detailed understanding of mantle plumes. Observations suggest plumes could be both thermal and chemical in nature, can attain complex and broad shapes, and that more than 18 plumes might be rooted in regions of the lowermost mantle. The case for a deep mantle origin is strengthened by the geochemistry of hotspot volcanoes that provide evidence for entrainment of deeply recycled subducted components, primordial mantle domains and, potentially, materials from Earth’s core. Deep mantle plumes often appear deflected by large-scale mantle flow, resulting in hotspot motions required to resolve past tectonic plate motions. Future research requires improvements in resolution of seismic tomography to better visualize deep mantle plume structures at smaller than 100-km scales. Concerted multi-proxy geochemical and dating efforts are also needed to better resolve spatiotemporal and chemical evolutions of long-lived mantle plumes.

Composition and Pressure Effects on Partitioning of Ferrous Iron in Iron-Rich Lower Mantle Heterogeneities

Both seismic observations of dense low shear velocity regions and models of magma ocean crystallization and mantle dynamics support enrichment of iron in Earth’s lowermost mantle. Physical properties of iron-rich lower mantle heterogeneities in the modern Earth depend on distribution of iron between coexisting lower mantle phases (Mg,Fe)O magnesiowüstite, (Mg,Fe)SiO3 bridgmanite, and (Mg,Fe)SiO3 post-perovskite. The partitioning of iron between these phases was investigated in synthetic ferrous-iron-rich olivine compositions (Mg0.55Fe0.45)2SiO4 and (Mg0.28Fe0.72)2SiO4 at lower mantle conditions ranging from 33–128 GPa and 1900–3000 K in the laser-heated diamond anvil cell. The resulting phase assemblages were characterized by a combination of in situ X-ray diffraction and ex situ transmission electron microscopy. The exchange coefficient between bridgmanite and magnesiowüstite decreases with pressure and bulk Fe# and increases with temperature. Thermodynamic modeling determines that incorporation and partitioning of iron in bridgmanite are explained well by excess volume associated with Mg-Fe exchange. Partitioning results are used to model compositions and densities of mantle phase assemblages as a function of pressure, FeO-content and SiO2-content. Unlike average mantle compositions, iron-rich compositions in the mantle exhibit negative dependence of density on SiO2-content at all mantle depths, an important finding for interpretation of deep lower mantle structures.

ScanArray—A Broadband Seismological Experiment in the Baltic Shield

Abstract The ScanArray international collaborative program acquired broadband seismological data at 192 locations in the Baltic Shield during the period between 2012 and 2017. The main objective of the program is to provide seismological constraints on the structure of the lithospheric crust and mantle as well as the sublithospheric upper mantle. The new information will be applied to studies of how the lithospheric and deep structure affect observed fast topographic change and geological-tectonic evolution of the region. The program also provides new information on local seismicity, focal mechanisms, and seismic noise. The recordings are generally of very high quality and are used for analysis by various seismological methods, including P- and S-wave receiver functions for the crust and upper mantle, surface wave and ambient noise inversion for seismic velocity, body-wave P- and S-wave tomography for upper mantle velocity structure using ray and finite frequency methods, and shear-wave splitting measurements for obtaining bulk anisotropy of the upper and lowermost mantle. Here, we provide a short overview of the data acquisition and initial analysis of the new data, together with an example of integrated seismological results obtained by the project group along a representative ∼1800-km-long profile across most of the tectonic provinces in the Baltic Shield between Denmark and the North Cape. The first models support a subdivision of the Paleoproterozoic Svecofennian province into three domains, where the highest topography of the Scandes mountain range in Norway along the Atlantic Coast has developed solely in the southern and northern domains, whereas the topography is more subdued in the central domain.

Elastic Properties of the Pyrite‐Type FeOOH‐AlOOH System From First‐Principles Calculations

AbstractThe stability, structure, and elastic properties of pyrite‐type (FeS2 structured) FeO2H were determined using density functional theory‐based computations with an internally consistent Coulombic self‐interaction term (Ueff). The properties of pyrite‐type FeO2H are compared to that of pyrite‐type AlO2H, with which it likely forms a solid solution at high temperature, as well as the respective lower pressure CaCl2‐type polymorphs of both endmembers: ϵ‐FeOOH and δ‐AlOOH. Due to substantial differences in the CaCl2‐type → pyrite‐type structural transition pressures of these endmembers, the stabilities of the (Al,Fe)O2H solid solution polymorphs are anticipated to be compositionally driven at lower mantle pressures. As the geophysical properties of (Al,Fe)OOH are structurally dependant, interpretations regarding the contribution of pyrite‐type FeO2H to seismically observed features must take into account the importance of this broad phase loop. With this in mind, Fe‐rich pyrite‐type (Al,Fe)OOH may coexist with Al‐dominant CaCl2‐type δ‐(Al,Fe)OOH in the deep Earth. Furthermore, pyrite‐type (Al0.5–0.6,Fe0.4–0.5)O2H can reproduce the reduced compressional and shear velocities characteristic of seismically observed ultra low velocity zones in the Earth's lowermost mantle while Al‐dominant but Fe‐bearing CaCl2‐type δ‐(Al,Fe)OOH may contribute to large low shear velocity provinces.

Ab initio lattice thermal conductivity of MgSiO3 across the perovskite-postperovskite phase transition

Lattice thermal conductivity ($\kappa_{lat}$) of MgSiO$_3$ postperovskite (MgPPv) under the Earth's lower mantle high pressure-temperature conditions is studied using the phonon quasiparticle approach by combing \textit{ab initio} molecular dynamics and lattice dynamics simulations. Phonon lifetimes are extracted from the phonon quasiparticle calculations, and the phonon group velocities are computed from the anharmonic phonon dispersions, which in principle capture full anharmonicity. It is found that throughout the lowermost mantle, including the D" region, $\kappa_{lat}$ of MgPPv is ~25% larger than that of MgSiO$_3$ perovskite (MgPv), mainly due to MgPPv's higher phonon velocities. Such a difference in phonon velocities between the two phases originates in the MgPPv's relatively smaller primitive cell. Systematic results of temperature and pressure dependences of both MgPPv's and MgPv's $\kappa_{lat}$ are demonstrated.

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% MgSiO3–10 mol% Fe2O3 system, rather than forming a homogeneous single phase. The Fe3+-rich bridgmanite has substantially lower velocities and a higher VP/VS ratio than MgSiO3 bridgmanite under lowermost-mantle conditions. Our modeling shows that the enrichment of Fe3+-rich bridgmanite in a pyrolitic composition can explain the observed features of the LLSVPs. The presence of Fe3+-rich materials within LLSVPs may have profound effects on the deep reservoirs of redox-sensitive elements and their isotopes.

A high-resolution map of Hawaiian ULVZ morphology from ScS phases

We use core reflected ScS waves sensitive to a broad region of the core-mantle boundary beneath Hawaii to create the first high-resolution map of the Hawaiian ultralow-velocity zone (ULVZ). Positive ScS-S differential times are used to identify regions of strong slow velocity anomalies in the lowermost mantle, and the presence of pre/post-cursors around the main ScS phase confirms the sharp top of a basal ULVZ layer. Pre/post-cursor arrivals are mapped into a volume across their region of sensitivity to produce a detailed image of ULVZ morphology. The variability observed across the ULVZ can be interpreted in terms of varying height or velocity reduction, but the large range of velocity variations required to explain observations suggests that most variability reflects varying layer thickness.The Hawaiian ULVZ is observed to be a large-scale regional feature of varying height (∼5-25 km) extending across a wide area along the edge of the Pacific large low-velocity province (LLVP). Variability in previous models of the Hawaiian ULVZ can be explained by studies imaging different parts of this strongly variable, large-scale feature. Maximum ULVZ thicknesses (no taller than 30 km) are found in a flat-topped, steep-sided region on the order of 1000 km in diameter located west of Hawaii. This feature is coincident with the previously identified Hawaiian mega-ULVZ, and is interpreted to represent the root of the Hawaiian plume, which is offset from the volcanic surface expression. A tall, asymmetric ridge of ULVZ material is observed to the east of Hawaii. Both regions of maximum ULVZ heights are bounded by the Pacific LLVP to the south, and fast seismic velocity anomalies interpreted as slab remnants to the north. This points to a potential geodynamical scenario where dense ULVZ material is pushed by subducted slab remnants into thicker piles against LLVP margins.

Reversal of carbonate-silicate cation exchange in cold slabs in Earth\u2019s lower mantle

The stable forms of carbon in Earth’s deep interior control storage and fluxes of carbon through the planet over geologic time, impacting the surface climate as well as carrying records of geologic processes in the form of diamond inclusions. However, current estimates of the distribution of carbon in Earth’s mantle are uncertain, due in part to limited understanding of the fate of carbonates through subduction, the main mechanism that transports carbon from Earth’s surface to its interior. Oxidized carbon carried by subduction has been found to reside in MgCO3 throughout much of the mantle. Experiments in this study demonstrate that at deep mantle conditions MgCO3 reacts with silicates to form CaCO3. In combination with previous work indicating that CaCO3 is more stable than MgCO3 under reducing conditions of Earth’s lowermost mantle, these observations allow us to predict that the signature of surface carbon reaching Earth’s lowermost mantle may include CaCO3.