Analysis and Classiﬁcation of Natural and Human-Induced Ground Deformations at Regional Scale (Campania, Italy) Detected by Satellite Synthetic-Aperture Radar Interferometry Archive Datasets

: The high levels of geo-hydrological, seismic, and volcanic hazards in the Campania region prompted full data collection from C-band satellites ERS-1 / 2, ENVISAT, and RADARSAT within regional (TELLUS) and national (PST-A) projects. The quantitative analysis, interpretation, and classiﬁcation of natural and human-induced slow-rate ground deformations across a span of two decades (1992–2010) was performed at regional scale (Campania, Italy) by using interferometric archive datasets, based on the Persistent Scatterer Interferometry approach. As radar satellite sensors have a side-looking view, the post-processing of the interferometric datasets allows for the evaluation of two spatial components (vertical and E-W horizontal ones) of ground deformation, while the N-S horizontal component cannot be detected. The ground deformation components have been analyzed across 89.5% of the Campania territory within a variety of environmental, topographical, and geological conditions. The main part (57%) of the regional territory was characterized during 1992–2010 by stable areas, where SAR signals do not have recorded signiﬁcant horizontal and vertical components of ground deformation with an average annual rate greater than + 1 mm / yr or lower than − 1 mm / yr. Within the deforming areas, the coastal plains are characterized by widespread and continuous strong subsidence signals due to sediment compaction locally enhanced by human activity, while the inner plain sectors show mainly scattered spots with locally high subsidence in correspondence of urban areas, sinkholes, and groundwater withdrawals. The volcanic sectors show interplaying horizontal and vertical trends due to volcano-tectonic processes, while in the hilly and mountain inner sectors the ground deformation is mainly controlled by large-scale tectonic activity and by local landslide activity. The groundwater-related deformation is the dominant cause of human-caused ground deformation. The results conﬁrm the importance of using Persistent Scatterer Interferometry data for a comprehensive understanding of rates and patterns of recent ground deformation at regional scale also within tectonically active areas as in Campania region.


