Investigating uplift in Lesina Marina (Southern Italy) with the aid of persistent scatterer SAR interferometry and in situ measurements

We apply persistent scatterer interferometry (PSI) techniques to synthetic aperture radar (SAR) data from ERS and ENVISAT satellites on the Lesina Marina area, a coastal tourist village in Apulia, Southern Italy, where the excavation of a canal exposed grey micro- and meso-crystalline gypsum which is now showing a high density of cavities and sinkholes due to gravitational collapse processes. We observe PS objects undergoing displacements, along the sensor line of sight, forming the same relatively smooth pattern in all the processed data stacks. Vertical displacement rates, derived through integration of ascending and descending geometries, reach about 4 mm/year on locations adjacent to the canal, gently decreasing towards the western end of the built-up area. High-precision leveling measurements, performed in 1999 and 2010, reveal a substantial agreement with the ENVISAT PSI data, taking into account a small bias due to the choice of the leveling reference point. The dataset, thus validated, suggests the presence of an uplift phenomenon going on steadily for the entire timespan covered by the SAR observations (1992--2009). These observations, supported by petrographic data and in situ investigations, seem only in part compatible with a residual diapirism, and hint instead to more complex processes, such as a combination of diapirism and the hydration of the residual anhydrite in the core of the gypsum mass. These results confirm the importance of the integration between in situ, geologic and geophysical, remotely sensed investigations, as the latter often represent an essential tool to infer whether a given phenomenon, which can be hypothesized by the former, is presently under development.

section_content: The village of Lesina Marina is sited between the Lesina Lagoon and the mouth of the Fortore River, on the coastal plain that extends along the northern Apulia Adriatic coastline (southern Italy, see Fig. 1a ). Standing only a few meters above the sea level, the village is a renowned tourist leisure area. It lies not far from the Pietre Nere (Black Stones) Point, a very interesting site and the only outcrop of magmatic rocks in the Apulia region and in the whole Adriatic coast of Italy. 

