Ground deformation modeling of flank dynamics prior to the 2002 eruption of Mt. Etna

On 22 September 2002, 1 month before the beginning of the flank eruption on the NE Rift, an M-3.7 earthquake struck the northeastern part of Mt. Etna, on the westernmost part of the Pernicana fault. In order to investigate the ground deformation pattern associated with this event, a multi-disciplinary approach is presented here. Just after the earthquake, specific GPS surveys were carried out on two small sub-networks, aimed at monitoring the eastern part of the Pernicana fault, and some baselines belonging to the northeastern EDM monitoring network of Mt. Etna were measured. The leveling route on the northeastern flank of the volcano was also surveyed. Furthermore, an investigation using SAR interferometry was performed and also the continuous tilt data recorded at a high precision sensor close to the epicenter were analyzed to constrain the coseismic deformation. The results of the geodetic surveys show a ground deformation pattern that affects the entire northeastern flank of the volcano, clearly shaped by the Pernicana fault, but too strong and wide to be related only to an M-3.7 earthquake. Leveling and DInSAR data highlight a local strong subsidence, up to 7 cm, close to the Pernicana fault. Significant displacements, up to 2 cm, were also detected on the upper part of the NE Rift and in the summit craters area, while the displacements decrease at lower altitude, suggesting that the dislocation did not continue further eastward. Three-dimensional GPS data inversions have been attempted in order to model the ground deformation source and its relationship with the volcano plumbing system. The model has also been constrained by vertical displacements measured by the leveling survey and by the deformation map obtained by SAR interferometry.     

Building vulnerability and seismic risk analysis in the urban area of Mt. Etna volcano (Italy)

The tectonic system of the eastern flank of Mt. Etna volcano (Sicily, Italy) is the source of most of the strongest earthquakes occurring in the area over the last 205 years. A total of 12 events with epicentre intensities ≥VIII EMS have occurred at Mt. Etna, 10 of which were located on the eastern flank. This indicates a mean recurrence time of about 20 years. This area is highly urbanised, with many villages around the volcano at altitudes up to 700 m a.s.l. The southern and eastern flanks are particularly highly populated areas, with numerous villages very close to each other. The probabilistic seismic hazard due to local faults for Mt. Etna was calculated by adopting a site approach to seismic hazard assessment. Only the site histories of local volcano-tectonic earthquakes were considered, leaving out the effects due to strong regional earthquakes that occurred in north-eastern and south-eastern Sicily. The inventory used in this application refers to residential buildings. These data were extracted from the 1991 census of the Italian National Institute of Statistics, and are grouped according to the census sections. The seismic vulnerability of the elements at risk belonging to a given building typology is described by a vulnerability index, in accordance with a damage model based on macroseismic intensities. For the estimation of economic losses due to physical damage to buildings, an integrated impact indicator was used, which is equivalent to the lost building volume. The expected annualised economic earthquake losses were evaluated both in absolute and in relative terms, and were compared with the geographical distribution of seismic hazard and with similar evaluations of losses for other regions.     

Variation ofb values before the Etnean eruption of March 1981

A fissural eruption occurred from the northern flank of Mt. Etna on March 17, 1981, and the associated earthquake activity, recorded by the seismic network operating on the volcano, was carefully examined to detect possible variations ofb values, according to the frequency-magnitude relationship, before the start of the eruption. The analysis of 4000B-type earthquakes was carried out according to the maximum likelihood method. A significant increase in theb value was noticed since the beginning of February 1981, followed by a sharp decrease in the days before the eruption started. The observed variation pattern was related to possible changes on the stress field acting on the volcano.     

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.     

Elastic model for the gravity and elevation changes before the 2001 eruption of Etna volcano

For 5 months before the 2001 Mt. Etna eruption, a progressive gravity decrease was measured along a profile of stations on the southern slope of the volcano. Between January and July 2001, the amplitude of the change reached 80 μGal, while the wavelength of the anomaly was of the order of 15 km. Elevation changes observed through GPS measurements during a period encompassing the 5-month gravity decrease, remained within 4–6 cm over the entire volcano and within 2–4 cm in the zone covered by the microgravity profile. We review both gravity and elevation changes by a model assuming the formation of new cracks, uniformly distributed in a rectangular prism. The inversion problem was formulated following a global optimization approach based on the use of Genetic Algorithms. Although it is possible to explain the observed gravity changes by means of the proposed analytical formulation, the results show that calculated elevation changes are significantly higher than those observed. Two alternative hypotheses are proposed to account for this apparent discrepancy: (1) that the assumptions behind the analytical formulation, used to invert the data, are fallacious at Etna, and thus, numerical models should be utilized; (2) that a second process, enabling a considerable mass decrease to occur without deformation, acted together with the formation of new cracks in the source volume.     

