Advantages and Limitations of Combined Diffusion-Phase Equilibrium Modelling for Pressure–Temperature–Time History of Metamorphic Rocks

Abstract

ABSTRACT This paper presents and discusses the results of phase diagram (Perple_X) and diffusion modelling (CZGM, or Compositional Zoning and its Modification by diffusion) to constrain the P–T path of metamorphism. The approach is based on the best fits between the zoning profile in measured garnet and that obtained by the intersections of garnet isopleths calculated by phase diagram modelling using whole rock bulk composition. The model was applied to garnets in natural rocks of various metamorphic grades, which were formed within different geotectonic environments. To compare the sequence of compositional change during Barrovian-type metamorphism, well-studied pelitic rocks from garnet–staurolite, kyanite–sillimanite, and sillimanite-K-feldspar metamorphic zones were selected. Garnets with two-stepped core and rim profiles that were formed during two different metamorphic stages or events were used for pressure–temperature (P–T) path constraint of each stage or event. For high-grade rocks, in which the original zoning profile in garnet was severely modified, the diffusion of the initial zoning profile was quantified to estimate the timescale of the metamorphic event. These rocks include high- to ultra-high-pressure rocks, which were subjected to thermal overprinting during collisional orogenesis. The results of the application of this approach allow for deciphering the reason why the calculated profile by phase diagram modelling does not fit with that of the measured garnet from low-grade rocks, in which garnet has preserved the original compositional zoning. This includes garnets whose nucleation was shifted from the garnet-in boundary to higher temperatures and pressures, as well as garnet crystallised during different metamorphic stages or events. Finally, the P–T paths in high-grade rocks were constrained after the multicomponent diffusion in garnet was quantified, and this was used for further P–T-time path constraint of metamorphism in the rocks.