A KINETIC MODEL OF METAMORPHISM: AN APPLICATION TO SILICEOUS DOLOMITES

Abstract

We use a kinetic model of a metamorphic system to study the effect of competing rates of reaction, fluid injection, and heating on the evolution of the reaction pathway in temperature/composition space at constant pressure. We show that for rocks in contact with mixed volatile (e.g., CO2-H2O) fluids the reaction path may be quite different from what is expected from equilibrium-based petrologic models. Equilibrium-based models, used to understand the development of rock systems undergoing mineral reactions during a metamorphic event, rely on the Gibbs phase rule and only consider stable phases. For constant pressure, the temperature-composition paths follow univariant curves and significant reactions may occur at invariant points. By contrast, the more general kinetic treatment is not constrained by equilibrium, although with the proper competing rates equilibrium is a possible endmember of the kinetic approach. The deviation from equilibrium depends on the competing rates of reaction, heating, and fluid injection. A key element required by the kinetic approach is the inclusion of metastable reactions in the formulation, whereas such reactions are irrelevant for equilibrium-based models. Metastable reactions are often involved in a complex interplay with common prograde stable metamorphic reactions. We present model results for the well-studied CaO-MgO-SiO2-CO2-H2O (CMS) system to show how the system evolves under kinetic control. Our simulations and discussion focus on the behavior of the CMS system under a number of closed and open system conditions. Special attention is paid to closed system behavior in the vicinity of the (first) isobaric invariant point (with Dol, Qtz, Tlc, Cal, and Tr). Also, for open systems with massive fluid infiltration we consider heating rates varying from contact to regional metamorphic conditions. For some geologically reasonable rates of reactions, heating, and fluid injection, our results demonstrate that equilibrium conditions may be significantly overstepped in metamorphic systems. We used overall mineral reactions in this model with rates based on experimental results. Future models could rely on more fundamental dissolution and precipitation reactions. Such an extension would require additional kinetic rate data, as well as mineral solubilities in mixed volatile fluids.