Introduction
Slow rate deformations of the ground surface can be caused by both natural dynamic processes (landslide, subsidence, groundwater oscillations; glacial, volcanic, and tectonic activity) and human activity (extraction of ground water and hydrocarbons, injection of fluids, exploitation of geothermal fields, mining, tunneling, excavation). These deformations can be measured by means of several geodetic tools, such as topographic leveling, GPS, aerial photogrammetry, and LIDAR surveys, but the Triassic to Early Miocene carbonate platform limestone and pelagic basin calcareous-pelitic sequences form the geological units in the Apenninic mountain chain sectors. These units are strongly deformed and thrusted mainly eastward [26][27][28], and are uncomfortably covered by Middle Miocene to Pliocene thrust-top basin fillings, formed by siliciclastic sequences (mainly clays, sandstones, and conglomerates) [29].
The Campania region territory is a tectonically active sector, currently experiencing extension in hinterland and axial sectors, causing opening of fluvio-lacustrine intramontane basins, and along the Tyrrhenian margin forming the Campanian plain [30].
The Quaternary extensional tectonic processes are also testified by the presence of active faults and seismogenic sources along the Apennine chain, which are the main cause of the strong seismic activity (M ≤ 7.0) occurring in Campania (http://emidius.mi.ingv.it/CPTI15-DBMI15/). The seismogenic sources are mainly related to NW-SE trending normal faults. Hypocentral depths are typically of about 10-12 km and focal mechanisms show extensional kinematics. In the XX century, Triassic to Early Miocene carbonate platform limestone and pelagic basin calcareous-pelitic sequences form the geological units in the Apenninic mountain chain sectors. These units are strongly deformed and thrusted mainly eastward [26][27][28], and are uncomfortably covered by Middle Miocene to Pliocene thrust-top basin fillings, formed by siliciclastic sequences (mainly clays, sandstones, and conglomerates) [29].
The Campania region territory is a tectonically active sector, currently experiencing extension in hinterland and axial sectors, causing opening of fluvio-lacustrine intramontane basins, and along the Tyrrhenian margin forming the Campanian plain [30].
The Quaternary extensional tectonic processes are also testified by the presence of active faults and seismogenic sources along the Apennine chain, which are the main cause of the strong seismic The documented vertical displacements evidence a Pleistocene uplift rate of Region of 0.24 mm/yr in the southwestern coastal sector, while the chain sector uplifted during Pleistocene at about 0.51 mm/yr [30]. Available leveling lines and GPS data evidence subsidence at Vesuvius and Ischia [35,36] and show current southwestward velocities between 2 and 4 mm/yr at the southwestern corner of the Campanian plain [37]. At Campi Flegrei, a general subsidence between 1985 and 2007 has been interrupted by a still growing uplift phases [38].

Persistent Scatterer Interferometry (PSI) Datasets Available for Campania Territory
InSAR is a geodetic tool widely used since the 90s to measure and monitor changes in the ground shape by measuring the difference in the straight line between radar sensor and ground target (line of sight, LOS) between space-borne radar acquisitions taken at different times [2].
In the area of the Campania Region, Persistent Scatterer Interferometry techniques (PSI), based on the analysis of radar images acquired by the C-band (5.6 cm wavelength) sensors onboard ERS-1/2, ENVISAT, and RADARSAT satellites, are able to supply a multi-temporal comprehensive view on the spatial distribution of ground deformation during the early 1990s to late 2000s. As these SAR satellites fly in the near-polar orbit, the N-S surface displacements have little contributions to the PSI derived LOS measurements. A simplified geometry is assumed for the ascending and descending SAR acquisitions, where the N-S horizontal displacement component is not considered in the LOS measurements. The vertical and E-W horizontal displacement components can thus be resolved from the ascending and descending LOS measurements by assuming the same incidence angles for the ascending and descending orbits [8].
Among MT-InSAR techniques, the Persistent Scatterer Interferometry (PSI) approach [2][3][4][5][6]9,12,39,40] is based on the use of a long series of co-registered, multitemporal SAR imagery allowing for precise and robust results [6,41]. The PSI techniques rely on the analysis of pixels that remain coherent over a sequence of interferograms and allow for the production of maps of displacement velocity, which are measured along the LOS. These maps show the temporal evolution of the ground displacement (time series) at each SAR acquisition epoch. The main results for each coherent PS include the annual average velocity of the displacement calculated by linear fitting of the time series of the variation in the sensor-target distance along the LOS over the entire acquisition period. The relative measurements of deformation and the zero in the time series are respectively referred to as a stable reference point, assumed motionless, and a master scene [42].
For the study area, six PSI datasets (Table 1) were collected and post-processed: (a) ENVISAT ascending orbit, (b) ENVISAT descending orbit, (c) RADARSAT ascending orbit, (d) RADARSAT descending orbit, (e) ERS-1/2 ascending orbit, and (f) ERS-1/2 descending orbit ( The PSI datasets include about 5.8 million PSs, which were identified within the study area with a coherence higher than 0.65. The used datasets are characterized by negative mean values of displacement velocities in both ascending and descending orbits (Table 2). Specific PSI techniques adopted for each used dataset are also listed ( Table 2). ERS-1/2 datasets were processed both with PS-InSAR [19][20][21] and PSP [22,23] methods, while the ENVISAT dataset is available only with the PSP method [22,23] and RADARSAT dataset only with the PS-InSAR method [19][20][21]. Therefore, for obtaining a greater consistency, the datasets processed with the same method (i.e., PSP [22,23]) were used for both European Agency satellite datasets (ERS and ENVISAT), while the only available dataset was used for RADARSAT (PS-InSAR method [19][20][21]). Each point of the PSI datasets is identified by a set of attributes (Table 2), including (a) identifier code; (b) coordinates (North, East); (c) time series of measurements (expressed in mm) of the displacement along the LOS; (d) average velocity expressed in mm/yr derived by the linear regression of (c); (e) standard deviation of the average velocity; (f) coherence, that is a normalized index of the local signal-to-noise ratio of the interferometric phase and reflects the accuracy of PS measurements [40,[50][51][52].
C-band PS time series are robustly constrained when measuring deformation at rates greater than 0.5-1.0 mm/yr [41,51,53,54]. Costantini et al. [10] confirmed the capabilities of PSI to measure the displacement of corner reflectors with millimetric accuracy, but also on typical PSs, corresponding to man-made structures, rocks, bare soil, and, in general, nonvegetated and noncultivated terrain. PSI processing usually allows to obtain very accurate measurements within 1.0 mm/year for the PS average velocity along the LOS by assuming a threshold value greater than 0.65 for coherence [51,53,55,56] when a large number of SAR images is used for the processing, as for used datasets in this study (Table 1). Previous studies [5,14] have demonstrated that SAR observations over the Campania territory result in high quality, strongly coherent datasets over time intervals of several years (from 6 to 10).
The validation process of PSI database produced by the Not-Ordinary Plan of Environmental Remote Sensing (PST-A) have been discussed in Di Martire et al. [11], while that referring to the Tellus project have been discussed in Vilardo et al. [5,14,21].
The six PSI datasets were spatially processed and georeferenced to the projection WGS-84 UTM Zone 33N using ArcGIS 10.2 software, producing a LOS velocity point map (i.e., vLOS asc and vLOS desc maps). In Figure 2, the data distribution for the used datasets is shown, referring to the administrative boundaries of the region to be studied. Specifically, the LOS velocity shows varying trends for the datasets referring to different satellites and orbits. There are some sectors with no data in the areas with high topographic relief and/or unfavorable land cover (vegetation or agricultural crops). Most of the PSs are in the plain sectors, where the density of buildings and infrastructures is very high. The LOS displacement rates vary from −24.5 to +15.2 mm/yr for ERS datasets, −51.1 to +30.7 mm/yr for RADARSAT datasets and −10.6 to +6.0 mm/yr for ENVISAT datasets ( Figure 2). The ENVISAT descending dataset is lacking for the southern sector of the regional territory.

PSI Datasets Post-Processing for Ground Deformation Components
Post-processing activities ( Figure 3) were finalized to the assessment and mapping of ground deformation processes that occurred during 1992-2010 in the Campania territory.
As SAR satellite sensors have a side-looking view, PSI displacements can be measured only along the LOS direction, which in the study area is in the range of 22-34 • (Table 1). Observed increases or decreases in the LOS can be due to ground subsidence and uplift, respectively, although horizontal motions do occur and can dominate in some locations [24]. For those areas common to both acquisition geometries, the overlapping availability of ascending and descending PSI datasets allows for the evaluation of two spatial components (vertical and E-W horizontal ones) of ground deformation, based on simple trigonometric rules [5,14], as the N-S horizontal component is not detected by the SAR satellite acquisition system.
Firstly, the regional scale annual average LOS velocity raster maps (vLOS asc -R and vLOS desc -R maps), referring to ascending and descending orbits of ERS-1/2, RADARSAT and ENVISAT PSI datasets, were produced by applying a spatial interpolation of the average velocity PS's point values with the Inverse Distance Interpolation Weighted (IDW) method. The IDW approach is commonly used to interpolate datasets formed by scattered points, thus allowing both reliable results and preservation of local variability of the data [17,18,[57][58][59]. In this study case, an interpolation method consisting of a quadratic weighting power of 2 within a 1000-m radius neighborhood was used to obtain 100-m regularly spaced grids, namely the ascending and descending LOS velocities raster maps derived from ERS-1/2, RADARSAT, and ENVISAT datasets with coherence > 0.65.   The components of the velocity were calculated by combining the vLOSasc and vLOSdesc raster maps on pixels common to both maps [5,12,14,35,49,60], by assuming that the ascending and descending LOS belong to the same East-Z plane and the look-angle is the same for both ascending and descending geometries. Based on these assumptions, the following equations [5,14] were used: vHori vLOSdesc vLOSasc /2 sin , where "v" is the displacement velocity vector of an investigated PS; vVert and vHori are, respectively, the v projection along the Cartesian vertical and horizontal axes; vLOSdesc and vLOSasc are the projections of velocity along the LOS of the two orbits; and q is the look-angle. A value of q given by the average of the two LOS incidence angles of the ascending and descending orbits was used for the calculations related to each satellite dataset, which was 22.5° for ERS-1/2, 33.25° for RADARSAT, and 23.5° for ENVISAT.
The obtained 100-m spaced grid maps ( Figure 4) show the distribution of the vertical and E-W horizontal components of ground deformation velocity in the Campania region for each dataset The components of the velocity were calculated by combining the vLOS asc and vLOS desc raster maps on pixels common to both maps [5,12,14,35,49,60], by assuming that the ascending and descending LOS belong to the same East-Z plane and the look-angle is the same for both ascending and descending geometries. Based on these assumptions, the following equations [5,14] were used: where "v" is the displacement velocity vector of an investigated PS; vVert and vHori are, respectively, the v projection along the Cartesian vertical and horizontal axes; vLOS desc and vLOS asc are the projections of velocity along the LOS of the two orbits; and q is the look-angle. A value of q given by the average of the two LOS incidence angles of the ascending and descending orbits was used for the calculations related to each satellite dataset, which was 22. satellite dataset is referred (i.e., 10 years for ERS-1/2, 2.5 years for RADARSAT, and 6.5 years for ENVISAT, as the overlapping years of RADARSAT and ENVISAT datasets count 0.5 year each). With reference to the southern sector of the regional territory, as the ENVISAT descending dataset is lacking, the average velocities were calculated by accounting only for ERS and RADARSAT grid maps.The use of SAR interferometry measurements to analyze ground deformation components has already been tested and validated in several sectors of the Campania Region (i.e., Campi Flegrei, Volturno plain, Sele plain) with millimetric levels of precision [5,13,14].  In vertical GD maps, red and blue colors, respectively, represent subsidence and uplift, while in horizontal GD maps, red and blue colors, respectively, represent westward and eastward movements.

PSI Datasets Post Processing for Ground Deformation Classification
After having classified vVertT and vHoriT raster maps (see Table 3 for used codes), the raster datasets were converted to two polygon features (VTF and HTF, respectively). Then a geometric union of the input features was computed for obtaining an output feature class assigning to each In vertical GD maps, red and blue colors, respectively, represent subsidence and uplift, while in horizontal GD maps, red and blue colors, respectively, represent westward and eastward movements.
In order to obtain a quantitative assessment of the two components of the ground deformation velocities (expressed in mm/yr) referring to the whole analyzed period (evaluated 19 years, as 1992 and 2010 partial acquisitions were considered representative of the full year), the average velocities were calculated with a pixel by pixel weighted average of the velocity of the different grid maps ( Figure 3). The weighted average was based on the number of years of the time interval to which each satellite dataset is referred (i.e., 10 years for ERS-1/2, 2.5 years for RADARSAT, and 6.5 years for ENVISAT, as the overlapping years of RADARSAT and ENVISAT datasets count 0.5 year each). With reference to the southern sector of the regional territory, as the ENVISAT descending dataset is lacking, the average velocities were calculated by accounting only for ERS and RADARSAT grid maps.
The results are two 100-m spaced grid maps showing, respectively, the cumulated distribution of vertical (vVert T ) and E-W horizontal (vHori T ) movements in the Campania territory (expressed in mm/yr) for the whole period of 1992-2010.
The use of SAR interferometry measurements to analyze ground deformation components has already been tested and validated in several sectors of the Campania Region (i.e., Campi Flegrei, Volturno plain, Sele plain) with millimetric levels of precision [5,13,14].

PSI Datasets Post Processing for Ground Deformation Classification
After having classified vVert T and vHori T raster maps (see Table 3 for used codes), the raster datasets were converted to two polygon features (VTF and HTF, respectively). Then a geometric union of the input features was computed for obtaining an output feature class assigning to each polygon feature the attributes related to both vertical and horizontal components of ground deformation within the nine possible combinations, i.e., U-W, U-N, U-E, S-W, S-N, S-E, D-W, D-N, and D-E classes. The class S-N means that no horizontal and vertical ground deformation have been recorded during 1992-2010, and so it designates stable areas. The remaining eight classes were used for areas affected by at least one component of ground deformation with an average annual rate greater than +1 mm/year or lower than −1 mm/year. The final Ground Deformation Classification map (GDC map) includes only features extending more than 0.1 km 2 in order to exclude too local deformations.

Ground Deformation Components
The interplay among different deformation processes (i.e., tectonic, volcanic, hydrothermal, hydrogeological, geomorphologic, and anthropic) causes a complex pattern of horizontal and vertical displacements during the 1992-2010 time range (Figures 5 and 6).
The PSI post-processing for the analysis of ground deformation components resulted in the assessment of both vertical and horizontal components for an area of 12,160 km 2 , representing the 89.5% of the whole regional territory (Table 4). No data areas (10.5%) are mainly located in steep mountainous inner sectors due to the dark gap generated by layover and shadow effects during SAR signals' acquisition. The almost total coverage of the regional territory allows having an overall and effective view of the effects of the various slow deformation processes acting on the ground at different scales. Table 4. Area distribution for vertical and E-W horizontal velocity value ranges. The main part (71%) of the classified regional territory results are stable with reference to the vertical component of ground deformation ( Figure 5, Table 4). Low subsidence rates (−3 to −1 mm/yr) are widespread (ca. 14.5%) in mountain and hilly Sannio and Irpinia sectors, in several coastal areas, and in southern Cilento; also, Roccamonfina and Vesuvius areas show subsidence trends. The areas affected by significant subsidence are presenting a minor extension (ca. 2.7% of the classified territory) and are characterized by high negative values (−23 to −3 mm/yr). These zones include ( Figure 5) the Volturno river plain, the Sele river coastal plain, Ischia Island, Campi Flegrei central sectors, and other sparse spots. The two sectors characterized by moderate uplift, testified by positive values (+1 to +5 mm/yr), are located in the central sector of the Campanian plain and in the eastern sector of the Cilento (i.e., Alburni Massif).

Velocity Range (mm/yr) Area (km 2 ) with Vertical Velocity Area (km 2 ) with Horizontal Velocity
Even if about half (52%) of the classified regional territory results are stable with reference to the E-W horizontal component of ground deformation ( Figure 6, Table 4), the regional territory is clearly divided into two main sectors. The western sector, encompassing the Campanian plain, the Roccamonfina relief, and the Sorrento peninsula, shows a slight eastward movement

Ground Deformation Classification
The PSI post-processing for the ground deformation classification allowed the classification for almost all (96.8%) of the regional territory ( Table 5). The S-N class characterizes the main part (ca. 57%) of the regional territory and designates stable areas during 1992-2010, where SAR signals do not have recorded significant horizontal and vertical component of ground deformation (Figure 7). The other eight classes characterize areas affected by at least one component of ground deformation with an average annual rate greater than +1 mm/yr or lower than −1 mm/yr.
The most common class (19.7%) is the W-S where only a westward horizontal component was recognized. This deformation type affects large inner hilly and mountainous sectors of Irpinia and Cilento (Figure 8a). The other main common classes are the N-D class (7.5%), where only a subsidence signal was recognized (Figure 8b), and the E-S class (6.7%), characterized by an eastward horizontal component (Figure 8a) mainly affecting coastal and intra-montane plains.

Ground Deformation Classification
The PSI post-processing for the ground deformation classification allowed the classification for almost all (96.8%) of the regional territory ( Table 5). The S-N class characterizes the main part (ca. 57%) of the regional territory and designates stable areas during 1992-2010, where SAR signals do not have recorded significant horizontal and vertical component of ground deformation (Figure 7). The other eight classes characterize areas affected by at least one component of ground deformation with an average annual rate greater than +1 mm/yr or lower than −1 mm/yr.
The most common class (19.7%) is the W-S where only a westward horizontal component was recognized. This deformation type affects large inner hilly and mountainous sectors of Irpinia and Cilento (Figure 8a). The other main common classes are the N-D class (7.5%), where only a subsidence signal was recognized (Figure 8b), and the E-S class (6.7%), characterized by an eastward horizontal component (Figure 8a) mainly affecting coastal and intra-montane plains.
The final ground deformation classification maps ( Figure 7) include only features extending more than 0.1 km 2 in order to exclude too local deformations. If the areal extension of the mapped classified areas is considered (Table 6), the main part of them falls within the local ranges 0.1-1 km 2 and 1-10 km 2 .
The largest classified area (>1000 km 2 range) and some others in the intermediate ranges (10-100 km 2 and 100-1000 km 2 ) are referred to as stable sectors (class N-S). The final ground deformation classification maps (Figure 7) include only features extending more than 0.1 km 2 in order to exclude too local deformations. If the areal extension of the mapped classified areas is considered (Table 6), the main part of them falls within the local ranges 0.1-1 km 2 and 1-10 km 2 . The largest classified area (>1000 km 2 range) and some others in the intermediate ranges (10-100 km 2 and 100-1000 km 2 ) are referred to as stable sectors (class N-S).

Analysis of Regional to Local Ground Deformation
In the following, the large scale ground deformation trends occurring in the different sectors (coastal plain, inner plain, Apennine sector, volcanic areas) of the Campania territory were analyzed in light of the available tectonic, volcanological, hydrogeological, geomorphological, and land use data. In addition, some examples of local scale ground deformations are described to show the potentiality of the used approach.

Coastal Plain Sectors
The Campania coastal plain sectors fall within the passive continental margin of Tyrrhenian Sea affected by a long-term tectonic subsidence process characterized by large-scale downward movements [27,28,61] with an average rate of about -0.2 mm/yr [62]. The analysis of the vertical component of ground deformation within the main coastal plains (Garigliano, Volturno and Sele rivers alluvial coastal plains) highlights the presence of wide low coast sectors characterized by strong subsidence movements with −2.5 mm/yr average rates, while Sarno and Alento rivers alluvial coastal plains show lower average rates (ca. −0.5 mm/yr).
The subsidence of alluvial coastal plain can be considered as a natural process mainly due to the compaction of the alluvial sediments infill under the lithostatic load, and the anthropic influences due to water pumping and urbanization loads are only additional factors that locally enhance the subsidence effects [17,18]. In Volturno, Sarno, and Sele river plains urban areas the subsidence-induced ground displacements due to groundwater withdrawals have relevant consequences to the exposed facilities [63,64]. The coastal subsidence phenomena are contributing to increasing the exposure of the subsiding low coastal areas to hazards related to inundation and erosion processes, and may cause a further increase in future relative sea level rates at local scale [17,65].
The vertical displacement values obtained for the Volturno river plain show a significant subsidence in the central axial sectors (−2.5 to −22 mm/yr) and in the river mouth area (−2.5 to −8 mm/yr). Moderate uplift is detected in the eastern part of the plain (+0.5 to +2.0 mm/yr), whereas other sectors of the study area are characterized by moderate subsidence and/or stability. A general eastward trend is also recognized (Figure 9). The subsidence recorded in the Volturno plain is mainly a consequence of a natural process related to the compaction under the lithostatic load of the fluvial and palustrine deposits that form the alluvial plain. The magnitude of the recorded subsidence has been found to be greater when thick peat layers occur in the subsoil [18]. The anthropic influences (e.g., water exploitation and urbanization) is able only to locally increase the subsidence effect [18], while in some sectors the negative balance between the water recharge rates of the Volturno hydrographic basin and the drainage operated by the artificial channeling, related to industrial and intense agricultural activities [5,66], may have larger effects on subsidence rates.
The uplift recognized in the eastern sector of the plain is related to the interplay between tectonic activity and hydrogeological conditions (Figure 9). The uplift can be due to the action of an E-W oriented active fault that bounds the area to the south [18]. Besides, the study area is characterized by a porous multi-layered aquifer system, formed by volcaniclastic and alluvial deposits, and the groundwater circulation can be considered unitary in the aquifer system. If the local hydrogeological conditions [67,68] are considered, the area is located on a minor groundwater divide (i.e., Caivano-Campi Flegrei, Figure 9) that is recharged by lateral groundwater flux coming from adjoining carbonatic aquifer sectors made by the calcareous massif bordering to the east the Campanian plain. In these conditions, the water fluxes may locally increase the soil pore pressure causing local soil oversaturation, which could be able to partly explain the local uplift observed in the Marcianise-Caivano sector.
The ground deformation components and classification ( Figure 10) show that the coastal sector of the Sele Plain is characterized by a general westward horizontal deformation and a complex vertical pattern, validated by GPS surveys [17]. The northern and central sectors of the plain are characterized by relative subsidence rates of about −3 to −7 mm/yr along the coastline and by stability in the hilly inland area; around the Sele river mouth, a narrow area with subsidence rates up to −8 mm/yr is present, while the southern sector is characterized by general condition of stability with minor subsiding areas.
The vertical ground deformation pattern is correlated with the non-uniform stratigraphy of the coastal Quaternary infill because of the presence of layers of clastic sediments with different thickness and degrees of compaction in the two sectors. In detail, the subsidence is higher in the northern sector of the plain, where the Quaternary alluvial-coastal deposits are thicker because of the structural asymmetry of the graben, while in the southernmost sector the rates are the lowest recorded in the plain due to the minor thickness of Holocene deposits and to the presence of very thick sedimentary bodies of travertine [17]. These differences in ground deformation rates expose the northern and central sectors of the plain to higher levels of inundation and erosion hazard than in the southern sector [17]. The contribution of the local tectonic deformation pattern could suggest prevailing dip-slip movements along the NNW-SSE striking faults and the E-W trending faults bounding the Sele coastal plain [5].
Sarno and Alento river plains show lower rates of subsidence in the coastal sectors ( Figure 5) [65]. The Sarno river coastal plain is characterized by moderate subsidence (up to −5 mm/yr) only near the river mouth, while the Alento river coastal plain shows low subsidence rates (up to −2 mm/yr) in the northern sector, while the coastal narrow strip displays stability. Along the Alento River course, a hot spot of subsidence (up to −3 mm/yr), developing 1-2 km inland, is more evident.

Volcanic Sectors
Ground deformation in active volcanic areas of Vesuvius, Campi Flegrei and Ischia has been extensively studied with PS-InSAR techniques (in [5,14] and references therein), while no ground deformation analysis has been done about the not active Roccamonfina volcanic system.

Volcanic Sectors
Ground deformation in active volcanic areas of Vesuvius, Campi Flegrei and Ischia has been extensively studied with PS-InSAR techniques (in [5,14] and references therein), while no ground deformation analysis has been done about the not active Roccamonfina volcanic system.
The ground deformation patterns referring to the whole 1992-2010 period ( Figure 11) show that the eastern sector of Campi Flegrei is characterized by westward velocity whereas the western sector is characterized by eastward velocity with maximum values in the caldera central area (around Pozzuoli), where also subsidence increases. This is only a cumulated deformation related to the different contribution experienced during the 1992-2010 period. In fact, at Campi Flegrei, after a strong uplift (about +3.  [38]). This ground deformation pattern is representative of the whole caldera area, even if with lower magnitudes in peripheral sectors. The largest ground deformation is localized within and around the structural border of the Campi Flegrei caldera and a systematic recurrence of opposite trends (uplift vs. subsidence) in the ground deformation of the inner caldera region with respect to the surrounding areas has been recognized [16]. The analysis of PSI velocity and acceleration annual variations also revealed intense yearly dynamics of the Campi Flegrei caldera collapse-resurgence system. This ground deformation field, combined with the re-activation of the caldera ring-faults, intra-caldera faults, and eruptive fissures, indicates a contraction (deflation) of the caldera due to a depressurization of the hydrothermal system and degassing from a magmatic reservoir (in [14] and references therein). The rates of vertical and horizontal components of ground deformation ( Figure 11) are in full agreement with deformation data based on levelling and GPS surveys made in the last decades [14,38].
Ischia Island shows a general vertical subsidence with velocities between −1 and −10 mm/yr, centered in the inner sector of the Mt. Epomeo volcanic complex, coupled with a general eastward horizontal displacement field characterized by a different trend in the eastern sector of the island (higher eastward velocity, from +3 to +25 mm/yr) respect to the central-western one (up to +3 mm/yr), suggesting an overall contraction movement ( Figure 11).
The vertical component velocities are comparable with those obtained by levelling lines (for example −12.7 mm/year at benchmark BM100; [36]) referred to the same time period [5]. The E-W striking, south dipping normal fault mechanism of the 21 August 2017 earthquake (Mw 3.9-Md 4.0; Imax EMS 8; hypocenter depth 800 m) that struck the northern sector of Ischia Island is very likely induced by the observed long-term subsidence phase, since the lithostatic load represents the principal vertical stress [69].
Ischia ground deformation pattern ( Figure 11) cannot be explained by a typical volcanic source [70,71]. The Ischia deformation pattern is related to a combination of endogenous and exogenous processes that include deflation of the island related to the depressurization of the local hydrothermal system, fault activity and landslides due to gravity instability on slopes [5,35,36], and to the coupling effects of crust rheology and the gravitational loading of the volcano [70]. The Somma-Vesuvius is an asymmetric, polygenic volcanic complex formed by the superimposition of two edifices, as the older Mt. Somma with a summit caldera and the younger Vesuvius cone. Its morphostructure results from the combined action of NW-SE faulting and large caldera collapses that occurred in 18 and 79 AD [72,73].
The ground deformation components ( Figure 12) evidence two separated zones of continuous subsidence within the Somma-Vesuvius volcano edifice. The first zone is the Vesuvius central cone and the southern flank, while the second one is represented by a discontinuous strip extending around the volcanic edifice at about 10 km of distance from the crater, where the outer flanks of the The Somma-Vesuvius is an asymmetric, polygenic volcanic complex formed by the superimposition of two edifices, as the older Mt. Somma with a summit caldera and the younger Vesuvius cone. Its morphostructure results from the combined action of NW-SE faulting and large caldera collapses that occurred in 18 and 79 AD [72,73].
The ground deformation components ( Figure 12) evidence two separated zones of continuous subsidence within the Somma-Vesuvius volcano edifice. The first zone is the Vesuvius central cone and the southern flank, while the second one is represented by a discontinuous strip extending around the volcanic edifice at about 10 km of distance from the crater, where the outer flanks of the volcanic edifice lay on the alluvial-marine sequences of the Campanian plain. The cumulated subsidence rates are mainly comprised between −1 and −10 mm/yr and are confirmed by the available GPS (up to −11.7 mm/yr in 2001-2012 [74]) and levelling data (up to −15.1 mm/yr in 1973-2009 [75]) for the flanks and crater. Somma-Vesuvius is affected by a small contraction phase, more marked in the areas with the greatest altitudes, according to the diffuse and modest subsidence observed in the central cone area [74].
This subsidence pattern is in agreement with the hypothesis that local seismicity, all clustered below the cone itself, is mainly driven by gravitational stress that should produce progressive subsidence of the highest relief zone [76][77][78]. The larger strip of annular, discontinuous subsidence is interpretable in terms of ring-like, shallow normal fault-like movements, occurring at the contact between the volcanic edifice and the rock basement, due to the high gravitational loading coupled with the embedding extensional tectonic stress field. At the main urban centers (Pomigliano D'Arco, Marigliano, and Saviano towns), located in the agricultural and industrial districts north of Somma-Vesuvius, several quasi-circular areas with subsidence rates up to −3 mm/yr can be observed, caused by the effects of water table changes [49,78]. In fact, the observed subsidence could reflect a decrease of pore pressure in the soil related to the intense artificial drainage from wells, causing a water deficit in the local water table [5].
The horizontal deformation ( Figure 12) shows a marked difference between the northern and eastern flanks of Mt. Somma, where reaches values up to +3 mm/yr (eastward movement), and the area encompassing the central cone sector and the southern and western flanks characterized by values up to −5 to −10 mm/yr (westward movement).
The available GPS data (up to 5 mm/yr in 2001-2012 [74]) confirm this pattern of horizontal deformation velocity. Vilardo et al. [5] considered the horizontal velocity field consistent with divergent movements between the western and eastern sectors of Vesuvius and Campanian plain due to a NNW-SSE tectonic structure crossing the volcano. Our data suggest that the horizontal pattern is more compatible with a gravitational deformation due to lateral SW-and W-directed collapses of Vesuvius Volcano driven by inherited tectonic faults (i.e., NW-SE normal faults and E-W strike-slip fault) as already hypothesized by Milia et al. [72].
In the sector located between Vesuvius and Campi Flegrei, our data suggest the presence of scattered ground deformation areas. In the north-eastern Naples urban area (Sebeto plain), a subsidence trend with rates of −1 to −5 mm/yr is present ( Figure 12).
This ground deformation trend can be caused by piezometric level lowering following groundwater pumping in the multi-layered pyroclastic-alluvial aquifer linked to a poro-elastic mechanism in the aquifer system [68]. Several small areas of subsidence (i.e., Vomero-Arenella and Scudillo-Stella districts) with rates up to −10 mm/yr induced by anthropic and natural processes are affecting residential districts of Naples urban areas.
A multiple association of triggering factors cause these subsidence processes [13] that are both of anthropic (subsoil excavations for the construction of transport infrastructures, filling/emptying cycles of large underground water reservoirs, rise of the water table due to the stop of ground water withdrawal) and natural (gravity slope instability related to local morphological factors, re-activation of the bradyseism phases in the Campi Flegrei caldera) origins.
The Roccamonfina volcanic complex is an extinct stratovolcano that was active from 550 to 150 ka in the Garigliano river rift valley. It and was affected by an intense plinian activity revealed by very large craters.
The central caldera is the result of the eruptive explosions at 353 ± 5 ka, while the latest stage of activity featured the edification of the central shoshonitic domes at 150 ka [79,80]. The ground deformation components ( Figure 13) evidence a large area of subsidence with vertical rates between −1 and −3 mm/yr centered on the eastern sector of the caldera rim; the E-W horizontal deformation is less relevant and shows discontinuous positive rates on the NE flanks of the volcanic complex and a localized westward deformation within the caldera southern sector (Figure 13).
Even if no interpretation about the residual volcano-tectonic activity geodynamics of Roccamonfina sector is available in literature, the generalized subsidence could be related to the structures produced by the caldera collapses, while the eastward horizontal movements on the eastern flanks could be also linked to slope instabilities.   [67,68], IFFI project landslides, and active fault datasets [79,80].

Inner Plain Sectors
The lowland area located in the inner Sarno plain is characterized by variable subsidence velocity patterns with values up to −6 mm/yr, coupled with eastward horizontal movement up to +3 mm/yr ( Figure 12). Several artificial wells for the water extraction occur in this area and are mainly responsible for the observed subsidence [81]. A strong correlation between the amounts of pumped groundwater and the groundwater level in the aquifer and, in turn, between this and the subsidence rates has been demonstrated in a study on ancient masonry buildings located in the urban area of Sarno. The available levelling data show that vertical negative movements up to −5 mm/yr (subsidence) affect the area confirming DInSAR results [63,82]. In addition, the historical and recent ground failures occurring in the densely populated urban center of Sarno can be ascribed either to the presence of underground cavities or to piping sinkholes [83].

Inner Plain Sectors
The lowland area located in the inner Sarno plain is characterized by variable subsidence velocity patterns with values up to −6 mm/yr, coupled with eastward horizontal movement up to +3 mm/yr (Figure 12). Several artificial wells for the water extraction occur in this area and are mainly responsible for the observed subsidence [81]. A strong correlation between the amounts of pumped groundwater and the groundwater level in the aquifer and, in turn, between this and the subsidence rates has been demonstrated in a study on ancient masonry buildings located in the urban area of Sarno. The available levelling data show that vertical negative movements up to −5 mm/yr (subsidence) affect the area confirming DInSAR results [63,82]. In addition, the historical and recent ground failures occurring in the densely populated urban center of Sarno can be ascribed either to the presence of underground cavities or to piping sinkholes [83].
Ground failures ascribed to piping sinkholes involved also the urban centers of Telese Terme (Figure 14), where collapses occurred in February 2002 and in August 2006 in a suburban area recently urbanized [21,83]. Our data show that signifi

Apennine Sectors
The Southern Apennines sector of the Campania Region appears to be a relatively stable area with reference to vertical deformation (Figure 8b) with scattered and restricted spots of negative vertical rates within −1 to −5 mm/yr. Only the large sector near Alburni massif shows positive vertical rates. Instead, large hilly and mountainous sectors in Irpinia, Sannio and Cilento are characterized by the westward horizontal component of ground deformation (Figure 8a).
The small scattered subsidence areas are probably related to landslide and sinkhole processes along the hilly and mountainous slopes. Most of the mapped landslides overlap the areas with subsidence and westward horizontal deformation (Table 7, Figure 15). Several studies used the PSI data for upgrading landslide maps and inventories in hilly and mountain sectors and defining their state of activity, especially for town-damaging landslides [12,21,84]. PSI data were also used to classify existing landforms susceptible to slow landslides along the Tammaro river valley in the central sector of Sannio area [85,86]. The detailed analysis of several case studies about towns affected by landslides in a mountain site, large landslides along fluvial valleys, and Deep Seated Gravitational Slope Deformation (DSGSD), confirm the results obtained in this study also at local scale.

Apennine Sectors
The Southern Apennines sector of the Campania Region appears to be a relatively stable area with reference to vertical deformation (Figure 8b) with scattered and restricted spots of negative vertical rates within −1 to −5 mm/yr. Only the large sector near Alburni massif shows positive vertical rates. Instead, large hilly and mountainous sectors in Irpinia, Sannio and Cilento are characterized by the westward horizontal component of ground deformation (Figure 8a).
The small scattered subsidence areas are probably related to landslide and sinkhole processes along the hilly and mountainous slopes. Most of the mapped landslides overlap the areas with subsidence and westward horizontal deformation (Table 7, Figure 15). Several studies used the PSI data for upgrading landslide maps and inventories in hilly and mountain sectors and defining their state of activity, especially for town-damaging landslides [12,21,84]. PSI data were also used to classify existing landforms susceptible to slow landslides along the Tammaro river valley in the central sector of Sannio area [85,86]. The detailed analysis of several case studies about towns affected by landslides in a mountain site, large landslides along fluvial valleys, and Deep Seated Gravitational Slope Deformation (DSGSD), confirm the results obtained in this study also at local scale.  Figure 15. Comparison of ground deformation classification with landslide regional inventory (IFFI project).   In the following sections, some relevant examples are described ( Figure 16). Castelpagano is a small city located in the high valley of the Fortore river, characterized by severe ground deformation affecting the urban area (Figure 16a) due to different types of landslide (earthflow, rotational slide, translational slide and complex). The ground deformation components show up to −18 mm/yr (westward) and up to −4 mm/yr (subsidence) indicating a general landslide direction toward SW compatible with the local geomorphological conditions. Field monitoring with GPS surveys and clinometer measurements in boreholes are coherent with the results of this study [15]. A large landslide affected the town of Calitri, located in eastern Irpinia in the Ofanto river valley, in the 80s to early 90s. The Calitri landslide is a complex slope movement, composed by a large deep-seated slide, reactivated after the 1980 earthquake, and a limited shallow mudslide showing impulsive activity linked with rainfall events [87]. No significant movements affected the area (Figure 16b) during the analyzed period (1992-2010), as confirmed by field monitoring activities [87]. In the following sections, some relevant examples are described ( Figure 16). Castelpagano is a small city located in the high valley of the Fortore river, characterized by severe ground deformation affecting the urban area (Figure 16a) due to different types of landslide (earthflow, rotational slide, translational slide and complex). The ground deformation components show up to −18 mm/yr (westward) and up to −4 mm/yr (subsidence) indicating a general landslide direction toward SW compatible with the local geomorphological conditions. Field monitoring with GPS surveys and clinometer measurements in boreholes are coherent with the results of this study [15]. A large landslide affected the town of Calitri, located in eastern Irpinia in the Ofanto river valley, in the 80s to early 90s. The Calitri landslide is a complex slope movement, composed by a large deep-seated slide, reactivated after the 1980 earthquake, and a limited shallow mudslide showing impulsive activity linked with rainfall events [87]. No significant movements affected the area (Figure 16b) during the analyzed period (1992-2010), as confirmed by field monitoring activities [87]. Moio della Civitella is a small town located in southern Cilento on a slope formed by mainly argillaceous-calcareous rocks. The deformation components along the slope (Figure 16c) are comprised up to −16 mm/yr (westward) and up to −2.5 mm/yr (subsidence) confirming the landslide direction toward SW. The ground deformation effects are related to large slow to very slow landslides, confirmed by GPS surveys, and are coupled with very high levels of damage in masonry old buildings. Actually, the damages are caused by other causes such as thermic effects, bad plans, Moio della Civitella is a small town located in southern Cilento on a slope formed by mainly argillaceous-calcareous rocks. The deformation components along the slope (Figure 16c) are comprised up to −16 mm/yr (westward) and up to −2.5 mm/yr (subsidence) confirming the landslide direction toward SW. The ground deformation effects are related to large slow to very slow landslides, confirmed by GPS surveys, and are coupled with very high levels of damage in masonry old buildings. Actually, the damages are caused by other causes such as thermic effects, bad plans, and inadequate foundations [12,21,88,89] and enhance those induced by slope movements, resulting in very high values of SAR ground deformation signals. Finally, a DSGSD affects Bisaccia, a small town located in eastern Irpinia in the Cervaro river high valley. The DSGSD involves a brittle lithotype (Pliocene conglomerates) resting over a structurally complex, mainly pelitic unit ("Argille Varicolori"). The deformation components measured in this study are comprised of up +8 mm/yr (eastward) and up to −13 mm/yr (subsidence) confirming the slope deformation kinematism. As a consequence of repeated seismic actions (M 6.7 in 1930, M 6.1 in 1962, M 6.9 in 1980 earthquakes), the top rigid plate made up of lithified conglomeratic layers resulted in being split in five portions, showing different rates of vertical and horizontal displacements (within some cm/yr), causing severe damage to the old village settlement [90].
The dams providing lakes forming water reservoirs for agricultural and town needs are another source of ground deformation that is present in the inner sectors. For example, the Conza della Campania earth dam (Figure 17), located in the Ofanto River valley in Pliocene clay, show relevant deformations. The deformation components measured in this study are comprised of up +4 mm/yr (eastward) and up to −2 mm/yr (subsidence) confirming the earth dam ground deformation pattern obtained by in-situ monitoring conventional instrumental data (levelling, extensometers, etc.) [91]. and inadequate foundations [12,21,88,89] and enhance those induced by slope movements, resulting in very high values of SAR ground deformation signals. Finally, a DSGSD affects Bisaccia, a small town located in eastern Irpinia in the Cervaro river high valley. The DSGSD involves a brittle lithotype (Pliocene conglomerates) resting over a structurally complex, mainly pelitic unit ("Argille Varicolori"). The deformation components measured in this study are comprised of up +8 mm/yr (eastward) and up to −13 mm/yr (subsidence) confirming the slope deformation kinematism. As a consequence of repeated seismic actions (M 6.7 in 1930, M 6.1 in 1962, M 6.9 in 1980 earthquakes), the top rigid plate made up of lithified conglomeratic layers resulted in being split in five portions, showing different rates of vertical and horizontal displacements (within some cm/yr), causing severe damage to the old village settlement [90].
The dams providing lakes forming water reservoirs for agricultural and town needs are another source of ground deformation that is present in the inner sectors. For example, the Conza della Campania earth dam (Figure 17), located in the Ofanto River valley in Pliocene clay, show relevant deformations. The deformation components measured in this study are comprised of up +4 mm/yr (eastward) and up to −2 mm/yr (subsidence) confirming the earth dam ground deformation pattern obtained by in-situ monitoring conventional instrumental data (levelling, extensometers, etc.) [91]. A large sector located south to the Alburni massif ( Figure 18) shows significant positive vertical rates (uplift up to +5 mm/yr) combined with negative horizontal rates (westward horizontal motion up to −10 mm/yr), while the Vallo di Diano valley is characterized by subsidence. The observed ground deformation pattern can be caused by tectonic activity. The Alburni Massif and and Vallo di Diano are crossed by NW-SE striking normal active faults, which could act with a strike-slip component of motion. In fact, the major earthquakes (M > 5) univocally indicate active extension, while the minor seismicity shows a more complex pattern with strike-slip, transtensive and occasional compressional events [92].
Very large sectors in Irpinia, Sannio, and Cilento are characterized by a prevailing westward horizontal component of ground deformation with rates within −1 to −5 mm/yr (Figure 7). Smaller areas with eastward velocities of +1 to +3 mm/yr or westward velocities of −5 to −10 mm/yr are also scattered within these sectors, and are likely linked to local movements due to landslides ( Figure 15). In the northern part of Campania, the combined vertical and horizontal velocity patterns at the boundary between the Apennines sectors of Sannio and Irpinia and the Campanian plain clearly show that the western sector is affected by subsidence or null vertical deformation and eastward deformations, whereas the eastern sector is characterized by subsidence and westward movements A large sector located south to the Alburni massif ( Figure 18) shows significant positive vertical rates (uplift up to +5 mm/yr) combined with negative horizontal rates (westward horizontal motion up to −10 mm/yr), while the Vallo di Diano valley is characterized by subsidence. The observed ground deformation pattern can be caused by tectonic activity. The Alburni Massif and and Vallo di Diano are crossed by NW-SE striking normal active faults, which could act with a strike-slip component of motion. In fact, the major earthquakes (M > 5) univocally indicate active extension, while the minor seismicity shows a more complex pattern with strike-slip, transtensive and occasional compressional events [92].
Very large sectors in Irpinia, Sannio, and Cilento are characterized by a prevailing westward horizontal component of ground deformation with rates within −1 to −5 mm/yr (Figure 7). Smaller areas with eastward velocities of +1 to +3 mm/yr or westward velocities of −5 to −10 mm/yr are also scattered within these sectors, and are likely linked to local movements due to landslides ( Figure 15). In the northern part of Campania, the combined vertical and horizontal velocity patterns at the boundary between the Apennines sectors of Sannio and Irpinia and the Campanian plain clearly show that the western sector is affected by subsidence or null vertical deformation and eastward deformations, whereas the eastern sector is characterized by subsidence and westward movements ( Figure 8). Additionally, in the southern part of Campania (Cilento sector) the ground deformation pattern is mainly characterized by subsidence and westward movements (Figure 8).  Figure 8). Additionally, in the southern part of Campania (Cilento sector) the ground deformation pattern is mainly characterized by subsidence and westward movements (Figure 8). The spatial discontinuity between these two horizontal velocity fields (eastward at the west side and westward at the east and south sides) roughly overlaps the Campanian plain-Apennines morphological boundary (Figure 19), which is marked by NW-SE to NNW-SSE striking faults, moving with prevailing strike-slip movements, as supported by the available seismic data, which show strike-slip to oblique-slip focal mechanisms for the seismic events [5].
The swarm of active seismogenic normal to strike-slip faults along the axis of the mountain belt along Sannio, Irpinia, and inner Cilento sectors has been associated to most of the large historical earthquakes, which concentrate in the extensional domain. The geodetically-estimated rate of regional extension varies from ~2-3.5 [37,93] to ~4-6 mm/yr [94][95][96] and can partly explain the observed ground deformation trends in the easternmost sectors of Sannio and Irpinia ( Figure 20).
The GNNS (Global Navigation Satellite System) values available at regional scale [97] are in the range of those estimated by PSI in the same areas and show very similar patterns. The velocity field derived from GNNS surveys spanning the 1995-2011 year interval in sites located throughout southern Italy [97] appears separated in two domains separated along the eastern axis of the mountain chain and characterized by different deformation trends. Velocity patterns evidence on the Tyrrhenian side of Campania region, a W-to SW-oriented motion at a rate of 3 to 6 mm/yr referring to the Apulia frame, while only the easternmost sectors of Sannio and Irpinia have a W-directed motion at a lower rate of 1-2 mm/yr with a differential rates of about 3 mm/yr between the two sectors. Sites located in the Campanian plain show a pattern that differs from the general trend due to local effects deriving from volcanic centers (Vesuvius, Campi Flegrei, Roccamonfina) and local subsiding effects (Volturno plain) on the western Campania margin (Figure 20). The spatial discontinuity between these two horizontal velocity fields (eastward at the west side and westward at the east and south sides) roughly overlaps the Campanian plain-Apennines morphological boundary (Figure 19), which is marked by NW-SE to NNW-SSE striking faults, moving with prevailing strike-slip movements, as supported by the available seismic data, which show strike-slip to oblique-slip focal mechanisms for the seismic events [5].
The swarm of active seismogenic normal to strike-slip faults along the axis of the mountain belt along Sannio, Irpinia, and inner Cilento sectors has been associated to most of the large historical earthquakes, which concentrate in the extensional domain. The geodetically-estimated rate of regional extension varies from~2-3.5 [37,93] to~4-6 mm/yr [94][95][96] and can partly explain the observed ground deformation trends in the easternmost sectors of Sannio and Irpinia ( Figure 20).
The GNNS (Global Navigation Satellite System) values available at regional scale [97] are in the range of those estimated by PSI in the same areas and show very similar patterns. The velocity field derived from GNNS surveys spanning the 1995-2011 year interval in sites located throughout southern Italy [97] appears separated in two domains separated along the eastern axis of the mountain chain and characterized by different deformation trends. Velocity patterns evidence on the Tyrrhenian side of Campania region, a W-to SW-oriented motion at a rate of 3 to 6 mm/yr referring to the Apulia frame, while only the easternmost sectors of Sannio and Irpinia have a W-directed motion at a lower rate of 1-2 mm/yr with a differential rates of about 3 mm/yr between the two sectors. Sites located in the Campanian plain show a pattern that differs from the general trend due to local effects deriving from volcanic centers (Vesuvius, Campi Flegrei, Roccamonfina) and local subsiding effects (Volturno plain) on the western Campania margin (Figure 20    Comparison among E-W horizontal ground deformation, Global Navigation Satellite System (GNNS) data [97], and seismogenic structures in Campania.

Figure 20.
Comparison among E-W horizontal ground deformation, Global Navigation Satellite System (GNNS) data [97], and seismogenic structures in Campania.

Discussion
Although the ground deformation datasets analyzed in this study are extensive, they are not fully complete at the regional scale and some limitations need to be considered, such as short temporal coverage, limited spatial resolution, deformation occurred at a rate that was below or over the threshold deformation rates of the adopted interferometric processing. The limitations in separating E-W horizontal and vertical ground deformation have been strongly dealt with the routine post-processing of ascending and descending interferometric datasets over the same locations. The comparison with available GPS and levelling datasets allowed validating and integrating the ground deformation spatial and temporal information.
The temporal coverage of the regional C-band PSI datasets is limited to 1992-2010 and for the southernmost sectors are limited to 1992-2007 because the spatial coverage is not at full regional scale for the descending ENVISAT dataset. Therefore, the deformation that might have occurred during other time periods (before of 1991 or post 2011 to date) were not detected.
Besides, many areas of anthropogenic or natural activity may be deforming at a level below the detection threshold (−/+1 mm/yr) in the analyzed periods, so in areas without deformation that can be detected by available PSI datasets human activities and natural processes may occur at lower or too fast rates to be detected by SAR processing. Besides, some types of human activity cause significant changes (>50−100 m) to the topography (e.g., quarrying, urbanization) resulting in a phase change in the interferometric processing and in a lack of interferometric data.
Another limitation is that deformation areas must be a few hundred meters in size to be detectable with the followed approach, so smaller areas are neglected, because deformation is too small to be detected.
In spite of these limitations, the deformation signals above the detection thresholds (i.e., area and rate) have been identified throughout the regional territory through the post-processing of the available archive interferometric datasets. Following the proposed approach, the main part of the slow-rate anthropogenic and natural ground deformations occurring in Campania territory were analyzed, classified, and interpreted. Even if the C-band interferometric data are not able to fully characterize all the possible deformations, they give an important contribution to the quantitative analysis of the regional geo-and morpho-dynamics with the integration of other geodetic technologies (i.e., levelling, GPS, LIDAR, and so on).
Attribution as anthropogenic or natural deformation in previous literature was firstly considered in the analysis of the classification results. Then, in order to distinguish an anthropogenic cause of deformation from natural deformation process, the visual assessment of the area with optical imagery was performed to recognize the presence of evidence of human activity (mining, groundwater extraction, infrastructures, urbanization, etc.) together with the geological, tectonic, geomorphological, and land use analysis of the area.
In spite of the limitations of our study, several key conclusions can be drawn. The ground deformation was detected across 89.5% of the Campania territory within a variety of environmental, topographical, and geological conditions. The dominant sources of ground deformation vary with location in Campania, and so individual characteristics such as the pattern and rates of deformation are not diagnostic of the source of deformation.
Classification and analysis of ground deformation components show that the coastal plains are characterized by widespread and continuous strong subsidence signals due to sediment compaction locally enhanced by human activity, while the inner plain sectors show mainly scattered spots with locally high subsidence in correspondence of urban areas, sinkholes, and groundwater withdrawals. The volcanic sectors show interplaying horizontal and vertical trends due to volcano-tectonic processes, while in the hilly and mountain inner sectors the ground deformation is mainly controlled by large scale tectonic activity and by local landslide activity. However, deformation can have complex spatial and temporal patterns that are not always obviously related to a single human activity or natural process.
Only detailed studies at the individual sites may allow to better understand the physical processes driving ground deformation. For example, the groundwater-related deformation is the dominant cause of human-caused ground deformation in our survey (see Volturno, Sele, and Sarno plains [17,18,57,65]), and alone produces signals varying from −1 to −10 mm/yr in the vertical component of ground deformation over areas ranging in diameter up to 2-3 km, while natural subsidence due to compaction of clay-rich sediments in the subsoil is able to cause vertical rates up to −20 mm/yr over larger areas up to 25-30 km in diameter. The interplay of these two subsidence inputs locally enhances the vertical component of ground deformation to larger negative peaks.

Conclusions
The results of this study are useful for the analysis both at local and regional scales of the current rates of ground deformations due to known (or unknown) natural processes of endogenous or exogenous origin and human activities. The monitoring of the intensifications of the effects of natural hazards, due for example to subsidence in flooding alluvial areas or to sea level rise in coastal areas, or the monitoring of hazards induced to infrastructure due to severe ground deformation are also relevant tasks.
The analysis of archive interferometric data has been partly limited by poor spatial detail and data quality from C-band ERS-1/2, ENVISAT, and RADARSAT. The new opportunities for wide area mapping using data from CosmoSkyMed, Sentinel-1a/b, and TerraSAR-X would improve the results of the used post-processing approach. On-going X-band SAR satellite missions should make routine observations over the entire region to develop a more complete and detailed understanding of anthropogenic and natural deformation. Besides, ground deformation can be variable in time, so analysis of time series of observations is needed to characterize these temporal variations.
The results of this study confirm the fundamental importance of using PSI data for a comprehensive understanding of rates and patterns of recent ground deformation at regional scale also within tectonically active areas as in the Campania region.
Funding: This research received no external funding.