3D seismic velocity models from local earthquake tomography furnish new insights on the Mount Etna volcano (Southern Italy)

We present a new seismotomography investigation providing a 3-D overall model of Vp, Vs and Vp/Vs for Mt. Etna, the largest and most active volcano in Europe. We estimated and jointly evaluated P- and S-wave velocity patterns together with the Vp/Vs ratio, particularly useful to discriminate the presence of groundwater, gas, and melts and thus very precious for volcano investigations. We applied the LOTOS software to ~ 4600 crustal earthquakes that occurred in the Etnean area during the last 26 years, the longest time-interval ever analysed for Mt. Etna. This wide dataset has allowed us to characterize the volcano velocity structure getting over possible singularities due to specific eruptive phases. Our results further refined the high velocity body widely recognized in the south-eastern sector of Mt. Etna by furnishing new clues on the possible former magma pathways. Moreover, the obtained 3D seismic velocity model depicted new anomalies revealing the presence of: (i) two shallow underground aquifers in the northern Etnean sector; (ii) a volume of strongly fractured rocks filled of fluids along the eastern flank; (iii) a quite deep region of probable fluid accumulation apparently not linked to the volcanic activity in the western sector.

Fibre optic distributed acoustic sensing of volcanic events

Understanding physical processes prior to and during volcanic eruptions has improved significantly in recent years. However, uncertainties about subsurface structures distorting observed signals and undetected processes within the volcano prevent volcanologists to infer subtle triggering mechanisms of volcanic phenomena. Here, we demonstrate that distributed acoustic sensing (DAS) with optical fibres allows us to identify volcanic events remotely and image hidden near-surface volcanic structural features. We detect and characterize strain signals associated with explosions and locate their origin using a 2D-template matching between picked and theoretical wave arrival times. We find evidence for non-linear grain interactions in a scoria layer of spatially variable thickness. We demonstrate that wavefield separation allows us to incrementally investigate the ground response to various excitation mechanisms. We identify very small volcanic events, which we relate to fluid migration and degassing. Those results provide the basis for improved volcano monitoring and hazard assessment using DAS.

Intense overpressurization at basaltic open-conduit volcanoes as inferred by geochemical signals: The case of the Mt. Etna December 2018 eruption

The balance between the amount of gas coexisting with mantle-derived magmas at depth and that emitted during intereruptive phases may play a key role in the eruptive potential of volcanoes. Taking the December 2018 eruption at Mt. Etna volcano as a case study, we discuss the geochemical data streams observed. The signals indicate a long-lasting prelude stage to eruption, starting in 2017 and involving magma-fluid accumulation in the deep (>7 km bsl) reservoir, followed by pressure buildup in the system at intermediate depth (5 to 2 km bsl), 6 to 7 months before the eruption. A brief preeruptive phase marks the pressurization at 2 to 3 km below the craters. By comparing the magma and fluid recharge at depth to the measured volcanic degassing from the plume, we provide evidence that Mt. Etna was in a state of extreme overpressurization in the weeks before the onset of the eruption.

2D Synthetic dataset of numerical simulations of long-period seismicity in a volcanic edifice and related sensitivity kernels

This work describes the data used in the EPSL research article “Quantifying strong seismic propagation effects in the upper volcanic edifice using sensitivity kernels”. The dataset is generated in order to investigate to what extent the seismic signals recorded on volcanoes are affected by near surface velocity structure. Data were calculated using the computational spectral elements scheme SPECFEM2D, where the wave propagation beneath Mount Etna volcano, Italy, was simulated in both homogeneous and heterogeneous models. The heterogeneous model comprises a low-velocity superficial structure (top several hundred meters) based on the previously published studies. Several different source mechanisms and locations were used in the simulations. The seismic wavefield was “recorded” by 15 surface receivers distributed along the surface of the volcano. The associated sensitivity kernels were also computed. These kernels highlight the region of the velocity model that affects the recorded seismogram within a desired time window. The text files describing the velocity models used in the simulations are also provided. The data may be of interest to volcano seismologists, as well as earthquake seismologists studying path effects and wave propagation through complex media.

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.

The Campi Flegrei caldera: unrest mechanisms and hazards

In the last four decades, Campi Flegrei caldera has been the world’s most active caldera characterized by intense unrest episodes involving huge ground deformation and seismicity, but, at the time of writing, has not culminated in an eruption. We present a careful review, with new analyses and interpretation, of all the data and recent research results. We deal with three main problems: the tentative reconstruction of the substructure; the modelling of unrest episodes to shed light on possible pre-eruptive scenarios; and the probabilistic estimation of the hazards from explosive pyroclastic products. The results show, for the first time at a volcano, that a very peculiar mechanism is generating episodes of unrest, involving mainly activation of the geothermal system from deeper magma reservoirs. The character and evolution of unrest episodes is strongly controlled by structural features, like the ring-fault system at the borders of the caldera collapse. The use of detailed volcanological, mathematical and statistical procedures also make it possible to obtain a detailed picture of eruptive hazards in the whole Neapolitan area. The complex behaviour of this caldera, involving interaction between magmatic and geothermal phenomena, sheds light on the dynamics of the most dangerous types of volcanoes in the world.

Magma ascent and the pressurization of Mount Etna's volcanic system

After a period of deflation during the 1991-1993 flank eruption, Mount Etna underwent a rapid inflation. Seismicity and ground deformation show that since 1994, a huge volume of magma intruded beneath the volcano, producing from 1998 onward a series of eruptions at the summit and on the flank of the volcano. The last of these, started on 27 October 2002, is still in progress and can be considered one of the most explosive eruptions of the volcano in recent times. Here we show how geodetic data and seismic deformation, between 1994 and 2001, indicate a radial compression around an axial intrusion, consistent with a repressurization of Mount Etna's plumbing system at a depth of 6 to 15 kilometers, which triggered most of the seismicity and provoked the dilatation of the volcano and the recent explosive eruptive activity.