High precision tilt observation at Mt. Etna Volcano,Italy

In 2007–2008, we installed on Mt. Etna two deep tilt stations using high resolution, self-leveling instruments. These installations are a result of accurate instrument tests, site selection, drilling and sensor positioning that has allowed detecting variations related to the principal diurnal and semidiurnal tides for first time on Mt. Etna using tilt data.     

Motif Discovery on Seismic Amplitude Time Series: The Case Study of Mt Etna 2011 Eruptive Activity

Algorithms searching for similar patterns are widely used in seismology both when the waveforms of the events of interest are known and when there is no a priori-knowledge. Such methods usually make use of the cross-correlation coefficient as a measure of similarity; if there is no a-priori knowledge, they behave as brute-force searching algorithms. The disadvantage of these methods, preventing or limiting their application to very large datasets, is computational complexity. The Mueen–Keogh (MK) algorithm overcomes this limitation by means of two optimization techniques—the early abandoning concept and space indexing. Here, we apply the MK algorithm to amplitude time series retrieved from seismic signals recorded during episodic eruptive activity of Mt Etna in 2011. By adequately tuning the input to the MK algorithm we found eight motif groups characterized by distinct seismic amplitude trends, each related to a different phenomenon. In particular, we observed that earthquakes are accompanied by sharp increases and decreases in seismic amplitude whereas lava fountains are accompanied by slower changes. These results demonstrate that the MK algorithm, because of its particular features, may have wide applicability in seismology.     

