Roberto Moretti

Nicole Di Métrich

Ilenia Arienzo

Valeria Di Renzo

Alessandro Aiuppa

Patrick Allard

Part I: Volatile stocking, gas fluxing, and the shift from low-energy to highly explosive basaltic eruptions

Keywords: Basaltic explosive volcanism, Mt. Etna, chemical thermodynamics, melt inclusions, total volatiles

HAL is

A C C E P T E D

M A N U S C R I P T

Introduction

Explosive volcanism is generally related to a large amount of volatiles, i.e., a high gas/melt mass ratio in magmatic system, which strongly determines the regime under which the magma ascends to surface in open conduits (e.g., Houghton and Gonnermann, 2008). Most commonly, explosive basaltic volcanism occurs in arc environments, where the mantle source is enriched in fluids derived from subducted slabs (Métrich and Wallace, 2008). Improved understanding of how volatiles are released from a basaltic volcano experiencing a range of eruptive styles can provide key insights into physico-chemical changes occurring in the plumbing system when activity evolves from effusive, or mild-explosive, to high-explosive, or even Plinian, eruptions. For this purpose, a scan of volatile contents dissolved in melt inclusions (MIs) is mandatory to provide a petrological imaging of plumbing systems (e.g., [START_REF] Marianelli | Shallow and deep reservoirs involved in magma supply of the 1944 eruption of Vesuvius[END_REF][START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF]Spilliaert et al., 2006a;[START_REF] Mangiacapra | The deep magmatic system of the Campi Flegrei caldera (Italy)[END_REF]Moretti et al., 2013b), and to quantify total volatiles and their variations due to fluid transfer (essentially CO 2 ) throughout the system (e.g., Barsanti et al., 2009).

Mt. Etna basaltic volcano, located in Sicily (southern Italy), shows such a wide eruptive spectrum, related to the emission of volatile-rich alkali-basalts to hawaiitic magmas. In this work, we studied melt inclusions (MIs) hosted in olivine crystals from products extruded during: i) Bocca

Nuova (BN) Strombolian activity in 1997-98, ii) the paroxysmal event of February 1999 from the South-East crater (SEC) (La [START_REF] Delfa | Activity of Mount Etna preceding the February 1999 fissure eruption: inferred mechanism from seismological and geochemical data[END_REF], and iii) the 122 BC Plinian eruption (Coltelli et al. 1998;2000). They have been analyzed in terms of major constituents, sulfur and chlorine concentrations and, for a subset of them, dissolved H 2 O and CO 2 contents and sulfur speciation. The relationships with the sulfur redox state are treated in a companion paper (Moretti et al., submitted).

The analysis of total (exsolved + dissolved) volatiles (H 2 O TOT and particularly CO 2

TOT

) allows us to infer gas/melt variations in the system and characterize volatile abundances of the deep primitive magma source. This allows obtaining precious information on the evolution of degassing patterns of other volatiles components, such as S and Cl. Because such an evolution is greatly controlled by processes of CO 2 addition or loss, already occurring in the deep source, the results are used to

A C C E P T E D M A N U S C R I P T

discriminate -and switch among-the different degassing behaviors characterizing mildly and highly explosive, up to Plinian, events. The proposed analysis and modeling, here developed for conditions relevant to Mt. Etna, can be extended to other basaltic explosive volcanoes worldwide.

Geological background

Volcanological and geochemical outline

Mt. Etna, located at the boundary of the colliding African and European plates, is the largest and most active volcano of Europe. Its buildup is related to the rollback of the slab subducting under the Aeolian Islands [START_REF] Doglioni | A geodynamic model of the southern Apennines accretionary prism[END_REF]). In the general classification of Mediterranean volcanism, Mt. Etna is considered as an example of anorogenic magmatism with hot-spot affinity [START_REF] A C C E P T E D M A N U S C R I P T Peccerillo | Plio-Quaternary volcanism in Italy[END_REF]. Mt. Etna alkali basalts to hawaiites are characterized by high contents of dissolved volatiles [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF]Spilliaert et al. 2006aSpilliaert et al. , 2006b)), as well as huge degassing rates at the surface, both as diffuse soil degassing from the volcano flanks and gas emissions from the summit crater plumes and fumarolic vents [START_REF] Allard | Eruptive and diffuse emissions of CO 2 from Mount Etna[END_REF][START_REF] Aiuppa | Rates of carbon dioxide plume degassing from Mount Etna volcano[END_REF][START_REF] Giammanco | Soil and fumarole gases of Mount Etna: geochemistry and relations with volcanic activity[END_REF]. Mt. Etna volcanic gases show compositions unusually enriched in carbon and sulfur (and poor in water) relative to volcanoes in convergent settings, and have compositions intermediate between arc and hot-spot and divergent-plate volcanoes [START_REF] Symonds | Volcanic-gas studies: methods, results, and applications[END_REF]. This peculiar gas signature suggests that the emitted alkaline magmas originate from a C-S-rich source, either of deep (e.g., a metasomatized mantle, likely related to the de-volatilization of the Ionian subducting slab; [START_REF] Tonarini | Subduction-like fluids in the genesis of Mt. Etna magmas: evidence from boron isotopes and fluid mobile elements[END_REF] or shallow (crustal) provenance (see for instance Chiodini et al. 2011).

The volcanic history in the Etnean area probably began at about 500 ka, when both submarine and subaerial tholeiitic lavas were erupted in the so-called «Pre-etnean Gulf» (Condomines et al. 1982;[START_REF] Gillot | The evolution of Mount Etna in the light of potassium-argon dating[END_REF]. Since then, the erupted lavas gradually changed composition from subalkaline towards Na-alkaline. Mt. Etna's historic and recent (up to 1970) activity principally consisted of continuous degassing from the summit craters Bocca Nuova (BN), North-East Crater (NEC), South-

ACCEPTED MANUSCRIPT

A C C E P T E D M A N U S C R I P T

East Crater (SEC), and Voragine (VOR) and frequent lava effusions and Strombolian eruptions of relatively low energy, both from the summit craters and from lateral vents that are chiefly concentrated in three rift zones [START_REF] Kieffer | Évolution structrurale et dynamique d'un grand volcan polygénique: stades d'édification et activité actuelle de l'Etna[END_REF][START_REF] Kieffer | Évolution structrurale et dynamique d'un grand volcan polygénique: stades d'édification et activité actuelle de l'Etna[END_REF][START_REF] Kieffer | Évolution structrurale et dynamique d'un grand volcan polygénique: stades d'édification et activité actuelle de l'Etna[END_REF].

Starting from the early 1970s, Mt. Etna has shown an escalation in eruptive frequency and magnitude (Allard et al. 2006a), associated with an increasingly alkaline and radiogenic signature of erupted products (e.g. Condomines et al. 1995;[START_REF] Métrich | Crystallization Driven by Decompression and Water Loss at Stromboli Volcano (Aeolian Islands, Italy)[END_REF]Ferlito and Lanzafame 2010 and references therein). Such an evolution tracks the inflow of a new magma type, associated with infiltration of CO 2 -rich deep gases in the upper magma storage regions (Spilliaert et al. 2006a;[START_REF] Ferlito | The role of supercritical fluids in the potassium enrichment of magmas at Mount Etna volcano (Italy)[END_REF]. This enhanced activity has consisted of very abundant gas emissions, persistent Strombolian activity, series of fire-fountains and several fissural eruptions with high effusion rate (e.g. 1991-93; 2001; 2002-2003). In particular, the last twenty years (i.e. after the 1991-93 flank eruption) have been characterized by several series of lava fountain episodes in late 1995 at the NEC and between 1998 and 2001 at SEC. The latter also include the February 1999 eruption, which was fed by a fissure on the SSE flank of the SEC (La [START_REF] Delfa | Activity of Mount Etna preceding the February 1999 fissure eruption: inferred mechanism from seismological and geochemical data[END_REF]. In particular, the SEC activity shows a periodicity of pre-eruptive magmatic processes, since its important role in the plumbing system (see [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF]Spilliaert et al., 2006a for details), likely volatile-and Mg-rich [START_REF] Kamenetsky | Arrival of extremely volatile-rich high-Mg magmas changes explosivity of Mount Etna[END_REF], is thought to have progressively replaced the trachybasaltic magma formerly residing in the volcano's upper plumbing system, and was first erupted "undiluted" during the 2001 eruption [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF]). At Mt. Etna, highenergy Plinian eruptions occurred in 3890 BP and 122 BC (Coltelli et al. 1998;2000). These events, at the same volcano, emphasize the critical role played by the exsolution, separation and expansion of the gas phase from the volatile-rich primitive Etnean magmas, and clearly show that at Mt. Etna fractional crystallization is only one of the processes contributing to the chemistry of magmas. In addition, magma mixing processes can explain both chemical and isotopic variations of the erupted products, as shown for the 1995-2005 activity at Mt. Etna summit craters (Corsaro et al. 2013;[START_REF] Kahl | Constraints on the nature and evolution of the magma plumbing system of Mt. Etna volcano (1991-2008) from a combined thermodynamic and kinetic modelling of the compositional record of minerals[END_REF], in which a radiogenic, more primitive but alkali-rich component can be recognized, the one marking the new alkali-rich trend observed at Etna in the last decades. Corsaro et al. (2013) defined the most primitive and radiogenic magma that erupted from SEC during the 17 May 2000 fire fountain, and an evolved, less radiogenic magma emitted on 6 June 1996 during Strombolian activity at NEC, as compositional endmembers 1 and 2 respectively. By choosing such endmembers, mixing proportions between end members 2 and 1 can be estimated to vary between 1:0.25 and 1:3 (Corsaro et al. 2013).

Volatiles ( H 2 O, CO 2 , S and Cl) in Mt. Etna magmas: state of the art.

Several studies based on MI chemistry and volatile contents for past and present Mt. Etna eruptions have been carried out to reconstruct structure and evolution of its plumbing system (e.g., [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF]Spilliaert et al. 2006aSpilliaert et al. , 2006b;;Corsaro et al. 2009;[START_REF] Ferlito | The role of supercritical fluids in the potassium enrichment of magmas at Mount Etna volcano (Italy)[END_REF] Figure 1). Generally, H 2 O and CO 2 determinations are fundamental to defining a) the total (dissolved plus exsolved) volatile content of the parental melt [START_REF] Papale | Determination of total H 2 O and CO 2 budgets in evolving magmas from melt inclusion data[END_REF], and b) the depth intervals of reservoirs in which magma evolution takes place (e.g. Moretti et al., 2013a,b;Spilliaert et al., 2006a). On the other hand, S content and speciation state can be used to understand the

A C C E P T E D M A N U S C R I P T

oxidation state controlling the evolving magma, whereas Cl represents an important tracer to understand the interplay of crystallization and degassing (e,g, [START_REF] Lesne | Experimental simulation of closed-system degassing in the system basalt-H 2 O-CO 2 -S-Cl[END_REF]Spilliaert et al., 2006b) and later stages of volcanic degassing (e.g., [START_REF] Alletti | Chlorine partitioning between a basaltic melt and H 2 O-CO 2 fluids at Mount Etna[END_REF].

Contents of dissolved H 2 O and CO 2 up to 3.6 wt% and ~ 3200 ppm, respectively, have been measured in Etna MIs (e.g., [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF]Spilliaert et al. 2006aSpilliaert et al. ,2006b)). The very different solubility behavior of these two volatiles indicates that Etnean magmas coexist at depth with a CO 2dominated fluid phase. [START_REF] Clocchiatti | Coexistence de «verres» alcalins et tholeiiliques saturés en CO 2 dans les olivines des hyaloclastites d'Aci Castello (Etna Sicile, Italie). Arguments en faveur d'un manteau anormal et d'un reservoir profond[END_REF] 1975[START_REF] Dixon | An experimental study of water and carbon dioxide solubilities in mid-ocean ridge basaltic liquids. part I: Calibration and solubility models[END_REF][START_REF] Allard | Eruptive and diffuse emissions of CO 2 from Mount Etna[END_REF]. CO 2 volcanic flux has dropped to 2 ktons/day since 2001 [START_REF] Aiuppa | Rates of carbon dioxide plume degassing from Mount Etna volcano[END_REF]), but rapidly increases upon replenishment of the shallow magmatic reservoir, yielding values up to ~ 30 times larger, such as in the case of the 2004-2005 effusive eruption [START_REF] Aiuppa | Rates of carbon dioxide plume degassing from Mount Etna volcano[END_REF].

High total (exsolved + dissolved) volatile contents have been assessed by [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF] and Spilliaert et al. (2006a), implying from 1.5 to 8 wt% of CO 2 -dominated free-fluid phase in the parental deep magma (calculations based on the [START_REF] Newman | VolatileCalc: a silicate melt-H2O-CO2 solution model written in Visual Basic for Excel[END_REF] model). This is in line with the amounts estimated by [START_REF] Scaillet | A model of sulphur solubility for hydrous mafic melts: application to the determination of magmatic fluid compositions of Italian volcanoes[END_REF] (up to 5% of free fluid phase) and [START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF] (H 2 O TOT =3.4 wt%, CO 2 TOT = 2wt%; S TOT = 0.32 wt%), based on either experimental data or modeling, the latter by using the Moretti et al. (2003) model updated for the H 2 O-CO 2 saturation model of [START_REF] Papale | The compositional dependence of the saturation surface of H 2 O+CO 2 fluids in silicate melts[END_REF]. These high total contents strongly suggest the addition of deep fluids from a devolatilizing slab, infiltrating into the astenospheric mantle that upraises through the slab window over which Mt. Etna has developed [START_REF] Tonarini | Subduction-like fluids in the genesis of Mt. Etna magmas: evidence from boron isotopes and fluid mobile elements[END_REF]).

However, recent simulations of the plagioclase stability fields indicate a large variability in water content within the magmatic system, such that primitive Etnean melts are water-poor

A C C E P T E D M A N U S C R I P T

(between 1 and 2 wt%) and water enrichment occurs on magma ascent to the surface through fluxing of water induced by decompression and water exsolution at relatively shallow depth, and only subordinately by fractional crystallization [START_REF] Ferlito | The volatile flushing triggers eruptions at open conduit volcanoes: Evidence from Mount Etna volcano (Italy)[END_REF]Giacomoni et al., 2014).