Ground and marine magnetic surveys of the lower eastern flank of Etna volcano (Italy)

In 1996 and 1997, two high-resolution magnetic surveys, one on land and the other at sea, were carried out on the lower eastern flank of Mount Etna. The magnetic surveys, covering an area of about 400 km², aimed to elucidate the relationships between the main tectonic and morphologic features of this flank of Mount Etna. Major features include widespread NNW- and NNE-trending active faults and the Valle del Bove, a depression considered to be the source area of the Chiancone deposit, the largest Etnean volcaniclastic sequence. Magnetic surveys show anomalies that roughly follow the trend of active main structures. Although few magnetization measurements are available for the most representative outcrops of the lower eastern side of Mount Etna, interpretation of the anomalies defines the underground geometry of the Chiancone deposit and its relationship with volcano stratigraphic units and the underlying sedimentary rocks. In particular, a volume of about 14 km³ was ascribed to the Chiancone deposit. Such a large amount of material was likely produced by a catastrophic event, and deposited at different periods at the exit of the Valle del Bove in an area produced by the interaction, on a regional scale, of the main tectonic structures affecting this flank of the 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.     

Hazards from pyroclastic density currents at Mt. Etna (Italy)

Despite the recent recognition of Mount Etna as a periodically violently explosive volcano, the hazards from various types of pyroclastic density currents (PDCs) have until now received virtually no attention at this volcano. Large-scale pyroclastic flows last occurred during the caldera-forming Ellittico eruptions, 15–16 ka ago, and the risk of them occurring in the near future is negligible. However, minor PDCs can affect much of the summit area and portions of the upper flanks of the volcano. During the past ~ 20 years, small pyroclastic flows or base-surge-like vapor and ash clouds have occurred in at least 8 cases during summit eruptions of Etna. Four different mechanisms of PDC generation have been identified during these events: (1) collapse of pyroclastic fountains (as in 2000 and possibly in 1986); (2) phreatomagmatic explosions resulting from mixing of lava with wet rock (2006); (3) phreatomagmatic explosions resulting from mixing of lava with thick snow (2007); (4) disintegration of the unstable flanks of a lava dome-like structure growing over the rim of one of the summit craters (1999). All of these recent PDCs were of a rather minor extent (maximum runout lengths were about 1.5 km in November 2006 and March 2007) and thus they represented no threat for populated areas and human property around the volcano. Yet, events of this type pose a significant threat to the lives of people visiting the summit area of Etna, and areas in a radius of 2 km from the summit craters should be off-limits anytime an event capable of producing similar PDCs occurs. The most likely source of further PDCs in the near future is the Southeast Crater, the youngest, most active and most unstable of the four summit craters of Etna, where 6 of the 8 documented recent PDCs originated. It is likely that similar hazards exist in a number of volcanic settings elsewhere, especially at snow- or glacier-covered volcanoes and on volcano slopes strongly affected by hydrothermal alteration.     

Evidence of multiple strain fields beneath the eastern flank of Mt. Etna volcano (Sicily,Italy) deduced from seismic and geodetic data during 2003–2004

We carried out a study of the seismicity and ground deformation occurring on Mt. Etna volcano after the end of the 2002–2003 eruption and before the onset of the 2004–2005 eruption. Data were recorded by the permanent local seismic network run by Istituto Nazionale di Geofisica e Vulcanologia – Sezione di Catania and by geodetic surveys carried out in July 2003 and July 2004 on the GPS network. Most earthquakes were grouped in two main clusters located in the northeastern and southeastern sectors of the volcano. The areal distribution of seismic energy associated with the recorded earthquakes allowed us to highlight the main seismogenic areas of Mt. Etna. In order to better understand the kinematic processes of the volcano, 3D seismic locations were used to compute fault plane solutions, and a selected dataset was inverted to determine stress and strain tensors. The focal mechanisms in the northeastern sector show clear left-lateral kinematics along an E-W fault plane, consistent with events occurring along the Pernicana Fault system. The fault plane solutions in the southeastern sector show mainly right-lateral kinematics along a NNE and ENE fault plane and left lateral-kinematics along NW fault planes that together suggest roughly E-W oriented compression. Surface ground deformation affecting Mt. Etna measured by GPS surveys highlighted a marked inflation during the same period and exceptionally strong seawards motion of its eastern flank. The 2D geodetic strain tensor distribution was calculated and the results show mainly ENE-WSW extension coupled with WNW-ESE contraction, indicating right-lateral shear along a NW-SE oriented fault plane. The different deformation of the eastern sector of the volcano, as measured by seismicity and ground deformation, must be interpreted by considering the different depths of the two signals. Seismic activity in the southeastern sector of volcano is located between 3 and 8 km b.s.l. and can be associated with a very strong additional E-W compression induced by a pressurizing source just westwards and at the same depth, located by inverting GPS data. Ground deformation, in contrast, is mainly affected by the shallower dynamics of the fast moving eastern flank which produces a shallower opposing E-W extension. The entire dataset shows that two different processes affect the eastern flank at the same time but at different depths; the boundary is clearly located at a depth of 3 km b.s.l. and could represent the décollement surface for the mobile flank.     

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.     

