Coupled dynamics and evolution of primordial and recycled heterogeneity in Earth’s lower mantle

The nature of compositional heterogeneity in Earth’s lower mantle remains a long-standing puzzle that can inform about the long-term thermochemical evolution and dynamics of our planet. Here, we use global-scale 2D models of thermo- chemical mantle convection to investigate the coupled evolution and mixing of (intrinsically-dense) recycled and (intrinsically- strong) primordial heterogeneity in the mantle. We explore the effects of ancient compositional layering of the mantle, as motivated by magma-ocean solidification studies, and of the physical parameters of primordial material. Depending on these physical parameters, our models predict various regimes of mantle evolution and heterogeneity preservation over 4.5 Gyrs. Over a wide parameter range, primordial and recycled heterogeneity are predicted to co-exist with each other in the lower mantle of Earth-like planets. Primordial material usually survives as mid-to-large scale blobs (or streaks) in the mid-mantle, around 1000-2000 km depth, and this preservation is largely independent on the initial primordial-material volume. In turn, recycled oceanic crust (ROC) persists as large piles at the base of the mantle and as small streaks everywhere else. In models with an additional dense FeO-rich layer initially present at the base of the mantle, the ancient dense material partially survives at the top of ROC piles, causing the piles to be compositionally stratified. Moreover, the addition of such an ancient FeO-rich basal layer significantly aids the preservation of the viscous domains in the mid-mantle. Finally, we find that primordial blobs are commonly directly underlain by thick ROC piles, and aid their longevity and stability. Based on our results, we propose an integrated style of mantle heterogeneity for the Earth, involving the preservation of primordial domains along with recycled piles. This style has important implications for early Earth evolution, and has the potential of reconciling geophysical and geochemical discrepancies on present-day lower-mantle heterogeneity.

Variable dynamic styles of primordial heterogeneity preservation in the Earth’s lower mantle

The evolution of the system Earth is critically influenced by the long-term dynamics, composition and structure of the mantle. While cosmochemical and geochemical constraints indicate that the lower mantle hosts an ancient primordial reservoir that may be enriched in SiO2 with respect to the upper mantle, geophysical observations and models point to efficient mass transfer and convective mixing across the entire mantle. Recent hypotheses of primordial-material preservation in a convecting mantle involve delayed mixing of intrinsically dense and/or intrinsically strong heterogeneity. Yet, the effects of composition-dependent rheology and density upon heterogeneity preservation and the dynamics of mantle mixing remain poorly understood. Here, we present two-dimensional numerical models in spherical geometry, investigating the preservation styles of primordial material as a function of its physical properties (i.e., viscosity and density contrasts). We establish multiple regimes of primordial-material preservation that can occur in terrestrial planets. These include (1) efficient mixing, (2) double-layered convection with or without topography, and (3) variable styles of partial heterogeneity preservation (e.g., as diffuse domains, piles or viscous blobs in the lower mantle). Some of these regimes are here characterised for the first time, and all regimes are put into context with each other as a function of model parameters. The viscous-blobs and diffuse-domains regimes can reconcile the preservation of primordial domains in a convecting mantle, potentially resolving the discrepancy between geochemical and geophysical constraints for planet Earth. Several, if not all, regimes characterised here may be relevant to understand the long-term evolution of terrestrial planets in general.

Temperature distribution in the three-layered upper mantle beneath a continent

Temperature distribution in the upper mantle underneath the continent, as well as temperature distribution in the lower mantle, is obtained. In the continental lithosphere, the solution to the heat transfer equation is obtained in the model of conduction heat transfer with inner heat within the crust. To calculate the temperature distribution in the upper and lower mantle, we use the results of laboratory and theoretical modeling of free convective heat transfer in a horizontal layer heated from below and cooled from above.

Incorporation mechanism of Fe and Al into bridgmanite in a subducting mid-ocean ridge basalt and its crystal chemistry

The compositional difference between subducting slabs and their surrounding lower-mantle can yield the difference in incorporation mechanism of Fe and Al into bridgmanite between both regions, which should cause heterogeneity in physical properties and rheology of the lower mantle. However, the precise cation-distribution has not been examined in bridgmanites with Fe- and Al-contents expected in a mid-ocean ridge basalt component of subducting slabs. Here we report on Mg0.662Fe0.338Si0.662Al0.338O3 bridgmanite single-crystal characterized by a combination of single-crystal X-ray diffraction, synchrotron 57Fe-Mössbauer spectroscopy and electron probe microanalysis. We find that the charge-coupled substitution AMg2+ + BSi4+ ↔ AFe3+(high-spin) + BAl3+ is predominant in the incorporation of Fe and Al into the practically eightfold-coordinated A-site and the sixfold-coordinated B-site in bridgmanite structure. The incorporation of both cations via this substitution enhances the structural distortion due to the tilting of BO6 octahedra, yielding the unusual expansion of mean <A–O> bond-length due to flexibility of A–O bonds for the structural distortion, in contrast to mean <B–O> bond-length depending reasonably on the ionic radius effect. Moreover, we imply the phase-transition behavior and the elasticity of bridgmanite in slabs subducting into deeper parts of the lower mantle, in terms of the relative compressibility of AO12 (practically AO8) and BO6 polyhedra.