# High CO2 content in magmas of the explosive andesitic Enco eruption of Mocho-Choshuenco volcano (Chile)

Jean-Guillaume Feignon, Nicolas Cluzel, Federica Schiavi, Séverine Moune, Olivier Roche, Jorge Clavero, Pierre Schiano, Marion Auxerre

► To cite this version:

Jean-Guillaume Feignon, Nicolas Cluzel, Federica Schiavi, Séverine Moune, Olivier Roche, et al. High CO2 content in magmas of the explosive andesitic Enco eruption of Mocho-Choshuenco volcano (Chile). Bulletin of Volcanology, 2022, 84 (4), pp.40. 10.1007/s00445-022-01550-y. insu-03643025

HAL Id: insu-03643025

https://insu.hal.science/insu-03643025

Submitted on 20 Apr 2022

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
# High CO₂ content in magmas of the explosive andesitic Enco eruption of Mocho-Choshuenco volcano (Chile)

Jean-Guillaume Feignon<fnref n="1" />,<fnref n="2" />✉ · Nicolas Cluzel<fnref n="1" /> · Federica Schiavi<fnref n="1" /> · Séverine Moune<fnref n="1" />,<fnref n="3" />,<fnref n="4" /> · Olivier Roche<fnref n="1" />✉ · Jorge Clavero<fnref n="5" /> · Pierre Schiano<fnref n="1" /> · Marion Auxerre<fnref n="1" />,<fnref n="6" />

Received: 7 February 2022 / Accepted: 2 March 2022 / Published online: 23 March 2022
© The Author(s) 2022

## Abstract

Mocho-Choshuenco volcano has produced several highly explosive eruptions during its history, which make it one of the most hazardous volcanoes in the southern volcanic zone of Chile, although it is still relatively little studied to date. We present a geochemical study of the products of the sub-Plinian, andesitic, Enco eruption that occurred about 1600 years ago. We determined the major and trace elements compositions, as well as the volatile (H₂O, CO₂, Cl, and S) contents of melt inclusions trapped in minerals (olivine, plagioclase, and pyroxene) using electron microprobe, ion microprobe (SIMS), and 3D confocal Raman mapping. Though the whole-rock composition of the Enco magma is andesitic (60.2 ± 1.1 wt.% SiO₂), the melt inclusions have SiO₂ contents ranging from 50.3 to 67.3 wt.%, following the magmatic series of Mocho-Choshuenco, and the compositions of the most mafic melt inclusions are close to those of the most mafic erupted magmas. Geochemical modeling indicates that mixing occurred between a mafic magma and an andesitic-to-dacitic magma. Glass analysis revealed typical parental arc magma values for H₂O (2.6–3.8 wt.%), S (116–1936 ppm), and Cl (620–1439 ppm). However, CO₂ contents are very high in some melt inclusions with concentrations above 4000 ppm (measured in the glass), suggesting trapping depths > ~17–22 km. Presence of solid carbonates inside inclusion-hosted bubbles clearly indicates that the CO₂ contents measured in the glass phase were minimum values. We conclude that a CO₂-rich basaltic magma ascended and mixed with a shallower andesitic magma. The magma cooled and exsolved high amounts of CO₂, which may have dramatically increased the pressure and triggered the highly explosive Enco eruption.

**Keywords** Carbon dioxide · Mocho-Choshuenco · Magma volatile content · Explosive eruption · 3D Raman mapping · Melt inclusions

Editorial responsibility: T.P. Fischer; Deputy Executive Editor: J. Tadeucci

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

## Introduction

Various parameters and/or properties may influence the explosivity of a volcanic system (e.g., Sparks et al. 1977; Gonnermann and Manga 2007; Scandone et al. 2007; Deegan et al. 2010; Williamson et al. 2010; Costa et al. 2013; Zimanowski et al. 2015; Cassidy et al. 2018). Differentiated magmas have high viscosities, which, combined with common high volatile contents (if magmas are not degassed), favor the development of high exsolved-gas pressure, and promote efficient magma fragmentation (Williamson et al. 2010; Costa et al. 2013). In contrast, poorly differentiated, low-viscosity magmas with (presumably) lower volatile contents commonly lead to eruptions dominated by effusive rather than explosive dynamics (Huppert and Woods 2002; Gardner 2009). However, examples of highly explosive eruptions involving basaltic to andesitic magmas (SiO₂ < ~60 wt. %) are reported in the literature, e.g., Fuego (Guatemala);

<footnote>Jean-Guillaume Feignon
jean-guillaume.feignon@univie.ac.at</footnote>

<footnote n="1">¹ Laboratoire Magmas Et Volcans, Université Clermont Auvergne, CNRS, IRD, OPGC, 63000 Clermont-Ferrand, France</footnote>

<footnote n="2">² Department of Lithospheric Research, University of Vienna, Althanstrasse 14, 1090 Vienna, Austria</footnote>

<footnote n="3">³ Institut de Physique du Globe de Paris, CNRS, Université de Paris, 75005 Paris, France</footnote>

<footnote n="4">⁴ Institut de Physique du Globe de Paris, Observatoire Volcanologique Et Sismologique de Guadeloupe, 97113 Gourbeyre, France</footnote>

<footnote n="5">⁵ Amawta Geoconsultores, Almirante Pastene 185 Of. 503, Providencia, Santiago, Chile</footnote>

<footnote n="6">⁶ Centre de Recherches Pétrographiques Et Géochimiques (CRPG), UMR CNRS 5873, Université de Lorraine, BP20, 54501 Vandœuvre les Nancy, France</footnote>
Pozzolane Rosse and Etna (Italy; Coltelli et al. 1998; Métrich et al. 2004; Spilliaert et al. 2006; Freda et al. 2010; Moretti et al. 2018); Masaya (Nicaragua; Szramek 2016); Llaima, Hudson, and Villarica (Chile; Stern et al. 2007); Tarawera (New Zealand; Houghton et al. 2004); the 1790 Keanakakoi eruption at Kilauea (Hawaii; McPhie 1990; Swanson et al. 2012); and Sunset Crater (Arizona; Allison et al. 2021). Possible causes of the high explosivity of these magmas are fast ascent and low degassing rates in the conduit (Zhang 1999), magma-water or ice interaction (Swanson et al. 2012), high magma viscosity due to high crystallinity (Lavallée et al. 2007), and unusually high volatile contents (Sable et al. 2006; Cross et al. 2014). An increasing number of studies have pointed out that highly explosive mafic magmas may also have high CO₂ (and to a lesser extent S and Cl) contents (e.g., Blundy et al. 2010; Aiuppa et al. 2011; Dallai et al. 2011; Freda et al. 2010; Allison et al. 2021). The low solubility of CO₂ (Bottinga and Javoy 1990; Le Gall and Pichavant 2016; Laporte et al. 2016), in particular, may be one of the key factors in promoting explosive dynamics of poorly differentiated magmas.

The relatively poorly studied, and dangerous Mocho-Choshuenco (Chile) represents another case of a volcano presenting explosive eruptions fed by basaltic to andesitic magmas. To better disentangle the origin of its explosive behavior and bring new data on the local magmatism, we focused our investigations on the volatile content of the products of a recent explosive eruption: the Enco-MC15 (as defined by Rawson et al. 2015) eruption. We considered this event because it generated the most mobile pyroclastic flows of the recent eruptive history of Mocho-Choshuenco. In addition, field observations do not provide evidence of phreato-magmatism (Rawson et al. 2015, and our own work), hence excluding a crucial role of magma-water or magma-ice interaction in the explosive character of this eruption. We characterized major, trace, and volatile element compositions of olivine-, pyroxene-, and plagioclase-hosted melt inclusions in order to determine the pre-eruptive conditions of the magmas that fed the Enco-MC15 eruption at Mocho-Choshuenco. Finally, we present a possible scenario of the magmatic storage conditions and the key factors controlling the explosive behavior of the eruption.

## Geological background and sampling

Mocho-Choshuenco volcano (2422 m a.s.l.) is located in Chile (39° 55'S 72° 2'W) in the Southern Volcanic Zone (SVZ) of the Andes (Moreno and Lara 2007; Fig. 1). This volcano has been characterized by 75 post-glacial (<18 ky) mainly andesitic to dacitic eruptions, including Plinian events that occurred during the Neltume (11 ka–5.3 km³ fallout deposit), the Pirehueico (9 ka–2.0 km³ fallout deposit), the Enco-MC15, and the 1864 AD eruptions (Rawson et al.

2015). The Neltume and Pirehueico eruptions probably led to a partial collapse of the edifice, forming a 4-km-wide caldera that was later filled by successive tephra falls and pyroclastic flow deposits forming the El Mocho cone (Rawson et al. 2015). The total magma volume emitted during the last 18 ky was estimated to ~20 km³, making this volcano the most productive of the SVZ (Rawson et al. 2015). It represents a major hazard as more than 400 000 people are living in a 100 km area from the volcano (Rawson et al. 2015). The Enco-MC15 sub-Plinian eruption (Rawson et al. 2015) emitted ~1.7 km³ of andesitic magma and generated pyroclastic fall and flow deposits. This event generated the most highly mobile pyroclastic flows of the recent eruptive history of Mocho-Choshuenco, with run-out distances of at least ~27 km from the vent in the eastern sector (Clavero et al. 2008).

We sampled an outcrop of a pyroclastic flow deposit of the Enco-MC15 eruption which belongs to unit UM4-5 as defined by Rawson et al. (2015). This outcrop, named VMJ32 in our study, is located on the Eastern flank of Mocho-Choshuenco volcano, 13 km away from the summit vent (Fig. 1). It is composed of several flow units whose tops are outlined by large (10–30 cm) scoria. The material is composed of coarse (mean grain size, $\phi = -3.35$ to $-2.54$) and well-sorted (sorting, $\sigma_{\phi} = 1.16$ to 1.24) juvenile pyroclasts, and rare (<5 vol. %) lithic clasts. The upper part of these flow units was sampled in the framework of this study. We found charcoal whose $^{14}$C dating gave an age of $1750 \pm 40$ BP (performed by Beta Analytic Inc., USA). After calibration using OxCal and the southern hemisphere calibration curve SH13 (Ramsey 2009), the age obtained is $1622 \pm 86$ BP. This is overlapping and consistent with calibrated ages of $1580 \pm 115$ BP given by Rawson et al. (2015) and Fontijn et al. (2016). Major element compositions of the whole rock samples were analyzed by ICP-AES at Laboratoire Magmas et Volcans, Clermont-Ferrand, France (Table 1, $\pm$10% relative error at the $2\sigma$ level). We focused our study on melt inclusion analysis, and the methods we used are presented below.

## Analytical methods

### Preparation of the melt inclusions

The lapilli-sized (main fraction between 8 and 16 mm with few larger (up to >32 mm) lapilli) pyroclasts yielding the most suitable phenocrysts to extract for investigation were selected. After crushing and sieving, the >300 µm fraction was retained to separate phenocrysts from the glass. Olivines and pyroxenes were separated by density using Bromoform (CHBr₃), whereas the remaining plagioclases and glass were separated magnetically using a FRANTZ magnetic separator.
<figure><img src="image_3.png" /><figcaption>**Fig. 1** a Map of the Southern Volcanic Zone (SVZ) with the main Holocene volcanic centers (modified from Rawson et al. 2015). b Simplified geological map of Mocho-Choshuenco volcano: units UM1 (> 350–200 ka), UM2 (< 200–130 ka), UC1 (> 100 ka), UC2 (80–20 ka), UM3 (> 60 ka), UM4-5 (18 ka to present). The outcrop of the pyroclastic flow deposit of the unit</figcaption></figure>

The selected crystal-hosted melt inclusions (hereafter denoted as MIs) were glassy, ellipsoidal, with shrinkage bubble(s) and no cracks. Inclusions aligned within fracture planes in the host minerals, due to late melt infiltration, were excluded.

Electron microprobe

Following the method of Moune et al. (2012), the major and volatile (S, Cl, and F) element compositions of 41 MIs,

Enco-MC15 considered in this study is located in unit UM4-5 and is called VMJ32 (modified from Rawson et al. 2015). c Photograph of outcrop VMJ32 showing a coarse-grained pyroclastic flow deposit of the Enco eruption between finely grained surge units. The material analyzed was sampled from juvenile material located in the upper part of the outcrop. The measuring tape is two meters long

