Elsevier

Lithos

Volumes 372–373, 1 November 2020, 105650
Lithos

Research Article
Revaluation of “equilibrium” P-T paths from zoned garnet in light of quartz inclusion in garnet (QuiG) barometry

https://doi.org/10.1016/j.lithos.2020.105650Get rights and content

Highlights

  • QuiG barometry potentially provides an estimate of garnet growth conditions

  • QuiG barometry in some localities is consistent with peak metamorphic estimates and in some localities not

  • QuiG barometry suggests that some garnets grow nearly isothermally and isobarically

  • QuiG barometry appears to preserve entrapment conditions through multiple metamorphic events

Abstract

Metamorphic P–T paths provide valuable constraints on the tectonics of collisional orogens and many such paths have been deduced from analysis of chemical zoning in garnet. Inclusion barometry (e.g. quartz-in-garnet or QuiG) provides complementary and sometimes contradictory information regarding the interpretation of garnet zoning and the P–T paths derived therefrom. For example, QuiG barometry on several generations of garnet from well-characterized samples from Fall Mountain, New Hampshire, is consistent with the P–T path inferred from other methods. However, QuiG barometry is inconsistent with P–T paths inferred from garnet zoning from the Orfordville Belt (Vermont), Townshend Dam (Vermont), the Connecticut Valley Trough (Vermont), and Sifnos (Greece). New data from the Perry Mountain Formation, southeastern New Hampshire, suggests that QuiG may preserve a record of polymetamorphic events. Specifically, QuiG barometry in coticule samples is interpreted to record a previously unrecognized early medium pressure (ca. 0.9 GPa) metamorphism and a later low pressure (ca. 0.3 GPa) metamorphism.

Introduction

Garnet has served as a central focus of metamorphic studies for the discernment of P–T paths because it has the capability of storing a record of its history through the crust through inclusion suites, chemical zoning, and age zoning (e.g. Ague and Carlson, 2013; Baxter et al., 2013; Baxter and Scherer, 2013; Caddick and Kohn, 2013). Although it is generally acknowledged that metamorphic rocks cannot react if in equilibrium and a finite amount of overstepping is required to drive metamorphic processes, it is generally assumed that the degree to which metamorphic rocks are out of equilibrium is sufficiently small that utilization of equilibrium constraints does not introduce a substantial error in the calculation of metamorphic P–T paths (e.g. Caddick and Kohn, 2013).

However, a number of studies have presented observations that challenge this assumption of a close approach to equilibrium throughout the duration of garnet growth. As early as 1969, Hollister (1969) demonstrated that overstepping of equilibrium reaction boundaries was required to explain textural relations of aluminosilicate minerals (staurolite, kyanite, andalusite and sillimanite) from the Kwoiek contact aureole in British Columbia. More recently, a pioneering study by Waters and Lovegrove (2002) provided unambiguous textural data documenting the sequence of porphyroblast growth in the contact aureole of the Bushveld complex and demonstrated that it was not the sequence that was predicted from equilibrium calculations, and drew the conclusion that considerable overstepping of the equilibrium phase boundary was required for porphyroblast nucleation. Pattison and Tinkham (2009) reach similar a similar conclusion for porphyroblast growth in the Nelson contact aureole by comparing the spacing of isograds mapped in the field with the spacing of isograds that would be predicted from equilibrium calculations and thermal modeling. Pattison and Spear (2018) presented compelling evidence that not only might garnet nucleation be overstepped during regional metamorphism, but also staurolite and aluminosilicates.

The introduction of a new tool into the petrologist's arsenal — inclusion barometry — has provided a technique for evaluating the extent of overstepping required for porphyroblast nucleation. In particular, quartz-in-garnet barometry, or QuiG, has been applied to rocks from collisional orogenic belts in New England (Spear et al., 2014; Wolfe and Spear, 2018, in review), the Cyclades subduction complex of Greece (Ashley et al., 2014; Castro and Spear, 2016), and to high pressure gneisses and eclogites (Alvaro et al., 2020; Gonzalez et al., 2019). In several of these studies, it was found that when the results of QuiG barometry were compared to the calculated P–T conditions of the garnet isograd based on equilibrium modeling, garnet was concluded to have nucleated not near the equilibrium isograd but only after considerable overstepping .

