A Model of Magmatic Crystallization

A computer model simulating fractional crystallization at oneatmosphere pressure incorporates nine broadly-defined minerals—magnetite,olivine, hypersthene, augite, quartz, plagioclase, orthoclase,leucite, and nepheline. The crystallization temperature of eachmineral is considered to be a smooth function of the compositionof the magmatic liquid. These mineral temperature equationsare obtained by multiple linear regression analysis of informationfrom published silicate systems and rock melting experiments.The nine equations are solved for any primary liquid, withinthe broad range of common magma types, to select the crystallizingmineral or minerals. Partition ratios from published experimentsand analyses of lavas and phenocrysts permit calculation ofthe composition of the crystallizing mineral assemblage. Subtractionof a small amount of that composition from the primary liquidyields a new liquid, which may be recycled to yield a sequenceof liquids during fractional crystallization. The crystallizationmodel handles assemblages of co-precipitating minerals, andcan trace progressive saturation in new minerals, substitutionof a new mineral for an old mineral, and cessation of crystallizationof a mineral. The sequences of minerals and liquids derivedfrom a broad set of primary liquids are geologically realistic,so the model is useful in predicting phenocrysts in volcanicrocks and events during crystallization of shallow intrusions.