Sr Isotope Zoning in Plagioclase from Parinacota Volcano (Northern Chile): Quantifying Magma Mixing and Crustal Contamination

Abstract

We present analyses of Sr isotope zoning by microdrilling and thermal ionization mass spectrometry in plagioclase crystals from Parinacota volcano (Central Volcanic Zone, northern Chile), which were analysed for major and minor element zoning in a previous study. Although the isotopic range of the bulk-rock samples is small at this volcano (0·7067–0·7070, except for one flow of mafic andesite at 0·7061), significant variations are seen (0·70649–0·70700) within and between plagioclase crystals. A general negative correlation is observed between Sr isotope composition and Sr concentration in the liquid in equilibrium with each plagioclase zone, as calculated from chemical zoning data and partition coefficients. Additional scatter is superimposed on this general trend, indicating a decoupling between isotopic and chemical variations for Sr. In one dacite sample a detailed isotopic profile shows increasing contamination during crystal growth, except for an abrupt decrease correlated with a dissolution surface and interpreted as a recharge event. We apply energy-constrained recharge, assimilation and fractional crystallization modelling to the melt evolution recorded in the chemical and isotopic zoning in this crystal. Results suggest 20% assimilation of the local wall-rock gneiss, at high initial temperatures. The isotopic data confirm the involvement of two contrasting mafic magmas, which are sampled at flank cinder cone vents. One (Lower Ajata) has a low Sr content with high 87Sr/86Sr, the other (Upper Ajata) has a high Sr content with lower 87Sr/86Sr. In some samples from Parinacota, the isotopic composition of plagioclase crystal rims or groundmass crystals is significantly higher than that of the high 87Sr/86Sr mafic magma. In others, where chemical zoning profiles suggest that recharge was from the low 87Sr/86Sr magma, the 87Sr/86Sr of the groundmass and crystal rims is higher than expected. This indicates either additional parent magmas to the two previously identified, or further crustal assimilation, either at lower crustal depths, before crystallization of plagioclase, or just after the last recharge. Our results illustrate the complexity of magma–crust interaction beneath Parinacota, which is likely to be representative of many other Central Andean volcanoes formed on thick crust. Such complex interactions can be revealed by combined study of chemical and isotopic zoning in plagioclase (in a textural petrographic context), despite a small whole-rock isotopic range. The distinct contamination patterns of various samples suggest an important role for the geometry, location and evolution of the magma plumbing system and, in general, variations of the thermal and compositional structure of the crust underneath the volcano.