# Magmatic gas percolation through the old lava dome of El Misti volcano

Yves Moussallam, Nial Peters, Pablo Masias, Fredy Apaza, Talfan Barnie, C.

Ian Schipper, Aaron Curtis, Giancarlo Tamburello, Alessandro Aiuppa,
Philipson Bani, et al.

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

To cite this version:

Yves Moussallam, Nial Peters, Pablo Masias, Fredy Apaza, Talfan Barnie, et al. Magmatic gas percolation through the old lava dome of El Misti volcano. Bulletin of Volcanology, 2017, 79, 10.1007/s00445-017-1129-5. insu-03708988

HAL Id: insu-03708988

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

Submitted on 30 Jun 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Magmatic gas percolation through the old lava dome of El Misti volcano

Yves Moussallam<fnref n="1" /> · Nial Peters<fnref n="1" /> · Pablo Masias<fnref n="2" /> · Fredy Apaza<fnref n="2" /> · Talfan Barnie<fnref n="3" /> ·

C. Ian Schipper<fnref n="4" /> · Aaron Curtis<fnref n="5" /> · Giancarlo Tamburello<fnref n="6" /> · Alessandro Aiuppa<fnref n="7" />,<fnref n="8" /> ·

Philipson Bani<fnref n="9" /> · Gaetano Giudice<fnref n="8" /> · David Pieri<fnref n="5" /> · Ashley Gerard Davies<fnref n="5" /> ·

Clive Oppenheimer<fnref n="1" />

Received: 2 February 2017 / Accepted: 16 May 2017 / Published online: 24 May 2017
© The Author(s) 2017. This article is an open access publication

**Abstract** The proximity of the major city of Arequipa to El Misti has focused attention on the hazards posed by the active volcano. Since its last major eruption in the fifteenth century, El Misti has experienced a series of modest phreatic eruptions and fluctuating fumarolic activity. Here, we present the first measurements of the compositions of gas emitted from the lava dome in the summit crater. The gas composition is found to be fairly dry with a H₂O/SO₂ molar ratio of 32 ± 3, a CO₂/ SO₂ molar ratio of 2.7 ± 0.2, a H₂S/SO₂ molar ratio of 0.23 ± 0.02 and a H₂/SO₂ molar ratio of 0.012 ± 0.002. This magmatic gas signature with minimal evidence of hydrothermal or wall rock interaction points to a shallow magma source that is efficiently outgassing through a permeable conduit and

lava dome. Field and satellite observations show no evolution of the lava dome over the last decade, indicating sustained outgassing through an established fracture network. This stability could be disrupted if dome permeability were to be reduced by annealing or occlusion of outgassing pathways. Continued monitoring of gas composition and flux at El Misti will be essential to determine the evolution of hazard potential at this dangerous volcano.

**Keywords** Volcanic hazard · Arequipa · Outgassing · ASTER · Multi-GAS · Trail by fire

**Highlights**
* First characterization of plume composition at El Misti volcano

* Gas chemistry reveals a clear magmatic signature

* Field and satellite observations suggest a stable outgassing pathway through the conduit and old lava dome

Editorial responsibility: P. Allard

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

Yves Moussallam
ym286@cam.ac.uk

<footnote n="1">¹ Department of Geography, University of Cambridge, Downing Place, Cambridge CB2 3EN, UK</footnote>

<footnote n="2">² Observatorio Vulcanológico del Ingemmet (OVI), Arequipa, Peru</footnote>

<footnote n="3">³ Nordic Volcanological Center, Institute of Earth Sciences, Sturlugata 7 – Askja, 101, Reykjavik, Iceland</footnote>

<footnote n="4">⁴ School of Geography, Environment and Earth Sciences, Victoria University of Wellington, PO Box 600, Wellington 6140, New Zealand</footnote>

<footnote n="5">⁵ Jet Propulsion Laboratory-California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA</footnote>

<footnote n="6">⁶ Istituto nazionale di geofisica e vulcanologia, sezione di Bologna, Bologna, Italy</footnote>

<footnote n="7">⁷ Dipartimento DiSTeM, Università di Palermo, Via archirafi 36, 90146 Palermo, Italy</footnote>

<footnote n="8">⁸ Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo Via La Malfa, 153, 90146 Palermo, Italy</footnote>

<footnote n="9">⁹ Université Clermont Auvergne, CNRS, IRD, OPGC, Laboratoire Magmas et Volcans, F-63000 Clermont-Ferrand, France</footnote>
## Introduction

Lava domes are the extruded part of a magma risen through a conduit. Typically associated with silicic (rhyolitic to andesitic) magmas, lava domes are often unstable and capable of generating pyroclastic density currents. Dome collapse events have been documented at volcanoes such as Soufrière Hills (Montserrat; e.g. Watts et al., 2002), Mount Unzen (Japan; e.g. Sato et al., 1992), Mount St. Helens (USA; e.g. Mellors et al., 1988) and Merapi (Indonesia; e.g. Komorowski et al., 2013). Critical in determining whether or not an otherwise gravitationally stable dome will become overpressured and collapse is the nature of magmatic outgassing through either the conduit wall or fracture networks (e.g. Jaupart and Allègre, 1991; Gonnermann and Manga, 2003; Boudon et al., 2015). While such considerations have been explored theoretically (e.g. Sparks, 1997; Melnik and Sparks, 1999, 2005; Hale and Mühlhaus, 2007) the inherent hazard associated with working on active lava domes has limited in situ collection of data on gas compositions, limiting data to that from remote sensing observations (e.g. Oppenheimer et al., 2002; Edmonds et al., 2003; Holland et al., 2011) with only a few direct measurements (e.g. Soufrière Hills; Hammouya et al., 1998).

