Phase equilibria modeling in igneous petrology: use of COMAGMAT model for simulating fractionation of ferro-basaltic magmas and the genesis of high-alumina basalt

Abstract

A new version of COMAGMAT-3.5 model designed for computer simulations of equilibrium and fractional crystallization of basaltic magmas at low to high pressures is presented. The most important modifications of COMAGMAT include an ability to calculate more accurately the crystallization of magnetite and ilmenite, allowing the user to study numerically the effect of oxygen fugacity on basalt magma fractionation trends. Methodological principles of the use of COMAGMAT were discussed based on its thermodynamical and empirical basis, including specific details of the model calibration. Using COMAGMAT-3.5 a set of phase equilibria calculations (called Geochemical Thermometry) has been conducted for six cumulative rocks from the Marginal Border Series of the Skaergaard intrusion. As a result, initial magma temperature (1165±10°C) and trapped melt composition proposed to be parental magma to the Skaergaard intrusion were determined. Computer simulations of perfect fractionation of this composition as well as another proposed parent produced petrochemical trends opposite to those followed from natural observations. This is interpreted as evidence for an initial Skaergaard magma containing a large amount of olivine and plagioclase crystals (about 40–45%), so that the proposed and calculated parents are related through the melt trapped in the crystal–liquid mixture. This promotes the conclusion that the Skaergaard magma fractionation process was intermediate between equilibrium and fractional crystallization. In this case the classic Wager's trend should be considered an exception rather than a rule for the differentiation of ferro-basaltic magmas. A polybaric version of COMAGMAT has been applied for the genetic interpretation of a volcanic suite from the Klyuchevskoi volcano, Kamchatka, Russia. To identify petrological processes responsible for the observed suite ranging from high-magnesia to high-alumina basalts, we used the model to simulate the Klyuchevskoi suite assuming isobaric crystallization of a parental HMB magma at a variety of pressures and a separate set of simulations assuming fractionation during continuous magma ascent from a depth of 60 km. These results indicate that the Klyuchevskoi trend can be produced by ∼40% fractionation of Ol–Aug–Sp±Opx assemblages during ascent of the parental HMB magma over the pressure range 19–7 kbar with the rate of decompression being 0.33 kbar/% crystallized (at 1350–1110°C), with ∼2 wt.% of H 2O in the initial melt and ∼3 wt.% of H 2O in the resultant high-Al basalt.