# Tracking slab surface temperatures with electrical conductivity of glaucophane

Geeth Manthilake, Ye Peng, Kenneth T. Koga, Mainak Mookherjee

► To cite this version:

Geeth Manthilake, Ye Peng, Kenneth T. Koga, Mainak Mookherjee. Tracking slab surface temperatures with electrical conductivity of glaucophane. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-97317-0. hal-03349098

HAL Id: hal-03349098

https://uca.hal.science/hal-03349098

Submitted on 20 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Tracking slab surface temperatures with electrical conductivity of glaucophane

Geeth Manthilake<fnref n="1" />✉, Ye Peng², Kenneth T. Koga<fnref n="1" /> & Mainak Mookherjee²

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.

Upon subduction, hydrous mineral phases in a subducting slab reach their limits of thermodynamic stability and produce a continuous flux of mobile phases, i.e. aqueous fluids, silicate melts, or supercritical fluids, into the overlying mantle through multi-phase, multi-variant dehydration reactions<fnref n="1" />. The migration of mobile phases into the overlying mantle metasomatizes the mantle wedge²,³, triggers partial melting in the peridotite mantle⁴, and plays a crucial role in recycling volatiles⁵⁻⁸, and trace elements in subduction zones⁹.

Among numerous prograde metamorphic dehydration reactions, the transition from blueschist facies to eclogite facies has been considered the most significant, as it characterizes the transformation of altered and less dense oceanic crust into an anhydrous and dense slab¹⁰. A consensus view of the geoscience community is that the blueschist-eclogite transition has a depth range that locates beneath arc volcanoes¹¹ and is considered as the principal driving force behind the fluid-assisted melting in the mantle wedge¹⁰,¹²,¹³.

The variations in slab temperature, composition, and subduction dynamics make each subduction zone unique and the differences are often reflected in the seismic velocity profiles, arc magmatism, and seismicity¹⁰,¹⁴. Slab temperature is one of the governing parameters that define the stability of hydrous minerals, and therefore magmatism of arc volcanoes. The subduction zone temperature profiles are estimated primarily based on numerical models, which use finite element analyses with a range of input parameters¹⁰,¹⁵,¹⁶. Geochemical constraints such as the H₂O/Ce ratio in slab fluid compositions have also been used to predict slab surface temperature¹⁷,¹⁸. The H₂O/Ce thermometer primarily depends on the chemical analyses of melt inclusion found in arc magmas. However, it has been observed that the degassing alters the H₂O/Ce ratio by losing the primitive H₂O even in inclusions¹⁹. While the existing models provide insights into the complex interplay of geophysical and geochemical parameters of the thermal structures of subduction zones, the lack of temperature-fixed points makes it harder to evaluate the accuracy of subduction zone temperature profiles. As a result, some thermal models report a geothermal gradient (dT/dz) < 5 K/km, which is unlikely for subduction zones²⁰. This suggests that although the present parametrizations in thermal models are internally consistent, they may not adequately estimate the realistic subduction zone conditions²¹. A temperature-fixed point in the subducting slab, if it can be constrained with an independent observation, is crucial for constructing geodynamic models and improving the understanding of mantle flow within subduction zones.

Glaucophane (Na₂(Mg₃Al₂)Si₈O₂₂(OH)₂) and lawsonite (CaAl₂Si₂O₇(OH)₂·H₂O) are the principal hydrous mineral phases associated with the blueschist facies lithology<fnref n="1" />,²,²²,²³. The presence of glaucophane and lawsonite is often associated with anomalous geophysical observations including a relatively low-velocity layer²⁴,²⁵ and enhanced electrical conductivity²⁶,²⁷ at the top of subducting slabs. The stability of lawsonite at high pressures

<footnote n="1">¹Laboratoire Magmas et Volcans CNRS, IRD, OPGC, Université Clermont Auvergne, 63000 Clermont-Ferrand, France. ²Earth Materials Laboratory, Department of Earth, Ocean and Atmospheric Sciences, Florida State University, Tallahassee, FL 32306, USA. ✉email: geeth.manthilake@uca.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The electrical conductivity of glaucophane as a function of reciprocal temperature. (a) At 1.5 GPa, electrical conductivity is characterized by three stages, separated by discontinuous increases of conductivity at ~825 K and 1000 K. (b) Electrical conductivity of glaucophane at 6 GPa. The dehydration of the sample occurs at 895 K. The activation enthalpies corresponding to the dominant conduction mechanism at each stage are shown next to the fitting lines. The experimental uncertainty can arise from the estimation of temperature, geometrical parameters. Fitting errors are estimated to be less than 5% and are smaller than the symbol size.</figcaption></figure>

