Chemical mass transfer in magmatic processes IV. A revised and internally consistent thermodynamic model for the interpolation and extrapolation of liquid-solid equilibria in magmatic systems at elevated temperatures and pressures

TL;DR: In this article, a regular solution-type thermodynamic model for twelve-component silicate liquids in the system SiO2-TiO 2-Al 2O3-Fe2O 3-Cr2O3 -FeO-MgO-CaO-Na2O-K 2O-P2O5-H2O is calibrated.

Abstract: A revised regular solution-type thermodynamic model for twelve-component silicate liquids in the system SiO2-TiO2-Al2O3-Fe2O3-Cr2O3-FeO-MgO-CaO-Na2O-K2O-P2O5-H2O is calibrated. The model is referenced to previously published standard state thermodynamic properties and is derived from a set of internally consistent thermodynamic models for solid solutions of the igneous rock forming minerals, including: (Mg,Fe2+,Ca)-olivines, (Na,Mg,Fe2+,Ca)M2 (Mg,Fe2+, Ti, Fe3+, Al)M1 (Fe3+, Al,Si)2TETO6-pyroxenes, (Na,Ca,K)-feldspars, (Mg,Fe2+) (Fe3+, Al, Cr)2O4-(Mg,Fe2+)2 TiO4 spinels and (Fe2+, Mg, Mn2+)TiO3-Fe2O3 rhombohedral oxides. The calibration utilizes over 2,500 experimentally determined compositions of silicate liquids coexisting at known temperatures, pressures and oxygen fugacities with apatite ±feldspar ±leucite ±olivine ±pyroxene ±quartz ±rhombohedral oxides ±spinel ±whitlockite ±water. The model is applicable to natural magmatic compositions (both hydrous and anhydrous), ranging from potash ankaratrites to rhyolites, over the temperature (T) range 900°–1700°C and pressures (P) up to 4 GPa. The model is implemented as a software package (MELTS) which may be used to simulate igneous processes such as (1) equilibrium or fractional crystallization, (2) isothermal, isenthalpic or isochoric assimilation, and (3) degassing of volatiles. Phase equilibria are predicted using the MELTS package by specifying bulk composition of the system and either (1) T and P, (2) enthalpy (H) and P, (3) entropy (S) and P, or (4) T and volume (V). Phase relations in systems open to oxygen are determined by directly specifying the fo2 or the T-P-fo2 (or equivalently H-P-fo2, S-P-fo2, T-V-fo2) evolution path. Calculations are performed by constrained minimization of the appropriate thermodynamic potential. Compositions and proportions of solids and liquids in the equilibrium assemblage are computed.