Magma Chamber Simulator
The Magma Chamber Simulator (MCS) is a numerical code that quantifies the effects of concurrent and/or serial magma Recharge (up to five recharge events), Assimilation of partial melts and/or stoped blocks, and Fractional Crystallization (RnAFC) on the compositional evolution of a magma body. The MCS integrates a user-defined sequence of open system interactions defined via the MCS IGOR (exec) interface with one of four user-selected MELTS engines (v 1.2.x, v 1.1.x, v 1.0.x or pMELTS) to generate self-consistent phase equilibria (composition and abundances of all crystals, melt, and supercritical H2O-CO2 fluid, temperature, pressure and redox conditions) in all subsystems [M, WR and R (n = 0 to 5)] throughout the complex series of user-defined RnAFC interactions. In addition to intensive and extensive thermodynamic and phase equilibria data, multiple trace element concentrations and isotopic ratios may also be computed based on user-input of partition coefficients and trace element bulk compositions. A comprehensive set of instructional tutorial videos that provide step-by- step instructions for code download and implementation are available by clicking on Video Tutorials from the Home page.
Details about the theory of MCS are provided in Bohrson et al. 2014, Journal of Petrology (v. 55, #9, pp.1685-1717). Examples of research studies in which MCS results figure prominently may be found by clicking MCS studies on the Home page. The MCS can be used to generate data for comparison with natural systems of whole-rock, melt inclusion, and phase compositions (major, minor, trace and isotopes) in order to establish the most significant RAFC processes defining the petrogenesis of suites of volcanic and plutonic rocks. The MCS can be used to model AFC processes where assimilation of partial melts of wallrock (WR) or wholesale stoping of blocks (of arbitrary composition and wholly solid, mushy or liquid) is the mechanism of contamination of resident magma (M). Additionally, RFC and FC-only simulations where assimilation is precluded by user-declared fiat can also be performed. Complex scenarios with wallrock partial melt assimilation, magma recharge, assimilation by stoping and fractional crystallization in arbitrary sequence are easily handled and computed using a user-defined exec input file (MES_filename).