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.

Behavior of Long-Chain Hydrocarbons at High Pressures and Temperatures.

Although long-chain aliphatic hydrocarbons are documented in meteorites, their origin is poorly understood. A key question is whether they are pristine or a byproduct of terrestrial alteration? To understand if these long-chain hydrocarbons are indigenous, it will be important to explore their thermodynamic and mechanical stability at conditions experienced by extraterrestrial objects during atmospheric entry and passage. Extreme pressures and temperatures experienced by meteorites are likely to alter the molecular organization of these long-chain hydrocarbons. These structural changes associated with extreme conditions are often documented via laboratory-based Raman spectroscopic measurements. So far, Raman spectroscopic measurements have investigated the effect of static compression on the aliphatic hydrocarbons. The effect of temperature on the structural changes remains poorly explored. To bridge this gap, in this study, we have explored the behavior of two aliphatic hydrocarbons at simultaneously high pressures and temperatures. We have used a resistively heated diamond anvil cell. On compression to moderate pressures, the appearance of new vibrational modes in the low-energy region confirms prior studies and is related to the bending of the linear chains. Upon heating to ∼220 °C, we note that the new low-energy mode undergoes softening. The mode softening is likely due to the combination of unbending of the alkane chain and mode anharmonicity.

Tracking slab surface temperatures with electrical conductivity of glaucophane.

Slab surface temperature is one of the key parameters that incur first-order changes in subduction dynamics. However, the current thermal models are based on empirical thermal parameters and do not accurately capture the complex pressure-temperature paths of the subducting slab, prompting significant uncertainties on slab temperature estimations. In this study, we investigate whether the dehydration-melting of glaucophane can be used to benchmark the temperature in the slab. We observe that dehydration and melting of glaucophane occur at relatively low temperatures compared to the principal hydrous phases in the slab and produce highly conductive Na-rich melt. The electrical properties of glaucophane and its dehydration products are notably different from the hydrous minerals and silicate melts. Hence, we conclude that the thermodynamic instability of glaucophane in the slab provides a unique petrological criterion for tracking temperature in the present-day subduction systems through magnetotelluric profiles.

Insights on the deep carbon cycle from the electrical conductivity of carbon-bearing aqueous fluids.

The dehydration and decarbonation in the subducting slab are intricately related and the knowledge of the physical properties of the resulting C-H-O fluid is crucial to interpret the petrological, geochemical, and geophysical processes associated with subduction zones. In this study, we investigate the C-H-O fluid released during the progressive devolatilization of carbonate-bearing serpentine-polymorph chrysotile, with in situ electrical conductivity measurements at high pressures and temperatures. The C-H-O fluid produced by carbonated chrysotile exhibits high electrical conductivity compared to carbon-free aqueous fluids and can be an excellent indicator of the migration of carbon in subduction zones. The crystallization of diamond and graphite indicates that the oxidized C-H-O fluids are responsible for the recycling of carbon in the wedge mantle. The carbonate and chrysotile bearing assemblages stabilize dolomite during the devolatilization process. This unique dolomite forming mechanism in chrysotile in subduction slabs may facilitate the transport of carbon into the deep mantle.

High-pressure elastic properties of dolomite melt supporting carbonate-induced melting in deep upper mantle.

Deeply subducted carbonates likely cause low-degree melting of the upper mantle and thus play an important role in the deep carbon cycle. However, direct seismic detection of carbonate-induced partial melts in the Earth's interior is hindered by our poor knowledge on the elastic properties of carbonate melts. Here we report the first experimentally determined sound velocity and density data on dolomite melt up to 5.9 GPa and 2046 K by in-situ ultrasonic and sink-float techniques, respectively, as well as first-principles molecular dynamics simulations of dolomite melt up to 16 GPa and 3000 K. Using our new elasticity data, the calculated VP/VS ratio of the deep upper mantle (∼180-330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated VP/VS ratio of the upper mantle from global seismic observations, supporting the pervasive presence of a low-degree carbonate-rich partial melt (∼0.05%) that is consistent with the volatile-induced or redox-regulated initial melting in the upper mantle as argued by petrologic studies. This carbonate-rich partial melt region implies a global average carbon (C) concentration of 80-140 ppm. by weight in the deep upper mantle source region, consistent with the mantle carbon content determined from geochemical studies.

Single crystal elasticity of natural topaz at high-temperatures.