Fluxing is typically due to continuous infiltration of CO 2 , which promotes magma dehydration through the plumbing system, and the volatile-induced differentiation with consequent enrichment of elements such as alkalis at shallow levels (Allard et al., 2006b;Spilliaert et al. 2006a;[START_REF] Ferlito | Volatile induced differentiation in the plumbing system of Mt. Etna volcano (Italy): evidence from glass tephra of the 2001 eruption[END_REF][START_REF] Collins | Melt inclusions track pre-eruption storage and dehydration of magmas at Etna[END_REF][START_REF] Ferlito | The role of supercritical fluids in the potassium enrichment of magmas at Mount Etna volcano (Italy)[END_REF]. CO 2 enrichment and H 2 O depletion through CO 2 fluxing occur stepwise, but are strongly recognizable in the 140-200 MPa range (Allard et al., 2006b;Spilliaert et al. 2006a;[START_REF] Collins | Melt inclusions track pre-eruption storage and dehydration of magmas at Etna[END_REF], i.e., at nearly isobaric conditions.

Fluxing allows the residing magma to re-equilibrate under different total volatile contents of H 2 O and CO 2 (Barsanti et al. 2009). The involved CO 2 -rich gas phase can be purely magmatic, separated at great depth (e.g., [START_REF] Clocchiatti | Coexistence de «verres» alcalins et tholeiiliques saturés en CO 2 dans les olivines des hyaloclastites d'Aci Castello (Etna Sicile, Italie). Arguments en faveur d'un manteau anormal et d'un reservoir profond[END_REF], or/and external and derived from decarbonation of the thick limestone platform present in the basement of Mt Etna at a depth > 4 km b.s.l. [START_REF] Lentini | The geology of the Mt. Etna basement[END_REF]Spilliaert et al. 2006a). The carbon isotope composition of the volcanic gas (Allard, 1986;Chiodini et al. 2011;[START_REF] Paonita | Geochemical evidence for mixing between fluids exsolved at different depths in the magmatic system of Mt Etna (Italy)[END_REF][START_REF] Rizzo | Real-time measurements of the concentration and isotope composition of atmospheric and volcanic CO 2 at Mount Etna (Italy)[END_REF] indicates that pre-1990 samples had somewhat lower  13 C (~ -3‰ vs PDB) than post-1990 samples (between -2 and -1 ‰ vs PDB), which may reflect a change in the source of gas supply. Such a change could be consistent with either (a) interaction with carbonate in the basement, or (b) increased mantle metasomatism by slabderived carbon-rich fluids [START_REF] Tonarini | Subduction-like fluids in the genesis of Mt. Etna magmas: evidence from boron isotopes and fluid mobile elements[END_REF]). However, important  13 C variations are intrinsic in the short-term magmatic activity of Mt. Etna. Measurements of CO 2 /SO 2 and of  13 C in CO 2 were carried during the 2014 activity and show that both indicators are concordant, revealing that early degassing of CO 2 may lead to peaks in both values and produce the highest  13 C values ever measured at Etna (1.5 ‰) [START_REF] Rizzo | Real-time measurements of δ 13 C, CO 2 concentration, and CO 2 /SO 2 in volcanic plume gases at Mount Etna, Italy, over 5 consecutive days[END_REF]. These results, in conjunction with the study of noble gas isotopes and abundances also confirm the occurrence of mixing between fluids exsolved at different depths in the magmatic system [START_REF] Paonita | Geochemical evidence for mixing between fluids exsolved at different depths in the magmatic system of Mt Etna (Italy)[END_REF].

A C C E P T E D M A N U S C R I P T

Sulfur is a sensitive tracer to investigate volcanic evolution to eruptions [START_REF] Oppenheimer | Mantle to surface gas trigger of the alkalic intraplate Erebus volcano, Antarctica[END_REF] and its speciation between SO 2 and H 2 S in the gas phase reflects the physico-chemical conditions of magmas feeding eruptions (e.g., [START_REF] A C C E P T E D M A N U S C R I P T Moretti | On the oxidation state and volatile behaviour in multicomponent gas-melt equilibria[END_REF][START_REF] Burgisser | Redox evolution of a degassing magma rising to the surface[END_REF]. Based on the composition of magmatic gases discharged from Etna volcano in July 1970 [START_REF] Huntingdon | The collection and analysis of volcanic gases from Mount Etna[END_REF][START_REF] Gerlach | Evaluation and restoration of the 1970 volcanic gas analyses from Mount Etna, Sicily[END_REF]), in 1975-1976(Allard 1986), in April 1983[START_REF] Le Guern | Ecoulements gazeux réactifs à hautes températures, mesures et modélisation[END_REF] and presently [START_REF] Aiuppa | H 2 S fluxes from Mt. Etna, Stromboli and Vulcano (Italy) and implications for the global volcanic sulfur budget[END_REF][START_REF] Aiuppa | Rates of carbon dioxide plume degassing from Mount Etna volcano[END_REF], SO 2 is the prevalent sulfur gaseous species at surface conditions. The SO 2 flux has traditionally been measured by using UV spectroscopy, yielding to a long-term time-average value of ~5,000 t/d with peaks to 25,000 t/d during eruptions and a large excess degassing rate of S compared to solid magma in the long term [START_REF] Allard | Endogenous magma degassing and storage at Mount Etna[END_REF][START_REF] Burton | FTIR remote sensing of fractional magma degassing at Mount Etna, Sicily[END_REF] for a review). Chlorine emissions, which are sourced by shallow (<50 MPa, Spilliaert et al. 2006a, b) degassing, average at ~1,000 to 2,000 t/d, although SO 2 /HCl proportions in the plume have been seen to vary substantially over time, and from one crater to another [START_REF] Aiuppa | Degassing of halogens from basaltic volcanism: Insights from volcanic gas observations[END_REF].

Samples, Analytical Procedures and Results

An olivine-hosted MIs study was performed on scorias extruded during the 1997-98 BN Strombolian eruption, during the paroxysmal event of February 1999 from the SEC (La [START_REF] Delfa | Activity of Mount Etna preceding the February 1999 fissure eruption: inferred mechanism from seismological and geochemical data[END_REF], and from products of the 122 BC Plinian eruption (Coltelli et al. 1998). Olivine crystals were handpicked under binocular microscope, from the 0.5-1 mm fraction, after gentle grinding and sieving. Because of rapid cooling of scoriae upon eruption, MIs were preserved as glass. Only in few cases they are-bubble free, and most often they display one bubble. The bubble/MI volume ratio does not appear constant and is highly variable, suggesting an heterogeneous trapping of a saturated magma (e.g., Spilliaert et al., 2006b). However, no attempt was made to determine if these are shrinkage bubbles (and therefore empty) or if they contain a CO 2 -or H 2 O-rich fluid phase.

Visually observed sulfide globules (Figure 5) are often associated with intergrowth of olivine and oxides, mostly (but not only) in the low-forsterite olivines.. They are typical of samples such that presented data and their processing were not affected by bubble decrepitation yielding CO 2 -loss from the inclusions during ascent (MacIennan 2017)

Melt inclusion chemistry and volatile contents

A total of 50 olivine phenocrysts (forsterite content ranging from 69 to 82 mol %; see Supplementary Material) were selected for studying MIs fully enclosed in their host crystal and embayments still in contact with their surroundings. Olivines in the analyzed products show forsterite (Fo) contents ranging from 69 to 83 mol.%, (Fig. 3). All the Fo range is recorded in the BN and SEC products, whereas the majority of olivine crystals from samples of the 122 BC eruption have Fo contents in the 74 to 83 mol % range. All these olivines entrap MIs characterized by highly variable volatile contents (H 2 O from 0.19 a 3.5 wt.%,, CO 2 from 81 a 1833 ppm, Cl from 0.05 a 0.46 wt.%, S from 0.017 a 0.39 wt.%). Because the highest Fo values are generally associated with the highest volatile contents (except for Cl), the data are consistent with degassing accompanied by crystallization, although the free exchange of crystals among different magma batches due to the occurrence of mixing/mingling processes cannot be discarded.

The occurrence of magmas differently enriched in potassium emerges from diagrams in Figure 4, in agreement with previous studies [START_REF] Metrich | The 1989-1990 activity of Etna: magma mingling and ascent of H 2 O-Cl-S-rich basaltic magma -evidence from melt inclusions[END_REF]Corsaro et al. 2009Corsaro et al. , 2013). An offset in K 2 O vs Na 2 O, K 2 O vs P 2 O 5 and K 2 O vs Cl (Fig. 4a-c) occurs between 122 BC and recent products. As a whole sulfur concentrations vary from  250 to > 4000 ppm, whereas chlorine contents vary from 500 to 4600 ppm over a wide interval of chemical differentiation (Fig. 3). MIs from the BN 1997-98 and SEC February 1999 eruptions encompass the same range of sulfur, chlorine and Being melt inclusions all from scoriae and explosive products, they were rapidly quenched upon eruption, preserving the volatile content of their pre-eruptive carrier melt [START_REF] Hartley | Diffusive over-hydration of olivine-hosted melt inclusions[END_REF]. As a matter of fact, MIs show a sharp contact with surrounding olivines. This indicates a very rapid cooling, which inhibited further exchanges between the melt droplet and the host crystal.

We discount full re-equilibration between melt and olivine after entrapment. If post-entrapment crystallization (PEC) occurred, it is negligible (see Appendix). Therefore, in the following elaborations we will use directly analytical data.

Similar considerations apply also to sulfides (Figure 5). Etnean magmas are generally FeS undersaturated, but may reach saturation under conditions in which sulfide separation would be driven by magma differentiation [START_REF] Scaillet | A model of sulphur solubility for hydrous mafic melts: application to the determination of magmatic fluid compositions of Italian volcanoes[END_REF][START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF]. Such an evolution is confirmed by Spilliaert et al. (2006b), who report sulfide globules in the samples from lava fountains of the 2002 flank eruptions. We discard the possibility of sulfide saturation as a postentrapment phenomenon within the inclusions. PEC being negligible (see Appendix), its eventual effect on S-concentration would be minor and unable to explain the presence of the episodically observed sulfides. Considering that stoichiometric pyrrhotite (density 2700 kg/m 3 ) contains 365,000 ppm of sulfide, if one restores S-rich ( > 2000 ppm) melt inclusions for the amount of sulfide contained in a globule occupying only 0.5% of the inclusion volume, unrealistic pre-entrapment concentrations (> 5,000 ppm) would be attained. This cannot be obtained via PEC at the scale of the single MI. Therefore, sulfides unmixed/precipitated as primary phases and were episodically trapped by MIs in virtue of their nugget effect concentrating sulfur, resulting in apparent extremely high local concentrations of sulfur.

Data Analysis and modeling

A C C E P T E D M A N U S C R I P T

Estimates from H 2 O-CO 2 data: pressure-depth and total volatiles

The new H 2 O-CO 2 data provided in this study (Supplementary Material) are plotted on saturation isobars (Fig. 6a) drawn with the [START_REF] Papale | The compositional dependence of the saturation surface of H 2 O+CO 2 fluids in silicate melts[END_REF] showing that volatile-rich Etna magmas can be found at depths > 15 km, at least) (Fig. 6a).

Particularly for the low soluble CO 2 component, dissolved contents are not representative of total abundances, given by the presence of an exsolved gas phase, in which H 2 O and CO 2 , represent more than 95 %. Especially at depth, these two components approximate the totality of such an excess gas phase (e.g., Moretti et al., 2003). On the saturation isobar diagram (Fig. 6a), each pair of This evidence of a CO 2 -rich gas phase already exsolved at large depth can be detailed by using the method described in [START_REF] Papale | Determination of total H 2 O and CO 2 budgets in evolving magmas from melt inclusion data[END_REF], focusing on the H 2 O-CO 2 pair dataset acquired in this study.

H 2 O-CO 2 dissolved
Following this author, we generated the total volatile (TV) lines, which represent the locus of total volatile contents consistent with the amounts dissolved in the liquid and measured in MIs (Fig. 6b).

Each TV line is generated by imposing different relative amounts of melt and gas phase and by calculating the H 2 O-CO 2 gas phase composition required to reproduce the (dissolved) H 2 O-CO 2 pair. Uncertainties in the procedure are then essentially related to the adopted melt-gas saturation model.

Along a TV line, the gas/melt ratio varies progressively, from zero (total contents = dissolved contents) to infinite (gas-dominated conditions) (Fig 6b (e.g., Spilliaert et al., 2006a;[START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF] and identify local reservoirs where CO 2 -rich magma is stocked [START_REF] Paonita | Geochemical evidence for mixing between fluids exsolved at different depths in the magmatic system of Mt Etna (Italy)[END_REF]. We then identify two possible fluxing arrays (arrows in Fig. 6c) qualitatively embracing all possible paths: a) CO 2 -gain (increasing CO 2

TOT

) and dehydration (left red arrow) and b) constant CO 2 TOT and dehydration (cyan arrow). Note that the reverse evolution of the magma system along such paths is not possible (unless reversing the pressure gradient) but could explain why some authors also refer to H 2 O fluxing throughout the Etna plumbing system [START_REF] Ferlito | The volatile flushing triggers eruptions at open conduit volcanoes: Evidence from Mount Etna volcano (Italy)[END_REF]. We must here remark that analytical uncertainties on H 2 O and CO 2 would make each TV backcalculation a field rather than a line, so intersections would correspond to vast polygonal areas spreading over a wide range of H 2 O TOT and CO 2 TOT values. However, our goal is not the

A C C E P T E D M A N U S C R I P T
identification of all populations of total volatiles needed to fit the whole MI database (e.g., Barsanti et al., 2009). Rather, we want to highlight the trends of either CO 2 -and H 2 O-fluxing and identify the total volatile contents that define a plausible initial pole for the parental magma. wt%) that are in line with those already used in [START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF]. However, the spread within IP region may also be determined by the infiltration of volatiles (essentially CO 2 ) released from the unsampled and deeper primitive magma.

