# Multi-parametric field experiment links explosive activity and persistent degassing at Stromboli

Simon Thivet, Andrew J. L. Harris, Lucia Gurioli, Philipson Bani

► To cite this version:

Simon Thivet, Andrew J. L. Harris, Lucia Gurioli, Philipson Bani. Multi-parametric field experiment links explosive activity and persistent degassing at Stromboli. Frontiers in Earth Science, 2021, 9, 10.3389/feart.2021.669661. insu-03708923

HAL Id: insu-03708923

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

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Multi-Parametric Field Experiment Links Explosive Activity and Persistent Degassing at Stromboli

Simon Thivet¹,²*, Andrew J. L. Harris¹, Lucia Gurioli¹, Philipson Bani¹, Talfan Barnie³, Maxime Bombrun⁴ and Emanuele Marchetti⁵

¹Laboratoire Magmas et Volcans, Université Clermont Auvergne-CNRS-IRD, OPGC, Clermont-Ferrand, France, ²Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Munchen, Germany, ³Nordic Volcanological Center, Institute of Earth Science, University of Iceland, Reykjavik, Iceland, ⁴Department of Information Technology, Division of Visual Information and Interaction, and Science for Life Laboratory, Uppsala University, Uppsala, Sweden, ⁵Dipartimento di Scienze Della Terra, Università Degli Studi di Firenze, Firenze, Italy

Visually unattainable magmatic processes in volcanic conduits, such as degassing, are closely linked to eruptive styles at the surface, but their roles are not completely identified and understood. To gain insights, a multi-parametric experiment at Stromboli volcano (Aeolian Islands, Italy) was installed in July 2016 focusing on the normal explosive activity and persistent degassing. During this experiment, gas-dominated (type 0) and particle-loaded (type 1) explosions, already defined by other studies, were clearly identified. A FLIR thermal camera, an Ultra-Violet SO₂ camera and a scanning Differential Optical Absorption Spectroscopy were deployed to record pyroclast and SO₂ masses emitted during individual explosions, as well as persistent SO₂ fluxes, respectively. An ASHER instrument was also deployed in order to collect ash fallouts and to measure the grain size distribution of the samples. SO₂ measurements confirm that persistent degassing was far greater than that emitted during the explosions. Further, we found that the data could be characterized by two periods. In the first period (25–27 July), activity was mainly characterized by type 0 explosions, characterized by high velocity jets. Pyroclast mass fluxes were relatively low (280 kg/event on average), while persistent SO₂ fluxes were high (274 t/d on average). In the second period (29–30 July), activity was mainly characterized by type 1 explosions, characterized by low velocity jets. Pyroclast mass fluxes were almost ten times higher (2,400 kg/event on average), while persistent gas fluxes were significantly lower (82 t/d on average). Ash characterization also indicates that type 0 explosions fragments were characterized by a larger proportion of non-juvenile material compared to type 1 explosions fragments. This week-long field experiment suggests that, at least within short time periods, Stromboli's type 1 explosions can be associated with low levels of degassing and the mass of particles accompanying such explosive events depends on the volume of a degassed magma cap sitting at the head of the magma column. This could make the classic particle-loaded explosions of Stromboli an aside from the true eruptive state of the volcano. Instead, gas-dominated explosions can be associated with high levels of degassing and are indicative of a highly charged (with gas) system. We thus suggest that relatively deep magmatic processes, such as persistent degassing and slug formation can

**OPEN ACCESS**

**Edited by:**
Rosa Anna Corsaro,
National Institute of Geophysics and
Volcanology, Italy

**Reviewed by:**
Sonia Calvari,
National Institute of Geophysics and
Volcanology, Italy
Jonathan Dehn,
University of Alaska Fairbanks,
United States

*Correspondence:
Simon Thivet
simon.thivet@lmu.de

**Specialty section:**
This article was submitted to
Volcanology,
a section of the journal
Frontiers in Earth Science

**Received:** 19 February 2021

**Accepted:** 14 May 2021

**Published:** 28 May 2021

**Citation:**
Thivet S, Harris AJL, Gurioli L, Bani P, Barnie T, Bombrun M and Marchetti E (2021) Multi-Parametric Field Experiment Links Explosive Activity and Persistent Degassing at Stromboli. *Front. Earth Sci.* 9:669661. doi: 10.3389/feart.2021.669661
rapidly influence the superficial behavior of the eruptive conduit, modulating the presence or absence of degassed magma at the explosion/fragmentation level.

**Keywords:** particle size distribution, persistent degassing, pyroclast mass flux, real-time monitoring, strombolian explosions, explosive degassing, ash characterization

INTRODUCTION

Characterizing physical processes that drive explosive eruptions, as well as post-fragmentation transport and deposition of the eruptive products, are mandatory if we are to understand the associated volcanic hazards and risks (e.g., Barberi, 1993; Nave et al., 2010; Rosi et al., 2013). A number of modern geophysical and sampling techniques allow to define key eruptive parameters that inform on those physical processes. These include gas flux and composition (Aiuppa et al., 2009; Burton et al., 2009; Barnie et al., 2015; Delle Donne et al., 2017; Pering et al., 2020), pyroclast dynamics and size distributions (e.g., Patrick et al., 2007; Harris et al., 2012; Taddeucci et al., 2012), and the physical and chemical characteristics of the eruptive products (e.g., D'Oriano et al., 2010; Pioli et al., 2014; Gurioli et al., 2015). When measured together, these parameters can allow well-constrained models to be built for the physical processes operating before and at fragmentation, as well as during particle transport and deposition. We therefore implemented an experimental multi-parametric monitoring network at Stromboli volcano (Italy) to synchronously measure all of these parameters so as to constrain the transport and deposition characteristics of Stromboli's explosive activity, and hence understand the magma ascent and fragmentation conditions that drive the activity.

Stromboli has been characterized by mildly explosive, “Strombolian” activity for at least 2000 years (Mercalli 1907; Barberi, 1993; Rosi et al., 2000), from a stable system of three summit craters (Washington 1917; Harris and Ripepe 2007a; **Figure 1**): North East (NE), Central (C) and South West (SW). This makes Stromboli an ideal “natural laboratory” (Harris and Ripepe 2007a) where we can study explosive volcanism and associated eruptive processes (Taddeucci et al., 2015) that can be tested at, transferred to and adapted for other systems characterized by “Strombolian” activity, such as Shishaldin, Erebus, Yasur (e.g., Vergniolle et al., 2004; Johnson et al., 2008; Gerst et al., 2013; Simons et al., 2020a; Simons et al., 2020b). The “normal” Strombolian activity on which we focus here is characterized by intermittent “mildly explosive” activity (Barberi, 1993) that typically involves 10²–10³ kg of material (Harris et al., 2013) in repeated emissions of gas, bombs, lapilli and ash in events that last a few to tens of seconds, (Rosi et al., 2013; Houghton et al., 2016). Normal Strombolian activity, and its associated geophysical (seismic and infrasonic) signals, is classically explained by the generation, ascent and bursting of large gas bubbles (slugs) at the free surface (e.g., Blackburn et al., 1976; Jaupart and Vergniolle 1988; Jaupart and Vergniolle 1989; Ntepe and Dorel 1990; Braun and Ripepe 1993; Vergniolle and Brandeis 1994; 1996; Parfitt and Wilson 1995; Chouet et al., 1999; Chouet et al., 2003). However, in-spite of the persistence of activity and commonality of source process, normal

Strombolian activity is highly variable in style, and three different types of normal explosions have been identified on the basis of plume type, componentry and ascent dynamics:

(i) Type 0 explosions are gas-dominated, involving high velocity jets moving at velocities of hundreds of meters per second, but only involving a few particles (Goto et al., 2014; Leduc et al., 2015; Gaudin et al., 2017).

(ii) Type 1 explosions are dominated by lapilli and bomb rich emissions (Patrick et al., 2007), dominated by particles that follow ballistic trajectories (Chouet et al., 1974) and involving masses of the order of 103 kg (Harris et al., 2013).

(iii) Type 2 explosions produce ash-rich plumes (Patrick et al., 2007), where plumes follow convective ascent dynamics (Figure 2).

These different plume types can be identified at most Strombolian systems, such as Yasur and Etna (Gaudin et al., 2017; Simons et al., 2020b), and have been increasingly associated

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | (A) Digital Elevation Model (DEM) of Stromboli Island, with the location of the Differential Optical Absorption Spectroscopy (DOAS) instrument. (B) DEM of the summit area of Stromboli volcano with the locations of the Forward Looking Infra-Red (FLIR) camera, the Ultra-Violet (UV) SO₂ camera, the ASH collector and sizER (ASHER) instrument (two different positions), and the DOAS instrument (two different positions). The three different active zones of Stromboli volcano are materialized: the South West (SW) crater, the Central (C) crater, and the North East (NE) crater.</figcaption></figure>
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | (A) Sequence of pictures of a Type 2a explosion at Stromboli volcano, representing the evolution of eruptive products (ballistics and plume) in the atmosphere since the explosion level. Orange numbers are described in (B). Terminology is based on Turner (1962), Turner (1969), Slawson and Csanady (1967). (B) Sketch of the different steps (marked in orange) of a normal explosion at Stromboli. Each step can be sampled and/or measured by different techniques and/or strategies, which are listed in the text. Instruments used in our field experiments are underlined in different colors following the one used in **Figure 1**.

with bubble bursting through a more (for type 1/2) or less (for type 0) “capped” surface; the cap comprising a layer of high-viscosity, degassed magma (Lautze and Houghton 2005, Lautze and Houghton 2007; Gurioli et al., 2014; Del Bello et al., 2015; Leduc et al., 2015; Capponi et al., 2016; Suckale et al., 2016; Oppenheimer et al., 2020) that can be involved also in more energetic activity (Caracciolo et al., 2021).

Between normal Strombolian explosions, activity is dominated by persistent degassing (Francis et al., 1993; Allard et al., 1994; Allard 1997; Harris and Stevenson 1997; Pering et al., 2016) and puffing, characterized by continuous, bursts of gas in a continuous train of “explosions” that are weakly (0.1 bar) over-pressured (Ripepe and Gordeev 1999; Harris and Ripepe 2007b). Again, this style of activity has been identified as common to many “Strombolian” systems, including Etna, Villarrica and Masaya (Gresta et al., 2004; Branan et al., 2008; Gurioli et al., 2008). Within this mixed (Strombolian explosions and persistent degassing) activity regime, outstanding questions are 1) related to the roles of, and relations between, these two co-existing activity