(>9 GPa) and high temperature up to 1300 K means that the lawsonite potentially remains stable to depths of more than 200 km²²,²³. In contrast, glaucophane is thermodynamically stable at lower P–T conditions and its dehydration marks the onset of the transition of hydrous oceanic crust to the dense anhydrous eclogite¹⁰. Notably, it has been observed that amphibole stability can be enhanced beyond 5 GPa by the substitution of fluorine for amphibole hydroxyl sites. Hence, depending on F contents, the blueschist-eclogite transition may shift to depths beyond 75 km²⁸,²⁹.

It is well established that the electrical conductivity of fluids or melts can be enhanced by the highly mobile Na⁺ as a charge carrier³⁰. Because glaucophane is one of the principal hosts of Na in subducting slabs, we focus on understanding whether the dehydration-melting of glaucophane can generate distinct geophysically detectable electrical anomalies, standing out from the hydrous minerals, fluid, and melts in subduction zones. We investigate the electrical conductivity of glaucophane in situ under high pressure and high temperature up to 6 GPa and 1258 K. We show that the dehydration melting of glaucophane produces highly conductive Na-rich fluids and melt that can serve as geophysical indicators for accurately benchmarking the slab temperatures.

## Results

The electrical conductivity of glaucophane increases discontinuously with increasing temperature (Fig. 1). At 1.5 GPa, two electrical conductivity discontinuities define three distinct conduction mechanisms in glaucophane corresponding to (i) the solid-state sample prior to dehydration (0.92 ± 0.08 eV), (ii) sample and aqueous fluid following dehydration (0.87 ± 0.08 eV), and (iii) partially molten sample (0.33 ± 0.16 eV) (Fig. 1a). The first discontinuous increase of conductivity to 0.1 S/m occurs at 830 K, which may correspond to glaucophane dehydration and the release of aqueous fluids. The second increase is observed at around 1000 K, possibly due to the melting of glaucophane and the production of a Na-rich melt. The sample exhibits electrical conductivity of 97 S/m at 1258 K. The electrical conductivity of glaucophane at 6 GPa exhibits a sharp increase of conductivity at around 895 K, which corresponds to the breakdown of glaucophane and the release of aqueous fluids in the sample. At 950 K, the sample exhibits electrical conductivity of 1 S/m (Fig. 1b). The impedance spectra of the samples display spectral features that are characteristic of dehydration and melting of a hydrous mineral⁷,²⁶ (Fig. 2).

The observed conductivity variations correspond to textural and chemical changes of the sample, which confirms the dehydration and melting of glaucophane (Fig. 3). The chemical analyses of the sample after the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Impedance spectra of the glaucophane at 1.5 GPa. (a) At 725 K, below dehydration temperature of glaucophane (830 K), the impedance arc resembles grain interior processes, which can be modeled with a resistor- constant phase element (R-CPE) circuit. (b) At 853 K, the impedance spectrum shows the development of an inductive loop in response to the dehydration and the presence of fluid in the sample. The equivalent circuit can be represented as an R-CPE with an inductive component (L). The induction at low frequencies can be explained by the adsorption of ionic species in the electrode surface or erosion of electrodes due to fluid phases. (c) Further increase in temperature to 1033 K results in a sudden decrease in sample resistance. The persisting induction loop in impedance spectra indicates a possible reaction of melt with Ni electrodes.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Back-scattered electrons (BSE) image of the recovered sample after the electrical conductivity measurements. (a) Energy-dispersive X-ray spectroscopy (EDS) compositional map overlapped with a back-scattered electrons (BSE) image of the recovered sample at 1.5 GPa. The phases present in the sample are olivine (green), enstatite (pink), plagioclase (light blue), and melt (indigo). (b) The BSE image showing the glaucophane starting sample recovered after electrical conductivity measurements at 6 GPa.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** A compilation of electrical properties of hydrous minerals, fluid, and melt in the slab. The comparison of the electrical properties of hydrous minerals and their dehydration product (Manthilake et al. 2021a, 2021b<sup>7,8</sup> and references therein) with glaucophane and its dehydration products. The dehydration and melting temperatures of glaucophane (gln) are lower than those observed for other major hydrous minerals in subduction zones.</figcaption></figure>