Two reasonable equilibrium approaches can be followed to reproduce analytical data. 1) Using selected total volatiles, degassing trends are all generated within the IP for different redox states ( -3  NNO  3; Figure 7a). In this case, measured datapoints are embraced by model-generated trends, whose curvature and varying slopes reflect variations of the redox state. 2) Processes other than degassing occur, such as gas (CO 2 )-fluxing and/or the precipitation of an immiscible sulfide-bearing phase (Figure 7b), in agreement with visual observations (e.g., Fig. 5). These additional processes alter the total amount of dissolved volatiles, shift the beginning of degassing trends outside the IP,

A C C E P T E D M A N U S C R I P T

to other starting conditions that are depleted in H 2 O or S, and allow generating the many possible degassing trends which match analytical data.

Although the possibility that the redox state of a magmatic system is highly variable, the former approach returns unrealistic results. First, it is not possible to generate with selected total volatiles a reduced degassing trend from the IP able to cross data at H 2 O  2.5 wt% and S  0.35 wt%. Second, and more important, oxidized trends crossing through high-S and low-H 2 O datapoints (e.g., NNO +2 or +3; Figure 7a) are characterized by sulfur dissolved as sulfate (SO 4 2-

), which is incompatible with both speciation data (Moretti et al., submitted) and the visual occurrence of sulfide globules (Figure 5).

We then choose the second approach, that is, a combination of degassing trends plus FeS separation and CO 2 -fluxing. Simulations in Figure 7b are presented for logfO 2 at NNO, in agreement with [START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF][START_REF] Aiuppa | Hydrogen in the gas plume of an open-vent volcano, Mount Etna, Italy[END_REF]. However, they are unaffected by logfO 2 variations of ± 0.5 log units, that is, within one log unit. In all cases, degassing trends are almost rectilinear, particularly Tr1, Tr2 and Tr3, suggesting that sulfur and H 2 O have similar solubilities under these fO 2 conditions (Fig. 7b).

The effect of CO 140 MPa, where a major discontinuity and magma ponding zone has already been inferred [START_REF] Métrich | 2001 flank eruption of the alkali-and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades[END_REF][START_REF] Ferlito | The volatile flushing triggers eruptions at open conduit volcanoes: Evidence from Mount Etna volcano (Italy)[END_REF]; see also Figure 6a).

Trends Tr4 and Tr4bis (Fig. 7) start from relatively H 2 O-rich and S-depleted contents (here taken at 3.2 wt % and 1300 ppm for Tr4, and 3.2 wt % and 2300 ppm for Tr4bis) which may be

A C C E P T E D M A N U S C R I P T

obtained only after an early FeS separation (Tr5) [START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF]. This may have taken place in the early 122 BC magma, causing depression of the dissolved S content from the original total value (~ 3500 ppm; Figs. 3 and4), as also discussed in a companion paper (Moretti et al., submitted).