of the host crystals, and of the groundmass glasses were analyzed using a SX-100 CAMECA electron microprobe at Laboratoire Magmas et Volcans, Clermont-Ferrand (France) with a 15 kV accelerating voltage. Analyses of glass samples were done using a probe current of 8 nA for major elements and 40 nA for sulfur, chlorine, and fluorine volatile species, with a defocused beam of 5, 10, or 20 µm according to the size of the MIs. For each sample, prior to sulfur concentration measurement, the S Kα peak maximum was
<table><caption>**Table 1** Whole rock and melt inclusion compositions including host composition. The olivine-hosted melt inclusion compositions were corrected from post entrapment crystallization, while plagioclase- and pyroxene-hosted MIs are uncorrected (i.e. measured value). Uncorrected data for olivine-hosted MIs can be found in Supplementary Material S4. Major oxides and H₂O given in wt.%, and S, Cl, and CO₂ measured in the MI glass given in ppm. Total Fe is expressed as FeO*. PEC post entrapment crystallization of the host mineral. The K_D expressed for plagioclase corresponds to the plagioclase (Albite-Anhortite)-liquid equilibrium as expressed by Putirka (2008). b.d.l. below detection limit</caption><thead><tr><th>Sample</th><th>Host composition</th><th>SiO<sub>2</sub></th><th>TiO<sub>2</sub></th><th>Al<sub>2</sub>O<sub>3</sub></th><th>FeO*</th><th>MnO</th><th>MgO</th><th>CaO</th><th>Na<sub>2</sub>O</th><th>K<sub>2</sub>O</th><th>P<sub>2</sub>O<sub>5</sub></th><th>Total</th><th>S</th><th>Cl</th><th>H<sub>2</sub>O</th><th>CO<sub>2</sub></th><th>K<sub>b</sub></th><th>PEC (%)</th></tr></thead><tbody><tr><td colspan="19" style="text-align: left; font-weight: bold;">Whole rock</td></tr><tr><td>VMJ32</td><td>60.2</td><td>1.19</td><td>16.3</td><td>7.57</td><td>0.17</td><td>2.22</td><td>5.35</td><td>4.79</td><td>1.29</td><td>0.40</td><td>99.74</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="19" style="text-align: left; font-weight: bold;">Groundmass (n=8)</td></tr><tr><td>VMJ32</td><td>61.5</td><td>1.26</td><td>15.6</td><td>6.78</td><td>0.25</td><td>2.09</td><td>5.06</td><td>4.68</td><td>1.56</td><td>0.38</td><td>99.08</td><td>138</td><td>986</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="19" style="text-align: left; font-weight: bold;">Olivine-hosted inclusions</td></tr><tr><td colspan="19" style="text-align: left; font-weight: bold;">Fo</td></tr><tr><td>O11_MII</td><td>80</td><td>51.0</td><td>1.15</td><td>19.9</td><td>6.62</td><td>0.25</td><td>3.61</td><td>10.0</td><td>3.27</td><td>0.59</td><td>0.11</td><td>96.54</td><td>1535</td><td>628</td><td>2.83</td><td>214</td><td>0.17</td><td>4.5</td></tr><tr><td>O12_MII</td><td>78</td><td>53.6</td><td>1.35</td><td>18.6</td><td>6.15</td><td>b.d.l</td><td>3.08</td><td>9.27</td><td>4.09</td><td>0.81</td><td>0.25</td><td>97.17</td><td>1610</td><td>675</td><td>2.99</td><td>541</td><td>0.16</td><td>4.0</td></tr><tr><td>O16-4_MII</td><td>81</td><td>51.5</td><td>0.99</td><td>18.9</td><td>6.84</td><td>0.15</td><td>4.41</td><td>9.62</td><td>3.22</td><td>0.58</td><td>0.08</td><td>96.21</td><td>1529</td><td>637</td><td>2.87</td><td>873</td><td>0.19</td><td>5.0</td></tr><tr><td>O16-5_MII</td><td>78</td><td>51.2</td><td>1.11</td><td>18.6</td><td>7.03</td><td>0.02</td><td>4.23</td><td>9.32</td><td>3.32</td><td>0.60</td><td>0.35</td><td>95.81</td><td>1516</td><td>684</td><td>2.62</td><td>2836</td><td>0.23</td><td>3.0</td></tr><tr><td>O16-9_MII</td><td>76</td><td>52.2</td><td>1.05</td><td>18.2</td><td>5.73</td><td>0.02</td><td>5.15</td><td>9.03</td><td>3.71</td><td>0.67</td><td>0.38</td><td>96.18</td><td>605</td><td>-</td><td>-</td><td>-</td><td>0.10</td><td>9.0</td></tr><tr><td>O16-12_MII</td><td>77</td><td>55.9</td><td>1.82</td><td>15.2</td><td>6.59</td><td>0.08</td><td>3.21</td><td>7.10</td><td>4.81</td><td>1.28</td><td>0.28</td><td>96.28</td><td>1341</td><td>1118</td><td>-</td><td>-</td><td>0.14</td><td>5.0</td></tr><tr><td>O16-13_MII</td><td>75</td><td>53.2</td><td>1.17</td><td>18.1</td><td>7.39</td><td>0.10</td><td>3.18</td><td>8.94</td><td>3.53</td><td>0.70</td><td>0.23</td><td>96.49</td><td>1936</td><td>692</td><td>0.48</td><td>225</td><td>0.23</td><td>2.5</td></tr><tr><td>O16-14_MII</td><td>81</td><td>50.6</td><td>1.11</td><td>19.1</td><td>6.80</td><td>0.15</td><td>4.28</td><td>9.69</td><td>3.22</td><td>0.60</td><td>0.28</td><td>95.87</td><td>1594</td><td>625</td><td>2.73</td><td>636</td><td>0.18</td><td>5.0</td></tr><tr><td>O16-14_MII</td><td>81</td><td>51.9</td><td>1.05</td><td>18.6</td><td>6.77</td><td>b.d.l</td><td>4.39</td><td>9.41</td><td>3.33</td><td>0.60</td><td>0.22</td><td>96.29</td><td>1494</td><td>638</td><td>-</td><td>-</td><td>0.17</td><td>5.5</td></tr><tr><td>O16-15_MII</td><td>81</td><td>52.2</td><td>1.07</td><td>18.7</td><td>6.84</td><td>0.13</td><td>4.05</td><td>9.38</td><td>3.39</td><td>0.62</td><td>0.20</td><td>96.48</td><td>1485</td><td>662</td><td>3.10</td><td>283</td><td>0.24</td><td>2.5</td></tr><tr><td>M16</td><td>76</td><td>50.3</td><td>1.00</td><td>18.0</td><td>8.13</td><td>0.12</td><td>5.05</td><td>9.71</td><td>3.01</td><td>0.51</td><td>0.21</td><td>96.01</td><td>1220</td><td>620</td><td>-</td><td>-</td><td>0.14</td><td>8.0</td></tr><tr><td colspan="19" style="text-align: left; font-weight: bold;">Plagioclase-hosted inclusions</td></tr><tr><td colspan="19" style="text-align: left; font-weight: bold;">An</td></tr><tr><td>P11_MII</td><td>50</td><td>62.3</td><td>0.83</td><td>14.0</td><td>4.66</td><td>0.43</td><td>1.37</td><td>3.62</td><td>3.64</td><td>1.77</td><td>0.35</td><td>92.95</td><td>200</td><td>1223</td><td>2.85</td><td>4213</td><td>0.14</td><td>-</td></tr><tr><td>P12_MII</td><td>76</td><td>56.3</td><td>1.49</td><td>14.3</td><td>8.09</td><td>0.26</td><td>3.39</td><td>5.33</td><td>4.03</td><td>1.46</td><td>0.33</td><td>94.99</td><td>1184</td><td>1045</td><td>3.84</td><td>1075</td><td>0.07</td><td>-</td></tr><tr><td>P13_MII</td><td>50</td><td>64.8</td><td>0.80</td><td>15.7</td><td>4.85</td><td>0.16</td><td>1.01</td><td>3.35</td><td>4.81</td><td>1.83</td><td>0.26</td><td>97.54</td><td>181</td><td>1223</td><td>2.78</td><td>212</td><td>0.11</td><td>-</td></tr><tr><td>P14_MII</td><td>51</td><td>65.9</td><td>0.78</td><td>15.4</td><td>4.39</td><td>0.03</td><td>1.15</td><td>3.27</td><td>5.04</td><td>1.80</td><td>0.31</td><td>98.00</td><td>120</td><td>1142</td><td>1.21</td><td>473</td><td>0.09</td><td>-</td></tr><tr><td>P15_MII</td><td>50</td><td>65.5</td><td>0.63</td><td>15.6</td><td>3.90</td><td>0.16</td><td>0.92</td><td>2.95</td><td>5.33</td><td>1.91</td><td>0.24</td><td>97.10</td><td>129</td><td>1019</td><td>-</td><td>-</td><td>0.08</td><td>-</td></tr><tr><td>P16_MII</td><td>48</td><td>64.4</td><td>0.71</td><td>15.9</td><td>4.12</td><td>0.18</td><td>0.99</td><td>3.26</td><td>4.93</td><td>1.79</td><td>0.22</td><td>96.51</td><td>114</td><td>1075</td><td>2.84</td><td>2710</td><td>0.11</td><td>-</td></tr><tr><td>P17_MII</td><td>50</td><td>60.8</td><td>0.97</td><td>15.5</td><td>5.93</td><td>0.17</td><td>2.11</td><td>4.36</td><td>4.95</td><td>1.64</td><td>0.37</td><td>96.79</td><td>273</td><td>969</td><td>3.16</td><td>321</td><td>0.14</td><td>-</td></tr><tr><td>P17_MII</td><td>50</td><td>64.4</td><td>0.77</td><td>15.9</td><td>4.62</td><td>0.11</td><td>1.21</td><td>3.33</td><td>4.67</td><td>1.82</td><td>0.07</td><td>96.96</td><td>212</td><td>1240</td><td>2.66</td><td>860</td><td>0.11</td><td>-</td></tr><tr><td>P18_MII</td><td>49</td><td>64.1</td><td>0.72</td><td>15.8</td><td>4.57</td><td>0.07</td><td>1.07</td><td>3.15</td><td>5.26</td><td>1.86</td><td>0.28</td><td>96.90</td><td>170</td><td>1124</td><td>2.48</td><td>4430</td><td>0.10</td><td>-</td></tr><tr><td>P19_MII</td><td>47</td><td>62.8</td><td>0.98</td><td>16.0</td><td>5.67</td><td>0.16</td><td>1.26</td><td>3.49</td><td>3.93</td><td>1.85</td><td>0.31</td><td>96.43</td><td>243</td><td>1309</td><td>1.85</td><td>689</td><td>0.16</td><td>-</td></tr><tr><td>M2</td><td>49</td><td>63.1</td><td>0.82</td><td>15.5</td><td>4.94</td><td>0.16</td><td>1.11</td><td>3.11</td><td>5.03</td><td>2.05</td><td>0.46</td><td>96.29</td><td>172</td><td>1373</td><td>-</td><td>-</td><td>0.10</td><td>-</td></tr><tr><td>M4</td><td>49</td><td>60.7</td><td>1.27</td><td>14.8</td><td>6.84</td><td>0.16</td><td>1.99</td><td>4.28</td><td>4.79</td><td>1.74</td><td>0.61</td><td>97.17</td><td>434</td><td>1174</td><td>-</td><td>-</td><td>0.11</td><td>-</td></tr><tr><td>M6</td><td>48</td><td>65.8</td><td>0.62</td><td>15.2</td><td>3.30</td><td>0.10</td><td>0.95</td><td>2.96</td><td>4.93</td><td>1.60</td><td>0.12</td><td>95.61</td><td>115</td><td>773</td><td>-</td><td>-</td><td>0.07</td><td>-</td></tr><tr><td>M8</td><td>44</td><td>66.0</td><td>0.69</td><td>15.2</td><td>4.42</td><td>0.18</td><td>0.97</td><td>2.95</td><td>5.10</td><td>1.72</td><td>0.12</td><td>97.33</td><td>117</td><td>1099</td><td>-</td><td>-</td><td>0.10</td><td>-</td></tr><tr><td>M15</td><td>45</td><td>64.4</td><td>0.58</td><td>14.7</td><td>4.74</td><td>0.09</td><td>1.13</td><td>3.08</td><td>5.06</td><td>1.73</td><td>0.35</td><td>95.85</td><td>160</td><td>1262</td><td>-</td><td>-</td><td>0.09</td><td>-</td></tr></tbody></table>

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 1 (continued)**</caption><thead><tr><th>Sample</th><th>Host com- position</th><th>SiO<sub>2</sub></th><th>TiO<sub>2</sub></th><th>Al<sub>2</sub>O<sub>3</sub></th><th>FeO*</th><th>MnO</th><th>MgO</th><th>CaO</th><th>Na<sub>2</sub>O</th><th>K<sub>2</sub>O</th><th>P<sub>2</sub>O<sub>5</sub></th><th>Total</th><th>S</th><th>Cl</th><th>H<sub>2</sub>O</th><th>CO<sub>2</sub></th><th>K<sub>b</sub></th><th>PEC (%)</th></tr></thead><tbody><tr><td>M20a</td><td>50</td><td>62.5</td><td>0.95</td><td>15.5</td><td>5.01</td><td>0.20</td><td>1.19</td><td>3.21</td><td>4.94</td><td>1.78</td><td>0.38</td><td>95.66</td><td>241</td><td>1364</td><td>-</td><td>-</td><td>0.12</td><td></td><td></td></tr><tr><td>M24</td><td>49</td><td>64.1</td><td>0.72</td><td>15.5</td><td>4.30</td><td>0.19</td><td>1.18</td><td>3.45</td><td>5.02</td><td>1.87</td><td>0.09</td><td>96.41</td><td>169</td><td>1310</td><td>-</td><td>-</td><td>0.12</td><td></td><td></td></tr><tr><td>M1</td><td>51</td><td>62.9</td><td>0.52</td><td>15.2</td><td>4.46</td><td>0.20</td><td>1.18</td><td>3.29</td><td>4.77</td><td>1.93</td><td>0.20</td><td>94.65</td><td>154</td><td>1221</td><td>-</td><td>-</td><td>0.09</td><td></td><td></td></tr><tr><td>M3</td><td>49</td><td>62.6</td><td>0.74</td><td>15.2</td><td>4.12</td><td>0.07</td><td>1.12</td><td>3.11</td><td>5.25</td><td>2.08</td><td>0.35</td><td>94.72</td><td>189</td><td>1393</td><td>-</td><td>-</td><td>0.09</td><td></td><td></td></tr><tr><td>M11</td><td>50</td><td>66.2</td><td>0.62</td><td>15.0</td><td>4.00</td><td>0.17</td><td>0.91</td><td>2.96</td><td>5.25</td><td>1.79</td><td>0.18</td><td>97.08</td><td>147</td><td>999</td><td>-</td><td>-</td><td>0.08</td><td></td><td></td></tr><tr><td>M12</td><td>50</td><td>63.7</td><td>0.79</td><td>15.1</td><td>4.41</td><td>0.16</td><td>1.05</td><td>3.11</td><td>5.08</td><td>1.95</td><td>0.26</td><td>95.60</td><td>174</td><td>1188</td><td>-</td><td>-</td><td>0.09</td><td></td><td></td></tr><tr><td>M14</td><td>48</td><td>63.8</td><td>0.40</td><td>15.5</td><td>4.08</td><td>0.21</td><td>1.01</td><td>3.08</td><td>5.45</td><td>1.65</td><td>0.12</td><td>95.25</td><td>141</td><td>1123</td><td>-</td><td>-</td><td>0.09</td><td></td><td></td></tr><tr><td>M25</td><td>49</td><td>59.3</td><td>1.22</td><td>14.8</td><td>6.52</td><td>0.24</td><td>2.34</td><td>4.55</td><td>4.43</td><td>1.47</td><td>0.26</td><td>95.10</td><td>242</td><td>1025</td><td>-</td><td>-</td><td>0.17</td><td></td><td></td></tr><tr><td colspan="20" style="text-align: left; font-weight: bold;">Pyroxene-hosted inclusions</td></tr><tr><td colspan="20" style="text-align: left; font-weight: bold;">En</td></tr><tr><td>M10</td><td>67</td><td>67.3</td><td>0.75</td><td>14.2</td><td>4.14</td><td>0.15</td><td>0.78</td><td>2.56</td><td>4.21</td><td>1.91</td><td>0.18</td><td>96.09</td><td>116</td><td>1182</td><td>-</td><td>-</td><td>0.14</td><td></td><td></td></tr><tr><td>M20b</td><td>59</td><td>64.3</td><td>1.22</td><td>14.0</td><td>4.75</td><td>0.21</td><td>0.51</td><td>2.53</td><td>3.69</td><td>1.87</td><td>0.35</td><td>93.44</td><td>156</td><td>1439</td><td>-</td><td>-</td><td>0.12</td><td></td><td></td></tr><tr><td>M7a</td><td>42</td><td>58.6</td><td>1.30</td><td>15.7</td><td>7.61</td><td>0.24</td><td>1.99</td><td>5.20</td><td>4.77</td><td>1.32</td><td>0.78</td><td>97.54</td><td>712</td><td>1058</td><td>-</td><td>-</td><td>0.20</td><td></td><td></td></tr><tr><td>M7b</td><td>42</td><td>61.5</td><td>1.51</td><td>15.7</td><td>5.34</td><td>0.14</td><td>1.02</td><td>3.21</td><td>5.43</td><td>1.90</td><td>0.53</td><td>96.22</td><td>177</td><td>1308</td><td>-</td><td>-</td><td>0.24</td><td></td><td></td></tr><tr><td>M7c</td><td>42</td><td>61.3</td><td>0.94</td><td>15.8</td><td>5.64</td><td>0.10</td><td>1.28</td><td>3.37</td><td>4.99</td><td>1.55</td><td>0.58</td><td>95.54</td><td>278</td><td>1230</td><td>-</td><td>-</td><td>0.25</td><td></td><td></td></tr><tr><td>M7d</td><td>42</td><td>63.1</td><td>1.16</td><td>16.0</td><td>5.38</td><td>0.23</td><td>1.55</td><td>3.89</td><td>4.53</td><td>1.43</td><td>0.41</td><td>97.66</td><td>303</td><td>1333</td><td>-</td><td>-</td><td>0.22</td><td></td><td></td></tr><tr><td>M7e</td><td>42</td><td>61.5</td><td>1.33</td><td>15.3</td><td>6.24</td><td>0.30</td><td>1.59</td><td>3.82</td><td>4.48</td><td>1.47</td><td>0.23</td><td>96.26</td><td>327</td><td>1151</td><td>-</td><td>-</td><td>0.20</td><td></td><td></td></tr></tbody></table>
first located by using the automated peak-scan routine, and the value of the peak position was recorded in the analysis setup (as described in Rose-Koga et al. 2020). Variations in the wavelength of sulfur Kα X-ray as a function of its oxidation state in silicate glasses (Wallace and Carmichael 1994; Métrich and Clocchiatti 1996) were taken into account during S analysis. For minerals analyses, the beam current was 15 nA, and the beam size was 1 µm. Spots were both acquired at the edge of the inclusions and at the mineral rim. The precision of the electron microprobe (2σ) was less than 10% for major elements, except for MnO, TiO₂, and P₂O₅, and it was 6% for S and Cl. Fluorine data were excluded as they were all below the detection limit (<370 ppm).

