skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: The Role of Crustal Contamination throughout the 1329–2005 CE Eruptive Record of Mt. Etna Volcano, Italy
Abstract 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 to 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.  more » « less
Award ID(s):
2151038
PAR ID:
10545691
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford Univeristy Press
Date Published:
Journal Name:
Journal of Petrology
Volume:
65
Issue:
4
ISSN:
0022-3530
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Geological Society of America)
    Several times in the past 60 ka, Mount Etna has erupted lavas with variable alkaline character. The most recent chemical excursion began in 1971, accompanied by an increase in explosivity and eruption frequency. The origin of the alkaline signature remains enigmatic, with endmember hypotheses involving dominant contributions from mantle vs. crust. For lavas that erupted between 1329 and 2016, we used thermodynamic modeling to test if post-1971 anomalous alkalinity is dominated by mantle processes. First, we assessed mantle melting conditions required to reproduce the chemistry of potential parental magmas. Second, we examined the differentiation of partial melts as they underwent closed-system crustal storage and ascent. The mantle melting conditions explored via ~300 models include source compositions (peridotite + pyroxenite ± metasomatic phases), pressure, extent of melting, and fO2. Best-fit models that reproduce the major element chemistry of volcanics interpreted as parental magmas involve 20 to 30% melting of a peridotite-pyroxenite mantle source that contained phlogopite, pargasite, and CO2 (~1 wt.%) between ~1-1.5 GPa (~30-45 km) along the QFM+1 buffer. Subsequent isentropic decompression (adiabatic and reversible) of mantle partial melts + crystallization at shallower pressures (0.8-0.2 GPa) were modeled to test the effects of closed system ascent and storage. Isentropic decompression models yield no crystallization although temperature decreases ~3°C/0.1 GPa. Decompression + closed system crystallization fail to replicate post-1971 glass samples and do not explain observed post-1971 alkali enrichments. We conclude that partial melting of a metasomatized source produced Etna parental magmas, but closed system crustal ascent and storage cannot fully account for alkali enrichment highlighted in post-1971 products at Etna. Open system modeling suggests that assimilation and crystallization (e.g., Takach et al. 2024) play a critical role, and ongoing modeling is testing the contributions of recharge (magma replenishment) and entrainment of previously formed mushes to the high alkalinity excursions. 
    more » « less
  2. ; ; ; (, AGU Fall Meeting 2024, held in Washington, D.C., 9-13 December 2024, Session: Volcanology, Geochemistry and Petrology / Experiment to Model to Interpretation: Insights into the Impact of Experimental Design on Models and Interpretation of Natural Systems Poster, Poster No. 3209, id. V11D-3209.)
    Documenting the processes and timescales of magma formation and diversification and defining the locations, shapes, volumes, and phase states of magma storage and transport zones rely on data produced by novel analytical techniques and state-of-the-art experimental methods. Computational modeling effectively links these critical tools. The Magma Chamber Simulator (MCS) is an internally consistent thermodynamic open system model that uses experimental constraints from rhyolite MELTS (Gualda et al. 2012, Ghiorso & Gualda, 2015) to compute paths of open system magmas that evolve via processes including crystallization, magma mixing, cumulate/mush entrainment, and host-rock assimilation. MCS results yield elemental, isotope, mass, and thermal characteristics of melt ± crystals ± volatiles in "resident" magma, crustal wallrock (melt and solids), recharge magma, and entrained material. To model the petrochemical evolution of igneous rocks related by open-system processes, one typically runs 200+ models that vary initial compositions, pressures, and temperatures of magma, wallrock, etc. Comparison of model results with whole rock, mineral, and melt inclusion chemical data and other constraints (e.g., thermobarometry) yield interpretations about igneous processes at a range of scales—from how crust forms and evolves to processes responsible for in situ geochemical records of crystals—and allows assessment of epistemic and aleatoric uncertainties. Two examples of computational studies will illustrate MCS's utility and flexibility. (1) Modeling of historical basalts at Mt. Etna (Italy) provides evidence for variable degrees of melting of metasomatized mantle, followed by magma recharge and assimilation of partial melts of carbonate-flysch crust. (2) MCS models reproduce whole rock and mineral data of plagioclase-rich basalts at Steens Mountain (USA) through entrainment of gabbroic mush that likely formed in early stages of Columbia River Basalt magmatism. To enhance understanding of trans-lithospheric magma systems, future work on MCS will prioritize (i) building a post-processing environment that utilizes select statistical methods to inform "best-fit" models and to quantitatively assess uncertainty, and (ii) increasing modeling efficiency by adding automated modeling capabilities. 
    more » « less
  3. ; ; ; ; ; (, Journal of Petrology)
    null (Ed.)
    Abstract The Okataina Volcanic Centre (OVC), located in the Taupo Volcanic Zone, New Zealand, is a dominantly rhyolitic magmatic system in an arc setting, where eruptions are thought to be driven by mafic recharge. Here, Sr–Pb isotopes, and compositional and textural variations in plagioclase phenocrysts from 10 rhyolitic deposits (two caldera, one immediately post-caldera, four intra-caldera, and three extra-caldera) are used to investigate the OVC magmatic system and identify the sources and assimilants within this diverse mush zone. Plagioclase interiors exhibit normal and reverse zoning, and are commonly in disequilibrium with their accompanying glass, melt inclusions, and whole-rock compositions. This indicates that the crystals nucleated in melts that differed from their carrier magma. In contrast, the outermost rims of crystals exhibit normal zoning that is compositionally consistent with growth in cooling and fractionating melts just prior to eruption. At the intra-crystal scale, the total suite of 87Sr/86Sr ratios are highly variable (0·7042–0·7065 ± 0·0004 average 2SE); however, the majority (95 %) of the crystals are internally homogeneous within error. At whole-crystal scale (where better precision is obtained), 87Sr/86Sr ratios are much more homogeneous (0·70512–0·70543 ± 0·00001 average 2SE) and overlap with their host whole-rock Sr isotopic ratios. Whole-crystal Pb isotopic ratios also largely overlap with whole-rock Pb ratios. The plagioclase and whole-rock isotopic compositions indicate significant crustal assimilation (≥20 %) of Torlesse-like metasediments (local basement rock) by a depleted mid-ocean ridge mantle magma source, and Pb isotopes require variable fluid-dominant subduction flux. The new data support previous petrogenetic models for OVC magmas that require crystal growth in compositionally and thermally distinct magmas within a complex of disconnected melt-and-mush reservoirs. These reservoirs were rejuvenated by underplating basaltic magmas that serve as an eruption trigger. However, the outermost rims of the plagioclase imply that interaction between silicic melts and eruption-triggering mafic influx is largely limited to heat and volatile transfer, and results in rapid mobilization and syn-eruption mixing of rhyolitic melts. Finally, relatively uniform isotopic compositions of plagioclase indicate balanced contributions from the crust and mantle over the lifespan of the OVC magmatic system. 
    more » « less
  4. ; ; (, Minerals)
    Costa, Simone; Caracciolo, Alberto (Ed.)
    The Tungnárhraun basalts in southern Iceland record a transcrustal magma system formed during Holocene deglaciation. These large-volume (>1 km3) Early through Mid-Holocene lavas contain ubiquitous plagioclase feldspar macrocrysts that are too primitive to have grown from the host lavas. Thermobarometry based on plagioclase melt and clinopyroxene melt equilibrium reveals a transcrustal structure with at least three distinct storage regions. A lower-crustal mush zone at ~14–30 km is fed by primitive, low 87Sr/86Sr magmas with diverse Ti/K and Al/Ti signatures. Plagioclase feldspar growth is controlled by an experimentally determined pseudoazeotrope where crystals develop inversely correlated An and Mg contents. The rapid ascent of magmas to mid-crustal levels (~8–9 km) allows the feldspar system to revert to conventional thermodynamic phase constraints. Continued plagioclase growth releases heat, causing olivine and pyroxene to be resorbed and giving the magmas their characteristic high CaO/Al2O3 values (~0.8–1.0) and Sc contents (~52 ppm in matrix material). Mid-Holocene MgO-rich lavas with abundant plagioclase feldspar macrocrysts erupted directly from this depth, but both older and younger magmas ascended to a shallow-crustal storage chamber (~5 km) where they crystallized olivine, clinopyroxene, and plagioclase feldspar and evolved to lower MgO contents. The Sr isotope differences between the plagioclase macrocrysts and their carrier melts suggest that the fractionation involves the minor assimilation of country rock. This model does not require the physical disruption of an established and long-lived gabbroic cumulate mush. The transcrustal structures documented here existed in south Iceland at least throughout the Holocene and likely influenced much of Icelandic magmatism. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; (, American Mineralogist)
    High-Sr/Y granitoids in continental settings are sometimes erroneously regarded as the products derived from partial melting of thickened/delaminated mafic lower curst under relatively higher pressures (1.5 GPa) in a collisional orogenic setting. In fact, multiple magmatic processes in the trans-crustal magma system, such as recycling of antecrysts, crustal assimilation, and fractional crystallization, can create or modify the primary “adakitic” signature. As a result, the generation of adakitic magmas in continental settings remains controversial from a bulk-rock perspective. Here, we address the origin of adakitic plutonic rocks through geochemical and textural characterization of rock-forming minerals in the pyroxene-bearing Zhuyuan granodiorite, West Qinling, China. The Zhuyuan granodiorite formed in a post-collisional setting and primarily consists of resorbed orthopyroxene, three types of clinopyroxene, amphibole, two types of plagioclases, K-feldspar, biotite, and quartz. Type-1 Cpx has high XMg (70.0–81.7). Type-2 Cpx displays normal zoning and decreasing XMg (80.9 to 71.5) from the core to rim. Type-3 Cpx is reversely zoned, where the rims have higher XMg (75.5–86.9), Ni, Cr, suggesting a recharge event. Orthopyroxene has high-Ni and -Cr contents, as well as high XMg (80.9–82.8), indicative of antecrysts that grew in mafic magma reservoirs. The injection of magmas from different sources is supported by sieve-textured plagioclase and crystal size distributions of non-poikilitic amphibole. Finally, non-sieve textured plagioclase, biotite, K-feldspar, and quartz are late-crystallized phases, indicative of an orthocrystic origin. The melts in equilibrium with these orthocrysts display significantly higher Sr/Y values than the magma batches that crystallized other mafic phases (i.e., amphibole, clinopyroxene, and orthopyroxene). Thus, we propose that the system involved an initial high-Sr/Y melts in equilibrium with the orthocryst assemblage was generated by water-fluxed melting of intermediate to felsic sources. The addition of low Sr/Y non-orthocrysts (e.g., amphibole and pyroxene) and associated melt diluted the original “adakitic signal” in the magma reservoir and drove the bulk composition to more mafic values. Consequently, the Zhuyuan pyroxene-bearing granodiorite represents a mixture of crystals with diverse origins and distinct magma batches of various compositions (from felsic to mafic compositions). Our study emphasizes that the origin of adakitic granitoids cannot be clearly deciphered without geochemical analysis of the constituent minerals. We also suggest that Sr/Y values in plutons should be cautiously used in paleo-crustal thickness estimates in collisional settings because of possible open system scenarios as described here. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email