Determination of zircon/melt trace element partition coefficients from SIMS analysis of melt inclusions in zircon

J. B. Thomas, R. J. Bodnar, N. Shimizu, A. K. Sinha

Research output: Contribution to journalArticle

120 Scopus citations

Abstract

Partition coefficients (zircon/meltDM) for rare earth elements (REE) (La, Ce, Nd, Sm, Dy, Er and Yb) and other trace elements (Ba, Rb, B, Sr, Ti, Y and Nb) between zircon and melt have been calculated from secondary ion mass spectrometric (SIMS) analyses of zircon/melt inclusion pairs. The melt inclusion-mineral (MIM) technique shows that DREE increase in compatibility with increasing atomic number, similar to results of previous studies. However, DREE determined using the MIM technique are, in general, lower than previously reported values. Calculated DREE indicate that light REE with atomic numbers less than Sm are incompatible in zircon and become more incompatible with decreasing atomic number. This behavior is in contrast to most previously published results which indicate D > 1 and define a flat partitioning pattern for elements from La through Sm. The partition coefficients for the heavy REE determined using the MIM technique are lower than previously published results by factors of ≈ 15 to 20 but follow a similar trend. These differences are thought to reflect the effects of mineral and/or glass contaminants in samples from earlier studies which employed bulk analysis techniques. DREE determined using the MIM technique agree well with values predicted using the equations of Brice (1975), which are based on the size and elasticity of crystallographic sites. The presence of Ce4+ in the melt results in elevated DCe compared to neighboring REE due to the similar valence and size of Ce4+ and Zr4+. Predicted zircon/meltD values for Ce4+ and Ce3+ indicate that the Ce4+/Ce3+ ratios of the melt ranged from about 10-3 to 10-2. Partition coefficients for other trace elements determined in this study increase in compatibility in the order Ba < Rb < B < Sr < Ti < Y < Nb, with Ba, Rb, B and Sr showing incompatible behavior (DM < 1.0), and Ti, Y and Nb showing compatible behavior (DM > 1.0). The effect of partition coefficients on melt evolution during petrogenetic modeling was examined using partition coefficients determined in this study and compared to trends obtained using published partition coefficients. The lower DREE determined in this study result in smaller REE bulk distribution coefficients, for a given mineral assemblage, compared to those calculated using previously reported values. As an example, fractional crystallization of an assemblage composed of 35% hornblende, 64.5% plagioclase and 0.5% zircon produces a melt that becomes increasingly more enriched in Yb using the DYb from this study. Using DYb from Fujimaki (1986) results in a melt that becomes progressively depleted in Yb during crystallization.

Original languageEnglish (US)
Pages (from-to)2887-2901
Number of pages15
JournalGeochimica et Cosmochimica Acta
Volume66
Issue number16
DOIs
StatePublished - Aug 1 2002
Externally publishedYes

    Fingerprint

ASJC Scopus subject areas

  • Geochemistry and Petrology

Cite this