## Mass spectrometry (SIMS and LA-ICPMS)

MIs larger than 40 µm in diameter were selected for mass spectrometry analysis. H₂O and CO₂ contents of 16 MIs (7 olivine- and 9 plagioclase-hosted) were analyzed using the CAMECA IMS 1270 Ion Probe (SIMS) at the Centre de Recherches Pétrographiques et Géochimiques, Nancy (France), with a 15-µm beam size for all analyses. Pyroxene-hosted MIs were not selected for SIMS and LA-ICPMS analysis being their size lower than the beam size. The selected crystals were transferred directly to indium sample holders, and no electron microprobe analyses were performed before using the SIMS in order to avoid carbon-coating contamination. Samples were gold coated and held in a gas chamber for degassing and cleaning (“pre-sputtering”). The glasses were sputtered with a 10 kV Cs⁺ primary beam of 10–15 nA. When a high CO₂ (> 1500 ppm) content was detected, the analysis was performed three times to check for any possible carbon contamination. Before and after each analytical session, a series of international reference standard basaltic and andesitic glasses were used for calibration of water and carbon contents (measurements on standard material and calibration curves are available in Supplementary Material S1). The error (2σ) was 1.5% for H₂O and 15% for CO₂.

In situ trace element analyses of 10 MIs (4 olivine- and 6 plagioclase-hosted) were performed by laser ablation inductively coupled-plasma mass-spectrometry (LA-ICPMS) at Laboratoire Magmas et Volcans using a Resonetics M50 Excimer Laser with a 193-nm wavelength coupled to an ICPMS Agilent 7500cs. Raw data were processed with Glitter Software using CaO content as an internal standard (see method details and analytical errors in Venugopal et al. 2016).

## Raman spectrometry

Raman measurements were carried out at the Laboratoire Magmas et Volcans using a Renishaw inVia confocal microspectrometer equipped with a 532-nm diode laser, a

Peltier-cooled CCD detector, and a Leica DM 2500 M optical microscope with a motorized XYZ stage. Scattered light was collected in a back-scattered geometry. Laser power on the sample was reduced to 8–16 mW to avoid damaging the samples. A 2400 lines/mm grating, a 20-µm slit aperture (high confocality setting), and a 100× microscope objective were used for the analyses. These analytical conditions resulted in spectral resolution of ~0.5 cm⁻¹ and lateral spatial resolution of ~1 µm. Vertical resolution was ~2–3 µm near the surface but deteriorated with depth mainly due to light refraction at the air/glass/bubble interfaces. Daily calibration of the spectrometer was performed based on the 520.5 cm⁻¹ peak of Si.

Two- or three-dimensional Raman mapping was performed on four MI-hosted bubbles. The 2D maps were acquired in the plane of the maximum bubble diameter with a step size of 1 µm. For 3D maps, we first performed depth profiles to define the vertical dimension of the map. The step size between acquisition points was 1 or 2 µm on the x- and y-axis and 2 µm on the z-axis. The acquisition time ranged from 20 to 60 s/point, and the spectra were continuously collected over 3 days in the ~200–1400 cm⁻¹ wavenumber range. Each 3D mapping produced a single hypermap of several thousand of spectra that was processed using WiRE™ 4.2 software. We applied a polynomial baseline correction to the entire dataset. For each phase present in the bubble, we selected a characteristic peak, whose intensity variations in the analyzed volume informed about the spatial distribution of that phase in the bubble. The intense peak located between 1085 and 1097 cm⁻¹ (Frezzotti et al. 2012) was used to map solid carbonates (Mg-calcite, magnesite, dolomite, and siderite), while a few grains of sulfides showing peaks near 320, 340, and 370 cm⁻¹ were identified as either pyrrhotite or pyrite. The Fermi doublet peaks of fluid CO₂ at ~1286 and 1388 cm⁻¹ were barely recognizable in some spectra as their intensities remained very weak, most likely due to low fluid density. The calibration by Lamadrid et al. (2017) was used to determine CO₂ density.

Several 2D sections along the z-axis were processed using ImageJ software to obtain the area fractions of each phase present inside the bubble. To estimate the volumetric fraction of each phase, the area fraction in each section was multiplied by the vol% of the associated spherical segment (i.e., the volume between adjacent x-y sections), and then the volumetric fractions of each phase in all spherical segments were summed (for further details, the reader is referred to Schiavi et al. 2020). The mass of carbonates was estimated from the measured volume and assuming a mean carbonate density value of 3.2 g/cm³. The total CO₂ concentration in a MI was given by the sum of the CO₂ dissolved in the melt and the CO₂ contained in the bubble. The inclusion volume was calculated assuming that inclusion thickness is equal to the arithmetic mean of the two
inclusion axes visible in the polished plane, as it seems to be the most accurate method (Tucker et al. 2019). The errors associated with the estimated inclusion volume and CO₂ concentration are relatively high because they include the uncertainties related to the spatial resolution of the technique, the decreasing intensity of the Raman signal with depth, the approximation of the applied geometric model, and the assumption about phase densities (more details on the method can be found in Supplementary Material S2 and in Schiavi et al. 2020).

# Results

## Petrography and textures of the melt inclusions

Samples contain plagioclase (mainly), clinopyroxene, olivine, rare orthopyroxene phenocrysts, and bubbles (see Supplementary Material S3). The phenocrysts represent ~5–10 vol.% of the whole rock. Finely crystallized plagioclase xenoliths are rare. Amphibole is absent, suggesting no excess of water. The microlites represent ~15–25 vol.% of the groundmass, whereas the vesicularity is ~30–40 vol.%. Figure 2 shows some of the textures commonly observed in MIs analyzed in this study. Three MI textures are identified in all studied host phases: (1) the typical rounded glassy case with one shrinkage bubble occupying less than 10% of the MI volume (Fig. 2a), (2) a shrinkage bubble occupying more than 50% of the inclusion volume (Fig. 2b), and (3) tens to hundreds of small bubbles ~1 to 10 µm in diameter, distributed along the rim or inside the whole MI volume, and occupying 10–50 vol. % of the MI, which are generally associated with a larger shrinkage bubble (Fig. 2c). The two last aforementioned textures represent 55% of the analyzed MIs.

