Formation of bridgmanite-enriched layer at the top lower-mantle during magma ocean solidification

Effect of bridgmanite-ferropericlase grain size evolution on Earth's average mantle viscosity: implications for mantle convection in early and present-day Earth

Recent experimental investigations of grain size evolution in bridgmanite-ferropericlase assemblages have suggested very slow growth for these bimodal phases. Despite numerous speculations on grain size-dependent viscosity, a comprehensive test with realistic grain size evolution parameters compatible with the lower mantle has been lacking. In this study, we develop self-consistent 2-D spherical half-annulus geodynamic models of Earth's evolution using the finite volume code StagYY to assess the role of grain size on lower mantle viscosity. We explore several models with and without grain size evolution to compare their effects on mantle viscosity. In models with grain size evolution, we consider three scenarios: (1) uniform grain growth throughout the entire mantle with a composite rheology, (2) different grain growth in the upper and lower mantle with a composite rheology, and (3) different grain growth in the upper and lower mantle with purely diffusion creep rheology. In the case of different grain size evolution, the upper mantle's grain size evolution law is controlled by forsterite-enstatite grain growth, while the lower mantle's grain size evolution law is controlled by bridgmanite-ferropericlase grain growth. Our results suggest that mantle viscosity is primarily controlled by temperature, whereas grain size has a minor effect compared to the effect of temperature. We attribute two primary reasons for this: First, the bridgmanite-ferropericlase growth is very slow in the lower mantle and the grain size variation is too small to significantly alter the mantle viscosity. Secondly, if grains grow too fast, thus the mantle deforms in the dislocation creep regime, making viscosity grain size-independent. To establish the robustness of this finding we vary several other model parameters, such as surface yield strength, phase transition grain size reset, different transitional stresses for creep mechanisms, pressure dependence on grain growth, and different grain damage parameters. For all our models, we consistently find that grain size has a very limited effect on controlling lower mantle viscosity in the present-day Earth. However, large grain size may have affected the lower mantle viscosity in the early Earth as larger grains of single phase bridgmanite could increase the viscosity of the early mantle delaying the onset of global convection.

Atomic and Electronic Structure in MgO-SiO2

Understanding disordered structure is difficult due to insufficient information in experimental data. Here, we overcome this issue by using a combination of diffraction and simulation to investigate oxygen packing and network topology in glassy (g-) and liquid (l-) MgO-SiO2 based on a comparison with the crystalline topology. We find that packing of oxygen atoms in Mg2SiO4 is larger than that in MgSiO3, and that of the glasses is larger than that of the liquids. Moreover, topological analysis suggests that topological similarity between crystalline (c)- and g-(l-) Mg2SiO4 is the signature of low glass-forming ability (GFA), and high GFA g-(l-) MgSiO3 shows a unique glass topology, which is different from c-MgSiO3. We also find that the lowest unoccupied molecular orbital (LUMO) is a free electron-like state at a void site of magnesium atom arising from decreased oxygen coordination, which is far away from crystalline oxides in which LUMO is occupied by oxygen's 3s orbital state in g- and l-MgO-SiO2, suggesting that electronic structure does not play an important role to determine GFA. We finally concluded the GFA of MgO-SiO2 binary is dominated by the atomic structure in terms of network topology.

Variation in bridgmanite grain size accounts for the mid-mantle viscosity jump

A viscosity jump of one to two orders of magnitude in the lower mantle of Earth at 800-1,200-km depth is inferred from geoid inversions and slab-subducting speeds. This jump is known as the mid-mantle viscosity jump1,2. The mid-mantle viscosity jump is a key component of lower-mantle dynamics and evolution because it decelerates slab subduction3, accelerates plume ascent4 and inhibits chemical mixing5. However, because phase transitions of the main lower-mantle minerals do not occur at this depth, the origin of the viscosity jump remains unknown. Here we show that bridgmanite-enriched rocks in the deep lower mantle have a grain size that is more than one order of magnitude larger and a viscosity that is at least one order of magnitude higher than those of the overlying pyrolitic rocks. This contrast is sufficient to explain the mid-mantle viscosity jump1,2. The rapid growth in bridgmanite-enriched rocks at the early stage of the history of Earth and the resulting high viscosity account for their preservation against mantle convection5-7. The high Mg:Si ratio of the upper mantle relative to chondrites8, the anomalous 142Nd:144Nd, 182W:184W and 3He:4He isotopic ratios in hot-spot magmas9,10, the plume deflection4 and slab stagnation in the mid-mantle3 as well as the sparse observations of seismic anisotropy11,12 can be explained by the long-term preservation of bridgmanite-enriched rocks in the deep lower mantle as promoted by their fast grain growth.

Insights into magma ocean dynamics from the transport properties of basaltic melt

