Petrology of lavas from the 2004–2005 flank eruption of Mt. Etna,Italy: inferences on the dynamics of magma in the shallow plumbing system

Following the 2001 and 2002–2003 flank eruptions, activity resumed at Mt. Etna on 7 September 2004 and lasted for about 6 months. This paper presents new petrographic, major and trace element, and Sr–Nd isotope data from sequential samples collected during the entire 2004–2005 eruption. The progressive change of lava composition allowed defining three phases that correspond to different processes controlling magma dynamics inside the central volcano conduits. The compositional variability of products erupted up to 24 September is well reproduced by a fractional crystallization model that involves magma already stored at shallow depth since the 2002–2003 eruption. The progressive mixing of this magma with a distinct new one rising within the central conduits is clearly revealed by the composition of the products erupted from 24 September to 15 October. After 15 October, the contribution from the new magma gradually becomes predominant, and the efficiency of the mixing process ensures the emission of homogeneous products up to the end of the eruption. Our results give insights into the complex conditions of magma storage and evolution in the shallow plumbing system of Mt. Etna during a flank eruption. Furthermore, they confirm that the 2004–2005 activity at Etna was triggered by regional movements of the eastern flank of the volcano. They caused the opening of a complex fracture zone extending ESE which drained a magma stored at shallow depth since the 2002–2003 eruption. This process favored the ascent of a different magma in the central conduits, which began to be erupted on 24 September without any significant change in eruptive style, deformation, and seismicity until the end of eruption.     

The role of the Pernicana Fault System in the spreading of Mt. Etna (Italy) during the 2002–2003 eruption

Flank instability and collapse are observed at many volcanoes. Among these, Mt. Etna is characterized by the spreading of its eastern and southern flanks. The eastern spreading area is bordered to the north by the E–W-trending Pernicana Fault System (PFS). During the 2002–2003 Etna eruption, ground fracturing along the PFS migrated eastward from the NE Rift, to as far as the 18 km distant coastline. The deformation consisted of dextral en-echelon segments, with sinistral and normal kinematics. Both of these components of displacement were one order of magnitude larger (~1 m) in the western, previously known, portion of the PFS with respect to the newly surveyed (~9 km long) eastern section (~0.1 m). This eastern section is located along a pre-existing, but previously unknown, fault, where displaced man-made structures give overall slip rates (1–1.9 cm/year), only slightly lower than those calculated for the western portion (1.4–2.3 cm/year). After an initial rapid motion during the first days of the 2002–2003 eruption, movement of the western portion of the PFS decreased dramatically, while parts of the eastern portion continued to move. These data suggest a model of spreading of the eastern flank of Etna along the PFS, characterized by eruptions along the NE Rift, instantaneous, short-lived, meter-scale displacements along the western PFS and more long-lived centimeter-scale displacements along the eastern PFS. The surface deformation then migrated southwards, reactivating, one after the other, the NNW–SSE-trending Timpe and Trecastagni faults, with displacements of ~0.1 and ~0.04 m, respectively. These structures, along with the PFS, mark the boundaries of two adjacent blocks, moving at different times and rates. The new extent of the PFS and previous activity over its full length indicate that the sliding eastern flank extends well below the Ionian Sea. The clustering of seismic activity above 4 km b.s.l. during the eruption suggests a deep décollement for the moving mass. The collected data thus suggests a significant movement (volume >1,100 km³) of the eastern flank of Etna, both on-shore and off-shore.Editorial responsibility: R. Cioni     

The July–August 2001 eruption of Mt. Etna (Sicily)

