# Chlorine as a geobarometer for alkaline magmas: Evidence from a systematic study of the eruptions of Mount Somma-Vesuvius

Hélène Balcone-Boissard, G Boudon, R Cioni, J.D. Webster, G Zdanowicz, G Orsi, L. Civetta

► To cite this version:

Hélène Balcone-Boissard, G Boudon, R Cioni, J.D. Webster, G Zdanowicz, et al. Chlorine as a geobarometer for alkaline magmas: Evidence from a systematic study of the eruptions of Mount Somma-Vesuvius. Scientific Reports, 2016, 6, pp.21726. 10.1038/srep21726. hal-01278918

HAL Id: hal-01278918

https://hal.sorbonne-universite.fr/hal-01278918

Submitted on 25 Feb 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Chlorine as a geobarometer for alkaline magmas: Evidence from a systematic study of the eruptions of Mount Somma-Vesuvius

Received: 29 September 2015

Accepted: 29 January 2016

Published: 18 February 2016

H. Balcone-Boissard<fnref n="1" />, G. Boudon², R. Cioni³, J. D. Webster⁴, G. Zdanowicz²,⁵, G. Orsi⁵,⁶ & L. Civetta⁵,⁷

Defining the magma storage conditions of a volcanic system is a major goal in modern volcanology due to its direct implications for the style of a possible eruption, and thus on the associated risk of any crisis and the necessary management and mitigation strategies. Below 200 MPa and at equivalent depths, the strongly non-ideal behaviour of the H-C-O-S-Cl-F system in the silicate melt causes unmixing of the fluid phase to form an H₂O-rich vapour and a hydrosaline phase in equilibrium with the silicate melt, both responsible for buffering the chlorine (Cl) concentration. Following this equilibrium, the Cl concentration in melts may be used as a geobarometer for alkaline magmas. Systematic application of this method to the main explosive eruptions of Mount Somma-Vesuvius highlights two main magma ponding zones, at ~180–200 and ~100 MPa. At these pressures, the maximum pre-eruptive H₂O contents for the different magma compositions can be estimated; the results obtained, largely in agreement with the current literature, therefore confirm the validity of the method. The Cl geobarometer may help scientists to define the variation of the magmatic reservoir location through time and thus provide strong constraints on pre-eruptive conditions, which are of utmost importance for volcanic crisis management.

Experimental studies on the solubility of volatiles (H₂O, CO₂, S, Cl, F) in fluid-saturated silicate melts show that Cl preferentially partitions into aqueous fluids rather than silicate melts<fnref n="1" />⁻⁵. The NaCl-H₂O system is representative of magmatic fluids and is characterized by fluid immiscibility under a wide range of P-T conditions, and NaCl contents (P < 200 MPa, T < 1000 °C⁶⁻⁹; Fig. 1), and NaCl-H₂O fluids in the presence of silicate melt similarly display immiscibility. The P-T conditions noted above are commonly encountered where magmas are evolving in shallow reservoirs, and the associated immiscibility strongly affects the evolution of magmatic melts and fluids. Numerous experimental studies on Cl solubility have highlighted its saturation conditions in silicate melts³,⁴,⁶⁻¹⁴ especially in alkaline magmas; however, few studies with natural samples¹⁵ are presently available. Experimental results shed light on the importance of considering the more general melt-H-C-O-S-Cl-F pseudo-system, in order to describe the behaviour of the fluid phase and its relationships with the silicate melt. Recent experimental data¹⁶,¹⁷ point out the possible complex interactions among volatiles in the fluid phase, depending on mixing processes and volatile speciation, which in turn influence volatile solubility in a silicate melt. Even if H₂O is the dominant species and plays a prominent role in magma behaviour, mainly through decompressional exsolution¹⁸, CO₂¹⁹, Cl⁴,⁵, S¹⁶,²⁰,²¹, and, to a lesser extent, F¹⁶ also play a significant role in determining the conditions of fluid saturation of a silicate melt. The solubility of Cl in silicate melt largely depends upon melt composition²⁻⁴,²²⁻²⁴