The extent to which the above results are true — that garnet only nucleates after considerable overstepping of the equilibrium phase boundary — raises a very important question as to the validity of previously published P–T paths that are based on the assumption of garnet having grown through a sequence of equilibrium states. If application of QuiG barometry reveals a P–T history for garnet growth that deviates from the P–T histories inferred from conventional chemical zoning analysis, then the nature and causes of this deviation urgently requires evaluation because of the numerous P–T paths that have been inferred from this latter approach. Inasmuch as metamorphic P–T paths, many of which have been constrained from garnet zoning studies, are a key constraint in the tectonic interpretation or orogenesis, evaluation of the accuracy of QuiG is of considerable importance.

The purpose of this paper is to summarize some of the recent results in which QuiG barometry has been applied to constrain metamorphic P–T conditions. In some of these studies, QuiG barometry reveals a different P–T history from that previously published whereas in others QuiG barometry reinforces published inferences about metamorphic P–T paths. Examples are also presented in which QuiG barometry provides new constraints on metamorphic histories in rocks for which traditional, equilibrium studies have seen only limited success.

Section snippets

Petrologic methods applied to garnet

Garnet has long enjoyed the status of an “uncommonly useful” mineral in petrologic studies (Baxter et al., 2013) and numerous approaches have been championed to deciphering P–T histories from garnet phase equilibria and compositional zoning (e.g. see Caddick and Kohn, 2013, for an excellent review). The earliest geothermometers utilized element partitioning between garnet and coexisting phases as a proxy for temperature of crystallization or equilibration (Kretz, 1959; Ramberg, 1952) and one of

The QuiG method

Inclusion barometry is based on the assumption that when a porphyroblast such as garnet overgrows another phase so that the phase becomes included in the host, the host-inclusion pair are initially in mechanical equilibrium such that the inclusion exactly fits inside of the hole in the host. Subsequent changes in P and T and eventual exhumation results in changes in the size of the host hole and inclusion so that at surface conditions the inclusion no longer fits perfectly inside the host and

Results

This section discusses the results of reevaluating P–T paths based on application of QuiG barometry to previously studied samples. Some examples reveal close consistencies between the originally published P–T path and the results of inclusion barometry whereas other examples show no consistency between the QuiG results and garnet zoning paths. Finally, two examples are discussed in which QuiG barometry places new constraints on the evolution of samples that were not previously accessible by

Discussion

The above examples demonstrate that in some cases QuiG barometry provides constraints on metamorphic P–T conditions that are consistent with previous estimates based on classical thermobarometry and thermodynamic calculations but in other cases there are major discrepancies. The question of whether QuiG, and other inclusion barometry, is providing an accurate assessment of the conditions of garnet nucleation and growth, is of considerable importance because many studies have utilized garnet

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

The authors are greatly indebted to Jennifer Thomson for providing the sample of coticule for this study. This work was supported by NSF grants 1447468 and 1750674 to Spear and The Edward P. Hamilton Distinguished Professor of Science Education chair.