In the area, an exotic rock body basically composed of Triassic gypsum crops out on a very gentle ridge, on its sides surrounded and covered by relatively thin Quaternary sandy deposits (e. g. Carella 1963 ) (Fig. 1b ). The area of Lesina Marina is characterized by a tectonic uplift since Middle Pleistocene (Ciaranfi et al. 1983 ). An average uplift rate of ca. 1. 5 mm/year was calculated by Mastronuzzi and Sanso `(2002) for the Holocene. 

The excavation of a canal through this area, performed in 1930, exposed grey micro-and meso-crystalline gypsum with intercalations of black limestones and marls of Upper Triassic age, mantled by loose sandy Quaternary deposits. Besides their uniqueness in the surface geology of southern Italy (e. g. Martinis and Pieri 1964), the importance of these gypsum rocks is also related to the large number of cavities, either dissolutional conduits or voids related to gravitational collapse processes, and cover suffosion sinkholes that have formed in the gypsum karst of Lesina Marina area, dangerously affecting the residential area during the last two decades (e. g. Fidelibus et al. 2011). The presence of the sinkholes and of the karst reactivation led to the demand of a monitoring system, for the detection of vertical displacements over a large area. To this end, SAR data acquired by the European ERS and ENVISAT satellites, spanning more than 10 years, were processed through persistent scatterer interferometry (PSI) techniques (Ferretti et al. 2001; Bovenga et al. 2004 ). These allow investigating millimetric displacements of the Earth surface occurring over long timespans, thanks to the regular revisit times and the uniform acquisition geometry of the used remote sensing satellites. Complete time series of displacement measurements can be derived over particular objects on the Earth surface, exhibiting radar backscattering signatures which stay coherent over time. 

The analysis was initially aimed at investigating the possible evidence of subtle movements connected to the known sinkholes in the area. In fact, applications of InSAR and/or PSI to sinkhole movement monitoring has been attempted in some instances, with success (Vaccari et al. 2013; Rucker et al. 2013; Gutie ´rrez et al. 2011; Carbonel et al. 2014). 

The most striking results of the processing showed instead the presence of steadily uplifting PS objects over the Lesina Marina urban area, in a rather smooth spatial pattern, with maximum values located along the Acquarotta canal, with average uplift velocities reaching about 4 mm/year, and decreasing gradually towards W-SW to about zero. The peculiar nature of this pattern, with its small spatial scale, and its position on such a geologically complex and hazard-prone area, prompted us to look for a comprehensive explanation of the phenomena possibly taking place at the Lesina Marina site. 

In the paper, after a description of the experimental data and methods, in ''Dataset'' section, we show and discuss in ''Results and discussion'' section the main results obtained from data processing and comparisons. Finally, in ''Conclusions'' section we draw some conclusions.
section_name: Introduction

section_content: During Mesozoic, in the Adria Plate (sensu Channell et al. 1979), the narrow Adriatic pelagic basin, flanked by the two carbonate Apulian, to the SW, and Adriatic Platform, to the NE, roughly occupied the same position of the present-day Adriatic Sea (Zappaterra 1990 (Zappaterra, 1994;; Bernoulli 2001 ), as sketched in Fig. 1a. Both the basinal and platform domains developed as a consequence of early Jurassic rifting of an epeiric area dominated by the deposition of carbonates. This carbonate platform-basin system, which developed mainly during the Mesozoic, was rooted on Norian anhydrites and shallow-water limestones and dolostones, i. e. the Burano Fm, which overlie Permian continental deposits draping the Hercynian basement (Ricchetti et al. 1988). Later, the basinal and platform domains were partially involved in the Tertiary shortening related to the Apennines and to the Dinarides-Albanides. Near and along the front of these orogens, the evaporites belonging to the Burano Fm promoted diapirism, mostly during Neogene (Underhill 1988; Zelilidis et al. 1998; Kamberis et al. 2000; Kokinou et al. 2005; Scrocca 2006; Alves et al. 2007; Geletti et al. 2008; Festa et al. 2014; Kokkalas et al. 2013). 

In northern Apulia (Fig. 1a ), in the Lesina Marina area, the cropping out of Triassic gypsum rocks, belonging to the Burano Fm (Fig. 1b ), rose up from the deep anhydrite source (Cotecchia and Canitano 1954; Bigazzi et al. 1996). According to Bigazzi et al. (1996), hydration of anhydrites, i. e. the formation of gypsum, possibly occurred after their uprise. Up to decametric blocks of Paleocenic mafic and ultramafic hypabissal rocks (De Fino et al. 1983; Bigazzi et al. 1996) and Triassic limestones (e. g. Posenato et al. 1994), that locally crop out embedded in the gypsum rocks (e. g. Amendolagine et al. 1964), would have been dragged, from the depth, during ascent by diapirism of the gypsum body (Cotecchia and Canitano 1954) or by tectonic squeezing of the anhydrites (Bigazzi et al. 1996). As b Fig. 1 a Schematic structural map of the Periadriatic region around the Adriatic Sea (modified from Zappaterra 1990 Zappaterra, 1994)). The Meso-Cenozoic paleogeographic position of the Adriatic Basin between the Apulian and Adriatic carbonate platforms is shown. The fronts of Apennines and Dinarides according to Scrocca (2006) and Fantoni and Franciosi (2010), respectively. The halokinetic structures are also indicated (modified from Geletti et al. 2008 ). b Schematic geological map of northern Apulia, around the area of Lesina Marina village and of the Tremiti Islands (modified from Boni et al. 1969; Cremonini et al. 1971 ) the location of the underlying halokinetic structure called Tremiti diapir is also indicated (modified from Festa et al. 2014). c Schematic geological map of the Lesina Marina village area. d Leveling measurements map observed in the exposed rocks by Selleri and Mastronuzzi (2003) and by Fidelibus et al. (2011), the gypsum bedrock in the subsurface of Lesina Marina, where its top gently dips toward West (Melidoro and Panaro 2000; Fidelibus et al. 2011; Caggiano et al. 2012 ), (Fig. 1c ), shows a high density of cavities, either dissolutional conduits or voids related to gravitational collapse processes. Most of these cavities, metric in size, are partially or totally filled basically with sand deposits derived from the overlying Quaternary sedimentary cover (Melidoro and Panaro 2000). According to Selleri and Mastronuzzi (2003), these characters may be related to a probable karst undercover genesis and evolution. 

Recently, Fidelibus et al. (2011) have demonstrated that the changes in the hydrogeological functioning have favoured both erosion of the filling of cavities and karstification processes. Furthermore, hydrochemical evidence reveals that gypsum dissolution is a currently active process favoured by fresh water and sea water mixing and cation exchange processes; however, the dominant evolution time of gypsum dissolution is much greater than human lifetime (Campana and Fidelibus 2015).
section_name: Geological setting

section_content: We used SAR data from the European Space Agency's ERS-1/2 and ENVISAT satellites, from both ascending (24 ERS, 33 ENVISAT acquisitions) and descending geometry (32 ERS, 28 ENVISAT acquisitions), covering a total time interval from July 1992 to October 2000 (ERS), and then from September 2003 to December 2009 (ENVISAT). The data were processed through the DORIS open source software (Kampes et al. 2003) to obtain stacks of co-registered differential interferograms, and then through the SPINUA software suite (Bovenga et al. 2004 ) to extract relevant information about displacements of stable points located on the area. 

Tables 1 and 2 summarize data relevant to the used SAR scenes, such as the date of acquisition and the perpendicular component of the sensor baselines, i. e. the distance separating satellite orbit tracks in a plane perpendicular to the satellite trajectory from that of the reference (master) acquisition. 

PSI techniques generally rely on the density of detected stable objects to extract information about spatially correlated phenomena, such as the atmospheric phase screen (APS) superimposed on the phase signal of each SAR image in a time series (Ferretti et al. 2001). When analysing scarcely urbanized areas, or, such as in this case, relatively isolated, small urban centers, it is useful to use robust approaches in extracting APS information to increase the number of detected PS. This is generally possible by restricting the processed area, so that APS variations can be modeled with smooth functions such as low-order polynomial surfaces. The approach has been demonstrated to give good results in several instances, e. g. in the study of small-scale phenomena such as landslides in low-urbanization settings (Bovenga et al. 2006 ). Here, we applied the same approach to increase the number of detected PS on the urban center of Lesina Marina, which has a size of about 1 km Â 1 km. 

To reduce the chance of processing artefacts such as spurious velocity gradients and/or offset velocity differences, the local results were compared to global, full-frame results obtained through the same or analogous software, and to GPS data from permanent stations distributed throughout the area. The comparison, not shown here for brevity, does not show significant artefacts, thus confirming the reliability of the measurements described below.
section_name: Dataset Satellite data and processing

section_content: High-precision leveling measurements were performed on several tens of topographic points on the study area, during two ground truth campaigns, in 2000 and in 2010. The leveling measurement positions can be seen in the map in Fig. 1d, as well as on the enlarged detail map in Fig. 5. The leveling stations are located mainly close to the canal, in a strip roughly oriented in a N-S direction, about 150 m wide. All leveling measurements are expressed with respect to a given reference, corresponding to the area of the ''Maddalena'' Hotel, shown as a magenta star in Fig. 1d.
section_name: Leveling

section_content: Field geological mapping, sampling and petrographic observations were carried out in the Lesina Marina area. Since the Quaternary sandy deposits are very poorly exposed, although their presence is well constrained by many holes drilled in the area (Melidoro and Panaro 2000; Fidelibus et al. 2011), the geological mapping focused on Upper Triassic gypsum rocks, well exposed especially along the northern part of the Acquarotta canal (Fig. 1c ). 

Here, the gypsum rocks appear from mainly coarse grained massive to locally finer grained layered. Within the massive portion, an intergranular lutite is often present between the larger gypsum crystals, which are euhedral to subhedral. The individual grains, that typically are colourless, vitreous and characterized by clear cleavage planes, reach a maximum size of 2-3 cm (Fig. 2a ); however, white gypsum is locally present as well (Fig. 2b ). Centimetric to metric in size clasts of both Paleocenic mafic to ultramafic hypabissal rocks and Triassic limestones, as well as of finer grained layered gypsum rocks, locally crop out embedded in the massive gypsum matrix (Fig. 2b-d, respectively ). In addition, a layering in the gypsum massive portion is sometimes highlighted by some centimeters-thick layers of black shales (Fig. 2e ). Observations of thin sections, under optical microscope in polarized transmitted light, reveal for all the samples (Fig. 2c ) that gypsum crystals, usually characterized by low birefringence with first-order grey interference colours are represented by poikilocrystals with small inclusions of anhydrite, typically showing bright, second-order interference colours (Fig. 2f ). It is worth to note that very fine anhydrite crystals can be also found as inclusions in the finer grained gypsum rocks. Therefore, these relationships between gypsum and anhydrite indicate that gypsum is secondary and that it formed at the expense of the original anhydrite crystals by hydration processes. In finer grained gypsum rocks the layering is well developed and characterized by an alternation of whitish to dark greyish millimeter-to centimeter-thick layers that exhibit a plane parallel geometry (Fig. 2g ); sometimes a chaotic convolute geometry is shown by this layers (Fig. 2h ). Finally, a network of centimeter-thick fractures filled by secondary gypsum (Fig. 2j ) crosscuts both the massive and the layered gypsum portions.
section_name: Field and petrographic observations

section_content: For all the processed SAR data stacks, we observe a good number of PS objects quite uniformly distributed over the village area, with a few more scattered PS in the outskirts. PS displacement are measured along the satellite line of sight (LOS), i. e. the sensor looking direction, which is determined by the nominal incidence angle, of about 23 degrees with respect to the vertical for all the images, and the orbit heading angle with respect to the N direction, which is about 8 and 172 degrees for ascending and descending geometries, respectively. The displacement LOS mean velocities, estimated over the ERS and ENVI-SAT observation periods, appear spatially homogeneous, gently decreasing from positive (uplift) to slightly negative (subsidence) values in a roughly WSW direction, with a few scattered points at the west of the urban area subsiding up to À7 mm/year. The PS mean velocity maps for the four processed data stacks are shown in Fig 3a-d. 

The described velocity spatial trend is present in the PS maps derived from both ascending and descending geometries, in ERS and ENVISAT data. In view of the above-mentioned controls, which allow to exclude processing artefacts, the presence of such a steady uplift signal in both time series periods appears striking, and claims for a geophysical interpretation. 

Relying on the relatively smooth nature of the observed phenomenon, average velocities from ascending and descending PS data were interpolated through Kriging (Isaaks and Srivastava 1989 ) over a common geographic grid, with 30 m posting, then combined to obtain vertical (up-down) and east-west velocity components, denoted by v U h i and v E h i, according to the following equations:

where v A h i and v D h i are the ascending and descending interpolated mean velocity fields, respectively, and h is the SAR sensor look angle. The third, north-south component of movement is hardly retrievable from satellite InSAR data, due to the geometric configuration of the spacecraft orbit, which causes the LOS to have very small north-south projection (e. g. Samiei-Esfahany et al. 2010 ). In our case, north-south motion can be assumed negligible, due to the geographic configuration of the test site, with the Acquarotta canal running roughly in a N-S direction. Therefore, this component of displacements is ignored in the following. 

Derived vertical displacement rates for ERS data (see Fig. 3e ) reach $ 4 mm/year on locations adjacent to the canal, gently decreasing to slightly negative values towards the western end of the built-up area (about 600 m away). Similar values are derived from ENVISAT data (Fig. 3f ), with somewhat lower maximum values on the canal, and a slightly more pronounced decrease towards the West. The east-west velocity component for the ERS data (Fig. 3g ) also show a W-SW trend, with rates going from about zero close to the canal, to about À1 to À2 mm/year (i. e. towards the West), on the westwards margin of the town area. Slightly lower values are reached in the ENVISAT data (Fig. 3g ).
section_name: Results and discussion

section_content: Mean displacement velocities derived from each leveling data point were compared with the nearest PSI measurements within a fixed search radius of 50 m. Mean ascending and descending PSI velocities within this search radius were combined to derive the vertical and horizontal motion components, as in (1). The v U h i values derived in this way were then compared to the leveling measurements. Plots of the two data series are shown in Fig. 4, as a function of the northings UTM coordinate. In the figure, PS velocities' error bars are derived from estimates of the phase dispersion (Bovenga et al. 2006 ). Uncertainties on leveling mean velocities can only be inferred from nominal accuracy values of each of the two measurements, and are negligible compared to the former. As can be seen, apart from a vertical offset due to the different chosen reference point, the two trends look similar. Differences are contained within AE2 mm/year, with a mean error of about 1. 5 mm/year and a standard deviation of about 0. 54 mm/ year. The correlation between the two measurements is about 0. 62, thus the hypothesis that PS vertical velocity is statistically independent of leveling velocity can be ruled out with a high level of confidence ([99 % using Fisher's test). These figures seem in agreement with those reported in several comparisons of PS-derived velocity measurements with ground-based leveling, see e. g. Perissin et al. (2007), or Crosetto et al. (2007). A few examples of time series of PS points close to leveling stations, including one possibly corresponding to a sinkhole collapse site, are reported in the following sections. This validation of the ENVISAT data allows to safely confirm also the trends observed from the ERS data, which show basically similar characteristics, thus depicting a regional uplift phenomenon going on in a rather steady way over the total *17-year timespan of SAR observations.
section_name: Comparison with leveling data

section_content: Figure 5 shows sample time series of PS close to leveling stations, represented as coloured triangles, with colour coded according to the measured average displacement rates. PS displacements are defined relative to the position in the master acquisition; the leveling deformation lines were also shifted vertically so as to refer them to the master are observed (see Fig. 3 ) we convert LOS-projected PS individual velocities to vertical velocities by assuming a null horizontal velocity component, then using the transformation v U ¼ v A=D = cos h, as can be deduced by (1) by assuming v E ¼ 0. 

As can be seen, in all cases the agreement of the PS time series with the leveling displacements looks satisfactory. The oscillations in the PS time series displacements, enhanced by the small scale of the plots' y-axis, are due in part to the absence of any filtering in the temporal dimension. Apart from this, no significant departures from a linear displacement behaviour are visible for all the points.
section_name: PS time series

section_content: In Fig. 6 we show a zoom map around a pair of detected PS close to a single leveling measurements which reported a decreasing height of up to 28 mm in 10 years. The PS points closest to these subsiding leveling benchmarks exhibit a slow downward movement, of about À1 mm/ year. Compared to the surrounding, positive uplift velocity values, around 2-3 mm/year, this means a relative subsidence (negative) rate of about 3-4 mm/year. 

Displacement time series of the mentioned two PS, as well as a relative mean displacement of the nearby leveling benchmark, are shown in Fig. 7. The time series seem to indicate a relatively fast subsidence during the years 2005 and 2006, followed by a relatively more stable behaviour. This seems compatible with the total displacement measured by the leveling. 

Due to the known widespread presence of several sinkholes in the Lesina Marina area, it is tempting to associate these localized, isolated displacement time series with the precursory movements associated with a nearby, possibly undetected sinkhole. However, no other clear instances of moving PS close to sinkholes could be reliably determined on the site. The most likely explanation is that in the period between 1999 and 2008, the distribution of newly detected sinkholes reveals a trend towards a clustering with old ones, rather than a random positioning (see e. g. fig. 8 in Caporale et al. 2013). This is in line with the general rock mass behaviour, as gathered by numerous direct and geophysical surveys carried out, in terms of rock fracturing and subsequent anisotropy of permeability (Spilotro et al. 2009), and the resulting groundwater flow (Campana and Fidelibus 2015). As a consequence, areas affected by sinkholes in the period covered by the SAR data were left free of man-made urban structures, and are therefore lacking PS. The only exception is the abovementioned single cluster of PS, belonging, as described, to a building sufficiently close to anomalous measurements of the leveling data, which however lack direct evidence of nearby sinkholes. 

In fact, sinkholes usually develop on small timescales, over limited spatial areas, and are therefore considered at the limit of applicability of spaceborne remote sensing techniques. Some recent examples of applications of SAR interferometry to sinkhole detection and/or monitoring can be found, e. g. in Bruno et al. (2008 ), or Gutie ´rrez et al. (2011). It is usually argued that sinkhole phenomena occur suddenly in landscapes characterized by rocks or other Fig. 6 ENVISAT ascending (coloured circles) and descending (coloured squares) PS, over a cadastral map of the area in the vicinity of a detected sinkhole. Larger circles in green are centered on points where leveling detected subsidence (circle radius proportional to subsidence rate) kinds of substrate with limited deformation capabilities. Nevertheless, the chance of detecting small precursory terrain movements in correspondence to sinkholes can be increased, when the sinkhole phenomenon occurs on a site with a cover of sand, clay or other compactable material. In such cases, in fact, small precursory episodic deformations, per se not detectable by remote sensing, could cause stronger compaction phenomena in the soft cover, which could be detected as small movements preceding the final catastrophic collapse.
section_name: Local analysis

section_content: From a geological point of view, these results are not of immediate interpretation. The investigated area affected by the uplift (Ciaranfi et al. 1983; Mastronuzzi and Sanso 2002 ) coincides with the higher part of a salt dome (Cotecchia and Canitano 1954). Therefore, according to the conclusions made by Festa et al. (2014) about the near Tremiti diapir (Fig. 1b ), a still active diapirism cannot be excluded also for the diapir of the Lesina Marina area. An uplift rate of ca. 1 mm/year can be estimated for the Tremiti diapir, since, according to Festa et al. (2014), a distance upwards of ca. 5000 m was covered by this diapir in the last ca. 5 Ma. This uplift rate shows some correlation with that inferred by Mastronuzzi and Sanso `(2002) for the Holocene, nearby Lesina Marina, i. e. 1. 5 mm/year. However, our data point to uplift rates well above that value, i. e. up to about 4 mm/year. Accordingly, a further contribution to uplift, in addition to that due to diapirism, may be related to the transformation of anhydrite in gypsum (Fig. 2f ). As demonstrated in several studies (Pettijohn 1957; Bell et al. 1986; Azam 2006; Rauh and Thuro 2007), this transformation involves a volume increase of up to ca. 60 % of the hydrating mass. Therefore, the hydration of the residual anhydrite surface could have determined, and could continue to determine, a volume increase of the gypsum mass in the area of Lesina Marina. This phenomenon may have been accelerated by an increase in circulation of the water in the fractured bedrock due to the Acquarotta canal, along which the maximum uplift rates have been recorded, especially in its northern part (Fig. 3e, f ). Here, the horizontal movements are practically absent, whereas they have been recorded, together with the subsidence, at the western side of the area with maximum uplift rates (Fig. 3g, h ), as would be expected. More quantitative evaluations or modeling about the volumetric expansion of the sulphate rock mass would have to rely on at least the amount of anhydrite/gypsum transformation, as well as on the geometrical 3D reconstruction of the sulphate rock body. Unfortunately, these quantitative evaluations would be highly speculative due to both the high heterogeneity of these sulphate rocks (i. e. random variations of volume percentages of anhydrite and gypsum) and the not currently available deep data. In this respect, as a final remark, it is worth noting that the anhydrite/gypsum transformation has been observed to cause consistent and sudden uplift phenomena on other sites, such as the city of Staufen, Germany (Lubitz et al. 2014), causing serious infrastructure damages.
section_name: A geological Interpretation