Finally, the arrow labeled Tr6 within the initial pole (IP) domain (Fig. 7) represents the H 2 O-S pattern followed by Etnean primitive melts that started crystallizing at depth, initially concentrating H 2 O and S when degassing is dominated by CO 2 exsolution (see also [START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF].

Therefore, early exsolution of a CO 2 -rich gas, due to crystallization at depth, fluxes the upward magma and allows retaining contents of soluble H 2 O and S up to 3.5 wt% and 0.45 wt%, respectively, embodying highest values measured in MIs (Supplementary Material and Figure 3).

Note that the adoption of H 2 O TOT and CO 2 TOT of 4 wt% and 2 wt%, respectively (see the bounding degassing trend in Fig. 6a), and S TOT of 0.41 wt% (highest measured dissolved S content;

Supplementary Material and Fig. 3), could be alternative to Tr6 within the error limits for S contents and would shift the IP toward the upper right corner of Figure 7b, that is at a higher starting pressure. This would not impact appreciably on trends departing from the IP (Tr1, Tr2 and Tr5).

The different behaviors displayed by H 2 O-S co-variations corroborate a view in which the magmatic system is parceled in different portions with variable total volatiles, reflecting the effect of CO 2 fluxing (Spilliaert et al., 2006a) and in agreement with the model of mixing of fluids exsolved at various depths developed by [START_REF] Paonita | Geochemical evidence for mixing between fluids exsolved at different depths in the magmatic system of Mt Etna (Italy)[END_REF] on conservative components such as CO 2 , Ar and He.

The role of crystallization and gas fluxing on chlorine release

The behavior of chlorine reflects the high solubility of this component in Etnean melts, which explains the constancy of chlorine concentration throughout chemical differentiation (Fig. 4c).

Although most chlorine tends to degas at shallow pressure, hence depth, its solubility behaviour depends on the relative H 2 O/CO 2 ratio in the gas phase [START_REF] Alletti | Chlorine partitioning between a basaltic melt and H 2 O-CO 2 fluids at Mount Etna[END_REF]. We then used equation ) in the frame of the Moretti et al. (2003) and [START_REF] A C C E P T E D M A N U S C R I P T Moretti | On the oxidation state and volatile behaviour in multicomponent gas-melt equilibria[END_REF] approach. Moreover, we also included the role of crystallization: a maximum crystallization of 67 wt% can in fact be evaluated on the basis of variations recorded by the P 2 O 5 component (from 0.5 to 1.5 %; Fig. 4b), considering P 2 O 5 as an incompatible component (apatite not on the liquidus; see also Spilliaert et al., 2006b) that saturation properties require chlorine to be already exsolved at large depth in the CO 2dominated gas phase. However, Cl-exsolution is very little, close to zero, for 122BC, such that the Cl-dissolved content corresponds to the total one at 400 MPa (Figure 8). In the case of the recent products, deep Cl-exsolution amounts to 500 ppm, that is, the difference between Cl TOT and the dissolved amount at ~350 MPa (1800 ppm; Figure 8). In agreement with the experimental evidences of [START_REF] Alletti | Chlorine partitioning between a basaltic melt and H 2 O-CO 2 fluids at Mount Etna[END_REF] on Etna samples, this deep exsolution is likely due to the high amount of exsolved CO 2 at depth. Differences in H 2 O TOT do not appear to appreciably affect the Cl-

A C C E P T E D M A N U S C R I P T

H 2 O evolution path, such that trends for H 2 O TOT at 3.5 wt% and 4 wt% overlap within uncertainty (see paragraph 4). Both systems exsolve chlorine progressively while depressurizing (solid lines).

When crystallization is considered, strong Cl-enrichment can be modeled. An important event of crystallization occurs in fact at shallow depth (P  20 MPa, i.e., H 2 O < 0.5 wt%) for both systems (122 BC and recent products), determining a change in slope in H 2 O-Cl trends (Fig. 8), and leading to increasing Cl contents. It is worth stressing here that, for each system, closed-system and opensystem degassing paths overlap until crystallization takes place. In case of recent products, low-P data are best embraced by a simple closed system degassing trend in absence of crystallization and an open-system degassing trend in presence of crystallization. In case of 122 BC, a simple closedsystem degassing trend, in presence of crystallization, matches well with available data.

The major feature emerging from simulations is then the difference in Cl TOT between 122 BC and recent products, which needs alternate explanations. One could be the involvement of different magma sources for 122 BC and recent activity, possibly complicated by mixing phenomena (e.g., Corsaro et al. 2009Corsaro et al. , 2013)). Some high chlorine concentrations can in fact be found in products from the recent activity that seem to pursue the 122BC trend at relatively high K 2 O values (Fig. 4). However, crystallization may easily account for the Cl-enrichment of these MIs, most of which cannot be reported in Fig. 8 as they lack the determination of H 2 O. Moreover, involvement of magmas from two different sources and their mixing does not afford the basic question: how to switch from the relatively Cl-rich, K-poor (or Na-enriched) 122 BC source to the relatively Cl-poor, K-rich (or Na-enriched) source of recent products ?

The alternate explanation is suggested by the experimental work of [START_REF] Alletti | Chlorine partitioning between a basaltic melt and H 2 O-CO 2 fluids at Mount Etna[END_REF], reporting that the introduction of CO 2 in the system tends to diminish D Cl g/m . Therefore, CO 2 -fluxing can operate a significant switch in chlorine contents, from a 122 BC-like parental melt to a term corresponding to that involved in recent activity. Two fluxing lines (one for a pure CO 2 gas, the other for a gas with H 2 O=5 wt% and CO 2 = 95 wt% ) are modeled in Figure 8, which allow the deep transition from the 122 BC trend (H 2 O TOT = 4 wt%; CO 2 TOT = 2 wt%; Cl TOT = 3500 ppm) to

A C C E P T E D M A N U S C R I P T

that of recent products (H 2 O TOT 3.5 wt%; CO 2 TOT  2 wt%; Cl TOT = 2300 ppm). Such a transition is modeled by extracting 20 wt% of gas phase and adding a corresponding amount of fluxing gas at each computational step corresponding to 5 MPa of depressurization. This produces a progressive CO 2 -enrichment and dehydration, and also depresses the dissolved chlorine content at concentration levels that at P < 350 MPa (Fig. 6a) are close to those typical of the basic trend for recent products, Cl TOT = 2300 ppm.

Discussion

Total volatiles and the shift of volatiles release patterns

Measured volatile abundances suggest a polybaric evolution of the Etnean system, in which simple paths of degassing and crystallization cannot be easily distinguished due to gas infiltration and magma mixing. Complex degassing behaviors are seen, particularly in case of highly soluble and reactive components, such as sulfur and chlorine (Fig 4). Such a complexity could reflect mixing processes, whose occurrence can partly explain the chemistry of BN and SEC products from recent activity. The observed variations (Fig. 4) can in fact result from self-mixing between i) the same volatile and K-rich end-member, which became substantial from 1999 to 2000 (Corsaro et al.

2013

) and is best represented by aphyric primitive products erupted during [START_REF] Newman | VolatileCalc: a silicate melt-H2O-CO2 solution model written in Visual Basic for Excel[END_REF]-2003(Spilliaert et al. 2006a;;Corsaro et al. 2009), and ii) its degassed and more evolved products. Self-mixing would also explain the presence of the whole spectrum of MIs in the Fo 80-70 olivine range (Fig. 3). Also the offset between 122 BC and recent products in K 2 O vs Na 2 O, K 2 O vs P 2 O 5 and K 2 O vs Cl plots (Fig. 4a-c) can be given by mixing processes superimposing to the normal trend of fractional crystallization [START_REF] Metrich | The 1989-1990 activity of Etna: magma mingling and ascent of H 2 O-Cl-S-rich basaltic magma -evidence from melt inclusions[END_REF]Corsaro et al. 2009Corsaro et al. , 2013)).

However, self-mixing alone cannot justify the existence of clearly fluxed magma portions (Fig. 6a) and the occurrence of variable intersection points in Figure 6c, characterized by variable amounts of H 2 O TOT and, particularly, CO 2 TOT which suggest a complex evolution of the Etnean

A C C E P T E D M A N U S C R I P T