## Compositions of melt inclusions

### Post-entrapment modifications

Following entrapment, MIs can undergo modifications that include bubble formation, host crystallization, chemical diffusion (e.g., Fe and Mg diffusion in olivine- and pyroxene-hosted MIs), formation of daughter minerals, and volatile leakage (Wallace et al. 2015; Le Voyer et al. 2017). All analyzed MIs lacked daughter minerals and were not fractured. Crystallization of the host mineral, however, commonly leads to the formation of an overgrowth rim which may also modify the MI composition (see also, Schiano 2003). The presence of an overgrowth at the MI rim is assumed based on the disequilibrium between the melt and host-crystal given by $K_D$ values. This rim must be accounted for in order to determine the original melt composition. We used the well-established Fe–Mg distribution coefficient between olivine and melt to correct MIs compositions ($K_{Dol-melt} = [X^{Fe}_{ol}/X^{Mg}_{ol}] / [X^{Fe}_{melt}/X^{Mg}_{melt}] = 0.30 \pm 0.03$; Roeder and Emslie 1970, Toplis 2005). The $Fe^{3+}/Fe_{tot}$ ratio was calculated to be ~0.2 based on the $fO_2$ (from NNO-0.5 to NNO + 0.5, Rawson et al. 2015). For the Enco-MC15 Mocho-Choshuenco eruption, the uncorrected $K_{Dol-melt}$ values obtained from measurements range between 0.10 and 0.24, indicative of post-entrapment crystallization (PEC). The oxide contents of the MIs were corrected by adding olivine (i.e., with the starting composition being the measured composition of olivine to which the MI belongs) incrementally to the inclusion composition until the equilibrium was reached, assuming that these corrected values represent the equilibrium compositions (PEC ranging from 2.5 to 9%). Uncorrected and corrected data for major, trace, and volatile elements composition, $K_D$ and PEC values are reported

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Textures of the MIs. The identification of inclusions refers to Table 1. **a** Inclusion Ol6-15_MI1 with one retraction bubble of small volume. **b** Inclusion Ol6-13_MI1 with a big bubble (~40% of the inclusion volume).</figcaption></figure>

**c** Inclusions Ol6-5_MI1 with a shrinkage bubble and numerous small bubbles inside the inclusion
in Table 1 and in Supplementary material S4. Taking into account the PEC, the corrected MI compositions hosted in olivines represent the melt trapped at the time of host crystallization.

Correction of pyroxene- and plagioclase-hosted inclusions from host crystallization is more challenging due to the uncertainty associated with their respective distribution coefficients. The Fe-Mg distribution coefficient between clinopyroxene and MIs measured in the studied samples is $K_D = 0.20–0.25$, which is slightly lower, but within the uncertainty of the equilibrium $K_D$ of $0.27 \pm 0.08$ provided in literature (Putirka 2008). This suggests that host crystallization of pyroxene within the MIs was minimal. The measured $K_D$ between orthopyroxene and MIs is equal to $0.12–0.14$. Unlike the clinopyroxene-hosted MIs, the calculated $K_D$ using Putirka (2008) equation ($K_D^{Fe–Mg}(opx-liq) = 0.4805–0.3733X_{Si}^{liq}$, where $X_{Si}$ is the fraction of silica cation in the melt) is ~0.30, indicating that the orthopyroxene-hosted MIs are not at equilibrium. Consequently, these data were not used for modeling. The $K_D$ for the plagioclase (albite-anorthite)-liquid equilibrium found in literature is $0.10 \pm 0.05$ for a magmatic liquid with a temperature less than 1050 °C (Putirka 2008). In our samples, the temperature of the plagioclase-hosted MIs (andesitic composition) is estimated at $1035 \pm 60$ °C (see the “Petrography and textures of the melt inclusions” section), and the observed $K_D$ is between 0.07 and 0.17, suggesting that the melt may be at equilibrium with the host mineral. Taken as a whole, the major element compositions of most of the pyroxene- and plagioclase-hosted MIs follow the trend defined by the bulk-rock samples (Fig. 3), suggesting negligible modification of the original melt composition by post-entrapment crystallization. As no correction was applied to the MIs hosted in pyroxene and plagioclase, the concentrations of volatile and incompatible elements in these inclusions could be slightly higher than the original concentrations in the trapped melt. In this regard, we point out that even a large amount of PEC (i.e., 10%) has a very low impact on concentration of incompatible elements.

## Major and trace element compositions

The whole rock composition is andesitic, with 60.2 wt.% SiO₂, 2.22 wt.% MgO, and 6.08 wt.% total alkalis (Table 1). The groundmass glass that represents the last stage in the melt evolution is andesitic and homogeneous, with 61.5 wt.% SiO₂ and 6.24 wt.% total alkalis representing the average of eight measurements. The MIs, however, have compositions ranging from basaltic to dacitic depending on the nature of the host mineral. The host olivines have a narrow forsterite range, $F_{O_{75–81}}$ (Table 1). Melt inclusions from the same olivine host do not vary considerably in major element contents. The olivine-hosted MIs have well constrained oxide

contents and are basaltic to basaltic-andesites in composition with: 50.3–55.9 wt.% SiO₂, 3.08–5.15 wt.% MgO, and 0.51–1.28 wt.% K₂O (Table 1, Fig. 3). These compositions are more primitive than the whole rock and groundmass glasses. The pyroxene (En₄₂–₆₇) and plagioclase (An₄₄–₅₁ and one An₇₆) hosted MIs are more evolved, being andesitic to dacitic in composition, with 56.3–67.3 wt.% SiO₂, 0.78–3.39 wt.% MgO, and 1.43–2.08 wt.% K₂O. Moreover, these compositions are generally more evolved than those of the groundmass and whole rock. These data can be compared with the eruptive glass compositions at Mocho-Choshuenco that range from ca. 55 wt.% to ca. 73 wt.% SiO₂ (Rawson et al. 2015) and the erupted whole rock compositions that range from ca. 50 wt.% to ca. 68 wt.% SiO₂ (Rawson et al. 2016). The most mafic MIs compositions, in particular, are close to that of the most mafic erupted magmas, whereas the most evolved MIs are following the general trend of Mocho-Choshuenco erupted magmas.

Trace element compositions in the inclusions show typical arc-like patterns with low concentrations of HFSEs such as Nb, and enrichment in U, Pb, and Ba, as shown in the spider diagram (Supplementary Material S4). The compositional variability of the observed trends is related to the composition of the host minerals, so that the plagioclase-hosted MIs are more enriched in trace elements, whereas more primitive MIs in olivines display lower trace element concentrations.

## Volatile element compositions

The different volatile species show contrasting patterns as a function of the degree of differentiation (i.e., increasing K₂O content, Fig. 4). The H₂O contents in olivine-hosted inclusions range from 2.62 to 3.10 wt.%, except for an inclusion with a large bubble (Ol6-13_MI1, 0.48 wt.%). The plagioclase-hosted MIs present a large variability in H₂O content, ranging from 1.21 up to 3.84 wt.%. The highest H₂O concentration is found in a MI with hundreds of small bubbles hosted in An₇₆ plagioclase (Pl2_MI1). The plagioclase-hosted MIs show a general trend of decreasing H₂O content with increasing K₂O (Fig. 4c), which is consistent with crystallization accompanied by degassing of a H₂O-saturated magma.

The CO₂ content measured in the MI glass is below 1100 ppm in most cases, independent of the host mineral, but four inclusions show much higher values of 2836 ppm in an olivine-hosted MI (Ol6-5_MI1), and 2710 ppm (Pl6_MI1), 4213 ppm (Pl1_MI1), and 4430 ppm (PIM1_MI1) in plagioclase-hosted MIs. Two of the inclusions mentioned above (Ol6-5_MI1 and PIM1_MI1) have several bubbles inside the inclusion (hundreds of them in Ol6-5_MI1, see Fig. 2c). The two other inclusions with high CO₂ content have a regular shape with a shrinkage
<figure><img src="image_6.png" /><figcaption>**Fig. 3** a–f Major oxide diagrams of investigated MIs. Olivine-hosted inclusions represent the primitive end-member and the plagioclase-hosted inclusions the most evolved one. The olivine-hosted MIs compositions are corrected from post-entrapment crystallization of host mineral at the rim of the MI. The two curves represent Petrolog3</figcaption></figure>

output of fractional crystallization models. The numbers represent the percentage of fractional crystallization. Diamonds on the fractional crystallization (FC) model curves represent increments with steps of: 1% from 1 to 5% FC, 5% from 5 to 30% FC, and 10% from 30% to the end of the FC. Error bars are 2σ. See Table 1 for values
bubble similar to the MI shown in Fig. 2a. Furthermore, it is worth noting that the two other MIs containing more than a hundred of small bubbles (Ol6-4_MI1 and Pl2-MI1) have relatively high CO₂ contents in the glass: 873 ppm and 1075 ppm, respectively. Raman analysis of olivine-hosted inclusions containing a shrinkage bubble revealed precipitation of solid carbonates at the rim of the bubble. In the analyzed bubbles (Ol6-4_MI1, Ol6-5_MI1, Ol6-14_MI1, and Ol6-15_MI1 MIs), the volume of the solid carbonates ranges from 5 ± 2 to 26 ± 7% of the bubble volume (supplementary material S2). Consequently, the quantity of CO₂ within the bubble represents 41 to 93% of the total amount, excepted for Ol6-5 with only 10 to 20% of the CO₂ inside the bubble (see supplementary material S2). Carbonates were observed only in the bubbles of olivine-hosted MIs (Fig. 5). When CO₂ stored in carbonates is taken into account, the corrected values of CO₂ for Ol6-4_MI1, Ol6-5_MI1, Ol6-14_MI1, and Ol6-15_MI1 (Fig. 5) are 2033 ± 187, 3155 ± 190, 1348 ± 266, and 3077 ± 744 ppm, respectively (873, 2836, 636, and 283 ppm in the MI glasses, respectively). The presence of carbonates indicates that the CO₂ contents in MIs (glass)

given in Table 1 are minima, underestimated values for olivine-hosted MIs. Our results are in agreement with previous estimates of CO₂ lost in the shrinkage bubble during post-entrapment (15–98%; Anderson and Brown 1993; Moore et al. 2015; Wallace et al. 2015; Robidoux et al. 2016; Tucker et al. 2019; Venugopal et al. 2020; Schiavi et al. 2020; Allison et al. 2021). These results suggest that the magma was originally very rich in CO₂, with contents possibly larger than 4000 ppm.

The S content decreases with increasing degree of differentiation (i.e., increasing K₂O; Fig. 4b). It ranges from 605 to 1936 ppm in olivine-hosted MIs. It is in general lower in pyroxene- and plagioclase-hosted MIs (excepted two inclusions with 712 and 1184 ppm). Considering oxidation conditions at Mocho-Choshuenco between NNO-0.5 and NNO+0.5 (Rawson et al. 2015), the S content at saturation in basaltic melts is estimated at 1500 ± 500 ppm at 200 MPa (Wallace and Carmichael 1992; Lesne et al. 2015) and is expected to vary from 1500 to 4500 ppm as fO₂ increases at ~1 GPa (Jugo 2009). The large variability in the S content of the olivine-hosted MIs (Fig. 4b) suggests saturation of the basaltic magma. In contrast, the Cl content generally

<figure><img src="image_7.png" /><figcaption>**Fig. 4** a–d Volatile contents relative to K₂O contents of the studied MIs. For four olivine-hosted MIs, both corrected (bubble + glass) and uncorrected (only glass) CO₂ contents are reported. The variability of S, H₂O, and CO₂ observed for a given K₂O value is an evidence of</figcaption></figure>

degassing. The high CO₂ contents in both basaltic and andesitic melts suggest that the two magmas that mixed and fed the Enco eruption were rich in CO₂