types in driving changes in activity level and style, and 2) related to the underlying (deep and shallow conduit system) physical processes. In setting up experiments to answer these questions, we can improve definition and design of monitoring networks, as well as processing, integration and interpretation of the data output, so as to better assess hazard types and risk levels at the associated system (e.g., Calvari et al., 2005; Ripepe et al., 2005; Bertolaso et al., 2009).

Hence, the aim of our study was to monitor and track, over a few days, the normal Strombolian and persistent degassing activity of Stromboli by 1) observing the explosive activity, 2) recording synchronous data concerning the solid (pyroclast) and gas (persistent and explosive) emissions of Stromboli volcano, and 3) collecting and characterizing ash fallout samples. Combining a carefully designed and deployed experiment, involving multiple field-based remote sensing measurements and samplings with post-eruptive laboratory analysis, is crucial to bridge the gap between visually unattainable magmatic processes (i.e., conduit conditions) and their superficial
manifestations (i.e., explosion style, pyroclast mass emission, particle size distribution, gas flux) as recordable using multispectral measurements spanning laser, ultraviolet and infrared. More specifically, we aim to answer the follow open questions:

(i) Is there any correlation between persistent and/or explosive gas and solid fluxes emitted by the normal Strombolian explosions?

(ii) Does ash characterization help in understanding the fragmentation mechanisms of the associated measured explosions?

(iii) What are the magmatic processes at the origin of the transitions observed in the eruptive styles, especially between type 0 and type 1 explosions?

Most importantly, we show how multi-parametric measurements, involving traditional and modern remote sensing technologies, can be deployed and coupled to monitor, track and understand activity levels at active persistently degassing Strombolian systems.

## METHODS

This multi-parametric field experiment involved deployment of three ground-based remote sensing instruments (Figure 1): a high speed thermal FLIR (Forward Looking Infra-Red) camera, a UV (Ultra-Violet) camera and a scanning DOAS (Differential Optical Absorption Spectroscopy); plus a new instrument capable of remotely sampling active fall out in real time using a laser barrier, an ASHER (ASH collector and sizER). These instruments were deployed around the active summit craters and on the lower flanks for a six-day-long period between July 25 and 30, 2016. All instruments were constantly geo-located and synchronized at the same local time in order to get a synchronized and exploitable dataset over the entire experiment duration. In addition, observational logs were kept during all summit visits detailing the time, location and type of each explosion (as summarized in Supplementary Table S1).

## Measurement Strategy and Instrument Layout

Normal explosions at Stromboli are generally understood to be sourced by the generation and ascent of a large gas bubble (or slug), generated at some depth, which then ascends to burst at the free-surface to cause fragmentation of the magma resident in the conduit (e.g., Chouet et al., 1999; Chouet et al., 2003; Burton et al., 2007a). After fragmentation, the mixture of gas and particles ascends the conduit to be released from the vent as an explosive emission (Ripepe et al., 2001; Ripepe et al., 2002). Thereafter, evolution of a normal explosive event at Stromboli follows five chronologically ordered steps (Figure 2A). The first is emission of the mixture of gas and particles from the vent as a jet (step 1 in Figure 2A), followed by initial plume ascent and development of a rooted thermal (step 2 in Figure 2A). In steps three and four,

bombs, blocks and larger lapilli following ballistic trajectories begin to land proximal to the vent (step 3 in Figure 2A), and a cloud of finer particles and gas detaches to drift in the wind (step 4 in Figure 2A). Finally, finer particles (ash) falls out from the cloud at distal locations (step 5 in Figure 2A). Each step can thus be sampled using different remote sensing and real-time sampling approaches to parameterize various components of the emission. These span the at-vent conditions in terms of exit velocity and particle size distribution to distal fall out conditions in terms of ash mass and componentry. We thus designed our experiment to constrain (in terms of gas/particle masses and particle sizes) the mixture comprising the plume, as well as the dynamics (velocities and fluxes) and nature (chemistry and texture) of the particles; all constraints that can be used to inform on the conduit conditions at the time of fragmentation. As sketched in Figure 2B, the network configuration was thus designed to use each of the steps to extract data to inform on the fullest possible range of parameters as follows:

(i) High speed thermal video of the mixture leaving the vent (step 1 in Figure 2B) allows to measure the pyroclast mass, the particle size distributions and the vent-leaving velocities of the particles (cf. Forward Looking Infra-red Camera: Insights Into Pyroclast Mass Emission).

(ii) The UV SO₂ camera then allows measurement of the SO₂ mass within the plume (step 2 in Figure 2B), and with the thermal camera data allows for pyroclast/gas mass ratio to be characterized (cf. Ultra-violet SO₂ Camera: Insights Into Explosive SO₂ Emission).

(iii) Vent-proximal fallout of ballistics (step 3 in Figure 2B) can be sampled to inform on the texture and chemistry of the conduit magma if conditions are ideal (e.g., Lautze and Houghton 2008; Gurioli et al., 2014; Leduc et al., 2015); however the geometry of the crater meant that this was not possible during our deployment.

(iv) The SO₂ flux associated with the gas-rich drifting cloud (step 4 in Figure 2B), as well as persistent degassing and puffing between explosions, can be measured with the DOAS instrument (cf. Differential Optical Absorption Spectroscopy: Insights Into Persistent SO₂ Degassing).

(v) Finally, sedimentation of particles from the drifting cloud (step 5 in Figure 2B) can be sampled and characterized in terms of grain size and terminal velocity by an ASHER instrument (cf. ASH Collector and sizER and Laboratory Ash Characterization: Insights Into Grain Size Distribution and Componentry). Sample return then allows textural, chemical and componentry analysis.

When tracked through time, this configuration allows to link changes of the emission properties to conduit conditions as witnessed by, for example, changes in particle mass, particle size distribution, componentry and/or gas flux, to changes in the eruptive style and the persistent gas flux. We next detail the new elements of our integrated, ground-based remote sensing array, i.e., the thermal camera, SO₂ camera, DOAS and ASHER, and the data processing steps.
Forward Looking Infra-Red Camera:
Insights into Pyroclast Mass Emission

A FLIR 655sc camera (Figure 3A) was used to derive the total pyroclast masses ejected during individual explosions using the particle tracking algorithm and particle parameter (velocity, dimensions, shape, volume and mass) calculation procedures described in Bombrun et al. (2014), Bombrun et al. (2015). This procedure allows detection of all emitted particles and removal of background pixels, which may contain hot particles lying on the ground as deposited during previous explosions. Pixel dimensions are computed, the long and short axis of each particle retrieved, and particle volume is calculated by assuming rotational symmetry about the long axis. The mass is calculated by multiplying the volume by an appropriate density: two densities from samples collected at Stromboli between 2008 and 2011 were defined so that the conversion

depending on particle size. While the density distribution for bombs (> 64 mm, following the Schmid, 1981 nomenclature) is tight with a mean and standard deviation of 1.8 ± 0.2 g/cm³, that for large lapilli (between 50 and 64 mm) forms a cluster with a mean and standard deviation of 1.0 ± 0.2 g/cm³ (Gurioli et al., 2013; Gurioli et al., 2014; Bombrun et al., 2015). These two values are used for the size-dependent conversion to mass, which gives an error on the final total mass of ± 10%.

The focal plane array of the FLIR camera is composed of uncooled microbolometers sensitive between 7 and 14 µm. The camera was installed viewing into the South West (SW) crater from the Pizzo Sopra la Fossa (PSF) over a distance from the vent of around 250 m (Figure 1B). The camera was equipped with a ×3.6 magnification lens and acquired 640 × 120 images at 200 frames per second for a rectangular zone immediately above the vent (see Bombrun et al. (2014), Bombrun et al. (2015) for targeting configuration). The

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

**FIGURE 3** | (A) FLIR camera aiming pyroclasts from the South West (SW) crater. (B) Ultra-Violet (UV) SO₂ camera aiming volcanic plumes from the SW crater. (C) ASH collector and sizER (ASHER) instrument measuring and collecting ash fallouts. In the background, panorama of the Central (on the left) and the North East (on the right) active craters of Stromboli volcano in July 2016. (D) Differential Optical Absorption Spectroscopy (DOAS) instrument with the summit area of Stromboli in the background.
focal length was 88.9 mm and the instantaneous field of view was 0.19 mrad, giving for the line of sight tilt angle of -30°, a pixel size of 5 cm.

## Ultra-Violet SO₂ Camera: Insights into Explosive SO₂ Emission

Explosive SO₂ masses (i.e., gas emitted by an explosion) were measured using an UV SO₂ camera system following Barnie et al. (2015). Our UV SO₂ camera system consists of two QSI-620s UV cameras and an Ocean Optics USB2000 + UV spectrometer (Figure 3B). The cameras are mounted with 25 mm f2.8 Universe Optics quartz lenses, giving a Field Of View (FOV) of 23° by 17.25°. Both lenses are fitted with Asahi Spectra 10 nm bandpass filters in front of the lens, one at 310 nm (SO₂ absorption) and one at 330 nm (no SO₂ absorption). The spectrometer is connected to an Ocean Optics 74-UV 10 mm f2 telescope oriented within the FOV of the cameras with an FOV of 5.7°.

The UV SO₂ camera was set up just next to the FLIR, also viewing into the SW crater (Figure 1B), and images and spectra were collected every 5 s. Each image was subject to a flat field correction using dark current and clear sky images and corrected for incidence angle wavelength shift effects (Smekens et al., 2015). Each corrected 330/310 nm image pair was then aligned using the phase correlation technique (Kuglin and Hines 1975; De Castro and Morandi 1987) and used to calculate an Apparent Absorption (AA) image. The location of the spectrometer FOV within the camera FOV was found by simulating camera AA for each spectrum, interpolating the spectral AA time series to the camera acquisition times, smoothing the camera AA time series spatially on a length scale of the spectrometer FOV, and finding the location within the camera FOV that gave the highest correlation coefficient between smoothed camera AA and interpolated spectral AA. The total column amount of SO₂ in ppmm (parts per million meters) was found using the iFit technique (Barnie et al., 2015; Burton and Sawyer 2016). The calibration relationship was found by fitting the camera AA time series at the center of the spectrometer FOV with the spectral ppmm interpolated to the image acquisition times. This relationship was then used to calibrate all the AA images to ppmm. Images containing explosions were identified, the explosion plumes digitized and the total area and mean ppmm values calculated, and the total SO₂ mass derived using the ideal gas law.