El Misti volcano, southern Peru (Fig. 1), is close to Arequipa, the second largest city in the country, with about one million inhabitants. El Misti is a composite stratovolcano composed of four stratocones (Thouret et al., 2001). The youngest, Misti 4, forms the current summit of the volcano and has erupted at least ten times since 11,000 years. B.P producing pyroclastic-surges and lahars that have travelled up to 13 km from the vent (Thouret et al., 2001). The last plinian eruption occurred ca 2050 years ago and produced a widely dispersed pumice-fall deposit extending ≥25 km from the vent (Thouret et al., 1995, 2001; Harpel et al., 2011; Cobeñas et al., 2012). Since then, eruptive activity has been mild, with minor events in 655–865 A.D., 1304–1398 A.D. (Thouret et al., 2001) and in 1440–1470 A.D. (Murúa, 1946, 1987); episodes of increased fumarolic or seismic activity, and small phreatic eruptions (Barriga, 1951; Zamácola and Jaureguí, 1958; Hantke and Parodi, 1966; Chávez Chávez, 1992; Simkin and Siebert, 1994; Thouret et al., 2001). Within the

crater of El Misti is an old lava dome of unknown age, potentially as old as the last magmatic eruption (fifteenth century), on which a fumarole field is located (Fig. 2). The dome is covered with sulphur sublimate previously characterised by Birnie and Hall, (1974). The current fumarolic activity has persisted since at least 1787 (Thouret et al., 2001), and a persistent thermal anomaly of ~+6 K has been identified at the summit in Advanced Spaceboune Thermal Emission and Reflection Radiometer (ASTER; Yamaguchi et al., 1998) thermal infrared images from 2000 to 2010 (Jay et al., 2013). ASTER is on-board NASA's *Terra* spacecraft.

Despite a number of studies emphasising the hazards and risks of future volcanic activity of El Misti (e.g. Sandri et al., 2014), very little work has focused on its ongoing outgassing. Here, we present the first compositional data for the gas plume emitted from the El Misti lava dome together with an analysis of photographic, thermal camera and satellite images. We evaluate the results in terms of outgassing behaviour and the implications for hazard evolution.

## Methods

### In situ gas measurements