<figure><img src="image_8.png" /></figure>
**Fig. 5** **a** Raman spectrum showing the occurrence of solid carbonates in the bubble of Ol6-14_MI1 MI (red spectrum). The olivine-host mineral is represented by the orange spectrum. The numbers in bold and the table indicate the characteristic peaks of calcite and Mg-calcite. A clear, intense peak is centered at 1085 cm⁻¹. A 3D-mapping was performed in the area with small bubbles, but CO₂ peaks could not be detected mainly due to the small size of the bubbles (< 1–2 µm in diameter) and probably low CO₂ density. **b** Typical Raman spectrum from the large bubble of Ol6-15_MI1 MI.

As for Ol6-14_MI1, the curve is mainly characterized by an intense and large peak centered at 1092 cm⁻¹, which indicates the occurrence of different solid carbonates (probably magnesite, siderite, and/or dolomite). The other characteristic peaks for these solid carbonates are indicated in the table.
**c** Raman 3D-map of the bubble in the Ol6-15 MI, showing the distribution of carbonates inside the bubble. The blue scale represents the number of counts: a high number of counts indicate a high content in solid carbonates (bright blue), whereas dark areas correspond to low carbonates content. In addition, solid carbonates and bubble volumes (V) are indicated.

The yellow area corresponds to sulfur; S in the bubble is not discussed here; however, it confirms that the bubble is also a reservoir for S, as suggested in Venugopal et al. (2020). A photograph of the bubble is shown below the map.

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

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

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

**Fig. 6** **a–d** Incompatible and compatible trace elements, with **a** Th/Ce versus Th, **b** Rb/Co versus Rb, **c** Sr/Y versus Sr, and **d** Zr/V versus Zr. The blue curve represents output of mixing calculation model (Günduz and Asan 2021) between the most primitive and the most evolved measured MIs. The diamonds represent increment of 10% mixing by addition of the evolved end-member. The modeling fits very well with the measurements, supporting magma mixing as a main process in melt evolution. In contrast, modeled fractional crystallization (FC) paths (black dashed line) obtained using Petrolog3

(Danyushevsky and Plechov 2011) do not fit with the entire range of measured compositions, excepted in the case of Sr/Y versus Sr content. The black diamonds represent increments with steps of: 1% from 1 to 5% FC, 5% from 5 to 30% FC, and 10% from 30% to the end of the fractional crystallization. Distribution coefficients of the selected trace elements were obtained from LaTourrette and Burnett (1992), Hart and Dunn (1993), Forsythe et al. (1994), Bindeman et al. (1998), and Zanetti et al. (2004)

increases with differentiation (Fig. 4a), which is consistent with its incompatible behavior during differentiation. It ranges from 620 to 1118 ppm in olivine-hosted MIs and from 773 to 1439 ppm in pyroxene- and plagioclase-hosted MIs. Considering the Cl solubility in a hydrated magma (less than 1 wt.% Cl; Balcone-Boissard et al. 2016), the basaltic and andesitic melts of the Mocho-Choshuenco Enco eruption were undersaturated in Cl.

# Discussion

## Magmatic processes and storage conditions

The MI compositions permitted us to estimate their trapping temperature and pressure. In order to estimate MIs trapping temperatures, olivine-liquid thermometry was used for basaltic MIs. Equation 22 from Putirka (2008), considering
<figure><img src="image_13.png" /><figcaption>**Fig. 7** CO₂ vs. H₂O contents diagram for the studied MIs. The green circles represent the corrected CO₂ content for the basaltic MIs Ol6-4_MI1, Ol6-5_MI1, Ol6-14_MI1, and Ol6-15_MI1. Closed-system (with and without vapor in excess) and open-system degassing paths starting from the composition of the CO₂- and H₂O-richest basaltic MI have been calculated using VolatileCalc (Newman and Lowenstern 2002). Most of data are well explained by degassing under closed-system conditions. In addition, the high CO₂ content measured in two plagioclase-hosted MIs suggests the occurrence of CO₂ flushing from depth. One sample (Pl2_MI1, 3.84 wt. % H₂O) is well outside this path and could be an indicator of a volatile-supersaturated melt (Pichavant et al. 2013)</figcaption></figure>

hydrous conditions, yields temperatures of ~1050–1150 °C (note that the standard error of estimate (SEE) is 29 °C). It is well-known that water loss can occur due to H⁺ diffusion through the host olivine (Chen et al. 2011). In the set of studied MIs, there is no clear evidence of post-entrapment water loss (except for Ol6-13_MI1); however, we cannot exclude its occurrence. In this regard, lapilli smaller than 29 mm in size are expected to have similar cooling behavior to ash, with limited water loss (Lloyd et al. 2013). Thus, for the olivine-liquid thermometry, the calculation was based on all MIs except Ol6-13_MI1. For andesitic and dacitic MIs, plagioclase-liquid thermometry, using Eq. 24a from Putirka (2008), yields temperatures of ~975–1095 °C (with SEE of 36 °C). Additionally, clinopyroxene-liquid thermometry, using the pressure-independent parametrization of Putirka et al. (1996), results in temperatures of ~1045–1105 °C (with SEE of 69 °C).

We estimated trapping pressure of the MIs using the models of Newman and Lowenstern (2002), Papale et al. (2006) and Witham et al. (2012) for basaltic MIs, and the experimental study of Botcharnikov et al. (2006) for andesitic melts. When considering this approach, loss of volatiles (H₂O, CO₂, S) from the melt into the bubble and/or to the host crystal must be taken into account. Diffusion of volatiles from the melt to the bubble can be important, as shown by Raman investigations (e.g., Esposito et al. 2016; Robidoux et al. 2018; Venugopal et al. 2020; Schiavi et al. 2020); therefore, the amounts measured in the glass are minimal

estimates. We calculated that a basaltic liquid containing 0.3 wt.% CO₂ and 3.0 wt.% H₂O is in equilibrium with a fluid phase containing a mole fraction of H₂O ≤ 0.1 at a temperature of 1100 °C and a pressure of 600 MPa. The very low amount of H₂O that is predicted to be present inside the bubble would most likely be in the form of a very thin fluid film on the bubble wall (Venugopal et al. 2020), which may not be detected during the analyses, in agreement with our observations. Moreover, significant amounts of CO₂ and S can be stored in the bubble as mineral precipitates (carbonates, sulfides, or sulfates; Venugopal et al. 2020; Schiavi et al. 2020); thus, MIs volatile contents can be underestimated if their presence is neglected. Diffusion of volatiles from the melt to the bubble can partially explain the variability of CO₂ contents in the glass (Fig. 4d). Assuming melt saturation in both H₂O and CO₂ and considering the inclusions with the highest volatile contents, we obtain the following trapping pressure range: 560–720 MPa (17–22 km depth) for the two most enriched basaltic olivine-hosted MIs (i.e., Ol6-15_MI1 with 3077 ± 744 ppm CO₂ and Ol6-5_MI1 with 3155 ± 190 ppm CO₂) based on the corrected CO₂ content (melt + bubble, see Supplementary Material S2 for details on the correction calculations); 490 ± 10 MPa (14.5–15 km depth) for the andesitic plagioclase-hosted MIs (i.e., PlM1_MI1 and Pl1_MI1 with 4430 and 4213 ppm CO₂, respectively). The latter is a minimal pressure estimate because the CO₂ content of the bubbles in plagioclase-hosted MIs is not known.

The major element compositions of MIs define two partially overlapping clusters: (1) with most olivine-hosted MIs being basaltic and representing parental magma and (2) plagioclase and pyroxene-hosted MIs (and olivine-hosted Ol6-12_MI1) exhibiting more evolved compositions, i.e., mainly andesitic. The MI major element contents seem to follow a differentiation trend from the parental to the more evolved magma, with a compositional gap between the two groups of MIs. The presence of this compositional trend suggests that the two groups of MIs share the same origin. To investigate the role played by fractional crystallization, we modeled the evolution of melt compositions using Petrolog3 software (Danyushevsky and Plechov 2011). Details on the models selected are provided in Supplementary Material S5. The starting composition is defined by the less evolved MI composition (i.e., M16 MI), and starting pressure is fixed at 720 MPa, corresponding to the maximal estimated trapping pressure, while the starting temperature is fixed at 1100 °C. The water was treated incompatibly as the effect of CO₂ on H₂O solubility is not taken into account in the models used. The two fractional crystallization models that better fit the major element data are shown in Fig. 3. In both models, the composition of the melt evolves by fractional crystallization of olivine, clinopyroxene, plagioclase, and magnetite in presence of H₂O. The first model (1) set at
<figure><img src="image_14.png" /><figcaption>**Fig. 8** Comparison of CO₂ contents in MIs of the Enco magma (red diamonds on the left side of the diagram, showing exceptionally high contents) with data of other arc volcanic systems. The diagram shows studies giving “uncorrected” CO₂ contents (measured only in the MI</figcaption></figure>

glass) and “corrected” CO₂ contents (i.e., taking into account the CO₂ content of both MI glass and bubble). The mean associated error with these data is ~15%

constant pressure (P = 720 MPa) partially reproduces the compositional variability of the mafic olivine-hosted MIs but stops at 56% crystallization, due to depletion of TiO₂, before reaching more evolved andesitic and dacitic compositions. Instead, the second model (2) of fractional crystallization associated to a fast decompression rate of the magma (–5 MPa·°C⁻¹) can reproduce the entire compositional range of major element data within about 75% of fractional crystallization. Despite the quite good fit obtained for the major element composition, fractional crystallization cannot be the only process responsible for the evolution from the basaltic end-member to the andesitic-dacitic one because of at least three main reasons. First, the final mineral compositions obtained in both models are not in agreement with the measured compositions in the investigated samples, i.e., Fo₆₄–An₆₁–En₅₄ and Fo₄₅–An₄₃–En₂₇ for models (1) and (2), respectively, against measured mineral compositions of Fo₇₅–八¹–An₄₄–五₁/An₇₆–En₄₂–六七. Possibly, as a significant CO₂ content in the melt and the presence (rare) of orthopyroxene in the investigated samples are not considered in the models used, the accuracy of the output may be low. Thus, some care should be taken with these modeling results. Second,

simple calculations demonstrate that variations in absolute trace element concentrations as well as in trace element ratios (Fig. 6) cannot be explained by a process of fractional crystallization alone. A good fit is only observed for Sr/Y versus Sr contents (Fig. 6c), probably due to the fact that Sr has a compatible behavior with plagioclase (Bindeman et al. 1998). Instead, fractional crystallization can partially account for the compositional variability (major and trace elements) of the olivine-hosted basaltic melts. Third, as the mafic and evolved end-members in the MIs were co-erupted, it is not likely that they were directly related by only fractional crystallization, with the host minerals being most probably phenocrysts. In fact, physical models taking into account the evolution of parameters such as magma density, viscosity, and crystal nucleation/growth rates (Brophy 2009) cannot explain the occurrence of an extreme degree of fractional crystallization (i.e., up to 75%) coupled to rapid magma ascent.

We consider trace elements to assess if the observed compositional trends are the result of a process of magma mixing. In the case of major elements, mixing lines between averaged compositions of olivine-hosted MIs and
plagioclase/pyroxene-hosted MIs were calculated (but not shown in Fig. 3 for clarity) and may explain the clusters observed. However, as fractional crystallization possibly occurs, especially in the case of basaltic melts, major elements do not show well-defined linear trends. PetroGram (Gündüz and Asan 2021) was used to model magma mixing for trace elements (Fig. 6). Highly to moderately incompatible elements like Th, Rb, Ce, Y, and Zr have low mineral-melt partition coefficients, while Sr is compatible in plagioclase, Co has strong affinity with olivine (and to a lesser extent, clinopyroxene), and V is slightly incompatible in clinopyroxene under slightly oxidizing conditions (e.g., see Laubier et al. 2014; Geochemical Earth Reference Model 2018). The magma mixing appears as a hyperbolic curve in incompatible/incompatible or incompatible/compatible versus incompatible trace element plots (Fig. 6). Incremental mixing from the most primitive melt measured (Ol6-14_MI1) to the most evolved end-member (represented by Pl5_MI1 sample) reproduces successfully trace element data. At the scale of the lapilli, mingling-like textures were observed with clearer portions in contact with darker portions. Thus, the compilation of sample observations, major, trace element data, and modeling suggests that the Enco eruption occurred as a result of mixing between a hot basaltic magma (represented by olivine-hosted MIs), which was evolving by fractional crystallization at depth, and a colder, more evolved, andesitic to dacitic magma (represented by plagioclase and pyroxene-hosted MIs).

## Volatiles in the Enco magma

