Electrical conductivity of melts: implications for conductivity anomalies in the Earth's mantle

Magmatic liquids, including silicate and carbonate melts, are principal agents of mass and heat transfer in the Earth and terrestrial planets, and they play a crucial role in various geodynamic processes and in Earth's evolution. Electrical conductivity data of these melts elucidate the cause of electrical anomalies in Earth's interior and shed light on the melt structure. With the improvement in high-pressure experimental techniques and theoretical simulations, major progress has been made on this front in the past several decades. This review aims to summarize recent advances in experimental and theoretical studies on the electrical conductivity of silicate and carbonate melts of different compositions and volatile contents under high temperature and pressure. The electrical conductivity of silicate melts depends strongly on temperature, pressure, water content and the ratio of non-bridging oxygens to tetrahedral cations (NBO/T). By contrast, the electrical conductivity of carbonate melts exhibits a weak dependence on temperature and pressure due to their fully depolymerized structure. The electrical conductivity of carbonate melts is higher than that of silicate melts by at least two orders of magnitude. Water can increase electrical conductivity significantly and reduce the activation energy of silicate melts. Conversely, this effect is weak for carbonate melts. In addition, the replacement of alkali-earth elements (Ca²⁺ or Mg²⁺) with alkali elements causes a significant decrease in the electrical conductivity of carbonate melts. A distinct compensation trend is revealed for the electrical conductivity of silicate and carbonate melts under anhydrous and hydrous conditions. Several important applications of laboratory-based melt conductivity are introduced in order to understand the origin of high-conductivity anomalies in the Earth's mantle. Perspectives for future studies are also provided.

Behavior and properties of water in silicate melts under deep mantle conditions

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-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 almost 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.

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

