Modeling of zoning patterns in garnet: Thermodynamic and kinetic aspects

Results from the modeling of compositional zoning patterns in garnet porphyroblasts from the medium-grade metapelitic schist of northern Ladoga area are considered. The P-T pseudosections in the model KMnFMASH system were calculated for this purpose using THERMOCALC software (Powell et al., 1998). Particular emphasis is placed upon the effect of garnet growth kinetics on the model zoning profiles for Mn (Gulbin, 2013). They fit the observed profiles if intergranular diffusion-controlled growth is assumed for porphyroblasts. Additionally, a model of metamorphic fractional crystallization is used to characterize the oscillations in both the garnet core and rim. Starting from the assumption that a reservoir, where garnet grows, consists of chlorite, and that this mineral is intensely replaced with biotite and staurolite at the onset of crystallization, a partial release of Mn from the chlorite structure and the concentration of this component in intergranular space is inferred. In terms of the model under consideration, the coefficient of the Mn partition between garnet and reservoir temporarily increases at the early stage of garnet growth, giving rise to the enrichment of the intermediate zone of porphyroblasts in Mn. In addition to the modeling of garnet growth zoning, its subsequent diffusion modification is estimated on the basis of intracrystalline diffusion profile simulation. The reverse zoned, Mn-rich and Mg-poor garnet rims are related to retrograde growth of garnet at the late stage of porphyroblast formation. The data obtained are used to constrain metamorphic evolution and the P-T-t path of staurolite-bearing rocks in the northern domain of the studied area.