Gas composition data were obtained using a portable Multi-GAS instrument (Shinohara, 2005) deployed directly inside the crater and on the dome of El Misti (S 16° 17' 57.22"; W 71° 24' 20.69"; 5600 m a.s.l.). The instrument incorporated SO₂, H₂S and H₂ electrochemical sensors. The SO₂ and H₂ sensors have calibration range of 0–200 ppmv while the H₂S sensor had a calibration range of 0–100 ppmv. An NDIR sensor was used for CO₂ and calibrated for 0–10,000 ppmv with an accuracy ±2%. A relative humidity (R.H.) sensor (Galltec) was used to measure H₂O, providing a measuring range of 0–100% R.H. with an accuracy of ±2%.

The conversion from relative humidity to water mixing ratio was made following Buck (1981) and using the following equation:

$$H_2O = \frac{6.1121 \times (1.0007 + 3.46 \times P^{-6}) \times \exp\left[\frac{17.502 \times T}{240.97 + T}\right] \times \frac{Rh}{100} \times 10^6}{P} \quad (1)$$

where H₂O is the absolute water concentration in parts per million by volume, T is the temperature in degrees Celsius, Rh is the relative humidity in percent and P is the atmospheric pressure in millibars. The gas temperature used in this equation is measured in real time by the Multi-GAS, the pressure is also measured by the Multi-GAS and assumed to remain constant during the

measurements. All sensors were housed inside a weatherproof box, with the ambient air sampled via Teflon tubing connected to a HEPA filter fed through an inlet in the box and circulated via a miniature 12-V rotary pump through the sensors. An on-board data-logger captured measurements at a rate of 1 Hz. The complete system is powered by a small (6 Ah) 12 V LiPo battery.
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Location map showing the location of El Misti volcanoes together with the location of all Peruvian Holocene volcanoes and the Nazca subduction trench. The centre of Arequipa (one million inhabitants) is located 18 km from El Misti's summit</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Visible (a) and infrared (b) images of the lava dome of El Misti taken on 1 December 2015 from within the crater, looking south. The dome is about 150 m across</figcaption></figure>

Similar systems have now been successfully deployed at many volcanoes (e.g. Aiuppa et al. 2011, 2012, 2014, 2015, Moussallam et al. 2012, 2014, 2016). All sensors were calibrated in the laboratory at INGV Palermo (October 2015), with target gases of known concentration. The differences in response time for the different sensors were corrected by finding the lag times from correlation analysis of the various time series. Ambient air composition was subtracted from the CO₂, H₂O and H₂ data, SO₂ interference on H₂S data was calibrated and corrected. Multi-GAS measurements were taken on 1 December 2015 for 90 min. Post processing was performed using the Ratiocalc software (Tamburello, 2015).

## Infrared and visible camera

Thermal infrared images were acquired from within the crater rim pointing at the dome. Images were taken using an OPTRIS PI400 camera (spectral range of 7.5–13 µm with an optical resolution of 382 × 288 pixels), fitted with an 8-mm lens providing a FOV of 62°×48°. Thermal images were acquired on 1 December 2015 using the temperature range setting of –20 to 100 °C.

Photographs of the lava dome were acquired from the crater rim by OVI personnel during repeated ascents in the period from 2007 to 2016.

## Satellite observations

We used land surface temperature maps derived from ASTER thermal images (AST08 Land Surface Temperature, or LST,
product at 90 m/pixel spatial resolution), short wavelength infrared Advanced Land Imager (ALI) images (nine bands from 0.4 to 2.4 µm at 30 m/pixel spatial resolution, with a panchromatic band (PAN) at 10 m/pixel spatial resolution—Ungar et al. 2003), Hyperion hyperspectral data (196 usable bands from 0.4 to 2.4 µm at 30 m/pixel spatial resolution) (Pearlman et al. 2003) and high resolution visible images (from Google Earth) to track potential changes in activity at the lava dome in the 15 years leading up to the fieldwork presented here. Hyperion and ALI are on-board NASA’s *Earth Observing 1* (EO-1) spacecraft. Observations by EO-1 were obtained via the NASA Volcano Sensor Web (VSW, e.g. Davies et al., 2015), which utilizes advanced spacecraft operations software developed to streamline the process to re-task spacecraft as quickly and efficiently as possible to effect rapid data acquisition of dynamic targets (Chien et al., 2005).

All ASTER granules covering El Misti were acquired from NASA Reverb at both level 1b (radiometrically and geometricaly corrected at sensor radiance) and at level 3 AST08, extending the coverage of Jay et al., (2013) from 2010 to late 2016. The AST08 LST dataset is produced by applying the temperature emisivity separation algorithm (TES) to atmospherically corrected images (Gillespie et al., 1998). We used DA White’s Aster Preprocessing Toolkit (APTK, White 2016) to extract the correct rational polynomial coefficients (RPCs) from the level 1b data, which were then used to orthorectify the LST product to the 30-m shuttle radar topography mission digital elevation model with bilinear resampling. All images were classified by inspection as either cloudy or cloud free, then pixels in the summit region were extracted from the cloud free images. Previous studies have manually selected thermal anomalies in ASTER images of Andean volcanoes (Jay et al. 2013); however in this study, we opted for an automated approach for expediency. Blank zero padding pixels that fell within the summit region were removed using a threshold and a simple measure of the thermal anomaly, $T_{anomaly}$, was calculated:

$$T_{anomaly} = T_{max} - T_{50\%}$$

Where $T_{max}$ is the maximum temperature in the summit area, and $T_{50\%}$ is the 50th percentile (or median) temperature. This gives a measure of the maximum temperature relative to a robust estimate of the surrounding temperature. We limited our study to night time images, when the contrast between cool ground and geothermally heated areas is highest. The LST maps have a resolution of 90 m and allowed us to track changes in temperature over a broad area of the summit of El Misti.

Time series images of the lava dome in the visible where taken from Google Earth from the period of 2002 to 2016 using images from Digital Globe, NASA, Landsat/Copernicus and CNES/Astrium.

Images from the ALI and Hyperion on EO-1 were acquired around the fieldwork period on 11th October and 4th

December 2015. Additional ALI images from 2002, 2003 and 2014 were also downloaded from the USGS Earth Explorer archive. All ALI images were processed to level 1T (precision orthorectified product). The short wavelength infrared bands of ALI have a resolution of 30 m, and allow small but very hot features to be identified, however all images were acquired during daytime so small thermal anomalies might be lost amid reflected sunlight. The VSW processes all Hyperion data searching for thermal anomalies (see Davies et al., 2006) using software originally developed to do this thermal emission detection on-board the spacecraft (the Autonomous Sciencecraft Experiment (ASE)—Chien et al., 2005; Davies et al., 2006), and has proved capable of detecting small thermal anomalies even in daylight. The Hyperion pixel brightness temperature detection limits of the ASE thermal classifier software are 426 K at 2.28 µm and 530 K at 1.65 µm (Davies et al., 2006).

## Results

### Gas composition

We obtained 30 min of very high quality measurements presented in Fig. 3, which shows four scatter plots of SO₂ vs CO₂; H₂O; H₂S and H₂ mixing ratios in the El Misti plume. H₂O and CO₂ concentrations and mixing ratios in the volcanic gas are shown after subtraction of their respective mean concentrations in ambient air (measured by the Multi-GAS in the crater but prior to entering the plume). H₂S mixing ratios are shown after correction for laboratory-determined interference with SO₂ gas (16%). The strong positive covariations observed between SO₂ and the other detected volatiles confirm a single, common, volcanic origin. The gas/SO₂ molar plume ratios were obtained by calculating the gradients of the best-fit regression lines. Scatter plots yield CO₂/SO₂ molar ratios of 2.7 ± 0.2, a H₂O/SO₂ molar ratios of 32 ± 3.4, a H₂S/SO₂ molar ratios of 0.23 ± 0.02 and a H₂/SO₂ molar ratios of 0.012 ± 0.002. Together these data yield molar proportions of H₂O, CO₂, SO₂, H₂S and H₂ gases of 89, 7.5, 2.8, 0.6 and 0.03 mol% (Table 1). The correlation between H₂ and SO₂ mixing ratios is much weaker than for other species ($R^2$ value of 0.21 compared to $R^2$ values of 0.93, 0.88 and 0.94 for CO₂, H₂O and H₂S vs SO₂). The two red lines in Fig. 3 show a conservative estimate of the range of H₂ to SO₂ molar ratio from 0.05 to 0.003 that could be derived from the data.

### Vent temperature

Parodi (1966) estimated the fumarole gas temperatures at over 250 °C (unknown date) stating that sulphur appeared melted in the sources’ emission gaps. Birnie and Hall (1974) estimated fumaroles temperatures of 100–125 °C while Thouret et al. (2001) indicate a maximum temperature of 220 °C measured
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 3** CO₂, H₂O, H₂S and H₂ vs SO₂ scatter plots of the mixing ratios in the El Misti plume. Measurements were taken on 1 December 2015 for 30 min at an acquisition frequency of 1 Hz. Least-square regression lines</figcaption></figure>

in December 1997. In 2013 OVI personnel measured fumaroles temperatures between 270 and 310 °C. In 2015, our thermal camera measurements indicate a temperature in excess of 125 °C (temperature at which the image saturated). Given the distance at which the image was acquired and the small size of the vents, fumarole temperatures in excess of 200 °C are expected. Inspection of short wave infrared images from the ALI on EO-1 acquired in 2002, 2003, 2014 and 2015 reveals the presence of a small plume above the dome, but no thermal anomaly. This observation, together with the absence of incandescence observed on the dome at night, suggests vent

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

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

are shown in dotted blue on each plot. Red lines on the H₂ vs SO₂ scatter plot show a conservative estimate of the range of potential gas ratio

temperatures below around 600 °C (basaltic rock is known to glow red at ~700 °C; Decker and Christiansen, 1984).

## Time series observations

Direct images of the dome taken from the crater rim from the period 2007 to 2016 are shown in Fig. 4. To our knowledge, the oldest published photograph of the dome dates from 1967 and is presented in Birnie and Hall, (1974). During the investigated period and since 1967, no changes in the dome morphology are apparent. The location of the fumarole field does

<table><caption>**Table 1** X/SO₂ molar and mass ratios measured by multi-gas and gas composition of the mixed plume at El Misti volcano. Error are expressed as the standard error of the regression analysis and subsequent error propagation</caption><thead><tr><th>Volcano</th><th>Gas</th><th>Mixed plume<br>molar ratio<br>(X/SO<sub>2</sub>)</th><th>Error<br>(1σ)</th><th>Mixed plume<br>mass ratio<br>(X/SO<sub>2</sub>)</th><th>Error<br>(1σ)</th><th>Mixed plume<br>composition<br>(mol%)</th><th>Error<br>(1σ)</th></tr></thead><tbody><tr><td rowspan="5">El Misti</td><td>H<sub>2</sub>O</td><td>32</td><td>3</td><td>9</td><td>1</td><td>89</td><td>9</td></tr><tr><td>CO<sub>2</sub></td><td>2.7</td><td>0.2</td><td>1.9</td><td>0.1</td><td>7.5</td><td>0.6</td></tr><tr><td>SO<sub>2</sub></td><td>1</td><td>0</td><td>1</td><td>1</td><td>2.8</td><td>0</td></tr><tr><td>H<sub>2</sub>S</td><td>0.23</td><td>0.02</td><td>0.12</td><td>0.01</td><td>0.64</td><td>0.06</td></tr><tr><td>H<sub>2</sub></td><td>0.012</td><td>0.002</td><td>0.00038</td><td>0.00005</td><td>0.033</td><td>0.005</td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Fig. 4** Still photograph of the El Misti lava dome taken from the crater rim during repeated ascent from 2007 to 2016 by OVI personnel. The pictures from 2008 were taken by Victor Aguilar (Universidad Nacional de San Agustin). Note the lack of changes in both the dome morphology and location of fumaroles</figcaption></figure>

not change, and while variations in the intensity of the outgassing can be seen, these may simply indicate variability

in atmospheric conditions (e.g. variations in temperature- and humidity-dependent steam condensation) rather than true
variability in gas flux. Periods of increased gas flux have, however, been documented at El Misti, such as in 2011 when a constant gas plume could be seen from Arequipa (OVI observations). Satellite observations at visible wavelengths spanning the period 2002 to 2016 are shown in Fig. 5. Consistent with field observations, these also show no variations in the dome morphology, area or situation of the fumarole field. ASTER observations of land surface temperature covering the period 2002 to 2016 are shown in Fig. 6. Compared with the time series from neighbouring Sabancaya volcano, which has shown varying levels of activity during the same period, the heat output from El Misti's lava dome has stayed relatively constant, consistent with a 'steady state' outgassing and heat output. Inspection of short wave infrared images from the ALI on EO-1 acquired in 2002, 2003, 2014 and 2015 reveal the presence of a small plume above the dome, but no detectable thermal anomaly. Neither was any anomalous thermal emission detected in Hyperion images.

## Discussion

### A magmatic gas signature

The composition of the gases emitted from the El Misti lava dome is indicative of derivation directly from outgassing magma. This is apparent in the high SO₂ concentration and the relatively low H₂O/SO₂ ratio, which together imply small to negligible contributions from a hydrothermal system. Given that we measured both H₂O-H₂ and SO₂-H₂S redox couples, the gas-melt equilibrium temperature and oxygen fugacity can be calculated. Following Giggenbach (1987, 1996) and using the thermodynamic data of Stull et al. (1969) provides two equations with two unknowns:

$$ \log \frac{H_2}{H_2O} = -\frac{12707}{T} + 2.548 - \frac{1}{2} \log f O_2 \quad (2) $$

and

$$ \log \frac{SO_2}{H_2S} = \frac{27377}{T} - 3.986 + \frac{3}{2} \log f O_2 - \log f H_2O \quad (3) $$

Solution yields an equilibrium temperature of 532 °C and a logfO₂ equivalent to ΔQFM = +2.8 (where QFM refers to the quartz-fayalite-magnetite buffer, and where ΔQFM = logfO₂-logfO₂ of QFM at corresponding temperature) or ΔNNO = +2.0 (where NNO refers to the nickel-nickel oxide buffer, and where ΔNNO = logfO₂-logfO₂ of NNO at corresponding temperature). Error on the measured gas ratios, especially on the H₂/SO₂ ratio results in non-symmetrical uncertainty of -96 and +133 °C on the equilibrium temperature and of -1.2 and +1.2 log units on the deviation from the QFM or NNO buffer. Equations for the QFM and NNO buffer used

here are from Frost, (1991). The value of fH₂O used here is 0.88 given that at 1 bar the fugacity of a gas is equal to its partial pressure and that P(H₂O) = (Ptot × nH₂O)/ntot = [(1 bar)(0.88ntot)]/(ntot) = 0.88 bar.

The equilibrium temperature of 532 °C is higher than the temperature at which scrubbing of magmatic gases by hydrothermal systems is expected to become significant (Symonds et al., 2001; Gerlach et al., 2008) giving confidence that the reported gas composition has not been affected by secondary processes other than cooling. The high oxidation state preserved by the gases further indicates limited interaction with low temperature rock (Giggenbach, 1987). The absence of contamination of the dome gas composition by hydrothermal fluids is consistent with the idea of Finizola et al., (2004) who suggested that the hydrothermal system at El Misti is sealed by hydrothermal alteration. The clear magmatic signature and high equilibrium temperature of the emitted gases, together with their high exit temperature (270–310 °C), suggest a relatively shallow magmatic source (see model from Stevenson, 1993). While the current composition of the magmatic source is unknown it may be similar to either the rhyolitic or andesitic magmas that interacted during the ca 2050 BP eruption (Tepley et al., 2013). Tepley et al. (2013) estimated the temperature of the andesitic magma at 940 ± 40 °C using pyroxenes pairs thermometry and the temperature of the rhyolitic magma at 816 ± 30 °C using groundmass ilmenite and magnetite thermometry. The temperature recorded by the gas composition (Eqs. 2–3) represent the temperature at which the gases were last in equilibrium and hence falls between the temperature of the vents and that of the magma.

### A permeable conduit and stable outgassing pathway

Both direct (Figs. 2 and 4) and satellite (Figs. 5 and 6) observations show a very stable dome structure with little to no variations in terms of geometry, distribution of fumaroles, or heat output over at least the last 15 years. We note that previous InSAR surveys at El Misti also found no deformation of the edifice between 1992 and 2002 (Pritchard and Simons, 2004) nor between 2006 and 2009 (Gonzales, 2009). Together these observations imply a stable structure with established percolation pathways for the gas and little to no build-up of pressure within the edifice (Fig. 7). Figure 2 shows a strong correspondence between the location of thermal hotspots on the dome and the location of gas discharge.

While the exact date at which the lava dome formed within the crater is unknown, reports of fumarolic activity date back to the last magmatic eruption in 1440–1470 A.D. (Murúa, 1946, 1987). Several periods of increased fumarolic activity have since been reported (Thouret et al., 2001 and references therein) but no major eruptions have occurred, and it can hence be assumed that the current lava dome dates from the fifteenth century. The high equilibrium temperature we
<figure><img src="image_11.png" /><figcaption>**Fig. 5** Satellite images of the El Misti lava dome from 2002 to 2016 taken from Google Earth and using images from Digital Globe, NASA, Landsat/Copernicus and CNES/Astrium. Note the lack of changes in both the dome morphology and location of the fumaroles field. Up is north on all images. The lava dome at the centre of each image is approximately 200 m in diameter</figcaption></figure>

derived from the gas composition may partly explain the long-term preservation of permeable pathway through the conduit as the gas flow maintains a temperature too high for extensive

precipitation of solids from the gas phase. Another possibility is that the relatively dry gas composition (unaffected by the hydrothermal system) is not prone to mineral precipitation

<figure><img src="image_12.png" /><figcaption>**Fig. 6** Time series of temperature anomalies (maximum temperature minus median temperature for all pixels in the summit area) derived from the ASTER AST08 land surface temperature product for El Misti and Sabancaya</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. 7** Schematic cross-section through the El Misti conduit highlighting the main conclusions from this study. Magmatic gases are released from a reservoir at unknown depth and quickly migrating to the surface through a network of established fracture with the conduit and lava dome. During ascent, the gas has very limited chemical interaction with the host rock and remains isolated from contamination by the surrounding hydrothermal system</figcaption></figure>

before extensive cooling. The flux of gas coming out of the dome is fairly low (see discussion below) and it would hence be surprising if the gas pressure was responsible for maintaining a permeable fracture network.

## Implications for hazards and volcanic monitoring

The current ‘stable’ activity suggests an opportunity in terms of monitoring and hazard assessment in the sense that the gas composition, being largely unaffected by hydrothermal and scrubbing processes, would quickly respond to any changes happening in the magmatic system or conduit. For instance, a change of the temperature or composition of the magma following a recharge event—such as inferred for the last Plinian eruption in ca 2050 BC (Tepley et al., 2013)—should produce a measurable change in the equilibrium temperature and oxidation state of the gases. A change in the established outgassing pathways by fracture healing (e.g. Heap et al., 2015) or pore network collapse (e.g. Kennedy et al., 2016), for instance, could result in a decrease of the total gas flux. Visual observations (Fig. 4) point toward a very low SO₂ flux at present, probably <50 t/day, rendering measurements by UV-based SO₂ flux measurements challenging but maybe not impossible.

The report of phreatic eruptions since the fifteenth century (11 possible events, occurring in 1542?, 1599?, 2 May 1677, 9

July 1784, 28 July and 10 October 1787, August 1826?, August 1830 and 1831?, September 1869, and March 1870? see Parodi, 1966; Simkin and Siebert, 1994 and references therein) begs the question of their origin. The measured gas composition indicates that the current outgassing is isolated from the surrounding hydrothermal system. One hypothesis would be that following periods of increased precipitation, groundwater could interact with the hot magmatic gases generating small eruptions. Confirmed events do not however seem to occur preferentially in the rainy season. In the case of El Misti another hypothesis is that most recorded phreatic eruptions rather refer to periods of increased magmatic gas release either following a transient sealing then opening of the outgassing pathway or an increased flux of magmatic gas from depth. Distinguishing the relative roles of groundwater and magmatic fluids based on these historical observations recorded from a distance is not possible, and neither hypothesis can be dismissed. Both scenarios however would follow a perturbation of the system that has not been seen over our observational period.

## Conclusion

We measured the composition of gases emitted by El Misti in December 2015 and examined the evolution of the lava dome using satellite and direct observations dating back to 2002. The gas composition is indicative of magmatic outgassing with negligible contamination from a hydrothermal system. Together with the apparent stability of the lava dome and fumarolic field through time this implies the efficient release of magmatic gases from the reservoir to the surface through an established fracture network. Future gas monitoring campaigns will be worthwhile in order to track any potential evolution of the magmatic system and conduit.

**Acknowledgements** This research was conducted as part of the ‘Trail By Fire’ expedition (PI: Y. Moussallam). The project was supported by the Royal Geographical Society (with the Institute of British Geographers) with the Land Rover Bursary; the Deep Carbon Observatory DECADE Initiative; Santander, Ocean Optics; Crowcon; Air Liquide; Thermo Fisher Scientific; Cactus Outdoor; Turbo Ace and Team Black Sheep. We thank Jean-loup Guyot, Sebastien Carretier, Rose-Marie Ojeda, Pablo Samaniego and Jean-Luc Lepennec together with IRD South-America personnel for all their logistical help. We are extremely grateful to Marco Rivera and all OVI personnel for their help and support. YM acknowledges support from the Scripps Institution of Oceanography Postdoctoral Fellowship program. A.A and G.T acknowledge the ERC grant no. 305377 (BRIDGE). CIS acknowledges a research start-up grant from Victoria University of Wellington. C.O. is supported by the NERC Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics. The Earth Observing-1 (EO-1) spacecraft is managed by NASA’s Goddard Space Flight Center, Greenbelt, Maryland, USA. We thank Daniel Tran, Steve Chien and Joshua Doubleday (all JPL) and Stuart Frye (GSFC) for their assistance in acquiring EO-1 data. We are very grateful to James White, Patrick Allard
and an anonymous reviewer for comments which significantly improved
the quality of this manuscript.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

References

Aiuppa A, Shinohara H, Tamburello G, Giudice G, Liuzzo M, Moretti R (2011) Hydrogen in the gas plume of an open-vent volcano, Mount Etna, Italy. J Geophys Res 116:8 PP

Aiuppa A, Giudice G, Liuzzo M, Tamburello G, Allard P, Calabrese S, Chaplygin I, McGonigle AJS, Taran Y (2012) First volatile inventory for Gorely volcano. Kamchatka Geophys Res Lett 39:L06307

Aiuppa A, Robidoux P, Tamburello G, Conde V, Galle B, Avard G, Bagnato E, De Moor JM, Martínez M, Muñóz A (2014) Gas measurements from the Costa Rica–Nicaragua volcanic segment suggest possible along-arc variations in volcanic gas chemistry. Earth Planet Sci Lett 407:134–147

Aiuppa A, Bani P, Moussallam Y, Di Napoli R, Allard P, Gunawan H, Hendrasto M, Tamburello G (2015) First determination of magma-derived gas emissions from Bromo volcano, eastern Java (Indonesia). J Volcanol Geotherm Res 304:206–213

Barriga, V. M., 1951. Los terremotos en Arequipa (1582–1868)

Birnie RW, Hall JH (1974) The geochemistry of El Misti volcano, Peru fumaroles. Bull Volcanol 38:1

Boudon G, Balcone-Boissard H, Villemant B, Morgan DJ (2015) What factors control superficial lava dome explosivity? Sci Rep 5:14551

Buck AL (1981) New equations for computing vapor pressure and enhancement factor. J Appl Meteorol 20:1527–1532

Chávez Chávez, J. , . 1992. La erupción del volcán Misti, Pasado, presente, futuro: Arequipa

Chien S, Sherwood R, Tran D, Cichy B, Rabideau G, Castano R, Davis A, Mandl D, Trout B, Shulman S, Boyer D (2005) Using autonomy flight software to improve science return on Earth Observing One. J Aerosp Comput Inf Commun 2:196–216

Cobeñas G, Thouret J-C, Bonadonna C, Boivin P (2012) The c.2030yr BP Plinian eruption of El Misti volcano, Peru: eruption dynamics and hazard implications. J Volcanol Geotherm Res 241–242:105–120

Davies AG, Chien S, Baker V, Doggett T, Dohm J, Greeley R, Ip F, Castano R, Cichy B, Rabideau G, Tran D, Sherwood R (2006) Monitoring active volcanism with the autonomous sciencecraft experiment on EO-1. Remote Sens Environ 101:427–446

Davies AG, Chien S, Tran D, Doubleday J (2015) The NASA volcano sensor web, advanced autonomy and the remote sensing of volcanic eruptions: a review. Geol Soc Lond Spec Publ 426:SP426.3

Decker R, Christiansen R (1984) Explosive eruptions of Kilauea volcano, Hawaii in: Explos. Volcanism inception Evol. Hazards. National Academy Press, Washington, pp 122–132

Edmonds M, Oppenheimer C, Pyle DM, Herd RA, Thompson G (2003) SO2 emissions from Soufrière Hills volcano and their relationship to conduit permeability, hydrothermal interaction and outgassing regime. J Volcanol Geotherm Res 124:23–43

Finizola A, Lénat J-F, Macedo O, Ramos D, Thouret J-C, Sortino F (2004) Fluid circulation and structural discontinuities inside Misti volcano (Peru) inferred from self-potential measurements. J Volcanol Geotherm Res 135:343–360

Frost BR (1991) Introduction to oxygen fugacity and its petrologic importance. Rev Mineral Geochem 25:1–9

Gerlach TM, McGee KA, Doukas MP, 2008 Emission Rates of CO2, SO2, and H2S, scrubbing, and preeruption excess volatiles at Mount St. Helens, 2004–2005, in: Volcano Re Kindled Renewed Erupt. Mt. St Helens 2004–2006. Geological Survey Professional Paper 1750, pp. 543–571

Giggenbach WF (1987) Redox processes governing the chemistry of fumarolic gas discharges from White Island, New Zealand. Appl Geochem 2:143–161

Giggenbach, W., 1996. Chemical composition of volcanic gases. Scarpa R Tilling RI Eds Monit. Mitig. Volcano hazards 202–226.

Gillespie A, Rokugawa S, Matsunaga T, Cothern JS, Hook S, Kahle AB (1998) A temperature and emissivity separation algorithm for Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images. IEEE Trans Geosci Remote Sens 36:1113–1126

Gonnermann HM, Manga M (2003) Explosive volcanism may not be an inevitable consequence of magma fragmentation. Nature 426:432–435

Gonzales, K., 2009. Monitoreo insar de los volcanes misti, ubinas y ticsani - 20098.

Hale AJ, Mühlhaus H-B (2007) Modelling shear bands in a volcanic conduit: implications for over-pressures and extrusion-rates. Earth Planet Sci Lett 263:74–87

Hammouya G, Allard P, Jean-Baptiste P, Parello F, Semet MP, Young SR (1998) Pre- and syn-eruptive geochemistry of volcanic gases from Soufrière Hills of Montserrat, West Indies. Geophys Res Lett 25: 3685–3688

Hantke, G., Parodi, A., 1966. Catalogue of the active volcanoes of the world

Harpel CJ, de Silva S, Salas G (2011) The 2 ka eruption of Misti volcano, Southern Peru—the most recent Plinian eruption of Arequipa’s iconic volcano. Geol Soc Am Spec Pap 484:1–72

Heap MJ, Farquharson JI, Wadsworth FB, Kolzenburg S, Russell JK (2015) Timescales for permeability reduction and strength recovery in densifying magma. Earth Planet Sci Lett 429:223–233

Holland ASP, Watson IM, Phillips JC, Caricchi L, Dalton MP (2011) Outgassing processes during lava dome growth: insights from Santiaguito lava dome. Guatemala J Volcanol Geotherm Res 202: 153–166

Jaupart C, Allègre CJ (1991) Gas content, eruption rate and instabilities of eruption regime in silicic volcanoes. Earth Planet Sci Lett 102:413–429

Jay JA, Welch M, Pritchard ME, Mares PJ, Mnich ME, Melkonian AK, Aguilera F, Naranjo JA, Sunagua M, Clavero J, 2013. Volcanic hotspots of the central and southern Andes as seen from space by ASTER and MODVOLC between the years 2000 and 2010. Geol. Soc. Lond. Spec. Publ. 380, SP380.1

Kennedy BM, Wadsworth FB, Vasseur J, Ian Schipper C, Mark Jellinek A, von Aulock FW, Hess K-U, Kelly Russell J, Lavallée Y, Nichols ARL, Dingwell DB (2016) Surface tension driven processes densify and retain permeability in magma and lava. Earth Planet Sci Lett 433:116–124

Komorowski J-C, Jenkins S, Baxter PJ, Picquout A, Lavigne F, Charbonnier S, Gertisser R, Preece K, Cholik N, Budi-Santoso A, Surono E (2013) Paroxysmal dome explosion during the Merapi 2010 eruption: processes and facies relationships of associated high-energy pyroclastic density currents. J Volcanol Geotherm Res 261:260–294

Mellors RA, Waitt RB, Swanson DA (1988) Generation of pyroclastic flows and surges by hot-rock avalanches from the dome of Mount St. Helens volcano. USA Bull Volcanol 50:14–25

Melnik O, Sparks RSJ (1999) Nonlinear dynamics of lava dome extrusion. Nature 402:37–41
Melnik O, Sparks RSJ (2005) Controls on conduit magma flow dynamics during lava dome building eruptions. *J Geophys Res Solid Earth* 110:B02209

Moussallam Y, Oppenheimer C, Aiuppa A, Giudice G, Moussallam M, Kyle P (2012) Hydrogen emissions from Erebus Volcano. *Antarctica Bull* Volcanol 74:2109–2120

Moussallam Y, Peters N, Ramírez C, Oppenheimer C, Aiuppa A, Giudice G (2014) Characterisation of the magmatic signature in gas emissions from Turrialba Volcano, Costa Rica. *Solid Earth* 5:1341–1350

Moussallam Y, Bani P, Curtis A, Barnie T, Moussallam M, Peters N, Schipper CI, Aiuppa A, Giudice G, Amigo Á, Velasquez G, Cardona C (2016) Sustaining persistent lava lakes: observations from high-resolution gas measurements at Villarrica volcano. *Chile Earth Planet Sci Lett* 454:237–247

de Murúa, F. M., 1946. Historia del origen, genealogía real de los Reyes Incas del Perú

de Murúa, F. M. 1987. Historia general del Perú

Oppenheimer C, Edmonds M, Francis P, Burton M (2002) Variation in HCl/SO2 gas ratios observed by Fourier transform spectroscopy at Soufrière Hills Volcano Montserrat. *Geol. Soc. Lond. Mem* 21:621–639

Parodi, A., 1966. The active volcanoes of Peru, in: Cat. Act. Volcanoes Solfatara Fields Colomb. Ecuad. Peru. International Association of Volcanology, Rome, pp. 65–73.

Pearlman JS, Barry PS, Segal CC, Shepanski J, Beiso D, Carman SL (2003) Hyperion, a space-based imaging spectrometer. *IEEE Trans Geosci Remote Sens* 41:1160–1173

Pritchard ME, Simons M (2004) An InSAR-based survey of volcanic deformation in the central Andes. *Geochem Geophys Geosystems* 5:Q02002

Sandri L, Thouret J-C, Constantinescu R, Biass S, Tonini R (2014) Long-term multi-hazard assessment for El Misti volcano (Peru). *Bull Volcanol* 76:771

Sato H, Fujii T, Nakada S (1992) Crumbling of dacite dome lava and generation of pyroclastic flows at Unzen volcano. *Nature* 360:664–666

Shinohara H (2005) A new technique to estimate volcanic gas composition: plume measurements with a portable multi-sensor system. *J Volcanol Geotherm Res* 143:319–333

Simkin, T., Siebert, L., 1994. Volcanoes of the world (2nd edition)

Sparks RSJ (1997) Causes and consequences of pressurisation in lava dome eruptions. *Earth Planet Sci Lett* 150:177–189

Stevenson DS (1993) Physical models of fumarolic flow. *J Volcanol Geotherm Res* 57:139–156

Stull, D. R., Westrum, E. F., Sinke, G. C., 1969. The chemical thermodynamics of organic compounds

Symonds RB, Gerlach TM, Reed MH (2001) Magmatic gas scrubbing: implications for volcano monitoring. *J Volcanol Geotherm Res* 108:303–341

Tamburello G (2015) Ratiocalc: software for processing data from multicolumnar volcanic gas analyzers. *Comput Geosci* 82:63–67

Tepley FJ, Silva SD, Salas G, 2013. Magma Dynamics and Petrological Evolution Leading to the VEI 5 2000 BP eruption of El Misti Volcano, Southern Peru. *J Petrol* egt040

Thouret J-C, Legros F, Gourgaud A, Salas G, Juvigne E, Gilot E, Uribe M, Rodriguez A (1995) Un exemple de prévision des risques volcaniques au Pérou méridional (région d'Arequipa), fondé sur l'étude de l'activité éruptive récente du strato-volcan El Misti. *Comptes Rendus Académie Sci. Sér. 2 Sci. Terre Planètes* 320:923–929

Thouret J-C, Finizola A, Fornari M, Legeley-Padovani A, Suni J, Frechen M (2001) Geology of the Misti volcano near the city of Arequipa. *Peru Geol Soc Am Bull* 113:1593–1610

Ungar SG, Pearlman JS, Mendenhall JA, Reuter D (2003) Overview of the Earth Observing One (EO-1) mission. *IEEE Trans Geosci Remote Sens* 41:1149–1159

Watts RB, Herd RA, Sparks RSJ, Young SR (2002) Growth patterns and emplacement of the andesitic lava dome at Soufrière Hills Volcano, Montserrat. *Geol Soc Lond Mem* 21:115–152

Yamaguchi Y, Kahle AB, Tsu H, Kawakami T, Pniel M (1998) Overview of advanced spaceborne thermal emission and reflection radiometer (ASTER). *IEEE Trans Geosci Remote Sens* 36:1062–1071

Zamácola Jaureguí, J. D. 1958. Apuntes para la historia de Arequipa. Primer festival del libro arequipeño