Petrologic studies suggest that carbonate-rich melts are present in the Earth’s upper mantle and play an important role in the deep carbon cycle. However, seismic detection of these melts is difficult due to the lack of understanding of the elastic properties of carbonate melts. Here we determined the sound velocity and density of dolomite melt at upper mantle conditions using high-pressure experiments and theoretical simulations. The calculated velocities of carbonate melt-bearing mantle using these new elasticity data were compared with global seismic observations. The result suggests that ∼0.05% carbonate-rich melt can be pervasively present in the deep upper mantle, implying a global average mantle carbon concentration of 80-140 ppm.

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 V_P/V_S ratio of the deep upper mantle (∼180–330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated V_P/V_S 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.

Transport properties of carbonated silicate melt at high pressure

Carbon dioxide, generally considered as the second most abundant volatile component in silicate magmas, is expected to significantly influence various melt properties. In particular, our knowledge about its dynamical effects is lacking over most of Earth's mantle pressure regime. Here, we report the first-principles molecular dynamics results on the transport properties of carbonated MgSiO3 liquid under conditions of mantle relevance. They show that dissolved CO2 systematically enhances the diffusion rates of all elements and lowers the melt viscosity on average by factors of 1.5 to 3 over the pressure range considered. It is remarkable that CO2 has very little or no influence on the electrical conductivity of the silicate melt under most conditions. Simulations also predict anomalous dynamical behavior, increasing diffusivity and conductivity and decreasing viscosity with compression in the low-pressure regime. This anomaly and the concomitant increase of pressure and temperature with depth together make these transport coefficients vary modestly over extended portions of the mantle regime. It is possible that the melt electrical conductivity under conditions corresponding to the 410- and 660-km seismic discontinuities is at a detectable level by electromagnetic sounding observation. In addition, the low melt viscosity values of 0.2 to 0.5 Pa⋅s at these depths and near the core-mantle boundary may imply high mobility of possible melts in these regions.

Heterogeneity in mantle carbon content from CO2-undersaturated basalts

The amount of carbon present in Earth's mantle affects the dynamics of melting, volcanic eruption style and the evolution of Earth's atmosphere via planetary outgassing. Mantle carbon concentrations are difficult to quantify because most magmas are strongly degassed upon eruption. Here we report undegassed carbon concentrations from a new set of olivine-hosted melt inclusions from the Mid-Atlantic Ridge. We use the correlations of CO₂ with trace elements to define an average carbon abundance for the upper mantle. Our results indicate that the upper mantle carbon content is highly heterogeneous, varying by almost two orders of magnitude globally, with the potential to produce large geographic variations in melt fraction below the volatile-free solidus. Such heterogeneity will manifest as variations in the depths at which melt becomes interconnected and detectable, the CO₂ fluxes at mid-ocean ridges, the depth of the lithosphere-asthenosphere boundary, and mantle conductivity.

Electrical conductivity during incipient melting in the oceanic low-velocity zone

A low viscosity layer in the upper mantle, the Asthenosphere, is a requirement for plate tectonics¹. The seismic low velocities and the high electrical conductivities of the Asthenosphere are attributed either to sub-solidus water-related defects in olivine minerals^2-4 or to a few volume percents of partial melt^5-8 but these two interpretations have shortcomings: (1) The amount of H₂O stored in olivine is not expected to be higher than 50 ppm due to partitioning with other mantle phases⁹, including pargasite amphibole at moderate temperatures¹⁰, and partial melting at high temperatures⁹; (2) elevated melt volume fractions are impeded by the too cold temperatures prevailing in the Asthenosphere and by the high melt mobility that can lead to gravitational segregation^11,12. Here we determined the electrical conductivity of CO₂-H₂O-rich melts, typically produced at the onset of mantle melting. Electrical conductivity modestly increases with moderate amounts of H₂O and CO₂ but it dramatically increases as CO₂ content exceeds 6 wt% in the melt. Incipient melts, long-expected to prevail in the asthenosphere^10,13-15, can therefore trigger its high electrical conductivities. Considering depleted and enriched mantle abundances in H₂O and CO₂ and their effect on the petrology of incipient melting, we calculated conductivity profiles across the Asthenosphere for various plate ages. Several electrical discontinuities are predicted and match geophysical observations in a consistent petrological and geochemical framework. In moderately aged plates (>5Ma), incipient melts most likely trigger both the seismic low velocities and the high electrical conductivities in the upper part of the asthenosphere, whereas for young plates⁴, where seamount volcanism occurs⁶, higher degree of melting is expected.

Global electromagnetic induction constraints on transition-zone water content variations

Small amounts of water can significantly affect the physical properties of mantle materials, including lowering of the solidus, and reducing effective viscosity and seismic velocity. The amount and distribution of water within the mantle thus has profound implications for the dynamics and geochemical evolution of the Earth. Electrical conductivity is also highly sensitive to the presence of hydrogen in mantle minerals. The mantle transition zone minerals wadsleyite and ringwoodite in particular have high water solubility, and recent high pressure experiments show that the electrical conductivity of these minerals is very sensitive to water content. Thus estimates of the electrical conductivity of the mantle transition zone derived from electromagnetic induction studies have the potential to constrain the water content of this region. Here we invert long period geomagnetic response functions to derive a global-scale three-dimensional model of electrical conductivity variations in the Earth's mantle, revealing variations in the electrical conductivity of the transition zone of approximately one order of magnitude. Conductivities are high in cold, seismically fast, areas where slabs have subducted into or through the transition zone. Significant variations in water content throughout the transition zone provide a plausible explanation for the observed patterns. Our results support the view that at least some of the water in the transition zone has been carried into that region by cold subducting slabs.

Dry mantle transition zone inferred from the conductivity of wadsleyite and ringwoodite

The Earth's mantle transition zone could potentially store a large amount of water, as the minerals wadsleyite and ringwoodite incorporate a significant amount of water in their crystal structure. The water content in the transition zone can be estimated from the electrical conductivities of hydrous wadsleyite and ringwoodite, although such estimates depend on accurate knowledge of the two conduction mechanisms in these minerals (small polaron and proton conductions), which early studies have failed to distinguish between. Here we report the electrical conductivity of these two minerals obtained by high-pressure multi-anvil experiments. We found that the small polaron conductions of these minerals are substantially lower than previously estimated. The contributions of proton conduction are small at temperatures corresponding to the mantle transition zone and the conductivity of wadsleyite is considerably lower than that of ringwoodite for both mechanisms. The dry model mantle shows considerable conductivity jumps associated with the olivine-wadsleyite, wadsleyite-ringwoodite and post-spinel transitions. Such a dry model explains well the currently available conductivity-depth profiles obtained from geoelectromagnetic studies. We therefore conclude that there is no need to introduce a significant amount of water in the mantle transition to satisfy electrical conductivity constraints.

Electromagnetic detection of a 410-km-deep melt layer in the southwestern United States

A deep-seated melt or fluid layer on top of the 410-km-deep seismic discontinuity in Earth's upper mantle, as proposed in the transition-zone 'water filter' hypothesis, may have significant bearing on mantle dynamics and chemical differentiation. The geophysical detection of such a layer has, however, proved difficult. Magnetotelluric and geomagnetic depth sounding are geophysical methods sensitive to mantle melt. Here we use these methods to search for a distinct structure near 410-km depth. We calculate one-dimensional forward models of the response of electrical conductivity depth profiles, based on mineral physics studies of the effect of incorporating hydrogen in upper-mantle and transition-zone minerals. These models indicate that a melt layer at 410-km depth is consistent with regional magnetotelluric and geomagnetic depth sounding data from the southwestern United States (Tucson). The 410-km-deep melt layer in this model has a conductance of 3.0 x 10(4) S and an estimated thickness of 5-30 km. This is the only regional data set that we have examined for which such a melt layer structure was found, consistent with regional seismic studies. We infer that the hypothesized transition-zone water filter occurs regionally, but that such a layer is unlikely to be a global feature.

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%.

Hydrous olivine unable to account for conductivity anomaly at the top of the asthenosphere

The oceanic asthenosphere is observed to have high electrical conductivity, which is highly anisotropic in some locations. In the directions parallel and normal to the plate motion, the conductivity is of the order of 10(-1) and 10(-2) S m(-1), respectively, which cannot be explained by the conductivity of anhydrous olivine. But because hydrogen can be incorporated in olivine at mantle pressures, this observation has been attributed to olivine hydration, which might cause anisotropically high conductivity by proton migration. To examine this hypothesis, here we report the effect of water on electrical conductivity and its anisotropy for hydrogen-doped and undoped olivine at 500-1,500 K and 3 GPa. The hydrous olivine has much higher conductivity and lower activation energy than anhydrous olivine in the investigated temperature range. Nevertheless, extrapolation of the experimental results suggests that conductivity of hydrous olivine at the top of the asthenosphere should be nearly isotropic and only of the order of 10(-2) S m(-1). Our data indicate that the hydration of olivine cannot account for the geophysical observations, which instead may be explained by the presence of partial melt elongated in the direction of plate motion.

Geophysical evidence from the MELT area for compositional controls on oceanic plates

Magnetotelluric and seismic data, collected during the MELT experiment at the southern East Pacific Rise, constrain the distribution of melt beneath this mid-ocean-ridge spreading centre and also the evolution of the oceanic lithosphere during its early cooling history. Here we focus on structures imaged at distances approximately 100 to 350 km east of the ridge crest, corresponding to seafloor ages of approximately 1.3 to 4.5 million years (Myr), where the seismic and electrical conductivity structure is nearly constant and independent of age. Beginning at a depth of about 60 km, we image a large increase in electrical conductivity and a change from isotropic to transversely anisotropic electrical structure, with higher conductivity in the direction of fast propagation for seismic waves. Conductive cooling models predict structure that increases in depth with age, extending to about 30 km at 4.5 Myr ago. We infer, however, that the structure of young oceanic plates is instead controlled by a decrease in water content above a depth of 60 km induced by the melting process beneath the spreading centre.

Asymmetric Electrical Structure in the Mantle Beneath the East Pacific Rise at 17 degrees S

The magnetotelluric component of the Mantle Electromagnetic and Tomography (MELT) Experiment measured the electrical resistivity structure of the mantle beneath the fast-spreading southern East Pacific Rise (EPR). The data reveal an asymmetric resistivity structure, with lower resistivity to the west of the ridge. The uppermost 100 kilometers of mantle immediately to the east of the ridge is consistent with a dry olivine resistivity structure indicating a mantle depleted of melt and volatiles. Mantle resistivities to the west of the ridge are consistent with a low-melt fraction (about 1 to 2 percent interconnected melt) distributed over a broad region and extending to depths of about 150 kilometers. The asymmetry in resistivity structure may be the result of asymmetric spreading rates and a westward migration of the ridge axis and suggests distinct styles of melt formation and delivery in the mantle beneath the two plates.