The July–August 2001 eruption of Mt. Etna stimulated widespread public and media interest, caused significant damage to tourist facilities, and for several days threatened the town of Nicolosi on the S flank of the volcano. Seven eruptive fissures were active, five on the S flank between 3,050 and 2,100 m altitude, and two on the NE flank between 3,080 and 2,600 m elevation. All produced lava flows over various periods during the eruption, the most voluminous of which reached a length of 6.9 km. Mineralogically, the 2001 lavas fall into two distinct groups, indicating that magma was supplied through two different and largely independent pathways, one extending laterally from the central conduit system through radial fissures, the other being a vertically ascending eccentric dike. Furthermore, one of the eccentric vents, at 2,570 m elevation, was the site of vigorous phreatomagmatic activity as the dike cut through a shallow aquifer, during both the initial and closing stages of the eruption. For 6 days the magma column feeding this vent was more or less effectively sealed from the aquifer, permitting powerful explosive and effusive magmatic activity. While the eruption was characterized by a highly dynamic evolution, complex interactions between some of the eruptive fissures, and changing eruptive styles, its total volume (~25×10 ⁶ m ³ of lava and 5–10×10 ⁶ m ³ of pyroclastics) was relatively small in comparison with other recent eruptions of Etna. Effusion rates were calculated on a daily basis and reached peaks of 14–16 m ³ s ^-1, while the average effusion rate at all fissures was about 11 m ³ s ^-1, which is not exceptionally high. The eruption showed a number of peculiar features, but none of these (except the contemporaneous lateral and eccentric activity) represented a significant deviation from Etna's eruptive behavior in the long term. However, the 2001 eruption could be but the first in a series of flank eruptions, some of which might be more voluminous and hazardous. Placed in a long-term context, the eruption confirms a distinct trend, initiated during the past 50 years, toward higher production rates and more frequent eruptions, which might bring Etna back to similar levels of activity as during the early to mid seventeenth century.     

Lava channel roofing,overflows, breaches and switching: insights from the 2008–2009 eruption of Mt. Etna

During long-lived basaltic eruptions, overflows from lava channels and breaching of channel levées are important processes in the development of extensive 'a'ā lava flow-fields. Short-lived breaches result in inundation of areas adjacent to the main channel. However, if a breach remains open, lava supply to the original flow front is significantly reduced, and flow-field widening is favoured over lengthening. The development of channel breaches and overflows can therefore exert strong control over the overall flow-field development, but the processes that determine their location and frequency are currently poorly understood. During the final month of the 2008–2009 eruption of Mt. Etna, Sicily, a remote time-lapse camera was deployed to monitor events in a proximal region of a small ephemeral lava flow. For over a period of ~10 h, the flow underwent changes in surface elevation and velocity, repeated overflows of varying vigour and the construction of a channel roof (a required prelude to lava tube formation). Quantitative interpretation of the image sequence was facilitated by a 3D model of the scene constructed using structure-from-motion computer vision techniques. As surface activity waned during the roofing process, overflow sites retreated up the flow towards the vent, and eventually, a new flow was initiated. Our observations and measurements indicate that flow surface stagnation and flow inflation propagated up-flow at an effective rate of ~6 m h⁻¹, and that these processes, rather than effusion rate variations, were ultimately responsible for the most vigorous overflow events. We discuss evidence for similar controls during levée breaching and channel switching events on much larger flows on Etna, such as during the 2001 eruption.     

Regimes of magma recharge and their control on the eruptive behaviour during the period 2001–2005 at Mt. Etna volcano

Mount Etna volcano (Italy) during the period 2001–2005 has undergone a period of intense eruptive activity marked by three large eruptions (2001, 2002–2003 and 2004–2005). These eruptions encompassed diverse eruptive styles and regimes: from intensely explosive, during 2001 and 2002–2003 eruptions, to exclusively effusive in the 2004–2005 event. In this work, we put forward the idea that these three eruptions are the response of the progressive arrival into the uppermost segment of the open-conduit system of a new magma, which was geochemically distinct in terms of trace element and Sr–Nd–Pb isotope signature from the products previously emitted by the Etnean volcano. The magma migrated upwards mainly through a peripheral tectonic system, which can be considered as eccentric in spite of its relative proximity to the main system. The ingress of the new magma and its gradual displacement from the eccentric system into the uppermost sector of the open-conduit gave rise to different eruptive behaviours. At the beginning, the ascent of the undegassed magma, able to exsolve a gas phase at depth, and its interaction with closed-system magma reservoirs less than 10 km deep gave rise to the explosive events of 2001 and 2002–2003. Later, when the same magma entered into the open-conduit system, it took part in the steady-state degassing and partially lost its volatile load, leading to a totally effusive eruption during the 2004–2005 event. One further consideration highlighted here is that in 2001–2005, migration of the feeding axis from an eccentric and peripheral position towards the main open-conduit has led to the development of a new vent (South East Crater 2) located at the eastern base of the South East Crater through which most of the subsequent Etnean activity occurred.     

Damage mechanics preceding the September–October 1989 flank eruption at Mount Etna volcano inferred by seismic scaling exponents