electrical conductivity measurements are reported in Supplementary Table S1. Previous studies find that the dehydration of glaucophane at high temperature and high pressure follows the reactions: (1) glaucophane = jadeite + enstatite + nyböite + H₂O<sup>31</sup> or (2) glaucophane = jadeite + enstatite + quartz + H₂O<sup>32</sup>. However, the chemical analyses of our sample after the electrical conductivity measurements at 1.5 GPa indicate enstatite, olivine, plagioclase, and melt (Fig. 3a). This disparity can be explained by the dissolution of quartz in the melt and olivine crystallization in consequence via reaction, enstatite = olivine + SiO₂ (in melt)<sup>33</sup>. The plagioclase found in our sample appears to crystallize from the melt during the cooling of the sample. To preserve the mineral assemblages after the dehydration, in our experiment at 6 GPa, we limit the maximum temperature to 950 K. The breakdown of starting glaucophane sample results in jadeite, enstatite, coesite, and secondary amphiboles (cummingtonite and barroisite) (Supplementary Table S2). Omphacite has been observed surrounding jadeite grains (Fig. 3b).

## Discussion

The electrical conductivity of glaucophane before the dehydration is slightly lower than the electrical conductivity of naturally occurring Cl- and/or F-bearing amphiboles, and broadly comparable to the Na-rich plagioclase (albite)<sup>34-36</sup>. This signifies that crystal structure-bound Na<sup>+</sup> may not be as mobile as the defect-bound volatiles or halogens. Na occupies both the alkali (A) and (B) sites in the crystal structure of amphiboles. Previous studies have shown that the preferential incorporation of Na in the large A sites may enhance the electrical conduction in amphiboles<sup>37</sup>. It appears that Na in glaucophane preferentially occupies the B-sites (0.02 mols in the A-site, 1.91 mols in the B-site), thus exerting negligible influence on the conductivity before dehydration.

The fluid released during the dehydration of glaucophane shows comparable electrical conductivity to that of aqueous fluids observed under subduction zone pressure and temperature conditions<sup>26,38-41</sup> (Fig. 4). While Na<sup>+</sup> can be mobile in fluids, the compatible nature of Na to glaucophane may restrict the partition of Na to the fluid phase during the dehydration<sup>42</sup>, thus limiting the role of Na<sup>+</sup> as a charge carrier in the fluid phase. Compared to the electrical conductivity of the fluid at 1.5 GPa, the increase of electrical conductivity associated with dehydration of glaucophane at 6 GPa appears to be larger (Fig. 1). We attribute this observation to a higher fluid fraction in the sample or the higher ion concentrations, such as ionization of H₂O in the supercritical fluid at the elevated pressure<sup>26</sup>.

Na appears to have a profound effect on the electrical conductivity of silicate melts. The melting of glaucophane increases the conductivity by a factor of ~100 compared to the conductivity of the aqueous fluid at 1.5 GPa (Fig. 1a). Our observation suggests that the electrical conduction in silicate melt, produced after dehydration of amphibole, is largely controlled by the mobility of Na<sup>+</sup> ions (Fig. 5). This observation agrees well with previous experimental results, which identified Na<sup>+</sup> as the principal charge carrier in basaltic and albite melts<sup>30</sup>. Compared to Na<sup>+</sup>, the effect of other major cations, such as Ca<sup>2+</sup>, Al<sup>3+</sup>, and Fe<sup>2+/3+</sup>, on the electrical conductivity is found to be negligible<sup>7</sup>. It has also been shown that the presence of H₂O promotes the mobility of Na<sup>+</sup> in silicate melts<sup>30</sup>. The Na-rich hydrous melt produced by the melting of glaucophane exhibits an electrical conductivity of 97 S/m at 1250 K. Extrapolation to high temperatures shows that the Na-rich melts occurring in the subduction zone settings are significantly more conductive than the silicate melts observed in other geological settings<sup>43-45</sup> (Fig. 4). In contrast, the hydrous silicate melts produced by melting of tremolite, actinolite, ferri-kaersutite, and hornblende amphibole compositions show a convergence of electrical conductivities around 1 S/m above 1250 K (Fig. 5). These values closely resemble the electrical conductivity values obtained for hydrous basalts that are commonly found in volcanic-arc settings<sup>43,46</sup> (Fig. 4).
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Electrical conductivity of silicate melt as a function of Na contents. The electrical conductivity of melts produced by dehydration and melting of amphiboles such as, ferri-kaersutite⁷, actinolite⁷, hastingsite⁷, tremolite⁷, and glaucophane (present study) at 1.5 GPa are plotted against the Na₂O contents in the melt. The uncertainties in composition and electrical conductivities are smaller than the symbol size. The highest temperatures of electrical measurements are shown next to individual symbols.</figcaption></figure>

