Francesco Guglielmino

Structural assessment of Mount Etna volcano from Permanent Scatterers analysis

A study of the deformation pattern of Mount Etna volcano based on the results from the Permanent Scatterers (PS) technique is reported. Ground motion data provided by the interferometric synthetic aperture radar (InSAR) PS technique from 1995 to 2000 are compared and validated by GPS data. An analysis of the ascending and descending line of sight (LOS) components of ground velocities has yielded detailed ground deformation maps and cross sections. This analysis allows detection and constraint of discontinuities in the surface velocity field. LOS velocities have then been combined to calculate the vertical and horizontal (E-W) ground velocities. A wide inflation of the edifice has been detected on the western and northern flanks (over an area of about 350 km2). A seaward motion of the eastern and southern flanks has also been measured. PS data allows the geometry and kinematics of the several blocks composing the unstable flanks to be defined even in the highly urbanized areas, and their displacement rates have been measured with millimeter precision. This analysis reveals the extension of some features beyond their field evidences and defines new important features. The results of this work depict a new comprehensive kinematic model of the volcano highlighting the gravitational reorganization of the unbuttressed volcanic pile on its slippery clay basement on the southern flank, but an additional drag force due to a strong subsidence of the continental margin facing the Etna volcano is necessary to explain the PS velocity field observed on the eastern flank.

Dynamics of Mount Etna before, during, and after the July–August 2001 eruption inferred from GPS and differential synthetic aperture radar interferometry data

We acknowledge the ‘‘Istituto Nazionale di Geofisica e Vulcanologia’’, the Italian ‘‘Dipartimento per la Protezione Civile’’ and the European Community (contract INGV-DPC UR V3_6/36 and VOLUME Project) for their economic contribution to this research. The SAR data are provided by ESA-ESRIN.

Simultaneous and Integrated Strain Tensor Estimation From Geodetic and Satellite Deformation Measurements to Obtain Three-Dimensional Displacement Maps

We propose a new technique, named SISTEM, based on the elastic theory, to efficiently estimate 3-D displacements for producing deformation maps by integrating sparse Global Positioning System (GPS) measurements of deformations and differential interferometric synthetic aperture radar (DInSAR) maps of movements of the Earths surface. Previous approaches in the literature to combine GPS and DInSAR data require two steps: a first step in which sparse GPS measurements are interpolated in order to fill in GPS displacements in the DInSAR grid and a second step to estimate the 3-D surface displacement maps by using a suitable optimization technique. One of the advantages of the proposed approach, compared to previous ones, is that it does not require the preliminary interpolation of the observed deformation pattern. Indeed, we propose a linear matrix equation which accounts for both the GPS and DInSAR data whose solution simultaneously provides the strain tensor, the displacement field, and the rigid body rotation tensor. The mentioned linear matrix equation is solved by using the weighted least square (WLS), thus assuring both numerical robustness and high computation efficiency. The methodology was tested on both synthetic and experimental data, these last from GPS and DInSAR measurements carried out on Mount Etna during the 2003-2004 period. In order to appreciate the accuracy of the results, the estimated standard errors computed by the WLS are provided. These tests also allow optimizing the choice of specific parameters of this algorithm. This method can be further exploited to account for other available data sets, such as additional interferograms or other geodetic data (e.g., leveling, tilt, etc.), in order to achieve higher accuracy.

Very shallow dyke intrusion and potential slope failure imaged by ground deformation: The 28 December 2014 eruption on Mount Etna

On 28 December 2014, eruptive activity resumed at Mount Etna with fire fountain activity feeding two lava flows on the eastern and southwestern upper flanks of the volcano. Unlike all previous summit activity, this eruption produced clear deformation at the summit of the volcano. GPS displacements and Sentinel-1A ascending interferograms were calculated in order to image the ground deformation pattern accompanying the eruption. The displacements observed by GPS depict a local ground deformation pattern, affecting only the upper part of the volcano. Despite snow cover on the summit, differential interferometry synthetic aperture radar (DInSAR) data allowed obtaining more detail on the ground deformation pattern on the upper eastern side of the volcano. Three-dimensional GPS displacements inversion located a very shallow NE-SW intrusion just beneath the New Southeast Crater. However, this model cannot justify all the deformation observed by DInSAR thus revealing a gravitational failure of the lava flow field.

Decomposing DInSAR Time-Series into 3-D in Combination with GPS in the Case of Low Strain Rates: An Application to the Hyblean Plateau, Sicily, Italy