## Differential Optical Absorption Spectroscopy: Insights Into Persistent SO₂ Degassing

During the 6 days of field recording we carried out DOAS measurements in a stationary scanning mode at three different locations. We used locations to capture best the total SO₂ persistent flux from the three summit craters (Figure 1B), so measurement sites were selected depending on wind direction and velocity so that the whole plume was scanned (Williams-Jones et al., 2006). On some days, measurements at the summit

were difficult due to wind conditions that caused the gas to pond as a “vog” (Businger et al., 2015) over the Pizzo Sopra la Fossa and the upper slopes of the *Sciara* del Fuoco. Accordingly, on July 25 and 27, combined vertical scanning measurements were performed at about 400 m from the active craters (870 m above sea level). On July 26 and 30, vertical scanning was carried out 2 km northeast of the crater (65 m above sea level). Finally, on July 28 and 29 the system was positioned at Pizzo Sopra la Fossa, around 150 m from the eastern crater rim (920 m above sea level), but performing sub-horizontal scanning. The spectrometer used was an Ocean Optics USB2000+ with a spectra range of 290–440 nm and a spectral resolution of 0.5 FWHM (Figure 3C). The SO₂ column amounts in ppmm were retrieved using DOAS calibration and standard analysis procedures (Platt and Stutz 2008). Reference spectra included in the non-linear fit were obtained by convolving high resolution SO₂ and O₃ cross sections with the instrument line shape (Voigt et al., 2001; Bogumil et al., 2003). A Fraunhofer reference spectrum and ring spectrum were also included in the fit. The total column amount of the plume cross section was then multiplied by the plume velocity to derive the SO₂ emission rates. In this work, for the vertical scanning, the plume dispersal speed was considered equivalent to the wind velocity, obtained from a hand-held anemometer. For the sub-horizon scanning, the plume ascent speeds were obtained using the thermal infrared camera and used as the velocity.

## ASH Collector and sizER and Laboratory Ash Characterization: Insights into Grain Size Distribution and Componenty

Characterizing the eruptive products, especially ash (particle diameter <2 mm following the Schmid 1981 nomenclature), is also a key to understand the eruptive processes involved in generating a Strombolian explosion. In order to collect and characterize the ash particles related to the measured explosions, an optical barrier (ASHER, Figure 3D) was deployed at two different locations around the summit craters (Figure 1B). Between July 25 and 28, it was installed at Pizzo Sopra la Fossa (at 915 m above sea level), while on July 29 it was installed closer to the SW crater (at 855 m above sea level).

The ASHER optical barrier consists of a laser and a receiver which detects particles when they cross the laser beam: a 90 CM-M140 laser (wavelength 635 ± 5 nm) is used to produce the linear (30 × 1 mm) laser beam. Operating at a sampling rate of 30 kHz, particles are detected from laser obscuration with the amplitude and duration of the absorption peak being representative of the particle diameter and terminal velocity, respectively. Ash particles are collected in a container below the laser beam. ASHER measurements are performed for 10 s every 30 s, allowing near-real time analysis of Grain Size Distribution (GSD). For this, the number of detected particles is provided and divided into 13 bins of diameter spanning from $Φ > = -1.5 (2,800 μm)$ to $Φ < = 3.5 (90 μm)$, with a step of 0.5 $Φ$, thus covering the coarse ash fractions.

Exploitable ash samples, emitted from the SW crater, were then analyzed. Laboratory sieving was performed to isolate the
different grain size bins with a step of 0.5 Φ and to compare with the syn-eruptive ASHER data. Selected ash fractions were prepared in polished epoxy grain mounts. Componentry analysis was performed on the 180–250 µm grain size bin, based on Back-Scattered Electron (BSE) mosaic images acquired with the Cameca SxFiveTactis Electon Probe Micro Analyzer (EPMA). In parallel, in-situ glass and crystal analysis were performed using the same EPMA instrument, with a 15 kV acceleration voltage and a focused spot size (1 µm) for crystals and a spot diameter of 20 µm for the glass. Ash morphology measurements were also carried out using the automatized morpho-grainsizer Malvern Morphologi G3, following the procedure described in Thivet et al. (2020a), Thivet et al. (2020b). This instrument allows measurement of the Apparent Projected shape of ASH (APASH) on hundreds of ash particles, listing several shape parameters (described in Leibrandt and Le Pennec 2015) for each analyzed particle. Roughness (solidity, which corresponds to the area of the projected ash shape divided by the convex hull area of the projected ash shape) and morphological (convexity, which corresponds to the convex hull perimeter of the projected ash shape divided by the perimeter of the projected ash shape) parameters are used to identify shape variations between the ash samples. These measurements were performed on the same grain size bin as the componentry analysis (i.e., 180–250 µm).

# RESULTS

## Field Observations

Summit-based observations made between 25 July and July 30, 2016 permitted to identify two significantly different periods in terms of explosion types (Figure 4A): the first between 25 and 27 July being dominated by Type 0 explosions, and the second between 29 and 30 July being dominated by Type 1 events (a transition period can also be identified on 28 July being dominated by Type 0 explosions but with the progressive occurrence of type 1 explosions). As shown by our field experiments results, these two main periods were also marked by systematic changes in particle velocities, numbers, and masses (Figure 4B), as well as SO₂ degassing (Figure 4C). All the field experiments results are displayed in Figure 4D and are further detailed in the following results sections. Weather conditions were relatively stable during all field measurements, with low (typically 3.5 m/s) wind speed, stable direction (generally SSW) and no precipitation, but sometimes with summit cloud cover (all observations are reported in Supplementary Table S1).

During the first period, between 25 and 27 July, an average of five explosions per hour was measured. These explosions were mainly from vents in the SW crater and were characterized by gas-dominated type 0 explosions, which accounted for 80% of the total recorded explosions. These explosions involved gas-roars, with no or very few particles heard landing. On the 28 July, loads of lapilli (particle diameter between 2 and 64 mm according to the Schmid, 1981 nomenclature) slightly increased. The 20% remaining events were silent type 2 (four 2a and four 2b) explosions, both within SW and North East (NE) craters.

These type 2 explosions produced thermals that rose up to around 300 m above the craters.

During the second period, on 29 and 30 July, a total of 27 events were counted over an observation period of five and a half hours, giving a frequency of around five explosions per hour on average. Again, these explosions were mainly occurring in the SW crater (with very scarce events in the C and NE craters), but were now characterized by particle-loaded type 1 explosions which accounted for 90% of the events. These explosions involved bomb-loaded ejections, which fell back into the crater or being sent out onto the Sciara del Fuoco to the north, prohibiting sampling. The 10% remaining events were silent type 2b explosions that occurred in the NE crater.

## Pyroclasts Mass Emission

Three type 0 explosions during the first period were fully recorded with the FLIR camera on the 27 July. The first recorded explosion yields a mass of just 5 kg of pyroclast distributed between 11 particles. The two other recorded explosions are represented through the FLIR image sequences shown in Figures 5A,B, corresponding, respectively, to the second (total pyroclast mass of 645 kg with 456 detected particles) and the third (total pyroclast mass of 190 kg with 64 detected particles) explosions recorded during this day. For all of these type 0 explosions, a relatively low pyroclast mass emission (280 kg per event on average) and a relatively low number of detected particles (180 per event on average) was measured (Figures 4B,D). On average, the mean particle size was 6 cm (with a maximum value of 46 cm) and the vent-leaving velocity was 50 m/s (with a minimum and maximum values of 8 and 140 m/s, respectively).

Instead, ten type 1 explosions were fully recorded during the second period on 29 July. The minimum of pyroclast mass recorded for one of these events was 434 kg distributed between 1,409 particles. In Figure 5C we give the FLIR image sequence for the explosion with the highest pyroclast mass emission recorded (2,590 kg distributed between 5,348 detected particles). Instead, in Figure 5D we give the FLIR image sequence of the explosion with the maximum number of detected particles (9,029, representing a total of 880 kg). For all recorded type 1 explosions in the second period, a relatively high pyroclast mass emission (2,400 kg per event on average, one order of magnitude higher than the type 0 explosions) and relatively high number of detected particles (5,260 per event on average, also one order of magnitude higher than the type 0 explosions) were measured (Figures 4B,D). On average, the mean particle size was the same as in the first period (6 cm; with a maximum measured value of 46 cm), but the vent-leaving velocity was lower (30 m/s, with a minimum and maximum values of 6 and 187 m/s, respectively). All parameters extracted from the FLIR camera data are available in Supplementary Table S2.

## Explosive SO₂ Mass Emission

During the entire survey period, two gas plume measurements associated with two different type 0 explosions within the SW crater, were exploitable; both collected on 27 July (Figure 6). Other measurements were not good enough due to heavy
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Field experiment results performed between July 25, 2016 and July 30, 2016. (A) Field observations. (B) Summary of the pyroclast mass emissions measured with the Forward Looking Infra-Red (FLIR) camera. (C) Summary of the persistent SO₂ fluxes measured with the Differential Optical Absorption Spectroscopy (DOAS) instrument. (D) Detailed timeline of the performed measurements. Red squares represent the pyroclast mass emissions during type 0 and 1 explosions, calculated from the FLIR acquisitions; the size of the red squares is proportional to the number of the detected particles. Green squares are the two explosive SO₂ measurements from the Ultra-Violet (UV) SO₂ camera. Yellow circles represent the explosions related to the three exploitable ash samples collected by the ASH collector and sizer (ASHER) instrument.</figcaption></figure>

condensation within the gas plumes or cloud formation. These two exploitable measurements yielded SO₂ masses (Figure 4D) of 2.61 ± 0.67 kg (explosion measured at 14:54; all times are local) and 2.07 ± 0.50 kg (explosion measured at 15:21). Values, assumptions and conversions applied to obtain these masses are given in Supplementary Table S3.

## Persistent SO₂ flux

The persistent SO₂ flux (Figures 4C,D) also clearly discriminates the two periods. During the first period, the 121 Scanning DOAS measurements emphasize a persistently high SO₂ flux that fluctuated between 73 ± 11 t/d (= 0.84 kg/s) and 693 ± 6 t/d (= 8.02 kg/s) with a mean value of 274 t/d (= 3.17 kg/s). The second period, that includes 130 Scanning DOAS profiles, highlights a drop in the SO₂ flux. This is more than three time less than in the first period, with an average of 82 t/d (= 0.95 kg/s), the minimum flux being 26 ± 14 t/d (= 0.30 kg/s) and the maximum being 206 ± 4 t/d (= 2.38 kg/s). All values and errors are given in Supplementary Table S4.