**Geological significance of dehydration melting of glaucophane in the slab.** The phase relations at the solidus of the subducted oceanic crust have been investigated by several studies¹,²,⁴,⁷. Typically, amphiboles in the slab are thought to break down at pressures greater than 2.5 GPa. The dehydration reaction of amphibole starts at around 50 km depth in hot subduction systems and continues up to 75 km depending on the slab surface temperature at 2.5 GPa. The stability of amphibole could be extended to a depth beyond 75 km by substituting OH with F in the amphibole crystal structure²,²⁹. In our 6 GPa experiment, glaucophane goes through a dehydration reaction forming jadeite + enstatite + quartz. We also observe the formation of secondary F-bearing amphibole (Supplementary Table S2) attesting to the stability of the amphibole phase at this pressure condition. The crystallization of F-bearing amphiboles has also been demonstrated by other amphibole dehydration experiments conducted under subduction zone pressure and temperature conditions⁷.

**Electrical conductivity of glaucophane as a geothermometer.** The dehydration temperature we observed for glaucophane using the in situ electrical conductivity measurements appears to occur at lower temperatures compared to the decomposition of OH⁻ and F⁻ bearing amphiboles²⁸,²⁹ (Fig. 6). We attribute this decrease in dehydration and melting temperature to the presence of Na. Albite, another principal host of Na in the slab, transforms to jadeite + quartz through a fluid absent reaction at high pressure and temperature⁴⁸,⁴⁹, thus, may not influence the electrical conductivity measurements. The low dehydration and melting temperature associated with the high electrical conductivity makes glaucophane an ideal indicator-mineral for tracking the slab surface temperature in subduction zones (Fig. 7).

As a case study, we investigated the 2D magnetotelluric (MT) profile at the Sodo ridge, where the old and cold Pacific plate subducts beneath North-East Japan⁵⁰,⁵¹. The MT profile of the Sodo ridge provides a case for the extreme scenario found in subduction systems, where the high electrical anomaly appears to occur at depths greater than 150 km. The conductivity-depth profile at the Sodo ridge displays an increase of electrical conductivity from 0.01 to 0.1 S/m at the slab surface at 150–170 km depths, which is the characteristic electrical conductivity for aqueous fluids⁸,²⁶. This apparent dehydration event is followed by a further increase of conductivity to more than 1 S/m at the slab surface at a depth of 185 km⁵⁰. Given that glaucophane breaks down at relatively low temperatures, and any other types of melting in the slab cannot occur before the melting of glaucophane, we consider that the extremely high conductivity displayed at 185 km depth may be related to the melting of F-bearing glaucophane and the subsequent release of Na-rich melt in the slab.

Assuming that incorporation of F stabilizes amphiboles below depths of 75 km and the high conductivity anomaly is caused by the melting of F-bearing glaucophane; we can estimate the depth where the blueschist-eclogite transition happens in the slab. In the case of Sodo ridge, the blueschist-eclogite transition can be fixed at a depth of 160 km (~5.5 GPa) and a temperature of ~887 K (Fig. 8). By setting the solidus of glaucophane parallel to the dehydration curve, we could estimate the onset of melting at 185 km depth (~6.5 GPa) and 1070 K. Our calculations suggest a thermal gradient of 7.3 K/km along the slab surface for the North-East Japan subduction zone. We conclude that the precise temperature for slab surfaces can be derived from the phase relations of glaucophane. These estimates can benchmark the slab-surface temperatures and provide a temperature-fixed point for thermal models of subduction zones.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** The stability of glaucophane at high-pressure and high-temperature. The solid red line indicates a linear fit through experimental data points obtained at 1.5 and 6 GPa. The red dashed line is drawn parallel to the dehydration curve through the melting temperature obtained at 1.5 GPa. The dark blue, light blue, and green lines are the dehydration and melting temperatures observed for F and Na-free amphibole²⁹, Na-free amphibole with 0.43 mol F replacing (OH⁻)²⁸, and Na-free amphibole with 1 mol F²⁹ at the pressure range from 3.0 to 5.5 GPa, respectively. *am* amphibole; *gln* glaucophane.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** The electrical conductivity based geothermometer. The locus of dehydration of glaucophane in the slab is shown in blue circles and the locus of melting of glaucophane is shown in pink circles. The characteristic electrical conductivity at dehydration of glaucophane is 0.01 S/m, the electrical conductivity at the onset of melting of glaucophane is 0.1 S/m. In comparison, the electrical conductivity of the hydrous slab is about 0.001 S/m. The corresponding slab surface temperatures are shown at different pressures.</figcaption></figure>

