skip to main content

Title: High-pressure single-crystal elasticity and thermal equation of state of omphacite and their implications for the seismic properties of eclogite in the Earth's interior
Omphacite is a major mineral phase of eclogite, which provides the main driving force for the slab subduction into the Earth's interior. We have measured the single‐crystal elastic moduli of omphacite at high pressures for the first time up to 18 GPa at ambient temperature using Brillouin spectroscopy. A least‐squares fit of the velocity‐pressure data to the third‐order finite strain equation of state yields KS0'=4.5(3), G0'=1.6(1) with ρ0=3.34(1) g/cm3, KS0=123(3) GPa, and G0=74(2) GPa. In addition, the synchrotron single‐crystal X‐ray diffraction data have been collected up to 18 GPa and 700 K. The fitting to Holland‐Powell thermal‐pressure equation of state yields KT0'=4.6(5) and α0=2.7(8) ×10‐5 K‐1. Based on the obtained thermoelastic parameters of omphacite, the anisotropic seismic velocities of eclogite are modeled and compared with pyrolite between 200 to 500 km. The largest contrast between the eclogite and pyrolite in terms of seismic properties is observed between ~310 to 410 km.
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of Geophysical Research: Solid Earth
Sponsoring Org:
National Science Foundation
More Like this
  1. Clinopyroxene (Cpx) is commonly believed to be the best structural water (hydrogen) carrier among all major upper mantle nominally anhydrous minerals (NAMs). In this study, we have measured the single-crystal elastic properties of a Cpx, a natural omphacite with ~710 ppm water at ambient pressure (P) and temperature (T) conditions. Utilizing the single-crystal X-ray diffraction (XRD) and electron microprobe data, the unit cell parameters and density were determined as a = 9.603(9) Å, b = 8.774(3) Å, c = 5.250(2) Å, β = 106.76(5)o, V = 255.1(4) Å3, and ρ = 3.340(6) g/cm3. We performed Brillouin spectroscopy experiments on fourmore »single crystals along a total of 52 different crystallographic directions. The best-fit single-crystal elastic moduli (Cijs), bulk and shear moduli were determined as: C11 = 245(1) GPa, C22 = 210(2) GPa, C33 = 249.6(9) GPa, C44 = 75.7(9) GPa, C55 = 71.2(5) GPa, C66 = 76(1) GPa, C12 = 85(2) GPa, C13 = 70(1) GPa, C23 = 66(2) GPa, C15 = 8.0(6) GPa, C25 = 6(1) GPa, C35 = 34.7(6) GPa, and C46 = 8.7(7) GPa, KS0 = 125(3) GPa, and G0 = 75(2) GPa, respectively. Compared with the anticipated elastic properties of an anhydrous omphacite with the same chemical composition, our results indicate that the incorporation of ~710 ppm structural water has no resolvable effect on the aggregate elastic properties of omphacite, although small differences (up to ~9 GPa) were observed in C13, C25, C44, and C66.« less
  2. Abstract The 13 single-crystal adiabatic elastic moduli (Cij) of a C2/c jadeite sample close to the ideal composition (NaAlSi2O6) and a natural P2/n diopside-rich omphacite sample have been measured at ambient condition by Brillouin spectroscopy. The obtained Cij values for the jadeite sample are: C11 = 265.4(9) GPa, C22 = 247(1) GPa, C33 = 274(1) GPa, C44 = 85.8(7) GPa, C55 = 69.3(5) GPa, C66 = 93.0(7) GPa, C12 = 84(1) GPa, C13 = 66(1) GPa, C23 = 87(2) GPa, C15 = 5.4(7) GPa, C25 = 17(1) GPa, C35 = 28.7(6) GPa, C46 = 14.6(6) GPa. Voigt-Reuss-Hill averaging of themore »Cij values yields aggregate bulk modulus KS = 138(3) GPa and shear modulus G = 84(2) GPa for jadeite. Systematic analysis combing previous single-crystal elasticity measurements within the diopside-jadeite solid solution indicates that the linear trends are valid for most Cij values. The νp and νs of omphacite decrease with diopside content, though the velocity changes are small as diopside component exceeds 70%. We also found that both the isotropic νp and νs, as well as the seismic anisotropy of eclogite, changed strongly with the bulk-chemical composition. The relationship between the anisotropic velocities of eclogite and the chemical composition can be a useful tool to trace the origin of the eclogitic materials in the Earth's mantle.« less
  3. Abstract As a major nominally anhydrous mineral (NAM) in the Earth’s upper mantle, orthopyroxene could host up to several hundred parts per million H2O in its crystal structure and transport the H2O to the deep Earth. To study the effect of structural H2O on the elasticity of orthopyroxene, we have measured the single-crystal elasticity of Mg1.991Al0.065Si1.951O6 with 842–900 ppm H2O and 1.64 ± 0.20 wt% Al2O3 at ambient conditions using Brillouin spectroscopy. The best-fit single-crystal elastic moduli (Cijs), bulk (KS0), and shear (G0) modulus of the hydrous Al-bearing orthopyroxene were determined as: C11 = 235(2) GPa, C22 = 173(2) GPa,more »C33 = 222(2) GPa, C44 = 86(1) GPa, C55 = 82(1) GPa, C66 = 82(1) GPa, C12 = 75(3) GPa, C13 = 67(2) GPa, and C23 = 49(2) GPa, KS0 = 111(2) GPa, and G0 = 78(1) GPa. Systematic analysis based on the results presented in this and previous studies suggests that the incorporation of 842–900 ppm H2O would increase C13 by 12.0(7)% and decrease C23 by 8.6(8)%. The effects on C11, C22, C33, C44, C66, KS0, and VP are subtle if not negligible when considering the uncertainties. The C55, C12, G0, and VS are not affected by the presence of structural H2O. Although laboratory experiments show that Fe,Al-bearing orthopyroxenes can host up to 0.8 wt% H2O in its structure, future high-pressure-temperature elasticity measurements on orthopyroxene with higher H2O content are needed to help better quantify this effect.« less
  4. Abstract Calcium carbonate (CaCO3) is one of the most abundant carbonates on Earth's surface and transports carbon to Earth's interior via subduction. Although some petrological observations support the preservation of CaCO3 in cold slabs to lower mantle depths, the geophysical properties and stability of CaCO3 at these depths are not known, due in part to complicated polymorphic phase transitions and lack of constraints on thermodynamic properties. Here we measured thermal equation of state of CaCO3-Pmmn, the stable polymorph of CaCO3 through much of the lower mantle, using synchrotron X-ray diffraction in a laser-heated diamond-anvil cell up to 75 GPa andmore »2200 K. The room-temperature compression data for CaCO3-Pmmn are fit with third-order Birch-Murnaghan equation of state, yielding KT0 = 146.7 (±1.9) GPa and K′0 = 3.4(±0.1) with V0 fixed to the value determined by ab initio calculation, 97.76 Å3. High-temperature compression data are consistent with zero-pressure thermal expansion αT = a0 + a1T with a0 = 4.3(±0.3)×10-5 K-1, a1 = 0.8(±0.2)×10-8 K-2, temperature derivative of the bulk modulus (∂KT/∂T)P = –0.021(±0.001) GPa/K; the Grüneisen parameter γ0 = 1.94(±0.02), and the volume independent constant q = 1.9(±0.3) at a fixed Debye temperature θ0 = 631 K predicted via ab initio calculation. Using these newly determined thermodynamic parameters, the density and bulk sound velocity of CaCO3-Pmmn and (Ca,Mg)-carbonate-bearing eclogite are quantitatively modeled from 30 to 80 GPa along a cold slab geotherm. With the assumption that carbonates are homogeneously mixed into the slab, the results indicate the presence of carbonates in the subducted slab is unlikely to be detected by seismic observations, and the buoyancy provided by carbonates has a negligible effect on slab dynamics.« less
  5. The high-pressure structure and stability of the calcic amphibole tremolite (Ca2Mg5Si8O22(OH)2) was investigated to ~40 GPa at 300 K by single-crystal X-ray diffraction using synchrotron radiation. C2/m symmetry tremolite displays a broader metastability range than previously studied clinoamphiboles, exhibiting no first-order phase transition up to 40 GPa. Axial parameter ratios a/b and a/c, in conjunction with finite strain versus normalized pressure trends, indicate that changes in compressional behavior occur at pressures of ~5 and ~20 GPa. An analysis of the finite strain trends, using third-order Birch-Murnaghan equations of state, resulted in bulk moduli (𝐾) of 72(7), 77(2), and 61(1) GPamore »for the compressional regimes from 0-5 GPa (regime I), 5-20 GPa (II), and above 20 GPa (III), respectively, and accompanying pressure-derivatives of the bulk moduli (𝐾′) of 8.6(42), 6.0(3), and 10.0(2). The results are consistent with first-principle theoretical calculations of tremolite elasticity. The axial compressibility ratios of tremolite, determined as 𝛽a : 𝛽b : 𝛽c = 2.22:1.0:0.78 (regime I), 2.12:1.0:0.96 (II), and 1.03:1.0:0.75 (III), demonstrate a substantial reduction of the compressional anisotropy of tremolite at high pressures, which is a notable contrast with the increasingly anisotropic compressibility observed in the high-pressure polymorphs of the clinoamphibole grunerite. The shift in compression-regime at 5 GPa (I-II) transition is ascribed to stiffening along the crystallographic a-axis corresponding to closure of the vacant A-site in the structure, and a shift in the topology of the a-oriented surfaces of the structural I-beam from concave to convex. The II-III regime shift at 20 GPa corresponds to an increasing rate of compaction of the Ca-polyhedra and increased distortion of the Mg-octahedral sites, processes which dictate compaction in both high-pressure compression-regimes. Bond-valence analyses of the tremolite structure under pressure show dramatic overbonding of the Ca-cations (75% at 30 GPa), with significant Mg-cation overbonding as well (40%). These imply that tremolite’s notable metastability range hinges on the calcium cation’s bonding environment. The 8-fold coordinated Ca-polyhedron accommodates significant compaction under pressure, while the geometry of the Ca-O polyhedron becomes increasingly regular and inhibits the reorientation of the tetrahedral chains that generate phase transitions observed in other clinoamphiboles. Peak/background ratio of diffraction data collected above 40 GPa and our equation of state determination of bulk moduli and compressibilities of tremolite in regime III, in concert with the results of our previous Raman study, suggest that C2/m tremolite may be approaching the limit of its metastability above 40 GPa. Our results have relevance for both the metastable compaction of tremolite during impact events, and for possible metastable persistence of tremolite within cold subduction zones within the Earth.« less