## Ash Fallout Characteristics

Two ash samples were collected from two distinct type 0 explosions on the 25 and 27 July, and an ash sample representative of six type 1 explosions was collected on 29 July (Figure 4D). Ash fallouts were systematically experienced around 2 min after the explosions occurred. Syn-eruptive GSDs directly measured by the ASHER were comparable with laboratory sieving data with a similar grain size range of between 63 and 710 µm for all samples (Figure 7). No specific GSD differences appear between the type 0 and type 1 ash particles, with GSD modes ranging between the 180–250 µm and 355–500 µm bins for laboratory

sieving (wt%) and between the 125–180 and 180–250 µm bins for ASHER syn-eruptive data (Figure 7).

Componentry analysis performed on the three ash samples (Supplementary Figures S1–S3), and on the 180–250 µm grain size bin, identifies four different types of particles within the bulk samples (Figure 8A), which are described hereafter in order of abundance. The first are sideromelane particles that are characterized by glassy matrixes with the presence of rounded vesicles and very scarce microlites. These particles either show fluidal, spongy or dense textures. They appear transparent to light brown under the binocular microscope. The second are tachylite particles that are characterized by fully microcrystalline matrixes with the presence of very scarce irregular-shaped vesicles. They appear dark and opaque under the binocular microscope. Thirdly, phenocryst-free particles are observed (olivine, pyroxene and plagioclase in order of decreasing abundance) and counted as such as soon as the phenocryst represents more than 50% of the entire particle. Finally, scarce altered particles are observed. These are characterized by highly altered features. They often appear red under the binocular microscope. The componentry analysis shows similar results between the two type 0 samples. However, slight differences are observed between the ash collected between the type 0 and type 1 explosions. Sideromelane particles are less abundant in the type 0 (49% on average) than in the type 1 (64%) samples. Tachylite particles are also more abundant in the type 0 samples (18% on average) than in the type 1 sample (8%). Similar proportions of phenocryst-free and altered particles are found between the two sample types (Figure 8A). APASH measurements reveal only very subtle differences between the
<figure><img src="image_7.png" /></figure>

**FIGURE 5** | Sequences of images of four explosions in the South West (SW) crater, recorded with the Forward Looking Infra-Red (FLIR) camera. (A) Type 0 explosion recorded the 27 July at 13:43 (all times quoted in this paper are local times). At *t* = 0 s, only intermittent gas puffing is observed; between *t* = 0.90 s and *t* = 1.85 s, the explosion begins, a gas jet is visible with temperature up to 200°C, with the synchronous emission of scarce particles; between *t* = 3.35 s and *t* = 4.51 s, the gas jet decreases, still with the emission of scarce particles; at *t* = 8.32 s, the explosion ends. (B) Type 0 explosion recorded the 27 July at 16:26. At *t* = 0 s, only intermittent gas puffing is observed; between *t* = 1.21 s and *t* = 1.46 s the explosion begins, a gas jet is visible with temperature more than 300°C; between *t* = 2.34 s and *t* = 2.69 s, the gas jet decreases in intensity with the emission of very scarce particles; at *t* = 3.88 s, the gas jet ends. (C) Type 1 explosion recorded the 29 July at 14:44. At *t* = 0 s, no particles and no gas puffing are visible; at *t* = 1.79 s, the explosion starts with the emission of a lot of particles; at *t* = 3.00 s, the explosions intensity reach its maximum in term of pyroclast mass emission; between *t* = 4.50 s and *t* = 8.24 s the explosion decreases; at *t* = 16.98 s, the explosion ends with still hot pyroclasts on the crater walls. (D) Type 1/2 explosion recorded the 29 July at 16:16. At *t* = 0 s, no particles and no gas puffing are visible; at *t* = 0.86 s, the explosion starts with the emission of a lot of particles with a distinguishable small ash and/or gas plume; between *t* = 1.77 s and *t* = 7.43 s the explosion oscillates in term of intensity with the emission of a lot of relatively fine particles (lapilli); at *t* = 15.01 s, the explosion ends with still hot pyroclasts on the crater walls.

type 0 and type 1 ash samples (Figure 8B) which can be attributed to
the slight differences in the component proportions between the
corresponding samples. In the type 1 ash sample, another particle
population (highlighted by the red circle in Figure 8B) appears,
which is mainly characterized by the occurrence of particles with a
higher morphological roughness range, mainly composed of
sideromelane particles. In-situ glass chemical analysis performed
on the sideromelane particles (Figure 8C) show that the emitted
magmas, both during type 0 and type 1 explosions, are part of the
Highly Porphyritic (HP) composition field, which is typical of the
normal Stromboli activity. Only scarce glass fragments fall in the
stagnant magma composition field as defined by Gurioli et al. (2014).
No significant chemical differences are observed between the two
types of explosions. Raw data on the ash particle characterization is
available in Supplementary Table S5.

DISCUSSION

Link Between the Pyroclast Mass and
Persistent Degassing

Our vent-leaving pyroclast emission parameters recorded
during 13 different Strombolian explosions within the SW

crater (Figures 4B,D; Supplementary Table S2) are
consistent with those already recorded during normal explosive
activity at Stromboli which show a typical range of 10–103 kg (e.g.,
Chouet et al., 1974; Blackburn et al., 1976; Weill et al., 1992; Ripepe
et al., 1993; Hort and Seyfried 1998; Patrick et al., 2007; Taddeucci
et al., 2012; Harris et al., 2013; Bombrun et al., 2015; Leduc et al.,
2015; Gaudin et al., 2017; Pioli and Harris 2019). Our
measurements confirm that gas-dominated type 0 explosions
(Figures 5A,B) produce relatively low total pyroclast masses
(between 5 and 645 kg for a single event) with low numbers of
particles (between 11 and 456 for a single event). Instead, type 1
explosions (Figures 5C,D) produce relatively high total pyroclast
masses (between 434 and 2,589 kg for a single event) and involve
higher numbers of particles (between 1,409 and 9,029 for a single
event).

Our two explosive SO₂ mass measurements show relatively
low values (between 2 and 3 kg per event; **Figure 4D,**
**Supplementary Table S3),** which fall within the lowest range
of previous measurements made at Stromboli (e.g., Mori and
Burton 2009; Tamburello et al., 2012; Barnie et al., 2015).
Likewise our 251 persistent SO₂ measurements of between 26
and 693 t/d; **Figures 4C,D,** **Supplementary Table S4)** overlap
with other studies, where Allard et al. (1994) found persistent SO₂
<figure><img src="image_8.png" /></figure>

**FIGURE 6** | Sequences of images of the two exploitable SO₂ plumes above the South West (SW) crater, 27 July, recorded with the Ultra-Violet (UV) SO₂ camera. (A, B) Plume recorded at 14:54. (C, D) Explosion recorded at 15:21. (A, C) are 330 nm images (no SO₂ absorption). (B, D) are 310 nm images (SO₂ absorption). The SO₂ plume is visible on the 310 nm images in black. Image pairs framed in red were chosen to calculate Apparent Absorption (AA) images and calibrated SO₂ images in ppm that are represented in (E, F) for the 14:54 and 15:21 explosions, respectively. Eruption plumes are outlined in red. Histograms to the right of both images show the image (gray curve) and plume (red curve) ppm distributions.

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

**FIGURE 7** | Grain size of the three collected ash samples measured both by the ASH collector and sizER (ASHER) instrument and by laboratory sieving.

degassing fluxes between 340 and 1,420 t/d. This highlights that persistent degassing represents a variable, but much higher flux, than the explosive SO₂ flux. If we take a typical event frequency of five explosions per hour we obtain a total explosive flux of just 0.3 t/d. This is in agreement with Mori and Burton (2008) who found that the explosive SO₂ flux at Stromboli represents a relatively small contribution (3–8%) of the total degassing budget compared to persistent flux (92–97%). Likewise,

Tamburello et al. (2012) showed that puffing, which is generally located within the Central crater, emits twice as much gas as Strombolian explosions. Thus, although Stromboli’s explosions make for a spectacular distraction, the system is really dominated by persistent degassing.

Previous studies at Stromboli, especially before or during uncommon eruptive activity (i.e., paroxysmal or effusive, see terminology in Rosi et al., 2013), already showed that explosive
<figure><img src="image_10.png" /><figcaption>**FIGURE 8** | (A) Componentry analysis of the three collected samples performed within the 180–250 µm grain size bin and using Back-Scattered Electron (BSE) images acquired with a Scanning Electronic Microscope (SEM). The percentage of each component is derived from the number of particles identified and counted manually (the numbers of counted particle are represented in parentheses). (B) Apparent Projected shape of ASH (APASH) for the STR160725_1 (type 0 explosion) and STR160729_12 (type 1 explosions) samples, represented within the Solidity vs. Convexity graphs. Linear trends are represented by the dashed lines. The red circle shows the same area in both graphs, evidencing the occurrence of another ash population for the STR160729_12 sample. (C) *In-situ* glass compositions of the sideromelane particles. The different composition fields were reconstructed from Metrich (2005), Gurioli et al. (2014), Pioli et al. (2014), Leduc et al. (2015) data.</figcaption></figure>
degassing and persistent degassing fluxes are correlated (Delle Donne et al., 2017), which are also correlated with the explosion amplitudes (Ripepe et al., 2007; Ripepe et al., 2009; Valade et al., 2016). Thus, these parameters have become critical monitoring and potential precursors to uncommon and potentially dangerous eruptive activity. This is in contrast to the case documented here for normal, very mild, activity: our specific multi-parametric field experiment, performed on the short period of 25–30 July 2016 highlights that the explosive dynamics (explosion type, particle mass and number) observed on the main active vent within the SW crater are not at-all coupled with the total persistent SO₂ degassing flux (Figure 4). This decoupling is indeed revealed by the two distinct periods of degassing and activity style. During the first period, the persistent SO₂ degassing flux was relatively high (274 t/d on average, with peaks near 700 t/d), but activity was dominated by type 0 explosions and emitted particle masses were low. During the second period, the persistent SO₂ degassing flux declined by more than three times (to 82 t/d on average, with lowest values near 25 t/d) compared to the first period, and mainly type 1 explosions occurred (without any type 0 explosions) and the mass and number of tephra particles increased by two and three orders of magnitude, respectively. To sum up, we observe a strong anti-correlation between the amount of magma being degassed and the mass of pyroclasts being ejected from the main active vent (SW crater), as already documented by Leduc et al. (2015). We note, though, that during our observation period, the NE crater was not at all active, where we witnessed just seven exceedingly weak type 2 eruptions emitting ash to no more than a few tens of meters above the crater rim during the observation period. The Central crater, though, was active with its characteristic puffing behavior (e.g., Harris and Ripepe 2007a; Ripepe et al., 2007; Tamburello et al., 2012), emitting around one gas puff every 2 s, but not being the source of significant normal explosive activity with only two events observed during the whole observation period; the primary source of explosive activity being our targeted vent: the SWC.

