An automated method for the calculation of P–T paths from garnet zoning, with application to metapelitic schist from the Kootenay Arc, British Columbia, Canada

An automated method for the calculation of P–T paths based on garnet zoning is presented and used to interpret zoning in metapelitic schist from the southern Canadian Cordillera. The approach adopted to reconstruct the P–T path is to match garnet compositions along a radial transect with predictions from thermodynamic forward models, while iteratively modifying the composition to account for fractional crystallization. The method is applied to a representative sample of garnet‐ and staurolite‐bearing schist from an amphibolite facies Barrovian belt in the southern Canadian Omineca belt. Garnet zoning in these schists is concentric and largely continuous from core to rim. Three zones are present, the first two of which coincide with sector‐zoned cores of garnet crystals. Similar zoning is developed in rocks that contain or lack staurolite, respectively, suggesting garnet growth was restricted to the initial part of the prograde P–T path prior to the development of staurolite. Growth zoning in large garnet crystals has not been significantly modified by diffusion. This interpretation is based on zoning characteristics of garnet crystals and is further supported by results of a forward model incorporating the effects of simultaneous fractional crystallization and intracrystalline diffusion. The P–T path calculated for this rock includes an initial, linear stage with a high d P / d T , and a later stage dominated by heating. The approach adopted in this study may have application to other garnet‐bearing rocks in which growth zoning is preserved.