The different textures in the MIs reflect various formation and/or post-entrapment processes. While the majority of the analyzed MIs show a regular shape and exhibit only one shrinkage bubble, the olivine-hosted inclusion Ol6-13_MI1 with a single large bubble (Fig. 2b) has the lowest water (0.48 wt.%) and CO₂ (225 ppm) contents, which suggest fast decompression of the inclusion (depressurization and then vesiculation or decrepitation) (Lowenstern 1995), and important volatiles loss. Among the studied inclusions, the olivine-hosted inclusion Ol6-5_MI1 (Fig. 7) has the highest measured CO₂ content in the glass (2836 ppm) and the highest corrected total (glass + bubble) CO₂ concentration of 3155 ± 190 ppm for this mineral phase. However, when taking into account the correction incertitude, the inclusion Ol6-15_MI1 may display the highest corrected CO₂ content with 3077 ± 744 ppm. Higher CO₂ contents have been measured in plagioclase-hosted MIs with up to 4430 ppm in PIM1_MI1 (uncorrected data, glass content only). A peculiar texture has been observed in this sample and also in Ol6-5_MI1 (Fig. 2c), Ol6-4_MI1, and Pl2_MI1, with numerous small bubbles evenly distributed that suggest fast exsolution of the volatiles and quenching of the

small bubbles hence formed. The presence of a large number of very small bubbles preserved in the glass of several inclusions indicates volatile supersaturation causing a high nucleation rate and is possible evidence of a very rapid syn-eruptive decompression followed by rapid quench, as shown by theory and experimental results (Mourtada-Bonnefoi and Laporte 2004; Toramaru 2006; Le Gall and Pichavant 2016). It is worth noting that the two olivine-hosted MIs showing this texture are those with the highest CO₂ contents in the glass, whereas high CO₂ concentration in the glass is not systematically associated to a large number of bubbles in plagioclase-hosted MIs. Moreover, both textures (i.e., bubble-rich, and single bubble) are observed in the same tephra unit. These observations can be explained by the high sensitivity of bubble nucleation to small variations in initial conditions (e.g., Mourtada-Bonnefoi and Laporte 2004; Toramaru 2006; Le Gall and Pichavant 2016), including the internal pressure and volatile concentration in the melt, but also to post-entrapment crystallization processes, thermal contraction, and elastic deformation (volumetric re-equilibration) (Schiavi et al. 2016; Drignon et al. 2020) that affect MIs before quenching.

An early degassing before the MIs entrapment could explain the large variability of CO₂, H₂O, and S contents observed at a given K₂O value (Fig. 4). Additionally, mixing between the mafic magma and the andesitic-dacitic one is able to affect the volatile contents, as observed in the case of S (Fig. 4b). Most data plotted in the CO₂ vs. H₂O contents diagram (Fig. 7) can be explained by a closed-system equilibrium degassing of the melt in presence of excess vapor. This is consistent with the strong explosivity behavior observed for the Enco eruption, as a closed-system equilibrium degassing would keep the volatiles under pressure (see also, Cassidy et al. 2018). The high CO₂ content measured in some MIs, in particular in two plagioclase-hosted MIs, also suggests the occurrence of CO₂ flushing from depth, as observed in several other volcanic systems (Blundy et al. 2010). We do not have a clear explanation for the high H₂O content (3.84 wt.%) of the inclusion Pl2_MI1 hosted in the Ca-rich plagioclase (An₇₆) that plots above the main trend, but we cannot exclude the possibility that it is a previously formed xenocryst or that it registered supersaturation of the magma in a deep portion of the system (Pichavant et al. 2013).

The H₂O, S, and Cl contents of the Enco magma are within the ranges given for other arc volcanic systems (Wallace 2005). The CO₂ contents in some of our MIs, however, are exceptionally high (up to 4430 ppm in the glass of PIM1_MI1, Table 1). The four most enriched MIs have CO₂ contents higher than those given by studies on other arc volcanoes in South and Central America (generally no more than 2000 ppm; see Fig. 8). In those studies, CO₂ in the bubbles was generally not taken into account. Two exceptions
are the MIs from San Cristobal (Nicaragua), for which corrected CO₂ contents are 718–1889 ppm versus 71–527 ppm for uncorrected data (Robidoux et al. 2016), and the MIs from Fuego (Guatemala) with corrected CO₂ contents of 513–2598 ppm versus 189–511 ppm for uncorrected data (Moore et al. 2015). Recent work by Rasmussen et al. (2020) gave corrected CO₂ contents for Fuego (Guatemala) ranging from 370 to 3409 ppm.

A possible mechanism of CO₂ enrichment at upper levels could be a gas flushing event (Caricchi et al. 2018), which seems in agreement with our observations (Fig. 7). According to the trace elements analysis, the magma of the Enco eruption of Mocho-Choshuenco volcano derived from a mixing between a deep and primitive basaltic magma rising from more than 22 km in depth and an evolved andesitic to dacitic magma stored in a hypothetic, shallower reservoir at a minimum estimated depth of 15 km. These estimates, as already stated previously, are minimal. However, only few studies in the literature report the composition of magmas erupted at Mocho-Choshuenco volcano. Moreover, no geophysical study, showing the position of a magmatic reservoir is, until now, available at Mocho-Choshuenco in order to confirm the minimal depth estimates obtained with MI investigations. Thus, our data are also important in providing this first vital information for future investigations that will better define magma storage conditions below Mocho-Choshuenco. The recharge of the magma chamber by this primitive CO₂-rich magma has possibly resulted in the enrichment in CO₂ and the triggering of the eruption. Hence, the high CO₂ content could have been the main cause of the explosivity of the Enco eruption event that generated highly mobile pyroclastic flows. In fact, other possible causes of explosivity seem not significant. The high temperature (> 1000 °C) and the low crystallinity (<30 vol.%, including macrocysts and microliths) of the erupted magma both suggest that the magma had a low viscosity and was not prone to an explosive regime (Lavallée et al. 2007). From the model of Vetere et al. (2006) for andesitic melts, with a water content of ~3.0 wt.%, the viscosity of the magma that fed the Enco eruption is estimated to be ~10²–10³ Pa.s, which is in the low viscosity range, while the high viscosity range is between 10⁸ and 10¹³ Pa.s (Vetere et al. 2006). As a comparison, the range of viscosities estimated at Mocho-Choshuenco for the Enco eruption is at least three orders of magnitude lower than the effusive eruption of Quizapu volcano (Chile), which emitted a 5 km³ andesitic-dacitic mingled lava flow in A.D. 1846–1847, with pre-eruptive H₂O contents of ~4–6 wt.% (Ruprecht and Bachmann 2010). The effusive behavior of this eruption, despite the relatively high viscosity of its magma, was explained by a late reheating following a mixing of dacite with andesite recharge;

leading to an acceleration of volatile diffusion; enhancing bubble nucleation, growth, and coalescence; and inhibiting brittle fragmentation (Ruprecht and Bachmann 2010). In the case of the Quizapu, the andesitic-dacitic mixing has prevented an explosive eruption to occur; as generally, silicic magma explosivity is strongly controlled by the viscosity (Cassidy et al. 2018). In contrast, the magma mixing which occurred at Mocho-Choshuenco seems to have rather favored the overpressure, possibly due to the high CO₂ contents, leading to the sub-Plinian Enco eruption. The CO₂ contents measured at Mocho-Choshuenco are also significantly higher than those recorded in some other Chilean volcanoes such as Llaima or Villarica, with CO₂ contents measured in MI glass below 500 ppm (see Fig. 8; Witter et al. 2004; Bouvet de Maisonneuve et al. 2012).

Finally, H₂O, S, and Cl contents are within the range of other arc volcanic systems around the world, with no exceptional concentrations measured (Wallace 2005). In this context, the low solubility of CO₂ (Bottinga and Javoy 1990; Le Gall and Pichavant 2016; Laporte et al. 2016) could have controlled the eruption dynamics. Moreover, CO₂ is known to significantly decrease the solubility of H₂O (Lowenstern 2001 and references therein). Water contents are also high; hence, high amounts of exsolved volatiles (CO₂ and H₂O) could have increased the magma buoyancy, generated high overpressures, and favored fast ascent of the magma and explosive fragmentation (Papale and Polacci 1999; Blundy et al. 2010; Freda et al. 2010; Aiuppa et al. 2011; Dallai et al. 2011; Allison et al. 2021). Owing to the exceptionally high CO₂ and fairly high H₂O contents we report, it appears that the magma of the Enco-MC15 eruption at Mocho-Choshuenco was particularly prone to explode violently despite its moderate silica content.

## Conclusions

We measured the volatile contents of MIs in minerals extracted from juveniles of a pyroclastic flow deposit formed during the andesitic Enco-MC15 eruption of Mocho-Choshuenco volcano. Though H₂O, S, and Cl contents are in ranges given for other arc volcanic systems, the CO₂ contents are exceptionally high compared to those given in published studies, with concentrations up to 4430 ppm. The initial CO₂ content was probably even higher, as suggested by the MI textures (i.e., presence of hundreds of small bubbles) as well as precipitated carbonates inside shrinkage bubbles. 3D Raman mapping allowed us to quantify the volume of solid carbonates trapped in inclusion bubbles and, thus, correct the CO₂ content in the MIs. As a whole, the volatile contents indicate that the basaltic MIs hosted in olivines were trapped at minimal depths of ~17–22 km, whereas trapping

<figure><img src="image_15.png" /></figure>
of andesitic to dacitic melts in plagioclase crystals occurred at ~15 km depth.

We conclude that the exceptionally high content and the low solubility of CO₂, combined with the fact that this volatile phase enhanced exsolution of H₂O, could have caused high overpressure in the magma reservoir and led the crystal-poor Enco andesitic magma of low viscosity to erupt explosively. A possible scenario considers that the deeper basaltic CO₂-rich magma ascended and recharged a shallower reservoir where the andesitic-dacitic magma was ponding. Once there, the mafic magma began to cool, exsolved volatiles, and mixed with the evolved magma, which resulted in rising of gas pressure (similar to gas sparging of Bachmann and Bergantz 2006); possible rapid bubble expansion associated with fast magma ascent finally triggered the explosive eruption. We note that the wide range of glass compositions and textures (i.e., scarce presence of phe-nocrysts and microliths, hundreds of bubbles preserved in the inclusion glass) of the studied samples seem to testify to a mild mixing process associated with a very rapid ascent of the recharging magma through the shallower storage levels, most likely caused by the strong enrichment in CO₂ of the basaltic melt. Further investigations on other explosive eruptions of Mocho-Choshuenco volcano are needed to confirm if the CO₂ enrichment is a typical feature throughout the life of the edifice or if it concerned only the Enco eruption. Our study supports the conclusions of earlier works that pointed out the possible correlation between high CO₂ contents and explosivity of poorly to moderately differentiated magmas. It suggests also that initial CO₂ contents in arc magmas may be commonly underestimated and that further work is required to better constrain the role that CO₂ may have in promoting explosive eruptions of basic to intermediate magmas.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00445-022-01550-y.

**Acknowledgements** The authors thank the Associate Editor, Tobias Fischer, as well as Tom Shea, and an anonymous reviewers for their constructive and helpful reviews. We thank Jean-Luc Devidal (LMV) for assistance with the electron microprobe and LA-ICP-MS and Etienne Deloule for his help with the ion probe at CRPG Nancy.

**Funding** Open access funding provided by University of Vienna. This research was financed by the French National Research Institute for sustainable Development (IRD), the French Government Laboratory of Excellence initiative No. ANR-10-LABX-0006, the Region Auvergne and the European Regional Development Fund, and the Chilean Science Commission (Fondecyt projects 1070162 and 1090387). This is Laboratory of Excellence ClerVolc contribution number 533.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are

included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

Aiuppa A, Burton M, Allard P, Caltabiano T, Giudice G, Gurrieri S, Liuzzo M, Salerno G (2011) First observational evidence for the CO₂-driven origin of Stromboli's major explosions. Solid Earth 2(2):135-142. https://doi.org/10.5194/se-2-135-2011

Allison CM, Roggensack K, Clarke AB (2021) Highly explosive basal- tic eruptions driven by CO₂ exsolution. Nat Commun 12:217. https://doi.org/10.1038/s41467-020-20354-2

Anderson AT, Brown GG (1993) CO₂ contents and formation pressures of some Kilauean melt inclusions. Am Mineral 78(7-8):794-803

Bachmann O, Bergantz GW (2006) Gas percolation in upper-crystal silicic crystal mushes as a mechanism for upward heat advection and rejuvenation of near-solidus magma bodies. J Volcanol Geotherm Res 149(1-2):85-102. https://doi.org/10.1016/j.jvolg-eores.2005.06.002