## Insights From ASHER Measurements and Ash Characterization

During this experiment, the ASHER was installed for the first time in an active fall out field. It was thus important to check and validate its operation. The number of counted ash particles by the ASHER (384, 33, and 381 for the STR160725_1, STR160727_7 and STR160728_12 samples, respectively) differ from the actual number of particles counted in the bulk sample (830, 650, and, 7,974, respectively). This is due to the sampling interval of the instrument. ASHER GSDs are slightly skewed toward finer grain size bins compared the laboratory sieving analyses (Figure 7). This is explained by the fact that ASHER GSD histograms represent the number of detected particles while laboratory sieving ones represent the total weight of the each grain size bins (the finer the particle, the lighter they are). We thus argue that GSDs acquired with these two different methods show coherent results. Moreover, no significant GSD differences are observed between the type 0 and type 1 explosions, both looking

at the vent-leaving particles (FLIR measurements) and the ash fallout from the associated drifting cloud (ASHER and laboratory measurements). At vent, lapilli and bombs (up to 46 cm in width) are observed with a mean particle size of 7 and 6 cm for type 0 and type 1 explosions, respectively. For the sampled ash particles, GSDs systematically comprise a range between 63 and 710 µm, which represent the size range of the coarsest particles carried by the drifting cloud, the finest particles being transported further. Terminal velocity values measured by the ASHER were between 0.1–0.5 and 1–2 m/s depending the grain size bins, and are coherent with those recorded by Freret-Lorgeril et al. (2019).

The different ash components (Figure 8A) identified in our samples were consistent with previous observations made for samples collected during normal activity at Stromboli (e.g., D’Oriano et al., 2010; Lautze et al., 2012; Lautze et al., 2013). Textures of sideromelane particles are typically formed by HP magmas during normal Strombolian activity. No highly vesicular particles are observed, which suggest the absence of Low Porphyritic (LP) magmas. The dominant contribution of HP magmas to ash formation, with scarce stagnant superficial magma, is confirmed by the *in situ* glass analyses (Figure 8C). The scarcity of a stagnant magma, microlite-rich component in the ash particles, is explained by the fact that the stagnant magma is not efficiently fragmented, only forming coarse pyroclasts that are mostly found in the vent-proximal zone (Gurioli et al., 2014; Leduc et al., 2015) and it only represents the “crust” at the top of the degassed magma, so volumetrically limited (as explained in Caracciolo et al., 2021). Tachylite particles are considered as fully crystallized stagnant magma or non-juvenile recycled particles. Phenocryst-free particles come from both juvenile and non-juvenile fraction. Altered particles are considered as non-juvenile lithic particles. A subtle but important ash component variation is observed: the proportion of juvenile (sideromelane) particles is greater for the type 1 explosion sample (64%) than the type 0 explosion samples (49%) looking at the 180–250 µm grain size bin. This difference is also identifiable in the APASH signature of the ash particles (occurrence of sideromelane, juvenile particles visible within the red circle in Figure 8B). To sum up, no LP magma has been identified in the type 0 and type 1 products. The only significant difference between the ash samples from both explosion types is that juvenile particles are more present within the type 1 products than the type 0 ones. This observation suggests that type 1 explosions fragment a large proportion of degassed HP magma than type 0, in agreement with magma accumulation. In contrast, Type 0 emits less magma, but a larger proportion of non-juvenile particle than type 1 explosions because of free from magma accumulation and therefore prone at cleaning the conduit (e.g., Gurioli et al., 2014; Leduc et al., 2015; Gaudin et al., 2017).

## An Accidental Driver for Eruption Magnitude: Implications for Eruptive Processes

The frequencies of the Strombolian explosions were the same in both observation periods. This indicates that for the whole investigated time, the magnitude, intensity, frequency and style
of the explosive activity was not driven by, or linked to, deep processes that feed the persistent degassing flux that is changing from a period to another. Thus, measurements of the plume properties cannot be systematically used to assess or infer how active, in terms of persistent degassing, the system is. In fact, the reverse holds for the normal mild activity (without coupling between gas and particles): the more abundant is the material emitted, the less intense the degassing activity is (Figure 4). This confirms that the magma erupted at such systems is, at least for some periods during normal Strombolian activity, in no way associated with the magma that is degassed to drive the persistent SO₂ flux at such systems. This is consistent with the idea that processes controlling fragmentation at Stromboli can be driven from relatively superficial parameters, where the explosivity is controlled by the presence, and thickness, of a thin and viscous cap on the magma column (e.g., Burton et al., 2007a; Gurioli et al., 2014; Barnie et al., 2015; Leduc et al., 2015; Gaudin et al., 2017; Simons et al., 2020). From these previous models, we can argue that type 0 explosions are characterized by a low content of particles because they occurred when the shallow part of the eruptive conduit contains no (or little) degassed magma. Instead, the presence of degassed magma adds particles, but poor fragmentation efficiency in the high viscosity caps result in bomb-dominated type 1 explosions. The presence of any recycled backfill material will then be enhancing type 2 explosions (e.g., Patrick et al., 2007; Capponi et al., 2016; Simons et al., 2020). Our data support this global model for type 0 and 1 explosions whereby the mass of the bomb load accompanying an explosive event depends on the volume of the degassed magma cap sitting at the head of the magma column. The thicker the degassed magma cap, the greater quantity of material is available for the slug to fragment so that the emission carries a greater accidental bomb and lapilli load. Lower vent-leaving particle velocities (33 m/s on average) for type 1 explosions compared to the type 0 velocities (46 m/s on average) are also consistent with the model presented in Gaudin et al. (2017).

According to our data, the association of type 1 explosions with relatively low persistent gas flux periods suggests that periods of reduced degassing allow caps to build, while relatively high persistent gas flux periods do not allow cap formation. Our model is also consistent with the observed spatial pattern of activity. That is, the Central crater is believed to be the main, central pathway for gas ascent, hosting the hottest magma (Harris and Ripepe 2007b; Landi et al., 2011) and being the source of around 70% of the active degassing, that is represented by the gas released by the explosions and puffing (Tamburello et al., 2012). The lack of normal activity at this location would thus be explained by lack of cap formation at this focused source of degassing. As in other studies, the switch between low and high degassing periods occurs as a step, the change happening at most over an hour or two (e.g., Ripepe et al., 2002; Ripepe et al., 2008). Here, we also see that the cap responds over a similar time scale, but may be slightly lagged as the cap become reamed out at the onset of high degassing phase: eruptions continue at the same frequency, but steadily exhaust the superficial reservoir of stagnant cap magma.

As shown by Burton et al. (2007a); Burton et al. (2007 M. R.), gas slug formation and persistent degassing likely has its origin in

relatively deep levels between sea level and near the volcano/crust interface (around 2 km below sea level, and 3 km below the summit craters). Moreover, Burton et al. (2007a); Burton et al. (2007 M. R.) showed that the depth of slug formation may highly modulate the explosion intensities at the surface, the smallest explosions being produced by relatively small shallow slugs and the largest explosions being produced by large deep slugs. We thus support that while the mass of gas in an explosive normal event is controlled by deep processes, the mass of pyroclasts is rather controlled by shallow processes so that the two masses do not, necessarily, correlate. We suggest that these deep magmatic processes, especially persistent degassing fluctuations, can rapidly influence the superficial structure of the eruptive conduit, modulating the presence or the absence of degassed magma at the explosion/fragmentation level. When deep-seated magma-gas-coupled, enter the system, then the correlation between explosivity and mass flux and gas flux is reestablished.

## CONCLUSION

In this study, we present a multi-parametric monitoring methodology, based on ground-based remote sensing, for tracking activity at persistently degassing Strombolian systems, and present a framework within which to assess coupled variations in tracked parameters. Our combined measurements, on a short and specific period during normal Strombolian activity at Stromboli, highlight that low levels of degassing can be associated with establishment of a degassed magma cap at the conduit head, which contributes to the bomb load for type 1 eruptions. Instead, high levels of persistent degassing can be associated with the occurrence of gas-dominated type 0 gas jets, with very few particles.

Essentially, if there were no degassed cap, Strombolian eruptions at Stromboli would not be characterized by the bomb emission with which they have so famously become associated (e.g., Mercalli 1907; Macdonald 1972; Francis et al., 1993). Instead, pure normal Strombolian activity in a system which is purely driven from the deeper source would be mainly characterized by gas jetting. However, the top down influence exerted by the presence of a cap or scoria vent-choking material means that bombs and ash are invariably involved, with greater bomb loads implying an increasing influence of the cap and greater ash loads indicating a greater role of scoria choking (e.g., Capponi et al., 2016; Simons et al., 2020). We suggest that relatively low persistent degassing fluxes allow preservation, formation and accumulation of viscous caps within the shallow eruptive conduits near the fragmentation level. Instead, relatively high degassing rates promote cap-free conditions, so that gas bursts typically provoke type 0 explosions. Drops in persistent degassing flux thus rapidly favor the formation of a degassed magma cap, enhancing the occurrence of type 1 explosions. Type 2 explosions are probably more independent with this process as this explosion type is more likely provoked by relatively thick and dense cap formed by recycled backfill material.

In parallel to our results and associated model, it is important to keep in mind that monitoring measurements also usually show that mass loads for normal Strombolian activity at the NE and SW craters increase during periods of enhanced tremor, degassing and puffing, allowing prediction of the arrival of uncommon and
larger explosive or effusive activity at Stromboli (e.g., Ripepe et al., 2007; Ripepe et al., 2009; Valade et al., 2016; Delle Donne et al., 2017). This suggests that the system switches between a decoupling magma-gas behavior vs. a coupled magma-gas system due to the arrival of deep-seated gas-rich magma and that the cap formation may be variable from crater to crater, and can persist at some craters/vents during increases in gas flux, especially if the central crater takes the main burden. We conclude that solid masses erupted during normal/mild Strombolian activity could inform little on the actual state of the system, but instead could be a measure of the quantity of degassed material sitting at the head of the eruptive conduit. Indeed, the gas can be completely decoupled from the magma it fragments: the gas having been sourced from several kilometers of depth, but the fragments being sourced from just a few tens of meters of depth (Harris et al., 2013). Therefore, persistent degassing flux is a better parameter to track deeper magmatic processes that determine the volumes of gas involved in explosive events. The link documented here between persistent degassing and eruptive dynamics can provide indirect information of the real magmatic state of the volcano: emissions characterized by bomb-poor gas jetting are indicative of a highly charged system in gas. Instead, systems with high and sustained bomb masses likely imply the presence, and persistence, of a cap (e.g., Leduc et al., 2015; Gaudin et al., 2017; Simons et al., 2020).

## DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

## REFERENCES

Aiuppa, A., Federico, C., Giudice, G., Giuffrida, G., Guida, R., Gurrieri, S., et al. (2009). The 2007 Eruption of Stromboli Volcano: Insights from Real-Time Measurement of the Volcanic Gas Plume CO2/SO2 Ratio. *J. Volcanology Geothermal Res.* 182 (3-4), 221–230. doi:10.1016/j.jvolgeores.2008.09.013

Allard, P., Carbonnelle, J., Métrich, N., Loyer, H., and Zettwoog, P. (1994). Sulphur Output and Magma Degassing Budget of Stromboli Volcano. *Nature* 368 (6469), 326–330. doi:10.1038/368326a0

Allard, P. (1997). Endogenous Magma Degassing and Storage at Mount Etna. *Geophys. Res. Lett.* 24 (17), 2219–2222. doi:10.1029/97gl02101

Barberi, F. (1993). Volcanic hazard Assessment at Stromboli Based on Review of Historical Data. *Acta Vulcanologica*. 3, 173–187.

Barnie, T., Bombrun, M., Burton, M. R., Harris, A., and Sawyer, G. (2015). Quantification of Gas and Solid Emissions during Strombolian Explosions Using Simultaneous sulphur Dioxide and Infrared Camera Observations. *J. Volcanology Geothermal Res.* 300, 167–174. doi:10.1016/j.jvolgeores.2014.10.003

Bertolaso, G., De Bernardinis, B., Bosi, V., Cardaci, C., Ciolli, S., Colozza, R., et al. (2009). Civil protection Preparedness and Response to the 2007 Eruptive Crisis of Stromboli Volcano, Italy. *J. Volcanology Geothermal Res.* 182 (3-4), 269–277. doi:10.1016/j.jvolgeores.2009.01.022

Blackburn, E. A., Wilson, L., and Sparks, R. S. J. (1976). Mechanisms and Dynamics of Strombolian Activity. *J. Geol. Soc.* 132, 429–440. doi:10.1144/gsjgs.132.4.0429

Bogumil, K., Orphal, J., Homann, T., Voigt, S., Spietz, P., Fleischmann, O., et al. (2003). Measurements of Molecular Absorption Spectra with the SCIAMACHY Pre-flight Model: Instrument Characterization and Reference Data for

## AUTHOR CONTRIBUTIONS

ST performed the ash analysis, participated to the discussion, wrote the manuscript and created the figures. AH, LG, PB, TB, MB, and EM performed the field experiments, participated to the discussion and reviewed the manuscript.

## FUNDING

This project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement N731070 and this research was financed by the French government IDEX-ISITE initiative 16-IDEX-0001 (CAP 20-25). The ASHER machine was purchased by the ANR-STRAP project funded by the Agence Nationale de la Recherche (ANR-14-CE03-0004-04). This is Laboratory of Excellence ClerVolc contribution number 483.

## ACKNOWLEDGMENTS

We thank J-L. Devidal for its invaluable help with the EPMA at LMV. This is a contribution among some VOs and VRIs of WP4 and WP8 of EUROVOLC.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.669661/full#supplementary-material

Atmospheric Remote-Sensing in the 230–2380 Nm Region. *J. Photochem. Photobiol. A: Chem.* 157 (2-3), 167–184. doi:10.1016/s1010-6030(03)00062-5

Bombrun, M., Barra, V., and Harris, A. (2014). Algorithm for Particle Detection and Parameterization in High-Frame-Rate thermal Video. *J. Appl. Remote Sens* 8 (1), 083549. doi:10.1117/1.jrs.8.083549

Bombrun, M., Harris, A., Gurioli, L., Battaglia, J., and Barra, V. (2015). Anatomy of a Strombolian Eruption: Inferences from Particle Data Recorded with thermal Video. *J. Geophys. Res. Solid Earth* 120 (4), 2367–2387. doi:10.1002/2014jb011556

Branan, Y. K., Harris, A., Watson, I. M., Phillips, J. C., Horton, K., Williams-Jones, G., et al. (2008). Investigation of At-Vent Dynamics and Dilution Using thermal Infrared Radiometers at Masaya Volcano, Nicaragua. *J. Volcanology Geothermal Res.* 169 (1-2), 34–47. doi:10.1016/j.jvolgeores.2007.07.021

Braun, T., and Ripepe, M. (1993). Interaction of Seismic and Air Waves Recorded at Stromboli Volcano. *Geophys. Res. Lett.* 20 (1), 65–68. doi:10.1029/92gl02543

Burton, M. R., and Sawyer, G. M. (2016). iFit: An Intensity-Based Retrieval for SO2 and BrO from Scattered Sunlight Ultraviolet Volcanic Plume Absorption Spectra. *Atmos. Meas. Tech. Discuss.* (380), 1–47. doi:10.5194/amt-2015-380

Burton, M., Allard, P., Mure, F., and La Spina, A. (2007a). Magmatic Gas Composition Reveals the Source Depth of Slug-Driven Strombolian Explosive Activity. *Science* 317 (5835), 227–230. doi:10.1126/science.1141900

Burton, M. R., Mader, H. M., and Polacci, M. (2007b). The Role of Gas Percolation in Quiescent Degassing of Persistently Active Basaltic Volcanoes. *Earth Planet. Sci. Lett.* 264 (1-2), 46–60. doi:10.1016/j.epsl.2007.08.028

Burton, M. R., Caltabiano, T., Murè, F., Salerno, G., and Randazzo, D. (2009). SO2 Flux from Stromboli during the 2007 Eruption: Results from the FLAME Network and Traverse Measurements. *J. Volcanology Geothermal Res.* 182 (3-4), 214–220. doi:10.1016/j.jvolgeores.2008.11.025
Businger, S., Huff, R., Pattantyus, A., Horton, K., Sutton, A. J., Elias, T., et al. (2015). Observing and Forecasting Vog Dispersion from Kilauea Volcano, Hawaii. *Bull. Am. Meteorol. Soc.* 96 (10), 1667–1686. doi:10.1175/BAMS-D-14-00150.1

Calvari, S., Spampinato, L., Lodato, L., Harris, A. J. L., Patrick, M. R., Dehn, J., et al. (2005). Chronology and Complex Volcanic Processes during the 2002-2003 Flank Eruption at Stromboli Volcano (Italy) Reconstructed from Direct Observations and Surveys with a Handheld thermal Camera. *J. Geophys. Res.* 110 (B2), 1–23. doi:10.1029/2004jb003129

Capponi, A., James, M. R., and Lane, S. J. (2016). Gas Slug Ascent in a Stratified Magma: Implications of Flow Organisation and Instability for Strombolian Eruption Dynamics. *Earth Planet. Sci. Lett.* 435, 159–170. doi:10.1016/j.epsl. 2015.12.028

Caracciolo, A., Gurioli, L., Marianelli, P., Bernard, J., and Harris, A. (2021). Textural and Chemical Features of a “Soft” Plug Emitted during Strombolian Explosions: A Case Study from Stromboli Volcano. *Earth Planet. Sci. Lett.* 559, 116761. doi:10.1016/j.epsl.2021.116761

Chouet, B., Hamisevicz, N., and McGetchin, T. R. (1974). Photoballistics of Volcanic Jet Activity at Stromboli, Italy. *J. Geophys. Res.* 79 (32), 4961–4976. doi:10.1029/jb079i032p04961

Chouet, B., Saccorotti, G., Dawson, P., Martini, M., Scarpa, R., De Luca, G., et al. (1999). Broadband Measurements of the Sources of Explosions at Stromboli Volcano, Italy. *Geophys. Res. Lett.* 26 (13), 1937–1940. doi:10.1029/ 1999gl900400

Chouet, B., Dawson, P., Ohminato, T., Martini, M., Saccorotti, G., Giudicepietro, F., et al. (2003). Source Mechanisms of Explosions at Stromboli Volcano, Italy, Determined from Moment-Tensor Inversions of Very-Long-Period Data. *J. Geophys. Res.* 108 (B1), ESE 7–1–ESE 7–25. doi:10.1029/2002jb001919

D'Oriano, C., Bertagnini, A., and Pompilio, M. (2010). Ash Erupted during normal Activity at Stromboli (Aeolian Islands, Italy) Raises Questions on How the Feeding System Works. *Bull. Volcanology* 73 (5), 471–477. doi:10.1007/s00445- 010-0425-0

De Castro, E., and Morandi, C. (1987). Registration of Translated and Rotated Images Using Finite Fourier Transforms. *IEEE Trans. Pattern Anal. Mach. Intell.* 9 (5), 700–703. doi:10.1109/tpami.1987.4767966

Del Bello, E., Lane, S. J., James, M. R., Llewelin, E. W., Taddeucci, J., Scarlato, P., et al. (2015). Viscous Plugging Can Enhance and Modulate Explosivity of Strombolian Eruptions. *Earth Planet. Sci. Lett.* 423, 210–218. doi:10.1016/j.epsl. 2015.04.034

Delle Donne, D., Tamburello, G., Aiuppa, A., Bitetto, M., Lacanna, G., D'Aleo, R., et al. (2017). Exploring the Explosive-Effusive Transition Using Permanent Ultraviolet Cameras. *J. Geophys. Res. Solid Earth* 122 (6), 4377–4394. doi:10. 1002/2017jb014027

Francis, P., Oppenheimer, C., and Stevenson, D. (1993). Endogenous Growth of Persistently Active Volcanoes. *Nature* 366 (6455), 554–557. doi:10.1038/ 366554a0

Freret-Lorgeril, V., Donnadieu, F., Eychenne, J., Soriaux, C., and Latchimy, T. (2019). *In Situ* terminal Settling Velocity Measurements at Stromboli Volcano: Input from Physical Characterization of Ash. *J. Volcanology Geothermal Res.* 374, 62–79. doi:10.1016/j.jvolgeores.2019.02.005

