Susi Pepe

Deformation and eruptions at Mt. Etna (Italy): A lesson from 15 years of observations

This work was partly funded by INGV and the Italian DPC and was supported by ASI, the Preview Project and CRdC-AMRA. DPC-INGV Flank project providing the funds for the publication fees.

Magma injection beneath the urban area of Naples: a new mechanism for the 2012–2013 volcanic unrest at Campi Flegrei caldera

We found the first evidence, in the last 30 years, of a renewed magmatic activity at Campi Flegrei caldera from January 2012 to June 2013. The ground deformation, observed through satellite interferometry and GPS measurements, have been interpreted as the effect of the intrusion at shallow depth (3090 ± 138 m) of 0.0042 ± 0.0002 km3 of magma within a sill. This interrupts about 28 years of dominant hydrothermal activity and occurs in the context of an unrest phase which began in 2005 and within a more general ground uplift that goes on since 1950. This discovery has implications on the evaluation of the volcanic risk and in the volcanic surveillance of this densely populated area.

Ground deformation and source geometry of the 24 August 2016 Amatrice earthquake (Central Italy) investigated through analytical and numerical modeling of DInSAR measurements and structural-geological data

We investigate the ground deformation and source geometry of the 2016 Amatrice earthquake (Central Italy) by exploiting ALOS2 and Sentinel-1 coseismic differential interferometric synthetic aperture radar (DInSAR) measurements. They reveal two NNW-SSE striking surface deformation lobes, which could be the effect of two distinct faults or the rupture propagation of a single fault. We examine both cases through a single and a double dislocation planar source. Subsequently, we extend our analysis by applying a 3-D finite elements approach jointly exploiting DInSAR measurements and an independent, structurally constrained, 3-D fault model. This model is based on a double fault system including the two northern Gorzano and Redentore-Vettoretto faults (NGF and RVF) which merge into a single WSW dipping fault surface at the hypocentral depth (8 km). The retrieved best fit coseismic surface deformation pattern well supports the exploited structural model. The maximum displacements occur at 5–7 km depth, reaching 90 cm on the RVF footwall and 80 cm on the NGF hanging wall. The von Mises stress field confirms the retrieved seismogenic scenario.

Detachment depth revealed by rollover deformation: An integrated approach at Mount Etna

This work was partially funded by INGV and the DPC‐INGV project “Flank”, and partially by the ASI (SRV project).

Capturing the fingerprint of Etna volcano activity in gravity and satellite radar data

Long-term and high temporal resolution gravity and deformation data move us toward a better understanding of the behavior of Mt Etna during the June 1995 – December 2011 period in which the volcano exhibited magma charging phases, flank eruptions and summit crater activity. Monthly repeated gravity measurements were coupled with deformation time series using the Differential Synthetic Aperture Radar Interferometry (DInSAR) technique on two sequences of interferograms from ERS/ENVISAT and COSMO-SkyMed satellites. Combining spatiotemporal gravity and DInSAR observations provides the signature of three underlying processes at Etna: (i) magma accumulation in intermediate storage zones, (ii) magmatic intrusions at shallow depth in the South Rift area, and (iii) the seaward sliding of the volcanos eastern flank. Here we demonstrate the strength of the complementary gravity and DInSAR analysis in discerning among different processes and, thus, in detecting deep magma uprising in months to years before the onset of a new Etna eruption.

Evidence of a shallow persistent magmatic reservoir from joint inversion of gravity and ground deformation data: The 25–26 October 2013 Etna lava fountaining event

To evaluate the volcanic processes leading to the 25–26 October 2013 lava fountain at Mount Etna, we jointly investigated gravity, GPS, and DInSARmeasurements covering the late-June to early-November time interval. We used finite element modeling to infer a shallow magmatic reservoir which (i) inflated since July 2013, (ii) fed the volcanic activity at the summit craters during 25–26 October, and (iii) deflated due to magma drainage related to this volcanic activity. We suggested that this reservoir belongs to a shallow volume, which is located beneath the summit area and is replenished by magma rising from deep reservoirs and fed the short-term volcanic activity, representing a persistent shallow magmatic plumbing system of Etna. In addition, the model results show that there is a large discrepancy between the erupted and shallow reservoir deflation volumes, which could be reasonably attributable to a highly compressible volatile-rich magma.

SBAS-DInSAR time series in the last eighteen years at Mt. Etna volcano (Italy)