Topaz is an aluminosilicate mineral phase stable in the hydrothermally altered pegmatitic rocks and also in subducted sedimentary lithologies. In nature, topaz often exhibits solid solution between fluorine and hydrous end members. We investigated elasticity of naturally occurring single crystal topaz (Al2SiO4F1.42(OH)0.58) using Resonant Ultrasound Spectroscopy. We also explored the temperature dependence of the full elastic constant tensor. We find that among the various minerals stable in the Al2O3-SiO2-H2O ternary system, topaz exhibits moderate elastic anisotropy. As a function of temperature, the sound velocity of topaz decreases with [Formula: see text] and [Formula: see text] being -3.10 and -2.30 × 10-4 km/s/K. The elasticity and sound velocity of topaz also vary as a function of OH and F content. The effect of composition ([Formula: see text]) on the velocity is equally important as that of the effect of temperature. We also note that the Debye temperature ([Formula: see text]) of topaz at room temperature condition is 910 K and decreases at higher temperature. The Debye temperature shows positive correlation with density of the mineral phases in the Al2O3-SiO2-H2O ternary system.

Carbon-bearing silicate melt at deep mantle conditions.

Knowledge about the incorporation and role of carbon in silicate magmas is crucial for our understanding of the deep mantle processes. CO2 bearing silicate melting and its relevance in the upper mantle regime have been extensively explored. Here we report first-principles molecular dynamics simulations of MgSiO3 melt containing carbon in three distinct oxidation states - CO2, CO, and C at conditions relevant for the whole mantle. Our results show that at low pressures up to 15 GPa, the carbon dioxide speciation is dominated by molecular form and carbonate ions. At higher pressures, the dominant species are silicon-polyhedral bound carbonates, tetrahedral coordination, and polymerized di-carbonates. Our results also indicate that CO2 component remains soluble in the melt at high pressures and the solution is nearly ideal. However, the elemental carbon and CO components show clustering of carbon atoms in the melt at high pressures, hinting towards possible exsolution of carbon from silicate melt at reduced oxygen contents. Although carbon lowers the melt density, the effect is modest at high pressures. Hence, it is likely that silicate melt above and below the mantle transition zone, and atop the core-mantle boundary could efficiently sequester significant amounts of carbon without being gravitationally unstable.

High-pressure phase of brucite stable at Earth's mantle transition zone and lower mantle conditions.

We investigate the high-pressure phase diagram of the hydrous mineral brucite, Mg(OH)2, using structure search algorithms and ab initio simulations. We predict a high-pressure phase stable at pressure and temperature conditions found in cold subducting slabs in Earth's mantle transition zone and lower mantle. This prediction implies that brucite can play a much more important role in water transport and storage in Earth's interior than hitherto thought. The predicted high-pressure phase, stable in calculations between 20 and 35 GPa and up to 800 K, features MgO6 octahedral units arranged in the anatase-TiO2 structure. Our findings suggest that brucite will transform from a layered to a compact 3D network structure before eventual decomposition into periclase and ice. We show that the high-pressure phase has unique spectroscopic fingerprints that should allow for straightforward detection in experiments. The phase also has distinct elastic properties that might make its direct detection in the deep Earth possible with geophysical methods.

Pressure induced elastic softening in framework aluminosilicate- albite (NaAlSi3O8).

Albite (NaAlSi3O8) is an aluminosilicate mineral. Its crystal structure consists of 3-D framework of Al and Si tetrahedral units. We have used Density Functional Theory to investigate the high-pressure behavior of the crystal structure and how it affects the elasticity of albite. Our results indicate elastic softening between 6-8 GPa. This is observed in all the individual elastic stiffness components. Our analysis indicates that the softening is due to the response of the three-dimensional tetrahedral framework, in particular by the pressure dependent changes in the tetrahedral tilts. At pressure <6 GPa, the PAW-GGA can be described by a Birch-Murnaghan equation of state with = 687.4 Å3, = 51.7 GPa, and = 4.7. The shear modulus and its pressure derivative are = 33.7 GPa, and = 2.9. At 1 bar, the azimuthal compressional and shear wave anisotropy = 42.8%, and = 50.1%. We also investigate the densification of albite to a mixture of jadeite and quartz. The transformation is likely to cause a discontinuity in density, compressional, and shear wave velocity across the crust and mantle. This could partially account for the Mohorovicic discontinuity in thickened continental crustal regions.

Dehydration of chlorite explains anomalously high electrical conductivity in the mantle wedges.

