A trace-element geochemical model for imperfect fractional crystallization associated with the development of crystal zoning

Although many petrological studies of volcanic rocks have suggested that crystallization proceeds within magma bodies, highly compatible trace elements do not display the marked variations and extreme depletions predicted to result from perfect fractional crystallization. Imperfect crystal–liquid separation is a key process in explaining this paradox. The presence of suspended crystals greatly affects variations in highly compatible elements, and has been quantitatively modeled by assuming perfect equilibrium between the suspended crystals and the liquid (equilibrium crystallization and imperfect separation; ECIS); however, volcanic rocks generally contain zoned phenocrysts that reflect the absence of solid-state equilibration. The present study develops a mass-balance model for zoned crystallization and imperfect separation (ZCIS). The ZCIS process is more efficient than the conventional ECIS process in generating depleted compatible elements. These two end-member models are able to explain the compositional range of igneous rocks that experienced imperfect fractional crystallization under natural conditions. The predicted compositional regions in bivariate trace-element diagrams successfully account for the sizes and shapes of the regions defined by whole-rock and melt-inclusion data from the Bishop Tuff, CA, USA.