## Methods

Finely ground powders of glaucophane (from Sesia-Lanzo Zone, Aosta Valley, Western Alps) were prepared from inclusion-free natural crystals, handpicked under a binocular microscope. The chemical composition of glaucophane was determined before the electrical conductivity measurements by electron probe microanalysis using a Cameca SxFiveTactis electron microprobe operating at an accelerating voltage of 15 kV and a beam current of 20 nA. The composition of the starting glaucophane is reported in Supplementary Table 1.
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Application of the glaucophane geothermometer to N–E Japan subduction zone. Two temperature-fixed points for the N–E Japan subduction slab are determined using the electrical conductivity model discussed in Fig. 6. The electrical conductivity structure of the subduction system is taken from Toh et al. 2006⁵⁰. The green, blue, and red zones were determined based on the electrical conductivity of the wedge mantle adjacent to the slab surface. The areas marked by dashedlines indicate high conductivity regions in the mantle wedge.</figcaption></figure>

High pressure–temperature experiments were performed using a 1500-ton multi-anvil apparatus installed at the Laboratoire Magmas et Volcans. The experiments were performed in an 18/11 assembly, with an 18 mm edge length Cr₂O₃-doped MgO octahedral pressure medium and 11 mm WC anvil truncation edge length (Supplementary Figure S1). Cylindrical sample specimens for electrical conductivity measurements were synthesized by hot pressing of amphibole power in rhenium (Re) capsule at 1.5 GPa and 700 K for 1 h. Pre-sintered cylindrical samples were placed in a polycrystalline MgO capsule, which electrically insulated the sample from the furnace during the measurements. Two Ni disks placed at the top and bottom of the sample served as electrodes for the electrical conductivity measurements. These Ni disks also served as an oxygen buffer, controlling the oxygen fugacity in the sample close to the Ni-NiO buffer. A W₉₅Re₅–W₇₄Re₂₆ thermocouple junction was placed at one side of the sample, which monitored the temperature. One cable formed the thermocouple, and a separate W₉₅Re₅ cable was placed at the opposite side of the sample, connected to the impedance spectroscopy for the electrical conductivity measurements. MgO ceramic sleeves insulated the electrode wires from the furnace. All ceramic assembly parts, including the pressure media, were baked at 1273 K for more than 12 h and stored at 400 K in high-vacuum furnaces (< 100 mTorr) before assembling. This step reduced the exposure of assembly components to atmospheric moisture and other impurities.

The electrical conductivity of glaucophane samples was determined at pressures of 1.5 and 6 GPa, based on the impedance spectroscopy method with a Modulab MTS Impedance gain-phase analyzer in the frequency range of 10⁶–10¹ Hz. Sample resistance was measured in 30–50 K temperature steps. The maximum temperatures attained were 1258 K at 1.5 GPa and 950 K at 6 GPa.

Polycrystalline samples can be characterized by a combination of resistor-capacitor (R–C) or resistor-constant phase element (R–CPE) circuits, and the resistance can be obtained by fitting the impedance spectra to appropriate equivalent circuits. Once the sample resistance has been determined, conductivity can be calculated using the sample diameter and length measured after each experiment, assuming the sample geometry remained unchanged during the experiment. The activation enthalpy (ΔH) of each conduction mechanism can be obtained by fitting the data to the Arrhenius equation, $σ = σ₀exp(-ΔH/RT)$, where σ is the electrical conductivity (S/m), T is the absolute temperature, σ₀ is the pre-exponential factor (S/m), and R is the gas constant (J/K/mol).