Seismicity preceding the September–October 1989 flank eruption at Mt. Etna volcano has been analysed by means of the fractal two-point correlation dimension D_c relative to epicentres distribution and the scaling exponent D inferred from the seismic b-value. Time evolution of D_c and D preceding the eruption onset revealed: a general decrease of D_c and increase of D (9–4 months prior the eruption onset), intermittent rates of D_c increase and D decrease (4 to 1 month), steady high values of D_c and low values of D (1 month) and sharp decreases of D_c and D marking the eruption onset. Results allowed to constrain the increment of seismic damage by recognising negative (mechanical hardening) and positive (mechanical softening) feedback prior the eruption onset.     

SO2 emissions at Mt. Etna with particular reference to the period 1993–1995

The variations in sulfur dioxide (SO₂) emission from the Summit Craters of Mt. Etna were determined, with particular reference to the period 1993–1995. Vehicle-based weekly measurements of SO₂ flux, using a correlation spectrometer (COSPEC), suggest new input of magma into the main feeder system of the volcano between 1993 and 1995. Minimal flux values (<1000 t/day) preceded the two eruptive events in the period 1987–1995. Only approximately 9.5% of the magma that contributed the SO₂ emission was erupted during the same period. Received: 3 November 1997 / Accepted: 21 September 1998     

Lava fountains during the episodic eruption of South–East Crater (Mt. Etna), 2000: insights into magma-gas dynamics within the shallow volcano plumbing system

Mt. Etna, in Sicily (Italy) is well known for frequent effusive and explosive eruptions from both its summit and flanks. South-East Crater (SE Crater), one of the four summit craters, has been the most active in the last 20 years and often produces episodic lava fountains over periods lasting from a few weeks to months. The most striking of such eruptive phases was in 2000. Sixty four lava fountains, separated by quiescent intervals and sometimes associated with lava overflows, occurred that year between January and June, a time period during which we consider the volcano to have been in episodic eruption. This paper presents mainly results of petrochemical investigations carried out on both tephra and lavas collected during a number of the lava fountain episodes in 2000. The new data have been integrated with volcanological and seismic information in order to correlate the features of the eruptive activity with magma-gas dynamics in the plumbing system of SE Crater. The main findings allow us to characterise the 2000 episodic eruption in the framework of the recent SE Crater activity. In particular, we infer that the onset of the 2000 eruption was triggered by the ascent of new, more primitive and volatile-rich magma that progressively intruded into the SE Crater reservoir, where it mixed with the resident, more evolved magma. Furthermore, we argue that the 2000 SE Crater lava fountains largely resulted from the instability of a foam layer accumulated at the top of the underlying reservoir and rebuilt prior to each episode, in agreement with the collapsing foam model for lava fountains.     

Representivity of incompletely sampled fall deposits in estimating eruption source parameters: a test using the 12–13 January 2011 lava fountain deposit from Mt. Etna volcano,Italy

The Southeast Crater (SEC) of Mt. Etna, Italy, is renowned for its high activity, mainly long-lived eruptions consisting of sequences of individual paroxysmal episodes which have produced more than 150 eruptive events since 1998. Each episode typically forms eruption columns followed by tephra fallout over distances of up to about 100 km from the vent. One of the last sequences consisted of 25 lava fountaining events, which took place between January 2011 and April 2012 from a pit-vent on the eastern flank of the SEC and built a new scoria cone renamed New Southeast Crater. The first episode on 12–13 January 2011 produced tephra fallout which was unusually dispersed toward to the South extending out over the Mediterranean Sea. The southerly deposition of tephra permitted an extensive survey at distances between ~1 and ~100 km, providing an excellent characterization of the tephra deposit. Here, we document the stratigraphy of the 12–13 January fallout deposit, draw its dispersal, and reconstruct its isopleth map. These data are then used to estimate the main eruption source parameters. The total erupted mass (TEM) was calculated by using four different methodologies which give a mean value of 1.5?±?0.4?×?10⁸ kg. The mass eruption rate (MER) is 2.5?±?0.7?×?10⁴ kg/s using eruption duration of 100 min. The total grain-size (TGS) distribution, peaked at ?3 phi, ranges between ?5 and 5 phi and has a median value of ?1.4 phi. Further, for the eruption column height, we obtained respective values of 6.8–13.8 km by using the method of Carey and Sparks (1986) and 3.4?±?0.3 km by using the methods of Wilson and Walker (1987), Mastin et al. (2009), and Pistolesi et al. (2011) and considering the mean value of MER from the deposit. We also evaluated the uncertainty and reliability of TEM and TGS for scenarios where the proximal and distal samples are not obtainable. This is achieved by only using a sector spanning the downwind distances between 6 and 23 km. This scenario is typical for Etna when the tephra plume is dispersed eastward, i.e., in the prevailing wind direction. Our results show that, if the analyzed deposit has poorer sample coverage than presented in this study, the TEM (3.4?×?10⁷ kg) is 22 % than the TEM obtained from the whole deposit. The lack of the proximal (<6 km) deposit may cause more significant differences in the TGS estimations.     