References (73)

  • R.J. Angel et al.

    EosFit-Pinc: a simple GUI for host-inclusion elastic barometry

    Am. Mineral.

    (2017)
  • R.J. Angel et al.

    Stress, strain and Raman shifts

    Zeitschrift für Kristallographie-Crystalline Materials

    (2019)
  • M. Aoya

    Extending the applicability of the Raman carbonaceous-material geothermometer using data from contact metamorphic rocks

    J. Metam. Geol.

    (2010)
  • K.T. Ashley et al.

    Geothermobarometric history of subduction recorded by quartz inclusions in garnet

    Geochem. Geophys. Geosyst.

    (2014)
  • E.F. Baxter et al.

    Garnet geochronology: Timekeeper of tectonometamorphic processes

    Elements

    (2013)
  • E.F. Baxter et al.

    Garnet: common mineral, uncommonly useful

    Elements

    (2013)
  • R.G. Berman

    Internally-consistent thermodynamic data for minerals in the system Na2O - K2O - CaO - MgO - FeO - Fe2O3 - Al2O3 - SiO2 - TiO2 - H2O - CO2

    J. Petrol.

    (1988)
  • O. Beyssac et al.

    Raman spectra of carbonaceous material in metasediments: a new geothermometer

    J. Metam. Geol.

    (2002)
  • M. Bonazzi et al.

    Assessment of the reliability of elastic geobarometry with quartz inclusions

    Lithos

    (2019)
  • M.J. Caddick et al.

    Garnet: Witness to the evolution of destructive plate boundaries

    Elements

    (2013)
  • C. de Capitani et al.

    The computation of equilibrium assemblage diagrams with Theriak/Domino software

    Am. Mineral.

    (2010)
  • A.E. Castro et al.

    Reaction overstepping and reevaluation of the peak P-T conditions of the blueschist unit Sifnos, Greece: Implications for the Cyclades subduction zone

    Int. Geol. Rev.

    (2016)
  • J.N. Christensen et al.

    Rates of tectonometamorphic processes from rubidium and strontium isotopes in garnet

    Science

    (1989)
  • J.A.D. Connolly et al.

    An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

    J. Metam. Geol.

    (2002)
  • B. Dragovic et al.

    Fluid production rate during the regional metamorphism of a pelitic schist

    Contrib. Mineral. Petrol.

    (2018)
  • P.C. England et al.

    The influence of erosion upon the mineral facies of rocks from different metamorphic environments

    J. Geol. Soc. Lond.

    (1977)
  • P.C. England et al.

    Pressure - temperature - time paths of regional metamorphism, part I: Heat transfer during the evolution of regions of thickened continental crust

    J. Petrol.

    (1984)
  • J.D. Eusden

    Stratigraphy, structure, and metamorphism across the “dorsal zone”, Central New Hampshire

  • J.D. Eusden et al.

    Silurian and Devonian rocks in the Alton and Berwick quadrangles, New Hampshire and Maine

  • J.M. Ferry et al.

    Experimental calibration of the partitioning of Fe and Mg between biotite and garnet

    Contrib. Mineral. Petrol.

    (1978)
  • P.J. Gonzalez et al.

    Quartz-in-garnet and Ti-in-quartz (QuiG-TiQ) thermobarometry: methodology and first application to a quartzofeldspathic gneiss from eastern Papua New Guinea

    J. Metam. Geol.

    (2019)
  • C.V. Guidotti

    Metamorphic petrology, mineralogy and polymetamorphism of a portion of NW Maine

  • C.V. Guidotti

    The classic high-T – low P metamorphism of west Central Maine. Is it post-tectonic or syntectonic? Evidence from porphyroblast-matrix relations: discussion

    Can. Mineral.

    (2000)
  • W.E. Hames et al.

    Fluid-assisted modification of garnet composition along rims, cracks, and mineral inclusion boundaries in samples of amphibolite facies schists

    Am. Mineral.

    (1993)
  • M. Hanfland et al.

    Graphite under pressure: equation of state and first-order Raman modes

    Phys. Rev. B

    (1989)
  • K.V. Hodges et al.

    Error estimation and empirical geothermobarometry for pelitic systems

    Am. Mineral.

    (1985)
  • Cited by (16)

    • The sedimentary record of ultrahigh-pressure metamorphism: a perspective review

      2022, Earth-Science Reviews
      Citation Excerpt :

      Such elastic approaches require (i) knowledge of the essential host–inclusion parameters, which includes the equation of state parameters, the elastic tensor, and the phonon-mode Grüneisen tensor (Angel et al., 2014, 2019); and (ii) determination of inclusion strains resulting from different thermoelastic properties of host–inclusion pairs, calculation of the mean stress (i.e. inclusion pressure), and calculation of the isomeke, that is a curve in P–T space along which both the host and inclusion have the same fractional volume change and therefore are in mechanical equilibrium (Angel et al., 2015, 2017, 2019; Murri et al., 2018; Mazzucchelli et al., 2021). These elastic models are successfully developed for elastically anisotropic mineral inclusions, such as quartz and zircon, entrapped in a quasi-elastically isotropic host minerals, such as garnet and diamond, and applied in several studies and in combination with other thermobarometric methods (e.g. Gonzalez et al., 2019, 2020; Zhong et al., 2019; Alvaro et al., 2020; Spear and Wolfe, 2020). The field of elastic thermobarometry is developing rapidly, and a recently introduced theoretical model for anisotropic mineral inclusions in anisotropic host minerals by Gonzalez et al. (2021) gives cause for optimism to study more complex host–inclusion systems in future.

    • Prolonged metamorphism of garnet-orthoamphibole gneisses from the Fuyun area: New insights into metamorphic evolution of the southern Chinese Altai orogen

      2021, Lithos
      Citation Excerpt :

      The operation of reaction (1) implies that the rocks may have experienced a marked decompression subsequently to the M1 stage. Garnet, when containing large inclusions (i.e. intergrowths of Ged + Pl + Bt), usually displays radial cracking (Fig. 5b), which was likely caused by overpressure of inclusions during decompression (Spear and Wolfe, 2020). M2 assemblages are interpreted to represent a granulite–facies overprint at relatively lower pressures as suggested by the formation of Spl + Crd after Sil, symplectitic intergrowths of Opx + Crd surrounding Grt, and the development of radial cracking in garnet while it contains large composite inclusions.

    View all citing articles on Scopus
    View full text