The viscosity of magma plays a crucial role in the dynamics of the Earth: from the crystallization of a magma ocean during its initial stages to modern-day volcanic processes. However, the pressure-dependence behavior of viscosity at high pressure remains controversial. In this study, we report the results of first-principles molecular dynamics simulations of basaltic melt to show that the melt viscosity increases upon compression along each isotherm for the entire lower mantle after showing minima at ~6 GPa. However, elevated temperatures of the magma ocean translate to a narrow range of viscosity, i.e., 0.01-0.03 Pa.s. This low viscosity implies that the crystallization of the magma ocean could be complete within a few million years. These results also suggest that the crystallization of the magma ocean is likely to be fractional, thus supporting the hypothesis that present-day mantle heterogeneities could have been generated during the early crystallization of the primitive mantle.

Deep mantle origin of large igneous provinces and komatiites

Large igneous provinces (LIPs) resulted from intraplate magmatic events mobilizing volumes of magma up to several million cubic kilometers. LIPs and lavas with deep mantle sources have compositions ranging from komatiites found in Archean greenstone belts to basalts and picrites in Phanerozoic flood basalt and recent oceanic islands. In this study, we identify the mantle conditions appropriate to each type of lava based on an experimental study of the melting of pyrolite. The depth of the mantle source decreases from 600 to 700 km for the oldest komatiites to 100 to 300 km for picrites and basalts, and the extent of mantle melting ranges from 10 to 50%. We develop a geodynamical model that explains the origin of the hot mantle plumes capable of generating these melting P-T conditions. Within a superadiabatic temperature gradient persisting in the deep mantle, the ascent of hot mantle plumes creates excess temperatures up to 250 to 300 K by adiabatic decompression.

Inefficient Water Degassing Inhibits Ocean Formation on Rocky Planets: An Insight from Self-Consistent Mantle Degassing Models

A sufficient amount of water is required at the surface to develop water oceans. A significant fraction of water, however, remains in the mantle during magma ocean solidification, and thus the existence of water oceans is not guaranteed even for exoplanets located in the habitable zone. To discuss the likelihood of ocean formation, we built two models to predict the rate of mantle degassing during the magma ocean stage and the subsequent solid-state convection stage. We find that planets with low H₂O/CO₂ ratios would not have a sufficient amount of surface water to develop water oceans immediately after magma ocean solidification, and the majority of the water inventory would be retained in the mantle during their subsequent evolution regardless of planetary size. This is because oceanless planets are likely to operate under stagnant lid convection, and for such planets, dehydration stiffening of the depleted lithospheric mantle would limit the rate of mantle degassing. In contrast, a significant fraction of CO₂ would already be degassed during magma ocean solidification. With a strong greenhouse effect, all surface water would exist as vapor, and water oceans may be absent throughout planetary evolution. Volatile concentrations in the bulk silicate Earth are close to the threshold amount for ocean formation, so if Venus shared similar concentrations, small differences in solar radiation may explain the divergent evolutionary paths of Earth and Venus.

Constraining composition and temperature variations in the mantle transition zone

The mantle transition zone connects two major layers of Earth’s interior that may be compositionally distinct: the upper mantle and the lower mantle. Wadsleyite is a major mineral in the upper mantle transition zone. Here, we measure the single-crystal elastic properties of hydrous Fe-bearing wadsleyite at high pressure-temperature conditions by Brillouin spectroscopy. Our results are then used to model the global distribution of wadsleyite proportion, temperature, and water content in the upper mantle transition zone by integrating mineral physics data with global seismic observations. Our models show that the upper mantle transition zone near subducted slabs is relatively cold, enriched in wadsleyite, and slightly more hydrated compared to regions where plumes are expected. This study provides direct evidence for the thermochemical heterogeneities in the upper mantle transition zone which is important for understanding the material exchange processes between the upper and lower mantle.

A new study by @JinZhang_MP models the global distribution of wadsleyite proportion, temperature and water content in the upper mantle transition zone.

Destabilization of deep oxidized mantle drove the Great Oxidation Event

The rise of Earth's atmospheric O₂ levels at ~2.4 Ga was driven by a shift between increasing sources and declining sinks of oxygen. Here, we compile recent evidence that the mantle shows a significant increase in oxidation state leading to the Great Oxidation Event (GOE), linked to sluggish upward mixing of a deep primordial oxidized layer. We simulate this scenario by implementing a new rheological model for this oxidized, bridgmanite-enriched viscous material and demonstrate slow mantle mixing in simulations of early Earth's mantle. The eventual homogenization of this layer may take ~2 Ga, in line with the timing of the observed mantle redox shift, and would result in the increase in upper mantle oxidation of >1 log(fO₂) unit. Such a shift would alter the redox state of volcanic degassing products to more oxidized species, removing a major sink of atmospheric O₂ and allowing oxygen levels to rise at ~2.4 Ga.

Simultaneous generation of ultrahigh pressure and temperature to 50 GPa and 3300 K in multi-anvil apparatus

