skip to main content


Title: Deformation and Strength of Mantle Relevant Garnets: Implications for the Subduction of Basaltic-rich Crust
Garnet is an important mineral phase in the upper mantle as it is both a key component in bulk mantle rocks, and a primary phase at high-pressure within subducted basalt. Here, we focus on the strength of garnet and the texture that develops within garnet during accommodation of differential deformational strain. We use X-ray diffraction in a radial geometry to analyze texture development in situ in three garnet compositions under pressure at 300 K: a natural garnet (Prp60Alm37) to 30 GPa, and two synthetic majorite-bearing compositions (Prp59Maj41 and Prp42Maj58) to 44 GPa. All three garnets develop a modest (100) texture at elevated pressure under axial compression. Elasto-viscoplastic self-consistent (EVPSC) modeling suggests that two slip systems are active in the three garnet compositions at all pressures studied: {110}<1-21 11> and {001}<110>. We determine a flow strength of ~5 GPa at pressures between 10 to 15 GPa for all three garnets; these values are higher than previously measured yield strengths measured on natural and majoritic garnets. Strengths calculated using the experimental lattice strain differ from the strength generated from those calculated using EVPSC. Prp67Alm33, Prp59Maj41 and Prp42Maj58 are of comparable strength to each other at room temperature, which indicates that majorite substitution does not greatly affect the strength of garnets. Additionally, all three garnets are of similar strength as lower mantle phases such as bridgmanite and ferropericlase, suggesting that garnet may not be notably stronger than the surrounding lower mantle/deep upper mantle phases at the base of the upper mantle.  more » « less
Award ID(s):
2017294 1620423 1855336 1654687
NSF-PAR ID:
10232777
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
American Mineralogist
Volume:
106
ISSN:
0003-004X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:log⁡Dm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlog⁡Dm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions. 
    more » « less
  2. SUMMARY We derive exact expressions for the thermal expansivity, heat capacity and bulk modulus for assemblages with arbitrarily large numbers of components and phases, including the influence of phase transformations and chemical exchange. We illustrate results in simple two-component, two-phase systems, including Mg–Fe olivine-wadsleyite and Ca–Mg clinopyroxene-orthopyroxene and for a multicompontent model of mantle composition in the form of pyrolite. For the latter we show results for the thermal expansivity and heat capacity over the entire mantle pressure–temperature regime to 40 GPa, or a depth of 1000 km. From the thermal expansivity, we derive a new expression for the phase buoyancy parameter that is valid for arbitrarily large numbers of phases and components and which is defined at every point in pressure–temperature space. Results reveal regions of the mantle where the magnitude of the phase buoyancy parameter is larger in magnitude than for those phase transitions that are most commonly included in mantle convection simulations. These regions include the wadsleyite to garnet and ferropericlase transition, which is encountered along hot isentropes (e.g. 2000 K potential temperature) in the transition zone, and the ferropericlase and stishovite to bridgmanite transition, which is encountered along cold isentropes (e.g. 1000 K potential temperature) in the shallow lower mantle. We also show the bulk modulus along a typical mantle isentrope and relate it to the Bullen inhomogeneity parameter. All results are computed with our code HeFESTo, updates and improvements to which we discuss, including the implementation of the exact expressions for the thermal expansivity, heat capacity and bulk modulus, generalization to allow for pressure dependence of non-ideal solution parameters and an improved numerical scheme for minimizing the Gibbs free energy. Finally, we present the results of a new global inversion of parameters updated to incorporate more recent results from experiment and first principles theory, as well as a new phase (nal phase), and new species: Na-majorite and the NaAlO2 end-member of ferropericlase. 
    more » « less
  3. Abstract

    Quartz‐in‐garnet elastic geobarometry (QuiG) pressures in rocks from two Barrovian metamorphic terranes in the western US Cordilleran hinterland exceed pressures determined using chemical thermodynamics by 3–4 kbar. For this study, 135 quartz inclusions from the Funeral Mountains, California, were analysed using QuiG in five garnets from three locations representing metamorphic grades of upper greenschist, lower amphibolite, and middle amphibolite facies. From a second Barrovian terrane, the Wood Hills in northeastern Nevada, 125 quartz inclusions were analysed using QuiG in 14 garnets from a single rock sample metamorphosed to middle amphibolite facies. Pressures determined for rocks in the Funeral Mountains using QuiG and methods rooted in equilibrium thermodynamics yielded consistent pressure differences between locations, but QuiG pressures are higher. Similarly, QuiG pressures determined for rocks in the Wood Hills are higher than pressures determined by equilibrium thermodynamic approaches. Possible explanations for the pressure differences include garnet compositions not reflecting equilibrium, sources of error in thermodynamic calculations such as thermodynamic data or a‐X models, or an unknown source of systematic error that causes QuiG to overestimate pressures of entrapment. To test Raman spectroscopy's ability to reproduce inclusion pressures, pressures were calculated using Raman spectroscopy and synchrotron X‐ray diffraction, which yielded consistent pressures and support the use of the single mode‐shift of the 464 cm−1band of quartz for geobarometry, which simplifies the method by assuming hydrostatic compression of quartz. These results are compared with pressures obtained using Grüneisen tensors and show consistency between these different approaches.

     
    more » « less
  4. null (Ed.)
    The degree to which the Earth’s mantle stores and cycles water in excess of the storage capacity of nominally anhydrous minerals is dependent upon the stability of hydrous phases under mantle-relevant pressures, temperatures, and compositions. Two hydrous phases, phase D and phase H, are stable to the pressures and temperatures of the Earth’s lower mantle, suggesting that the Earth’s lower mantle may participate in the cycling of water. We build on our prior work of density functional theory calculations on phase H with the stability, structure, and bonding of hydrous phases D, and we predict the aluminum partitioning with H in the Al 2 O 3 -SiO 2 -MgO-H 2 O system. We address the solid solutions through a statistical sampling of site occupancy and calculation of the partition function from the grand canonical ensemble. We show that each phase has a wide solid solution series between MgSi 2 O 6 H 2 -Al 2 SiO 6 H 2 and MgSiO 4 H 2 -2 δ AlOOH + SiO 2 , in which phase H is more aluminum rich than phase D at a given bulk composition. We predict that the addition of Al to both phases D and H stabilizes each phase to higher temperatures through additional configurational entropy. While we have shown that phase H does not exhibit symmetric hydrogen bonding at high pressure, we report here that phase D undergoes a gradual increase in the number of symmetric H-bonds beginning at ∼30 GPa, and it is only ∼50% complete at 60 GPa. 
    more » « less
  5. The first in situ quantitative synchrotron X-ray diffraction (XRD) study of plastic strain-induced phase transformation (PT) has been performed on $\alpha-\omega$ PT in ultra-pure, strongly plastically predeformed Zr as an example, under different compression-shear pathways in rotational diamond anvil cell (RDAC). Radial distributions of pressure in each phase and in the mixture, and concentration of $\omega$-Zr, all averaged over the sample thickness, as well as thickness profile were measured. The minimum pressure for the strain-induced $\alpha-\omega$ PT, $p^d_{\varepsilon}$=1.2 GPa, is smaller than under hydrostatic loading by a factor of 4.5 and smaller than the phase equilibrium pressure by a factor of 3; it is independent of the compression-shear straining path. The theoretically predicted plastic strain-controlled kinetic equation was verified and quantified; it is independent of the pressure-plastic strain loading path and plastic deformation at pressures below $p^d_{\varepsilon}$. Thus, strain-induced PTs under compression in DAC and torsion in RDAC do not fundamentally differ. The yield strength of both phases is estimated using hardness and x-ray peak broadening; the yield strength in shear is not reached by the contact friction stress and cannot be evaluated using the pressure gradient. Obtained results open a new opportunity for quantitative study of strain-induced PTs and reactions with applications to material synthesis and processing, mechanochemistry, and geophysics. 
    more » « less