Gaudin, D., Taddeucci, J., Scarlato, P., del Bello, E., Ricci, T., Orr, T., et al. (2017). Integrating Puffing and Explosions in a General Scheme for Strombolian-Style Activity. *J. Geophys. Res. Solid Earth* 122, 1860. doi:10.1002/2016jb013707

Gerst, A., Hort, M., Aster, R. C., Johnson, J. B., and Kyle, P. R. (2013). The First Second of Volcanic Eruptions from the Erebus Volcano Lava lake, Antarctica- Energies, Pressures, Seismology, and Infrasound. *J. Geophys. Res. Solid Earth* 118 (7), 3318–3340. doi:10.1002/jgrb.50234

Goto, A., Ripepe, M., and Lacanna, G. (2014). Wideband Acoustic Records of Explosive Volcanic Eruptions at Stromboli: New Insights on the Explosive Process and the Acoustic Source. *Geophys. Res. Lett.* 41 (11), 3851–3857. doi:10. 1002/2014gl060143

Gresta, S., Ripepe, M., Marchetti, E., D'Amico, S., Coltelli, M., Harris, A. J. L., et al. (2004). Seismoacoustic Measurements during the July-August 2001 Eruption of Mt. Etna Volcano, Italy. *J. Volcanology Geothermal Res.* 137 (1-3), 219–230. doi:10.1016/j.jvolgeores.2004.05.017

Gurioli, L., Harris, A. J. L., Houghton, B. F., Polacci, M., and Ripepe, M. (2008). Textural and Geophysical Characterization of Explosive Basaltic Activity at Villarrica Volcano. *J. Geophys. Res.* 113 (B8), 1–16. doi:10.1029/2007jb005328

Gurioli, L., Harris, A. J. L., Colò, L., Bernard, J., Favalli, M., Ripepe, M., et al. (2013). Classification, landing Distribution, and Associated Flight Parameters for a Bomb Field Emplaced during a Single Major Explosion at Stromboli, Italy. *Geology* 41 (5), 559–562. doi:10.1130/g33967.1

Gurioli, L., Colo', L., Bollasina, A. J., Harris, A. J. L., Whittington, A., and Ripepe, M. (2014). Dynamics of Strombolian Explosions: Inferences from Field and Laboratory Studies of Erupted Bombs from Stromboli Volcanofield and Labora-Tory Studies of Erupted Bombs from Stromboli Volcano. *J. Geophys. Res. Solid Earth* 119, 319–345. doi:10.1002/2013JB010355

Gurioli, L., Andronico, D., Bachelery, P., Balcone-Boissard, H., Battaglia, J., Boudon, G., et al. (2015). MeMoVolc Consensual Document: A Review of Cross-Disciplinary Approaches to Characterizing Small Explosive Magmatic Eruptions. *Bull. Volcanol* 77 (6), 1–33. doi:10.1007/s00445-015-0935-x

Harris, A., and Ripepe, M. (2007a). Synergy of Multiple Geophysical Approaches to Unravel Explosive Eruption Conduit and Source Dynamics - A Case Study from Stromboli. *Geochemistry* 67 (1), 1–35. doi:10.1016/j.chemer.2007.01.003

Harris, A., and Ripepe, M. (2007b). Temperature and Dynamics of Dégassing at Stromboli. *J. Geophys. Res.* 112 (B3), 1–18. doi:10.1029/2006jb004393

Harris, A. J. L., and Stevenson, D. S. (1997). Magma Budgets and Steady-State Activity of Volcano and Stromboli. *Geophys. Res. Lett.* 24 (9), 1043–1046. doi:10.1029/97gl00861

Harris, A. J. L., Ripepe, M., and Hughes, E. A. (2012). Detailed Analysis of Particle Launch Velocities, Size Distributions and Gas Densities during normal Explosions at Stromboli. *J. Volcanology Geothermal Res.* 231-232, 109–131. doi:10.1016/j.jvolgeores.2012.02.012

Harris, A. J. L., Delle Donne, D., Dehn, J., Ripepe, M., and Worden, A. K. (2013). Volcanic Plume and Bomb Field Masses from thermal Infrared Camera Imagery. *Earth Planet. Sci. Lett.* 365, 77–85. doi:10.1016/j.epsl.2013.01.004

Hort, M., and Seyfried, R. (1998). Volcanic Eruption Velocities Measured with a Micro Radar. *Geophys. Res. Lett.* 25 (1), 113–116. doi:10.1029/97gl03482

Houghton, B. F., Taddeucci, J., Andronico, D., Gonnermann, H. M., Pistolesi, M., Patrick, M. R., et al. (2016). Stronger or Longer: Discriminating between Hawaiian and Strombolian Eruption Styles. *Geology* 44, 163–166. doi:10. 1130/G37423.1

Jaupart, C., and Vergniolle, S. (1988). Laboratory Models of Hawaiian and Strombolian Eruptions. *Nature* 331 (6151), 58–60. doi:10.1038/331058a0

Jaupart, C., and Vergniolle, S. (1989). The Generation and Collapse of a Foam Layer at the Roof of a Basaltic Magma Chamber. *J. Fluid Mech.* 203, 347–380. doi:10.1017/s0022112089001497

Johnson, J., Aster, R., Jones, K. R., Kyle, P., and McIntosh, B. (2008). Acoustic Source Characterization of Impulsive Strombolian Eruptions from the Mount Erebus Lava lake. *J. Volcanology Geothermal Res.* 177 (3), 673–686. doi:10.1016/ j.jvolgeores.2008.06.028

Kuglin, C. D., and Hines, D. C. (1975). The Phase Correlation Image Alignment Method. *Proc. 1975 Int. Conf. Cybernetics Soc.*, 163–165.

Landi, P., Marchetti, E., La Felice, S., Ripepe, M., and Rosi, M. (2011). Integrated Petrochemical and Geophysical Data Reveals thermal Distribution of the Feeding Conduits from Stromboli Volcano, Italy. *Geophys. Res. Lett.* 38 (8), 1–6. doi:10.1029/2010GL046296

Lautze, N. C., and Houghton, B. F. (2005). Physical Mingling of Magma and Complex Eruption Dynamics in the Shallow Conduit at Stromboli Volcano, Italy. *Geology* 33, 425–428. doi:10.1130/G21325.1

Lautze, N. C., and Houghton, B. F. (2007). Linking Variable Explosion Style and Magma Textures during 2002 at Stromboli Volcano, Italy. *Bull. Volcanol.* 69, 445–460. doi:10.1007/s00445-006-0086-1

Lautze, N. C., and Houghton, B. F. (2008). Single Explosions at Stromboli in 2002: Use of Clast Microtextures to Map Physical Diversity across a Fragmentation Zone. *J. Volcanology Geothermal Res.* 170 (3-4), 262–268. doi:10.1016/j. jvolgeores.2007.10.011

Lautze, N. C., Taddeucci, J., Andronico, D., Cannata, C., Tornetta, L., Scarlato, P., et al. (2012). SEM-based Methods for the Analysis of Basaltic Ash from Weak Explosive Activity at Etna in 2006 and the 2007 Eruptive Crisis at Stromboli. *Phys. Chem. Earth, Parts A/B/C* 45-46, 113–127. doi:10.1016/j.pce.2011.02.001

Lautze, N., Taddeucci, J., Andronico, D., Houghton, B., Niemeijer, A., and Scarlato, P. (2013). "Insights into Explosion Dynamics and the Production of Ash at Stromboli from Samples Collected in Real-Time," in *Understanding Open-Vent Volcanism and Related Hazards*. Editors W. I. Rose, J. L. Palma, H. D. Granados,
and N. Varley, Boulder Colorado: Geological Society of America Special Papers, Boulder Colorado (498), 125-139. doi:10.1130/2013.2498(08)

Leduc, L., Gurioli, L., Harris, A., Colò, L., and Rose-Koga, E. F. (2015). Types and Mechanisms of Strombolian Explosions: Characterization of a Gas-Dominated Explosion at Stromboli. *Bull. Volcanol* 77 (1), 1-15. doi:10.1007/s00445-014-0888-5

Leibrandt, S., and Le Pennec, J.-L. (2015). Towards Fast and Routine Analyses of Volcanic Ash Morphometry for Eruption Surveillance Applications. *J. Volcanology Geothermal Res.* 297, 11-27. doi:10.1016/j.jvolgeores.2015.03.014

Mercalli, G. (1907). Vulcani Attivi Della Terra Morfologia-Dinamismo-Prodotti Distribuzione Geografica-Cause. *Nature* 76, 291. doi:10.1038/076291a0

Metrich, N. (2005). Triggering Mechanism at the Origin of Paroxysms at Stromboli (Aeolian Archipelago, Italy): The 5 April 2003 Eruption. *Geophys. Res. Lett.* 32 (10), 1-4. doi:10.1029/2004gl022257

Nave, R., Isaia, R., Vilardo, G., and Barclay, J. (2010). Re-assessing Volcanic hazard Maps for Improving Volcanic Risk Communication: Application to Stromboli Island, Italy. *J. Maps* 6 (1), 260-269. doi:10.4113/jom.2010.1061

Ntepe, N., and Dorel, J. (1990). Observations of Seismic Volcanic Signals at Stromboli Volcano (Italy). *J. Volcanology Geothermal Res.* 43 (1-4), 235-251. doi:10.1016/0377-0273(90)90054-j

Oppenheimer, J., Capponi, A., Cashman, K. V., Lane, S. J., Rust, A. C., and James, M. R. (2020). Analogue Experiments on the Rise of Large Bubbles through a Solids-Rich Suspension: A “weak Plug” Model for Strombolian Eruptions. *Earth Planet. Sci. Lett.* 531, 115931. doi:10.1016/j.epsl.2019.115931

Parfitt, E. A., and Wilson, L. (1995). Explosive Volcanic Eruptions-IX. The Transition between Hawaiian-style Lava Fountaining and Strombolian Explosive Activity. *Geophys. J. Int.* 121 (1), 226-232. doi:10.1111/j.1365-246x.1995.tb03523.x

Patrick, M. R., Harris, A. J. L., Ripepe, M., Dehn, J., Rothery, D. A., and Calvari, S. (2007). Strombolian Explosive Styles and Source Conditions: Insights from thermal (FLIR) Video. *Bull. Volcanol* 69 (7), 769-784. doi:10.1007/s00445-006-0107-0

