Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp71.0Alm12.6Sps0.7Grs3.5Uvr12.2) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus (
Synchrotron‐based high‐pressure/high‐temperature single‐crystal X‐ray diffraction experiments to ~24 GPa and 700 K were conducted on eclogitic garnets (low‐Fe: Prp28Alm38Grs33Sps1and high‐Fe: Prp14Alm62Grs19Adr3Sps2) and omphacites (low‐Fe: Quad57Jd42Ae1and high‐Fe: Quad53Jd27Ae20), using an externally heated diamond anvil cell. Fitting the pressure‐volume‐temperature data to a third‐order Birch‐Murnaghan equation of state yields the thermoelastic parameters including bulk modulus (
- Award ID(s):
- 1722969
- NSF-PAR ID:
- 10460791
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 46
- Issue:
- 1
- ISSN:
- 0094-8276
- Page Range / eLocation ID:
- p. 179-188
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract K S 0), shear modulus (G 0), their pressure ((∂K S /∂P )T and (∂G /∂P )T ) and temperature ((∂K S /∂T )Pand (∂G/∂T )P) derivatives,K S 0 = 167.7(8) GPa,G 0 = 91.5(5) GPa, (∂K S /∂P )T = 4.3(1), (∂G /∂P )T = 1.4(1), (∂K S /∂T )P =− 0.0175(1) GPa/K and (∂G/∂T )P =− 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity,ρ /V P andρ /V S are much more sensitive to the compositional change of the SCLM. -
Abstract Eclogite thermobarometry is crucial for constraining the depths and temperatures to which oceanic and continental crust subduct. However, obtaining the pressure and temperature (
P–T ) conditions of eclogites is complex as they commonly display high‐variance mineral assemblages, and the mineral compositions only vary slightly withP–T . In this contribution, we present a comparison between two independent and commonly used thermobarometric approaches for eclogites: conventional thermobarometry and forward phase‐equilibrium modelling. We assess how consistent the thermobarometric calculations are using the garnet–clinopyroxene–phengite barometer and garnet–clinopyroxene thermometer with predictions from forward modelling (i.e. comparing the relative differences between approaches). Our results show that the overall mismatch in methods is typically ±0.2–0.3 GPa and ±29–42°C although differences as large as 80°C and 0.7 GPa are possible for a few narrow ranges ofP–T conditions in the forward models. Such mismatch is interpreted as the relative differences among methods, and not as absolute uncertainties or accuracies for either method. For most of the investigatedP–T conditions, the relatively minor differences between methods means that the choice in thermobarometric method itself is less important for geological interpretation than careful sample characterization and petrographic interpretation for derivingP–T from eclogites. Although thermobarometry is known to be sensitive to the assumedX Fe3+of a rock (or mineral), therelative differences between methods are not particularly sensitive to the choice of bulk‐rockX Fe3+, except at high temperatures (>650°C, amphibole absent) and for very large differences in assumedX Fe3+(0–0.5). We find that the most important difference between approaches is the activity–composition (a–x ) relations, as opposed to the end‐member thermodynamic data or other aspects of experimental calibration. When equivalenta–x relations are used in the conventional barometer,P calculations are nearly identical to phase‐equilibrium models (ΔP < 0.1). To further assess the implications of these results for real rocks, we also evaluate common mathematical optimizations of reaction constants used for obtaining the maximumP–T with conventional thermobarometric approaches (e.g. using the highesta Grs2 ×a Prp in garnet and Si content in phengite, and the lowesta Di in clinopyroxene). These approaches should be used with caution, because they may not represent the compositions of equilibrium mineral assemblages at eclogite facies conditions and therefore systematically biasP–T calculations. Assuming method accuracy, geological meaningfulP maxat a typical eclogite facies temperature of ~660°C will be obtained by using the greatesta Di,a Cel, anda Prp and lowesta Grs anda Ms; garnet and clinopyroxene with the lowest Fe2+/Mg ratios may yield geological meaningfulT maxat a typical eclogite facies pressure of 2.5 GPa. -
Abstract Identifying and locating the geochemical and geophysical heterogeneities in the Earth’s interior is one of the most important and challenging tasks for the deep Earth scientists. Subducted oceanic crust metamorphizes into the dense eclogite in the upper mantle and is considered as a major cause of geochemical and geophysical heterogeneities in the deep Earth. In order to detect eclogitic materials inside the Earth, precise measurements of the high pressure‐temperature single‐crystal elasticity of major minerals in eclogite are thus exceedingly important. Omphacite, a Na,Al‐bearing clinopyroxene, constitutes up to 75 vol% of eclogite. In the present study, we performed the first high pressure‐temperature single‐crystal elasticity measurements of omphacite using Brillouin spectroscopy. Utilizing the finite‐strain approach, we obtained the following thermoelastic parameters for omphacite:
K S0’ = 4.5(1),G 0’ = 1.53(5), ∂KS0/∂T = −0.029(5) GPa/K, ∂G0/∂T = −0.013(5) GPa/K, withK S0 = 123(3) GPa,G 0 = 74(2) GPa, andρ 0 = 3.34(1) g/cm3. We found that the seismic velocities of undeformed eclogite are similar to pyrolite at the depths of 200–300 and 410–500 km, thus eclogite is seismically invisible at these depths. Combined with the lattice‐preferred orientations of the omphacite in naturally deformed eclogites, we also modeled seismic anisotropy of eclogite at various pressure‐temperature conditions. A 10 km thick subducted eclogitic crust can result in ∼0.2 s shear wave splitting in the Earth’s upper mantle. -
Abstract Fe‐Al‐bearing bridgmanite may be the dominant host for ferric iron in Earth's lower mantle. Here we report the synthesis of (Mg0.5Fe3+0.5)(Al0.5Si0.5)O3bridgmanite (FA50) with the highest Fe3+‐Al3+coupled substitution known to date. X‐ray diffraction measurements showed that at ambient conditions, the FA50 adopted the LiNbO3structure. Upon compression at room temperature to 18 GPa, it transformed back into the bridgmanite structure, which remained stable up to 102 GPa and 2,600 K. Fitting Birch‐Murnaghan equation of state of FA50 bridgmanite yields
V 0 = 172.1(4) Å3,K 0 = 229(4) GPa withK 0′ = 4(fixed). The calculated bulk sound velocity of the FA50 bridgmanite is ~7.7% lower than MgSiO3bridgmanite, mainly because the presence of ferric iron increases the unit‐cell mass by 15.5%. This difference likely represents the upper limit of sound velocity anomaly introduced by Fe3+‐Al3+substitution. X‐ray emission and synchrotron Mössbauer spectroscopy measurements showed that after laser annealing, ~6% of Fe3+cations exchanged with Al3+and underwent the high‐ to low‐spin transition at 59 GPa. The low‐spin proportion of Fe3+increased gradually with pressure and reached 17–31% at 80 GPa. Since the cation exchange and spin transition in this Fe3+‐Al3+‐enriched bridgmanite do not cause resolvable unit‐cell volume reduction, and the increase of low‐spin Fe3+fraction with pressure occurs gradually, the spin transition would not produce a distinct seismic signature in the lower mantle. However, it may influence iron partitioning and isotopic fractionation, thus introducing chemical heterogeneity in the lower mantle. -
Abstract We present ab initio (LDA + U
s c ) studies of high‐temperature and high‐pressure elastic properties of pure as well as iron‐bearing (ferrous, Fe2+, and ferric, Fe3+) and aluminum‐bearing MgSiO3postperovskite, the likely dominant phase in the deep lower mantle of the Earth. Thermal effects are addressed within the quasiharmonic approximation by combining vibrational density of states and static elastic coefficients. Aggregate elastic moduli and sound velocities for the Mg end members are successfully compared with scarce experimental data available. Effects of iron (Fe) and aluminum (Al) substitutions on elastic properties and their pressure and temperature dependence have been thoroughly investigated. At the observed perovskite to postperovskite transition (P = 125 GPa andT = 2,500 K), compressional and shear velocities increase by 0–1% and 1.5–3.75%, respectively. This observation is consistent with some seismic studies of the D′ ′ discontinuity beneath the Caribbean, which suggests that our robust estimates of elastic properties of the postperovskite phase will be very helpful to understand lateral velocity variations in the deep lower mantle region and to constrain its composition and thermal structure.