Balcone-Boissard H, Boudon G, Cioni R, Webster JD, Zdanowicz G, Orsi G, Civetta L (2016) Chlorine as a geobarometer for alkaline magmas: evidence from a systematic study of the eruptions of Mount Somma-Vesuvius. Sci Rep 6:21726. https://doi.org/10.1038/srep21726

Bindeman IN, Davis AM, Drake MJ (1998) Ion microprobe study of plagioclase-basalt partition experiments at natural concentration levels of trace elements. Geochim Cosmochim Acta 62(7):1175-1193. https://doi.org/10.1016/S0016-7037(98)00047-7

Blundy J, Cashman KV, Rust A, Witham F (2010) A case for CO₂-rich arc magmas. Earth Planet Sci Lett 290(3-4):289-301. https://doi.org/10.1016/j.epsl.2009.12.013

Botcharnikov RE, Behrens H, Holtz F (2006) Solubility and speciation of C-O-H fluids in andesitic melt at T = 1100-1300 °C and P = 200 and 500 MPa. Chem Geol 229(1-3):125-143. https://doi.org/10.1016/j.chemgeo.2006.01.016

Bottinga Y, Javoy M (1990) MORB degassing: bubble growth and ascent. Chem Geol 81(4):255-270. https://doi.org/10.1016/0009-2541(90)90050-H

Bouvet de Maisonneuve C, Dungan MA, Bachmann O, Burgisser A (2012) Insights into shallow magma storage and crystallization at Volcan Llaima (Andean Southern Volcanic Zone, Chile). J Volcanol Geotherm Res 211-212:76-91. https://doi.org/10.1016/j.jvolgeores.2011.09.010

Brophy JG (2009) Decompression and H₂O exsolution driven crys- tallization and fractionation: development of a new model for low-pressure fractional crystallization in calc-alkaline magmatic systems. Contrib Mineral Petrol 157:797-811. https://doi.org/10.1007/s00410-008-0365-2

Caricchi L, Sheldrake TE, Blundy J (2018) Modulation of magmatic processes by CO₂ flushing. Earth Planet Sci Lett 491:160-171. https://doi.org/10.1016/j.epsl.2018.03.042

Cassidy M, Manga M, Cashman K, Bachmann O (2018) Controls on explosive-effusive volcanic eruption styles. Nat Commun 9:2839. https://doi.org/10.1038/s41467-018-05293-3

Chen Y, Provost A, Schiano P, Cluzel N (2011) The rate of water loss from olivine-hosted melt inclusions. Contrib Mineral Petrol 162:625-636. https://doi.org/10.1007/s00410-011-0616-5

Clavero J, Roche O, Rivera A (2008) Recent unusual fines- depleted and highly-mobile andesitic pyroclastic flows from
Mocho-Choshuenco volcano, southern Andes. IAVCEI General Assembly, August 18–25, Reykjavik, Iceland.

Coltelli M, Del Carlo P, Vezzoli L (1998) Discovery of a Plinian basaltic eruption of Roman age at Etna volcano. Italy Geology 26(12):1095–1098. https://doi.org/10.1130/0091-7613(1998)026%3c1095:DOAPBE%3e2.3.CO;2

Costa F, Andreastuti S, Bouvet de Maisonneuve C, Pallister JS (2013) Petrological insights into the storage conditions, and magmatic processes that yielded the centennial 2010 Merapi explosive eruption. J Volcanol Geotherm Res 261:209–235. https://doi.org/10.1016/j.jvolgeores.2012.12.025

Cross JK, Tomlinson EI, Giordano G, Smith VC, De Benedetti AA, Roberge J, Manning CJ, Wulf S, Menzies MA (2014) High level triggers for explosive mafic volcanism: Albano Maar, Italy. Lithos 190–191:137–153. https://doi.org/10.1016/j.lithos. 2013.11.001

Dallai L, Cioni R, Boschi C, D'Oriano C (2011) Carbonate-derived CO₂ purging magma at depth: influence on the eruptive activity of Somma-Vesuvius. Italy Earth Planet Sci Lett 310(1-2):84-95. https://doi.org/10.1016/j.epsl.2011.07.013

Danyushevsky LV, Plechov P (2011) Petrolog3: integrated software for modeling crystallization processes. Geochemi Geophys Geosys 12(7):Q07021. https://doi.org/10.1029/2011GC003516

Deegan FM, Troll VR, Freda C, Misiti V, Chadwick JP, McLeod CL, Davidson JP (2010) Magma-carbonate interaction processes and associated CO₂ release at Merapi volcano, Indonesia: insights from experimental petrology. J Petrol 51:1027–1051. https://doi. org/10.1093/petrology/egg010

Drignon MJ, Arbaret L, Cluzel N, Nielsen RL, Bodnar RJ (2020) Experimentally induced volumetric re-equilibration of plagiocles-hosted melt inclusions. Geochem Geophys Geosyst 22(1):2020GC009357. https://doi.org/10.1029/2020GC009357

Edmonds M, Humphreys MC, Hauri EH, Herd RA, Wadge G, Rawson H, Ledden R, Plail M, Barclay J, Aiuppa A (2014) Pre-eruptive vapour and its role in controlling eruption style and longevity at Soufrière Hills Volcano. In: Wadge G, Robertson REA, Voight B (eds.) The Eruption of Soufrière Hills Volcano, Montserrat from 2000 to 2010, Geological Society, London, Memoirs, 39, pp. 291–315. https://doi.org/10.1144/M39.16

Esposito R, Lamadrid HM, Redi D, Steele-MacInnis M, Bodnar RJ, Manning CE, De Vivo B, Cannatelli C, Lima A (2016) Detection of liquid H₂O in vapour bubbles in reheated melt inclusions: implications for magmatic fluid composition and volatile budgets of magmas? Am Mineral 101:1691–1695. https://doi.org/10.2138/ am-2016-5689

Fontijn K, Rawson H, Van Daele M, Moernaut J, Abarzúa AM, Heirman K, Bertrand S, Pyle DM, Mather TA, De Batist M, Naranjo JA, Moreno H (2016) Synchronisation of sedimentary records using tephra: a postglacial tephrochronological model for the Chilcan Lake District. Quat Sci Rev 137:234–254. https://doi.org/10. 1016/j.quascirev.2016.02.015

Forsythe LM, Nielsen RL, Fisk MR (1994) High-field-strength element partitioning between pyroxene and basaltic to dacitic magmas. Chem Geol 117:107–125. https://doi.org/10.1016/0009-2541(94) 09124-4

Freda C, Gaeta M, Giaccio B, Marra F, Palladino DM, Scarlato P, Sottili G (2010) CO₂-driven large mafic explosive eruptions: the Pozzolane Rosse case study from the Colli Albani Volcanic District (Italy). Bull Volcanol 73(3):241–256. https://doi.org/10.1007/ s00445-010-0406-3

Frezzotti ML, Tecce F, Casagli A (2012) Raman spectroscopy for fluid inclusion analysis. J Geochem Explor 112:1–20. https://doi.org/ 10.1016/j.gexpol.2011.09.009

Gardner JE (2009) The impact of pre-existing gas on the ascent of explosively erupted magma. Bull Volcanol 71:835–844. https:// doi.org/10.1007/s00445-009-0276-8

Geochemical Earth Reference Model (GERM) Partition Coefficient (Kd) Database (2018) https://earthref.org/KDD/ (accessed 20th March 2018)

Gonnermann HM, Manga M (2007) The fluid mechanics inside a volcano. Annu Rev Fluid Mech 39:321–356. https://doi.org/10.1146/ annurev.fluid.39.050905.110207

Gündüz M, Asan K (2021) PetroGram: an Excel-based petrology program for modeling of magmatic processes. Geosci Front 12(1):81–92. https://doi.org/10.1016/j.gsf.2020.06.010

Hart SR, Dunn T (1993) Experimental cpx/melt partitioning of 24 trace elements. Contrib Mineral Petrol 113:1–8. https://doi.org/ 10.1007/BF00320827

Houghton BF, Wilson CJN, Del Carlo P, Catelli M, Sable JE, Carey JR (2004) The influence of conduit processes on changes in style of basaltic Plinian eruptions: Tarawera 1886 and Etna 122 BC. J Volcanol Geotherm Res 137(1–3):1–14. https://doi.org/10.1016/j. jvolgeores.2004.05.009

Huppert HE, Woods AW (2002) The role of volatiles in magma chamber dynamics. Nature 420(6915):493–495. https://doi.org/10. 1038/nature01211

Jugo P (2009) Sulfur content at sulfide saturation in oxidized magmas. Geology 37 (5):415–418. https://doi.org/10.1130/G25527A.1

Lamadrid HM, Moore LR, Moncada D, Rimstidt JD, Burruss RC, Bodnar RJ (2017) Reassessment of the Raman CO₂ densimeter. Chem Geol 450:210–222. https://doi.org/10.1016/j.chemgeo.2016.12. 034

Laporte D, Hardiagon M, Provost A (2016) Ascent and degassing of magmas: the conditions of homogeneous bubble nucleation in the system basalt-CO₂. Experimental Mineralogy, Petrology and Geochemistry, EMPG XV - International Conference, Zürich, Switzerland, June 5–8, 2016.

LaTourrette TZ, Burnett DS (1992) Experimental determination of U and Th partitioning between clinopyroxene and natural and synthetic basaltic liquid. Earth Planet Sci Lett 110:227–244. https:// doi.org/10.1016/0012-821X(92)90050-6

Laubier M, Grove TL, Langmuir CH (2014) Trace element mineral/ melt partitioning for basaltic and basaltic andesitic melts: an experimental and laser ICP-MS study with application to the oxidation state of mantle source regions. Earth Planet Sci Lett 392:265–278. https://doi.org/10.1016/j.epsl.2014.01.053

Lavallée Y, Hess KU, Cordonnier B, Dingwell DB (2007) Non-Newtonian law for highly crystalline dome lavas. Geology 35(9):843– 846. https://doi.org/10.1130/G23594A.1

Le Gall N, Pichavant M (2016) Homogeneous bubble nucleation in H₂O- and H₂O-CO₂-bearing basaltic melts: results of high temperature decompression experiments. J Volcanol Geotherm Res 327:604–621. https://doi.org/10.1016/j.jvolgeores.2016.10.004

Le Voyer M (2009) Rôle des fluides dans la genèse des magmas d'arc: analyses in-situ des éléments volatils et des isotopes du bore dans les inclusions magmatiques des olivines primitives. PhD Thesis, Geochemistry, Université Blaise Pascal, Clermont-Ferrand II, France, 276p. https://tel.archives-ouvertes.fr/tel-00453889

Le Voyer M, Kelley KA, Cottrell E, Hauri EH (2017) Heterogeneity in mantle carbon content from CO₂-undersaturated basalts. Nat Commu 8:14062. https://doi.org/10.1038/ncomms14062

Lesne P, Scaillet B, Pichavant M (2015) The solubility of sulfur in hydrous basaltic melts. Chem Geol 418:104–116. https://doi.org/ 10.1016/j.chemgeo.2015.03.025

Lloyd AS, Plank T, Ruprecht P, Hauri E, Rose W (2013) Volatile loss from melt inclusions in pyroclasts of differing sizes. Contrib Mineral Petrol 165(1):129–153. https://doi.org/10.1007/ s00410-012-0800-2

Lowenstern JB (1995) Application of silicate-melt inclusions to the study of magmatic volatiles. In: Thompson JFH (ed.) Magmas, Fluids and Ore Deposits. Mineralogical Association of Canada, pp. 71–100.
Lowenstern JB (2001) Carbon dioxide in magmas and implications for hydrothermal systems. Miner Deposita 36(6):490-502. https://doi.org/10.1007/s001260100185

McPhie J (1990) Phreatomagmatic and phreatic fall and surge deposits from explosions at Kilauea Volcano, Hawaii, 1790 A.D. Keanakakoi Ash Member. Bull Volcanol. 52(5):334-354. https://doi.org/10.1007/BF00302047

Métrich N, Clocchiatti R (1996) Sulfur abundance and its speciation in oxidized alkaline melts. Geochim Cosmochim Acta 60:4151-4160. https://doi.org/10.1016/S0016-7037(96)00229-3

Métrich N, Allard P, Spilliaert N, Andronico D, Burton M (2004) 2001 flank eruption of the alkali- and volatile-rich primitive basalt responsible for Mount Etna's evolution in the last three decades. Earth Planet Sci Lett 228:1-17. https://doi.org/10.1016/j.epsl.2004.09.036

