The Role of Crustal Contamination Throughout the 1329–2005 CE Eruptive Record of Mt. Etna Volcano, Italy

The nearly continuous volcanic eruption record at Mt. Etna dating back ~700 years provides an excellent opportunity to investigate the geochemical evolution of a highly active volcano. Of particular interest is elucidating the cause of selective enrichment in alkali elements (K and Rb) and 87Sr/86Sr observed in various episodes of past activity. More recently, this alkali enrichment trend started to manifest in the 17th century and accelerated after 1971, and was accompanied by an increase in the volume, frequency, and explosivity of eruptions. Explanations for this signature include recharge of alkali-enriched magmas and/or crustal contamination from the subvolcanic basement. This study quantitatively examines the role of crustal contamination in post-1971 Etnean magma compositions via hundreds of open-system phase equilibria and trace element calculations based upon whole-rock major oxides, trace elements, 87Sr/86Sr ratios, and mineral compositional data. Available pre-1971 petrochemical data are satisfactorily reproduced by fractional crystallization of a high whole-rock MgO (12–17 wt.%), Ni (135–285 ppm), and Cr (920–1330 ppm) parental magma composition that is documented in Etna’s ~4-ka fall-stratified deposit. Observed post-1971 whole-rock and glass trends and phase equilibria are reproduced via modeled assimilation of a skarn and flysch mixture, lithologies that represent the uppermost 10–15 km of sedimentary rocks beneath Etna. Notably, models show that K2O (wt.%) and Rb (ppm) behave incompatibly during partial melting of skarn/flysch. Additionally, the observed elevation of 87Sr/86Sr in post-1971 samples is consistent with the addition of radiogenic Sr from wallrock partial melts. In best-fit models, which yield observed post-1971 K2O, Rb, and 87Sr/86Sr trends, ~17% anatectic melt is assimilated and there may be a subordinate stoped wallrock component of ≤2% (percentage is relative to the starting mass of pristine magma). Previous work has shown that metasomatized spinel lherzolite and garnet pyroxenite can be melted in different proportions to reproduce long- and short-term changes observed in Etna’s geochemical products. We propose that the alkali enrichment signature observed after 1971 can be fully explained through the combination of mantle heterogeneity and crustal contamination. In particular, up to ~20% crustal input coupled with mantle heterogeneity of primitive melts explains the geochemical signals quite well. The influence of crustal contamination on post-1971 lavas is, in part, the result of frequent recharge of magmas that thermally primed the middle to upper crust and enhanced its partial melting.