Lower crustal control in the iron isotope systematics of plutonic xenoliths from Adak Island, Central Aleutians, with implications for arc magma geochemistry

Abstract

We present bulk-rock and mineral Fe isotope data of ultramafic to mafic xenoliths and basaltic to andesitic lavas from Adagdak Volcano (Adak Island, Central Aleutians) to study the effects of early differentiation on the Fe isotopic evolution of island arc basalts and their crystallization products. The Fe isotope composition of ultramafic cumulate xenoliths increases from dunite (δ56Fe = –0.09 to –0.02 ‰) to clinopyroxenite (δ56Fe = +0.06 to +0.09 ‰), consistent with higher modal proportions of clinopyroxene (δ56Fe = –0.05 to +0.11 ‰) relative to olivine (δ56Fe = –0.10 to +0.06 ‰) in the latter. Mid-crustal cumulate amphibole gabbro and hornblendite cumulates also record heavier Fe isotope compositions (δ56Fe = +0.04 to +0.08 ‰) due to the abundance of isotopically heavy amphibole (δ56Fe = +0.07 to +0.09 ‰) and magnetite (δ56Fe = +0.11 to +0.13 ‰) in these rocks. High inter-mineral fractionations observed in spinel-olivine and spinel-clinopyroxene pairs (Δ56Fespl-ol = +0.12 to +0.28 and Δ56Fespl-cpx = +0.06 to +0.19) suggest that spinel is not recording equilibrium crystallization conditions for the ultramafic assemblages, likely due to subsolidus Fe-Mg exchange. Our data also include Fe isotope measurements of one mantle dunite (δ56Fe = +0.03 ± 0.05 ‰). Five Adagdak lavas, spanning from basalts to andesites, yield a narrow range of δ56Fe between +0.03 and +0.06 ‰. Our results highlight the potential of amphibole in driving the Fe isotope depletion trends observed in many erupted arc lavas, as amphibole hosts 28–99 % of the FeOT budget in the amphibole gabbro and hornblendite cumulates. This is also supported by single-crystal synchrotron Mössbauer spectroscopy of two amphibole grains, the first from an amphibole gabbro and the second from a hornblendite, which yield Fe3+/ΣFe ratios of 0.55 ± 0.06 and 0.58 ± 0.02, respectively. Water content and hydrogen isotope compositions determined by secondary-ion mass spectrometry from the same amphibole grains indicate partial dehydrogenation. Using Rayleigh fractionation modeling to account for oxidation during post-crystallization dehydrogenation, we calculate magmatic Fe3+/ΣFe ratios of 0.41 ± 0.04 for the amphibole gabbro and 0.30 ± 0.05 for the hornblendite. These data are then used to estimate an appropriate Fe force constant for Adagdak amphibole and quantitatively evaluate the effects of amphibole fractionation. Through a fractional crystallization model, we show how arc melts may experience periods of increasing δ56Fe during olivine-dominated fractionation, followed by decreasing δ56Fe once magnetite and amphibole saturate as cumulate phases. Notably, this dichotomy between fractionation of isotopically light versus heavy cumulate assemblages and its effects on the Fe isotope evolution of arc magmas is not captured by the Adagdak lava record, highlighting the utility of cumulates in chronicling the early isotopic evolution of magmatic systems.