We investigate the deformation of Mt. Etna volcano (Italy) by exploiting the advanced Differential Synthetic Aperture Radar Interferometry (DInSAR) technique referred to as the Small BAseline Subset (SBAS) algorithm. In particular, we take advantage of the multi-sensor data processing capability of the SBAS algorithm which allows us to generate Mt. Etna mean deformation velocity maps and the corresponding time series in the last eighteen years. To achieve this task we exploit different set of SAR data collected by the European (ERS-1/2, ENVISAT) satellites in the 1992–2010 time interval, and by the Italian COSMO-SkyMed constellation during 2009–2010 period. We also benefit from the availability of ERS-ENVISAT multi-orbit (ascending and descending) data in order to discriminate the vertical and East-West components of the volcano edifice displacements and generate the relevant time series. Finally, we evidence how the higher spatial resolution and denser temporal sampling of the COSMO-SkyMed data, with respect to the European satellites, permit to follow with more details the complex deformative pattern which has characterized the volcano in the last two years.

A quantitative assessment of DInSAR Time series accuracy in volcanic areas: From the first to second generation SAR sensors

We perform a quantitative assessment of the accuracy of Differential SAR Interferometry (DInSAR) time series in volcanic areas, retrieved through “first” and “second generation” SAR data. In particular, we analyze the impact that the wavelengths and looking geometries may have in the DInSAR measurement retrieval depending on the radar system. To this aim, we focus on the DInSAR algorithm referred to as Small BAseline Subset (SBAS) to generate mean deformation velocity maps and corresponding time series starting from sequences of SAR images. Moreover, we consider collections of SAR data acquired by the ERS-1/2 and ENVISAT (C-band), and COSMO-SkyMed (Xband) sensors over the volcanic area of the Campi Flegrei caldera, Southern Italy. We invert these SAR data sequences through the SBAS-DInSAR technique, thus obtaining C- and X- band deformation time series that we compare to continuous GPS measurements, the latter assumed as reference. The achieved results provide, in addition to a clear picture of the surface deformation phenomena already occurred and occurring in the selected case study, relevant indications for the analysis of the SBAS-DInSAR time series accuracies in volcanic areas passing from the first to second generation SAR sensors.

The SBAS-DInSAR technique as a tool for the observation of active volcanic areas: Results and future perspectives

In this work we describe the application of the basic small BAseline subset (SBAS) technique, which exploits multilook interferograms, to a number of active volcanic areas. The use of such a technique reduces the amount of data to be processed and simplifies the analysis of geophysical phenomena occurring in extended areas (up to 100 km by 100 km). Moreover, it can be trivially extended in order to combine data acquired by different sensors with similar geometrical and electromagnetic characteristics, i.e., ERS and ENVISAT IS-2 mode, thus allowing an easy extension of the temporal observation window. The selected test sites for this study are Long Valley caldera, Mt. Etna and the Neapolitan Volcanic district (Mt. Vesuvio and Campi Flegrei caldera). Finally, we shortly analyze the implications of the use of forthcoming SAR sensors that operates in L-and X-bands.

The Use of Massive Deformation Datasets for the Analysis of Spatial and Temporal Evolution of Mauna Loa Volcano (Hawai’i)

In this work, we exploited large DInSAR and GPS datasets to create a 4D image of the magma transfer processes at Mauna Loa Volcano (Island of Hawai’i) from 2005 to 2015. The datasets consist of 23 continuous GPS time series and 307 SAR images acquired from ascending and descending orbits by ENVISAT (ENV) and COSMO-SkyMed (CSK) satellites. Our results highlight how the joint use of SAR data acquired from different orbits (thus with different look angles and wavelengths), together with deformation data from GPS networks and geological information can significantly improve the constraints on the geometry and location of the sources responsible for the observed deformation. The analysis of these datasets has been performed by using an innovative method that allows building a complex source configuration. The results suggest that the deformation pattern observed from 2005 to 2015 has been controlled by three deformation sources: the ascent of magma along a conduit, the opening of a dike and the slip along the basal decollement. This confirms that the intrusion of the magma within a tabular system (rift dikes) may trigger the sliding of the SE portion of the volcanic edifice along the basal decollement. This case study confirms that it is now possible to exploit large geodetic datasets to improve our knowledge of volcano dynamics. The same approach could also be easily applied in other geodynamical contexts such as geothermal reservoirs and regions with complex tectonics.