<footnote n="1">¹Sorbonne Universités, UPMC Univ. Paris 06, CNRS, Institut des Sciences de la Terre de Paris (iSTeP), 4 place Jussieu 75005 Paris, France. ²Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Univ Paris Diderot, CNRS, F-75005 Paris, France. ³Dip.to Scienze della Terra, Universita' degli Studi di Firenze, Via La Pira 4, 50121 Firenze, Italy. ⁴Department of Earth and Planetary Sciences, American Museum of Natural History, Central Park West at 79th St., NY, NY 10024-5192 USA. ⁵Dipartimento di Scienze della Terra, dell'Ambiente e delle Risorse, Università degli Studi di Napoli Federico II, Largo S. Marcellino 10, 80138 Napoli, Italy. ⁶Dipartimento di Fisica "E. R. Caianiello" Università degli Studi di Salerno, Via Giovanni Paolo II 132, 84084 Fisciano, Salerno, Italy. ⁷Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo, Via Ugo La Malfa 153, 90146 Palermo, Italy. Correspondence and requests for materials should be addressed to H.B.-B. (email: helene.balcone_boissard@upmc.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Link between Cl behaviour and the geometry of a plumbing system. (A) the H₂O-NaCl-silicate melt pseudo-system, redrawn from phase equilibrium data⁶,⁷. Blue and purple domains: the silicate melt is in equilibrium with fluids, made of a vapour (H₂O and/or CO₂ and or S-rich) and a brine; the two colours correspond to two different melt temperatures. Within this subsolvus domain, drawn here for a given temperature, the Cl concentration in the silicate melt is buffered. Grey domain: the domain at higher P (T is fixed) than the solvus; the silicate melt is in equilibrium with a single fluid phase (see text for details).</figcaption></figure>

(B) Cross section of a plumbing system. Question marks: the position of the shallow crustal reservoir(s) will be constrained by Cl concentration in the silicate melt.

(being higher in alkali-rich magmas)¹³, pressure¹¹ and the co-occurrence of other volatile species, such as S²⁰,²¹. Temperature has little effect on Cl solubility⁴. In H₂O-bearing systems, the Cl concentration in the exsolved vapour may increase with increasing Cl concentration in the silicate melt. Yet, these systems exhibit strongly non-ideal mixing properties between H₂O, CO₂, and Cl. If, depending on P-T conditions, the crystallizing magma is at solvus or sub-solvus conditions (in the system H₂O-NaCl), the exsolved fluid phase unmixes to form a vapour and a brine, and so what was a pseudo-binary system becomes pseudo-ternary (Fig. 1). In such cases, the Cl partitioning behaviour is complex, as it is ruled by the coexistence of three phases in the volcanic system (vapour, silicate melt and brine). Because of this three-way partitioning, Cl exerts a strong control on fluid phase equilibria at the magmatic P-T conditions⁸,²⁵⁻²⁷, with the stability of the vapour and brine at sub-solvus conditions dependent on Cl activity in the fluids¹⁴,²⁷. In such a sub-critical domain at constant P and T, the Cl concentration in the silicate melt is invariant as the composition of both vapour phase and brine is fixed, due to Gibb's phase rule. This non-ideal behaviour of the fluid is expressed in the magmatic system by a buffered value of the Cl concentration in the silicate melt. The precise Cl buffer value varies with pressure: it increases at low pressure, as the immiscibility gap of the H₂O-salt system becomes wider, and decreases at higher pressure (Fig. 1). At magmatic T and shallow crustal P, the pseudo-ternary H₂O-NaCl-melt system predicts the presence of both vapour and brine, even for low NaCl contents¹⁷. During magma differentiation at hyper-solvus conditions, an increasing activity of Cl in the fluid may be ultimately responsible for the increase of the non-ideality of the mixing process among volatiles. Conversely, during differentiation at sub-solvus conditions, H₂O and Cl concentrations in the melt, vapour and brine will be fixed, while the vapour/brine ratio decreases with increasing Cl in the system¹⁷ (Fig. 2). The solubility relationships may be integrated to interpret the Cl content in terms of the pressure conditions at which the silicate melt is stored.

The aim of this paper is to demonstrate the power of Cl as a geobarometer for alkaline magmas, once the buffering effect of its concentration has been identified. The case study we use is based on the detailed study of material from several explosive eruptions, of variable magnitude, of Mount Somma-Vesuvius (MSV). For the purposes of geobarometry, we propose to use the systematic Cl buffering effect on the silicate melt to estimate the pressure, and hence depth, of the magmatic reservoir feeding a given eruption from solubility data (Fig. 3). The existence of such a buffering effect for the magmas feeding the studied eruptions has been systematically detected on the basis of the Cl content in the residual matrix glass. The Cl geobarometer may be used for a single eruption or, if sampling and chemistry permits, to reconstruct the behaviour of a volcanic system through time, as we do here for MSV. In addition, the maximum pre-eruptive H₂O content dissolved in a magmatic melt can be estimated using the corresponding experimental H₂O solubility law²⁸⁻³¹, based on the values of pressure of the
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

**Figure 2.** Cl behaviour and evidence of the buffering effect in silicate melt (melt inclusion and residual matrix glass composition). (A) Theoretical Cl behaviour. From the initial melt composition at a given Cl content (orange domain), melt differentiation under fluid-undersaturated conditions explains the Cl increase that is linearly correlated with the differentiation index (such as CaO or F) through the behaviour of an incompatible non-volatile component (solid grey line). When fluid saturation occurs (star), the pressure-temperature-composition of the silicate melt requires that the exsolved fluid phase unmixes into a low density, vapour phase and a Cl-rich brine or both phases exsolve directly from the melt simultaneously. Consequently, Cl partitions (dashed yellow line) between fluid phases and silicate melt to respect the Gibb's phase law. Thus, the Cl concentrations in the silicate melt, vapour and brine are fixed at a buffered concentration corresponding to the thermodynamic equilibrium of the system (solid violet line); (B) Example data from melt from the Pollena eruption of MSV: melt inclusions (green diamonds) and residual glass (pink triangles) are reported and interpreted respect to the model as in 2a.

<figure><img src="image_6.png" /></figure>

**Figure 3.** Location of the pressure domain of the magma reservoir estimated from the Cl buffer value.
Experimental Cl solubility law in Na-Phonolitic¹¹, K-Phonolitic¹¹ and Trachytic¹² melt saturated with aqueous fluid + brine are represented versus pressure. From the Cl buffering concentration in silicate melt, by using the relevant experimental Cl solubility law, it's possible to estimate the pressure domain or domains at which the magma is stored at depth. The uncertainty of the pressure estimates depends on the Cl buffering concentration (mean values obtained on several pumice clasts by EPMA) and on the experimental Cl solubility law. Arrow (1): from the Cl buffer value until the experimental Cl solubility law; Arrow (2): from the experimental Cl solubility law to the pressure. Solid line: the corresponding conversion of the Cl buffer value. Dashed lines: the uncertainty from the Cl concentration determined by EPMA, that constrains the precision on the pressure domain estimate.

magma storage derived from the Cl buffer (Fig. 4). As for Cl, the H₂O solubility model should take into account the influence of other volatiles, such as CO₂¹⁷,³²,³³, as proposed in Cl-H₂O solubility models¹⁷.

**Terminology**

Hereafter, we shall take *silicate melt* to refer to a molten silicate magmatic body that evolves by fractional crystallization in a reservoir within the crust. *Fluid* is a general term to mean a multicomponent, non-crystalline phase. In volcanic environments the volatile components exsolve as a *vapour phase*, which is a low-density fluid phase dominated by H₂O and/or CO₂ with more or less S species (SO₂, H₂S), and a *brine*, or *hypersaline liquid*, being a halogen (Cl and F as halide) species-rich, water-bearing fluid of magmatic origin. Silicate melt may be in equilibrium with either a vapour or a brine, or both. *Volatiles* is a general term referring to all species that may strongly partition in favour of a fluid phase. The *solubility* of a given volatile species in a silicate melt is the amount of that
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Pre-eruptive H₂O content estimates from Cl buffer value. How to convert a Cl buffer value into a pre-eruptive dissolved H₂O content of a melt? (A) From the domain of pressure deduced from the Cl buffer value, the H₂O concentration in the silicate melt is deduced by reading the H₂O content at this pressure using the relevant experimental H₂O solubility law. Solid line: the corresponding conversion of the mean Cl or H₂O value. Dashed lines: the uncertainty for the estimated Cl or H₂O concentration. Theoretical Cl (blue) and H₂O (red) experimental solubility laws are reported. (B) Experimental H₂O solubility law for Na-phonolite²⁹, K-phonolite²⁸, trachyte³⁰ and tephrite³¹ are reported. As for Cl, the key point is to have the H₂O solubility laws that correspond to the studied melt composition; this constrains the uncertainty of the method as the Cl content is accurately measured by EPMA and the analytical contribution to uncertainty is relatively minor.</figcaption></figure>

species dissolvable into that melt, at given P-T conditions, when the melt is saturated in either vapour (for H₂O solubility), or brine (for Cl solubility) or both (for both H₂O and Cl solubility). It also depends upon the mechanism through which it is incorporated in the silicate network, i.e. its speciation in the melt phase.

## How to use Cl as a geobarometer?

To infer magma storage conditions from Cl concentrations detected in silicate glass, it is necessary to take into account (1) the influence of the silicate melt composition and (2) the effect of the fluid composition on Cl behaviour.

**The composition of the melt.** The influence of the evolutionary path of a silicate melt on Cl solubility was first suspected on the basis of the results of experimental studies that shed light on the role of the (Na + K)/Al and Al/Si molar ratios³,¹¹, and the alkali⁷,⁹ and CaO contents²³. Modelling of the predicted Cl content in silicate melt, extrapolated from the results of solubility experiments, demonstrates the major role played by the silicate melt composition on Cl solubility¹³. Modelling also highlights the influence of each major cation on Cl solubility. Such multi-component considerations allow us to take into account the significant changes in Cl behaviour resulting from fractional crystallisation of a silicate melt, even for an apparently subtle variation of its composition. During crystallisation, the Cl concentration in a fluid-free silicate melt increases as it behaves as an incompatible element (Fig. 2). Simultaneous decreases in Mg, Ca, Fe, and Al content may be responsible for the direct exsolution of a brine¹⁷. Thus, melt differentiation with accompanying increase of the Cl concentration in residual melt, leads to brine exsolution, with or without vapour, which may be seen as a particular form of second boiling¹⁷.

**The effect of other volatile species.** The behaviour of the volatiles in a fluid composed of only one phase (vapour or brine), differs from that in a fluid composed of coexisting immiscible vapour and brine. At sub-solvus conditions, when applied to the melt composition of interest, Cl modelling can provide the T-P conditions that best reproduce the Cl content measured in the same glass. Importantly, Cl solubility is also strongly affected by the presence of other volatile species, such as S. In fact, Cl solubility decreases by up to 30% in S-rich silicate melts¹³,¹⁶,²⁰,²¹,³⁴, so that Cl is observed to be most abundant in S-poor melts¹⁶. As S substantially reduces the solubility of Cl in silicate melts under oxidizing conditions¹⁷, and given the comparatively high oxygen fugacity and S contents of MSV magmas, the S effect in these magmas is likely to be significant. For this reason, the S effect has to be taken into account during Cl solubility modelling. The role played by F in terms of Cl solubility is less clear, although Cl solubility likely increases with increasing F¹³. Variation of the characteristics (i.e., compositions and phase relations) of the fluids due to the addition of other volatile components may also affect the Cl concentration in silicate melts (Fig. 1). For instance, an increase in CO₂ will lead to an increase in the non-ideality of the volatile mixing processes¹⁶. In the presence of both vapour and brine, CO₂ expands the stability field of the two fluids³⁵, whereas S reduces this field, thus favouring mixing¹⁶. Such competing effects need to therefore be modelled to accurately determine Cl solubility.

**Solubility and equilibrium.** The equilibrium between silicate melt and fluids is generally inherited from conditions established in the reservoir rather than during magma ascent¹⁸. The rate of magma ascent during an eruption controls the degree to which equilibrium conditions are attained between silicate melt, vapour phase and brine at a given P. In the case of a sufficient fast ascent rate, a silicate melt may have insufficient time to re-equilibrate before fragmentation, thus preserving a record of the pre-eruptive conditions. Therefore, the
<figure><img src="image_8.png" /></figure>

**Figure 5.** The main explosive eruptions of Mount Somma-Vesuvius (A) Schematic chronogram of Mount Somma-Vesuvius activity for studied explosive eruptions (modified from reference³⁷). Length and colour of arrows refer to the estimated VEI. Breaks in the chronogram mark changes of timescale. (B) Residual glass compositions are reported in a total alkalis-silica diagram. These data refer to mean glass compositions of samples (at least 10 analyses per pumice), from the density mode of each eruptive unit described for each eruption, in which the Cl buffering effect has been identified (Table S1); except for Pompeii where the mean composition for each eruptive unit is represented.

Cl-buffering effect can be assessed through the analysis of residual matrix glass, representing the melt quenched at fragmentation (Fig. 2). Due to the lower diffusivity of Cl with respect to that of H₂O in melts³⁶, syn-eruptive, degassing-induced crystallisation of microlites during magma ascent³⁷ could increase the Cl concentration in the residual melt. Therefore, the Cl concentration measured in the residual glass has to be corrected for the micro-lite effect in order to determine the pre-eruptive Cl buffer value. An alternative process which could induce Cl buffering is the saturation in a Cl-bearing crystalline phase (e.g. feldspathoids of the cancrinite-sodalite group, commonly present in the K-alkaline magmas of MSV). If Cl-bearing minerals do not represent a significant proportion of the crystallizing mineral assemblage, it is unlikely that they can induce a buffering effect on Cl concentration in the silicate melt at non-equilibrium conditions.

Using the Cl solubility law for either a given silicate melt composition, or the law for the most similar melt composition available, the pressure conditions of magma storage can be deduced by means of the determined Cl buffer value (Fig. 3). Unfortunately, such experimental data are not available for all alkaline melt compositions. As Cl solubility largely depends upon melt composition, the pressure of the magma reservoir estimated by using a Cl experimental solubility law determined on a silicate melt of different composition (in terms of alkali ratio, or Ca, Mg, Fe content) yields an additional and significant uncertainty. In order to avoid such uncertainties, it is preferable to use a more general Cl solubility model which takes into account the whole silicate melt composition, such as the one developed by Webster and collaborators¹³,¹⁷.

Melt inclusions (MIs) can also provide interesting information as they represent droplets of silicate melt trapped at different stages during the differentiation process³⁸. MIs trapped in phenocrysts highlight Cl incompatible behaviour during fractional crystallisation until the point of fluid saturation (Fig. 2). Those trapped after fluid saturation, in the subcritical domain in which vapour and brine are in equilibrium with the silicate melt, record the pre-eruptive achievement of Cl saturation in the magma reservoir (Fig. 2).

## Results

We analysed the products of 12 explosive eruptions of MSV³⁷,³⁹,⁴⁰ (Fig. 5), including four Plinian (Pomici di Base, Mercato, Avellino, Pompeii), five sub-Plinian (Greenish, AP1, AP2, Pollena, 1631 AD) and three violent Strombolian events of the recent series (1822, 1906, 1944)³⁹,⁴¹. All magmas emitted during the MSV activity were alkali-rich melts that had undergone complex processes of differentiation, including mixing, assimilation, and fluid exsolution⁴⁰,⁴². The residual matrix glass compositions range from tephritic to phonolitic and trachytic (Fig. 5). Many of these deposits show a variable compositional zoning, with a general decrease in melt evolution from the earliest to the latest erupted products. Such compositional zoning has been interpreted as the result of progressive tapping of a compositionally zoned reservoir⁴⁰,⁴²,⁴³, accompanied by complex syn-eruptive magma
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Cl vs. F contents in residual glass from the first eruptive units of each explosive eruption of Mount Somma-Vesuvius in which the Cl buffer effect has been identified. One point represents a mean Cl content of one single fragment describing the density mode (at least 10 analyses); except for Pompeii (melt inclusions + mean Cl content for each eruptive unit). Data for Pompeii⁵⁸ and Avellino⁵⁹ are also reported. Melt inclusion data are not represented except for Pompeii (Table S2; see discussion in the text). (A) Pomici di Base, Greenish, Pollena and 1631; (B) Mercato, Avellino, AP1, 2 and Pompeii; (C) 1822, 1906, 1944. Cl uncertainty is within symbol size, contrary to F (see symbol).</figcaption></figure>

mixing processes⁴². We have focused our research on the earliest erupted, most-evolved products of each event, likely representing the uppermost, fluid-saturated portion of the reservoir. Both residual matrix glass (Table S1, S2) and MI (Table S3) compositions were studied. MIs register a more complex history by exhibiting a wide range of Cl content, with values from within fluid-undersaturated conditions to within fluid-saturated conditions similar to those measured in the residual glass (Fig. 2). Here, we focus only on the residual matrix glass composition, as our interest is in the shallow plumbing system. In such glasses of the products of all eruptions, the F concentration positively correlates with the degree of melt differentiation, thus suggesting its overall behaviour as an incompatible element with minimal partitioning into a fluid. We note that CaO may also be used as a differentiation index in alkali magmas. On the contrary, the Cl content is constant, regardless of the evolving melt composition (Fig. 6; Tables S1, S2). As Cl-rich minerals are rare in the investigated rocks, the buffered Cl values, corrected for microlite crystallization, must be a result of the silicate melt being in equilibrium with a two-phase fluid. This is because the presence of a two phase fluids is the only other phase assemblage that can fix the Cl content in the system. Microlite crystallization induced by syn-eruptive degassing clearly occurred and varied from 0 to 30 vol% (Table S1), depending on the eruption.

The trachytic melt involved during the earliest-erupted Pomici di Base (22,030 yr BP) and Greenish (19,065 yr BP) eruptions, show buffered Cl contents of 6,764 ± 180 and 5,534 ± 290 ppm, respectively (Fig. 6a). The phonolitic magma of the Mercato (8,890 yr BP) eruption displays a value of 6,306 ± 320 ppm, slightly higher than those of the 5,200 ± 400 and 5,368 ± 380 ppm displayed by the other large Plinian eruptions of Avellino (4,365 yr BP) and Pompeii (AD 79), respectively (Fig. 6b). The phonolitic tephrite to phonolite melts feeding the two main AP eruptions at 3,500 yr BP (AP1 and AP2³⁷), occurring shortly after the major Plinian Avellino eruption, share close Cl values of 5,475 ± 120 and 5,196 ± 490 ppm, respectively (Fig. 6b). The phonolitic tephrite melt from the
Pollena (AD 472) and 1631 eruptions, the two latest sub-Plinian events, display Cl values restricted to $6,907 \pm 270$ and $6,688 \pm 280$ ppm, respectively (Fig. 6b). The tephritic melt of the more recent and violent 1822, 1906 and 1944 Strombolian eruptions also highlight the same Cl buffering effect, at values of $7,524 \pm 230$, $6,226 \pm 190$, and $5,022 \pm 218$ ppm, respectively (Fig. 6c).

Discussion

The Cl solubility in silicate melt⁵, mainly dependent upon composition and P (negligible T effect), is a key parameter to convert the Cl buffer value detected in residual matrix glass into a P value indicative of the depth of the shallow magma reservoir feeding the eruption. The Cl solubility in a given melt composition can be estimated through two different methods. One method is based on the use of the experimentally determined solubility law defined on the same magma composition¹¹,¹², hereafter referred to as “experimental Cl solubility”, while the other on the use of modelling¹³,¹⁷, hereafter referred to as “modelled Cl solubility”. In the following we will discuss and compare the results for the studied MSV eruptions by using both methods, which will demonstrate the very strong dependence of Cl solubility on melt composition and will reveal the pre-eruptive pressure conditions at which Cl buffer was reached for the different cases analysed.

As the pre-eruptive Cl content in the different MSV melts varies in a bracket of $\pm 200-300$ ppm, equilibrium pre-eruptive pressures estimated with the experimental Cl solubility method for each eruption, yield a small uncertainty (around 10 MPa; Fig. 3). Among the different MSV erupted magmas, Cl solubility has been experimentally measured only for the trachytic Pomici di Base¹² and the K-phonolitic Pompeii¹¹, melts. For these eruptions, the pressure domain estimated from the measured Cl buffering value in the natural samples is highly precise, and gives values of $100 \pm 10$ and $185 \pm 10$ MPa, respectively. Cl solubility curves are not available for the compositions of the magmas feeding the other studied eruptions. Using the results of experiments on the Na-phonolitic rocks of Montaña Blanca (Canary Islands)¹², that are compositionally similar to the residual glasses of the Mercato, Avellino, AP 1, and AP 2 products, we have defined pressure domains for the magma reservoirs feeding these four eruptions of $170 \pm 10$, $200 \pm 10$, $195 \pm 10$ and $200 \pm 10$ MPa, respectively. For the Greenish eruption glass, that straddles the phonolite and trachyte fields (Fig. 5), we have estimated a large pressure domain of $150-200$ MPa on the basis of the Cl solubility laws available for both trachytic and K-phonolitic melts. Cl solubility for the Pompeii K-phonolite melt has also been used for the mildly evolved glasses of both Pollena and 1631 eruptions, estimating equilibrium pressures of $105 \pm 10$ and $115 \pm 15$ MPa, respectively. For the tephritic glasses of the 1822, 1906 and 1944 eruptions, the closest available Cl solubility law, although poorly constrained, derives from experiments performed on a basaltic composition⁴⁴. The magmas of these three eruptions show different Cl buffering values, but the same equilibrium pressure value of ~105 MPa, as imposed by the very few experimental data.

A modelled Cl solubility method, such as the one developed by Webster and collaborators¹⁶,¹⁷ (Fig. 7) takes into account composition (in particular including S content), P, and T of the investigated melt. In our study we have used the composition of each MI and the mean of the compositions of the residual matrix glass detected in variable fragments (Tables S2, S3; Fig. 7). No systematic S measurements have been performed, but literature data provide considerable constraints on S concentrations in melt. The residual matrix glasses of the large plinian eruptions (Mercato, Avellino, Pompeii) exhibit S concentrations between $200-300$ ppm⁴⁵, and those of the sub-plinian eruptions of Pollena⁴⁶ and 1631⁴⁵ between 1800-2000 and 1000-1200 ppm, respectively (EPMA data). Poor constraints exist for the Pomice di Base and Greenish eruptions that give S estimates of 500-1300 and 200-2800 ppm, respectively (EDS data³⁹). For the last eruptions of 1822, 1906 and 1944, that involve less-differentiated melts, melt inclusions in olivines provide S contents of 1800-1900 ppm for 1906 and 1000-2400 ppm for 1944⁴⁷. In addition, the K-tephritic to K-basanite magma that is thought to supply the superficial reservoir has S contents between 1500-3000 ppm⁴²,⁴⁷. The Cl saturation pressure values for the Pomici di Base and Pompeii melts, estimated by using the modelled Cl solubility taking into account S effect and the experimental Cl solubility methods (Fig. 7), are in very good agreement (Figs 3 and 6a). For the Mercato, Avellino, AP 1 and AP 2 melts, the pressure values estimated using the Cl solubility law for Na-phonolite and those derived from Cl solubility modelling taking into account S effect are in the same range. For the magmas feeding the Greenish, Pollena, 1631, 1822, 1906, and 1944 eruptions, for which experimental solubility data are not available, Cl solubility modelling (corrected for the S effect) provides better constrained pressure estimates. For the Greenish feeding reservoir, a pressure of 100 MPa gives the best fit (Fig. 7). The latest two sub-Plinian eruptions of Pollena and 1631 also highlight an equilibrium pressure for their magma reservoirs of ~100 MPa. Cl modelling requires lower equilibrium pressures for the magmas extruded during the latest three studied eruptions, of ~80 MPa for those occurred in 1822 and 1906, and ~50 MPa for the 1944 event.

Whatever the method used to estimate the Cl solubility, experimental or modelled, the pressure values we have obtained by using the buffered Cl contents of melts as a geobarometer, are in very good agreement with published data using other pressure constraints for the same MSV eruptions. The pressure domain estimated from the Cl buffering effect for the reservoirs feeding the Mercato, Avellino and Pompeii Plinian eruptions is similar to that arising from experimental petrology ($200 \pm 20$ MPa)⁴⁸. The value of ~100 MPa that we have estimated for the magma feeding the Pollena eruption is in agreement with the results of both experimental petrology ($100 \pm 20$ MPa)⁴⁸ and MI studies (~95 MPa)¹⁵,⁴⁶. Also, our pressure estimates for the reservoirs feeding the post-Pompeii eruptions are in good agreement with the results of both experimental petrology and MI studies⁴⁵⁻⁴⁸. The pressure value of the magmatic reservoir of the Pomici di Base eruption deduced from the Cl-buffering effect ($100 \pm 10$ MPa) strongly contrasts with that obtained on feldspar geobarometers ($300-400$ MPa)⁴⁹ while it is coherent with the estimates resulting from water fugacity-based studies (below 165 MPa)⁵⁰. Modelling Cl concentrations in the silicate melt at 300-400 MPa gives unrealistic values. The similar pressure value obtained for the reservoirs feeding the Pomici di Base and the Greenish eruptions³⁷,³⁹ (Table S1), that are very close in magma composition and dynamics, indirectly confirms the suitability of our
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Modelled and measured Cl concentrations vs. CaO in glass (residual glass and melt inclusions) from two eruptions of MSV, Pomice di Base and Greenish eruptions. (A) Measured Cl concentration in residual glass (mean values of at least 10 analyses) and melt inclusions (single measurement). (B) Measured Cl concentration compared to modelled Cl content without S correction, at a pressure of 100 MPa. (C) Measured Cl concentration compared to modelled Cl content with S correction (reduction of 30% in relative of Cl solubility<sup>17</sup>), at the same pressure of B (100 MPa).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 8.** The evolution through time of the shallow plumbing system for the last 22 ky period of activity of Mount Somma-Vesuvius. The uncertainty in the pressure determination is within symbol size (see text for explanation; Table S1 and S4).</figcaption></figure>

estimate. Neither direct experimental petrology nor field (e.g. nature of accidental lithic fragments) data that could contribute to assessments of the depth of the reservoirs feeding these two eruptions are available.

The volatile solubility data presented here depict a time and space evolution of the shallow MSV plumbing system (Fig. 8) that is slightly different for that previously discussed by other authors<sup>37,39,45,48</sup>. This evolution summarizes the discussion on the choice of the best method to estimate the pressure of the shallow reservoir to establish this time and space evolution (see above and Table S4). In fact, basing on the new P estimated, the magmas emitted during the earliest high-magnitude eruptions (Pomici di Base and Greenish) were stored in a shallow reservoir located at a depth corresponding to a pressure of ~100 MPa, in contrast with what previously suggested<sup>49</sup>. A deepening of the magma storage to pressure conditions between ~180 and 200 MPa began with the formation of the reservoir of the Mercato eruption (between 19,265 and 8,890 yrs BP) and lasted at least until the Pompeii eruption (79 AD). From the Pollena (472 AD) until the 1631 eruption, the magma storage depth
switched again to a shallower level (corresponding to a pressure of ~100 MPa), similar to that of the Pomici di Base event. Subsequently, it progressively raised towards more and more shallow depth, up to the one corresponding to ~50 MPa, that fed the less intense, violent Strombolian eruptions of 1822, 1906 and 1944. The post-9 ka evolution of the magmatic system confirms the previous models for this period of activity⁴⁸. On the basis of our results, we suggest that the MSV volcanic system has been characterised through time by two main zones of magma storage, with the roof of the reservoirs located at ~180–200 (~7–8 km depth) and ~100 MPa (~3–4 km depth). The lower pressure values are consistent with petrologic constraints³⁸,⁵¹,⁵². The higher pressure values, in good agreement with experimental petrology data⁴⁸, are also consistent with the results of seismic surveys, suggesting the presence of a still active large reservoir with the top at ~8 km depth⁵³. Finally, recent electrical resistivity investigations suggest the presence of residual cooling magma batches at 3–4 km depth, in addition to the deeper one⁵⁴.

As a corollary to these new and robust constraints on the depth of the shallow plumbing system, determined by the recognition of the Cl buffering effect on the MSV products, a first estimate of the maximum pre-eruptive H₂O content in all these silicate melts can also be derived (Fig. 4; Table S1). We estimate H₂O concentrations using a H₂O solubility law in the corresponding melt at the pressure deduced from the Cl-buffering effect. These first estimates give between 4.8 and 7 wt% H₂O for the phonolitic and trachytic magmas, and 3.5 and 4.8 wt% for the less differentiated tephritic and phonotephritic melts. This approach requires that the H₂O solubility law for the composition of the erupted magma is known, and it neglects the influences of other volatiles on H₂O solubility (mainly CO₂). Since CO₂ is not considered into the H₂O solubility law, and as CO₂ may influence H₂O solubility³⁵, these estimates are first estimate (maximum values), as in particular the presence of CO₂ could possibly depress the amount of water dissolved in glass at saturation, by dilution of the fluid at low bulk Cl concentration³⁵, introducing consequently large error on H₂O estimates. However, volatile studies on the products of the AD 79 Pompeii eruption confirm such values, also revealing the absence of an important CO₂ component⁵⁵. It is also likely that CO₂ fluxing through the shallow magma chambers represents a significant process at MSV, as suggested by oxygen isotope data on olivine and diopsidic clinopyroxene⁵⁶. This process could be responsible for an early establishment of volatile saturation in the magma, due to the decrease in H₂O solubility induced by the changes in fugacity of different volatile species in the fluid phase as CO₂ increases. This effect could be relevant especially in small magma chambers, such as those feeding the post-Pompeii mafic activity³¹,⁴⁵. Furthermore, it could also explain the achievement of Cl buffering for such mafic compositions, as H₂O would be forced to exsolve in response to CO₂ fluxing¹⁶,⁵⁶. Compared to H₂O contents measured directly in melt inclusions, the minimum exsolved H₂O content during magma ascent may be estimated. In several cases, the H₂O contents measured in melt inclusions are consistent with our estimates (Table S1).

## Conclusion

A reliable assessment of volcanic hazard requires the time, place, and scenario of the next eruption to be defined. Facing this issue, the Cl content of evolved, Cl-enriched alkaline melts represented by matrix glasses may be used as (1) a geobarometer for locating a crustal reservoir and (2) an indicator of H₂O dissolved content in silicate melt prior to eruption. The Cl geobarometer here presented may be useful for a large composition domain of alkaline magmas, such as rhyolite, although this was not included in our study¹⁵. When systematically applied to the products of the past activity of a given volcano, this geobarometer can also constrain the temporal evolution of its plumbing system. All these parameters are of interest for alkali magmas as they may allow a better understanding of storage conditions (pressure, depth and pre-eruptive volatile content) that are linked to eruptive style as highlighted for phonolitic-trachytic magmas⁵⁷. Our results also highlight that fluid phase relations are key parameters to explain the H₂O-CO₂-S-Cl dissolution in a silicate melt at relevant magmatic pressures. Unfortunately, not all these relationships, such as the influence of CO₂, S or F, are yet clearly defined. This also has implications for acquiring a deeper understanding of halogen fractionation in fluid-saturated magmas.

## Methodology

For the products of each investigated eruption, we selected one representative outcrop including the basal fallout sequence of the entire pyroclastic succession and collected samples from each major eruptive unit. Samples from the deposits of pyroclastic flow phases were not included in our study. The sampling strategy was based on the collection of at least 100 pumice clasts per eruptive unit. Each clast was then sawed into three parts, and a density measurement performed on one third of each of them by the "three weights" method⁵⁸,⁵⁹. Several pumice clasts representative of the eruptive unit were chosen in the mode of the density distribution, that is among those with the most representative density value (generally ~80% of the total erupted magma). Electron microprobe (EPMA) analyses (Cameca SX 100, Camparis - France) were performed on one third of the selected pumice⁵⁸,⁵⁹, with an acceleration voltage of 15 kV and a beam current of 4 nA. The dwell time was 10 s, except for Si (5 s) and Na (5 s). F and Cl were measured with a dwell time of 60 s, an acceleration voltage of 15 kV and a beam current of 40 nA. Analytical conditions (in particular dwell times) were established to minimize Na and F diffusion during EPMA analyses using glass standards. For inter-calibration of EPMA sessions, three natural obsidians were analysed⁵⁸,⁵⁹. Microlite contents of residual glass may be estimated by chemical mapping⁵⁹. MI studies were performed on phenocrysts separated from several clasts per each studied eruptive unit. In order to have a consistent dataset, only MIs hosted in phenocrysts with diameter >125 µm, were chosen.

The Cl modelling was based on the model published by Webster and collaborators¹³ expressing the influence of major elements on Cl solubility through the definition of association coefficients linking Cl to cations in melt. Their Excel™ spreadsheet allows one to introduce and account for each matrix glass and MI composition determined by EPMA measurements (Si, Ti, Al, Fe (as FeO), Mg, Mn, Ca, Na, K, F) at a given temperature-pressure condition. The result of the computation is the modelled Cl concentration in the silicate melt. The S effect is introduced manually by reducing the Cl contents of the melts by 30 relative percent, after calculating the Cl solubility
for the correct bulk melt composition, P, and T. This option is not routine in the version of the Excel™ spreadsheet because it would need more systematically run experiments with S for other melts than phonolitic in composition. This estimate of S effect is based on experiments performed by Webster and collaborators¹⁶,¹⁷ on phonolitic melt with S content (2000 ppm) in the range of that of MSV magmas.

## References

1. Kilinc, I. a. & Burnham, C. W. Partitioning of chloride between a silicate melt and coexisting aqueous phase from 2 to 8 kilobars 1. *Econ. Geol.* **67**, 231–235 (1972).

2. Kravchuk, I. F. & Keppler, H. Distribution of Chloride Between Aqueous Fluids and Felsic Melts At 2-Kbar and 800-Degrees-C. *Eur. J. Mineral.* **6**, 913–923 (1994).

3. Metrich, N. & Rutherford, M. Experimental study of chlorine behavior in hydrous silicic melts. *Geochim. Cosmochim. Acta* **56**, 607–616 (1992).

4. Webster, J. D. Chloride Solubility in Felsic Melts and the Role of Chloride in Magmatic Degassing. *J. Petrol.* **38**, 1793–1807 (1997).

5. Webster, J. D. The exsolution of magmatic hydrosaline chloride liquids. *Chem. Geol.* **210**, 33–48 (2004).

6. Anderko, A. & Pitzer, K. S. Equation-of-state representation of phase equilibria and volumetric properties of the system NaCl-H₂O above 573 K. *Geochim. Cosmochim. Acta* **57**, 1657–1680 (1993).

7. Anderko, A. & Pitzer, K. S. Phase-Equilibria and Volumetric Properties of the Systems KCl-H₂O and NaCl-KCl-H₂O above 573-K - Equation of State Representation. *Geochim. Cosmochim. Acta* **57**, 4885–4897 (1993).

8. Sourirajan, S. & Kennedy, G. C. The system H₂O-NaCl at elevated temperatures and pressures. *Am. J. Sci.* **260**, 115–141 (1962).

9. Chou, I. Phase relations in the system NaCl-KCl-H₂O. III: Solubilities of halite in vapor-saturated liquids above 445 °C and redetermination of phase equilibrium properties in the system NaCl-H₂O to 1000 °C and 1500 bars. *Geochimica Cosmochimica Acta* **51**, 1965–1975 (1987).

10. Shinohara, H. Exsolution of immiscible vapor and liquid phases from a crystallizing silicate melt: Implications for chlorine and metal transport. *Geochim. Cosmochim. Acta* **58**, 5215–5221 (1994).

11. Signorelli, S. & Carroll, M. R. Solubility and fluid-melt partitioning of Cl in hydrous phonolitic melts. *Geochim. Cosmochim. Acta* **64**, 2851–2862 (2000).

12. Signorelli, S. & Carroll, M. R. Experimental study of Cl solubility in hydrous alkaline melts: Constraints on the theoretical maximum amount of Cl in trachytic and phonolitic melts. *Contrib. to Mineral. Petrol.* **143**, 209–218 (2002).

13. Webster, J. D. & De Vivo, B. Experimental and modeled solubilities of chlorine in aluminosilicate melts, consequences of magma evolution, and implications for exsolution of hydrous chloride melt at Mt. Somma-Vesuvius. *Am. Mineral.* **87**, 1046–1061 (2002).

14. Webster, J. D. & Mandeville, C. W. Fluid immiscibility in volcanic systems. In: Leibscher, A., Heinrich, C. (Eds.), Fluid-Fluid Equilibria in the Crust. *Rev. Mineral. Geochemistry* **65**, 313–362 (2007).

15. Lowenstern, J. B. Chlorine, fluid immiscibility, and degassing in peralkaline magmas from Pantelleria, Italy. *Am. Mineral.* **79**, 353–369 (1994).

16. Webster, J. D., Goldoff, B., Sintoni, M. F., Shimizu, N. & De Vivo, B. C-O-H-Cl-S-F Volatile Solubilities, Partitioning, and Mixing in Phonolitic-Trachytic Melts and Aqueous-Carbonic Vapor Saline Liquid at 200 MPa. *J. Petrol.* **55**, 2217–2248 (2014).

17. Webster, J. D. et al. Experimental and modeled chlorine solubilities in aluminosilicate melts at 1 to 7000 bars and 700 to 1250 °C: Applications to magmas of Augustine Volcano, Alaska. *Am. Mineral.* **100**, 522–535 (2015).

18. Sparks, S. R. J. Dynamics of magma degassing. *Geol. Soc. London, Spec. Publ.* **213**, 5–22 (2003).

19. Joyce, D. B. & Holloway, J. R. An experimental determination of the thermodynamic properties of H₂O-CO₂-NaCl fluids at high pressures and temperatures. *Geochim. Cosmochim. Acta* **57**, 733–746 ST-An experimental determination of the (1993).

20. Botcharnikov, R. E., Behrens, H., Holtz, F., Koepke, J. & Sato, H. Sulfur and chlorine solubility in Mt. Unzen rhyodacitic melt at 850 °C and 200 MPa. *Chem. Geol.* **213**, 207–225 (2004).

21. Webster, J. D., Sintoni, M. F. & De Vivo, B. The partitioning behavior of Cl, S, and H₂O in aqueous vapor-±saline-liquid saturated phonolitic and trachytic melts at 200 MPa. *Chem. Geol.* **263**, 19–36 (2009).

22. Webster, J. D. Fluid-melt interactions involving Cl-rich granites: Experimental study from 2 to 8 kbar. *Geochim. Cosmochim. Acta* **56**, 659–678 (1992).

23. Chevychelov, V. Y. Chlorine dissolution in fluid-rich granitic melts: The effect of calcium addition. *Geochemistry Int.* **37**, 522–535 (1999).

24. Carroll, M. R. Chlorine solubility in evolved alkaline magmas. *Ann. Geophys.* **48**, 619–631 (2005).

25. Bodnar, R. J. & Sterner, S. M. Synthetic Fluid Inclusions. *Hydrothermal Exp. Tech.* (eds Ulmer G.C., Barnes H.L.) Wiley Interscience. (1987).

26. Driesner, T. & Heinrich, C. A. The system H₂O-NaCl. Part I: Correlation formulae for phase relations in temperature-pressure-composition space from 0 to 1000 °C, 0 to 5000 bar, and 0 to 1 XNaCl. *Geochim. Cosmochim. Acta* **71**, 4880–4901 (2007).

27. Liebscher, A. & Heinrich, C. A. Fluid-fluid interactions in the Earth's Lithosphere. *Rev. Mineral. Geochemistry* **65**, 1–13 (2007).

28. Iacono Marziano, G., Schmidt, B. C. & Dolfi, D. Equilibrium and disequilibrium degassing of a phonolitic melt (Vesuvius AD 79 'white pumice') simulated by decompression experiments. *J. Volcanol. Geotherm. Res.* **161**, 151–164 (2007).

29. Schmidt, B. C. & Behrens, H. Water solubility in phonolite melts: Influence of melt composition and temperature. *Chem. Geol.* **256**, 258–267 (2008).

30. Di Matteo, V., Carroll, M. R., Behrens, H., Vetere, F. & Brooker, R. A. Water solubility in trachytic melts. *Chem. Geol.* **213**, 187–196 (2004).

31. Lesne, P., Scaillet, B., Pichavant, M., Iacono-Marziano, G. & Beny, J.-M. The H₂O solubility of alkali basaltic melts: an experimental study. *Contrib. to Mineral. Petrol.* **162**, 133–151 (2011).

32. Witham, F. et al. SolEx: A model for mixed COHSCI-volatile solubilities and exsolved gas compositions in basalt. *Comput. Geosci.* **45**, 87–97 (2011).

33. Newman, S. & Lowenstern, J. B. Volatile Calc : a silicate melt-H₂O-CO₂ solution model written in Visual Basic for excel. *Comput. Geosci.* **28**, 597–604 (2002).

34. Beermann, O. The solubility of sulfur and chlorine in H₂O-bearing dacites of Krakatau and basalts of Mt. Etna. *Etna. PhD Thesis* (2010).

35. Botcharnikov, R. E., Holtz, F. & Behrens, H. The effect of CO₂ on the solubility of H₂O-Cl fluids in andesitic melt. *Eur. J. Mineral.* **19**, 671–680 (2007).

36. Balcone-Boissard, H., Baker, D. R., Villemant, B. & Boudon, G. F and Cl diffusion in phonolitic melts: Influence of the Na/K ratio. *Chem. Geol.* **263**, 89–98 (2008).

37. Cioni, R., Bertagnini, A., Santacroce, R. & Andronico, D. Explosive activity and eruption scenarios at Somma-Vesuvius (Italy): Towards a new classification scheme. *J. Volcanol. Geotherm. Res.* **178**, 331–346 (2008).

38. Lima, A., Danyushevsky, L. V., De Vivo, B. & Fedele, L. A model for the evolution of the Mt. Somma-Vesuvius magmatic system based on fluid and melt inclusion investigations. *Dev. Volcanol. J.* **5**, 227–249 (2003).

39. Santacroce, R. et al. Age and whole rock-glass compositions of proximal pyroclastics from the major explosive eruptions of Somma-Vesuvius: A review as a tool for distal tephrostratigraphy. *J. Volcanol. Geotherm. Res.* **177**, 1–18 (2008).
40. Arrighi, S., Rosi, M. & Principe, C. Violent strombolian and subplinian eruptions at Vesuvius during post-1631 activity. *Bull. Volcanol.* **63**, 126–150 (2001).

41. Di Renzo, V. et al. Magmatic history of Somma-Vesuvius on the basis of new geochemical and isotopic data from a deep borehole (Camaldoli della Torre). *J. Petrol.* **48**, 753–784 (2007).

42. Cioni, R. et al. Compositional Layering and Syn-eruptive Mixing of a Periodically Refilled Shallow Magma Chamber: the AD 79 Plinian Eruption of Vesuvius. *J. Petrol.* **36**, 739–776 (1995).

43. Webster, J. D., De Vivo, B. & Tappen, C. Volatiles, magmatic degassing and eruptions of Mt. Somma-Vesuvius: Constraints from silicate melt inclusions, Cl and H₂O solubility experiments and modeling. *Melt Inclusions Volcan. Syst. Methods, Appl. Probl.* **5**, 207–226 vr258 (2003).

44. Alletti, M. et al. Chlorine partitioning between a basaltic melt and H₂O-CO₂ fluids at Mount Etna. *Chem. Geol.* **263**, 37–50 (2009).

45. Pappalardo, L. & Mastrolorenzo, G. Short residence times for alkaline Vesuvius magmas in a multi-depth supply system: Evidence from geochemical and textural studies. *Earth Planet. Sci. Lett.* **296**, 133–143 (2010).

46. Fulignati, P. & Marianelli, P. Tracing volatile exsolution within the 472 AD 'Pollena' magma chamber of Vesuvius (Italy) from melt inclusion investigation. *J. Volcanol. Geotherm. Res.* **161**, 289–302 (2007).

47. Marianelli, P., Sbrana, A., Metrich, N. & Cecchetti, A. The deep feeding system of Vesuvius involved in recent violent strombolian eruptions. *Geophys. Res. Lett.* **32**, 2–5 (2005).

48. Scaillet, B., Pichavant, M. & Cioni, R. Upward migration of Vesuvius magma chamber over the past 20,000 years. *Nature* **455**, 216–219 (2008).

49. Landi, P., Bertagnini, A. & Rosi, M. Chemical zoning and crystallization mechanisms in the magma chamber of the Pomici di Base plinian eruption of Somma-Vesuvius (Italy). *Contrib. to Mineral. Petrol.* **135**, 179–197 (1999).

50. Fabbrizio, A., Scaillet, B. & Carroll, M. R. Estimation of pre-eruptive magmatic water fugacity in the Phlegrean Fields, Naples, Italy. *Eur. J. Mineral.* **1**, 107–116 (2009).

51. De Natale, G., Troise, C., Trigila, R., Dolfi, D. & Chiarabba, C. Seismicity and 3-D substructure at Somma-Vesuvius volcano: Evidence for magma quenching. *Earth Planet. Sci. Lett.* **221**, 181–196 (2004).

52. Pichavant, M., Scaillet, B., Pommier, A., Iacono-Marziano, G. & Cioni, R. Nature and Evolution of Primitive Vesuvius Magmas: an Experimental Study. *J. Petrol.* **55**, 2281–2310 (2014).

53. Auger, E. Seismic Evidence of an Extended Magmatic Sill Under Mt. Vesuvius. *Science* (80–). **294**, 1510–1512 (2001).

54. Pommier, a. et al. A new petrological and geophysical investigation of the present-day plumbing system of Mount Vesuvius. *Geochemistry, Geophys. Geosystems* **11** (2010).

55. Cioni, R. Volatile content and degassing processes in the AD 79 magma chamber at Vesuvius (Italy). *Contrib. to Mineral. Petrol.* **140**, 40–54 (2000).

56. Dallai, L., Cioni, R., Boschi, C. & D'Oriano, C. Carbonate-derived CO₂ purging magma at depth: Influence on the eruptive activity of Somma-Vesuvius, Italy. *Earth Planet. Sci. Lett.* **310**, 84–95 (2011).

57. Andújar, J. & Scaillet, B. Relationships between pre-eruptive conditions and eruptive styles of phonolite–trachyte magmas. *Lithos* **152**, 122–131 (2012).

58. Balcone-Boissard, H., Villemant, B., Boudon, G. & Michel, A. Non-volatile vs volatile behaviours of halogens during the AD 79 plinian eruption of Mt. Vesuvius, Italy. *Earth Planet. Sci. Lett.* **269**, 66–79 (2008).

59. Balcone-Boissard, H. et al. Magma degassing and eruption dynamics of the Avellino pumice Plinian eruption of Somma-Vesuvius (Italy). Comparison with the Pompeii eruption. *Earth Planet. Sci. Lett.* **331–332**, 257–268 (2012).

## Acknowledgements

This research was supported by the Grant INSU-CNRS Aléas 2014. We thank M. Fialin, F. Couffignal, N. Rividi for help support at the electronic microprobe (Camparis, Paris, France) and O. Boudouma for textural analyses by SEM (Paris, France). A. Carandente and P. Belviso helped in sampling fallout deposits and S. Hidalgo carried some density measurements. D. Morgan is acknowledged for critical comments of the manuscript.

## Author Contributions

H.B.B. and G.B. defined and performed data acquisition and interpretation, wrote the manuscript and prepared the figures. R.C. participated to data interpretation and the writing of the manuscript text. J.D.W. was involved in the chlorine modelling, interpretation and the writing of the manuscript text. G.Z. participated to some data acquisition. G.O. and L.C. were involved in the interpretation and the writing of the manuscript text. All authors reviewed the manuscript.

## Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing financial interests: The authors declare no competing financial interests.

**How to cite this article:** Balcone-Boissard, H. et al. Chlorine as a geobarometer for alkaline magmas: Evidence from a systematic study of the eruptions of Mount Somma - Vesuvius. *Sci. Rep.* **6**, 21726; doi: 10.1038/srep21726 (2016).

<figure><img src="image_12.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/