Reconstructing the eruption magnitude and energy budgets for the pre-historic eruption of the monogenetic ∼5 ka Mt. Gambier Volcanic Complex, south-eastern Australia

Understanding explosive volcanic eruptions, especially phreatomagmatic eruptions, their intensities and energy budgets is of major importance when it comes to risk and hazard studies. With only a few historic occurrences of phreatomagmatic activity, a large amount of our understanding comes from the study of pre-historic volcanic centres, which causes issues when it comes to preservation and vegetation. In this research, we show that using 3D geometrical modelling it is possible to obtain volume estimates for different deposits of a pre-historic, complex, monogenetic centre, the Mt. Gambier Volcanic Complex, south-eastern Australia. Using these volumes, we further explore the energy budgets and the magnitude of this eruption (VEI 4), including dispersal patterns (eruption columns varying between 5 and 10 km, dispersed towards north-east to south), to further our understanding of intraplate, monogenetic eruptions involving phreatomagmatic activity. We also compare which thermodynamic model fits best in the creation of the maar crater of Mt. Gambier: the major-explosion-dominated model or the incremental growth model. In this case, the formation of most of the craters can best be explained by the latter model.     

The November 2002 eruption of Piton de la Fournaise,Réunion: tracking the pre-eruptive thermal evolution of magma using melt inclusions

The November 2002 eruption of Piton de la Fournaise in the Indian Ocean was typical of the activity of the volcano from 1999 to 2006 in terms of duration and volume of magma ejected. The first magma erupted was a basaltic liquid with a small proportion of olivine phenocrysts (Fo₈₁) that contain small numbers of melt inclusions. In subsequent flows, olivine crystals were more abundant and richer in Mg (Fo_83–84). These crystals contain numerous melt and fluid inclusions, healed fractures, and dislocation features such as kink bands. The major element composition of melt inclusions in this later olivine (Fo_83–84) is out of equilibrium with that of its host as a result of extensive post-entrapment crystallization and Fe²⁺ loss by diffusion during cooling. Melt inclusions in Fo₈₁ olivine are also chemically out of equilibrium with their hosts but to a lesser degree. Using olivine–melt geothermometry, we determined that melt inclusions in Fo₈₁ olivine were trapped at lower temperature (1,182 ± 1°C) than inclusions in Fo_83–84 olivine (1,199–1,227°C). This methodology was also used to estimate eruption temperatures. The November 2002 melt inclusion compositions suggest that they were at temperatures between 1,070°C and 1,133°C immediately before eruption and quenching. This relatively wide temperature range may reflect the fact that most of the melt inclusions were from olivine in lava samples and therefore likely underwent minor but variable amounts of post-eruptive crystallization and Fe²⁺ loss by diffusion due to their relatively slow cooling on the surface. In contrast, melt inclusions in tephra samples from past major eruptions yielded a narrower range of higher eruption temperatures (1,163–1,181°C). The melt inclusion data presented here and in earlier publications are consistent with a model of magma recharge from depth during major eruptions, followed by storage, cooling, and crystallization at shallow levels prior to expulsion during events similar in magnitude to the relatively small November 2002 eruption.     

Magmatic plumbing systems of the Koryakskii–Avacha Volcanic Cluster as inferred from observations of local seismicity and from the regime of adjacent thermal springs