magmatic system, due to gas addition or loss (Spilliaert et al. 2006a;[START_REF] Collins | Melt inclusions track pre-eruption storage and dehydration of magmas at Etna[END_REF][START_REF] Paonita | Geochemical evidence for mixing between fluids exsolved at different depths in the magmatic system of Mt Etna (Italy)[END_REF]. The set of total volatile contents in the parental magma (P > 400 MPa) have been identified at about 4 wt% and 2 wt% of H 2 O TOT and CO 2 TOT on products from the 3890 Plinian BP eruption and recent products (Fig. 6a and 6b; [START_REF] Kamenetsky | Arrival of extremely volatile-rich high-Mg magmas changes explosivity of Mount Etna[END_REF]. Magma parcels with total volatile contents even exceeding those of the parental melt have been recognized from products of the recent activity (Fig. 6b), revealing that a large role is played by CO 2 addition throughout the plumbing system, which is characterized by intervals in which large amount of volatiles are stocked. This evidence, which does not exclude self-mixing processes of the magma batches, is Our simulations (Figs. 6,7,and 8) have shown that CO 2 fluxing greatly determines the evolution of the other measured dissolved volatiles, in combination to crystallization (e.g., Fig. 7).

surely
In case of sulfur, modeling has demonstrated that CO 2 -fluxing may shift total volatile abundances producing the initial trends that on the H 2 O-S diagrams run through measured datapoints (e.g. trend Tr2 in Fig. 7b). However, our simulations were conducted by investigating melt-gas equilibria, independently of the presence of a FeS-bearing phase, either solid or liquid. Thus, it is not immediate to reconcile the presence of sulfide globules with high-S concentrations in such melts fluxed by CO 2 , because the S-budget of the whole trapped system (MI + sulfide globule) largely exceeds the value of S TOT . The effect is particularly evident for CSE 1999and, particularly, BN 1997-98 MIs that fall on Tr2 and Tr3 trends, and show the largest extent of CO 2 -fluxing and consequent dehydration (Fig. 6b). In these dehydrated MIs, sulfur, contrary to H 2 O, is retained and concentrates instead within the MIs, reaching high values despite the presence of sulfide globules

A C C E P T E D M A N U S C R I P T

(e.g., sample d9 from the BN 1997-1998 dataset; Supplementary Material) that concentrate large amounts of S. However, the presence of sulfide globules determines a "nugget" effect, such that the amount of sulfur in the sulfide globule plus the sulfur dissolved in the melt may not be representative of the actual S-budget within the volume of the trapped MIs, but of total sulfur in a much larger equilibrium volume. In such a case, trapped MIs carrying a sulfide globule have an apparently extremely high S-budget (dissolved sulfur in melt, plus sulfur in the sulfide globule) which instead should be distributed over the much larger but unconstrained equilibrium magma volume which includes also MIs not displaying sulfide globules.

Consequently, high S-budgets of some MIs may either reflect 1) unexpectedly high sulfur abundance inherited from the deep parental source or, 2) combination of multiple crystallization, increasing the S-content of the residual melt, and nugget effect by precipitating sulfide globules in a magma whose patchy distribution of sulfur is being recorded by MIs entrapment. CO 2 -fluxing, yielding dehydration and consequent crystallization (e.g., [START_REF] Blundy | Petrologic reconstruction of magmatic system variables and processes[END_REF], clearly relates to the second case, as testified by the sulfide-bearing BN 97-98 MIs compositions, more evolved than the sulfide-lacking 2001 ones.

Also [START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF], in their modeling of sulfur isotopes ( 34 S) from Etnean volcanics, highlighted the two above possibilities. Starting from an initial mantle value ( 34 S = 0 ‰ vs CDT), they were in fact able to model measured  34 S values assuming an early trend of FeS separation, and showed that that the same initial condition for fractionation patterns of  34 S can be obtained either by starting from a parental melt with S  1.1 wt% from which only FeS precipitates, or by allowing the CO 2 -rich parental melt (with S TOT = 0.32 wt% ) to degas and crystallize. Based on their independent findings, early degassing at depth is likely to lead to S (and also H 2 O) enrichment, similarly to Tr6 in Figure 7, unlike FeS separation [START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF]. CO 2 exsolved in such an early degassing stage fluxes the magma ponding upward at lower pressures, and according to our simulations (e.g., Tr2, Fig. 7b), makes the melt a sort of repository of the highly soluble S component, thus preventing its loss.

A C C E P T E D M A N U S C R I P T

On the other hand, trend Tr5 (Fig. 7b), along which degassing trends such as Tr4 and Tr4-bis start taking place, can be explained by simple separation of FeS from IP, which in absence of CO 2fluxing, determines S-depletion in the melt. Sulfide separation would require the onset of more oxidized conditions along Tr5 proportional with the removal of S 2-and Fe 2+ , which make up the FeS separated phase, thus yielding high residual proportions of Fe 3+ and SO 4 2-in the residual melt.

This is confirmed by the large amount of S as sulfate from XANES analyses of hydrous ( H 2 O  2.5 wt%) MIs from 2001 eruption. (Moretti et al., submitted).

Summarizing, earlier separation of sulfide globules can affect the IP and its sulfur total content.

Sulfide separation may occur in absence of CO 2 fluxing (Tr5), substantially resulting in the depletion of the melt S-content (Tr5). On the contrary, in presence of CO 2 fluxing (Tr2), the melt retains high concentration of sulfur. Subsequent degassing trends will take place from such arrays explaining the spread of data observed on the S-H 2 O diagram (Fig. 7). Direct degassing from IP (Tr1), without sulfide separation, may also occur and appears as an upper limit of the former process. Finally, the IP may experience joint crystallization and degassing (Tr6) and consequent CO 2 -fluxing of upward magma parcels, leading to enrichment of highly soluble sulfur and water components.

Consistent with such an early degassing, also probed on the basis of sulfur isotopes [START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF], the IP itself can represent an apparent initial condition, being affected by volatile infiltration (CO 2 -addition) from a deeper, unsampled, exsolving primitive magma likely characterized by a lower CO 2 TOT /H 2 O TOT than the IP. This possibility is enlightened by our simulations of how CO 2 -fluxing can determine a negative shift in chlorine dissolution from the 122 BC to the K-enriched MIs from recent eruptive products (Fig. 8). Such a shift of Cl TOT and Cl degassing patterns is associated to a change of the Na/K ratio that may reflect vapor transport of alkalis, particularly K (e.g. [START_REF] Ferlito | The role of supercritical fluids in the potassium enrichment of magmas at Mount Etna volcano (Italy)[END_REF]Moretti et al., 2013a). Extraction of metal cations seems in fact to increase with increasing Cl content of the gas/fluid phase (Webster et al., 2011 and references therein). Therefore, the occurrence of fluxing may partly explain the so-called

A C C E P T E D M A N U S C R I P T

low-K and high-K differentiation trends of Etnean magmas, as their concomitant occurrence in same products likely reflects self-mixing within degassed and undegassed, or even fluxed, parcels of the same plumbing system. This phenomenon is likely to sum up to the contribution from slabderived fluids, which appears to be dominat in light of the inferred decoupling between alkalis and Cl (Corsaro and Métrich, 2016).

However, we cannot push further these arguments, also due to the intrinsic limits of the reported C-H-O-S-Cl-melt modeling. In fact, model results on the effects of CO 2 -fluxing (Fig. 8) are based on the gas-melt chlorine partition coefficient given in [START_REF] Alletti | Chlorine partitioning between a basaltic melt and H 2 O-CO 2 fluids at Mount Etna[END_REF], for which on one side the authors suggest an expression only dependent on pressure, but on the other side they report that the introduction of CO 2 in the system tends to diminish the partition coefficient. At present the latter evidences cannot lead to a H 2 O/CO 2 -dependent partition coefficient of Cl, in light of the uncertainty of data, which are derived by mass balance and do not allow a full characterization of the speciation state of the gaseous phase [START_REF] Alletti | Chlorine partitioning between a basaltic melt and H 2 O-CO 2 fluids at Mount Etna[END_REF]. A definite assessment of the relative roles of H 2 O and CO 2 on Cl dissolution requires new experimental data as well as a solubility model for chlorine accounting for chemical interactions with the other melt and gas species.

Relations with eruptive styles

By considering the ensemble of data and simulations reported in Figures 678, the following observations and deductions emerge:

1) An early degassing stage characterizes the deepest and least differentiated magma that is potentially common to all studied products. Such a magma is identified from MIs involved in February 1999 and 2001 activities (Fig. 6c), for which the early degassing stage coupled to deep crystallization is recognized around the IP (Tr6 in Figure 7), also according to the S-isotope modeling of [START_REF] Marini | Sulfur Isotopes in Magmatic-Hydrothermal Systems, Melts, and Magmas[END_REF]. Such an early degassing stage may produce H 2 O and S-enrichment and, in particular, produces the huge flux of

A C C E P T E D M A N U S C R I P T

CO 2 infiltrating upward and characterizing the Etna activity [START_REF] Allard | Eruptive and diffuse emissions of CO 2 from Mount Etna[END_REF]. The early degassing stage is inferred to be a stable and major feature of the Etnean plumbing system, associated with the growth of the deep crystal mush [START_REF] Allard | Endogenous magma degassing and storage at Mount Etna[END_REF]). An early degassing stage, and consequent CO 2 -fluxing, is then inferred to produce in time the switch from the 122 BC primitive magma to that involved in the recent activity (Fig. 8).

Alkali transport in the fluxing gas phase may explain differences in the Na/K ratio (Fig 4), which is lower for the 122 BC and higher for products from recent activities.

2) All products record evidences of CO 2 -fluxing, especially at P  100 MPa, as shown by the TV lines (Fig. 6c,d). The upper limit of CO2-fluxing is given by the gas-buffered trend at CO 2(gas) /H 2 O (gas)  5, which takes place for P < 200 MPa (Fig. 6a) and is associated with low-density and buoyant magma batches within the plumbing system.

3) Recent eruptive products show that CO 2 -fluxing also characterizes deep reservoirs (P up to 200 MPa and even more), as revealed by the shift of the IP in the H 2 O-S covariation of figure 7, and the position of TV lines of 1999 and 2001 products, which tend to identify more than one stocking zone of CO 2 -rich volatiles (Fig. 6c). The same cannot be inferred for the 122 BC system: TV lines show the occurrence of open system conditions (Fig. 6d), suggesting the absence of important stocking volumes of volatiles within the plumbing system and prior to erupt. Furthermore, 122 BC TV lines point to H 2 O-fluxing superimposing to an earlier event of CO 2 -fluxing.

4)

Simulations of chlorine dissolution (Fig. 8) reveal that 122 BC plinian products experience a final stage in closed-system degassing conditions, thus retaining the whole gas phase. Therefore, despite CO 2 -fluxing took place throughout the 122 BC plumbing system, gas stocking and accumulation characterizes only the very late and low-P preeruptive stage. On the other hand, recent products are better explained by a final evolution under open-system degassing conditions, which also characterizes latest preeruptive conditions.

A C C E P T E D M A N U S C R I P T

Summarizing, CO 2 -fluxing is a common phenomenon at Mt. Etna, likely related to important CO 2 release during the early stages of degassing and crystallization in the roots of the plumbing system. Nevertheless, in case of 122 BC, the continuous structure of the infiltrating flow allowed open-system degassing conditions to take place throughout the entire magmatic system, with the gas accumulating uniquely at the top of the system prior to eruption (figs. 8 and9). This favored the subsequent H 2 O-fluxing, associated to important water exsolution in the shallow system, very likely in response to crystallization. On the contrary, the plumbing system of recent activities is characterized by multiple storage levels, where the gas exsolved at large depth accumulates in magma batches evolving under closed-system conditions (Figure 9). These magma batches are well represented by products emitted in low to intermediate explosive eruptions, such as the BN 1997-98

and February 1999 episodes (Fig. 6b,7): they rise under volatile rich conditions and experience open-system degassing conditions prior to eruption (Figs. 8 and9). Even among the products from recent activities distinctions are possible. In fact, MIs from products relative to the February 1999 event, characterized by highly explosive phases, suggest that the erupted magma experienced lower CO 2 -fluxing and dehydration compared to the BN1997-98 eruption: the Tr1 trend, representing the closed system degassing trend generated from the initial IP condition, crosses many February 1999 data, which are not particularly affected by CO 2 -infiltration and which also do not represent remobilized crystals that experienced different ascent and degassing pathways (Figure 7).

We then suggest that CO 2 -fluxing, hence dehydration, may be inversely related to explosivity, because produces relatively small magma batches characterized by high gas/melt ratios (up to ~ 0.3 in mass; Fig. 6a) which are mobilized upward in gas-dominated conditions. Such magma batches, being already dehydrated by CO 2 -fluxing, during late decompression and pre-eruptive stages cannot exsolve high water amounts that significantly increase the explosive power: CO 2 -fluxing keeps the magma supersaturated under high gas/melt ratios and inhibits the role of decompression, thus favoring the buoyancy of water-depleted small magma batches, erupted in effusive or mildlyexplosive conditions. Such low-density and buoyant magma batches fall on a gas-buffered trend at

A C C E P T E D M A N U S C R I P T

CO 2(gas) /H 2 O (gas)  5, a value which appears to be typical for all volcanic systems in the western Mediterranean region [START_REF] Mangiacapra | The deep magmatic system of the Campi Flegrei caldera (Italy)[END_REF][START_REF] Arienzo | The feeding system of Agnano-Monte Spina eruption (Campi Flegrei, Italy): Dragging the past into present activity and future scenarios[END_REF][START_REF] Pino | Seismic precursors of a basaltic paroxysmal explosion track deep gas accumulation and slug upraise[END_REF]Moretti et al., 2013) and which likely marks the density and rheological thresholds for buoyancy.