After an electrical conductivity experiment, cross-sections of the experimental run products were investigated using an electron probe micro analyzer and energy dispersive X-ray spectroscopy (EDS) chemical mapping with a JEOL JSM-5910LV scanning electron microscope, to identify mineral assemblages and the dehydration reaction products. The compositions of the fluid and melt phases were determined by mass-balance calculations based on the mineral proportions and their chemical compositions of recovered samples after the dehydration/melting events⁷.

## Data availability

All data generated or analyzed during this study are included in this article. The raw electrical conductivity data are available from the corresponding author on reasonable request.

Received: 18 May 2021; Accepted: 17 August 2021

Published online: 09 September 2021
## References

1. Schmidt, M. W. & Poli, S. Experimentally based water budgets for dehydrating slabs and consequences for arc magma generation. *Earth Planet. Sci. Lett.* **163**, 361–379 (1998).

2. Scambelluri, M., Müntener, O., Ottolini, L., Pettke, T. T. & Vannucci, R. The fate of B, Cl and Li in the subducted oceanic mantle and in the antigorite breakdown fluids. *Earth Planet. Sci. Lett.* **222**, 217–234 (2004).

3. Ishimaru, S., Arai, S., Ishida, Y., Shirasaka, M. & Okrugin, V. M. Melting and multi-stage metasomatism in the mantle wedge beneath a frontal arc inferred from highly depleted peridotite xenoliths from the avacha volcano, Southern Kamchatka. *J. Petrol.* **48**, 395–433 (2007).

4. Hirschmann, M. M. Mantle solidus: Experimental constraints and the effects of peridotite composition. *Geochem. Geophys. Geosyst.* **1**, 10 (2000).

5. Symonds, R. B., Reed, M. H. & Rose, W. I. Origin, speciation, and fluxes of trace-element gases at Augustine volcano, Alaska: Insights into magma degassing and fumarolic processes. *Geochem. Cosmochim. Acta* **56**, 633–657 (1992).

6. Liebscher, A. Aqueous fluids at elevated pressure and temperature. *Geofluids* **10**, 3–19 (2010).

7. Manthilake, G., Koga, K. T., Peng, Y. & Mookherjee, M. Halogen bearing amphiboles, aqueous fluids, and melts in subduction zones: Insights on halogen cycle from electrical conductivity. *J. Geophys. Res. Earth* **126**, e2020jb021339 (2021).

8. Manthilake, G., Mookherjee, M. & Miyajima, N. Insights on the deep carbon cycle from the electrical conductivity of carbon-bearing aqueous fluids. *Sci. Rep.* **1**, 1–10 (2021).

9. Spandler, C., Hermann, J., Arculus, R. & Mavrogenes, J. Redistribution of trace elements during prograde metamorphism from lawsonite blueschist to eclogite facies; implications for deep subduction-zone processes. *Contrib. Miner. Petrol.* **146**, 205–222 (2003).

10. Peacock, S. M. Thermal structure and metamorphic evolution of subducting slabs. *Geophys. Monogr. Ser.* **138**, 7–22 (2004).

11. Tatsumi, Y. Formation of the volcanic front in subduction zones. *Geophys. Res. Lett.* **12**, 717–720 (1986).

12. Pawley, A. R. & Holloway, J. R. Water sources for subduction zone volcanism: New experimental constraints. *Science (80-)* **260**, 664–667 (1993).

13. Liu, J., Bohlen, S. R. & Ernst, W. G. Stability of hydrous phases in subducting oceanic crust. *Earth Planet. Sci. Lett.* **143**, 161–171 (1996).

14. Zhao, D. Global tomographic images of mantle plumes and subducting slabs: Insight into deep Earth dynamics. *Phys. Earth Planet. Inter.* **146**, 3–34 (2004).

15. Syracuse, E. M. et al. The global range of subduction zone thermal models. *Phys. Earth Planet. Inter.* **183**, 73–90 (2010).

16. van Keken, P. E. et al. A community benchmark for subduction zone modeling. *Phys. Earth Planet. Inter.* **171**, 187–197 (2008).

17. Cooper, L. B. et al. Global variations in H2O/Ce: 1. Slab surface temperatures beneath volcanic arcs. *Geochem. Geophys. Geosyst.* **13**, 3024 (2012).