An analysis of local seismicity within the Avacha–Koryakskii Volcanic Cluster during the 2000–2016 period revealed a sequence of plane-oriented earthquake clusters that we interpret as a process of dike and sill emplacement. The highest magmatic activity occurred in timing with the 2008–2009 steam–gas eruption of Koryakskii Volcano, with magma injection moving afterwards into the cone of Avacha Volcano (2010–2016). The geometry of the magma bodies reflects the NF geomechanical conditions (tension and normal faults, \(S_V>S_{H_{\text{max}}}>S_{h_{\text{min}}}\)) at the basement of Koryakskii Volcano dominated by vertical stresses S _v , with the maximum horizontal stress \(S_{h_{\text{max}}}\) pointing north. A CFRAC simulation of magma injection into a fissure under conditions that are typical of those in the basement of Koryakskii Volcano (the angle of dip is 60°, the size is 2 × 2 km², and the depth is –4 km abs.) showed that when the magma discharge is maintained at the level of 20000 kg/s during 24 hours the fissure separation increases to reach 0.3 m and the magma injection is accompanied by shear movements that occur at a rate as high as 2 × 10^–3 m/s, thus corresponding to the conditions of local seismic events with Mw below 4.5. We are thus able to conclude that the use of planeoriented clusters of earthquakes for identification of magma emplacement events is a physically sound procedure. The August 2, 2011 seismicity increase in the area of the Izotovskii hot spring (7 km from the summit of Koryakskii Volcano), which is interpreted as the emplacement of a dike, has been confirmed by an increase in the spring temperature by 10–12°C during the period from October 2011 to July 2012.     

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.     

Comparison of Integrated Geodetic Data Models and Satellite Radar Interferograms to Infer Magma Storage During the 1991–1993 Mt. Etna Eruption

We compare the results obtained from the modelling of EDM, GPS, levelling and tilt data measured in the first part of the 1991–1993 eruption at Etna to the InSAR data acquired during the second part. The geodetic changes are very marked in the first half of the eruption and constrain a deflation source located at a few kilometers of depth ( 3 km b.s.l.), in agreement with other independent geophysical evidence. SAR data, available during the second part of the eruption, were analysed for different time intervals in the second part of the eruption. The interpretation of SAR interferograms reveals a large-scale but less marked deflation of the volcano that could be caused by a deeper source. This second source is in accord with a second deeper anomaly revealed by recent seismic investigations. The combination of geodetic data modelling and SAR images suggests a complex plumbing system composed at least of two possible storage regions located at different depths.     

Pre-eruptive physical conditions of El Reventador volcano (Ecuador) inferred from the petrology of the 2002 and 2004–05 eruptions

A petrological study of the eruptive products of El Reventador allowed us to infer the magmatic processes related to the 2002 and 2004–05 eruptions of this andesitic stratovolcano. On November 3, 2002, El Reventador experienced a highly explosive event, which was followed by emplacement of two lava flows in November–December 2002. Silica contents range from 62 to 58 wt.% SiO₂ for the November 3 pyroclastic deposits to 58–56 and 54–53 wt.% SiO₂ for the successive lava flows. In November 2004 eruptive activity resumed supplying four new lava flows (56–54 wt.% SiO₂) between November 2004 and August 2005.     

New Insights on Mt. Etna’s Crust and Relationship with the Regional Tectonic Framework from Joint Active and Passive P-Wave Seismic Tomography

In the Central Mediterranean region, the production of chemically diverse volcanic products (e.g., those from Mt. Etna and the Aeolian Islands archipelago) testifies to the complexity of the tectonic and geodynamic setting. Despite the large number of studies that have focused on this area, the relationships among volcanism, tectonics, magma ascent, and geodynamic processes remain poorly understood. We present a tomographic inversion of P-wave velocity using active and passive sources. Seismic signals were recorded using both temporary on-land and ocean bottom seismometers and data from a permanent local seismic network consisting of 267 seismic stations. Active seismic signals were generated using air gun shots mounted on the Spanish Oceanographic Vessel ‘Sarmiento de Gamboa’. Passive seismic sources were obtained from 452 local earthquakes recorded over a 4-month period. In total, 184,797 active P-phase and 11,802 passive P-phase first arrivals were inverted to provide three different velocity models. Our results include the first crustal seismic active tomography for the northern Sicily area, including the Peloritan–southern Calabria region and both the Mt. Etna and Aeolian volcanic environments. The tomographic images provide a detailed and complete regional seismotectonic framework and highlight a spatially heterogeneous tectonic regime, which is consistent with and extends the findings of previous models. One of our most significant results was a tomographic map extending to 14 km depth showing a discontinuity striking roughly NW–SE, extending from the Gulf of Patti to the Ionian Sea, south-east of Capo Taormina, corresponding to the Aeolian–Tindari–Letojanni fault system, a regional deformation belt. Moreover, for the first time, we observed a high-velocity anomaly located in the south-eastern sector of the Mt. Etna region, offshore of the Timpe area, which is compatible with the plumbing system of an ancient shield volcano located offshore of Mt. Etna.