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.

Electronic, Structural, and Mechanical Properties of SiO_{2} Glass at High Pressure Inferred from its Refractive Index

We report the first direct measurements of the refractive index of silica glass up to 145 GPa that allowed quantifying its density, bulk modulus, Lorenz-Lorentz polarizability, and band gap. These properties show two major anomalies at ∼10 and ∼40  GPa. The anomaly at ∼10  GPa signals the onset of the increase in Si coordination, and the anomaly at ∼40  GPa corresponds to a nearly complete vanishing of fourfold Si. More generally, we show that the compressibility and density of noncrystalline solids can be accurately measured in simple optical experiments up to at least 110 GPa.

The Critical Point and the Supercritical State of Alkali Feldspars: Implications for the Behavior of the Crust During Impacts

The position of the vapor-liquid dome and of the critical point determine the evolution of the outermost parts of the protolunar disk during cooling and condensation after the Giant Impact. The parts of the disk in supercritical or liquid state evolve as a single thermodynamic phase; when the thermal trajectory of the disk reaches the liquid-vapor dome, gas and melt separate leading to heterogeneous convection and phase separation due to friction. Different layers of the proto-Earth behaved differently during the Giant Impact depending on their constituent materials and initial thermodynamic conditions. Here we use first-principles molecular dynamics to determine the position of the critical point for NaAlSi3O8 and KAlSi3O8 feldspars, major minerals of the Earth and Moon crusts. The variations of the pressure calculated at various volumes along isotherms yield the position of the critical points: 0.5-0.8 g cm-3 and 5500-6000 K range for the Na-feldspar, 0.5-0.9 g cm-3 and 5000-5500 K range for the K-feldspar. The simulations suggest that the vaporization is incongruent, with a degassing of O2 starting at 4000 K and gas component made mostly of free Na and K cations, O2, SiO and SiO2 species for densities below 1.5 g cm-3. The Hugoniot equations of state imply that low-velocity impactors (<8.3 km s-1) would at most melt a cold feldspathic crust, whereas large impacts in molten crust would see temperatures raise up to 30000 K.

Density of NaAlSi2O6 Melt at High Pressure and Temperature Measured by In-Situ X-ray Microtomography

In this study, the volumetric compression of jadeite (NaAlSi2O6) melt at high pressures was determined by three-dimensional volume imaging using the synchrotron-based X-ray microtomography technique in a rotation-anvil device. Combined with the sample mass, measured using a high-precision analytical balance prior to the high-pressure experiment, the density of jadeite melt was obtained at high pressures and high temperatures up to 4.8 GPa and 1955 K. The density data were fitted to a third-order Birch-Murnaghan equation of state, resulting in a best-fit isothermal bulk modulus K T 0 of 10.8 − 5.3 + 1.9 GPa and its pressure derivative K T 0 ′ of 3.4 − 0.4 + 6.6 . Comparison with data for silicate melts of various compositions from the literature shows that alkali-rich, polymerized melts are generally more compressible than alkali-poor, depolymerized ones. The high compressibility of jadeite melt at high pressures implies that polymerized sodium aluminosilicate melts, if generated by low-degree partial melting of mantle peridotite at ~250–400 km depth in the deep upper mantle, are likely denser than surrounding mantle materials, and thus gravitationally stable.

High-Pressure Sound Velocity Measurements of Liquids Using In Situ Ultrasonic Techniques in a Multianvil Apparatus

Sound velocity and equation of state of liquids provide important constraints on the generation, presence, and transport of silicate and metallic melts in the Earth’s interior. Unlike their solid counterparts, these properties of liquids pose great technical challenges to high-pressure measurements and are poorly constrained. Here we present the technical developments that have been made at the GSECARS beamline 13-ID-D of the Advanced Photon Source for the past several years for determination of sound velocity of liquids using the ultrasonic techniques in a 1000-ton Kawai-type multianvil apparatus. Temperature of the sound velocity measurements has been extended to ~2400 K at 4 GPa and ~2000 K at 8 GPa to enable studies of liquids with very high melting temperatures, such as the silicate liquids.

Structure and density of basaltic melts at mantle conditions from first-principles simulations

The origin and stability of deep-mantle melts, and the magmatic processes at different times of Earth's history are controlled by the physical properties of constituent silicate liquids. Here we report density functional theory-based simulations of model basalt, hydrous model basalt and near-MORB to assess the effects of iron and water on the melt structure and density, respectively. Our results suggest that as pressure increases, all types of coordination between major cations and anions strongly increase, and the water speciation changes from isolated species to extended forms. These structural changes are responsible for rapid initial melt densification on compression thereby making these basaltic melts possibly buoyantly stable at one or more depths. Our finding that the melt-water system is ideal (nearly zero volume of mixing) and miscible (negative enthalpy of mixing) over most of the mantle conditions strengthens the idea of potential water enrichment of deep-mantle melts and early magma ocean.

Development of chemical and topological structure in aluminosilicate liquids and glasses at high pressure

The high pressure structure of liquid and glassy anorthite (CaAl(2)Si(2)O(8)) and calcium aluminate (CaAl(2)O(4)) glass was measured by using in situ synchrotron x-ray diffraction in a diamond anvil cell up to 32.4(2) GPa. The results, combined with ab initio molecular dynamics and classical molecular dynamics simulations using a polarizable ion model, reveal a continuous increase in Al coordination by oxygen, with 5-fold coordinated Al dominating at 15 GPa and a preponderance of 6-fold coordinated Al at higher pressures. The development of a peak in the measured total structure factors at 3.1 Å(-1) is interpreted as a signature of changes in topological order. During compression, cation-centred polyhedra develop edge- and face- sharing networks. Above 10 GPa, following the pressure-induced breakdown of the network structure, the anions adopt a structure similar to a random close packing of hard spheres.

Structural change in molten basalt at deep mantle conditions

Silicate liquids play a key part at all stages of deep Earth evolution, ranging from core and crust formation billions of years ago to present-day volcanic activity. Quantitative models of these processes require knowledge of the structural changes and compression mechanisms that take place in liquid silicates at the high pressures and temperatures in the Earth's interior. However, obtaining such knowledge has long been impeded by the challenging nature of the experiments. In recent years, structural and density information for silica glass was obtained at record pressures of up to 100 GPa (ref. 1), a major step towards obtaining data on the molten state. Here we report the structure of molten basalt up to 60 GPa by means of in situ X-ray diffraction. The coordination of silicon increases from four under ambient conditions to six at 35 GPa, similar to what has been reported in silica glass. The compressibility of the melt after the completion of the coordination change is lower than at lower pressure, implying that only a high-order equation of state can accurately describe the density evolution of silicate melts over the pressure range of the whole mantle. The transition pressure coincides with a marked change in the pressure-evolution of nickel partitioning between molten iron and molten silicates, indicating that melt compressibility controls siderophile-element partitioning.