18. Plank, T., Cooper, L. B. & Manning, C. E. Emerging geothermometers for estimating slab surface temperatures. *Nat. Geosci.* **2**, 611–615 (2009).

19. Rose-Koga, E. F. et al. Volatile (F and Cl) concentrations in Iwate olivine-hosted melt inclusions indicating low-temperature subduction 5. Volcanology geofluid processes in subduction zones and mantle dynamics. *Earth Planets Sp.* **66**, 1–12 (2014).

20. Hacker, B. R., Abers, G. A. & Peacock, S. M. Subduction factory 1. Theoretical mineralogy, densities, seismic wave speeds, and H2O contents. *J. Geophys. Res.* **108**, 1–26 (2003).

21. Penniston-Dorland, S. C. & Ferry, J. M. Coupled dichotomies of apatite and fluid composition during contact metamorphism of siliceous carbonate rocks. *Am. Miner.* **90**, 1606–1618 (2005).

22. Schmidt, M. W. & Poli, S. The stability of lawsonite and zoisite at high pressures: Experiments in CASH to 92 kbar and implications for the presence of hydrous phases in subducted lithosphere. *Earth Planet. Sci. Lett.* **124**, 105–118 (1994).

23. Pawley, A. R. The pressure and temperature stability limits of lawsonite: Implications for H2O recycling in subduction zones. *Contrib. Miner. Petrol.* **118**, 99–108 (1994).

24. Chantel, J., Mookherjee, M. & Frost, D. J. The elasticity of lawsonite at high pressure and the origin of low velocity layers in subduction zones. *Earth Planet. Sci. Lett.* **349–350**, 116–125 (2012).

25. Mookherjee, M. & Bezacier, L. The low velocity layer in subduction zone: Structure and elasticity of glaucophane at high pressures. *Phys. Earth Planet. Inter.* **208–209**, 50–58 (2012).

26. Manthilake, G., Mookherjee, M., Bolfan-Casanova, N. & Andrault, D. Electrical conductivity of lawsonite and dehydrating fluids at high pressures and temperatures. *Geophys. Res. Lett.* **42**, 7398–7405 (2015).

27. Pommier, A., Williams, Q., Evans, R. L., Pal, I. & Zhang, Z. Electrical investigation of natural lawsonite and application to subduction contexts. *J. Geophys. Res. Solid Earth* **124**, 1430–1442 (2019).

28. Holloway, J. R. & Ford, C. E. Fluid-absent melting of the fluoro-hydroxy amphibole pargasite to 35 kilobars. *Earth Planet. Sci. Lett.* **25**, 44–48 (1975).

29. Foley, S. High-pressure stability of the fluor- and hydroxy-endmembers of pargasite and K-richterite. *Geochem. Cosmochim. Acta* **55**, 2689–2694 (1991).

30. Ni, H., Keppler, H., Manthilake, M. A. G. M. & Katsura, T. Electrical conductivity of dry and hydrous NaAlSi3O8 glasses and liquids at high pressures. *Contrib. Miner. Petrol.* **162**, 501–513 (2011).

31. Pawley, A. R. Experimental study of the compositions and stabilities of synthetic nyböite and nyböite-glaucophane amphiboles. *Eur. J. Miner.* **4**, 171–192 (1992).

32. Cheng, N., Jenkins, D. D. M. D. & Huang, F. Dehydration of glaucophane in the system Na2O-MgO-Al2O3-SiO2-H2O and the effects of NaCl-, CO2- and silicate-bearing aqueous fluids. *J. Petrol.* **60**, 2369–2386 (2019).

33. Carman, J. H. & Gilbert, M. C. Experimental studies on glaucophane stability. *Am. J. Sci.* **283**, 414–437 (1983).

34. Manthilake, G. et al. The electrical conductivity of liebermannite: Implications for water transport into the Earth's lower mantle. *J. Geophys. Res. Solid Earth* **125**, e2020jb020094 (2020).

35. Yang, X., Keppler, H., McCammon, C. & Ni, H. Electrical conductivity of orthopyroxene and plagioclase in the lower crust. *Contrib. Miner. Petrol.* **163**, 33–48 (2012).

36. Hu, H., Li, H., Dai, L., Shan, S. & Zhu, C. Electrical conductivity of alkali feldspar solid solutions at high temperatures and high pressures. *Phys. Chem. Miner.* **40**, 51–62 (2013).