section_content: We show results of the monitoring of the Lesina Marina area in Southern Italy, a tourist village plagued by the formation of sinkholes in the vicinity of an artificial canal cutting the area since the thirties. ERS and ENVISAT ascending and descending SAR data, processed through PSI techniques, show displacement velocities degrading spatially from a maximum, positive value in a strip close to the canal, towards null values at the western part of the urban area, with a few scattered subsiding points. Highprecision leveling, performed in two ground campaigns in 2000 and 2010, was compared to the PS-derived vertical velocity maps from ENVISAT data. 

Comparison between leveling mean velocities and PS mean vertical velocities reveals a good overall agreement, compatible with current state of the art PSI accuracy assessments. The velocity differences exhibit a small-dispersion distribution of the residuals, apart from a bias due to the different reference points in the two techniques. The similarity of the ERS data series thus points to a steady uplift phenomenon covering the entire $ 17-year timespan of SAR observations. The available dataset has been tentatively interpreted as the result of anhydrite/gypsum transformation that is still provoking a volume increase of the hydrated gypsum mass and, accordingly, uplift rates well above those attributable to the solely diapiric phenomenon. 

It should be noticed that the available leveling stations are located mostly close to the canal path, disposed in a roughly N-S direction, thus perpendicular to the direction of main variation of the uplift phenomenon. Further measurements along this latter direction would probably add information to this study. Moreover, recent and higher resolution SAR data will likely help in shedding light on the actual processes taking place in this interesting site. The petrographic observations cannot alone reveal if the anhydrite/gypsum transformation is still active. As demonstrated in the present case study, SAR interferometry and in situ measurements may represent an essential tool to infer if the hydration of anhydrite is under development. Therefore, the integration between both petrographic and geophysical investigations should be used to influence the engineering design choices for the risk mitigation in those urbanized areas characterized by the presence of sulphate rocks.
section_name: Conclusions

section_content: Acknowledgments Research funded through the Italian Ministry of Research-PRIN 2008 research grant ''Advanced technologies in the assessment and mitigation of the landslide risk: precursors detection, previsional models and thematic mapping''. ENVISAT and ERS data are copyright of European Space Agency -provided under Category-1 project N. 5367, ''Subsidence monitoring in Daunia and Capitanata (Puglia Region, Italy) through multi-temporal point-target DInSAR techniques''. GAP srl personnel is acknowledged for aid in InSAR processing. 

Open Access This article is distributed under the terms of the Creative Commons Attribution 4. 0 International License (http://crea tivecommons. org/licenses/by/4. 0/),which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.