Differential Interferometric SAR (DInSAR) time-series techniques can be used to derive surface displacement rates with accuracies of 1 mm/year, by measuring the one-dimensional distance change between a satellite and the surface over time. However, the slanted direction of the measurements complicates interpretation of the signal, especially in regions that are subject to multiple deformation processes. The Simultaneous and Integrated Strain Tensor Estimation from Geodetic and Satellite Deformation Measurements (SISTEM) algorithm enables decomposition into a three-dimensional velocity field through joint inversion with GNSS measurements, but has never been applied to interseismic deformation where strain rates are low. Here, we apply SISTEM for the first time to detect tectonic deformation on the Hyblean Foreland Plateau in South-East Sicily. In order to increase the signal-to-noise ratio of the DInSAR data beforehand, we reduce atmospheric InSAR noise using a weather model and combine it with a multi-directional spatial filtering technique. The resultant three-dimensional velocity field allows identification of anthropogenic, as well as tectonic deformation, with sub-centimeter accuracies in areas of sufficient GPS coverage. Our enhanced method allows for a more detailed view of ongoing deformation processes as compared to the single use of either GNSS or DInSAR only and thus is suited to improve assessments of regional seismic hazard.

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.

Gravitational collapse of Mount Etna’s southeastern flank

Gravitational collapse of Mount Etna’s SE flank: New seafloor geodetic data capture active displacement of underwater volcanic flank. The southeastern flank of Etna volcano slides into the Ionian Sea at rates of centimeters per year. The prevailing understanding is that pressurization of the magmatic system, and not gravitational forces, controls flank movement, although this has also been proposed. So far, it has not been possible to separate between these processes, because no data on offshore deformation were available until we conducted the first long-term seafloor displacement monitoring campaign from April 2016 until July 2017. Unprecedented seafloor geodetic data reveal a >4-cm slip along the offshore extension of a fault related to flank kinematics during one 8-day-long event in May 2017, while displacement on land peaked at ~4 cm at the coast. As deformation increases away from the magmatic system, the bulk of Mount Etna’s present continuous deformation must be driven by gravity while being further destabilized by magma dynamics. We cannot exclude flank movement to evolve into catastrophic collapse, implying that Etna’s flank movement poses a much greater hazard than previously thought. The hazard of flank collapse might be underestimated at other coastal and ocean island volcanoes, where the dynamics of submerged flanks are unknown.

Corrigendum: Global positioning system survey data for active seismic and volcanic areas of eastern Sicily, 1994 to 2013

Scientific Data 3:160062 doi: 10.1038/sdata.2016.62 (2016); Published 1 August 2016; Updated 13 September 2016 The authors regret that financial support for this Data Descriptor was not fully acknowledged. The Acknowledgements should have read:

GPS and DInSAR timeseries SISTEM integration for interseismic motion detection — A case study from the Hyblean Plateau in South-East Sicily

Within this study we tested a combination of processing tools in order to overcome certain limitations of Advanced Differential Interferometrie SAR (A-DInSAR) for the detection of intraplate surface motion. The study area is located over the Hyblean Plateau in south-east Sicily, Italy. Instead of usually applied spatio-temporal smoothing of the stacked dataset, we calculated the atmospheric contribution within the interferometrie signal by external data from the global meteorological model ERA-INTERIM. Furthermore we enhanced the coverage of the mean line-of-sight velocity by gap filling and applied a multi-directional filter. For this we used ENVISAT ASAR imagery (49 scenes descending and 58 scenes ascending orbit). The output was then integrated with local GPS data into the SISTEM framework to decompose the measurements into a 3-dimensional velocity field. While in areas of no GPS stations, detected deformation is not precise enough, the results in areas of good GPS coverage look promising with accuracies down to 1 mm/y and various zones of surface deformation could be detected successfully.

Analysis of the SBAS-DInSAR displacement time-series accuracies retrieved in volcanic areas through the first and second generation sensor SAR data

We carry out a quantitative assessment of the accuracy of the advanced Differential SAR Interferometry (DInSAR) approach referred to as Small BAseline Subset (SBAS) in mapping deformation phenomena affecting active volcanic areas, by exploiting SAR data acquired by the “first” and “second” generation SAR sensors, and by comparing the achieved results with independent geodetic measurements, the latter assumed as reference. In particular, we analyze the impact that the different wavelengths and looking geometries may have in the SBAS-DInSAR measurement retrieval depending on the radar system. To this aim, we consider SAR images acquired by the ERS-1/2 and ENVISAT (C-band) as well as COSMO-SkyMed (X-band) sensors (i.e., by the “first” and “second” generation SAR sensors, respectively) over the Campi Flegrei caldera (Southern Italy) and the in-situ displacement measurements provided by the continuous GPS stations deployed in the area. The performed comparative analysis reveals a general good agreement between DInSAR and geodetic data; it also shows that the single displacement measurements have a sub-centimetric accuracy with a standard deviation of about 0.5 cm.