A thermodynamic model for silicate melt in CaO–MgO–Al2O3–SiO2 to 50 kbar and 1800 °C

A fully thermodynamic model for mafic melt in CaO–MgO–Al2O3–SiO2 (CMAS) has been calibrated, for calculation of melting equilibria in the pressure range 0–50 kbar. It is intended as a preliminary step towards a large-system melt model, suitable for exploring melting, melt loss and crystallization processes in a wide range of natural rock compositions. Calibration was performed with attention to the model's behaviour in its compositional subsystems, as a rigorous test of model structure and parameterization. The model is consistent with the latest Holland & Powell thermodynamic data set, and can therefore be used to calculate phase relations in conjunction with the many solid-phase activity–composition models written for the data set. Model calculations successfully reproduce experimental melting reactions in CMAS spinel lherzolite and garnet lherzolite assemblages, as well as sapphirine- and kyanite-bearing assemblages, at moderate to high pressure. Thermodynamically sensitive features, such as thermal divides are also recovered. However, some changes to the model structure will be required before the model can describe the full range of mafic and ultramafic melt compositions known from experiment at low pressures