Origin of heavy Fe isotope compositions in high-silica igneous rocks: A rhyolite perspective

The origin of heavy Fe isotope compositions in high-silica (>70 wt% SiO2) igneous rocks remains a highly controversial topic. Considering that fluid exsolution in eruptive rocks is more straight-forward to constrain than in plutonic rocks, this study addresses the problem of Fe isotope fractionation in high-silica igneous rocks by measuring Fe isotope compositions of representative rhyolitic samples from the Neoproterozoic volcanic-sedimentary basins in southern China and the Triassic Tu Le Basin in northern Vietnam. The samples show remarkably varied δ56FeIRMM014 values ranging from 0.05 ± 0.05‰ to 0.55 ± 0.05‰, which is among the highest values reported from felsic rocks. The extensional tectonic setting and short melt residence time in magma chambers for the studied rhyolites rule out Soret diffusion and thermal migration processes as causes of the high δ56Fe values. Effects of volcanic degassing and fluid exsolution on bulk rock δ56Fe values for the rhyolites are also assessed using bulk rock geochemical indicators and Rayleigh fractionation models, and these processes are found to be insufficient to produce resolvable changes in Fe isotope compositions of the residual melt. The most probable mechanism accounting for heavy Fe isotope compositions in the high-silica rhyolites is narrowed down to fractional crystallization processes in the magma before rhyolite eruption. Removal of isotopically light Fe-bearing minerals (i.e. ulvöspinel-rich titanomagnetite, ilmenite and biotite) is proposed as the main cause of Fe isotope variation in silicic melts during magmatic evolution. This study implies that crystal fractionation is the dominant mechanism that controls Fe isotope fractionation in eruptive rocks and Fe isotopes could be used to study magmatic differentiation of high-silica magmas.