An official website of the United States government
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
MANTLE MELTING AT MOUNT ETNA: A THERMODYNAMIC MODELING APPROACH TO UNDERSTANDING THE POST-1971 ALKALI ENRICHMENT
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
- Award ID(s):
- 2151038
- PAR ID:
- 10629393
- Publisher / Repository:
- Geological Society of America
- Date Published:
- Format(s):
- Medium: X
- Location:
- Anaheim CA
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Arc magmatism is thought to be driven by peridotite melting in the mantle wedge. Yet pyroxenites are ubiquitous in the melting region beneath magmatic arcs. Because they typically have lower solidi temperatures and higher melt productivities compared to peridotite, pyroxenites likely play a significant role in magma generation. Here, we use the Zn/Fe ratios of a global database of Pliocene–Holocene primitive arc magmas to show that, as the crustal thickness of the overlying plate increases, so does the proportion of pyroxenite-derived melts relative to peridotite-derived melts. In fact, at arcs with crustal thicknesses >40 km, the majority of magmas are sourced from pyroxenite. Major and trace element geochemistry of pyroxenite melts is consistent with derivation from mafic magmas frozen in the mantle en route to the surface. We hypothesize that, as the thickness of the continental crust increases, the mantle wedge is displaced toward higher pressures and cooler temperatures, thereby lowering the extent of peridotite melting and allowing magmas sourced from the pyroxeniteveined mantle to dominate the arc budget.more » « less
-
Abstract Spongy clinopyroxene is common in most mantle-derived xenoliths and megacrysts of eclogitic and peridotitic parageneses. Its formation is commonly attributed to the partial melting of a primary clinopyroxene in response to various factors, including changes in pressure and temperature or infiltration of external melts or fluids. In order to study the mechanism of spongy clinopyroxene formation in detail, we selected six eclogitic clinopyroxene inclusions in diamonds with varying amounts of spongy clinopyroxene (from ~10 to 100%). We employed computed tomography, electron microprobe analysis, and Raman spectroscopy to study the textural characteristics, major element concentrations, and the types of volatiles present in both phases. We also used pMELTS to model the compositions of spongy clinopyroxene and associated melts produced by the melting of primary clinopyroxene over a range of pressures and temperatures. We compare these results with estimates from major element thermobarometry of the spongy clinopyroxene. We conclude that the studied spongy clinopyroxene is the solid product of partial melting that occurs upon decompression of the primary clinopyroxene within the diamond in a near-closed system. Melting of the primary clinopyroxene occurred continuously or in pulses at different depths during the diamond’s ascent to Earth’s surface and produced variable spongy clinopyroxene and melt compositions even within the same inclusion. This is possible due to relatively rapid kimberlite ascent. The degrees of melting are various and unexpectedly high for mantle melting (between <10 and 60% with an average of ~20–30%). The produced melts are highly silicic, phonolitic, and alkali-rich. pMELTS modelling shows the spongy clinopyroxene compositions can be reproduced at pressures between 0.5–2.7 GPa and temperatures of 850–1300°C, with the majority of them satisfying the P–T conditions of 1–2 GPa and 1100–1300°C. This indicates decompression melting of primary clinopyroxene at shallow upper mantle or lower crustal conditions.more » « less
-
Abstract Ignimbrite flare-ups are rare periods of intense silicic volcanism during which the pyroclastic volume and eruptive frequency is more than an order of magnitude higher than background activity. Investigating the compositional differences between flare-up and steady-state magmas provides critical constraints on the petrogenetic causes for the event and can offer unique opportunities to investigate the role of large-scale tectonic or geodynamic processes in arc magmatism. In this study, we focus on the bimodal Deschutes Formation ignimbrite flare-up of Central Oregon, which erupted unusually high volumes of pyroclastic material 6.25–5.45 Ma from a new axis of volcanism in the Cascades arc. This episode is marked by increased eruption rates and eruption of more silicic compositions relative to the Quaternary Cascade arc, which rarely erupts rhyolites. Ignimbrites are crystal-poor (<10%) dacite to rhyolites (mostly 65–77 wt.% SiO2) with anhydrous mineral assemblages and higher FeO/MgO, Y, Eu/Eu*, MREE and Zr/Sr, indicating drier magmatic evolution compared to the Quaternary arc, and are more similar to those from the rear-arc High Lava Plains (HLP) province that lies to the east. Magnetite-ilmenite oxybarometry indicates that Deschutes Formation felsic magmas tend to be hotter and more reduced (NNO-1 to NNO) than the Quaternary arc (NNO to NNO + 1.5). Rhyolite-MELTS geobarometry suggests complex storage of diverse Deschutes Formation magmas within the shallow crust (50–250 MPa), and the common co-eruption of multiple plagioclase populations, pumice compositions, and compositionally banded pumice suggest variable degrees of mixing and mingling of distinct magmas. Deschutes magmas also have low δ18Oplagioclase values that indicate partial melting and assimilation of hydrothermally altered shallow crust. Trace element systematics and rhyolite-MELTS modeling suggests that felsic pumice cannot be produced by simple fractionation of co-erupted mafic pumice or basaltic lavas, and requires a crustal melting origin, and trace elements and Pb isotopes suggest that young mafic crust may have been the primary protolith. We suggest that partial melting produced low-Si rhyolite melt (~72 wt.%) that acted as both a parent for the most evolved rhyolites, and as a mixing endmember to create the dacite to rhyodacite magmas with heterogenous plagioclase populations. Unlike the predominantly calc-alkaline basalts erupted in the Quaternary Cascade arc, Deschutes Formation primary basalts are mostly low-K tholeiites, indicative of decompression melting. These are similar to the compositions erupted during a contemporaneous pulse of low-K tholeiite volcanism across the whole HLP that reached into the Cascades rear-arc. We suggest that intra-arc extension focused decompression melts from the back-arc into the arc and that tensional stresses allowed this high flux of hot-dry-reduced basalt throughout the crustal column, causing partial melting of mafic protoliths and the production of hot-dry-reduced rhyolite melts. Depletion of incompatible elements in successive rhyolites implies progressive depletion in fertility of the protolith. Extension also allowed for the establishment of a robust hydrothermal system, and assimilation of hydrothermally-altered rocks by magmas residing in a shallow, complex storage network lead to low δ18O melts. Our findings suggest the integral role that extensional tectonics played in producing an unusual ignimbrite flare-up of hot-dry-reduced rhyolite magmas that are atypical of the Cascades arc and may be an important contributor to flare-ups at arcs worldwide.more » « less
-
We present computational modeling outcomes for bilithologic (peridotite and pyroxenite) mantle melting in divergent environments, considering equilibrium and disequilibrium porous flow melting of 0–50 % pyroxenite in thermal equilibrium with peridotite, potential temperatures of 1300 and 1400 °C, upwelling rates from 1–50 cm yr−1, maximum porosities of 0.1–2.0 %, and four compositions that span pyroxenite melting behavior. Basalt-like pyroxenites (G2) uniquely produce low (226Ra/230Th) and (231Pa/235U) with high (230Th/238U), but quantities greater than ~10 % produce anomalously thick crust, restricting their global abundance. Silica-deficient pyroxenite (M7-16 and MIX1G) melts are more moderate, but require chemical re-equilibration during transport to resemble global basalts, while hybrid lithologies (KG1) produce melts similar to those of peridotites. Uranium-series disequilibria in partial melts can also be decoupled from trace elements by radioactive decay in two-dimensional regimes. The mantle must thus contain multiple types of pyroxenite on a global scale, with melts traveling by complex networks and experiencing heterogeneous extents of chemical re-equilibration.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.1130/abs/2024AM-404497