New formulae for estimating lava flow volumes at Mt. Etna Volcano, Sicily

 A new data set of Etna lava flows erupted since 1868 has been compiled from eight topographic maps of the volcano published at intervals since then. Volumes of 59 flows or groups of flows were measured from topographic difference maps. Most of these volumes are likely to be considerably more accurate than those published previously. We cut the number of flow volumes down to 25 by selecting those examples for which the volume of an individual eruption could be derived with the highest accuracy. This refined data set was searched for high correlations between flow volume and more directly measurable parameters. Only two parameters showed a correlation coefficient of 70% or greater: planimetric flow area A (70%) and duration of the eruption D (79%). If only short duration (<18 days) flows were used, flow length cubed, L³, had a correlation coefficient of 98%. Using combinations of measured parameters, much more significant correlations with volume were found. Dh had a correlation coefficient of 90% (h is the hydrostatic head of magma above the vent), and  , 92% (where W is mean width and E is the degree of topographic enclosure), and a combination of the two , 97%. These latter formulae were used to derive volumes of all eruptions back to 1868 to compare with those from the complete data set. Values determined from the formulae were, on average, lower by 16% (Dh), 7% (, and 19% . Received: 30 November 1998 / Accepted: 20 June 1999     

Monitoring Seismo-volcanic and Infrasonic Signals at Volcanoes: Mt. Etna Case Study

Volcanoes generate a broad range of seismo-volcanic and infrasonic signals, whose features and variations are often closely related to volcanic activity. The study of these signals is hence very useful in the monitoring and investigation of volcano dynamics. The analysis of seismo-volcanic and infrasonic signals requires specifically developed techniques due to their unique characteristics, which are generally quite distinct compared with tectonic and volcano-tectonic earthquakes. In this work, we describe analysis methods used to detect and locate seismo-volcanic and infrasonic signals at Mt. Etna. Volcanic tremor sources are located using a method based on spatial seismic amplitude distribution, assuming propagation in a homogeneous medium. The tremor source is found by calculating the goodness of the linear regression fit (R ²) of the log-linearized equation of the seismic amplitude decay with distance. The location method for long-period events is based on the joint computation of semblance and R ² values, and the location method of very long-period events is based on the application of radial semblance. Infrasonic events and tremor are located by semblance–brightness- and semblance-based methods, respectively. The techniques described here can also be applied to other volcanoes and do not require particular network geometries (such as arrays) but rather simple sparse networks. Using the source locations of all the considered signals, we were able to reconstruct the shallow plumbing system (above sea level) during 2011.     

Time Correlation Analysis of the Microseismicity of the Low Eastern Flank of Mt. Etna Volcano (Italy)

—An analysis in terms of time correlation functions has been applied to the time distribution of microseismicity. We considered the single station detections (M < 2) recorded by each one of the four working stations of the network of the Osservatorio Sismologico di Protezione Civile, Acireale (Catania), located on the active faults systems (Timpe) of the low eastern flank of Mt. Etna. Information obtained on the time regularity of seismic crisis and on the time correlation between the activities recorded by different stations allowed a better understanding of the role of the two main structures. The Timpa di Santa Tecla is the most seismogenic structure in the area and probably played a role in the magma transfer process for the 1989 eruption. Modifications of the local stress field by actions of this feature may have produced gravitational instability along a second structural system (the Timpa di S. Leonardello) with a time delay of about 100 days. By contrast, magma penetration associated with the larger eruption of 1991–1993 had a deeper rise system and had no effect on local stress conditions, consequently there was no correlation between volcanism and seismicity.     

Volcanic tremor at Mt. Etna: State-of-the-art and perspectives

Volcanic tremor on Etna seems to have its origin within the main magma feeding system. On the basis of both spectral analyses at two permanent seismic stations and periodical measurements along the slopes of the volcano, two distinct sources are proposed. The former, characterized by low frequency contents (f<1.5 Hz), is located in a 2 km deep flat magma chamber, whereas the latter source seems to be linked to the upper part of the active vents.Turbulent motions in the magma-gas mixture, induced by escaping gases within the conduits, is one proposed cause of volcanic tremor on Etna (Seidl et al., 1981).From spectral analyses we propose approximate models of the feeding system of the main summit craters.Time variations of tremor energy were also investigated, and no regular patterns have been observed for the studied eruptions.More systematic information seems to be needed for a better knowledge of both the source model and location, and correlation between tremor features and volcanic activity.     

Intrusive Mechanisms at Mt. Etna Forerunning the July-August 2001 Eruption from Seismic and Ground Deformation Data

In this work we present seismological and ground deformation evidence for the phase preparing the July 18 to August 9, 2001 flank eruption at Etna. The analysis performed, through data from the permanent seismic and ground deformation networks, highlighted a strong relationship between seismic strain release at depth and surface deformation. This joint analysis provided strong constraints on the magma rising mechanisms. We show that in the last ten years, after the 1991–1993 eruption, an overall accumulation of tension has affected the volcano. Then we investigate the months preceding the 2001 eruption. In particular, we analyse the strong seismic swarm on April 20–24, 2001, comprising more than 200 events (M_max = 3.6) with prevalent dextral shear fault mechanisms in the western flank. The swarm showed a ca. NE-SW earthquake alignment which, in agreement with previous cases, can be interpreted as the response of the medium to an intrusive process along the approximately NNW-SSE volcano-genetic trend. These mechanisms, leading to the July 18 to August 9, 2001 flank eruption, are analogous to ones observed some months before the 1991–1993 flank eruption and, more recently, in January 1998 before the February-November 1999 summit eruption.     

Towards an Automatic Monitoring System of Infrasonic Events at Mt. Etna: Strategies for Source Location and Modeling

Active volcanoes characterized by open conduit conditions generate sonic and infrasonic signals, whose investigation provides useful information for both monitoring purposes and studying the dynamics of explosive processes. In this work, we discuss the automatic procedures implemented for a real-time application to the data acquired by a permanent network of five infrasound stations running at Mt. Etna volcano. The infrasound signals at Mt. Etna consist in amplitude transients, called infrasound events. The adopted procedure uses a multi-algorithm approach for event detection, counting, characterization and location. It is designed for an efficient and accurate processing of infrasound records provided by single-site and array stations. Moreover, the source mechanism of these events can be investigated off-line or in near real-time by using three different models: (1) Strombolian bubble; (2) resonating conduit and (3) Helmholtz resonator. The infrasound waveforms allow us to choose the most suitable model, to get quantitative information about the source and to follow the time evolution of the source parameters.     

Volcanic tremor at Mt. Etna: Inferences on magma dynamics during effusive and explosive activity

During 1999, the volcanic activity at Mt. Etna was both explosive and effusive at the summit craters: Strombolian activity, lava fountains and lava flows affected different areas of the volcano, involving three of the four summit craters. Results from analysis of the 1999 volcanic tremor features are shown at two different time scales. First, the long-term time variation of the features of the volcanic tremor (including spectral and polarization parameters), during the entire year, was compared with the evolution of the eruptive activity. This approach demonstrated the good agreement between tremor data and observed eruptive activity; the activation of different tremor sources was suggested. Then, a more refined analysis of the volcanic tremor, recorded during 14 lava fountain eruptions, was performed. In particular, a shift of the dominant frequencies towards lower values was noted which corresponds with increasing explosive activity. Similar behaviour in the frequency content has already been observed in other explosive eruptions at Mt. Etna as well as on other volcanoes. This behaviour has been explained in terms of either an increase in the tremor source dimension or a decrease in the sound speed in the magma within the conduit. These results confirm that the volcanic tremor is a powerful tool for better understanding the physical processes controlling explosive eruptions at Mt. Etna volcano.     

Review of Microgravity Observations at Mt. Etna: A Powerful Tool to Monitor and Study Active Volcanoes

Microgravity observations at Mt. Etna have been routinely performed as both discrete (since 1986) and continuous (since 1998) measurements. In addition to describing the methodology for acquiring and reducing gravity data from Mt. Etna, this paper provides a collection of case studies aimed at demonstrating the potential of microgravity to investigate the plumbing system of an active volcano and detect forerunners to paroxysmal volcanic events. For discrete gravity measurements, results from 1994–1996 and 2001 are reported. During the first period, the observed gravity changes are interpreted within the framework of the Strombolian activity which occurred from the summit craters. Gravity changes observed during the first nine months of 2001 are directly related to subsurface mass redistributions which preceded, accompanied and followed the July-August 2001 flank eruption of Mt. Etna. Two continuous gravity records are discussed: a 16-month (October 1998 to February 2000) sequence and a 48-hour (26–28 October, 2002) sequence, both from a station within a few kilometers of the volcano's summit. The 16-month record may be the longest continuous gravity sequence ever acquired at a station very close to the summit zone of an active volcano. By cross analyzing it with contemporaneous discrete observations along a summit profile of stations, both the geometry of a buried source and its time evolution can be investigated. The shorter continuous sequence encompasses the onset of an eruption from a location only 1.5 km from the gravity station. This gravity record is useful for establishing constraints on the characteristics of the intrusive mechanism leading to the eruption. In particular, the observed gravity anomaly indicates that the magma intrusion occurred “passively” within a fracture system opened by external forces.     

Geochemistry of gases and waters discharged by the mud volcanoes at Paternò, Mt. Etna (Italy)

 Approximately 20 km south of Mt. Etna craters, at the contact between volcanic and sedimentary formations, three mud volcanoes discharge CO₂-rich gases and Na–Cl brines. The compositions of gas and liquid phases indicate that they are fed by a hydrothermal system for which temperatures of 100–150  °C were estimated by means of both gas and solute geothermometry. The hydrothermal system may be associated with CO₂-rich groundwaters over a large area extending from the central part of Etna to the mud volcanoes. Numerous data on the He, CH₄, CO₂ composition of the gases of the three manifestations, sampled over the past 5 years, indicate clearly that variations are due to separation processes of a CO₂-rich gas phase from the liquid. The effects of these processes have to be taken into account in the interpretation of the monitoring data collected for the geochemical surveillance of Etna volcano. Received: 4 September 1995 / Accepted: 14 February 1996     

Continuous Hidden Markov Models: Application to automatic earthquake detection and classification at Las Canãdas caldera,Tenerife

A possible interaction of (volcano-) tectonic earthquakes with the continuous seismic noise recorded in the volcanic island of Tenerife was recently suggested. Also recently the zone close to Las Canadas caldera shows unusual high number of near (< 25 km), possibly volcano-tectonic, earthquakes indicating signs of reawakening of the volcano putting high pressure on the risk analyst. Certainly for both tasks consistent earthquake catalogues provide valuable information and thus there is a strong demand for automatic detection and classification methodologies generating such catalogues. Therefore we adopt methodologies of speech recognition where statistical models, called Hidden Markov Models (HMMs), are widely used for spotting words in continuous audio data. In this study HMMs are used to detect and classify volcano-tectonic and/or tectonic earthquakes in continuous seismic data. Further the HMM detection and classification is evaluated and discussed for a one month period of continuous seismic data at a single seismic station. Being a stochastic process, HMMs provide the possibility to add a confidence measure to each classification made, basically evaluating how “sure” the algorithm is when classifying a certain earthquake. Moreover, this provides helpful information for the seismological analyst when cataloguing earthquakes. Combined with the confidence measure the HMM detection and classification can provide precise enough earthquake statistics, both for further evidence on the interaction between seismic noise and (volcano-) tectonic earthquakes as well as for incorporation in an automatic early warning system.     

Seismic hazard assessment in the volcanic region of Mt. Etna (Italy): a probabilistic approach based on macroseismic data applied to volcano-tectonic seismicity

In the framework of the UPStrat-MAFA project, a seismic hazard assessment has been undertaken in the volcanic region of Mt. Etna as a first step in studies aimed at evaluating the risk on an urban scale. The analysis has been carried out with the SASHA code which uses macroseismic data in order to calculate, starting from the site seismic history, the maximum intensity value expected in a given site with a probability of exceedance of 10 % (I_ref), for a fixed exposure time. Depending on the aims of the project, hazard is estimated for local volcano-tectonic seismicity and short exposure times (10 and 30 years), without taking into account the contribution of “regional” events characterized by much longer recurrence times. Results from tasks A, B and D of the project have produced an updated macroseismic dataset, better performing attenuation models and new tools for SASHA, respectively. The maps obtained indicate that the eastern flank of Etna, the most urbanized sector of the volcano, is characterized by a high level of hazard with I_ref values up to degree VIII EMS, and even IX EMS locally. The disaggregated data analysis allows recognizing the “design earthquake” and the seismogenic fault which most contribute to the hazard at a site-scale. The latter analysis is the starting point to select the scenario earthquake to be used in the analyses of tasks C and F of the project dealing with, respectively, synthetic ground motion simulations and the evaluation of the Disruption Index.     

First study of fault trench stratigraphy at Mt. Etna volcano, Southern Italy: understanding Holocene surface faulting along the Moscarello fault

Paleoseismology, the study of past earthquakes based on their geological record in the stratigraphy and landscape, is a successful newly developing field of research. The application of fault trench studies in volcanic environments is one of the youngest branches of paleoseismology. In this paper, we present the results of the first exploratory trenches excavated at Mt. Etna in Sicily, the largest European volcano. Modern surface faulting at Etna is a very well known feature, which poses significant hazard to the local community, both in terms of ground displacement of essential lifelines and ground shacking from frequent damaging earthquakes. However, while the geomorphology and the seismicity of the active fault in the Etna region consistently show very high rates of tectonic activity, the Holocene cumulative throw and slip-rates, along with the nature (coseismic vs. creeping fault slip), dimension and timing of the displacement events, are still poorly constrained. For this purpose, we selected as a sample area the Moscarello fault, one of the most outstanding segment of the Timpe system of active normal faults in the volcano’s lower eastern flank. Displaced landforms and volcanic units at the Fondo Macchia basin, in the central sector of this fault, indicate some hundreds of meters of vertical offset in the last ca. 80 kyr, with a long-term slip-rate substantially higher than 1.5–2.0 mm/yr. According to the historical sources and instrumental observations, the Moscarello fault ruptured four times in the last 150 years during shallow (H < 5 km) and moderate magnitude (M < 4.8) earthquakes. These events were associated with severe damage in a narrow epicentral area (macroseismic intensities up to the IX–X grade of the MSK scale) and extensive surface faulting (end-to-end rupture length up to 6 km, vertical offsets up to 90 cm). This clearly indicates very high modern rates of deformation along this fault. We conducted trench investigations at the Fondo Macchia site, in a point where eyewitnesses observed ca. 20 cm of coseismic vertical displacement after the April 21, 1971, M_s=3.7, earthquake. The excavated sections provided direct stratigraphic evidence for a vertical slip-rate of 1.4–2.7 mm/yr in the last ca. 6 kyr. This should be regarded as a minimum slip-rate for the central section of the fault. We explored a single scarp at a single site, while we know from recent historical observations that several parallel scarps may rupture coseismically at Fondo Macchia. Thus, the relevant deformation rate documented for the modern period might be likely extended back in the past to a time-span of some thousands of years at least. As expected, for such a volcanic environment, the activity rates of the Moscarello fault are also significantly higher than for the Apennines normal faults, typically showing slip-rates lower than 1 mm/yr. The agriculturally reworked trench hangingwall stratigraphy did not allow to recognize individual displacement events. Nevertheless, the sedimentary structures observed in the trench footwall strongly suggest that, as for the last 150–200 years of detailed historical record, fault behavior at Fondo Macchia is governed by coseismic surface displacement rather than fault creep. This research confirms that paleoseismology techniques can be effectively applied also in active volcanic environments, typically characterized by rheology and, consequently, seismicity and fault dynamics very different from those of other tectonic environments in which paleoseismology has been firstly developed and is today extensively applied.     

Unusual magma storage conditions at Mt. Etna (Southern Italy) as evidenced by plagioclase megacryst-bearing lavas: implications for the plumbing system geometry and summit caldera collapse

At Mt. Etna volcano, the emission of plagioclase megacryst-bearing lavas, known locally as “cicirara”, has occurred rarely and generally in association with unusual volcanological phenomena. In this work, we interpret the magma chamber processes and the structural features of the plumbing system that led to the production of these peculiar volcanic rocks, based on a detailed study of plagioclase megacrysts, including their oscillatory zoning, sieve textures, and fluid inclusions. Patchy zoning suggests limited ascent in the deep levels of the plumbing system, based on the plagioclase nucleation threshold and the volatile saturation depth. At intermediate, water-undersaturated levels of the plumbing system ascent is faster, as indicated by crystals with coarse sieve textures. Storage at shallow, water-saturated levels (less than 6 km deep) is associated with oscillatory zoning with very small changes in An. Slightly larger An variations coupled with different wavelengths provide evidence of convection of crystals across distinct zones of the chamber. Stripes of melt inclusions formed at steps of magma ascent and volatile loss, whereas layers of fluid inclusions may be related to episodes of volatile flushing into the magma chamber. In contrast, strongly sieve-textured envelopes with An increase and constant FeO may be related to mixing with more volatile-rich magmas of similar composition. We interpret the repeated occurrence of “cicirara” lavas as evidence that the shallow portion of the plumbing system underwent a progressive coalescence of a complex network of dykes and sills in response to increasing rates of magma supply from depth. Major magma withdrawals from this larger reservoir may be linked to episodes of summit instability associated with major caldera collapses.