skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Anomalous Behavior of Viscosity and Electrical Conductivity of MgSiO 3 Melt at Mantle Conditions
Abstract Silicate melts have served as transport agents in the chemical and thermal evolution of Earth. Molecular dynamics simulations based on a deep neural network potential trained byab initiodata show that the viscosity of MgSiO3melt decreases with increasing pressure at low pressures (up to ∼6 GPa) before it starts to increase with further compression. The melt electrical conductivity also behaves anomalously; first increasing and then decreasing with pressure. The melt accumulation implied by the viscosity turnover at ∼23 GPa along mantle liquidus offers an explanation for the low‐velocity zone at the 660‐km discontinuity. The increase in electrical conductivity up to ∼50 GPa may contribute to the steep rise of Earth's electrical conductivity profiles derived from magnetotelluric observations. Our results also suggest that small fraction of melts could give rise to detectable bulk conductivity in deeper parts of the mantle.  more » « less
Award ID(s):
2001074
PAR ID:
10359970
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
48
Issue:
13
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Water (H2O) as one of the most abundant fluids present in Earth plays crucial role in the generation and transport of magmas in the interior. Though hydrous silicate melts have been studied extensively, the experimental data are confined to relatively low pressures and the computational results are still rare. Moreover, these studies imply large differences in the way water influences the physical properties of silicate magmas, such as density and electrical conductivity. Here, we investigate the equation of state, speciation, and transport properties of water dissolved in Mg1-xFexSiO3 and Mg2(1-x)Fe2xSiO4 melts (for x = 0 and 0.25) as well as in its bulk (pure) fluid state over the entire mantle pressure regime at 2000 to 4000 K using first-principles molecular dynamics. The simulation results allow us to constrain the partial molar volume of the water component in melts along with the molar volume of pure water. The predicted volume of silicate melt+water solution is negative at low pressures and becomes zero above 15 GPa. Consequently, the hydrous component tends to lower the melt density to similar extent over much of the mantle pressure regime irrespective of composition. Our results also show that hydrogen diffuses fast in silicate melts and enhances the melt electrical conductivity in a way that differs from electrical conduction in the bulk water. The speciation of the water component varies considerably from the bulk water structure as well. Water is dissolved in melts mostly as hydroxyls at low pressure and as -O-H-O-, -O-H-O-H- and other extended species with increasing pressure. On the other hand, the pure water behaves as a molecular fluid below 15 GPa, gradually becoming a dissociated fluid with further compression. On the basis of modeled density and conductivity results, we suggest that partial melts containing a few percent of water may be gravitationally trapped both above and below the upper mantle-transition region. Moreover, such hydrous melts can give rise to detectable electrical conductance by means of electromagnetic sounding observations. 
    more » « less
  2. Abstract The continental crust is produced by the solidification of aluminosilicate‐rich magmas which are sourced from deep below the surface. Migration of the magma depends on the density (ρ) contrast to source rocks and the melt viscosity (η). At the surface, these silica‐rich melts are typically sluggish due to highη > 1,000 Pa s. Yet at their source regions, the melt properties are complexly influenced by pressure (P), temperature (T), and water contents (). In this study, we examined the combinedP‐T‐ effects on the behavior of melts with an albite stoichiometry (NaAlSi3O8). We usedfirst‐principlesmolecular dynamics simulations to examine anhydrous (0 wt % H2O) and hydrous (5 wt % H2O) melts. To constrain thePandTeffects, we exploredP ≤ 25 GPa across several isotherms between 2500 and 4000 K. The melts show anomalousP‐ρrelationships at lowP ∼ 0 GPa and highT ≥ 2500 K, consistent with vaporization. At lithospheric conditions, meltρincreases with compression and is well described by a finite‐strain formalism. Water lowers the melt density (ρhydrous < ρanhydrous) but increases the compressibility, that is, 1/Khydrous>1/KanhydrousorKhydrous < Kanhydrous. We also find that the meltηdecreases with pressure and then increases with further compression. Water decreases the viscosity (ηhydrous < ηanhydrous) by depolymerizing the melt structure. The ionic self‐diffusivities are increased by the presence of water. The decreasedρandηby H2O increase the mobility of magma at crustal conditions, which could explain the rapid eruption and migration timescales for rhyolitic magmas as observed in the Chaitén volcano in Chile. 
    more » « less
  3. In this study, we use f irst-principles molecular dynamics simulations to explore the behavior of anhydrous aluminosilicate melt with a stoichiometry of NaAlSi2O6 up to pressures of ∼30 GPa and temperatures between 2500 and 4000 K. We also examine the effect of water (∼4 wt % H2O) on the equation of state and transport properties of the aluminosilicate melt and relate them to atomistic scale changes in the melt structure. Our results show that water reduces the density and bulk modulus of the anhydrous melt. However, the pressure derivative of the bulk modulus of the hydrous melt is larger than that of the anhydrous melt. The pressure dependence of the transport property exhibits an anomalous behavior. At a pressure of ∼12 GPa, anhydrous aluminosilicate melts exhibit maxima in diffusion and minima in viscosity. Dissolved water in melts also affects both diffusion and viscosity. In hydrous aluminosilicate melts, the maxima in diffusion and the minima in viscosity occur at ∼14 GPa. The anomalous behavior of transport properties is related to the pressure-induced changes in the melt structure. At shallower depths, i.e., up to 100 km, relevant for subduction zone settings, the lower density compounded by the lower viscosity of hydrous aluminosilicate melts is likely to provide buoyancy for upward migration. At greater depths of ∼180−200 km, greater compressibility of the hydrous aluminosilicate melts together with the minimum viscosity could hinder magma migration and may explain the presence of a partial melt layer at the lithosphere−asthenosphere boundary. 
    more » « less
  4. The continental crust is rich in aluminosilicates and formed by the crystallization of arc magmas. However, the magma produced at sub-arc depths is often silica-poor. The chemical evolution of sub-arc magma from silica-poor to aluminosilicate-rich is perplexing. Magnetotelluric (MT) observations in subduction zones and complementary laboratory-based constraints of electrical conductivity (σ) are crucial to understanding this chemical evolution. The σ of a magma is sensitive to pressure (P), temperature (T), and chemistry (X). To date, laboratory-based measurements on the σ of silicate melts have helped to interpret MT observations at P ≤ 2 GPa. Yet, the melting in subduction zones could occur deeper, at P ≤ 6−7 GPa. The σ of melt at such pressures is poorly constrained. To address this, we performed experiments at P ≤ 6 GPa to examine the σ of basaltic to andesitic melts, which are common in subduction zones. We constrained the effects of silica, alumina, alkali, alkaline, and water (H2O) contents on the σ of melt. The activation volume of σ increases with silica contents. Hence, the σ of basaltic melt is overall greater than that of an andesitic counterpart. The σ of basaltic magma is also less sensitive to P than andesitic magma. Water lowers the activation energy and enhances σ for all melt compositions. Our results help constrain how the electrical properties of a magma change with an evolving composition in a subduction zone. 
    more » « less
  5. Abstract We present phase‐equilibria experiments of a K‐bearing, depleted peridotite (Mg# 92) fluxed with a mixed CO2‐H2O fluid (0.5 wt.% CO2and 0.94 wt.% H2O in the bulk) to gain insight into the stability of volatile‐bearing partial melts versus volatile‐bearing mineral phases in a depleted peridotite system. Experiments were performed at 850–1150 °C and 2–4 GPa using a piston‐cylinder and a multianvil apparatus. Olivine, orthopyroxene, clinopyroxene, and spinel/garnet are present at all experimental conditions. Textural confirmation of partial melt is made at temperatures as low as 1000 °C at 2 GPa, 950 °C at 3 GPa, and 1000 °C at 4 GPa marking the onset of melting at 900–1000 °C at 2 GPa, 850–950 °C at 3 GPa, and 950–1000 °C at 3 GPa. Phlogopite and magnesite breakdown at 900–1000 °C at 2 GPa, 950–1000 °C at 3 GPa, and 1000–1050 °C at 4 GPa. Comparison with previously published experiments in depleted peridotite system with identical CO2‐H2O content introduced via a silicic melt show that introduction of CO2‐H2O as fluid lowers the temperature of phlogopite breakdown by 150–200 °C at 2–4 GPa and stabilizes partial melts at lower temperatures. Our study thus, shows that the volatile‐bearing phase present in the cratonic mantle is controlled by bulk composition and is affected by the process of volatile addition during craton formation in a subduction zone. In addition, volatile introduction via melt versus aqueous fluid, leads to different proportion of anhydrous phases such as olivine and orthopyroxene. Considering the agent of metasomatism is thus critical to evaluate how the bulk composition of depleted peridotite is modified, leading to potential stability of volatile‐bearing phases as the cause of anomalously low shear wave velocity in mantle domains such as mid lithospheric discontinuities beneath continents. 
    more » « less