We present an advanced differential synthetic aperture radar (SAR) interferometry (DInSAR) processing chain, based on the Parallel Small BAseline Subset (P-SBAS) technique, for the efficient generation of deformation time series from Sentinel-1 (S-1) interferometric wide (IW) swath SAR data sets. We first discuss an effective solution for the generation of high-quality interferograms, which properly accounts for the peculiarities of the terrain observation with progressive scans (TOPS) acquisition mode used to collect S-1 IW SAR data. These data characteristics are also properly accounted within the developed processing chain, taking full advantage from the burst partitioning. Indeed, such data structure represents a key element in the proposed P-SBAS implementation of the S-1 IW processing chain, whose migration into a cloud computing (CC) environment is also envisaged. An extensive experimental analysis, which allows us to assess the quality of the obtained interferometric products, is presented. To do this, we apply the developed S-1 IW P-SBAS processing chain to the overall archive acquired from descending orbits during the March 2015-April 2017 time span over the whole Italian territory, consisting in 2740 S-1 slices. In particular, the quality of the final results is assessed through a large-scale comparison with the GPS measurements relevant to nearly 500 stations. The mean standard deviation value of the differences between the DInSAR and the GPS time series (projected in the radar line of sight) is less than 0.5 cm, thus confirming the effectiveness of the implemented solution. Finally, a discussion about the performance achieved by migrating the developed processing chain within the Amazon Web Services CC environment is addressed, highlighting that a two-year data set relevant to a standard S-1 IW slice can be reliably processed in about 30 h.The presented results demonstrate the capability of the implemented P-SBAS approach to efficiently and effectively process large S-1 IW data sets relevant to extended portions of the earth surface, paving the way to the systematic generation of advanced DInSAR products to monitor ground displacements at a very wide spatial scale.

For the first time a physical model, constrained by monitoring data, is used to derive a quantitative estimate of the evolution in time of magmatic gases that enter a hydrothermal system of an active volcano. The site is Campi Flegrei (west of Naples, in Italy), a caldera that had a large ground inflation in 1982–1984 followed by 20 yr of subsidence. More recently the behavior of the system has changed: the magmatic component of fumaroles has increased, swarms of earthquakes are more frequent, and the ground has started a general uplifting trend, indicating that the hydrothermal system undergoes repeated injections of magmatic fluid. Physical simulations of the process show that total injected fluid masses are the same order of magnitude as those emitted during small to medium size volcanic eruptions, and their cumulative curve highlights a current period of increasing activity. Gas emission studies coupled with physical modeling can be extremely effective in predicting magmatic evolution and eruptive activity at volcanoes.

The Neapolitan volcanic area includes three active and high-risk volcanoes: Campi Flegrei caldera, Somma–Vesuvius, and Ischia island. The Campi Flegrei volcanic area is a typical example of a resurgent caldera, characterized by intense uplift periods followed by subsidence phases (bradyseism). After about 21 years of subsidence following the 1982–1984 unrest, a new inflation period started in 2005 and, with increasing rates over time, is ongoing. The overall uplift from 2005 to December 2019 is about 65 cm. This paper provides the history of the recent Campi Flegrei caldera unrest and an overview of the ground deformation patterns of the Somma–Vesuvius and Ischia volcanoes from continuous GPS observations. In the 2000–2019 time span, the GPS time series allowed the continuous and accurate tracking of ground and seafloor deformation of the whole volcanic area. With the aim of improving the research on volcano dynamics and hazard assessment, the full dataset of the GPS time series from the Neapolitan volcanic area from January 2000 to December 2019 is presented and made available to the scientific community.

From 2006 to spring 2013, Campi Flegrei (CF) caldera, Italy, was mostly uplifting at an increasing rate, particularly high from 2011. We show that the 2011–2013 accelerated uplift and 1980–2010 inflation and deflation phases can be explained by a two-source conceptual model similar to that proposed by Amoruso et al. (2014) (reference model). However, pressurization of the sole thin quasi-horizontal ∼4000 m deep source, responsible for large-scale 1980–2010 deformation, can explain the whole 2011–2013 deformation, while activity of the shallower Solfatara hydrothermal source, responsible for residual 1980–2010 deformation, appears constant. These results suggest a predominantly magmatic unrest in 2011–2013. Near-real-time comparison of observations and reference model predictions can provide additional information for short-term eruption forecasting at CF; a similar approach could be followed also in other volcanic environments.