Mantle wedge regions in subduction zone settings show anomalously high electrical conductivity (~1 S/m) that has often been attributed to the presence of aqueous fluids released by slab dehydration. Laboratory-based measurements of the electrical conductivity of hydrous phases and aqueous fluids are significantly lower and cannot readily explain the geophysically observed anomalously high electrical conductivity. The released aqueous fluid also rehydrates the mantle wedge and stabilizes a suite of hydrous phases, including serpentine and chlorite. In this present study, we have measured the electrical conductivity of a natural chlorite at pressures and temperatures relevant for the subduction zone setting. In our experiment, we observe two distinct conductivity enhancements when chlorite is heated to temperatures beyond its thermodynamic stability field. The initial increase in electrical conductivity to ~3 × 10(-3) S/m can be attributed to chlorite dehydration and the release of aqueous fluids. This is followed by a unique, subsequent enhancement of electrical conductivity of up to 7 × 10(-1) S/m. This is related to the growth of an interconnected network of a highly conductive and chemically impure magnetite mineral phase. Thus, the dehydration of chlorite and associated processes are likely to be crucial in explaining the anomalously high electrical conductivity observed in mantle wedges. Chlorite dehydration in the mantle wedge provides an additional source of aqueous fluid above the slab and could also be responsible for the fixed depth (120 ± 40 km) of melting at the top of the subducting slab beneath the subduction-related volcanic arc front.

Crystal structure, equation of state, and elasticity of phase H (MgSiO4H2) at Earth's lower mantle pressures.

Dense hydrous magnesium silicate (DHMS) phases play a crucial role in transporting water in to the Earth's interior. A newly discovered DHMS, phase H (MgSiO4H2), is stable at Earth's lower mantle, i.e., at pressures greater than 30 GPa. Here we report the crystal structure and elasticity of phase H and its evolution upon compression. Using first principles simulations, we have explored the relative energetics of the candidate crystal structures with ordered and disordered configurations of magnesium and silicon atoms in the octahedral sites. At conditions relevant to Earth's lower mantle, it is likely that phase H is able to incorporate a significant amount of aluminum, which may enhance the thermodynamic stability of phase H. The sound wave velocities of phase H are ~2-4% smaller than those of isostructural δ-AlOOH. The shear wave impedance contrast due to the transformation of phase D to a mixture of phase H and stishovite at pressures relevant to the upper part of the lower mantle could partly explain the geophysical observations. The calculated elastic wave velocities and anisotropies indicate that phase H can be a source of significant seismic anisotropy in the lower mantle.

Hydrous silicate melt at high pressure.

The structure and physical properties of hydrous silicate melts and the solubility of water in melts over most of the pressure regime of Earth's mantle (up to 136 GPa) remain unknown. At low pressure (up to a few gigapascals) the solubility of water increases rapidly with increasing pressure, and water has a large influence on the solidus temperature, density, viscosity and electrical conductivity. Here we report the results of first-principles molecular dynamics simulations of hydrous MgSiO3 melt. These show that pressure has a profound influence on speciation of the water component, which changes from being dominated by hydroxyls and water molecules at low pressure to extended structures at high pressure. We link this change in structure to our finding that the water-silicate system becomes increasingly ideal at high pressure: we find complete miscibility of water and silicate melt throughout almost the entire mantle pressure regime. On the basis of our results, we argue that a buoyantly stable melt at the base of the upper mantle would contain approximately 3 wt% water and have an electrical conductivity of 18 S m(-1), and should therefore be detectable by means of electromagnetic sounding.

The effect of water on the electrical conductivity of olivine.

It is well known that water (as a source of hydrogen) affects the physical and chemical properties of minerals--for example, plastic deformation and melting temperature--and accordingly plays an important role in the dynamics and geochemical evolution of the Earth. Estimating the water content of the Earth's mantle by direct sampling provides only a limited data set from shallow regions (<200 km depth). Geophysical observations such as electrical conductivity are considered to be sensitive to water content, but there has been no experimental study to determine the effect of water on the electrical conductivity of olivine, the most abundant mineral in the Earth's mantle. Here we report a laboratory study of the dependence of the electrical conductivity of olivine aggregates on water content at high temperature and pressure. The electrical conductivity of synthetic polycrystalline olivine was determined from a.c. impedance measurements at a pressure of 4 GPa for a temperature range of 873-1,273 K for water contents of 0.01-0.08 wt%. The results show that the electrical conductivity is strongly dependent on water content but depends only modestly on temperature. The water content dependence of conductivity is best explained by a model in which electrical conduction is due to the motion of free protons. A comparison of the laboratory data with geophysical observations suggests that the typical oceanic asthenosphere contains approximately 10(-2) wt% water, whereas the water content in the continental upper mantle is less than approximately 10(-3) wt%.