We attempted to generate ultrahigh pressure and temperature simultaneously using a multi-anvil apparatus by combining the technologies of ultrahigh-pressure generation using sintered diamond (SD) anvils, which can reach 120 GPa, and ultrahigh-temperature generation using a boron-doped diamond (BDD) heater, which can reach 4000 K. Along with this strategy, we successfully generated a temperature of 3300 K and a pressure of above 50 GPa simultaneously. Although the high hardness of BDD significantly prevents high-pressure generation at low temperatures, its high-temperature softening allows for effective pressure generation at temperatures above 1200 K. High temperature also enhances high-pressure generation because of the thermal pressure. We expect to generate even higher pressure in the future by combining SD anvils and a BDD heater with advanced multi-anvil technology.

Low-spin ferric iron in primordial bridgmanite crystallized from a deep magma ocean

The crystallization of the magma ocean resulted in the present layered structure of the Earth's mantle. An open question is the electronic spin state of iron in bridgmanite (the most abundant mineral on Earth) crystallized from a deep magma ocean, which has been neglected in the crystallization history of the entire magma ocean. Here, we performed energy-domain synchrotron Mössbauer spectroscopy measurements on two bridgmanite samples synthesized at different pressures using the same starting material (Mg_0.78Fe_0.13Al_0.11Si_0.94O₃). The obtained Mössbauer spectra showed no evidence of low-spin ferric iron (Fe³⁺) from the bridgmanite sample synthesized at relatively low pressure of 25 gigapascals, while that directly synthesized at a higher pressure of 80 gigapascals contained a relatively large amount. This difference ought to derive from the large kinetic barrier of Fe³⁺ rearranging from pseudo-dodecahedral to octahedral sites with the high-spin to low-spin transition in experiments. Our results indicate a certain amount of low-spin Fe³⁺ in the lower mantle bridgmanite crystallized from an ancient magma ocean. We therefore conclude that primordial bridgmanite with low-spin Fe³⁺ dominated the deeper part of an ancient lower mantle, which would contribute to lower mantle heterogeneity preservation and call for modification of the terrestrial mantle thermal evolution scenarios.

Machinable boron-doped diamond as a practical heating element in multi-anvil apparatuses

Being refractory and X-ray transparent, a boron-doped diamond (BDD) heater is considered an ideal heating element in a multi-anvil apparatus under diamond-stable pressures. However, the extremely high hardness of diamond makes it difficult to manufacture a BDD tube, which, in turn, hinders the wide application of BDD heaters in multi-anvil apparatuses. Here, I sintered a machinable BDD (MBDD) from a mixture of BDD powder and pitch (C_nH_2n+2) by its annealing in Ar at 1273 K for 5 h. The BDD powder was bound by a small amount of graphite (<10 wt. %) during the sintering process. Tubes (such as 1.2/0.7/4.0 mm in outer/inner diameter/length) can be manufactured from the MBDD block using a lathe or a computer numerical control machine. Due to the low content of graphite in MBDD, the graphite-diamond conversion has a small effect on heating performance. The MBDD heater shows a comparable performance in ultrahigh temperature generation with a high-pressure synthesized BDD heater by generating a temperature higher than 3300 K and melted Al₂O₃ under a pressure of 15 GPa. With good heating performance and excellent machinability, MBDD is a practical heating element in multi-anvil apparatuses. The achievement of stable temperature generation over 3300 K by the MBDD heater enables various measurements on the physicochemical properties of melts under the Earth's mantle conditions.

A strip-type boron-doped diamond heater synthesized by chemical vapor deposition for large-volume presses

We have developed a high-pressure furnace assembly with a commercially available chemical-vapor-deposition synthesized boron-doped diamond heater consisting of four strips for large-volume multi-anvil presses (LVPs). This assembly successfully generated temperatures up to 2990 K at 15 GPa. It also has highly reproducible power-temperature relations, enabling us to estimate temperature from power reliably. It can be used for experiments above 9 GPa and is particularly useful for synchrotron x-ray experiments because of the x-ray transparency. It is also competitive in price. This technique is, thus, practical in various LVP experiments in the diamond-stability field.

Boron-MgO composite as an X-ray transparent pressure medium in the multi-anvil apparatus

X-ray transparent materials are very beneficial for in situ X-ray experiments in the multi-anvil apparatus. We sintered machinable blocks of boron-MgO composites at 800-1000 °C under atmospheric pressure from a mixture of amorphous boron and brucite or Mg(OH)₂. The machinability of composite blocks improved with an increase in the brucite content in the starting material; a brucite content higher than 15 wt. % showed reasonable machinability in forming various shapes such as octahedron, cylinder, and sleeve. We confirmed the feasibility of the boron-MgO pressure medium by successfully generating lower mantle pressure (>23 GPa); its pressure generation efficiency is comparable to that of a Cr₂O₃ doped MgO pressure medium. The boron-MgO composite is expected to be an excellent thermal insulator owing to the extremely low thermal conductivity of amorphous boron; we confirmed its better thermal insulation performance through a comparative heating test with a zirconia sleeve in a Cr₂O₃ doped MgO pressure medium. Constituting light elements, the boron-MgO composite has high X-ray transparency, which enables us to conduct various cutting edge X-ray measurements in the large volume multi-anvil apparatus.