Implications of overstepping of garnet nucleation for geothermometry, geobarometry and P–T path calculations

Abstract

The composition of a porphyroblast such as garnet is not governed by equilibrium relations when it nucleates and grows after considerable overstepping. The approach followed here is to assume the composition is controlled by the maximum driving force (MDF) or parallel tangent model and to examine the potential pitfalls for assuming equilibrium. The MDF model ensures Fe-Mg partitioning equilibrium between the matrix phases and growing porphyroblast, but the net transfer equilibrium controlling grossular content has a different stoichiometry from the equilibrium relations. An important result of this study is that garnet zoning profiles generated using this model and assuming chemical fractionation are nearly identical to zoning profiles grown assuming continuous equilibrium and it is impossible to discern whether garnet grew under near-equilibrium conditions or only after considerable overstepping by examination of the zoning profile alone.Pressure–temperature paths calculated using the method of intersecting isopleths from simulations in which garnet was grown after overstepping and assuming the MDF model typically display nearly isothermal loading or loading with minor heating paths, even when the garnet was grown under isothermal and isobaric conditions. Garnet-plagioclase barometry on the garnet core and rim also suggests loading during garnet growth, even for crystals grown isothermally and isobarically. Garnet-biotite Fe-Mg exchange thermometry does, however, recover the correct temperature to within 10 °C for both the garnet core and rim.Comparison of natural samples to the theoretical calculations suggests that the local effective bulk composition in the vicinity of garnet is depleted in Na2O, in addition to FeO, MnO and CaO, thus rendering calculation schemes that rely on knowledge of the bulk composition difficult to administer. To the extent that garnet has nucleated and grown after considerable overstepping, application of methods of extracting P–T paths from the chemical zoning in garnet will yield incorrect results if those methods assume continuous equilibration with the matrix. The major discrepancy comes from calculation of the P–T conditions of the garnet core. Simulations and natural samples suggest that under typical scenarios the garnet rim is in near equilibrium with the matrix assemblage so equilibrium approaches (thermobarometry or intersecting isopleths) should recover the peak metamorphic conditions with reasonable assuredness.