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.

We have modeled the fast dike intrusion that started on 24 December 2018 at Mount Etna. The intrusion was accompanied by an intense seismicity swarm that also continued the following day. Since previous studies did not detail the overall chain of events in time during the magma ascent, here we propose a combined analytical and FEM modeling of all available continuous deformation data, focusing on the signals over 2 days (24–25 December) when the continuous deformation networks recorded clear variations directly related to the dike ascent. High‐rate GPS enabled obtaining an early and reliable source model. Borehole instruments (strainmeters and tiltmeters) highlighted clear variations, starting about 1 hr before those recorded by GPS, and moreover, made it possible to improve the dike ascent modeling. In particular, our continuous deformation data clearly revealed not one but two dikes and showed how they evolved in time. We inferred one dike that began propagating in the first kilometer above sea level, continued to rise with a maximum opening of 1.9 m and increased its horizontal dimension until reaching the ground surface. Soon after, continuous deformation networks revealed a new elongated intrusion in the southern flank, matching the south shifted position of the seismic swarm. This second dike, with a thicker opening of 4.9 m, started from a depth of about −3 km (below sea level) but did not reach the ground surface. This proposed multiparametric modeling of continuous deformation data has therefore enabled disentangling the complexity of the real volcanic processes.

We demonstrate the capability of distributed acoustic sensing (DAS) to record volcano-related dynamic strain at Etna (Italy). In summer 2019, we gathered DAS measurements from a 1.5 km long fibre in a shallow trench and seismic records from a conventional dense array comprised of 26 broadband sensors that was deployed in Piano delle Concazze close to the summit area. Etna activity during the acquisition period gives the extraordinary opportunity to record dynamic strain changes (∼ 10⁻⁸ strain) in correspondence with volcanic events. To validate the DAS strain measurements, we explore array-derived methods to estimate strain changes from the seismic signals and to compare with strain DAS signals. A general good agreement is found between array-derived strain and DAS measurements along the fibre optic cable. Short wavelength discrepancies correspond with fault zones, showing the potential of DAS for mapping local perturbations of the strain field and thus site effect due to small-scale heterogeneities in volcanic settings.

Between 13 December 2020 and 21 February 2022, Etna volcano produced a sequence of 66 paroxysmal explosive eruptions, with Strombolian activity at the summit craters climaxing in lava fountains and eruption columns extending several kilometers above the craters, accompanied by minor and short-lasting lava flows from the crater rim. We selected three of these episodes that occurred within a short space of time, between 13 December 2020 and 12 March 2021, of different magnitude (i.e., erupted volume) and intensity (i.e., mass eruption rate or instantaneous eruption rate), and analyzed them from a multidisciplinary perspective. The aim was to gain insights into those parameters that mostly reveal the eruptive process for hazard assessment purposes. The multidisciplinary data consist of calibrated visible images, thermal images, seismic and infrasound data, ground deformation detected from the strainmeters, as well as satellite SEVIRI images. From these data, we obtained the timing of each paroxysmal event, the erupted volume in terms of tephra and lava flows, and the corresponding deflation of the source region, together with the development of the lava fountains and eruption columns with time. The results enabled determining that the smallest episode was that of 13 December 2020, which comprised three distinctive pulses but did not produce an eruptive column detectable from either monitoring cameras or satellites. The 28 February 2021 episode was remarkable for the short amount of time required to reach the climax, and was the most intense, whereas the 12 March 2021 event showed the longest duration but with an intensity between that of the previous two. Our results show that these three paroxysmal events display a typical trend, with the first event also being the smallest in terms of both erupted volume and intensity, the second being the most intense, and the third the one of greatest magnitude but less intense than the second. This is coherent with the end of the first paroxysmal phase on 1 April 2021, which was followed by 48 days of eruptive pause before starting again. In this context, the end of the paroxysmal phase was anticipated by a more effusive episode, thus heralding a temporary decline in the gas content within the feeding magma batch.