火山地区二维自然电位异常反演

将线性反演理论与正则化法相结合，提出一个火山地区二维自然电位（ＳＰ）异常反演方法．地下岩浆活动（对流体）是引起自电异常的极化源，反演的目的是计算地下视极化强度分布和极化源中心位置．可以通过用已知的各种资料改变地下极化源区域（简称源区）的方式来进行多次反演，以便获得一个比较合理的地下视极化强度分布模型．最后，分别给出意大利和墨西哥的两个火山地区的自然电位异常分析结果．     

火山地区二维自然电位异常反演

将线性反演理论与正则化法相结合,提出一个火山地区二维自然电位（ＳＰ）异常反演方法．地下岩浆活动（对流体）是引起自电异常的极化源,反演的目的是计算地下视极化强度分布和极化源中心位置．可以通过用已知的各种资料改变地下极化源区域（简称源区）的方式来进行多次反演,以便获得一个比较合理的地下视极化强度分布模型．最后,分别给出意大利和墨西哥的两个火山地区的自然电位异常分析结果．     

Pyroclastic density currents resulting from the interaction of basaltic magma with hydrothermally altered rock: an example from the 2006 summit eruptions of Mount Etna,Italy

After 16 months of quiescence, Mount Etna began to erupt again in mid-July 2006. The activity was concentrated at and around the Southeast Crater (SEC), one of the four craters on the summit of Etna, and eruptive activity continued intermittently for 5 months. During this period, numerous vents displayed a wide range of eruptive styles at different times. Virtually all explosive activities took place at vents at the summit of the SEC and on its flanks. Eruptive episodes, which lasted from 1 day to 2 weeks, became shorter and more violent with time. Volcanic activity at these vents was often accompanied by dramatic mass-wasting processes such as collapse of parts of the cone, highly unusual flowage processes involving both old rocks and fresh magmatic material, and magma–water interaction. The most dramatic events took place on 16 November, when numerous rockfalls and pyroclastic density currents (PDCs) were generated during the opening of a large fracture on the SE flank of the SEC cone. The largest PDCs were clearly triggered explosively, and there is evidence that much of the energy was generated during the interaction of intruding magma with wet rocks on the cone’s flanks. The most mobile PDCs traveled up to 1 km from their source. This previously unknown process on Etna may not be unique on this volcano and is likely to have taken place on other volcanoes. It represents a newly recognized hazard to those who visit and work in the vicinity of the summit of Etna.     

Evolution of the etna volcano: Information from the southern wall of the Valle Del Bove Caldera

The greater part of Etna can be regarded as a complex strato-shield volcano constructed from the overlapping products of several centres of trachy-basaltic activity. The Valle del Bove is a horse-shoe-shaped caldera, 8 km long and 5 km wide, cut into the eastern flanks of Etna. The caldera is one of the few areas on the volcano where historic eruptions have not obscured the products of pre-historic centres of activity and these are well exposed in the cliff walls surrounding the caldera. Examination of these older volcanics provides important information on the eruptive style and internal plumbing of the Etna volcano during pre-historic times, and suggests that both were significantly different from the present day.Much of the southern wall of the Valle del Bove represents a surviving portion of the Trifoglietto II volcano, the largest pre-historic centre of activity. A stratigraphy is constructed for the southern wall, the Trifoglietto II lavas and pyroclastics rest unconformably upon the eroded remnants of an older centre, and are themselves overlain by the products of younger centres. All the lavas exposed in the southern wall are of alkalic affinity and comprise a trachybasaltic suite ranging from hawaiite to benmoreite. Variation in the chemistry of the lavas can be explained by their differentiation at high levels in the crust from a more basic magma of alkali olivine-basalt/hawaiite composition. An anomalous trend in the TiO₂ content of the Trifoglietto II lavas may be explained by the fractionation of kaersutite (Ti-rich amphibole).A study has been made of the numerous dykes exposed in the walls of the Valle del Bove, the alignments of which parallel trends which are important on Etna at the present time.It is proposed that the initial opening of the Valle del Bove occurred sometime between 20,000 and 10,000 y. B.P., as a result of a phreatic or phreato-magmatic explosion near the base of the eastern flank of Trifoglietto II. This is visualised as triggering a slope failure and resulting in the destruction of much of the centre by a catastrophic landslide. This mechanism has much in common with the explosive eruptions which produced both the Bandai-san (Japan) caldera in 1888, and the Mount St Helens caldera in May, 1980.     

