Zircon saturation revisited: temperature and composition effects in a variety of crustal magma types

Hydrothermal experiments in the temperature range 750–1020°C have defined the saturation behavior of zircon in crustal anatectic melts as a function of both temperature and composition. The results provide a model of zircon solubility given by: In D_Zr^zircon/melt= ?3.80?[0.85(M?1)]+12900/T where D_Zr^zircon/melt is the concentration ratio of Zr in the stoichiometric zircon to that in the melt, T is the absolute temperature, and M is the cation ratio (Na + K + 2Ca)/(Al · Si). This solubility model is based principally upon experiments at 860°, 930°, and 1020°C, but has also been confirmed at temperatures up to 1500°C for M = 1.3. The lowest temperature experiments (750° and 800°C) yielded relatively imprecise, low solubilities, but the measured values (with assigned errors) are nevertheless in agreement with the predictions of the model.For M = 1.3 (a normal peraluminous granite), these results predict zircon solubilities ranging from ～ 100 ppm dissolved Zr at 750°C to 1330 ppm at 1020°C. Thus, in view of the substantial range of bulk Zr concentrations observed in crustal granitoids (～ 50–350 ppm), it is clear that anatectic magmas can show contrasting behavior toward zircon in the source rock. Those melts containing insufficient Zr for saturation in zircon during melting can have achieved that condition only by consuming all zircon in the source. On the other hand, melts with higher Zr contents (appropriate to saturation in zircon) must be regarded as incapable of dissolving additional zircon, whether it be located in the residual rocks or as crystals entrained in the departing melt fraction. This latter possibility is particularly interesting, inasmuch as the inability of a melt to consume zircon means that critical geochemical “indicators” contained in the undissolved zircon (e.g. heavy rare earths, Hf, U, Th, and radiogenic Pb) can equilibrate with the contacting melt only by solid-state diffusion, which may be slow relative to the time scale of the melting event.