Tidal tomography constrains Earth's deep-mantle buoyancy

Volcanic arc rigidity variations illuminated by coseismic deformation of the 2011 Tohoku-oki M9

Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales.

GIA imaging of 3-D mantle viscosity based on palaeo sea level observations - Part I: Sensitivity kernels for an Earth with laterally varying viscosity

A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth's viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data.

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-I[Formula: see text]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-I[Formula: see text]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.

Vestiges of impact-driven three-phase mixing in the chemistry and structure of Earth's mantle

Highly siderophile elements (HSEs; namely Ru, Rh, Pd, Re, Os, Ir, Pt, and Au) in Earth's mantle require the addition of metals after the formation of Earth's core. Early, large collisions have the potential to deliver metals, but the details of their mixing with Earth's mantle remain unresolved. As a large projectile disrupts and penetrates Earth's mantle, a fraction of its metallic core may directly merge with Earth's core. Ensuing gravitational instabilities remove the remaining projectile's core stranded in Earth's mantle, leaving the latter deprived of HSEs. Here, we propose a framework that can efficiently retain the metallic components during large impacts. The mechanism is based on the ubiquitous presence of a partially molten region in the mantle beneath an impact-generated magma ocean, and it involves rapid three-phase flow with solid silicate, molten silicate, and liquid metal as well as long-term mixing by mantle convection. In addition, large low-shear-velocity provinces in the lower mantle may originate from compositional heterogeneities resulting from the proposed three-phase flow during high-energy collisions.

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.

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

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

Assembly of the basal mantle structure beneath Africa

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

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

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

The influence of δ-(Al,Fe)OOH on seismic heterogeneities in Earth's lower mantle

The high-pressure phases of oxyhydroxides (δ-AlOOH, ε-FeOOH, and their solid solution), candidate components of subducted slabs, have wide stability fields, thus potentially influencing volatile circulation and dynamics in the Earth's lower mantle. Here, we report the elastic wave velocities of δ-(Al,Fe)OOH (Fe/(Al + Fe) = 0.13, δ-Fe13) to 79 GPa, determined by nuclear resonant inelastic X-ray scattering. At pressures below 20 GPa, a softening of the phonon spectra is observed. With increasing pressure up to the Fe3+ spin crossover (~ 45 GPa), the Debye sound velocity (vD) increases. At higher pressures, the low spin δ-Fe13 is characterized by a pressure-invariant vD. Using the equation of state for the same sample, the shear-, compressional-, and bulk-velocities (vS, vP, and vΦ) are calculated and extrapolated to deep mantle conditions. The obtained velocity data show that δ-(Al,Fe)OOH may cause low-vΦ and low-vP anomalies in the shallow lower mantle. At deeper depths, we find that this hydrous phase reproduces the anti-correlation between vS and vΦ reported for the large low seismic velocity provinces, thus serving as a potential seismic signature of hydrous circulation in the lower mantle.

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

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

Ancient helium and tungsten isotopic signatures preserved in mantle domains least modified by crustal recycling

Rare high-³He/⁴He signatures in ocean island basalts (OIB) erupted at volcanic hotspots derive from deep-seated domains preserved in Earth's interior. Only high-³He/⁴He OIB exhibit anomalous ¹⁸²W-an isotopic signature inherited during the earliest history of Earth-supporting an ancient origin of high ³He/⁴He. However, it is not understood why some OIB host anomalous ¹⁸²W while others do not. We provide geochemical data for the highest-³He/⁴He lavas from Iceland (up to 42.9 times atmospheric) with anomalous ¹⁸²W and examine how Sr-Nd-Hf-Pb isotopic variations-useful for tracing subducted, recycled crust-relate to high ³He/⁴He and anomalous ¹⁸²W. These data, together with data on global OIB, show that the highest-³He/⁴He and the largest-magnitude ¹⁸²W anomalies are found only in geochemically depleted mantle domains-with high ¹⁴³Nd/¹⁴⁴Nd and low ²⁰⁶Pb/²⁰⁴Pb-lacking strong signatures of recycled materials. In contrast, OIB with the strongest signatures associated with recycled materials have low ³He/⁴He and lack anomalous ¹⁸²W. These observations provide important clues regarding the survival of the ancient He and W signatures in Earth's mantle. We show that high-³He/⁴He mantle domains with anomalous ¹⁸²W have low W and ⁴He concentrations compared to recycled materials and are therefore highly susceptible to being overprinted with low ³He/⁴He and normal (not anomalous) ¹⁸²W characteristic of subducted crust. Thus, high ³He/⁴He and anomalous ¹⁸²W are preserved exclusively in mantle domains least modified by recycled crust. This model places the long-term preservation of ancient high ³He/⁴He and anomalous ¹⁸²W in the geodynamic context of crustal subduction and recycling and informs on survival of other early-formed heterogeneities in Earth's interior.