On the other side, data on the 122 BC pre-eruptive volatile content, along with micro-textural evidence [START_REF] Sable | Changing conditions of magma ascent and fragmentation during the Etna 122BC basaltic Plinian eruption: evidence from clast microtextures[END_REF], indicate that basaltic Plinian volcanism at Mt Etna in 122 BC was enabled by degassing associated with phenomena of H 2 O-fluxing and high rates of bubble nucleation, decompression and magma ascent, enhanced by abundant microlite crystallization and consequent increase of magma effective viscosity [START_REF] Sable | Changing conditions of magma ascent and fragmentation during the Etna 122BC basaltic Plinian eruption: evidence from clast microtextures[END_REF]. Such a decompression and consequent magma ascent could have been triggered by the initial emission of the volatile-rich magma evolved in closed-system degassing conditions on top of the system.

In general, we tentatively propose that explosivity at Mt. We suggest, on these grounds, that it is possible to build up a criterion to understand the evolution toward either low-to mildly explosive or highly explosive -up to Plinian-eruptions at basaltic volcanoes. An important step in the application of this criterion would be the joint monitoring of ratios involving CO 2 , sulfur and chlorine, in discharged gases and MIs from emitted scorias and lavas. In fact, the different roles of CO 2 -fluxing and CO 2 -separation could be evaluated by their contrasting effect on the solubility of other volatiles such as sulfur and halogens. For example, CO 2 addition to the system and consequent melt dehydration are expected to keep

A C C E P T E D M A N U S C R I P T

dissolved S as high as the IP region, until shallow degassing takes place (Figure 7). At same time, CO 2 addition and dehydration soon depress Cl dissolution in the melt (Figure 8). On the other hand, products of highly-explosive to Plinian activity should be characterized by decreasing amounts of dissolved S, probably due to early sulfide separation, but also by high contents of dissolved Cl.

However, further studies on sulfur and particularly chlorine dissolution in presence of variable amounts of CO 2 are needed to confirm this internal mechanism.

Are Etnean MIs really representative of the melt-volatile equilibrium ?

By simulating melt-plagioclase equilibrium constrained by plagioclase occurrence in the alkaline lavas of Mt. Etna, [START_REF] Ferlito | The volatile flushing triggers eruptions at open conduit volcanoes: Evidence from Mount Etna volcano (Italy)[END_REF][START_REF] Mollo | Reconstruction of magmatic variables governing recent Etnean eruptions: Constraints from mineral chemistry and P-T-fO 2 -H 2 O modeling[END_REF] suggest that the "high" water content found in MIs is not representative of the deep magma source, but results from H 2 Ofluxing, i.e, from a process of upward water migration that would enrich the originally waterundersaturated ( 1.4wt%) magma ponding at various levels within the feeding system (Ferlito et al., 2014). Water fluxing within the mushy plumbing system is consistent with the late degassing of the residual gas phase formed by a continuously crystallizing magma which has already exsolved CO 2 in its early fractionation stages. However, from MIs of recent Etna products we do not have direct evidences of such a process: the possible fluxing paths that we identified involve increasing CO 2 TOT /H 2 O TOT with decreasing pressure (i.e., CO 2 -fluxing), rather than a decreasing ratio as required by H 2 O-fluxing (Fig. 6c). We also remark that Giacomoni et al. ( 2015) report results based on MELTS simulations in which CO 2 cannot be considered, such that the modeled system is necessarily underaturated at depths of interest. The role of a CO 2 -rich gas phase in such simulations, together with a better appraisal of the role of fO 2 , should be better investigated by experimental petrology in order to appropriately recognize the extent to which MI-based and melt-plagiocase geothermobarometries agree. For example, experimental phase equilibria on the volatile-rich Kbasaltic magmas of Stromboli [START_REF] A C C E P T E D M A N U S C R I P T Di Carlo | Experimental crystallization of a high-K arc basalt: the golden pumice, Stromboli volcano (Italy)[END_REF] as Etna most of CO 2 must be stored into primary bubbles associated with the magma trapped as MIs. Further loss to bubble is however possible during post-entrapment cooling and contraction.

This would contribute to the observed dispersion and of course to underestimating pressure [START_REF] Moore | Bubbles matter: An assessment of the contribution of vapor bubbles to melt inclusion volatile budgets[END_REF][START_REF] Wallace | Melt inclusion CO2 contents, pressures of olivine crystallization, and the problem of shrinkage bubbles[END_REF]MacIennan, 2017). However, following [START_REF] Wallace | Melt inclusion CO2 contents, pressures of olivine crystallization, and the problem of shrinkage bubbles[END_REF], the generally high amounts of water in our inclusions are expected to decrease the CO 2 -loss to vapor bubble due to shrinkage. However, even if our recomputed pressures were underestimated, this would not diminish appreciably the representativeness of investigated MIs and our reconstruction of the Etna system. In fact the CO 2 -fluxing phenomenon here recognized would be even larger, then reinforcing our findings as an explanation to the well known gas/melt imbalance of Etna, which results in too much gas compared to produced magma (e.g., [START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF]. Noteworthy, our pressure estimates agree with the independent petrologic geobarometer of [START_REF] Armienti | A new model to estimate deep-level magma ascent rates, with applications to Mt. Etna (Sicily, Italy)[END_REF].

Conclusions

The trigger processes of basaltic Plinian eruptions, observed at Etna and likely many other open-vent volcanoes, are still far from being completely understood (e.g., Houghton and

A C C E P T E D M A N U S C R I P T

Gonnermann 2008). However, explosive basaltic volcanoes such as Etna are characterized by very abundant volatiles in the system, as testified by the contrast between the huge gas but low magma emission. How such gases are released can thus be related to the ample spectrum of activity at this volcano. In this study we have then showed that a careful analysis and modeling of degassing patterns within the H 2 O-CO 2 -S-Cl-melt system can give insights into the shift from low-explosive to Plinian activity. MIs in fact record the variable distribution of the amount of gas at equilibrium with melts throughout the Etna plumbing system and allow discriminating the various processes which affect its activity. These are:

1) self-mixing processes superimposing to the normal trend of fractional crystallization;

2) early degassing and crystallization of a deep (P > 200 MPa) magma, contributing to the endogenous growth of the mushy system, and releasing a huge amount of CO 2 fluxing the upward system 3) CO 2 -infiltration throughout the plumbing system, which can either result in the preeruptive accumulation of gas at various levels within the plumbing system (recent activities), or in a direct stocking at the top of system (122 BC).

4) Subsequent H 2 O-fluxing, which can be inferred for 122 BC, due to large water exsolution in the shallow magmatic system.

Gas accumulation at various levels yield high gas-melt (up to ~ 0.3 in mass, with molar CO 2,gas /H 2Ogas  5) magma batches, which experience strong dehydration and crystallization, and which are mobilized upward in low-explosive eruption characterized by little magma emission (e.g., BN 1997-98 activity). When the huge amount of fluxing gas accumulates essentially at the top of the system, MIs show open-system degassing conditions throughout the plumbing system associated with H 2 O-fluxing, but highlight the occurrence of a last pre-eruptive step in closedsystem degassing conditions, likely favoured by abundant microlite crystallization and consequent increase of magma effective viscosity. This determines the Plinian activity, like in the 122 BC case.

Continuous geochemical and petrologic monitoring within the H 2 O-CO 2 -S-Cl-melt system can thus

A C C E P T E D M A N U S C R I P T

help recognizing whether the magmatic system is tending to this latter scenario or not. However, further experimental studies about S and Cl dissolution at different H 2 O (gas) /CO 2(gas) and gas/melt are required to attain a highly reliable parameterization.

  reactivation in 1996 up to the onset of 2001 flank eruption, with the alternation of fractional crystallization and mixing between residing and new intruding magma batches, controlled by geometry of the conduit and the balance between refilling magma and erupted lavas. On the contrary, an upper storage system, corresponding to the central conduits region (i.e. depths <100 MPa, Collins et al. 2009), with different size and/or supply rate is hypothesized for the other three summit craters (Corsaro et al. 2013). The post-1970 activity of Mt. Etna has been featured by short but intense periods of explosive activity (e.g. the 2001 July-August eruption), and even of sub-Plinian scale, such as the paroxysmal event on 22 July 1998 at VOR crater. Such a change in the eruptive intensity has to be related to the change in magma composition, initiated during the 1974 flank eruption (Corsaro et al. 2009), with intrusion of a new volatile-rich and more alkali-rich basaltic magma component. This magma component, which has come to play an increasing A C C E P T E D M A N U S C R I P T

  Figure 2 reports the Total-Alkali-Silica classification diagram (TAS; Le Maitre 1979). Only for this diagram, compositions have been normalized to 100 wt% on a volatile free basis, in order to compare them with whole rocks. MIs from SEC and BN (1999 and 1997-98 respectively), are K-trachybasalts and K-basaltic trachyandesites displaying a strong tephritic affinity. Those from products of the 122 BC eruption fall in the field of intermediate terms of alkaline series, which are ascribed to trachybasalts and A C C E P T E D M A N U S C R I P T trachybasaltic andesites. Because of their sodic affinity they may be classified as hawaiites or as hawaiites and mugearites.

