Insights Into the Mechanisms of Phreatic Eruptions From Continuous High Frequency Volcanic Gas Monitoring: Rincón de la Vieja Volcano, Costa Rica

Understanding the trigger mechanisms of phreatic eruptions is key to mitigating the effects of these hazardous but poorly forecastable volcanic events. It has recently been established that high-rate volcanic gas observations are potentially very suitable to identifying the source processes driving phreatic eruptions, and to eventually detecting precursory changes prior to individual phreatic blasts. In February-May 2017, we deployed a Multi-GAS instrument to continuously monitor gas concentrations in the crater lake plume of Rincón de la Vieja, a remote and poorly monitored active volcano in Costa Rica, site of frequent phreatic/phreatomagmatic eruptions. Forty-two phreatic/phreatomagmatic eruptions were seismically recorded during our investigated period, 9 of which were also recorded for gas by the Multi-GAS. To the best of our knowledge, these represent the ﬁrst instrumentally measured gas compositions during individual phreatic/phreatomagmatic explosions at an active volcano. Our results show that during background quiescent degassing the Rincón de la Vieja crater lake plume was characterized by high CO 2 /SO 2 ratios of 64 ± 59 and H 2 S/SO 2 ratios of 0.57 ± 0.20. This composition is interpreted as reﬂecting hydrothermal (re)processing of magma-sourced gas in the sub-limnic environment. Phreatic blasts were recorded by the Multi-GAS as brief (1–2 min long) pulses of elevated gas mixing ratios (up to ∼ 52 ppmv SO 2 and > 3,000 ppmv CO 2 ), or more than an order of magnitude higher than during background degassing ( ∼ 1 ppmv SO 2 and ∼ 450 ppmv CO 2 ). During the phreatic eruption(s), the H 2 S/SO 2 ratio was systematically lower ( < 0.18) than during background degassing, but the CO 2 /SO 2 ratio remained high (and variable), ranging from 37 to 390. These S-poor compositions for the eruptive gas imply extensive processing of the source magmatic gas during pre-eruptive hydrothermal storage, likely by deposition of native S and/or triggered by accumulation of magmatic-hydrothermal gases beneath a hydrothermal seal. We claim that real-time Multi-GAS monitoring is urgently needed at other crater lake-hosting volcanoes (e.g., Ruapehu, Aso), where phreatic eruptions may similarly be preceded by phases of reduced S degassing at the surface.