37. Peng, Y., Manthilake, G. & Mookherjee, M. Electrical conductivity of metasomatized lithology in subcontinental lithosphere. *Am. Miner.* https://doi.org/10.2138/am-2021-7942 (2021).

38. Manthilake, G., Bolfan-Casanova, N., Novella, D., Mookherjee, M. & Andrault, D. Dehydration of chlorite explains anomalously high electrical conductivity in the mantle wedges. *Sci. Adv.* **2**, e1501631 (2016).

39. Hu, H., Dai, L., Li, H., Hui, K. & Sun, W. Influence of dehydration on the electrical conductivity of epidote and implications for high-conductivity anomalies in subduction zones. *J. Geophys. Res. Solid Earth* **122**, 2751–2762 (2017).

40. Hu, H., Dai, L., Li, H., Sun, W. & Li, B. Effect of dehydrogenation on the electrical conductivity of Fe-bearing amphibole: Implications for high conductivity anomalies in subduction zones and continental crust. *Earth Planet. Sci. Lett.* **498**, 27–37 (2018).

41. Sun, W. et al. Effect of dehydration on the electrical conductivity of phyllite at high temperatures and pressures. *Mineral. Petrol.* https://doi.org/10.1007/s00710-017-0494-2 (2017).

42. Manning, C. E. Fluid composition at the blueschist-eclogite transition in the model system Na2O-MgO-Al2O3-SiO2-H2O-HCl. *Swiss Bull. Miner. Petrol.* **78**, 225–242 (1998).
43. Ni, H., Keppler, H. & Behrens, H. Electrical conductivity of hydrous basaltic melts: Implications for partial melting in the upper mantle. *Contrib. Miner. Petrol.* **162**, 637–650 (2011).

44. Yoshino, T., McIsaac, E., Laumonier, M. & Katsura, T. Electrical conductivity of partial molten carbonate peridotite. *Phys. Earth Planet. Inter.* **194–195**, 1–9 (2012).

45. Gaillard, F., Malki, M., Iacono-Marziano, G., Pichavant, M. & Scaillet, B. Carbonatite melts and electrical conductivity in the asthenosphere. *Science* **322**, 1363–1365 (2008).

46. Gailler, L. et al. 3D electrical conductivity imaging of Halema'uma'u lava lake (Kilauea volcano). *J. Volcanol. Geotherm. Res.* https://doi.org/10.1016/j.jvolgeores.2019.06.001 (2019).

47. Poli, S. & Schmidt, M. W. Petrology of subducted slabs. *Annu. Rev. Earth Planet. Sci.* **30**, 207–235 (2002).

48. Liu, L. G. High-pressure phase transformations of albite, jadeite and nepheline. *Earth Planet. Sci. Lett.* **37**, 438–444 (1978).

49. Mookherjee, M. et al. Pressure induced elastic softening in framework aluminosilicate-albite (NaAlSi3O8). *Sci. Rep.* **6**, 1–10 (2016).

50. Toh, H. et al. Two-dimensional electrical section beneath the eastern margin of Japan sea. *Geophys. Res. Lett.* **33**, 1–5 (2006).

51. Ichiki, M., Baba, K., Toh, H. & Fuji-ta, K. An overview of electrical conductivity structures of the crust and upper mantle beneath the northwestern Pacific, the Japanese Islands, and continental East Asia. *Gondwana Res.* **16**, 545–562 (2009).

## Acknowledgements

We thank E. Voyer for helping with the SEM analyses, J-L Devidal for the electron microprobe analyses, Antoine Mathieu, Jean-Louis Fruquière and Cyrille Guillot for the technical assistance. Christian Nicollet is thanked for providing the glaucophane sample. The constructive comments from Lidong Dai and an anonymous reviewer help improve the quality of the manuscript. GM acknowledges funding from the INSU-CNRS. M.M. is supported by the US National Science Foundation grants, EAR 1753125, 1763215. This research was financed by the French Government Laboratory of Excellence initiative no. ANR-10-LABX-0006, the Région Auvergne, and the European Regional Development Fund (ClerVolc contribution number 498).

## Author contributions

G.M. and Y.P. performed the experiments. G.M., K.T.K. analyzed the data. G.M., K.T.K., Y.P., and M.M. wrote the paper. All authors discussed the results and commented on the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-97317-0.

**Correspondence** and requests for materials should be addressed to G.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021