H 2 O

 2 , CO 2 and S contents versus K 2 O have been plotted in Figure4d-f showing the peculiar enrichment in H 2 O and CO 2 , but not S, of low-K products. H 2 O and CO 2 are generally negatively correlated with K 2 O, indicating degassing coupled to crystallization. However, this occurs along multiple paths and is complicated by variable degrees of potassium enrichment, particularly evident in the K 2 O vs H 2 O diagram (Fig.4d). The negative correlation with K 2 O is less evident for sulfur, especially because low-K 2 O MIs from the 122 BC Plinian event show lower sulfur contents (S < 2000 ppm) than those from the other analyzed products (Fig.4f).

  H 2 O contents (S  ~2500 ppm, 1000 < Cl < 4600 ppm and H 2 O  ~ 1.2 wt%), but few MIs from BN 1997-98, although hosted in a Fo 80 olivine, are characterized by the highest observed values of A C C E P T E D M A N U S C R I P T dissolved sulfur (S~ 4200 ppm wt). The low K 2 O (2 wt%), P 2 O 5 (0.59 wt%) and high CaO/Al 2 O 3 ratio (0.7) of these samples reflect the primitive nature of these MIs.

  Despite the availability of only two H 2 O-CO 2 pairs, conditions of important CO 2 -enrichment are recorded by BN 1997-98 activity, yielding very high CO 2 TOT /H 2 O TOT (mass) ratios and crossing at CO 2 TOT /H 2 O TOT ≈ 20 (dashed lines; Fig. 6c), a value consistent with the trend bordering datapoints on the left side of Fig.6a. TV lines relative to MIs in crystals from the 122 BC Plinian samples do not allow inferred total volatile contents consistent with those of February 1999 and 2001 products. In fact, we miss intersection points, and TV lines depict a divergent fan, especially if we include the data from Del A C C E P T E D M A N U S C R I P T Carlo and Pompilio (2004) (dashed lines in Fig 6d), which reveals the onset of complex patterns of either gas loss or gain in open system conditions. With respect to MIs from recent activity, TV lines generated by 122 BC MIs are shifted toward lower CO 2 TOT /H 2 O TOT ratios, with four of them collapsing along the abscissa, which suggests important CO 2 -loss. On the other hand, a few of them (five) are steep and reveal dehydrated and CO 2 -rich conditions: they approach the highest CO 2 TOT /H 2 O TOT ratios observed for BN 1997-1998 MIs, reflecting extensive CO 2 fluxing and addition producing the gas-buffered trend in Fig. 6a. The never crossing TV lines generated from 122BC MIs could be compatible with H 2 O-fluxing (Ferlito et al., 2014), superimposing to previous CO 2 -fluxing. Water fluxing, may have occurred in the shallow (P ≤ 150 MPa; Figure 6a) part of the 122 BC plumbing system, in response to exsolution of the highly soluble water component. This could have produced the flattening of the H 2 O-rich side of Figure 6c. Therefore, 122 BC samples display peculiar features of the extensive open-system degassing occurred at large depth: the magmatic system, first, got rid of its CO 2 , fluxing the shallower system but without accumulating in intermediate magma storage zones like those observed throughout the plumbing system of recent activity; secondly, subsequent water exsolution and fluxing in the shallow system (P ≤ 150 MPa) occurred, as recorded by TV lines.This may suggest that the set of total volatiles, particularly CO 2 TOT , inferred for recent products may not be fully representative of the initial volatile content of the primitive magma: instead, it may represent an intermediate accumulation level resulting from fluxing of CO 2 released from an even deeper but unsampled pristine source at the roots of the Etnean system. However, for the purposes of this study, and in agreement with[START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF], we will retain H 2 O TOT  3.5 wt% and CO 2 TOT  2wt%.

4. 2

 2 Figures 7a,b report the S-H 2 O covariation and degassing trends from a region (labeled IP) of

  2 addition, due to fluxing, on the shift of H 2 O-S degassing patterns is accounted for by trend Tr2, which shifts Tr1 to Tr3 and well explains the largest extent of dehydration displayed by BN 97-98 and CSE 1999 MIs (Fig. 6b,7b). Along Tr2, sulfur, contrary to H 2 O, is retained and concentrates instead within the MIs. This array interpolates the starting poles of the many possible degassing paths (such as Tr3) that commence from progressively dehydrated melts because of fluxing or a combination of fluxing plus self-mixing between deep volatile-rich and degassed magma portions. The fluxing array (Tr2) ends and the degassing trend Tr3 starts-at

  et al. (2009) [ D Cl f/m = 12.732 x P -0:0853 ; P in MPa ] to introduce a P-dependent chlorine partition coefficient between gas and melt (D Cl g/m

  reminiscent of the mixing model of fluids released at different depths proposed by Paonita et al. (2012). On the other hand, products from the 122 BC Plinian eruption depict a somehow different scenario, in which open-system degassing and continuous gas extraction takes place since the deep parental magma region. The consequent lack of crossing TV lines does not allow to definitely constrain total volatiles for the 122 BC (Fig. 6c Papale, 2005), which are considered to be the same as for recent activity and the analogous 3890 BP eruption (Fig. 6a,b).

  Etna may be governed by a balance involving either CO 2 -fluxing (recent products) or H 2 O-fluxing (122 BC) and crystallization throughout the plumbing system. If magma batches remain hydrous, also because of H 2 O-fluxing, the huge and rapid exsolution of water makes them prone to erupt explosively: considering the molar weights of water and anhydrous melt, rapid exsolution of ~ 3 wt% in the last 100 MPa of decompression (Figure 6a) corresponds to the transfer of ~ 10 mol% of H 2 O from the liquid to the gas phase. Fast decompression can effectively determine high H 2 O super-saturation in the melt, and trigger very rapid nucleation of H 2 O-dominated bubbles, leading to a strong acceleration of the magma.

Figure 1 .

 1 Figure Captions

Figure 2 .

 2 Figure 2. TAS diagram of samples object of this study. Only in this case, data (see Supplementary

Figure 3 .

 3 Figure 3. Concentration of dissolved H 2 O (panel a), CO 2 (panel b), S (panel c) and Cl (panel d)in

Figure 4 .

 4 Figure 4. Panel a) Plots of Na 2 O vs K 2 O, allows discriminating among magmas characterized by

Figure 5 .

 5 Figure 5. Etna MI hosted in olivines coexisting with sulfide globules and multiple oxide minerals.

Figure 2 A

 2 Figure 2

Figure 3 A

 3 Figure 3

Figure 4 A

 4 Figure 4

Figure 5 A

 5 Figure 5

Figure 6 A

 6 Figure 6

Figure 8 A

 8 Figure 8

Figure 9 A

 9 Figure 9

  measured CO 2 concentrations of ~ 98-99 wt% in monophasic fluid inclusions hosted in high-Mg olivines (Fo 92 mol%) from the Mt. Maletto suite (7

ka BP;

Condomines et al. 1995)

, implying that CO 2 is already exsolved at a pressure of 800 MPa, i.e., at a lithostatic depth of about 21-24 km, corresponding to the base of the crust. Large production of CO 2 is testified by the background CO 2 volcanic flux of 4 ktons/day, which was assessed in

see Appendix and Supplementary Material), Infrared

	characterized by low sulfate proportions (Moretti et al., submitted) and were also found in
	clinopyroxene and plagioclase crystals..
	We just provide here a summary of the conditions used for minerals and MIs analyses,
	whereas the details of the analytical procedures are provided in Appendix. Major elements,
	chlorine, phosphorous and sulfur concentrations in MIs were measured by electron microprobe
	(EMP) CAMECA-SX 50 (CAMPARIS Jussieu, France). Major elements in melt inclusions were analyzed with 10 nA beam current, 10 mm spot size, and counting time of 10 to 25 s. Sulfur, chlorine and phosphorus analyses were carried out with 30 nA beam current, 15 mm spot size, counting times of 100 to 200 s. Analytical uncertainties were 3.5% for Si, Al, Mg, Fe, Ca, K and Na, and 5% for Cl, S and P. On a subset of doubly polished MIs (Spectroscopy measurements of dissolved H 2 O and CO 2 were carried out prior to any other analysis, by using a NICOLET ® Magna 550, associated to a microscope Spectra-Tech ® IR-plan. Errors are M A N U S C R I P T 10% for H 2 O and at most 25% for CO 2 . We have of course discarded inclusions lying on cracks,
	A C C E P T E D

  model, for a temperature of 1200°C, which is close to the 1220°C value of[START_REF] Kamenetsky | Arrival of extremely volatile-rich high-Mg magmas changes explosivity of Mount Etna[END_REF] for deep volatile-rich Etnean magmas. Note that this temperature is only 10% greater than temperatures recorded by lavas extruded at surface vents (close to 1090°C) during the February 1999 eruption (La Delfa et al.

2001

). Maximum entrapment pressure was computed at ~300 MPa for MIs from crystals hand picked from the February 1999 sample, and agrees well with previous results available in the literature (Fig.

6a

).

Although we did not analyze H 2 O and CO 2 contents by FTIR on the majority of MIs entrapped in primitive olivine crystals from products of the 122 BC eruption, literature data on products from this and from similar Plinian eruptions, such as the 3930 BP event

(Coltelli et al. 2000; Del Carlo and Pompilio 2004;[START_REF] Kamenetsky | Arrival of extremely volatile-rich high-Mg magmas changes explosivity of Mount Etna[END_REF]

, are characterized by H 2 O and CO 2 contents up to 4 wt% and 4000 ppm, respectively, yielding the highest MI entrapment pressures ( 500 MPa) and

  in MIs can be produced by a virtually infinite number of degassing trends,Spilliaert et al., 2006a;[START_REF] Arienzo | The feeding system of Agnano-Monte Spina eruption (Campi Flegrei, Italy): Dragging the past into present activity and future scenarios[END_REF] Moretti et al., 2013b). The curvature of such a trend around the 200-100 MPa region, where the contrast of H 2 O and CO 2 solubilities is such that H 2 O starts exsolving significantly and little CO 2 stays dissolved above typical detection limits for

	2004; such a component, is greatly controlled by H 2 O TOT and especially CO 2	TOT .
	In Figure 6a, we see that the bulk of data can be bounded by a closed-system degassing trend
	characterized by high total volatile contents (H 2 O TOT = 4 wt%; CO 2 is up to ten times more abundant than the largest dissolved amount measured in a SEC 1999 MI, TOT TOT = 2wt%), such that CO 2 M A N U S C R I P T and ~5 times larger than highest dissolved contents from literature data. These values are in line
	A C C E P T E D
	depending 1) on the presence and relative amount of the initial gas phase, hence total
	(exsolved+dissolved) H 2 O and CO 2 volatile abundances (H 2 O TOT and CO 2	TOT ), and 2) the degassing