川滇地区7级大震前中强震震源机制变化

分析了 70年代以来 ,川滇地区发生的 8次 7级大震前 5年内 ,发生在大震孕震区和震源区内的中强震震源机制解时空分布。结果表明 ,最早中强震发生在大震震源区或其附近 ,其发震应力场与区域构造应力场一致 ,与大多数大震发震应力场一致或接近。大多数中强震震源破裂特征与大震明显不同。之后有多次中强震发生在距大震震源区较远的大震孕震区内其他地方 ,它们的发震应力场往往经历了与区域构造应力场和大震应力场一致与不一致的多次交替变化。大震前最后 1个中强震也发生在距大震震源区较近的地方 ,其发震应力场与大震发震应力场明显不一致 ,偏转了 30°～ 5 0° ,或更多 ,大多数也与区域构造应力场不一致 ,有的中强震发震断裂破裂特征与大震不一致。大震前中强震震源机制的变化 ,反映了大震孕育过程的不同阶段 ,区域构造应力场的时空调整变化和增强过程 ,以及由此引发的构造断裂异常活动 ,揭示出与大震发生有关的应力场和震源破裂特征信息     

Coseismic and Aseismic Tilt Variations on Mount Etna

?—?In the last ten years (1990–1999), 21 discrete variations of continuous tilt signal have been recorded on Mount Etna, among which one episode was caused by the opening of the eruptive fracture. The remaining 20 anomalies can be classified into two categories: the first comprises 5 “instantaneous” tilt variations recorded in correspondence to the most energetic seismic events (M _L ?≥?3.3) localized on the high western part of the volcano; the second consists of 15 transient anomalies ranging from some hours to 1–2 days, observed at different times at the various tilt stations, with no correlation to seismic events or other evident volcanic episodes. The aseismic variations propagate through the volcanic edifice with a velocity between 4.5–6.0?km/day. Modeling studies suggest that the deformation is generated by a tensile source located 3–6?km SW from Etna volcano summit and 5–10?km depth.     

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.     

A model for radial dike emplacement in composite cones based on observations from Summer Coon volcano,Colorado, USA

We mapped the geometry of 13 silicic dikes at Summer Coon, an eroded Oligocene stratovolcano in southern Colorado, to investigate various characteristics of radial dike emplacement in composite volcanoes. Exposed dikes are up to about 7 km in length and have numerous offset segments along their upper peripheries. Surprisingly, most dikes at Summer Coon increase in thickness with distance from the center of the volcano. Magma pressure in a dike is expected to lessen away from the pressurized source region, which would encourage a blade-like dike to decrease in thickness with distance from the center of the volcano. We attribute the observed thickness pattern as evidence of a driving pressure gradient, which is caused by decreasing host rock shear modulus and horizontal stress, both due to decreasing emplacement depths beneath the sloping flanks of the volcano. Based on data from Summer Coon, we propose that radial dikes originate at depth below the summit of a host volcano and follow steeply inclined paths towards the surface. Near the interface between volcanic cone and basement, which may represent a neutral buoyancy surface or stress barrier, magma is transported subhorizontally and radially away from the center of the volcano in blade-like dikes. The dikes thicken with increasing radial distance, and offset segments and fingers form along the upper peripheries of the intrusions. Eruptions may occur anywhere along the length of the dikes, but the erupted volume will generally be greater for dike-fed eruptions far from the center of the host volcano owing to the increase in driving pressure with distance from the source. Observed eruptive volumes, vent locations, and vent-area intrusions from inferred post-glacial dike-fed eruptions at Mount Adams, Washington, USA, support the proposed model. Hazards associated with radial dike emplacement are therefore greater for longer dikes that propagate to the outer flanks of a volcano.     

Modeling unusual eruptive behavior of Mt. Etna,Italy, by means of event bush

One of the best-studied volcanoes of the world, Mt. Etna in Sicily, repeatedly exhibits eruptive scenarios that depart from the behavior commonly considered typical for this volcano. Episodes of intense explosive activity, pyroclastic flows, dome growth and cone collapse pose a variety of previously underestimated threats to human lives in the summit area of the volcano. However, retrospective analysis of these events shows that they were likely caused by the same very sets of premises and starting conditions as “normal” eruptions, yet combined in an unexpected, probably unique, way. To cope with such unexpected consequences, we involve an approach of artificial intelligence developed specially for needs of the geosciences, the event bush. Scenarios inferred from the event bush fit the observed ones and allow to foresee other low-probability events that may occur at the volcano. Application of the event bush provides a more impartial vision of volcanic phenomena and may serve as an intermediary between expert knowledge and numerical assessment, e.g., by means of Bayesian Belief Networks.