Zircon saturation in silicate melts: a new and improved model for aluminous and alkaline melts

The importance of zircon in geochemical and geochronological studies, and its presence not only in aluminous but also in alkaline rocks, prompted us to think about a new zircon saturation model that can be applied in a wide range of compositions. Therefore, we performed zircon crystallization experiments in a range of compositions and at high temperatures, extending the original zircon saturation model proposed by Watson and Harrison (Earth Planet Sci Lett 64: 295-304, 1983) and Boehnke et al. (Chem Geol 351: 324-334, 2013). We used our new data and the data from previous studies in peraluminous melts, to describe the solubility of zircon in alkaline and aluminous melts. To this effect, we devised a new compositional parameter called G[(3 center dot Al2O3 + SiO2)/(Na2O + K2O + CaO + MgO + FeO)] (molar proportions), which enables to describe the zircon saturation behaviour in a wide range of rock compositions. Furthermore, we propose a new zircon saturation model, which depends basically on temperature and melt composition, given by (with 1 sigma errors): ln [Zr] = (4.29 +/- 0.34) - (1.35 +/- 0.10) center dot ln G+ (0.0056 +/- 0.0002) center dot T(degrees C) where [Zr] is the Zr concentration of the melt in mu g/g, G is the new parameter representing melt composition and T is the temperature in degrees Celsius. The advantages of the new model are its straightforward use, with the G parameter being calculated directly from the molar proportions converted from electron microprobe measurements, the temperature calculated given in degrees Celsius and its applicability in a wider range of rocks compositions. Our results confirm the high zircon solubility in peralkaline rocks and its dependence on composition and temperature. Our new model may be applied in all intermediate to felsic melts from peraluminous to peralkaline compositions.