Abstract
Magmatic degassing, typically measured as SO2 flux, plays a fundamental role in controlling volcanic eruption style and is one of the key parameters used by volcano observatories to assess volcanic unrest and detect eruption precursors. Volcanic tremor, the integrated amplitude of seismic energy release over a range of frequencies, is also a key parameter in volcano monitoring. A connection between volcanic degassing and tremor has been inferred through correlations between the signals which are often, but not always, observed during periods of unrest or eruption. However, data are often equivocal and our understanding of the physical processes, which couple degassing with tremor are still evolving. New insights into degassing-tremor coupling can be made by investigation of the long-term relationship between degassing and tremor, focusing on the frequency-dependence of tremor and passive degassing behavior. In this study, we examine how long-term SO2 emission rates and volcanic tremor on Mt. Etna, track rapid variability in eruptive dynamics. Correlations between SO2 flux and tremor are explored in both quiescent and eruptive periods, comparing the two parameters at both long and short time-scales (< < 1 day) for ∼2 years. Our analysis reveals that over ∼month-long timescales passive degassing of SO2 and tremor tend to be well-correlated, but these correlations are lost over shorter timescales. This reflects a coupling process between passive degassing and tremor, produced by a combination of gas flow through permeable magma and the convective flow of magma within the conduit. Short-term correlations are lost because variations in the continuous degassing process are relatively small compared with the overall degassing rate and fall below measurement noise. During eruptive periods strong correlations are observed between degassing and tremor, with a significant contribution of higher frequency signal in tremor, controlled by eruptive style. These observations suggest that in syn-eruptive periods the tremor source is dominated by the coupling between the eruption column and the ground through infrasonic waves, rather than conduit processes. Our results demonstrate the importance of high quality long-term observations and offer new insights into the physical mechanisms which couple degassing and volcanic tremor at active volcanoes.
Introduction
Over the last decades, technological advances have allowed volcanic activity to be monitored at ever-increasing spatial and temporal resolutions (e.g., Heliker et al., 2003; Calvari et al., 2008; Johnson and Poland, 2013). Particularly, in the case of volcanic ground-based gas measurements the development of automated networks of spectrometer gas sensors (e.g., Edmonds et al., 2003a; Salerno et al., 2009b), ultraviolet and thermal cameras (e.g., Mori and Burton, 2006; Burton et al., 2015a; Lopez et al., 2015), and FTIR and multigas sensors (e.g., Burton et al., 2003; Shinohara, 2005; Taquet et al., 2017), have improved the temporal resolution of magmatic gas composition and flux observations, e.g., SO2 flux, from typically of order hours to ∼1 Hz. This has permitted comparison with geophysical measurements, thus allowing better integration of both geochemical and geophysical parameters for refining models and to identify eruptive anomalies and unrest (e.g., Aiuppa et al., 2010; Bonaccorso et al., 2011b; Poland et al., 2012; Patanè et al., 2013; Burton et al., 2015b; Hibert et al., 2015; Nadeau et al., 2015). In particular, volcanic SO2 emissions are important indicators of subsurface processes, and the study of their temporal evolution provides inferences on processes occurring at shallow depth (∼4–5 km from crater top). Measurements of SO2 emissions made over almost 40 years (e.g., Williams-Jones et al., 2008), have shown remarkable transitions from quiescence to unrest at both silicic magmatic systems (e.g., Fischer et al., 1994; Williams-Jones et al., 2001; Nadeau et al., 2011) and mafic volcanoes (e.g., Malinconico, 1979; Voight et al., 1999; Sutton et al., 2001; Caltabiano et al., 2004; Kazahaya et al., 2004). Measurement of SO2 outgassing also provide constraints on magma-degassing budgets and mass balance (e.g., Wallace and Gerlach, 1994; Allard, 1997; Shinohara, 2008; Steffke et al., 2011).
Volcanic tremor is observed at volcanoes as background seismic radiation in quiescent stages and as peaks in amplitudes during eruptive episodes such as explosive eruptions (e.g., Alparone et al., 2003; Patanè et al., 2013). Dominant frequency ranges between 0.1 and 10 Hz and episodes of high amplitude tremor may persist for months (e.g., Kubotera, 1974; McNutt, 1992; Zobin, 2003). Physical processes generating volcanic tremor are thought to be associated with unsteady mass transport-flow of magma dynamically coupled with the surrounding rocks (Steinberg and Steinberg, 1975; Schick and Mugiono, 1991; Neuberg and Pointer, 2000; Battaglia et al., 2005). However, several other mechanisms (e.g., Gordeev, 1993; Benoit and McNutt, 1997) and models have been proposed as potential sources of tremor (e.g., Aki et al., 1977; Chouet et al., 1987) depending on individual volcanoes and eruption style (Konstantinou and Schlindwein, 2002; Matoza and Fee, 2014). In particular, a review of the engineering literature associated with studies of two phase fluid flow induced vibration in pipes was recently published (Miwa et al., 2015), focusing on the hydrodynamic force produced by flows that generate potentially destructive vibrations in industrial machines and infrastructure. The review highlighted the impact of different flow patterns, turbulence and pipe geometry, producing a fluctuation force magnitude spectrum with a frequency range similar to that observed for tremor in volcanic settings when velocities were in the range of magma flow during quiescent degassing ∼0.6 ms−1 (Burton et al., 2007).
Persistent degassing from active volcanoes is widely associated with volcanic tremor (e.g., Williams-Jones et al., 2001; Konstantinou and Schlindwein, 2002; McNutt, 2002). Evidence of coupling between the two parameters has been discussed at different volcanic systems (e.g., Mt. Etna: Gresta et al., 1991, Bruno et al., 1995; Patanè et al., 2013; Zuccarello et al., 2013; Soufrière Hills: Miller et al., 1998; Edmonds et al., 2003b; Piton de la Fournaise: Battaglia et al., 2005; Colima: Vargas-Bracamontes et al., 2009). Increases in tremor amplitude were observed prior to and during effusive eruptions and synchronous with explosive activity (e.g., Usu: Omori, 1911; Pavlof: McNutt, 1986; Hekla: Brandsdóttir and Einarsson, 1992; Galeras: Fischer et al., 1994; Mt. Etna: Cannata et al., 2008; Bonaccorso et al., 2011b; Kilauea: Nadeau et al., 2015). Changes in the dominant tremor frequency have been commonly associated with changes in the regime and style of eruptive activity (e.g., Ereditato and Luongo, 1994; Thompson et al., 2002; Alparone et al., 2003; Bryan and Sherburn, 2003; Cannata et al., 2018). At Mt. Etna, Leonardi et al. (2000a) studied the relationship between volcanic tremor and SO2 flux in the period between 1987 and 1992 by cross-correlation analysis. Their results indicated that, in the case of eruptive activity, the two signals strongly correlated. The authors proposed that such as behavior might have resulted from common physical mechanisms related to magma dynamics. Similarly, at Soufrière Hills, Montserrat, Young et al. (2003) found systematic and direct correlation between SO2 flux and tremor analyzing data between December 1999 and January 2000, with gas flux lagging behind the tremor signal. Similar correlations between SO2 outgassing and seismic amplitude were observed also at Villarrica (Palma et al., 2008), Yasur (Bani and Lardy, 2007), Fuego (Nadeau et al., 2011), and Kilauea (Nadeau et al., 2015). Volcanic tremor has also shown good correlation with other geochemical parameters. For instance, Alparone et al. (2005) carried out a statistical analysis on radon emissions from the soil and reduced displacement of volcanic tremor both recorded during paroxysmal explosive phases of Mt. Etna’s summit craters. Their studies revealed increase in radon concentrations ∼58 ± 12 h prior to changes in volcanic tremor. More recently, investigating the relationship between CO2 flux and volcanic tremor at Mt. Etna, Cannata et al. (2009a) found that variations in the geochemical signal preceded those of volcanic tremor of ∼50 days. Nevertheless, in some cases contradictory or/and lack of relationship between eruptive activity and tremor amplitude has been observed. Doukas and Gerlach (1995), observed an episode of inverse correlation between SO2 flux and volcanic tremor at Mount Spurr, Alaska during the 1991–1993 eruptive activity, interpreting the gas declined as a result of SO2 absorption by the hydrothermal system. Similarly, inverse correlations between CO2 flux and volcanic tremor were observed at Stromboli (Aiuppa et al., 2009). At Soufrière hills volcano, Watson et al. (2000) reported that high rates of SO2 were associated with enhanced seismicity and ground deformation a month before the 1997 dome collapse. However, previous observations by Young et al. (1998) showed that the enhanced long-period seismicity at Soufrière Hills did not relate to increases in SO2 flux. This indicates that, though tremor and degassing are somewhat coupled, the nature of their relationship, as well as the source mechanism for tremor, are still poorly understood. It is likely that several processes are involved in the generation of volcanic tremor from gas and magma flows, making the unraveling of the tremor and degassing relationship challenging.
Here, we focus on three mechanisms by which volcanic tremor may be coupled with degassing: (1) flow of gas through permeable magma (e.g., Burton et al., 2007; La Spina et al., 2017), (2) magma flow within a conduit (e.g., Kazahaya et al., 1994; Beckett et al., 2014), and (3) coupling of eruptive processes to ground seismicity through infrasound during explosive activity (Matoza and Fee, 2014). The links between these processes and magmatic degassing are explored on Mt. Etna by comparing seismic tremor measured with the INGV seismic network with SO2 flux measurements collected with INGV FLAME network. Long- and short-timescale comparisons of the two parameters for six case studies selected between 2007 and 2008 were carried out by correlation analysis to investigate whether the correlation holds across different timescales.
Eruptive Activity Between 2007 and 2008
Mt. Etna is the most active volcano in Europe characterized by extensive quiescent and active degassing activity that occurs at the main craters (North-East Crater: NEC, South-East Crater: SEC, and central craters Voragine: VOR, and Bocca Nuova: BN, Figure 1; e.g., Allard, 1997; Aiuppa et al., 2008). Between 2007 and 2008, eruptive activity at Mt. Etna resumed after a quiescent period following the 2006 eruption (Bonaccorso et al., 2011a). Activity was vigorous and characterized by a series of relatively short explosive episodes and long-lasting lava effusion both occurring from the eastern flank of SEC. (Andronico et al., 2008; Corsaro and Miraglia, 2009; Bonaccorso et al., 2011a). This period is divided into two main phases, firstly the period January 2007 to early May 2008 characterized by intermittent lava fountains and the period May–December 2008, characterized by effusive activity (Figure 2). The first stage, from January 2007 to early May 2008, consisted of sporadic violent strombolian and lava fountaining episodes accompanied by short-lasting lava effusion (e.g., Andronico et al., 2008; Di Grazia et al., 2009; Bonaccorso et al., 2011a,b; Behncke et al., 2016). Paroxysms exhibited recurrent features consisting on increasing strombolian activity at first, lava flow output, and transition from strombolian to lava fountaining (Aloisi et al., 2009; Di Grazia et al., 2009; Langer et al., 2010). This activity was similar to the lava fountain sequences observed at SEC in 2000 (e.g., Alparone et al., 2003; Allard et al., 2005; La Spina et al., 2015), and between 2011 and 2015 (e.g., Calvari et al., 2011; Patanè et al., 2013; Corsaro et al., 2017). The effusion phase, between 13 May and December 2008, started 3 days after the 10 May paroxysm (Aloisi et al., 2009; James et al., 2011), and was dominated by lava effusion occasionally accompanied between May and early September 2008 with strombolian activity from the eruptive fissure (Cannata et al., 2009b; Bonaccorso et al., 2011a; Currenti et al., 2011).
FIGURE 1
FIGURE 2
Data Acquisition and Analysis
All data presented here were acquired by the continuous geochemical and geophysical monitoring system of Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo (INGV-OE).
SO2 Flux
The bulk SO2 flux from the summit craters and eruptive fissure of Mt. Etna was measured automatically by the FLAME scanning spectrometer network (Salerno et al., 2009b; Calvari et al., 2011). The network consists of ten ultraviolet scanning spectrometer stations spaced ∼7 km apart and installed at an altitude of ∼900 m above sea level (a.s.l.) on the flanks of Mt. Etna (Calvari et al., 2011). Recorded UV spectra were retrieved in SO2 column amounts using the DOAS method (e.g., Platt and Stutz, 2008), and applying a modeled background reference spectrum (Burton et al., 2009; Salerno et al., 2009a). Retrieved SO2 data were transmitted to INGV–OE where they were converted into mass flux rate. Each flux datum was time corrected to account for the travel time of the plume from the summit craters to the scanning plane of each scanner of the network (∼14 km). Uncertainty in computed SO2 flux depends to an extent on plume velocity, because at very low velocities the absolute velocity error becomes a larger proportion of the relative error in calculated flux. Assumptions on plume height are made based on observations of the relationship between plume velocity and height, and light scattering effects also contribute to the flux error budget (e.g., Mori et al., 2006; Campion et al., 2015). Salerno et al. (2009b) estimates uncertainty in SO2 flux by stationary automatic scanning array between −22 + 36%, but error can vary to greater than 100% depending on conditions.
Seismic Data
The permanent seismic network of INGV comprises 45 three-component broadband (40 s) digital stations with continuous data acquisition and transmission at a sampling frequency of 100 Hz (Figure 1; e.g., Di Grazia et al., 2006; Patanè et al., 2008). In this work, we use data from the Etna Cratere del Piano (ECPN) station, which is set up on the southern flank of the volcano at an altitude of 2900 m a.s.l. and ∼1 km from the summit craters (Figure 1). The ECPN station was chosen as it offers the longest data continuity provided during the studied period, the best signal to noise ratio and its proximity to the summit craters makes it especially sensitive to eruptive activity. Volcanic tremor spectral amplitude was calculated using the root mean square (RMS expressed in arbitrary unit) of the seismic signal recorded on the vertical component. In this study, the daily RMS at overall spectral amplitude (OSA) was used for long-term characterization of volcanic tremor and to explore its relationship with daily SO2 flux. In order to inspect the short-term relationship between SO2 flux and tremor and identify if any and which frequency components correlated most strongly with degassing rates, seismic amplitude was decomposed into separate frequency bands (e.g., Gresta et al., 1987; Thompson et al., 2002; Di Grazia et al., 2009). RMS spectral amplitude was calculated for several frequency bands of 1 Hz width in the frequency range 0.025–10.5 Hz, and tremor amplitude was obtained by averaging the RMS values with a 5 min time window. This time window was chosen according Cannata et al. (2009a), who report that time windows of this length are a good compromise between the stability of the signal and the opportunity of observing details of the evolution of the parameter during eruptions. Uncertainties affecting RMS depends on the dispersion of the averaged RMS time window; it ranges between 10% for the short-5 min intraday study and up to mean 25%, enveloped between a minimum of 6% and maximum of 77%, for the long-temporal daily averaged RMS.
Data Analysis
To investigate the relationship between the SO2 flux and volcanic tremor, we apply the correlation analysis, which expresses the strength of linkage or co-occurrence between two variables in a single value ranging from −1 to +1, i.e.,−1 ≤ r ≤ +1 (e.g., Davis, 1986; McKillup and Dyar, 2010). Usually for evaluating dependences between two parameters, the conventional Pearson’s correlation analysis is used. However, this method requires normal distribution of the parameter samples. Since both SO2 flux and volcanic tremor data were characterized by non-Gaussian distribution (mean skewness and kurtosis are 1.1 and 0.9, and 1.2 and 1.1, respectively; Table 1 and Supplementary Figure S1) and the SO2 flux sample sizes were small (maximum 120 daily-light observations), the non-parametric Spearman’s Rank correlation analysis was applied (e.g., Zar, 1972; Davis, 1986; Swan and Sandilands, 1995). Compared to the Pearson’s correlation analysis, the Spearman’s correlation method does not require continuous-level data (interval or ratio), as it uses ranks instead of assumptions on the distributions of two variables. This allows analysis of the association between variables of ordinal measurement levels. Mathematically, Spearman’s and Pearson’s correlations are very similar in the way that they use different measurements to calculate the strength of association of two parameters. Pearson’s correlation applies standard deviations, while Spearman’s the difference in ranks (Davis, 1986; Swan and Sandilands, 1995).
Table 1
| Phase | Parameter | Min | Max | Mean | δ | |
|---|---|---|---|---|---|---|
| Intermittent paroxysmal | ΦS02 | 400 | 13,500 | 2400 | 1500 | |
| Vt | −6.0 | −4.6 | −5.5 | 0.2 | ||
| Eruption 2008–2009 | ΦS02 | 450 | 20,000 | 2700 | 650 | |
| Vt | −5.7 | -4.5 | -5.4 | 0.2 | ||
| Case study | 1 | 2 | 3 | 4 | 5 | 6 |
| Date | April 29, 2007 | May 07, 2007 | September 04, 2007 | July 07, 2008 | July 10, 2008 | July 24, 2008 |
| Eruptive style | explosive | explosive | explosive | effusive – explosive | effusive – mild explosive | effusive-mild explosive |
| Frequency band Vt (Hz) | 5.5–6.5 | 5.5–6.5 | 2.5–3.5 | 0.025–0.5 | 1.5–2.5 | 0.5–1.5 |
| ρ | 0.0.7 | 0.9 | 0.7 | 0.6 | 0.8 | 0.6 |
| Min ΦS02 | 700 | 1300 | 1300 | 2100 | 1100 | 750 |
| Max ΦS02 | 20,000 | 13,000 | 10,000 | 11,000 | 17,000 | 10,300 |
| Mean ΦS02 | 5600 | 6000 | 3500 | 5500 | 7500 | 2600 |
| δ ΦSO2 | 4500 | 3700 | 1500 | 2500 | 4100 | 1950 |
| Min Vt | 4.0 × 10−7 | 7.2 × 10−7 | 4.2 × 10−6 | 1.9 × 10−7 | 1.4 × 10−6 | 8.8 × 10−7 |
| Max Vt | 1.8 × 10−5 | 2.8 × 10−5 | 1.5 × 10−5 | 3.3 × 10−7 | 3.6 × 10−6 | 1.9 × 10−6 |
| Mean Vt | 5.3 × 10−6 | 9.2 × 10−6 | 7.1 × 10−6 | 2.4 × 10−7 | 2.1 × 10−6 | 1.1 × 10−6 |
| d Vt | 3.0 × 10−6 | 1.1 × 10−5 | 2.5 × 10−6 | 2.6 × 10−8 | 4.9 × 10−7 | 1.8 × 10−7 |
| ΦSO2 = SO2 flux in t/d; Vt = Volcanic tremor amplitude (arbitrary units); δ = Standard deviation; ρ = Coefficient of correlation | ||||||
Details of the main features and parameters of the intermittent paroxysmal and eruptive 2008–2009 phases and of the six case studies investigated.
Results of the tremor frequency bands with maximum correlation used in this study together with statistical details on the SO2 emission rate and volcanic tremor.
Results
Figure 2 reports the weekly averaged SO2 flux and the daily OSA of the volcanic tremor RMS for long-term observation between 2007 and 2008. Over the investigated period, tremor and SO2 flux show common changes at different temporal and magnitude scales associated in both stages of quiescent-passive degassing and eruptive activity (Figure 2). During the intermittent paroxysmal phase, preceding the opening of the eruptive fissure on May 13, 2008, the SO2 rates and tremor behaved in a similar manner showing marked oscillations in correspondence to the explosive activity. These waxing-waning trends have higher amplitudes starting from the end of July 2007 until the opening of the 2008–2009 eruptive fissure (Table 1). Note that in Figure 2, volcanic tremor RMS is plotted in logarithmic scale to allow for better comparison with SO2 flux, and values are expressed in arbitrary units. Short-term intraday comparison between SO2 flux-volcanic tremor patterns were carried out to explore potential correlation between the two signals on a scale of minutes (Figure 3). Of the seven lava fountain episodes of the intermittent paroxysmal phase and of the strombolian activity fed by the eruptive fissure on 14 May, only 6 days could be analyzed (Table 1). The other episodes had to be excluded due to eruptive events at night (when UV spectroscopy was not possible) or on days when the plume was blown to the northern and western flanks of the volcano beyond the coverage of the FLAME network. In addition, some of the eruptive episodes were also excluded from the analysis due to the limited number of SO2 scans on the given days.
FIGURE 3
Figure 3 shows the intraday SO2 flux and volcanic tremor data for the six case studies. In each of the six graphs, SO2 flux is plotted together with RMS at a dominant frequency, i.e., that has shown the highest SO2 flux-volcanic tremor correlation coefficient between 0.025 and 10.5 Hz; likewise correlation analysis was also performed considering the OSA (Table 1; Figure 4). Signals display simultaneous changes in both magnitude and temporal scale repeatedly over the course of the measurements and persist over time scales ranging from 6 to 12 h. This correlation occurs during both stage of onset and waning phase of strombolian activity (case 1 and 2) and transition from strombolian to lava fountaining (case 3), both eruptive styles superimposed on short-lived and persistent lava flow output (4–6; Table 1). SO2 flux was strongly correlated with seismic tremor between 2.5 and 6.5 Hz during explosive eruptions and between 0.025 and 2.5 Hz (case studies 1–3, and in 4–6, respectively; Table 1), isolated anticorrelation were also observed for cases 1, 2, and 4, when SO2 emission rates increased while RMS decreased.
FIGURE 4
Discussion
From these observations, we see that in both the initial period between January 2007–May 2008 and the effusive eruption from May 2008 onwards, a strong correlation is observed between tremor and the weekly averaged SO2 flux (Figure 2). We highlight in particular the quiet eruptive period between May and December 2007, when the relationship between tremor and SO2 flux is marked, and a further period of close correlation is observed between May and December 2008, associated with the effusive eruption. We propose that this arises from a driving mechanism of magma flow, in which seismic energy is produced through friction between the flowing magma and conduit walls. Analysis of the source depth of seismic tremor demonstrates that the entire shallow conduit is a tremor source (e.g., Patanè et al., 2008; Cannata et al., 2013), consistent with our frictional flow hypothesis, with a peak in energy at a depth of 2–3 km, perhaps due to the lower viscosity of magma at this depth permitting a maximum ascent rate (e.g., Burton et al., 2007). There are two main regimes of magma flow, firstly as a result of effusive eruption, in which magma degasses and crystallizes during ascent before erupting, and secondly due to the convective overturn of magma during quiescent degassing (Kazahaya et al., 1994; Beckett et al., 2014) which is required to sustain the persistent degassing flux of Mt. Etna (e.g., Allard, 1997). This flow transports magma from below 3 km depth (from the volcano summit) to near-surface, which is the principle pressure regime where sulfur and water exsolution takes place at Mt. Etna (Metrich et al., 2004; Spilliaert et al., 2006; Wallace and Edmonds, 2011). Convective overturn leads to permanent burial of degassed magma in the roots of the plumbing system, revealed as a voluminous plutonic body in seismic tomography observations (Patanè et al., 2006; Díaz-Moreno et al., 2018). We tentatively propose that periods of weak correlation between SO2 flux and tremor, such as that observed in early 2007, are produced when magmatic overturn slows down, and degassing becomes dominated by permeable gas flow from depth. A further related process, which could play a key role, is the level and condition of magma in the conduit, as a lower magma level together with magma viscosity, gas and crystal content may attenuate the tremor amplitude (e.g., Collier et al., 2006). These processes will be investigated in future work through examination of CO2/SO2 and ground deformation time series. Short-term correlations on the scale of a day between gas flux and tremor during quiescent periods are not observed clearly, perhaps because variations in the continuous degassing process are relatively small compared with the overall degassing rate and fall below measurement noise. Future improvements in the precision and accuracy of SO2 flux quantifications are required to reveal any short-term correlations during passive degassing.
In order to characterize the main oscillations of both intermittent paroxysmal and 2008–2009 eruptive phases, corresponding to the more intense explosive episodes of the study period, an inspection at daylight intraday scale was carried out. Tremor was decomposed to its spectral components between 0.025 and 10.5 Hz signal, and then each volcanic tremor frequency was compared with SO2 flux using correlation analysis. This decomposition allowed statistical identification on the best correlation between SO2 flux and volcanic tremor frequency, and highlights that in different eruptive contexts the geochemical signal correlated with tremor at different dominant frequencies (Figure 4). We found that SO2 flux correlated strongly with volcanic tremor with high frequencies during the most explosive activity. This is consistent with the observation of coupling via infrasound between explosive volcanic activity and the ground, producing a high frequency tremor source (Matoza and Fee, 2014). In the six case studies investigated here, the first three of them, i.e., those falling within the intermittent paroxysmal phase (Figure 2) showed the highest SO2 flux/tremor correlation coefficients at frequencies of 2.5–6.5 Hz (Table 1; Figure 4). Conversely, the case studies pertaining to the 2008–2009 eruptive fissure activity displayed the highest correlations at lower frequencies, i.e., from 0.025 to 2.5 Hz (Table 1; Figure 4). The fact that the two parameters correlated at different tremor frequencies, and that the frequency depended on the intensity of the associated eruptive phenomenon, suggests that on long time scales, some of the signal features could be masked or missed if the RMS is calculated for all frequencies. Falsaperla et al. (2005) analyzed amplitude and frequency content of the seismic signal of the 2001 Mt. Etna’s flank eruption finding considerable changes in the volcanic tremor associated with different styles of eruptive activity. In particular, they observed that the dominant frequency of the signal decreased from ∼5 Hz to 3 Hz during lava fountaining, and further decreased to ∼2 Hz during intense lava emissions, supporting the infrasound-coupling hypothesis as the source of seismic tremor during more explosive activity. There were also periods where SO2 flux and tremor appear anticorrelated on short timescale, resulting from time shifts of the two parameters. Shifts were identified in cases 1, 2, and 4, with gas lagging behind seismic energy release. This behavior, which has also been observed on long-time scales by Leonardi et al., 2000b, has been interpreted as due to increasing pressurization of the volcano’s shallow feeder system, which simultaneously increases tremor but with a lag time before gas release at the onset of eruptive activity (e.g., Young et al., 2003; Nadeau et al., 2011). A further process might rely on the turbulent magma-flow rate in the upper conduit triggered by gas-slug dynamics during ongoing eruptive events (e.g., Parfitt, 2004) coupled with instability of magma column (Bercovici et al., 2013).
These results underline that though the strong association between magmatic degassing and tremor, several process are involved in generating seismic energy. Their mutual behavior might change depending on the physical mechanism of magma flow regime and gas/melt separation in the conduit, revealing the gas-tremor study puzzling. Although further efforts are required to advance our understating in magmatic degassing – seismic tremor release relationship, the strong correlation observed in this study between the two parameters, indicates that degassing generates seismic tremor, and that gas flux might thereby provide a proxy for eruptive style and intensity and short-term warning for impeding eruptions.
Conclusion
Both explosive and effusive eruptions are believed to be largely controlled by volatile content and magma flow rate (e.g., Woods and Cardoso, 1997), and by their mutual modulation within the shallow conduit (e.g., Jaupart and Vergniolle, 1988). Likewise, seismic tremor varies with volcanic activity and is considered originated by fluids dynamic process in volcanic conduit (e.g., Chouet, 1996). Results achieved in this study provide strong evidence that for extended periods the volcanic tremor and SO2 flux signals on Mt. Etna are strongly correlated. The mutual relationship rely on a physical mechanism of magma flow, which through friction between magma and conduit walls produces a tremor signal with frequency dominated between 0.025 and 2.5 Hz. This magma flow provides the source of SO2 emitted persistently at the summit craters during passive degassing and the eruptive vents during effusive activity, through exsolution and transport of SO2 during magma ascent. We tentatively attribute periods of low correlation between tremor and SO2 flux to a different degassing regime, dominated by fluxing of gas from depth, and this hypothesis will be tested in future work. Our examination of explosive activity demonstrates that the best correlation is achieved between high frequency tremor (2.5–6.5 Hz) and SO2 flux, which we attribute to a process of coupling through infrasound between the explosive activity in the atmosphere and the surface, creating high frequency tremor. Our results underline the clear link between the geochemical and the geophysical signals. Their common behavior during both quiescent and eruptive stages, and during transitions between eruptive regimes, emphasizes that magmatic degassing produces volcanic tremor, and that both originate from a common physical mechanism of magma dynamics in the shallow conduit. This study provides a framework for the interpretation of tremor and SO2 degassing to other persistently active basaltic systems worldwide, which will assist in the understanding of the processes and mechanisms controlling unrest and pre-eruptive activity.
Statements
Author contributions
GS and MB coordinated the research and mainly wrote the manuscript with the substantial, direct, and intellectual contribution to the work from all authors.
Acknowledgments
We greatly acknowledge F. Murè and V. Longo for their technical assistance in the FLAME network. We thank Dr. L. Spampinato for useful discussions. We are also grateful to the two reviewers for helpful comments and suggestions.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2018.00157/full#supplementary-material
References
1
AiuppaA.BurtonM.CaltabianoT.RipepeM. (2009). “Anti-correlation between gas flux and volcanic tremor on Stromboli volcano,” inProceedings of the Geophysical Research Abstracts 11, EGU2009-13671 EGU General Assembly 2009, Göttingen.
2
AiuppaA.CannataA.CannavòF.Di GraziaG.FerrariF.GiudiceG.et al (2010). Patterns in the recent 2007–2008 activity of Mount Etna volcano investigated by integrated geophysical and geochemical observations.Geochem. Geophys. Geosyst.11:Q09008. 10.1029/2010GC003168
3
AiuppaA.GiudiceG.GurrieriS.LiuzzoM.BurtonM.CaltabianoT.et al (2008). Total volatile flux from Mount Etna.Geophys. Res. Lett.35:L24302. 10.1029/2008GL035871
4
AkiK.FehlerM.DasS. (1977). Source mechanism of volcanic tremor: fluid-driven crack models and their application to the 1963 Kilauea eruption.J. Volcanol. Geotherm. Res.2259–287. 10.1016/0377-0273(77)90003-8
5
AllardP. (1997). Endogenous magma degassing and storage at Mount Etna.Geophys. Res. Lett.242219–2222. 10.1029/97GL02101
6
AllardP.BurtonM.MurèF. (2005). Spectroscopic evidence for a lava fountain driven by previously accumulated magmatic gas.Nature433407–410. 10.1038/nature03246
7
AloisiM.BonaccorsoA.CannavòF.GambinoS.MattiaM.PuglisiG.et al (2009). A new dyke intrusion style for the Mount Etna May 2008 eruption modelled through continuous tilt and GPS data.Terra Nova21316–321. 10.1111/j.1365-3121.2009.00889.x
8
AlparoneS.AndronicoD.LodatoL.SgroiT. (2003). Relationship between tremor and volcanic activity during the Southeast Crater eruption on Mount Etna in early 2000.J. Geophys. Res.108:2241. 10.1029/2002JB001866
9
AlparoneS.BehnckeB.GiammancoS.NeriM.PriviteraE. (2005). Paroxysmal summit activity at Mt. Etna (Italy) monitored through continuous soil radon measurements.Geophys. Res. Lett.32:L16307. 10.1029/2005GL023352
10
AndronicoD.CristaldiA.ScolloS. (2008). The 4–5 September 2007 lava fountain at South–East Crater of Mt. Etna, Italy.J. Volcanol. Geotherm. Res.173325–328. 10.1016/j.jvolgeores.2008.02.004
11
BaniP.LardyM. (2007). Sulphur dioxide emission rates from Yasur volcano, Vanuatu archipelago.Geophys. Res. Lett.34:L20309. 10.1029/2007GL030411
12
BattagliaJ.AkiK.FerrazziniV. (2005). Location of tremor sources and estimation of lava output using tremor source amplitude on the Piton de la Fournaise volcano: 2. Estimation of lava output.J. Volcanol. Geotherm. Res.147291–308. 10.1016/j.jvolgeores.2005.04.006
13
BeckettF. M.BurtonM.MaderH. M.PhillipsJ. C.PolacciM.RustA. C.et al (2014). Conduit convection driving persistent degassing at basaltic volcanoes.J. Volcanol. Geotherm. Res.28319–35. 10.1016/j.jvolgeores.2014.06.006
14
BehnckeB.FornaciaiA.NeriM.FavalliM.GanciG.MazzariniF. (2016). Lidar surveys reveal eruptive volumes and rates at Etna, 2007–2010.Geophys. Res. Lett.434270–4278. 10.1002/2016GL068495
15
BenoitJ. P.McNuttS. R. (1997). New constraints on source processes of volcanic tremor at Arenal Volcano, Costa Rica, using broadband seismic data.Geophys. Res. Lett.24449–452. 10.1029/97GL00179
16
BercoviciD.JellinekA. M.MichautC.RomanD. C.MorseR. (2013). Volcanic tremors and magma wagging: gas flux interactions and forcing mechanism.Geophys. J. Int.1951001–1022. 10.1093/gji/ggt277
17
BonaccorsoA.BonforteA.CalvariS.Del NegroC.Di GraziaG.GanciG.et al (2011a). The initial phases of the 2008–2009 Mount Etna eruption: a multidisciplinary approach for hazard assessment.J. Geophys. Res.116:B03203. 10.1029/2010JB007906
18
BonaccorsoA.CaltabianoT.CurrentiG.Del NegroC.GambinoS.GanciG.et al (2011b). Dynamics of a lava fountain revealed by geophysical, geochemical and thermal satellite measurements: the case of 10 April 2011 Mt. Etna eruption.Geophys. Res. Lett.38:L24307. 10.1029/2011GL049637
19
BrandsdóttirB.EinarssonP. (1992). “Volcanic tremor and low frequency earthquakes in Iceland,” inVolcanic Seismology, edsGaspariniP.ScarpaR.AkiK. (Canberra: IAVCEI), 212–222.
20
BrunoN.CaltabianoT.GrassoM. F.PortoM.RomanoR. (1995). “Studio dell ’emissione di anidride solforosa (SO2) all’ Etna - 1993/1994,” inProgetto-Etna 1993-1995, edsFerrucciF.InnocentiF. (Pisa: Giardini Editore), 29–33.
21
BryanC. J.SherburnS. (2003). Eruption-induced modifications to volcanic seismicity at Ruapehu, New Zealand, and its implications for eruption forecasting.Bull. Volcanol.6530–42. 10.1007/s00445-002-0236-z
22
BurtonM.AllardP.MurèF.OppenheimerC. (2003). “FTIR remote sensing of fractional magma degassing at Mount Etna, Sicily,” inVolcanic DegassingVol. 213edsOppenheimerC.PyleD. M.BarclayJ. (London: Geological Society Special Publication), 281–293. 10.1144/GSL.SP.2003.213.01.17
23
BurtonM.MaderH. M.PolacciM. (2007). The role of gas percolation in quiescent degassing o of persistently active basaltic volcanoes.Earth Planet. Sci. Lett.26446–60. 10.1016/j.epsl.2007.08.028
24
BurtonM. R.CaltabianoT.MurèF.SalernoG.RandazzoD. (2009). SO2 flux from Stromboli during the 2007 eruption: results from the FLAME network and traverse measurements.J. Volcanol. Geotherm. Res.182214–220. 10.1016/j.jvolgeores.2008.11.025
25
BurtonM. R.PrataF.PlattU. (2015a). Volcanological applications of SO2 cameras.J. Volcanol. Geotherm. Res.3002–6. 10.1016/j.jvolgeores.2014.09.008
26
BurtonM. R.SalernoG. G.D’AuriaL.CaltabianoT.MurèF.MaugeriR. (2015b). SO2 flux monitoring at Stromboli with the new permanent INGV SO2 Camera system: a comparison with the FLAME network and seismological data.J. Volcanol. Geotherm. Res.30095–102. 10.1016/j.jvolgeores.2015.02.006
27
CaltabianoT.BurtonM.GiammancoS.AllardP.BrunoN.MurèF.et al (2004). “Volcanic gas emissions from the summit craters and flanks of Mt. Etna, 1987-2000,” inMt. Etna: Volcano LaboratoryVol. 143edsBonaccorsoA.CalvariS.ColtelliM.Del NegroC.FalsaperlaS. (Washington, DC: American Geophysical Union), 111–128. 10.1029/143GM08
28
CalvariS.InguaggiatoS.PuglisiG.RipepeM.RosiM. (2008). The Stromboli Volcano: An Integrated Study of the 2002-2003 Eruption. Geophysical Monograph Series182. Washington, DC: American Geophysical Union, 10.1029/GM182
29
CalvariS.SalernoG. G.SpampinatoL.GouhierM.La SpinaA.PecoraE.et al (2011). An unloading foam model to constrain Etna’s 11–13 January 2011 lava fountaining episode.J. Geophys. Res.116:B11207. 10.1029/2011JB008407
30
CampionR.Delgado-GranadosH.MoriT. (2015). Image-based correction of the light dilution effect for SO2 camera measurements.J. Volcanol. Geotherm. Res.30048–57. 10.1016/j.jvolgeores.2015.01.004
31
CannataA.CataniaA.AlparoneS.GrestaS. (2008). Volcanic tremor at Mt. Etna: inferences on magma dynamics during effusive and explosive activity.J. Volcanol. Geotherm. Res.17819–31. 10.1016/j.jvolgeores.2007.11.027
32
CannataA.Di GraziaG.AliottaM.CassisiC.MontaltoP.PatanèD. (2013). Monitoring Seismo-volcanic and infrasonic signals at volcanoes: Mt. Etna case study.Pure Appl. Geophys.1701751–1771. 10.1007/s00024-012-0634-x
33
CannataA.Di GraziaG.GiuffridaM.GrestaS.PalanoM.SciottoM.et al (2018). Space- time evolution of magma storage and transfer at Mt. Etna volcano (Italy): the 20152016 reawakening of Voragine crater.Geochem. Geophys. Geosys.19471–495. 10.1002/2017GC007296
34
CannataA.GiudiceG.GuerrieriS.MontaltoP.AlparoneS.Di GraziaG.et al (2009a). Relationship between soil CO2 flux and volcanic tremor at Mt. Etna: implications for magma dynamics.Environ. Earth Sci.61477–489. 10.1007/s12665-009-0359-z
35
CannataA.MontaltoP.PriviteraE.RussoG.GrestaS. (2009b). Tracking eruptive phenomena by infrasound: May 13, 2008 eruption at Mt. Etna.Geophys. Res. Lett.36:L05304. 10.1029/2008GL036738
36
ChouetB. (1996). Long-Period volcano seismicity: its source and use in eruption forecasting.Nature380309–316. 10.1038/380309a0
37
ChouetB.KoyanagiR. Y.AkiK. (1987). “Origin of volcanic tremor in Hawaii, Part II: theory and discussion,” inVolcanism in Hawaii, U.S. Geological Survey Professional PaperVol. 1350edsDeckerR. W.WrightT. L.StaufferP. H. (Washington, DC: US Government Publishing office), 1259–1280.
38
CollierL.NeubergJ.LenskyN.LyakhovskyV. (2006). Attenuation in gas-charged magma.J. Volcanol. Geotherm. Res.15321–36. 10.1016/j.jvolgeores.2005.08.009
39
CorsaroR. A.AndronicoD.BehnckeB.BrancaS.CaltabianoT.CiancittoF.et al (2017). Monitoring the December 2015 summit eruptions of Mt. Etna (Italy): implications on eruptive dynamics.J. Volcanol. Geotherm. Res.34153–69. 10.1016/j.jvolgeores.2017.04.018
40
CorsaroR. A.MiragliaL. (2009). Dynamics of magma in the plumbing system of Mt. Etna volcano, Sicily, Italy: a contribution from petrologic data of volcanic erupted from 2007 to 2009.EOS Trans. AGU90:1690.
41
CurrentiG.NapoliR.Di StefanoA.GrecoF.Del NegroC. (2011). 3D integrated geophysical modeling for the 2008 magma intrusion at Etna: constraints on rheology and dike overpressure.Phys. Earth Planet. Int.18544–52. 10.1016/j.pepi.2011.01.002
42
DavisJ. C. (1986). Statistics and Data analysis in Geology, 2nd Edn.New York, NY: John Wiley & Sons.
43
Di GraziaG.CannataA.MontaltoP.PatanèD.PriviteraE.ZuccarelloL.et al (2009). A multiparameter approach to volcano monitoring based on 4D analyses of seismo-volcanic and acoustic signals: the 2008 Mt. Etna eruption.Geophys. Res. Lett.36:L18307. 10.1029/2009GL039567
44
Di GraziaG.FalsaperlaS.LangerH. (2006). Volcanic tremor location during the 2004 Mount Etna lava effusion.Geophys. Res. Lett.33:L04304. 10.1029/2005GL025177
45
Díaz-MorenoA.BarberiG.CocinaO.KoulakovI.ScarfìL.ZuccarelloL.et al (2018). New insights on Mt. Etna’s crust and relationship with the regional tectonic framework from joint active and passive P-wave seismic tomography.Surv. Geophys.3957–97. 10.1007/s10712-017-9425-3
46
DoukasM. P.GerlachT. M. (1995). “Sulfur dioxide scrubbing during the 1992 eruptions of Crater Peak, Mount Spurr Volcano, Alaska,” inThe 1992 Eruptions of Crater Peak Vent, Mount Spurr Volcano, AlaskaVol. 2139ed.KeithT. E. C. (Reston, VA: U.S. Geological Survey), 47–57.
47
EdmondsM.HerdR. A.GalleB.OppenheimerC. M. (2003a). Automated, high time- resolution measurements of SO2 flux at Soufrière Hills Volcano, Montserrat.Bull. Volcanol.65578–586. 10.1007/s00445-003-0286-x
48
EdmondsM.OppenheimerC.PyleD. M.HerdR. A.ThompsonG. (2003b). SO2 emissions from Soufrière Hills Volcano and their relationship to conduit permeability, hydrothermal interaction and degassing regime.J. Volcanol. Geotherm. Res.12423–43. 10.1016/S0377-0273(03)00041-6
49
EreditatoD.LuongoG. (1994). Volcanic tremor wave field during quiescent and eruptive activity at Mt. Etna (Sicily).J. Volcanol. Geotherm. Res.61239–251. 10.1016/0377-0273(94)90006-X
50
FalsaperlaS.AlparoneS.D’AmicoS.Di GraziaG.FerrariF.LangerH.et al (2005). Volcanic tremor at Mt. Etna, Italy, preceding and accompanying the eruption of July– August, 2001.Pure Appl. Geophys.1621–22. 10.1126/sciadv.1700219
51
FischerT. P.MorrisseyM. M.CalvacheM. L. V.GòmezD. M.TorresR. C.StixJ.et al (1994). Correlations between SO2 flux and long-period seismicity at Galeras volcano.Nature368135–137. 10.1038/368135a0
52
GordeevE. (1993). Modeling of volcanic tremor as explosive point sources in a singled-layered, elastic half-space.J. Geophys. Res.9819687–19703. 10.1029/93JB00348
53
GrestaS.ImposaS.PatanèD.PatanèG. (1987). Volcanic tremor at Mt. Etna: state-of- the-art and perspectives.Pure Appl. Geophys.1251079–1095. 10.1007/BF00879369
54
GrestaS.MontaltoA.PatanèG. (1991). Volcanic tremor at Mount Etna (January 1984–March 1985): its relationship to the eruptive activity and modelling of the summit feeding system.Bull. Volcanol.53309–320. 10.1007/BF00414527
55
HelikerC.KauahikauaJ.SherrodD. R.LisowskiM.CervelliP. (2003). “The rise and fall of Pu‘u ‘Ô‘ô cone, 1983–2002,” inThe Pu‘u ‘Ô‘ô -Kûpaianaha Eruption of Kîlauea Volcano, Hawai‘i: The First 20 YearsVol. 1676edsHelikerC.SwansonD. A.TakahashiT. J. (Denver, CO: U.S. Geological Survey Professional Paper), 29–51.
56
HibertC.MangeneyA.PolacciM.Di MuroA.VergniolleS.FerrazziniV.et al (2015). Toward continuous quantification of lava extrusion rate: results from the multidisciplinary analysis of the 2 January 2010 eruption of Piton de la Fournaise volcano, La Réunion.J. Geophys. Res. Solid Earth1203026–3047. 10.1002/2014JB011769
57
JamesM. R.ApplegarthL. J.PinkertonH. (2011). Lava channel roofing, overflows, breaches and switching: insights from the 2008–2009 eruption of Mt. Etna.Bull. Volcanol.74107–117. 10.1007/s00445-011-0513-9
58
JaupartC.VergniolleS. (1988). Laboratory models of Hawaiian and Strombolian eruptions.Nature33158–60. 10.1038/331058a0
59
JohnsonJ. H.PolandM. P. (2013). Seismic detection of increased degassing before Kîlauea’s 2008 summit explosion.Nat. Commun.4:1668. 10.1038/ncomms2703
60
KazahayaK.ShinoharaH.SaitoG. (1994). Excessive degassing of Izu-Oshima volcano: magma convection in a conduit.Bull. Volcanol.56207–216. 10.1007/BF00279605
61
KazahayaK.ShinoharaH.UtoK.OdaiM.NalkahoriY.MoriH.et al (2004). Gigantic SO2 emission from Miyakejima volcano, Japan, caused by caldera collapse.Geology32425–428. 10.1130/G20399.1
62
KonstantinouK. I.SchlindweinV. (2002). Nature, wavefield properties and source mechanism of volcanic tremor: a review.J. Volcanol. Geotherm. Res.119161–187. 10.1016/S0377-0273(02)00311-6
63
KuboteraA. (1974). “Volcanic tremors at Aso volcano,” inPhysical VolcanologyVol. 6edsCivettaL.GaspariniP.LuongoG.RapollaA. (Amsterdam: Elsevier), 29–47. 10.1016/B978-0-444-41141-9.50008-5
64
La SpinaA.BurtonM.AllardP.AlparoneA.MuréF. (2015). Open-path FTIR spectroscopy of magma degassing processes during eight lava fountains on Mount Etna.Earth Planet. Sci. Lett.413123–134. 10.1016/j.epsl.2014.12.038
65
La SpinaG.PolacciM.BurtonM.de’ Michieli VitturiM. (2017). Numerical investigation of permeability models for low viscosity magmas: application to the 2007 Stromboli effusive eruption.Earth Planet. Sci. Lett.473279–290. 10.1016/j.epsl.2017.06.013
66
LangerH.FalsaperlaS.MessinaA.SpampinatoS.BehnckeB. (2010). Detecting imminent eruptive activity at Mt. Etna, Italy, in 2007–2008 through pattern classification of volcanic tremor data.J. Volcanol. Geotherm. Res.2001–17. 10.1016/j.jvolgeores.2010.11.019
67
LeonardiS.GrestaS.MulargiaF. (2000a). Cross-correlation between tremor and SO2 flux data from Mount Etna volcano, 1987-1992.Phys. Chem. Earth25737–740. 10.1016/S1464-1895(00)00114-9
68
LeonardiS.GrestaS.MulargiaF. (2000b). Searching for a significant correlation between volcanic tremor amplitude and SO2 emissions at Mount Etna volcano, Sicily.Geophys. J. Int.141832–834. 10.1046/j.1365-246X.2000.00114.x
69
LopezT.ThomasH. E.PrataA. J.AmigoA.FeeD.MorianoD. (2015). Volcanic plume characteristics determined using an infrared imaging camera.J. Volcanol. Geotherm. Res.300148–166. 10.1016/j.jvolgeores.2014.12.009
70
MalinconicoL. L. (1979). Fluctuations in SO2 emission during recent eruptions of Etna.Nature27843–45. 10.1038/278043a0
71
MatozaR. S.FeeD. (2014). Infrasonic component of volcano-seismic eruption tremor.Geophys. Res. Lett.411964–1970. 10.1002/2014GL059301
72
McKillupS.DyarM. D. (2010). Geostatistics Explained. An Introductory Guide for Earth Scientists.Cambridge: Cambridge University Press, 396. 10.1017/CBO9780511807558
73
McNuttS. R. (1986). Observations and analysis of B-type earthquakes, explosions, and volcanic tremor at Pavlof Volcano, Alaska.Bull. Seismol. Soc. Am.76153–175.
74
McNuttS. R. (1992). “Volcanic tremor,” inEncyclopedia of Earth System ScienceVol. 4ed.NierenbergW. A. (Cambridge, MA: Academic Press),417–425.
75
McNuttS. R. (2002). “Volcano seismology and monitoring for eruptions,” inInternational Geophysics, International Handbook of Earthquake & Engineering Seismology, Part BVol. 81edsLeeW.KanamoriH.JenningsP.KisslingerC. (New York, NY: Elsevier), 383–406. 10.1016/S0074-6142(02)80228-5
76
MetrichN.AllardP.SpilliaertN.AndronicoD.BurtonM. (2004). 2001 flank eruption of the alkali- and volatile-rich primitive basalt responsible for Mount Etna’s evolution in the last three decades.Earth Planet. Sci. Lett.2281–17. 10.1016/j.epsl.2004.09.036
77
MillerA. D.StewartR. C.WhiteR. A.LuckettR.BaptieB. J.AspinallW. P.et al (1998). Seismicity associated with dome growth and collapse at the Soufriere Hills Volcano, Montserrat.Geophys. Res. Lett.253401–3404. 10.1029/98GL01778
78
MiwaS.MoriM.HibikiT. (2015). Two-phase flow induced vibration in piping systems.Prog. Nucl. Energy78270–284. 10.1016/j.pnucene.2014.10.003
79
MoriT.BurtonM. (2006). The SO2 camera: a simple, fast and cheap method for ground- based imaging of SO2 in volcanic plumes.Geophys. Res. Lett.33:L24804. 10.1029/2006GL027916
80
MoriT.MoriT.KazahayaK.OhwadaM.HirabayashiJ.YoshikawaS. (2006). Effect of UV scattering on SO2 emission rate measurements.Geophys. Res. Lett.33:L17315. 10.1029/2006GL026285
81
NadeauP. A.PalmaJ. L.WaiteG. P. (2011). Linking volcanic tremor, degassing, and eruption dynamics via SO2 imaging.Geophys. Res. Lett.38:L01304. 10.1029/2010GL045820
82
NadeauP. A.WernerC. A.WaiteG. P.CarnS. A.BrewerI. D.EliasT.et al (2015). Using SO2 camera imagery and seismicity to examine degassing and gas accumulation at Kîlauea Volcano, May 2010.J. Volcanol. Geotherm. Res.30070–80. 10.1016/j.jvolgeores.2014.12.005
83
NeubergJ.PointerT. (2000). Effects of volcano-topography on seismic broadband waveforms.Geophys. J. Int.143239–248. 10.1046/j.1365-246x.2000.00251.x
84
OmoriF. (1911). The Usu-san eruption and earthquake and elevation phenomena.Bull. Imp. Earthq. Investig. Comm.51–38.
85
PalmaJ. L.CalderE. S.BasualtoD.BlakeS.RotheryD. A. (2008). Correlations between SO2 flux, seismicity, and outgassing activity at the open vent of Villarrica volcano, Chile.J. Geophys. Res.113:B10201. 10.1029/2008JB005577
86
ParfittE. A. (2004). A discussion of the mechanisms of explosive basaltic eruptions.J. Volcanol. Geotherm. Res.13477–107. 10.1016/j.jvolgeores.2004.01.002
87
PatanèD.AiuppaA.AloisiM.BehnckeB.CannataA.ColtelliM.et al (2013). Insights into magma and fluid transfer at Mount Etna by a multiparametric approach: a model of the events leading to the 2011 eruptive cycle.J. Geophys. Res. Solid Earth1183519–3539. 10.1002/jgrb.50248
88
PatanèD.BarberiG.CocinaO.De GoriP.ChiarabbaC. (2006). Time-resolved seismic tomography detects magma intrusions at Mount Etna.Science313821–823. 10.1126/science.1127724
89
PatanèD.Di GraziaG.CannataA.MontaltoP.BoschiE. (2008). Shallow magma pathway geometry at Mt. Etna volcano.Geochem. Geophys. Geosys.9:Q12021.
90
PlattA.StutzU. (2008). Differential Optical Absorption Spectroscopy - Principles and Applications.Heidelberg: Springer, 10.1007/978-3-540-75776-4
91
PolandM. P.MikliusA.SuttonJ. A.ThornberC. R. (2012). A mantle-driven surge in magma supply to Kîlauea volcano during 2003–2007.Nat. Geosci.5295–300. 10.1038/NGEO1426
92
SalernoG. G.BurtonM. R.OppenheimerC.CaltabianoT.RandazzoD.BrunoN.et al (2009a). Three-years of SO2 flux measurements of Mt. Etna using an automated UV scanner array: comparison with conventional traverses and uncertainties in flux retrieval.J. Volcanol. Geotherm. Res.18376–83. 10.1016/j.jvolgeores.2009.02.013
93
SalernoG. G.BurtonM. R.OppenheimerC.CaltabianoT.TsanevV.BrunoN. (2009b). Novel retrieval of volcanic SO2 abundance from ultraviolet spectra.J. Volcanol. Geotherm. Res.181141–153. 10.1016/j.jvolgeores.2009.01.009
94
SchickR.MugionoR. (1991). Volcanic Tremor and Magma Flow.Jülich: Forschungszentrum Jülich GmbH, 200.
95
ShinoharaH. (2005). New technique to estimate volcanic gas composition: plume measurements with a portable multi-sensor system.J. Volcanol. Geotherm. Res.143319–333. 10.1016/j.jvolgeores.2004.12.004
96
ShinoharaH. (2008). Excess degassing from volcanoes and its role on eruptive and intrusive activity.Rev. Geophys.46:RG4005. 10.1029/2007RG000244
97
SpilliaertN.AllardP.MètrichN.SobolevA. V. (2006). Melt inclusion record of the conditions of ascent, degassing, and extrusion of volatile-rich alkali basalt during the powerful 2002 flank eruption of Mount Etna (Italy).J. Geophys. Res.111:B04203. 10.1029/2005JB003934
98
SteffkeA. M.HarrisA. J. L.BurtonM.CaltabianoT.SalernoG. G. (2011). Coupled use of COSPEC and satellite measurements to define the volumetric balance during effusive eruptions at Mt. Etna, Italy.J. Volcanol. Geotherm. Res.20547–53. 10.1016/j.jvolgeores.2010.06.004
99
SteinbergG. S.SteinbergA. S. (1975). On possible causes of volcanic tremor.J. Geophys. Res.801600–1604. 10.1029/JB080i011p01600
100
SuttonA. J.EliasT.GerlachT. M.StokesJ. B. (2001). Implications for eruptive processes as indicated by sulfur dioxide emission from Kîlauea Volcano, Hawai’i, 1979-1997.J. Volcanol. Geotherm. Res.108283–302. 10.1016/S0377-0273(00)00291-2
101
SwanA. R. H.SandilandsM. (1995). Introduction to Geological Data analysis.Hoboken, NJ: Blackwell Science.
102
TaquetN.Meza HernándezI.StremmeW.BezanillaA.GrutterM.CampionR.et al (2017). Continuous measurements of SiF4 and SO2 by thermal emission spectroscopy: insight from a 6-month survey at the Popocatépetl volcano.J. Volcanol. Geotherm. Res.341255–268. 10.1016/j.jvolgeores.2017.05.009
103
ThompsonG.McNuttS. R.TytgatG. (2002). Three distinct regimes of volcanic tremor associated with the eruption of Shishaldin Volcano, Alaska 1999.Bull. Volcanol.64535–547. 10.1007/s00445-002-0228-z
104
Vargas-BracamontesD. M.NavaA. F.Reyes-DávilaG. A. (2009). Time-scale wavelet patterns related to the 1998–1999 eruptions of the Colima volcano, and their possible implications for eruption forecasting.J. Volcanol. Geotherm. Res.184271–284. 10.1016/j.jvolgeores.2009.01.025
105
VoightB.SparksR. S. J.MillerA. D.StewartR. C.HoblittR. P.ClarkeA.et al (1999). Magma flow instability and cyclic activity at Soufriere Hills volcano, Montserrat, British West Indies.Science2831138–1142. 10.1126/science.283.5405.1138
106
WallaceP. J.EdmondsM. (2011). “The sulfur budget in magmas: evidence from melt inclusions, submarine glasses, and volcanic gas emissions,” inReviews in Mineralogy and GeochemistryVol. 73edsBehrensH.WebsterJ. D. (Chantilly, VA: Mineralogical Society of America), 215–246.
107
WallaceP. J.GerlachT. M. (1994). Magmatic vapor source for sulfur dioxide released during volcanic eruptions: evidence from Mount Pinatubo.Science265497–499. 10.1126/science.265.5171.497
108
WatsonI. M.OppenheimerC.VoightB.FrancisP. W.ClarkeA.StixJ.et al (2000). The relationship between degassing and ground deformation at Soufriere Hills Volcano, Montserrat.J. Volcanol. Geotherm. Res.98117–126. 10.1016/S0377-0273(99)00187-0
109
Williams-JonesG.StixJ.HeiligmannM.BarqueroJ.FernandezE.GonzalezE. D. (2001). A model of degassing and seismicity at Arenal volcano, Costa Rica.J. Volcanol. Geotherm. Res.108121–139. 10.1016/S0377-0273(00)00281-X
110
Williams-JonesG.StixJ.HicksonC. (2008). The COSPEC Cookbook: Making SO2 Measurements at Active Volcanoes.Rome: IAVCEI, 1. 10.13140/RG.2.2.13728.99845
111
WoodsA. W.CardosoS. S. S. (1997). Triggering basaltic volcanic eruptions by bubble-melt separation.Nature385518–520. 10.1038/385518a0
112
YoungS. R.FrancisP. W.BarclayJ.CasadevallT. J.GardnerC. A.DarrouxB.et al (1998). Monitoring SO2 emissions at the Soufriere Hills Volcano: implications for changes in eruptive conditions.Geophys. Res. Lett.253681–3684. 10.1029/98GL01406
113
YoungS. R.VoightB.DuffellH. J. (2003). Magma extrusion dynamics revealed by high- frequency gas monitoring at Soufrière Hills volcano, Montserrat.Geol. Soc. Spec. Publ.213219–230. 10.1144/GSL.SP.2003.213.01.13
114
ZarJ. H. (1972). Significance testing of the spearman rank correlation coefficient.J. Am. Stat. Assoc.67578–580. 10.1080/01621459.1972.10481251
115
ZobinV. M. (2003). Introduction to Volcanic Seismology (Developments in Volcanology; 6).Amsterdam: Elsevier Science.
116
ZuccarelloL.BurtonM. R.SaccorottiG.BeanC. J.PatanèD. (2013). The coupling between very long period seismic events, volcanic tremor, and degassing rates at Mount Etna volcano.J. Geophys. Res. Solid Earth1184910–4921. 10.1002/jgrb.50363
Summary
Keywords
Mt. Etna, SO2 flux, volcanic tremor, eruptive and quiescent degassing, volcano monitoring
Citation
Salerno GG, Burton M, Di Grazia G, Caltabiano T and Oppenheimer C (2018) Coupling Between Magmatic Degassing and Volcanic Tremor in Basaltic Volcanism. Front. Earth Sci. 6:157. doi: 10.3389/feart.2018.00157
Received
29 June 2018
Accepted
24 September 2018
Published
15 October 2018
Volume
6 - 2018
Edited by
Fátima Viveiros, Instituto de Investigação em Vulcanologia e Avaliação de Riscos (IVAR), Portugal
Reviewed by
Taryn Lopez, University of Alaska Fairbanks, United States; Oleg E. Melnik, Lomonosov Moscow State University, Russia
Updates
Copyright
© 2018 Salerno, Burton, Di Grazia, Caltabiano and Oppenheimer.
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.
*Correspondence: Giuseppe G. Salerno, giuseppe.salerno@ingv.it
This article was submitted to Volcanology, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.