style (closed-system, open-system, or a combination of both). As a rule of thumb, it is considered that the closed-system degassing trend crossing the rightmost datapoints on a H 2 O-CO 2 saturation isobar diagram well approximates the total volatile abundance of the system (e.g.,

Métrich et al.,

with earlier determination of

[START_REF] Scaillet | A model of sulphur solubility for hydrous mafic melts: application to the determination of magmatic fluid compositions of Italian volcanoes[END_REF]

and consistent with those of

[START_REF] Aiuppa | Forecasting Etna eruptions by real-time observation of volcanic gas composition[END_REF]

. On their left side, datapoints are limited by a gas buffered trend (gas/melt mass ratio ~ 0.3 ) for a molar CO 2(g) /H 2 O (g) ratio of 5. This trend, which involves low-density magmas basically unaffected by degassing because buffered by the exsolved volatiles, is in good agreement with same trends observed from Stromboli K-basaltic MIs (molar CO 2(g) /H 2 O (g) >3;

[START_REF] Pino | Seismic precursors of a basaltic paroxysmal explosion track deep gas accumulation and slug upraise[END_REF]

, latitic and trachytic MIs from Campi Flegrei (molar CO 2(g) /H 2 O (g)  4;

[START_REF] Mangiacapra | The deep magmatic system of the Campi Flegrei caldera (Italy)[END_REF][START_REF] Arienzo | The feeding system of Agnano-Monte Spina eruption (Campi Flegrei, Italy): Dragging the past into present activity and future scenarios[END_REF]

, and particularly K-trachybasaltic MIs from Ischia (molar CO 2(g) /H 2 O (g)  5;

Moretti et al., 2013b)

.

  , all TV lines generated from dissolved H 2 O-CO 2 pairs must converge to the same point, representing the total volatile content of the system (Fig.6binset). In case of gas separation during the degassing process (open-system conditions) TV lines will not converge to one common point.

		Therefore, in Figure 6c, points where TV lines cross identify H 2 O TOT -CO 2	TOT pairs that
	represent the total volatile contents for the considered samples. In particular, primitive MIs from the
	February 1999 and 2001 activity generate TV lines (solid lines) crossing at several H 2 O TOT -CO 2	TOT
	pairs, revealing that a quite complex degassing history and style also affects the deep magmatic
	system. Line intersections give values from 3.4 to 3.8 wt % (for H 2 O TOT ) and 1.9 wt% (for
	CO 2	TOT ), which best represent the initial concentrations in February 1999 parental melts. In the
	shallower portions of the plumbing system, the magma is characterized by larger CO 2	TOT /H 2 O TOT

). In case of an ideally closed degassing A C C E P T E D M A N U S C R I P T systemratios: note, for example, the two intersections at [H 2 O TOT  2%; CO 2 TOT  4.7%] and [H 2 O TOT  2.3%; CO 2 TOT  1%], which clearly show CO 2 addition due to fluxing through the plumbing system

  et al., 2011), shows that plagioclase crystallization, for T  1100°C and 400 MPa, occurs under hydrous conditions (dissolved H 2 O  2.5 wt% ). Furthermore, dissolved H 2 O-CO 2 pairs from phase equilibria investigations agree with dissolved MI contents(Pichavant et al., 2009 and references therein), with H 2 O TOT -CO 2 TOT contents in the experimental charges close to those inferred from MIs. Such evidence corroborates the validity of information about volatile contents that can be extracted by MIs, at least those from the very volatile-rich active volcanoes of Southern Italy. The consistency of our water determinations, as well those of[START_REF] Spilliaert | Melt inclusion record of the conditions of ascent, degassing, and extrusion of volatile-rich alkali basalt during the powerful 2002 flank eruption of Mount Etna (Italy)[END_REF] with available independent petrologic hygrometers over the same pressure range[START_REF] Armienti | A new model to estimate deep-level magma ascent rates, with applications to Mt. Etna (Sicily, Italy)[END_REF][START_REF] Mollo | Reconstruction of magmatic variables governing recent Etnean eruptions: Constraints from mineral chemistry and P-T-fO 2 -H 2 O modeling[END_REF] further supports our determinations.On the contrary, 122 BC melt inclusion suggest the entrapment of a H 2 O-fluxed magma at shallow depth (P ≤ 150 MPa). Therefore, it is possible that the resulting H 2 O-enrichment of magma may result in the subsequent onset of diffusive phenomena throughout olivines and secondary water enrichment of hosted MIs[START_REF] Hartley | Diffusive over-hydration of olivine-hosted melt inclusions[END_REF]. However, we have no means to discriminate H 2 O-rich MIs, but of H 2 O-poor MIs. Post-entrapment phenomena would then be antithetic to the above mentioned results from melt-plagioclase equilibrium modeling, which instead criticize the representativeness of high-H 2 O MIs. In addition, if due to post-entrapment phenomena, the attainment of the buoyancy threshold at CO 2(g) /H 2 O (g)  5 (Fig.6a) would then be a fortuity, rather than the consequence of addition of CO 2 and dehydration, involving all active volcanic systems of the western Mediterranean region. On another side, recent studies have shown that CO 2 may be transferred from the melt or glass to a coexisting vapor bubble within the MI (e.g., Esposito et al. 2011), during post-entrapment crystallization (Steele-MacInnis et al. 2011) or differential thermal contraction (Moore et al. 2015 and references therein). First of all, CO 2 corrections to measured values are not necessary as PEC can be neglected (see Appendix). Second, it is not disputed that in a CO 2 -rich (fluxed) system such

	M A N U S C R I P T M A N U S C R I P T
	A C C E P T E D between the primary process (water fluxing itself) and the secondary one (post-entrapment diffusive H 2 O gain), anyway associated to water fluxing. Diffusive H 2 O-enrichment is likely to play a role in 122BC samples, subordinate to the occurrence of H 2 O-fluxing. On the other hand, we can safely discount such a process for recent products as the same samples cannot simultaneously experience and provide evidence of both CO 2 -and H 2 O-fluxing. Another argument against MIs representativeness involves post-entrapment disequilibrium processes (e.g., H-loss), which may offer an explanation alternative to ours to explain the observed A C C E P T E D
	spread of datapoints on H 2 O-CO 2 , H 2 O-S and H 2 O-Cl diagrams (Figures 6-8). For example, post-
	entrapment diffusive phenomena occurring in olivines may cause a fast and preferential loss of H 2 O
	(Danyushevsky et al. 2002; Gaetani et al. 2011), which then would particularly affect the low-H 2 O,

, which can be considered analogous to the Etnean magma for its total volatiles (H 2 O TOT  3.4 wt%, CO 2 TOT  2 wt%; Aiuppa et al., 2010; Pino high-CO 2 and high-S BN 97-98 MIs. Despite these phenomena can contribute to the dispersion displayed by the data (Figs. 3, 4, 6a, 7), they are not expected to cause dissolved H 2 O to span ~3 wt%. Moreover, water loss due to post-entrapment diffusion would not question about the validity of

Acknowledgements -All acquired data are given as supplementary material. Funding by PRIN-MIUR 2009 is acknowledged. The paper greatly benefited of the detailed and highly constructive review by Margaret Hartley (Univ. Manchester).

yellow ellipses) since a depth corresponding to ~400 MPa (at least). This produces a CO 2 cap on top of the 122BC system but gas dominated magma batches (large two-colours ellipses) in recent activity. In case of 122BC water fluxing (small cyan ellipses) in the shallow system overprints the previous CO-fluxing and is accompanied by diffuses crystallization (black dispersed chips).

Endogenus growth [START_REF] Allard | Endogenous magma degassing and storage at Mount Etna[END_REF] explains why the plumbing system for recent activity was pictured bigger than the 122BC one.

Appendix

A C C E P T E D M A N U S C R I P T

Etnean glasses is given in Figure A1 and shows that equation A1 reduces to

Carbon is present in studied inclusions only as carbonates CO 3 2- [START_REF] Holloway | Application of experimental results to C-O-H species in natural melts[END_REF][START_REF] Blank | Experimental studies of carbon dioxide in silicate melts: solubility, speciation and stable isotope behavior[END_REF]. The measurement is made through the peak at 1515 cm -1 , whose absorption coefficient,  1515 , has been determined according to [START_REF] Dixon | An experimental study of water and carbon dioxide solubilities in mid-ocean ridge basaltic liquids. part I: Calibration and solubility models[END_REF] as a function of the atom ratio Na/(Na+Ca). Because the carbonate doublet including the peak of interest is located on the shoulder of the large absorption domain of the silicate network, absorbance was estimated only after the subtraction of the spectrum of an identical but volatile-free glass. Thus, we subtracted to each spectrum both an oxidised and a reduced spectrum of reference, thus averaging the two absorbance values obtained through peak deconvolution. FTIR-estimated errors are 10% for H 2 O and at most 25% for CO 2 , by considering the intervening factors as independent variables.

A2. About Post Entrapment Crystallization (PEC)

PEC determined a compositional shift tending obviously to override the original melt-crystal equilibrium. Its effect is to magnify the differences in Mg, Na, and Si abundances between MIs and glass-matrix. Also note that the PEC shift does not explain the variability within a single olivine grain (see

Supplementary Data).

As a matter of fact, the compositional variations of major and also volatile elements of the Etna are quite independent of the host crystal chemistry, as evidenced by appreciable chemical differences recorded by MIs within a single homogeneous crystal grain.

A C C E P T E D M A N U S C R I P T

We have also computed Fe 3+ /Fe tot values on the basis of [START_REF] Roeder | Olivine-liquid equilibrium[END_REF], who defined a temperature independent coefficient for the Fe/Mg exchange between olivine and coexisting liquid basalt (K D =0.3). Fe 3+ /Fe tot values range from zero up to 0.62 in one BC sample (Figure A3a). This wide range is in line with the two log units of variations in fO 2 , related to water exsolution during magma ascent in the conduit and magma emplacement near the surface [START_REF] A C C E P T E D M A N U S C R I P T Moretti | On the oxidation state and volatile behaviour in multicomponent gas-melt equilibria[END_REF][START_REF] Burgisser | Redox evolution of a degassing magma rising to the surface[END_REF][START_REF] Mollo | Reconstruction of magmatic variables governing recent Etnean eruptions: Constraints from mineral chemistry and P-T-fO 2 -H 2 O modeling[END_REF].

We however considered the role that PEC may have occurred and thus affected the original redox state of iron. If PEC occurred, it should have increased the Fe 3+ /Fe tot ratio because of FeO removal from the liquid. We then assumed a low Fe 3+ /Fe tot value of 0.2, which is a lowermost value in our dataset common to the whole compositional range (Figure A3a) and which corresponds to XANES measurements by [START_REF] Bonnin-Mosbah | Iron oxidation states in silicate glass fragments and glass inclusions with a XANES micro-probe[END_REF]. We then back-calculated MI compositions by adding to each MI composition FeO, MgO, SiO 2 aliquots corresponding to the surrounding host olivine until the 0.2 iron redox ratio was attained. Computed values (Figure A3b) show that even under this -quite extreme-redox assumption, crystallization would be at most 7%, peaking to 11%

for only one MI from the 122BC. In our exercise the larger the difference in redox state, the larger the extent of PEC, which explains why the computed PEC shows the same trend increasing with olivine Fo content as for the iron redox state from Fe/Mg olivine-liquid partitioning (Figure A3).

Given the small extent of PEC, encompassed by EMPA errors, we discard post-entrapment effects and safely use raw analytical data. Considering the focus of our study (volatile contents and their interpretation) we can state that the PEC effect does not alter at all the solubility properties of the melt and its mass balance for volatiles.

Appendix figure caption

ACCEPTED MANUSCRIPT

A C C E P T E D M A N U S C R I P T