LoadDef: A Python-Based Toolkit to Model Elastic Deformation Caused by Surface Mass Loading on Spherically Symmetric Bodies

Temporal variations of surface masses, such as the hydrosphere and atmosphere of the Earth, load the surfaces of planetary bodies causing temporal variations in deformation. Surface shear forces and gravitational fields also drive ongoing planetary deformation. Characterizing the spatiotemporal patterns of planetary deformation can constrain allowable models for the interior structure of a planetary body as well as for the distribution of surface and body forces. Pertinent applications include hydrology, glaciology, geodynamics, atmospheric science, and climatology. To address the diversity of emerging applications, we introduce a software suite called LoadDef that provides a collection of modular functions for modeling planetary deformation within a self-consistent, Python-based computational framework. Key features of LoadDef include computation of real-valued potential, load, and shear Love numbers for self-gravitating and spherically symmetric planetary models; computation of Love-number partial derivatives with respect to planetary density and elastic structure; computation of displacement, gravity, tilt, and strain load Green's functions; and computation of three-component surface displacements induced by surface mass loading. At a most basic level, only a planetary-structure model and a mass-load model must be supplied as input to LoadDef to utilize all the main features of the software. The end-to-end forward-modeling capabilities for mass-loading applications lay the foundation for sensitivity studies and geodetic tomography. LoadDef results have been validated with Global Navigation Satellite System observations and verified against independent software and published results. As a case study, we use LoadDef to predict the solid Earth's elastic response to ocean tidal loading across the western United States.

The adjoint equations for thermochemical compressible mantle convection: derivation and verification by twin experiments

The adjoint method is an efficient way to obtain gradient information in a mantle convection model relative to past flow structure, allowing one to retrodict mantle flow from observations of the present-day mantle state. While adjoint equations for isochemical mantle flow have been derived for both incompressible and compressible flows, here we extend the method to thermochemical mantle flow models, and present thermochemical adjoint equations in the elastic-liquid approximation. We verify the method with twin experiments, and retrodict the flow history of a thermochemical reference model (reference twin) assuming for the final state, either a consistent thermochemical interpretation, using the thermochemical adjoint equations, or an inconsistent purely thermal interpretation, using the isochemical adjoint equations. The consistent simulation correctly retrodicts the flow evolution of the reference twin. The inconsistent case, instead, restores a false flow history whereby internal buoyancy forces and convectively maintained topography are overestimated. Because the cost function is reduced in either case, our results suggest that the adjoint method can be used to link assumptions on the role of chemical mantle heterogeneity to geologic inferences of dynamic topography, thus providing additional means to test hypotheses on mantle composition and dynamics.

Valence and spin states of iron are invisible in Earth's lower mantle

Heterogeneity in Earth's mantle is a record of chemical and dynamic processes over Earth's history. The geophysical signatures of heterogeneity can only be interpreted with quantitative constraints on effects of major elements such as iron on physical properties including density, compressibility, and electrical conductivity. However, deconvolution of the effects of multiple valence and spin states of iron in bridgmanite (Bdg), the most abundant mineral in the lower mantle, has been challenging. Here we show through a study of a ferric-iron-only (Mg_0.46Fe³⁺_0.53)(Si_0.49Fe³⁺_0.51)O₃ Bdg that Fe³⁺ in the octahedral site undergoes a spin transition between 43 and 53 GPa at 300 K. The resolved effects of the spin transition on density, bulk sound velocity, and electrical conductivity are smaller than previous estimations, consistent with the smooth depth profiles from geophysical observations. For likely mantle compositions, the valence state of iron has minor effects on density and sound velocities relative to major cation composition.