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: Electrical conductivity of metasomatized lithology in subcontinental lithosphere
Abstract A plausible origin of the seismically observed mid-lithospheric discontinuity (MLD) in the subcontinental lithosphere is mantle metasomatism. The metasomatized mantle is likely to stabilize hydrous phases such as amphiboles. The existing electrical conductivity data on amphiboles vary significantly. The electrical conductivity of hornblendite is much higher than that of tremolite. Thus, if hornblendite truly represents the amphibole varieties in MLD regions, then it is likely that amphibole will cause high electrical conductivity anomalies at MLD depths. However, this is inconsistent with the magnetotelluric observations across MLD depths. Hence, to better understand this discrepancy in electrical conductivity data of amphiboles and to evaluate whether MLD could be caused by metasomatism, we determined the electrical conductivity of a natural metasomatized rock sample. The metasomatized rock sample consists of ~87% diopside pyroxene, ~9% sodium-bearing tremolite amphibole, and ~3% albite feldspar. We collected the electrical conductivity data at ~3.0 GPa, i.e., the depth relevant to MLD. We also spanned a temperature range between 400 to 1000 K. We found that the electrical conductivity of this metasomatized rock sample increases with temperature. The temperature dependence of the electrical conductivity exhibits two distinct regimes. At low temperatures <700 K, the electrical conductivity is dominated by the conduction in the solid state. At temperatures >775 K, the conductivity increases, and it is likely to be dominated by the conduction of aqueous fluids due to partial dehydration. The main distinction between the current study and the prior studies on the electrical conductivity of amphiboles or amphibole-bearing rocks is the sodium (Na) content in amphiboles of the assemblage. Moreover, it is likely that the higher Na content in amphiboles leads to higher electrical conductivity. Pargasite and edenite amphiboles are the most common amphibole varieties in the metasomatized mantle, and our study on Na-bearing tremolite is the closest analog of these amphiboles. Comparison of the electrical conductivity results with the magnetotelluric observations constrains the amphibole abundance at MLD depths to <1.5%. Such a low-modal proportion of amphiboles could only reduce the seismic shear wave velocity by 0.4–0.5%, which is significantly lower than the observed velocity reduction of 2–6%. Thus, it might be challenging to explain both seismic and magnetotelluric observations at MLD simultaneously.  more » « less
Award ID(s):
1763215 1753125 1638752
PAR ID:
10318509
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
American Mineralogist
Volume:
107
Issue:
3
ISSN:
0003-004X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Metasomatized mantle xenoliths containing hydrous minerals, such as amphiboles, serpentine, and phlogopite, likely represent the potential mineralogical compositions of the metasomatized upper mantle, where low seismic velocities are commonly observed. This study presents the first experimentally determined single‐crystal elasticity model of an Fe‐free near Ca, Mg‐endmember amphibole tremolite at high pressure and/or temperature conditions (maximum pressure 7.3(1) GPa, maximum temperature 700 K) using Brillouin spectroscopy. We found that sound velocities of amphiboles strongly depend on the Fe content. We then calculated the sound velocities of 441 hydrous‐mineral‐bearing mantle xenoliths collected around the globe, and quantitatively evaluated the roles that amphiboles, phlogopite and serpentine played in producing the low velocity anomalies in the metasomatized upper mantle. 
    more » « less
  2. Abstract Amphiboles are hydrous minerals that are formed in the oceanic crust via hydrothermal alteration. The partial substitution of halogens for OHmakes amphibole one of the principal hosts of Cl and F in the subducting slab. In this study, we investigated the electrical conductivity of a suite of halogen bearing amphibole minerals at 1.5 GPa up to 1,400 K. The discontinuous electrical behavior indicates dehydration of amphibole at ∼915 K. This is followed by dehydration induced hydrous melting at temperatures above 1,070 K. We find that the released aqueous fluids have an electrical conductivity of ∼0.1 S/m. This high electrical conductivity is likely to explain anomalously high electrical conductivity observed in certain subduction zone settings. This high electrical conductivity of an order of magnitude greater than the electrical conductivity of pure aqueous fluids at similar conditions is likely due to the partitioning of the F and Cl into the aqueous fluids. We also noted that subsequent to the dehydration, secondary phases form due to the breakdown of the primary halogen bearing amphibole. Chemical analyses of these secondary phases indicate that they are repositories of F and Cl. Hence, we infer that upon dehydration of the primary halogen bearing amphibole, first the F and Cl are partitioned into the aqueous fluids and then the halogens are partitioned back to the secondary mineral phases. These secondary minerals are likely to transport the halogen to the deep Earth and may in part explain the halogen concentration observed in ocean island basalt. 
    more » « less
  3. null (Ed.)
    Abstract We investigated the structure, equation of state, thermodynamics, and elastic properties of tremolite amphibole [Ca2Mg5Si8O22(OH)2] up to 10 GPa and 2000 K, using first principles simulations based on density functional perturbation theory. We found that at 300 K, the pressure-volume results can be adequately described by a third-order Birch-Murnaghan equation of state with bulk moduli K0 of 78.5 and 66.3 GPa based on local density approximation (LDA) and generalized gradient approximation (GGA), respectively. We also derived its coefficients of the elastic tensor based on LDA and GGA and found that the LDA result is in good agreement with the experimental results. At 300 K, the shear modulus G0 is 58.0 GPa based on LDA. The pressure derivative of the bulk modulus K′ is 5.9, while that of the shear modulus G′ is 1.3. The second Grüneisen parameter, or δT = [–1/(αKT)](∂KT/∂T)P, is 3.3 based on LDA. We found that at ambient conditions, tremolite is elastically anisotropic with the compressional wave velocity anisotropy AVP being 34.6% and the shear wave velocity anisotropy AVS being 27.5%. At higher pressure corresponding to the thermodynamic stability of tremolite, i.e., ~3 GPa, the AVP reduces to 29.5%, whereas AVS increases to 30.8%. To evaluate whether the presence of hydrous phases such as amphibole and phlogopite could account for the observed shear wave velocity (VS) anomaly at the mid-lithospheric discontinuity (MLD), we used the thermoelasticities of tremolite (as a proxy for other amphiboles), phlogopite, and major mantle minerals to construct synthetic velocity profiles. We noted that at depths corresponding to the mid-lithosphere, the presence of 25 vol% amphibole and 1 vol% phlogopite could account for a VS reduction of 2.3%. Thus based on our thermoelasticity results on tremolite amphibole, it seems that mantle metasomatism could partly explain the MLD. 
    more » « less
  4. Abstract Mid‐lithosphere discontinuities are seismic interfaces likely located within the lithospheric mantle of stable cratons, which typically represent velocities decreasing with depth. The origins of these interfaces are poorly understood due to the difficulties in both characterizing them seismically and reconciling the observations with thermal‐chemical models of cratons. Metasomatism of the cratonic lithosphere has been reported by numerous geochemical and petrological studies worldwide, yet its seismic signature remains elusive. Here, we identify two distinct mid‐lithosphere discontinuities at ∼87 and ∼117 km depth beneath the eastern Wyoming craton and the southwestern Superior craton by analyzing seismic data recorded by two longstanding stations. Our waveform modeling shows that the shallow and deep interfaces represent isotropic velocity drops of 2%–8% and 4%–9%, respectively, depending on the contributions from changes in radial anisotropy and density. By building a thermal‐chemical model including the regional xenolith thermobarometry constraints and the experimental phase‐equilibrium data of mantle metasomatism, we show that the shallow interface probably represents the metasomatic front, below which hydrous minerals such as amphibole and phlogopite are present, whereas the deep interface may be caused by the onset of carbonated partial melting. The hydrous minerals and melts are products of mantle metasomatism, with CO2‐H2O‐rich siliceous melt as a probable metasomatic reagent. Our results suggest that mantle metasomatism is probably an important cause of mid‐lithosphere discontinuities worldwide, especially near craton boundaries, where the mantle lithosphere may be intensely metasomatized by fluids and melts released by subducting slabs. 
    more » « less
  5. Abstract The goal of this study is to constrain the origins of layering in the seismic velocity structure within the cratonic mantle lithosphere (i.e. mid‐lithospheric discontinuities [MLDs]). For long‐lived stations in cratons worldwide, we calculated S‐to‐P converted phase receiver function stacks using time domain deconvolution and a k‐means algorithm to select robust, consistent receiver functions. Negative MLDs appear in only 50% of the receiver function stacks, indicating that negative MLDs are common but intermittent. The negative MLDs correspond to shear velocity drops of 1%–4%, which could be caused by layers of minerals created by metasomatism, although vertical layering in seismic anisotropy cannot be ruled out. In craton interiors, negative MLDs have a lower amplitude (<3% velocity drops) and can be explained by metasomatism of the original Archean mantle. Negative MLD amplitudes increase with decreasing upper mantle shear velocity (toward the outer margins of the cratons), but do not depend on the age of the craton. Thus, negative MLD amplitudes are not dominated by age‐related variations in the cratonic mantle composition, and, instead, are more strongly correlated with proximity to tectonic and metasomatic activity that occurred long after craton formation. Negative MLDs are less numerous among stations that have Paleoproterozoic and Archean thermotectonic ages, consistent with the view that shallow release of slab‐derived fluids during early “warm” subduction was less favorable for negative MLD formation. We also observe velocity gradients below 150 km at stations in craton boundaries and interiors, indicating the presence of seismic velocity changes at the cratonic lithosphere‐asthenosphere boundary and/or Lehmann discontinuity. 
    more » « less