Understanding the trigger mechanisms of phreatic eruptions is key to mitigating the effects of these hazardous but poorly forecastable volcanic events. It has recently been established that high-rate volcanic gas observations are potentially very suitable to identifying the source processes driving phreatic eruptions, and to eventually detecting precursory changes prior to individual phreatic blasts. In February-May 2017, we deployed a Multi-GAS instrument to continuously monitor gas concentrations in the crater lake plume of Rincón de la Vieja, a remote and poorly monitored active volcano in Costa Rica, site of frequent phreatic/phreatomagmatic eruptions. Forty-two phreatic/phreatomagmatic eruptions were seismically recorded during our investigated period, 9 of which were also recorded for gas by the Multi-GAS. To the best of our knowledge, these represent the first instrumentally measured gas compositions during individual phreatic/phreatomagmatic explosions at an active volcano. Our results show that during background quiescent degassing the Rincón de la Vieja crater lake plume was characterized by high CO 2 /SO 2 ratios of 64 ± 59 and H 2 S/SO 2 ratios of 0.57 ± 0.20. This composition is interpreted as reflecting hydrothermal (re)processing of magma-sourced gas in the sub-limnic environment. Phreatic blasts were recorded by the Multi-GAS as brief (1-2 min long) pulses of elevated gas mixing ratios (up to ∼52 ppmv SO 2 and INTRODUCTION Phreatic eruptions are among the most unpredictable and hazardous volcanic phenomena (Mastin and Witter, 2000;Browne and Lawless, 2001). These blasts involve the violent, explosive discharge of pressurized pockets of external (nonvolcanic, mostly meteoric) water, and are therefore particularly common at "wet" volcanoes whose summits are topped by crater lakes (Rouwet and Morrissey, 2015;Stix and de Moor, 2018). In addition to the large availability of exogenous water, active crater lakes are especially prone to developing phreatic eruptions owing to the presence of a persistent heat source (conductive/convective heating from shallow magma, and/or rising magmatic volatiles), and the frequent formation of permeability barriers at the lake bottom (e.g., impermeable layers of precipitated native sulfur or alteration minerals) that favor gas accumulation (Christenson et al., 2010;Christenson and Tassi, 2015;Delmelle and Bernard, 2015, and references cited therein). The exact mechanisms driving crater lake breaching eruptions are still not entirely understood, and the respective roles of the magmatic and sublimnic hydrothermal systems in triggering the eruptions are still a matter of debate (Takano et al., 1994;Christenson and Tassi, 2015;Rouwet and Morrissey, 2015).
Phreatic eruptions have typically occurred in the past without being preceded by any obvious precursor, as recently dramatically demonstrated by the deadly Ontake (in 2014, Oikawa et al., 2016) and Kusatsu-Shirane (in 2018) eruptions in Japan. The recent technical advances in real-time observation of lake plume gas compositions (Di Napoli et al., 2013;Shinohara et al., 2015;Tamburello et al., 2015;de Moor et al., 2016a;Gunawan et al., 2016) bring a new perspective on eruption forecasting. Using measurements from a permanently installed "lake" Multi-GAS (Multi-component Gas Analyser System; Aiuppa et al., 2006Aiuppa et al., , 2010Aiuppa, 2015), de Moor et al. (2016a) identified for the first time systematic short-term (days to weeks) variations in plume gas compositions prior to individual phreatic explosions at Laguna Caliente crater lake, Poás volcano (Costa Rica). These observations indicated an increase of magmatic volatiles input prior to individual phreatic blasts, which demonstrates the potential of high-frequency real-time gas monitoring.
Here, we report on the results of 3 months of instrumental monitoring of volcanic gas composition at Rincón de la Vieja volcano (10.49 N, 85.19 W), in the Costa Rican segment of the Central American volcanic arc (CAVA). Rincón de la Vieja is one of the most active and remote volcanoes in Costa Rica (Barquero and Segura, 1983;Alvarado et al., 1992), and hosts a highly acidic crater lake (Tassi et al., 2005(Tassi et al., , 2009, the source of recurrent phreatomagmatic to phreatic eruptions (Boudon et al., 1996). The intense (SO 2 > 60 tons/day; de Moor et al., 2017) gas emissions have only occasionally been studied in the past (Tassi et al., 2005(Tassi et al., , 2009Aiuppa et al., 2014;de Moor et al., 2017), due to the limited accessibility of the volcano summit. Our study here provides the first near-continuous gas dataset taken during a period of recurrent phreatic activity at Rincon de la Vieja, including the first measurement of the syn-explosive gas phase. The aim is to use our novel gas observations to derive new insights into the mechanisms driving crater lake phreatic explosions and assess the potential for forecasting eruptions at Rincón de la Vieja.

RINCÓN DE LA VIEJA VOLCANO
Rincón de la Vieja volcano is part of the Guanacaste volcanic range (Figure 1), a Quaternary magmatic range related to subduction of the Cocos plate underneath the Caribbean plate (Carr, 1984;Carr et al., 1990;DeMets, 2001;DeMets et al., 2010). The Guanacaste volcanic range consists of four andesitic central edifices (Orosí-Cacao, Rincón de la Vieja-Santa María, Miravalles, and Tenorio-Montezuma). Before construction of the Quaternary andesitic chain, intense explosive silicic volcanism generated a series of collapse events and associated calderas (Molina et al., 2014). The Rincón de la Vieja-Santa María volcanic complex (maximum elevation 1,916 m) was constructed within one of these calderas, the 120 km 2 Cañas Dulces Caldera (Molina and Martí, 2016). The Cañas Dulce caldera hosts several hydrothermal fields (Giggenbach and Correales, 1992), including several thermal manifestations on the Rincón de la Vieja massif itself (Tassi et al., 2005), aligned along a NW-SE trend running roughly parallel to the volcano's axis.
Phreatic/phreatomagmatic activity resumed in September 2011, after 13 years of quiescence. During 2012-2013, Rincón de la Vieja exhibited low seismicity and very infrequent eruptions. From September 2014, eruptive activity escalated to an average of 25 phreatic-phreatomagmatic events per month, peaking in October 2015 and March 2016 when 43 and 220 eruptions were recorded, respectively. Visual binocular microscope observations (made by G. Avard at OVSICORI-UNA, following the procedure described in Alvarado et al., 2016) revealed that erupted ash fragments (250-500 µm portion) contain well-preserved greenish minerals and a small portion (∼5%) of fresh-looking glassy and vesicular shards. These observations were interpreted as indicative of the involvement of shallow magma in 2015-2016, as is the case at many volcanoes with phreatic eruptions (e.g., Stix and de Moor, 2018). Native sulfur fragments were also systematically observed in the erupted products of 2015-2017, especially in the 1-2 mm ash grain-size fraction (G. Avard, pers. comm.). Sulfur-rich minerals were also repeatedly observed as veins and fillings in fractures of large ejecta and blocks dispersed throughout the crater. The frequency and seismic energy of eruptions increased in early 2017. During our gas monitoring interval (February 3-May 9, 2017) 42 seismic signals associated to phreatic/phreatomagmatic eruptions were recognized, and 9 of these eruptions were also identified geochemically (Table 1; see below). Some of these 9 eruptions (eruptions 1-3 and 9; Figure 3) were preceded by a long period (LP) seismic signal several seconds before the explosion. In almost all cases (eruptions 1, 3, 4-7), pulses of spasmodic tremor occurred a few hours before the eruption. Two major phreatic/phreatomagmatic eruptions (generating hot lahars that traveled outside the crater to the north) took place on May 23 and on June 11 (Global Volcanism Program, 2017), just a few weeks after our Multi-GAS stopped acquiring. The explosion on June 11 generated a 1-2 km high plume, and ejected coarse materials to the W and NW onto the upper N flank. This event, the largest (based on seismic energy) registered during 2011-2018, was preceded by a LP swarm 2 h before the explosion. Material collected from the May 23 lahar down the river contained only ∼1% juvenile component, but fragments erupted from the June 11 eruption (also collected from a lahar) contained ∼44% juvenile material (OVSICORI-UNA). Phreatic activity has intermittently continued until the time of writing (late 2018).

MATERIALS AND METHODS
The gas dataset we report on in this study was recorded using two distinct fully autonomous Multi-GAS instruments, one FIGURE 2 | (a) Satellite image of Rincón de la Vieja volcano showing location of the Multi-GAS measurement site, the VORI seismic station, and the Radio Bridge, map data: Google, Landsat/Copernicus image, SIO, NOAA, U.S. Navy, NGA, GEBCO; (b) The Rincón de la Vieja active crater, with convection cells and abundant gray sulfur spherules floating on the lake surface; (c) Water and sulfur deposits in a canyon on the northern slope of Rincón de la Vieja. designed by the Instituto Nazionale di Geofisica e Vulcanologia (INGV) and one designed by the U.S. Geological Survey (USGS). The two instruments (INGV and USGS-campaign) operated sequentially at the same site, located on the northern inner wall of the active crater ∼1 m downwind the lake shore (Lat. 10.8327 Long. −85.3355; see Figure 2). This site was selected because its position ensured regular fumigation by the lake gas plume, while minimizing any potential gas contribution from low-temperature crater fumaroles (not observed at the measurement site during our study). The Multi-GAS instruments were powered by an external (12 V, 40 Ah) battery and 3 solar panels in order to ensure proper operation in all climatic conditions. A PVC tube with two filters served as gas inlet, and allowed gas to be pumped (at 1.2 L/m) inside the Multi-GAS instruments. The first (INGV-type) Multi-GAS was installed on February 3, 2017 and operated until March 17, 2017, when it was dismantled (for use at a different volcano) and replaced by a second (USGS-campaign type; Gunawan et al., 2016) Multi-GAS. This latter operated until May 09, 2017, when it stopped transmitting data before being finally destroyed during the May 23 lahar-producing phreatic blast (Global Volcanism Program, 2017).
The two Multi-GAS instruments used very similar sensor kits (see below) and were re-calibrated at OVSICORI with the same standard gases before the installation, showing very similar response. The INGV-type Multi-GAS was recalibrated after returning from the field (in late March 2017), and the original calibration was reproduced within sensor precision. Both Multi-GAS instruments measured SO 2 and H 2 S mixing ratios (precision within ± 15% at 2σ; Lewicki et al., 2017) with the same specific electrochemical sensors (models TD2G-1A and TC4E-1A, respectively; all from City Technology and with ± 5% repeatability). Interference of SO 2 gas on the H 2 S sensor (15 %) was determined during calibration procedure and corrected with the Ratiocalc software. CO 2 mixing ratios were measured using an on-board spectrometer (INGV-type: Gascard EDI030102NG, measurement range = 0-3,000 ppmv, precision, ± 3% at 2σ; USGS-type LI-COR LI-840A, measurement range = 0-5,000 ppmv, precision, ± 1.5% at 2σ). The INGV-type Multi-GAS instrument, controlled by an arduino2 datalogger, acquired data at 0.1 Hz rate during 4 sampling periods per day, of 30 min duration each. In the USGS-type Multi-GAS, all data acquisition and scheduling was controlled by a Campbell Scientific data logger (CR1000 with NL115 module) outfitted with 2 GB of onboard memory. The station completed four data acquisition cycles per day, in which data were acquired for 1 h at a 1 Hz rate. Before and after each measurement cycle, the station performed sensor baseline checks by activating two miniature 3-way teflon solenoid valves (Cole-Parmer WU-01540-11) that formed a closed-loop and recirculated trapped sample gases through soda lime and desiccant for 3 min to remove acid gases and water vapor in order to provide a measure of within-run baseline drift and to clean and dry the instrument prior to shutting down. For both Multi-GAS instruments, data were telemetered (at the end of each acquisition cycle) via radio modem (Xetawave) to OVSICORI using a radio bridge installed on the northern base of the volcano (Figure 2). The acquired gas mixing ratios data were post-processed using the Ratiocalc software (Tamburello,  Nine of these eruptions were also recorded by the Multi-GAS (the "Multi-GAS peak Time" identifies the temporal window of Multi-GAS over which the syn-explosive gas was detected). The nine eruptions are numbered as in Figure 3.

2015)
to derive gas ratios between volatile couples (CO 2 /SO 2 and H 2 S/SO 2 ) using linear regression. In order to derive accurate gas ratios, the different response time (∼2-4 s) between the spectroscopic sensor and the electrochemical sensors was corrected for during data processing with the Ratiocalc software.
No pressure correction was applied to the electrochemical sensors owing to very similar altitude between calibration site (Heredia, ∼1,200 m) and Rincon's summit (∼1,500 m). Based on laboratory tests with standard gases, errors in the derived ratios are typically ≤15%. Unfortunately, the plume was condensing at all conditions encountered during observations, so it was not possible to determine a robust volcanic H 2 O signal. Seismic activity was characterized using the reference station VORI, operated by the Observatorio Sismológico y Vulcanológico de Arenal y Miravalles (OSIVAM) of the Instituto Costarricense de Electricidad (ICE). This station is located 1.8 km FIGURE 3 | Time-series of (a-c) gas mixing ratios (in ppmv); (d) RSEM, (e) eruptive energy. Note that Multi-GAS data (CO 2 , SO 2, and H 2 S mixing ratios) are not continuous but obtained during 4 daily cycles, each 30-60 min long. The shaded yellow bands identify the 9 explosions (numbered #1 to #9) during which the Multi-GAS was operating, allowing the syn-explosive gas composition to be determined. The five light gray areas identify all periods with no MultiGAS data. The small red insets in the CO 2 plot identify the expanded portions shown in Figures 4, 6.
SW of the active crater, and is equipped with a digital 3component broadband TRIMBLE REFTEK 151B sensor with flat response from 0.016 (60 s) to 50 Hz. Signals were sampled at 100 Hz. Close inspection of seismic records was carried out to extract discrete events and, in the particular case of eruptive signals, to estimate the associated seismic energy (E seismic ). Seismic energy was calculated using the relation for bodywaves generated by an isotropic source at the top of a homogeneous half space (Boatwright, 1980;Arámbula-Mendoza et al., 2011): where r is the distance from the source to the seismic station, ρ is rock density, c is the P wave velocity, S is the site effect, U is seismic amplitude and A is attenuation. We used r = 1,800 m, c = 1,500 m/s, ρ = 2,600 kg/m 3 . The site and attenuation effects were neglected for simplicity. The temporal evolution of seismic activity was explored using Real-time Seismic Energy Measurement, or RSEM (De la Cruz-Reyna and Reyes-Davila, 2001), which considers the square of the amplitude (directly related to seismic energy). The seismic signal was filtered in the band of 1-10 Hz in order to avoid strong ambient noise present in the records. Figure 3 is a temporal record of CO 2 , SO 2 , and H 2 S concentrations measured by the two Multi-GAS instruments.

RESULTS
The RSEM time-series is shown in Figure 3d for comparison, along with Figure 3e the timing of each of the seismicallyidentified phreatic explosions.
The temporal plots exhibit a sequence of gas peaks (numbered #1 to #9 in Figure 3) that clearly emerge above a persistent background (characterized by low concentration of all gases). By comparison with seismic data (Figures 3d,e, 4a), we find that each of the gas peaks corresponds to a Multi-GAS recording of a phreatic explosion, i.e., the peaks are records of the composition of the gas released during a phreatic blast (referred as synexplosive gas). CO 2 and SO 2 (Figures 3a,b, 4a) exhibit the most pronounced peaks in the syn-explosive gas, with peak concentrations of >3,000 ppmv and ∼52 ppmv, respectively. H 2 S peaks are more moderate (typically < 5 ppmv). Nine (Figure 3) of the 42 phreatic events that occurred during February 3-May 9, 2017 are captured by our Multi-GAS record (the remaining explosions took place outside the four daily Multi-GAS acquisition cycles). The background gas mixing ratios thus correspond to continuous passive degassing between the explosions, in what we refer to as the quiescent gas. The quiescent gas is typically characterized by low SO 2 and H 2 S mixing ratios (<2 ppmv), and CO 2 (<500 ppmv) slightly above atmospheric background values.
As in other recent work de Moor et al., 2017), the Multi-GAS-derived gas mixing ratios data are post-processed to calculate the ratios between the various volatiles using the procedure of Tamburello (2015).  Figure 3, see inset) occurred at ∼21:24 GMT. Mixing ratios increase after the blast, starting at ∼21:24 GMT. The yellow-colored area corresponds to the sub-interval (identified with Ratiocalc) where CO 2 , SO 2 , and H 2 S are positively correlated at statistically significant level (R 2 > 0.6). These subsets of CO 2 , SO 2 , and H 2 S mixing ratios are used to draw (b) CO 2 vs. SO 2 and (c) H 2 S vs. SO 2 correlation plots, and to calculate the time-averaged CO 2 /SO 2 and H 2 S/SO 2 ratios (within the yellow-colored temporal windows) from the slopes of the best-fit regression lines (listed in Table 2). In (a) the red and blue lines are high-resolution (1 Hz) records of CO 2 /SO 2 and H 2 S/SO 2 ratios, respectively, calculated by taking the ratio of individual couples of co-acquired mixing ratios. This point-to-point ratios (e.g., Pering et al., 2014) have large (∼ 50 %) associated errors at the ppm level [e.g., before the explosion, (a)], so they should only be viewed as semi-quantitative estimates. Nevertheless, the contrasting gas signatures for pre-explosive and syn-explosive gas are clearly observed. Example of CO 2 vs. SO 2 and H 2 S vs. SO 2 correlation plots for the pre-explosive gas [built from data within the gray-colored temporal windows in (a)], shown in (b') and (c'), confirm a very CO 2 -rich (or SO 2 -poor) gas composition prior to the blast. Gas ratios progressively return to the pre-eruptive levels toward the end of the Multi-GAS acquisition. (d) Example of gas and seismic records taken during a quiescent phase (March 21, 2017 between ∼17:55 and ∼18:30 GMT). Very low gas mixing ratios are observed throughout, as characteristic of the quiescent gas. The yellow-colored area corresponds to the only sub-interval (identified with Ratiocalc) where CO 2 , SO 2 , and H 2 S are positively correlated at statistically significant level (R 2 > 0.6). These subsets of CO 2 , SO 2 , and H 2 S mixing ratios are used to construct (e) CO 2 vs. SO 2 and (f) H 2 S vs. SO 2 correlation plots, from which the time-averaged CO 2 /SO 2 and H 2 S/SO 2 ratios (within the yellow-colored temporal windows) are calculated (listed in Table 2).
To this aim, the gas mixing ratios time-series are sequentially examined with Ratiocalc to identify individual temporal windows (of ≥250 s, corresponding to subsets of ≥ 30 data points) with high correlation coefficients (R 2 > 0.6) between gas couples.
The procedure is illustrated in Figures 4a,d are temporal plots of mixing ratios recorded during a single explosion (Figures 4a-c), and during a typical quiescent degassing phase (Figures 4d-f) (see insets in Figure 4). By sequentially scanning the two sub-datasets, we identify the temporal intervals where the best correlations (R 2 > 0.6) between gas mixing ratio pairs are observed (yellow bands in Figures 4a,d). These subsets of CO 2 , SO 2 , and H 2 S mixing ratios are then used to build CO 2 vs. SO 2 (Figures 4b,e) and H 2 S vs. SO 2 (Figures 4c,f) correlation plots, and to calculate the time-averaged CO 2 /SO 2 and H 2 S/SO 2 ratios (within the yellow-colored temporal windows) from the slopes of the best-fit regression lines. The procedure is repeated for the entire dataset, and all the obtained gas ratio pairs are listed in Table 2.
Our gas ratios, listed in Table 2 and illustrated in Figure 5, show contrasting H 2 S/SO 2 ratios composition for the quiescent gas and the syn-explosive gas. The quiescent gas composition is characterized by H 2 S/SO 2 ratios of 0.57 ± 0.2 (1 σ) (range, 0.01-1.5), or well above the H 2 S/SO 2 of the gas released during the phreatic eruptions. This syn-explosive gas exhibits H 2 S/SO 2 ratios of 0.04 ± 0.06 (range, 0.0003-0.18) (Figure 5). The quiescent gas H 2 S/SO 2 ratio exhibits a weak but appreciable declining trend, from February (∼1) to early May (∼0.5) ( Figure 5B), but remains systematically above the syn-explosive gas range (which shows no systematic trend). Table 2 shows that quiescent gas and syn-explosive gas released by the lake have overlapping CO 2 /SO 2 ratio compositions (64 ± 59 and 136 ± 110 at 1σ; Figure 5A). Both show a tendency of decreasing CO 2 /SO 2 ratios from February (∼120 ± 51) to early May (∼70 ± 20).
Prior to each individual explosion, the CO 2 , SO 2 , and H 2 S mixing ratios are typically very low (see Figures 4a, 6a). Thus, the composition of the quiescent gas released in the minutes/seconds before an explosion, here refereed as pre-explosive gas, is not well constrained. One important aspect is, however, that the compositional change from pre-explosive gas to syn-explosive gas is large and abrupt at the onset of eruption (Figures 4a,  6a). To highlight this observation, we show in Figure 4a the CO 2 /SO 2 and H 2 S/SO 2 ratios calculated by simply taking the ratio of individual couples of co-acquired mixing ratios. This point-to-point ratio approach (e.g., Pering et al., 2014) provides high-resolution (1 Hz) records of CO 2 /SO 2 and H 2 S/SO 2 ratios, see Figures 4a, 6a. The errors associated with this methodology are potentially very large (∼ 50%, as based on laboratory tests) at the few ppm level (e.g., before the explosions), so the ratios displayed in Figures 4a, 6a should only be viewed as semi-quantitative estimates. However, for the specific example shown in Figure 4a, the ratios calculated with the point-to-point ratio technique (respectively of 640 and 1.1 for CO 2 /SO 2 and H 2 S/SO 2 ) are reasonably close to those obtained with the scatter plot methodology described above (respectively of 628 and 0.9 for CO 2 /SO 2 and H 2 S/SO 2 ). A comparable similarity of pointto-point-derived and scatter plot-derived ratios is obtained for the explosion detailed in Figure 6a (see Table 2). We thus argue that the abrupt changes in CO 2 /SO 2 and H 2 S/SO 2 ratios at eruption onsets (Figures 4a,e, 6a) are real, and point to very distinct gas signatures for pre-explosive gas and syn-explosive gas, with the former being typically more SO 2 -poor. Very similar abrupt variations in gas ratios are observed for all 9 recorded explosions, followed by a gradual return to quiescent gas compositions over timescales of minutes (Figures 4a, 6a). The point-to-point-derived compositions of the pre-explosive gas are listed in Table 2.

DISCUSSION
Our Multi-GAS results, in tandem with previous results de Moor et al., 2017), constrain the composition of the gas plume released by the Rincón de la Vieja crater lake (Figure 7). Our 2017 results confirm that SO 2 is the prevalent S gas species in the Rincón de la Vieja crater lake plume during both quiescent and explosive degassing, with H 2 S/SO 2 <1 in all except six cases ( Table 2) (Figure 7). In contrast, low-temperature (<80 • C) fumaroles, scattered on the inner crater walls, and around the lake shore, typically exhibit high H 2 S/SO 2 ratios of ∼4 (Tassi et al., 2005). In these weaker hydrothermal manifestations, S re-equilibration (e.g., SO 2 conversion to H 2 S) during cooling and reaction with wall-rocks in the fumarole's feeding conduits is very likely. Similarly, de Moor et al. (2016b) showed that Poas fumaroles have higher H 2 S/SO 2 than the plume from the acid crater lake.
The 2017 crater lake plume is richer in CO 2 (CO 2 /SO 2 ratios of 64 ± 59 to 136 ± 110 for quiescent gas and synexplosive gas, respectively) than in 2013 (CO 2 /SO 2 = 27 ± 15; Aiuppa et al., 2014) and 2014-2016 (CO 2 /SO 2 of 4.3-9.5; de Moor et al., 2017) (a period over which Rincón's activity was progressively intensifying) (Figure 7). Thus, degassing at Rincón de la Vieja appears to be highly dynamic on yearly time scales in terms of gas compositions. Our continuous record from 2017 has much improved temporal resolution than past (campaign) surveys, and thus allows evaluation of the extent to which gas composition responds to highfrequency changes in volcanic activity. Three key observations emerge: i) the plume H 2 S/SO 2 ratio varied little during quiescent lake degassing in 2017: 0.57 ± 0.20 (Figures 5B, 7); this limited variability is suggestive of S speciation buffering by fluid + solid phase reaction(s) (at least partially) (Giggenbach, 1987(Giggenbach, , 1996; ii) during phreatic blasts, the gas emissions shift to more oxidized conditions (H 2 S/SO 2 <0.18) than during quiescent degassing (Figures 5B, 6); the compositional change at the eruption    The R 2 values (≥ 0.6 in all cases) identify the regression coefficients in the CO 2 vs. SO 2 and H 2 S vs. SO 2 scatter plots that concurred to determine each ratio. SO 2 max identifies the peak SO 2 concentration within each measurement interval. All data refer to quiescent gas except for those associated with an eruption number (which correspond to the syn-explosive gas). The quiescent gas measured prior to each explosion is referred as pre-explosive gas (PRE-EXPL in the table). The quoted compositions for the pre-explosive gas are 20-min averages of point-to-point ratios, e.g., calculated by taking the ratios of individual co-acquired gas concentration couples. Based on laboratory tests, we estimate error in these point-to-point ratios at ∼50%, so the quoted values should only be viewed as semi-quantitative estimates. Only in 2 cases (out of 9 explosions), gas mixing rates were correlated to significant levels (R 2 > 0.6) to allow ratios to be derived with the correlation plot methodology also (see Figures 4a, 6a). The scatter plot-derived gas ratios are listed in parenthesis for comparison in the table.
onset is rapid and abrupt, as indicated by the gas ratio contrast between the pre-explosive gas and explosive gas (Figures 4a, 6a); iii) both the quiescent gas and syn-explosive gas are similarly CO 2rich (relative to S;) compared to magmatic gases in Costa Rica and Nicaragua (see Figure 7 and reference therein). iv) The quiescent gas vented prior to an explosion (the preexplosive gas) has systematically higher CO 2 /SO 2 ratios than the corresponding explosive gas (Figures 4, 6). This difference can only partially reflect the different data processing technique (point-to-point vs. scatter plot, see Figures 4, 6). Thus, also in view of the H 2 S/SO 2 records discussed above (see point ii), we conclude that the explosive gas is especially SO 2 -rich.
These peculiar features of the Rincón de la Vieja crater lake plume have implications for the nature of the magmatic-hydrothermal system producing the emissions, and for the trigger mechanisms of the recurrent explosions, as detailed below.
Rincon de la Vieja Gas Signature: Magmatic or Hydrothermal?
In contrast to the less active (e.g., less acidic) volcanic lakes fed by hydrothermal H 2 S only (e.g., Hasselle et al., 2018), the low (typically <1) H 2 S/SO 2 ratios in the Rincón de la Vieja crater lake gas suggest supply of oxidized (SO 2 -dominated) magmatic fluids into the lake system (Christenson et al., 2010;Christenson and Tassi, 2015). Magmatic gases typically display equilibrium H 2 S/SO 2 ratios buffered by redox conditions in the coexisting FIGURE 5 | Time-series of gas (A) CO 2 /SO 2 and (B) H 2 S/SO 2 (molar) ratios. Quiescent gases, pre-explosive gases and syn-explosive gases are plotted with different colors. Symbol size is proportional to SO 2 mixing ratios. In (A), all pre-explosive gases plot at higher CO 2 /SO 2 ratios (as shown by the arrows).
silicate melt, according to Giggenbach (1987): This magma-inherited H 2 S/SO 2 ratio is then generally preserved during rapid ascent and cooling of magmatic gases, unless hydrothermal storage and re-equilibration occurs (Giggenbach, 1987). Resolving Equation (2) over a range of magmatic temperatures, and at redox conditions (e.g., H 2 /H 2 O ratios) buffered by the silicate melt at QFM (quartz-fayalite-magnetite) and NNO (Nickel-Nickel oxide) buffers (Carmichael and Ghiorso, 1986), a range of equilibrium H 2 S/SO 2 ratios can be obtained (as illustrated in Figure 8), well encompassing the Rincón de la Vieja crater lake gas range. For example, gases separating from magma at ∼900 • C and QFM will have equilibrium H 2 S/SO 2 of ∼1 (constituting the so-called magmatic gas H 2 S/SO 2 buffer of Giggenbach, 1987), similar to gases from the Rincón de la Vieja crater lake (Figures 7, 8).
A continuous supply of SO 2 -rich magmatic gases to the lake is also consistent with the dissolved ion composition of the crater lake water, placing Rincón de la Vieja in the field of high-activity volcanic lakes according to the classification of Pasternack and Varekamp (1997) and Varekamp et al. (2000). Tassi et al. (2005Tassi et al. ( , 2009 argue that the high Cl − /Na + ratio in the Rincón de la Vieja crater lake requires supply of magmatic HCl into the lake, a process also responsible for the high SO 2− 4 and F − concentrations (from the dissolution of magmatic SO 2, H 2 S, and HF). The hyper-acidic chloride-sulfate brine filling the crater lake is also rich in aluminum, iron, zinc, copper, and boron (because of elemental input via both magmatic gases and enhanced rock dissolution), and is thus compositionally similar to the magmatic gas-fed Laguna Caliente Crater Lake, Poás volcano (Tassi et al., 2009).
Whereas, S speciation and lake chemistry clearly support magmatic gas feeding the lake, the high CO 2 /SO 2 signature of the Rincón de la Vieja crater lake gas is inconsistent with a "pure magmatic" hypothesis. Figure 7 compares the Rincón de la Vieja gas compositions (this work and previous studies) with gases emitted by other recently active Costa Rican volcanoes (Poás and Turrialba;de Moor et al., 2016ade Moor et al., ,b, 2017. Both Poás and Turrialba datasets show a spread of volcanic gas compositions, which have been interpreted as reflecting variable extents of mixing between magmatic and hydrothermal end-members (Fischer et al., 2015;de Moor et al., 2016a,b). The magmatic CO 2 /SO 2 endmembers (see stars in Figure 7) have been estimated at ∼0.3-0.5 (Poás) and ∼1-2 (Turrialba), well below the Rincón de la Vieja crater lake gas range for both quiescent and syn-explosive gas types. Similarly, gas observations at nearby Nicaraguan volcanoes Masaya and Momotombo imply a magmatic gas CO 2 /SO 2 ratio of ∼3-6 de Moor et al., 2017), again implying that the Rincón de la Vieja crater lake gas in 2017 has too high CO 2 /SO 2 to be interpreted as pure magmatic (see Figure 7).
The CO 2 -rich composition of Rincón de la Vieja crater lake gases also does not agree with magmatic degassing models (Moretti et al., 2003) using input conditions relevant to Costa Rica-like magma (see de Moor et al., 2016b for details).  Figure 3, see inset) occurred at ∼6.21 GMT. Mixing ratios increase after the blast. The yellow-colored area corresponds to the sub-interval (identified with Ratiocalc) where CO 2 , SO 2 , and H 2 S are positively correlated at statistically significant level (R 2 > 0.6). These subsets of CO 2 , SO 2 , and H 2 S mixing ratios are used to draw (b) CO 2 vs. SO 2 and (c) H 2 S vs. SO 2 correlation plots, and to calculate the time-averaged CO 2 /SO 2 and H 2 S/SO 2 ratios (within the yellow-colored temporal windows) from the slopes of the best-fit regression lines (listed in Table 2). In (a) the red and blue lines are high-resolution (1 Hz) records of CO 2 /SO 2 and H 2 S/SO 2 ratios, respectively, calculated by taking the ratio of individual couples of co-acquired mixing ratios. This point-to-point ratios suggest contrasting gas signatures for pre-explosive and syn-explosive gas. Example of CO 2 vs. SO 2 and H 2 S vs. SO 2 correlation plots for the pre-explosive gas [built from data within the gray-colored temporal windows in (a)], shown in (b') and (c').
These numerical simulations (see Figure 7) predict equilibrium CO 2 /SO 2 ratios of 20 to 0.07 for gases exsolved from a Turrialbalike andesitic magma (at redox condition of QFM +1 to +3 and temperatures of 900-1,100 • C) decompressed from 250 MPa to 0.1 MPa pressure. The measured CO 2 /SO 2 ratios in Rincón de la Vieja crater lake gas (64±59 and 136±110 for the quiescent and syn-explosive gas, respectively) are clearly higher than the typical magmatic range, unless very high gas-melt separation pressures (>> 250 MPa, equivalent to > 9 km depth) are assumed. Such a deep magma source is unlikely, however, because heat pipes typically sustaining intense degassing and convective overturning at crater lakes require shallow magma (Figure 9; Christenson and Tassi, 2015). Also, magma involvement in the recent Rincón de la Vieja eruptions (implicated by the small but ubiquitous juvenile fragments in the eruption deposits) is indicative of shallow storage (Figure 9).
To summarize, the CO 2 -rich gas compositions with low H 2 S/SO 2 ratios (Figure 7) are incompatible with a direct shallow magmatic source for the gas vented by the Rincón de la Vieja crater lake. We propose that the pristine chemical composition of deep magmatic gases is altered during upward migration by gas-water-rock reactions within the (sub)limnic hydrothermal system.

The Role of the Sublimnic Hydrothermal System
The interactions of magmatic gases with either the crater lake or the subjacent hydrothermal system are the most obvious candidates for generating the observed CO 2 -rich (S-depleted) Rincón de la Vieja gas. Magmatic gas-lake and/or magmatic-hydrothermal reactions, if occurring, could also effectively control S partitioning between reduced and oxidized forms, thus modifying the original magmatic H 2 S/SO 2 ratio (Christenson et al., 2010;Christenson and Tassi, 2015).
One often-invoked mechanism of magmatic SO 2 dissolution into active volcanic lakes is via disproportionation reactions such FIGURE 7 | H 2 S/SO 2 vs.CO 2 /SO 2 scatter plot of Costa Rica and Nicaraguan volcanic gases. The Rincón de la Vieja gas (both quiescent gas and syn-explosive gas; same symbols as in Figure 5) is contrasted with the composition of volcanic gases from other volcanoes of the Costa Rica-Nicaragua arc segment. Momotombo and Masaya: de Moor et al., 2017); Poas: de Moor et al., 2017). Estimated compositional fields for hydrothermal (H) and magmatic (M) end-members are indicated (this study; based on averaging of magmatic and hydrothermal compositions quoted in the previous studies). The hypothetical composition of the Rincón hydrothermal pole (H) was estimated by averaging composition of the two measurements of quiescent gases with the lowest SO 2 concentrations (e.g., the most hydrothermal endmember). Numbers alongside the syn-explosive gas data-points identify the explosions (see Figure 3). as (Kusakabe et al., 2000): The two relationships can be combined into the following equilibrium: To assess the role played by magmatic gas-lake water reactions, we test if the measured H 2 S/SO 2 ratios in the Rincón de la Vieja gas (Figures 4-8) are consistent with the equilibrium ratio predicted by Equation (4) at lake water conditions, i.e., if the measured gas S speciation is controlled/buffered by dissolved S species in the crater lake water. In the Rincón de la Vieja crater lake, oxidized (exavalent) dissolved S species prevail (Tassi et al., 2005), and sulfur spherules are always observed on the lake surface (Figure 2), implying that all of the ingredients for reaction 4 are available in abundance. Solving Equation (4) for the H 2 S/SO 2 gas ratio, at crater lake water conditions (T = 30-50 • C; H 2 SO 4 = 8.16 × 10 −12 M at pH = 0.77 and SO 2− 4 = 0.146 M), and using the thermodynamic data from HSC thermochemical software (http://www.hsc-chemistry. com/), we obtain the curve labeled "Lake" in Figure 8. The H 2 S/SO 2 gas ratios predicted by reaction (4) are, however, H 2 S-dominated, pointing against a lake-buffered H 2 S/SO 2 ratio hypothesis (Figure 8). Our calculations thus suggest that either a kinetic process prevails (e.g., that Equation 4 does not go to completion in the lake, leading to only partial conversion of magmatic SO 2 into H 2 S) (Ohmoto and Lasaga, 1982;Rye, 2005), or that magmatic gas-water interactions occur deeper in the system, such as in the hotter sublimnic hydrothermal system (Figure 9).
A variety of magmatic gas-water-rock reactions have been proposed to occur in sublimnic hydrothermal systems (Christenson and Tassi, 2015;Shinohara et al., 2015). These hydrothermal reactions, typically occurring in the upper portion of a heat pipe (Figure 9), invariably enrich the gas phase in CO 2 relative to more reactive SO 2 and H 2 S (Symonds et al., 2001(Symonds et al., , 2003, and may thus generate the high CO 2 /SO 2 ratios in Rincón de la Vieja crater lake gas. In addition to leading to S depletion, these reactions can also act as controls on the H 2 S/SO 2 ratio (see below).
In addition to the reaction described by Equation (4), S scrubbing can occur via (Giggenbach, 1987;Christenson and Tassi, 2015): This reaction scavenges H 2 S and SO 2 in 2:1 proportions, and can thus (at least partially) justify the oxidized nature (SO 2 >H 2 S) of the Rincón de la Vieja gas. In fact, the situation is more complicated, because as pointed out by Shinohara et al. (2011), S deposition via Equation (5) can either increase or decrease FIGURE 8 | Temperature dependence of the SO 2 /H 2 S ratio, as predicted by the common magmatic and hydrothermal redox buffers (modified from Aiuppa et al., 2005). The solid lines represent the H 2 S/SO 2 ratios at equilibrium with redox conditions (oxygen fugacity, fO 2 ) fixed by QFM (Quartz-Fayalite-Magnetite), NNO (Nickel-Nickel oxide), and HM (Hematite-magnetite) redox buffers (Giggenbach, 1987(Giggenbach, , 1996. The magmatic gas range is from Aiuppa et al. (2005). The curve labeled GB and Lake illustrate the equilibrium H 2 S/SO 2 ratios at redox conditions fixed by the Giggenbach (1987) FeO-FeO 1.5 hydrothermal buffer (Equation 6) and dissolved S species in the lake buffer water (Equation 2), respectively. Pressure is 1 bar except where specified. The measured H 2 S/SO 2 ranges for both the quiescent and syn-explosive gas are indicated by gray and red bands, respectively.
the gas H 2 S/SO 2 ratio, depending on the original S speciation (SO 2 -dominated or H 2 S-dominated) in the feeding magmatic gas entering the hydrothermal system (unconstrained for Rincón). Also, concurrently with hydrothermal S deposition, iron minerals in the hydrothermal rock matrix can contribute buffering redox conditions (e.g., the gas H 2 /H 2 O ratio), via (Giggenbach, 1987): These reactions can then buffer, via Equation (2), the residual (after S deposition) hydrothermal gas H 2 S/SO 2 ratio. For example, the H 2 S/SO 2 gas ratios predicted (as a function of temperature) from Equation (2), and redox conditions buffered by either the Fe(II)-Fe(III) hydrothermal buffer (Equation 6; Giggenbach, 1987) or the Hematite-Magnetite buffer (HM) hydrothermal (Equation 7), are graphically illustrated in Figure 8. The figure shows that the measured H 2 S/SO 2 ratios in the Rincón de la Vieja gas would be consistent with those imposed by hydrothermal buffering at HM redox conditions, 290-400 • C, and 1 bar (Figure 8). Coexistence of oxidized (hematite) and reduced (magnetite, pyrite) iron forms is suggested at both Rincón de la Vieja (Tassi et al., 2005) and in the nearby Borinquen and Miravalles hydrothermal fields (Gherardi et al., 2002;Molina and Martí, 2016). However, since a second independent redox couple (e.g., the H 2 /H 2 O ratio; Aiuppa et al., 2011) is not measured in the Rincón gas, evidence for hydrothermal (HM) buffering of S speciation remains speculative.
In summary, available gas information at Rincón de la Vieja is suggestive of extensive hydrothermal processing of magmatic gases within the sublimnic hydrothermal system. These hydrothermal reactions consume S (compared to relatively inert C, thus explaining the high observed CO 2 /S ratio), but the exact S deposition mechanisms (Equations 3-5), and perhaps the role of hydrothermal redox buffering (Equation 7), cannot be quantitatively constrained with the present data. However the observed difference in H 2 S/SO 2 between the quiescent gas (0.57 ± 0.20) and the syn-explosive gas (0.04 ± 0.06), and the abrupt gas ratio changes at eruption onsets (Figures 4a, 6a), clearly imply less hydrothermal interaction (e.g., less SO 2 deposition) during the phreatic blasts (Figure 9). This can be explained by a combination of faster gas transit, transient disruption of the sublimnic hydrothermal envelop, and/or deeper (hotter) gas source during the explosions (Figure 9).

The Eruption Trigger
Whatever the exact process (Equations 3-6), S scrubbing reactions are likely to lead to rapid and effective formation of native S seals (Christenson and Tassi, 2015). The presence of S spherules in the Rincón de la Vieja crater lake points to the existence of a native sulfur layer at the lake bottom (Hurst et al., 1991;Takano et al., 1994;Christenson et al., 2010) (Figure 9).
The physical state of this S pool is known to be dependent on temperature, with a ∼2000-times viscosity increase upon heating in the 150-200 • C temperature range (Hurst et al., 1991;Oppenheimer, 1992;Takano et al., 1994). It is possible that magmatic-hydrothermal gases, while persistently fluxing through the Rincón de la Vieja sublimnic hydrothermal system, progressively heat the S pool, thus ultimately creating the conditions for the development of a viscous impermeable seal ( Figure 9A).
We envision a mechanism in which, during quiescence ( Figure 9A), the feeding magmatic gas interacts with the sublimnic hydrothermal system, undergoes S deposition (± buffering by hydrothermal minerals) and develops a CO 2 -rich gas. A fraction of this gas, perhaps after further S deposition into the lake, is emitted to the surface as quiescent gas, but the remaining fraction is accumulated at depth as the impermeable S seal develops (Figure 8a). The quiescent gas released shortly prior to explosions, e.g., the pre-explosive gas, is characterized by low gas mixing ratios and especially high CO 2 /SO 2 ratios, consistent with a reduction of gas transfer from underneath the seal, and extensive S loss to hydrothermal minerals (Figures 4a, 6a).
Ultimately, gas accumulation underneath the seal leads to pressure buildup and seal failure to trigger a phreatic explosion ( Figure 9B). This mechanism is also supported by the recurrent observation in the erupted products, especially in the 1-2 mm ash grain-size fraction, of native sulfur fragments. We argue that, during the phreatic blasts, rapid gas ascent reduces hydrothermal interactions as supported by the observed short-term variations in gas chemistry. In fact, our high-resolution (1 Hz) of CO 2 /SO 2 and H 2 S/SO 2 records (see Figures 4a, 6a) show that both ratios suddenly decrease at the eruption onset, implying the each syneruptive gas is SO 2 -richer than its corresponding pre-explosive gas. Rapid gas ascent, perhaps combined with gas ascent from FIGURE 9 | A schematic model for the Rincón de la Vieja volcano shallow feeding system, detailing the structure of the sublimnic hydrothermal/magmatic system (modified from Christenson et al., 2010). (A) During quiescence, the feeding magmatic gas undergoes S deposition (± buffering by hydrothermal minerals) within the sublimnic hydrothermal system. The resulting CO 2 -rich gas is partially surface-vented as quiescent gas, and partially accumulated beneath an impermeable S seal; (B) pressure buildup below the S seal triggers a phreatic explosion, characterized by rapid gas ascent, limited hydrothermal interaction, and deeper gas source. This syn-explosive gas retains a C-rich signature, implying a magmatic-hydrothermal nature.
deeper/hotter portions of the magmatic-hydrothermal vapor zone (Figure 9B), are likely implicated in producing the relatively SO 2 -richer syn-explosive gas.
We caution that our interpretation is based on only ∼3 months of observations and the capturing of only 9 out of 42 explosions in that time interval. Thus, longer observations are required to draw more concrete conclusions. However, our results confirm that continuous instrumental geochemical monitoring, in tandem with seismic monitoring, can contribute to understanding the mechanisms that drive the Rincón de la Vieja magmatic system toward potentially hazardous critical states. It is interesting to note that both the H 2 S/SO 2 and CO 2 /SO 2 ratios in the Rincón de la Vieja quiescent gas decreased from February to early May (Figures 5, 7), implying a decreasing extent of hydrothermal re-equilibration of the feeding magmatic gases. This trend toward more magmatic-like gas composition may tentatively be interpreted as reflecting an increase in the magmatic gas supply to the lake conduit, and thus a more unstable magmatic system leading to the major lahar-generating phreatic eruption of May 23 (Global Volcanism Program, 2017), which occurred just a few weeks after our Multi-GAS dataset ended (and actually destroyed the instrument). As stated above, this event (and the following June 11 event) erupted a sizeable juvenile fragment component, claiming for escalating magmatic activity.

Implications for Monitoring of Active Crater Lakes
Our measurements of the syn-explosive gas, the first of their nature at an active volcano, provide observational evidence for that formation of hydrothermal S seals may strongly be implicated in the generation of phreatic eruptions, as originally proposed at Ruapehu volcano in New Zeeland (Hurst et al., 1991;Takano et al., 1994;Christenson et al., 2010). It is arguable that hydrothermal seal formation may fuel phreatic/phreatomagmatic activity at other recurrently erupting volcanic lakes worldwide (Christenson and Tassi, 2015;Stix and de Moor, 2018). As our Rincón de la Vieja dataset shows, continuous gas records may help tracking such seal formation events, as marked by phases of reduced surface S degassing (see the low SO 2 and H 2 S mixing ratios, and high CO 2 /SO 2 ratios, in the pre-explosive gas; Figures 4a, 6a, 7). To test this hypothesis further, our results thus claim for the need of augmenting the number of volcanic crater lakes that are monitored by permanent gas instrumentation.
It is also interesting to note that the rapid and abrupt transition between pre-explosive and explosive gas at Rincón de la Vieja (Figures 4a, 6a), pointing to reduced surface S outgassing prior to each blast, is in stark contrast with what recurrently reported at Poás, where precursory trends toward more magmatic (more S-rich) gas composition are observed instead (de Moor et al., 2016a). We argue this difference may imply an overall lower magmatic gas input, or a more effective seal (or both), at Rincón de la Vieja, compared to the more active Poás volcano (where the magmatic system appears more implicated as eruption trigger). This comparison suggests that distinct trigger mechanisms, with different relative roles played by hydrothermal sealing vs. magmatic as influx, may in fact be operating at active crater lakes. This diversity in potential trigger mechanisms (Stix and de Moor, 2018) reinforces the need for more robust (and temporally continuous) gas records at volcanic crater lakes.

CONCLUSIONS
This study investigates degassing dynamics of Rincón de la Vieja during a period of intense phreatic activity. For the first time at any volcano, the composition of the gas released during discrete phreatic events (confirmed seismically) was resolved using Multi-GAS. Our results demonstrate chemically distinct gas compositions during quiescent degassing vs. explosive eruptive degassing. The quiescent gas is characterized by very low concentrations of sulfur gas species (SO 2 and H 2 S < 2 ppmv) and relatively high H 2 S/SO 2 ratios (mean 0.57), whereas the synexplosive gas shows much higher gas mixing ratios of SO 2 and CO 2 and a systematic decrease of H 2 S/SO 2 ratio (mean, 0.04). Both the quiescent and syn-explosive gases exhibit high CO 2 /S ratios relative to regional magmatic gases, but each explosive gas has lower CO 2 /S ratio than its corresponding pre-explosive gas.
These C-rich gases, if interpreted as magmatic, would require a very deep (P > 250 MPa) magmatic source, which is unlikely in view of the recent magma involvement in the 2017 Rincón phreatomagmatic eruptions. We thus favor a mechanism in which magmatic gas-water-rock reactions in the sublimnic hydrothermal system lead to deposition of native S and sulfate formation, and thus enrich the gas phase in C relative to SO 2 and H 2 S. The presence of abundant native sulfur spherules in the lake is consistent with this scenario, as the presence of S spherules suggests a native sulfur layer at the lake bottom. We propose that continuous gas fluxing may progressively lead to development of an impermeable seal near the lake/hydrothermal system interface. Gas accumulates underneath the seal, and ultimately ruptures this seal resulting in phreatic eruptions. The trend of decreasing quiescent CO 2 /SO 2 and H 2 S/SO 2 ratios from February to May 2017 can be interpreted as an increase in the magmatic gas supply, providing a precursor to the major eruption of 23 May.
These results, the first volcano reports for the gas phase released by discrete phreatic eruptions, confirm (Christenson and Tassi, 2015) that the complex interplay between rising magmatic gases and the sublimnic hydrothermal system likely plays a decisive role in triggering the explosions.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

AUTHOR CONTRIBUTIONS
AB, JMdM, and AA conceived the idea. AB, JMdM, GA, HB, HV, MB, GG, and PK prepared the instrumentation and conducted the field experiment. AB, MM, HB, and DD processed the data. AB and AA drafted the article with help from JMdM, PK, GA, HB, and MM.