Pering, T. D., McGonigle, A. J. S., James, M. R., Tamburello, G., Aiuppa, A., Delle Donne, D., et al. (2016). Conduit Dynamics and post Explosion Degassing on Stromboli: A Combined UV Camera and Numerical Modeling Treatment. *Geophys. Res. Lett.* 43, 5009-5016. doi:10.1002/2016GL069001

Pering, T. D., Liu, E. J., Wood, K., Wilkes, T. C., Aiuppa, A., Tamburello, G., et al. (2020). Combined Ground and Aerial Measurements Resolve Vent-specific Gas Fluxes from a Multi-Vent Volcano. *Nat. Commun.* 11, 3039. doi:10.1038/s41467-020-16862-w

Pioli, L., and Harris, A. J. L. (2019). Real-Time Geophysical Monitoring of Particle Size Distribution During Volcanic Explosions at Stromboli Volcano (Italy). *Front. Earth Sci.* 7, 52, 1-13. doi:10.3389/feart.2019.00052

Pioli, L., Pistolesi, M., and Rosi, M. (2014). Transient Explosions at Open-Vent Volcanoes: The Case of Stromboli (Italy). *Geology* 42 (10), 863-866. doi:10.1130/g35844.1

Platt, U., and Stutz, J. (2008). Differential Optical Absorption Spectroscopy. *Phys. Earth Space Environments* 135-174. doi:10.1007/978-3-540-75776-4

Ripepe, M., and Gordeev, E. (1999). Gas Bubble Dynamics Model for Shallow Volcanic Tremor at Stromboli. *J. Geophys. Res.* 104 (B5), 10639-10654. doi:10.1029/98jb02734

Ripepe, M., Rossi, M., and Saccorotti, G. (1993). Image Processing of Explosive Activity at Stromboli. *J. Volcanology Geothermal Res.* 54 (3-4), 335-351. doi:10.1016/0377-0273(93)90071-x

Ripepe, M., Ciliberto, S., and Della Schiava, M. (2001). Time Constraints for Modeling Source Dynamics of Volcanic Explosions at Stromboli. *J. Geophys. Res.* 106 (B5), 8713-8727. doi:10.1029/2000jb900374

Ripepe, M., Harris, A. J. L., and Carniel, R. (2002). Thermal, Seismic and Infrasonic Evidences of Variable Degassing Rates at Stromboli Volcano. *J. Volcanology Geothermal Res.* 118 (3-4), 285-297. doi:10.1016/s0377-0273(02)00298-6

Ripepe, M., Marchetti, E., Ulivieri, G., Harris, A., Dehn, J., Burton, M., et al. (2005). Effusive to Explosive Transition during the 2003 Eruption of Stromboli Volcano. *Geology* 33 (5), 341. doi:10.1130/g21173.1

Ripepe, M., Marchetti, E., and Ulivieri, G. (2007). Infrasonic Monitoring at Stromboli Volcano during the 2003 Effusive Eruption: Insights on the Explosive and Degassing Process of an Open Conduit System. *J. Geophys. Res.* 112 (B9), 1-13. doi:10.1029/2006jb004613

Ripepe, M., Donne, D. D., Harris, A., Marchetti, E., and Ulivieri, G. (2008). Dynamics of Strombolian Activity. *Geophys. Monogr. Ser.* (182), 39-48. doi:10.1029/182gm05

Ripepe, M., Delle Donne, D., Lacanna, G., Marchetti, E., and Ulivieri, G. (2009). The Onset of the 2007 Stromboli Effusive Eruption Recorded by an Integrated Geophysical Network. *J. Volcanology Geothermal Res.* 182 (3-4), 131-136. doi:10.1016/j.jvolgeores.2009.02.011

Rosi, M., Bertagnini, A., and Landi, P. (2000). Onset of the Persistent Activity at Stromboli Volcano (Italy). *Bull. Volcanology* 62 (4-5), 294-300. doi:10.1007/s004450000098

Rosi, M., Pistolesi, M., Bertagnini, A., Landi, P., Pompilio, M., and Di Roberto, A. (2013). Chapter 14 Stromboli Volcano, Aeolian Islands (Italy): Present Eruptive Activity and Hazards. *Geol. Soc. Lond. Mem.* 37 (1), 473-490. doi:10.1144/m37.14

Schmid, R. (1981). Descriptive Nomenclature and Classification of Pyroclastic Deposits and Fragments. *Geol. Rundsch* 70 (2), 794-799. doi:10.1007/bf01822152

Simons, B. C., Cronin, S. J., Eccles, J. D., Bebbington, M. S., and Jolly, A. D. (2020a). Spatiotemporal Variations in Eruption Style, Magnitude and Vent Morphology at Yasur Volcano, Vanuatu: Insights into the Conduit System. *Bull. Volcanol* 82 (8), 1-22. doi:10.1007/s00445-020-01394-4

Simons, B. C., Cronin, S. J., Eccles, J. D., Jolly, A. D., Garaebiti, E., and Cevuard, S. (2020b). Spatiotemporal Variations in Eruption Style and Magnitude at Yasur Volcano, Vanuatu: Part 2. Extending Strombolian Eruption Classifications. *Bull. Volcanology* 82:70 (11). doi:10.1007/s00445-020-01404-5

Slawson, P. R., and Csanady, G. T. (1967). On the Mean Path of Buoyant, Bent-Over Chimney Plumes. *J. Fluid Mech.* 28 (02), 311. doi:10.1017/s0022112067002095

Smekens, J.-F., Burton, M. R., and Clarke, A. B. (2015). Validation of the SO 2 Camera for High Temporal and Spatial Resolution Monitoring of SO 2 Emissions. *J. Volcanology Geothermal Res.* 300, 37-47. doi:10.1016/j.jvolgeores.2014.10.014

Suckale, J., Keller, T., Cashman, K. V., and Persson, P.-O. (2016). Flow-to-fracture Transition in a Volcanic Mush Plug May Govern normal Eruptions at Stromboli. *Geophys. Res. Lett.* 43 (12), 12071-12081. doi:10.1002/2016GL071501

Taddeucci, J., Scarlato, P., Capponi, A., Del Bello, E., Cimarelli, C., Palladino, D. M., et al. (2012). High-speed Imaging of Strombolian Explosions: The Ejection Velocity of Pyroclasts. *Geophys. Res. Lett.* 39 (2), 1-6. doi:10.1029/2011gl050404

Taddeucci, J., Edmonds, M., Houghton, B., James, M. R., and Vergniolle, S. (2015). "Hawaiian and Strombolian Eruptions," in *The Encyclopedia of Volcanoes*. Editors H. Sigurdsson. Amsterdam: Academic Press, 485-503. doi:10.1016/B978-0-12-385938-9.00027-4

Tamburello, G., Aiuppa, A., Kantzas, E. P., McGonigle, A. J. S., and Ripepe, M. (2012). Passive vs. Active Degassing Modes at an Open-Vent Volcano (Stromboli, Italy). *Earth Planet. Sci. Lett.* 359-360, 106-116. doi:10.1016/j.epsl.2012.09.050

Thivet, S., Gurioli, L., Di Muro, A., Derrien, A., Ferrazzini, V., Gouhier, M., et al. (2020a). Evidences of Plug Pressurization Enhancing Magma Fragmentation During the September 2016 Basaltic Eruption at Piton de la Fournaise (La Réunion Island, France). *Geochem. Geophys. Geosyst.* 21, e2019GC008611. doi:10.1029/2019GC008611

Thivet, S., Gurioli, L., Di Muro, A., Eychenne, J., Besson, P., and Nedelec, J.-M. (2020b). Variability of ash deposits at Piton de la Fournaise (La Reunion Island): insights into fragmentation processes at basaltic shield volcanoes. *Bull. Volcanol* 82 (9), 1-20. doi:10.1007/s00445-020-01398-0

Turner, J. S. (1962). The 'starting Plume' in Neutral Surroundings. *J. Fluid Mech.* 13 (03), 356-368. doi:10.1017/s0022112062000762

Turner, J. S. (1969). Buoyant Plumes and Thermals. *Annu. Rev. Fluid Mech.* 1 (1), 29-44. doi:10.1146/annurev.fl.01.010169.000333

Valade, S., Lacanna, G., Coppola, D., Laiolo, M., Pistolesi, M., Donne, D. D., et al. (2016). Tracking Dynamics of Magma Migration in Open-Conduit Systems. *Bull. Volcanol* 78 (11), 1-12. doi:10.1007/s00445-016-1072-x

Vergniolle, S., and Brandeis, G. (1994). Origin of the Sound Generated by Strombolian Explosions. *Geophys. Res. Lett.* 21 (18), 1959-1962. doi:10.1029/94gl01286
Vergniolle, S., and Brandeis, G. (1996). Strombolian Explosions: 1. A Large Bubble Breaking at the Surface of a Lava Column as a Source of Sound. *J. Geophys. Res.* 101 (B9), 20433–20447. doi:10.1029/96jb01178

Vergniolle, S., Boichu, M., and Caplan-Auerbach, J. (2004). Acoustic Measurements of the 1999 Basaltic Eruption of Shishaldin Volcano, Alaska. *J. Volcanology Geothermal Res.* 137 (1-3), 109–134. doi:10.1016/j.jvolgeores.2004.05.003

Voigt, S., Orphal, J., Bogumil, K., and Burrows, J. P. (2001). The Temperature Dependence (203-293 K) of the Absorption Cross Sections of O3 in the 230-850 Nm Region Measured by Fourier-Transform Spectroscopy. *J. Photochem. Photobiol. A: Chem.* 143 (1), 1–9. doi:10.1016/s1010-6030(01)00480-4

Washington, H. S. (1917). Persistence of Vents at Stromboli and its Bearing on Volcanic Mechanism. *Geol. Soc. America Bull.* 28 (1), 249–278. doi:10.1130/gsab-28-249

Weill, A., Brandeis, G., Vergniolle, S., Baudin, F., Bilbille, J., Fèvre, J.-F., et al. (1992). Acoustic Sounder Measurements of the Vertical Velocity of Volcanic Jets at Stromboli Volcano. *Geophys. Res. Lett.* 19 (23), 2357–2360. doi:10.1029/92gl02502

Williams-Jones, G., Horton, K. A., Elias, T., Garbeil, H., Mouginis-Mark, P. J., Sutton, A. J., et al. (2006). Accurately Measuring Volcanic Plume Velocity with Multiple UV Spectrometers. *Bull. Volcanol* 68 (4), 328–332. doi:10.1007/s00445-005-0013-x

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Thivet, Harris, Gurioli, Bani, Barnie, Bombrun and Marchetti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.