Moore L, Gazel E, Tuohy R, Lloyd A, Esposito R, Steele-MacInnis M, Hauri EH, Wallace PJ, Plank T, Bodnar RJ (2015) Bubbles matter: an assessment of the contribution of vapor bubbles to melt inclusion volatile budgets. Am Mineral 100(4):806-823. https://doi.org/10.2138/am-2015-5036

Moreno H, Lara L (2007) Geologia del Complejo Volcánico Mocho-Choshuenco, Carta Geológica de Chile. Serie Geología Básica 107:1-27

Moretti R, Métrich N, Arienzo I, Di Renzo V, Aiuppa A, Allard P (2018) Degassing vs. eruptive styles at Mt. Etna volcano (Sicily, Italy). Part I: Volatile stocking, gas fluxing, and the shift from low-energy to highly explosive basaltic eruptions. Chem Geol 482:1-17. https://doi.org/10.1016/j.chemgeo.2017.09.017

Moune S, Sigmarsson O, Schiano P, Thordarson T, Keiding JK (2012) Melt inclusion constraints on the magma source of Eyjafjallajökull 2010 flank eruption. J Geophys Res-Solid Earth 117:B00C07. https://doi.org/10.1029/2011JB008718

Mourtada-Bonnefoi CC, Laporte D (2004) Kinetics of bubble nucleation in a rhyolitic melt: an experimental study of the effect of ascent rate. Earth Planet Sci Lett 218(3-4):521-537. https://doi.org/10.1016/S0012-821X(03)00684-8

Myers ML, Geist DJ, Rowe MC, Harpp KS, Wallace PJ, Dufek J (2014) Replenishment of volatile-rich mafic magma into a degassed chamber drives mixing and eruption of Tungurahua Volcano. Bull Volcanol 76(11):1-17. https://doi.org/10.1007/s00445-014-0872-0

Newman S, Lowenstern JB (2002) VolatileCalc: a silicate melt-H₂O-CO₂ solution model written in visual basic for Excel. Comput and Geosci 28(5):597-604. https://doi.org/10.1016/S0098-3004(01)00081-4

Papale P, Polacci M (1999) Role of carbon dioxide in the dynamics of magma ascent in explosive eruptions. Bull Volcanol 60(8):583-594. https://doi.org/10.1007/s004450050253

Papale P, Moretti R, Barbato D (2006) The compositional dependence of the saturation surface of H₂O+CO₂ fluids in silicate melts. Chem Geol 229(1):78-95. https://doi.org/10.1016/j.chemgeo.2006.01.013

Pichavant M, Di Carlo I, Rotolo SG, Scaillet B, Burgisser A, Le Gall N, Martel C (2013) Generation of CO₂-rich melts during basalt magma ascent and degassing. Contrib Mineral Petrol 166:545-561. https://doi.org/10.1007/s00410-013-0890-5

Putirka KD (2008) Thermometers and barometers for volcanic systems. Rev Mineral Geochem 69(1):61-120. https://doi.org/10.2138/rmg.2008.69.3

Putirka K, Johnson M, Kinzler R, Walker D (1996) Thermobarometry of mafic igneous rocks based on clinopyroxene-liquid equilibria, 0-30 kbar. Contrib Mineral Petrol 123:92-108. https://doi.org/10.1007/s004100050145

Ramsey CB (2009) Dealing with outliers and offsets in radiocarbon dating. Radiocarbon 51(3):1023-1045. https://doi.org/10.1017/S0033822200034093

Rasmussen DJ, Plank TA, Wallace PJ, Newcombe ME, Lowenstern JB (2020) Vapor-bubble growth in olivine-hosted melt inclusions. Am Min 105:1898-1919. https://doi.org/10.2138/am-2020-7377

Rawson H, Naranjo JA, Smith V, Fontijn K, Pyle DM, Mather TA, Moreno H (2015) The frequency and magnitude of post-glacial explosive eruptions at Volcán Mocho-Choshuenco, southern Chile. J Volcanol Geotherm Res 299:103-129. https://doi.org/10.1016/j.jvolgeores.2015.04.003

Rawson H, Keller TM, Fontijn K, Pyle DM, Mather TA, Smith VC, Naranjo JA (2016) Compositional variability in mafic arc magmas over short spatial and temporal scales: evidence for the signature of mantle reactive melt channels. Earth Planet Sci Lett 456:66-77. https://doi.org/10.1016/j.epsl.2016.09.056

Robidoux P, Aiuppa A, Rotolo SG, Rizzo AL, Hauri EH, Frezzotti ML (2016) Volatile contents of mafic-to-intermediate magmas at San Cristóbal Volcano in Nicaragua. Lithos 272-273:147-163. https://doi.org/10.1016/j.lithos.2016.12.002

Robidoux P, Frezzotti ML, Hauri EH, Aiuppa A (2018) Shrinkage Bubbles: The C-O-H-S Magmatic Fluid System at San Cristóbal Volcano. J Petrol 59(11):2093-2122. https://doi.org/10.1093/petrology/egy092

Roeder PL, Emslie RF (1970) Olivine-liquid equilibrium. Contrib Mineral Petrol 29(4):275-289. https://doi.org/10.1007/BF00371276

Rose-Koga EF, Koga KT, Devidal J-L, Shimizu N, Le Voyer M, Dalou C, Döbeli M (2020) In-situ measurements of magmatic volatile elements, F, S, and Cl, by electron microprobe, secondary ion mass spectrometry, and heavy ion elastic recoil detection analysis. Am Min 105(5):616-626. https://doi.org/10.2138/am-2020-7221

Ruprecht P, Bachmann O (2010) Pre-eruptive reheating during magma mixing at Quizapu volcano and the implications for the explosiveness of silicic arc volcanoes. Geology 38(10):919-922. https://doi.org/10.1130/G31110.1

Sable JE, Houghton BF, Del Carlo P, Coltelli M (2006) Changing conditions of magma ascent and fragmentation during the 122 BC basaltic Plinian eruption: evidence from cast micro textures. J Volcanol Geotherm Res 158:333-354. https://doi.org/10.1016/j.jvolgeores.2006.07.006

Scandone R, Cashman KV, Malone SD (2007) Magma supply, magma ascent and the style of volcanic eruptions. Earth Planet Sci Lett 253:513-529. https://doi.org/10.1016/j.epsl.2006.11.016

Schiano P (2003) Primitive mantle magmas recorded as silicate melt inclusions in igneous minerals. Earth-Sci Rev 63(1-2):122-144. https://doi.org/10.1016/S0012-8252(03)00034-5

Schiavi F, Provost A, Schiano P, Cluzel N (2016) P-V-T-X evolution of olivine-hosted melt inclusions during high-temperature homogenization treatment. Geochim Cosmochim Acta 172:1-21. https://doi.org/10.1016/j.gca.2015.09.025

Schiavi F, Bolfan-Casanova N, Buso R, Laumonier M, Laporte D, Medjoubi K, Venugopal S, Gómez-Ulla A, Cluzel N, Hardiagon M (2020) Quantification of magmatic volatiles by Raman microscopy of glass inclusion-hosted bubbles. Geochem Perspect Lett 16:17-24. https://doi.org/10.7185/geochemlet.2038

Sparks SRJ, Sigurdsson H, Wilson L (1977) Magma mixing: a mechanism for triggering acid explosive eruptions. Nature 267:315-318. https://doi.org/10.1038/267315a0

Spilliaert N, Allard P, Metrich N, Sobolev AV (2006) 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). J Geophys Res 111:B04203. https://doi.org/10.1029/2005JB003934

Stern CR, Moreno H, López-Escobar L, Clavero JE, Lara LE, Naranjo JA, Parada MA, Skewes MA (2007) Chilean Volcanoes. In: Moreno T, Gibbons W (eds.) The geology of Chile, Geological Society, London, pp 147-178. https://doi.org/10.1144/GOCH.5

Swanson DA, Rose TR, Fiske RS, McGeehin JP (2012) Keanakāko'i Tephra produced by 300 years of explosive eruptions following
collapse of Kīlauea's caldera in about 1500 CE. J Volcanol Geotherm Res 215-216:8-25. https://doi.org/10.1016/j.jvolgeores.2011.11.009

Szramek LA (2016) Mafic plinian eruptions: is fast ascent required? J Geophys Res-Solid Earth 121(10):7119-7136. https://doi.org/10.1002/2016JB013208

Toplis MJ (2005) The thermodynamics of iron and magnesium partitioning between olivine and liquid: criteria for assessing and predicting equilibrium in natural and experimental systems. Contrib Mineral Petrol 149:22-39. https://doi.org/10.1007/s00410-004-0629-4

Toramaru A (2006) BND (bubble number density) decompression rate meter for explosive volcanic eruptions. J Volcanol Geotherm Res 154(3-4):303-316. https://doi.org/10.1016/j.jvolgeores.2006.03.027

Tucker JM, Hauri EH, Pietruszka AJ, Garcia MO, Marske JP, Trusdell FA (2019) A high carbon content of the Hawaiian mantle from olivine-hosted melt inclusions. Geochim Cosmochim Acta 254:156-172. https://doi.org/10.1016/j.gca.2019.04.001

Venugopal S, Moune S, Williams-Jones G (2016) Investigating the subsurface connection beneath Cerro Negro Volcano and the El Hoyo complex, Nicaragua. J Volcanol Geotherm Rese 325:211-224. https://doi.org/10.1016/j.jvolgeores.2016.06.001

Venugopal S, Schiavi F, Moune S, Bolfan-Casanova N, Druitt T, Williams-Jones G (2020) Melt inclusion vapour bubbles: the hidden reservoir for major and volatile elements. Sci Rep 10:9034. https://doi.org/10.1038/s41598-020-65226-3

Vetere F, Behrens H, Holtz F, Neuville DR (2006) Viscosity of andesitic melts—new experimental data and a revised calculation model. Chem Geol 228:233-245. https://doi.org/10.1016/j.chemgeo.2005.10.009

Wallace PJ (2005) Volatiles in subduction zone magmas: concentration and fluxes based on melt inclusion and volcanic gas data. J Volcanol Geotherm Res 140(1-3):217-240. https://doi.org/10.1016/j.jvolgeores.2004.07.023

Wallace PJ, Carmichael ISE (1992) Sulfur in basaltic magmas. Geochim Cosmochim Acta 56(5):1863-1874. https://doi.org/10.1016/0016-7037(92)90316-B

Wallace PJ, Carmichael ISE (1994) S speciation in submarine basaltic glasses as determined by measurements of S K-α X-ray wavelength shifts. Am Mineral 79(1-2):161-167

Wallace PJ, Kamenetsky VS, Cervantes P (2015) Melt inclusion CO₂ contents, pressures of olivine crystallization, and the problem of shrinkage bubbles. Am Mineral 100(4):787-794. https://doi.org/10.2138/am-2015-5029

Williamson BJ, Di Muro A, Horwell CJ, Spieler O, Llewellin EW (2010) Injection of vesicular magma into an andesitic dome at the effusive-explosive transition. Earth Planet Sci Lett 295(1-2):83-90. https://doi.org/10.1016/j.epsl.2010.03.027

Witham F, Blundy J, Kohn SC, Lesne P, Dixon J, Churakov SV, Botcharnikov RE (2012) SolEx: a model for mixed COHSCI-volatile solubilities and exsolved gas compositions in basalt. Comput and Geosci 45:87-97. https://doi.org/10.1016/j.cageo.2011.09.021

Witter J, Kress V, Delmelle P, Stix J (2004) Volatile degassing, petrology, and magma dynamics of the Villarrica lava lake Southern Chile. J Volcanol Geotherm Res 134(4):303-337. https://doi.org/10.1016/j.jvolgeores.2004.03.002

Zanetti A, Tiepolo M, Oberti R, Vannucci R (2004) Trace-element partitioning in olivine: modelling of a complete data set from a synthetic hydrous basanite melt. Lithos 75:39-54. https://doi.org/10.1016/j.lithos.2003.12.022

Zimanowski B, Büttner R, Dellino P, White JDL, Wohletz KH (2015) Magma-Water Interaction and Phreatomagmatic Fragmentation. In: Sigurdsson H, Hughton B, McNutt SR, Rymer H, Stix J (eds) The Encyclopedia of Volcanoes Second Edition, Academic Press, San Diego, pp 473-484. https://doi.org/10.1016/B978-0-12-385938-9.00026-2

Zhang Y (1999) A criterion for the fragmentation of bubbly magma based on brittle failure theory. Nature 402:648-650. https://doi.org/10.1038/45210