Vol49_1_2006def 253 ANNALS OF GEOPHYSICS, VOL. 49, N. 1, February 2006 Key words landslide movements – damage evalua- tion – differential SAR interferometry – remote sens- ing of land surface – Maratea 1. Rationale and methods Natural disasters caused by meteorological phenomena (such as hail, storms, floods, land- slides), are affecting countries of the European Union with increasing frequency, and are produc- ing heavier and heavier damage, partly hamper- ing the economic development of these countries. Timely assessment of damage caused by natural disasters and, above all, their continuous monitoring, are of fundamental importance in a modern well organised society to understand the effects of climatic changes on the social and economic fabric. In fact, a fast and synoptic da- ta acquisition is useful to guarantee rapid ade- quate intervention, and, on the other hand, it is fundamental in planning proper long-period countermeasures. In addition, loss pricing is highly interesting for institutions for socio-eco- nomic analysis, and also for insurance compa- nies. Most of such institutions and companies produce their analysis or set their annual premi- um on the basis of estimates of damage in the last period. Currently, these estimates are based on very rough damage assessment performed immediately after the disaster by central or local administrations of the various European coun- Satellite differential SAR interferometry for the evaluation of effects of hydrogeological disasters: definition of a scale for damage evaluation Vincenzo Rizzo (1) and Antonio Iodice (2) (1) Istituto di Ricerca per la Protezione Idrogeologica (IRPI), CNR, Rende (CS), Italy (2) Dipartimento di Ingegneria Elettronica e delle Telecomunicazioni, Università degli Studi di Napoli «Federico II», Napoli, Italy Abstract In this paper we present the results of monitoring soil movements over an about 10 km2 area around the border be- tween the Calabria and Basilicata regions in Italy. Monitoring has been performed using the satellite differential SAR interferometry measurements integrated with GPS measurements. In particular, we used ERS data acquired at time interval of several months (about two acquisitions per year), and after particularly strong pluvial events. Terrain dis- placement spatial and temporal analysis has been performed by employing the original method described in Berardi- no et al. (2002, 2003). Obtained results allow us to characterise unstable areas, and can be used within projects aimed at territory classification and characterisation, and at performing damage evaluation. Accordingly, this technique has been included in the framework of a proposed UE project (WEDELOP project) aimed at developing an integrated methodology to devise a damage scale for hydrogeological disasters. This scale is of great interest in many fields; in particular, definition of a damage scale would be highly desirable from the viewpoint of insurance companies. Mailing address: Dr. Antonio Iodice, Dipartimento di Ingegneria Elettronica e delle Telecomunicazioni, Univer- sità degli Studi di Napoli «Federico II», Via Claudio 21, 80125 Napoli, Italy; e-mail: iodice@unina.it 254 Vincenzo Rizzo and Antonio Iodice tries. This assessment is largely approximated, and usually direct and indirect costs for damage fixing and for productive activity re-starting tend to significantly increase with respect to the initial estimate, due to necessarily unpredictable factors. In general, objective and uniform crite- ria for damage evaluation are still lacking, also in the case of low intensity, but more frequent, events. Direct on site evaluation of damaged ob- jects is sometimes impossible, and is often very costly. The consequence of such a situation is that large-scale reliable estimates are not avail- able. By using modern well-established satellite re- mote sensing techniques, integrated with on site studies and measurements, an ad-hoc research can be conducted, aimed at defining damage scales. They will be functions of the event typol- ogy, its intensity, the urbanization degree, and so on. In order to do that, researches are needed on «sample events» used as a reference and to be in- vestigated in detail. Damaging events can be defined and coded first of all in terms of physical intensity using re- mote sensed data related to the considered event (e.g., height or duration of rain, wind strength, etc.). To these elements, some others can be added, relative to the intensity of effects on the physical environment. In fact, nowadays suffi- ciently experimented tools and methods are available for the analysis of soil changes and for the evaluation of the physical intensity of events. Physical effects at the soil produced by strong meteorological events consist of changes that are partly temporary (such as soil saturation, flood- ing, sea strength, etc.) and partly permanent (such as debris cover, terrain displacements, variation of coastline, etc.). Such effects can be detected and measured using remote sensing techniques (spectral and multispectral sensors, differential SAR interferometry, etc.), and then coded in terms of their intensity. The coding of the intensity of physical effects measured on the territory (for instance in terms of percentage of affected areas or number of damaged sites per in- vestigated area) will be then related to the coding of the event intensity. In order to obtain a better homogeneity of data, it is useful to consider a first distinction among different homogeneous physical and morphological environments (Alpi- ne or Appenninic areas, coastal area, basins ex- tension, kinds of rocks, etc.). Accordingly, a number of possibly recent and well documented «sample events» that cor- respond to the above cited coding will be chosen. This set of events will be the basis for a detailed damage estimate, limited to areas that are repre- sentative of different urbanization degrees. Ob- tained data will allow the availability of a large survey to establish a scale of relations among events, physical effects produced, and damage that allows an estimate relying on objective cod- ed criteria. This general procedure was recently sug- gested within the 6th Framework Programme («Weather Damage Evaluation and Loss Pric- ing», or WEDELOP proposal) for a large num- ber of weather risks. In particular with regard to landslide activity, permanent effects should be monitored in terms of slope movements, and a specific methodology could be assessed using remote sensing data. Fig. 1. Block scheme for the damage scale defini- tion procedure. 255 DiffSAR for the evaluation of effects of hydrogeological disasters: a scale for damage evaluation 2. Application to landslide damage evaluation Following the procedures proposed in WEDELOP we could also define a scale of damage for landslides scattered by rainy events whose elements are shown in fig. 1. The combi- nation of these elements produces a large num- ber of damage situations (and sampled events) to price: each combination of values of the ge- ological settings, event magnitude, and urban- ization degree indexes leads to one class of damage. However, to assess damage on large areas, some simplification can be made because only some factors are dominant. For example with regard to geological setting the «landslid- ing» and «morphodynamic» indexes are gener- ally related: geodynamic areas showing a con- siderable amount of landslide are usually also characterized by a percentage increase in active landslide. More properly in certain cases (where available) the indexes of morphodynam- ic activity could be represented by the percent- age of landslide remobilisation scattered by Fig. 2. The investigated area. Boxes A, B, C, and D show the locations of areas in the maps of fig. 5. 256 Vincenzo Rizzo and Antonio Iodice Table I. Interferometric pairs combined to generate the map of fig. 3. Dates of SAR acquisitions Spatial baseline Temporal baseline Percentage of selected points B⊥ [m] [days] (coherence > 0.3) 13/08/1997 31/12/1997 30 140 33.11% 31/12/1997 15/04/1998 15 105 46.74% 13/08/1997 15/04/1998 45 245 30.34% 15/04/1998 02/09/1998 110 140 30.85% 02/09/1998 09/06/1999 2 280 33.17% 02/09/1998 14/07/1999 6 315 35.28% 09/06/1999 14/07/1999 4 35 52.46% 09/06/1999 05/01/2000 6 210 32.05% 14/07/1999 05/01/2000 10 175 29.53% destabilising events. Finally, the degree of ur- banisation is very important in relation to the potential damage, and should be immediately related to the events magnitude. Geological setting (landslide and morpho- dynamic indexes for various lithological com- plexes) could be fixed on sampling areas; while differential SAR interferometry could help in assessing and setting up landslide activity, so that the morphodinamic index (D, E, F) can be updated. 3. Results of DiffSAR application on Earth slope movements performed at Maratea Differential Synthetic Aperture Radar (SAR) interferometry (DiffSAR) is able, in principle, to measure very small movements of the ground and to cover in continuity large areas (Godstein et al., 1994; Franceschetti and Lanari, 1999; Fer- retti et al., 2001), so that it can be considered an ideal tool to investigate landslides and other slope instability (Acache et al., 1995; Carnec et al., 1996). In our implementation of the DiffSAR technique, the problem of decorrelation noise is faced using a phase unwrapping approach that al- lows us to process sparse data, and the impact of atmospheric artefacts is reduced by performing a temporal analysis of the deformations observed in successive interferograms. As an example of the potentiality of this technique, we here show some results obtained by using it to improve our knowledge of the slope instability of a well investigated area (the Maratea Valley, see fig. 2) affected by continu- ous slow movements (Berardino et al., 2002, 2003). In particular, using this technique, and em- ploying the SAR acquisitions reported in table I, we analysed the time evolution of ground movements from 1997 to 2000. This time inter- val had also been explored in the past using oth- er techniques, such as distancemeter (EDM) and GPS measurements: GPS and DiffSAR re- sults were then compared. Obtained results are summarised in figs. 3 to 5. In the map of fig. 3, by visual inspection it is possible to verify the agreement between data ob- tained by the SAR technique (represented by the color scale defined in the bottom left of fig. 3) and those obtained by in situ GPS measurements (blue arrows). A 3D view of the map is reported in fig. 4, and a table of SAR versus GPS results is shown in table II. A full discussion of the com- parison between SAR and GPS measurements can be found in (Berardino et al., 2003). We here want to emphasize that the employed technique is also able to detect small scale events, possibly with antropic origin, as shown in fig. 5. Results of fig. 3 have shown that it is possi- ble to perform a temporal analysis of continuous slow landslide movements using a limited num- ber of ERS SAR data sets and low precision topographic information. All the acquired data (EDM, GPS and DiffSAR) are consistent and al- low a kinematic model of instability within the investigated time interval to be sketched. A map 257 DiffSAR for the evaluation of effects of hydrogeological disasters: a scale for damage evaluation Fig. 3. DiffSAR displacement map relative to the period from July 1997 to January 2000, and GPS data ac- quired in almost the same time interval. Displacements must be understood as projected along the SAR line of sight (about 23° from vertical and 104° from north direction, descending orbit). Fig. 4. DEM of the investigated area, showing DiffSAR displacements. The same color scale of fig. 3 is used, except for blue areas, showing slight positive vertical displacements: those are often low coherence area. 258 Vincenzo Rizzo and Antonio Iodice Fig. 5. Examples of geological validation for small phenomena, also of antropic nature, in the Maratea area (parts of maps, in scale 1:25000, of the vertical soil movements in the interval 1996-1998, only for pixels with coherence coefficient greater than 0.3). Top left (A): the Maratea dump (a few pixels large) is perfectly visible and well delim- ited, and turns out to be moving towards the sensor (i.e., upwards). Bottom left (B): last part of the Noce River, where detected upward and downward movements are coherent with erosion and deposition phenomena, with dumps, quar- ries, and so on. Top right (C): central part of the Noce Valley, where detected upward and downward movements are in reasonable agreement with low and progressive landslide movements. Bottom right (D): movements to be verified by GPS measurements, probably related to deep gravitative movements or tectonic alignments. 259 DiffSAR for the evaluation of effects of hydrogeological disasters: a scale for damage evaluation Table II. Comparison of GPS and DiffSAR measurements. Displacements are both relative to the period from June 1997 to March 2000 and are projected along the SAR line of sight. The obtained values show good agree- ment in the pixels at higher coherence. However, by considering all the points, the difference between GPS and DiffSAR measurements has a mean value equal to 2.3 cm (DiffSAR underestimates the movements) and a stan- dard deviation equal to 3.3 cm. These discrepancies can be explained by considering that in the low-coherence points the DiffSAR measurements are obtained by interpolating results obtained in surrounding high-coherence points, so that they can be considered as average displacement over a rather wide area, whereas GPS measure- ments are strictly referred to each single ground point. GPS GPS DiffSAR Coherence GPS-SAR benchmarck displacement displacement difference number (mm) (mm) (mm) 1 − 52 −54 0.3 2 2 −49 −63 0.6 14 3 −65 −63 0.5 −2 4 −52 −54 0.7 2 5 −69 −22 0.2 −47 6 −75 −24 0.6 −51 7 −18 −15 0.2 −3 8 −51 −8 0.3 −43 9 no −25 0.6 - 10 −28 −6 0.2 −22 11 −84 −56 0.2 −28 12 −13 −36 0.2 23 13 −105 −83 0.7 −22 14 −73 −52 0.3 −21 15 no −47 0.3 - 16 −32 −65 0.3 33 17 −17 −38 0.7 21 18 −112 −46 0.6 −66 19 no −14 0.4 - 20 −41 −15 0.4 −26 21 −102 −27 0.6 −80 22 no −58 0.3 - 23 −97 −51 0.3 −46 24 no −69 0.3 - 25 no −69 0.3 - 26 −104 −56 0.6 −48 27 −24 −75 0.2 51 28 −31 −18 0.3 −13 29 no 2 0.2 - 30 −47 −41 0.3 −6 31 −24 −18 0.2 −6 32 −120 −33 0.1 −87 33 −97 −13 0.2 −84 34 −83 −68 0.4 −15 35 −66 −61 0.5 −5 260 Vincenzo Rizzo and Antonio Iodice of slopes subject to different velocities and ver- tical displacements was delineated, modifying previous knowledge (Berardino et al., 2003). Within the valley a progressive and almost lin- ear displacement over time was confirmed. Dis- placement maps in figs. 3 and 5 detect areas subjected to movements, and hence the active landslides. This information can be used to set the morphodynamic index. 4. Conclusions Results of research on DiffSAR application to the monitoring of Western Basilicata (Ma- ratea and Noce Valley) show the possibility to use algorithms and operative conditions (inte- grated use of reference GPS networks) such that slow and progressive soil deformations can be detected, even in small areas and in medium-low coherence areas. However, the obtained final maps have some limitations: they do not cover the entire monitored area (in the considered case-study for about 40% of the territory the co- herence coefficient was below 0.3) and show some visual artefacts due to the sensor viewing geometry. Obtained data, however, can be inte- grated by ground based observations and by oth- er remote sensing images. These results encour- age applications aimed at studying landslide movements and erosion and deposition process- es close to rivers. Therefore, the outlined proce- dures are suitable to develop a project aimed at monitoring permanent changes produced by hy- drogeological events and at devising criteria for the definition of a relative damage scale. REFERENCES ACHACHE, J., B. FRUNEAU and C. DELACOURT (1995). Ap- plicability of SAR interferometry for operational mon- itoring of landslides, in Proceedings of the 2nd ERS Applications Workshop, London, 165-168. BERARDINO, P., M. COSTANTINI, G. FRANCESCHETTI, A. IODICE, L. PIETRANERA and V. RIZZO (2002): Differen- tial SAR interferometry for the study of slope instabil- ity at Maratea, Italy, in Proceedings of the Internation- al Geoscience and Remote Sensing Symposium, Toron- to, Canada, 2693-2695. BERARDINO, P., M. COSTANTINI, G. FRANCESCHETTI, A. IODICE, L. PIETRANERA and V. RIZZO (2003): Use of dif- ferential SAR interferometry in monitoring and model- ling large slope instability at Maratea (Basilicata, Italy), Eng. Geol., 68, 31-51. CARNEC, C., D. MASSONNET and C. KING (1996): Two ex- amples of the use of SAR interferometry on displace- ment-fields of small spatial extent, Geophys. Res. Lett., 23 (24), 3579-3582. FERRETTI, A., C. PRATI and F. ROCCA (2001): Permanent scatterers in SAR interferometry, IEEE Trans. Geosci. Remote Sensing, 39 (1), 8-20. FRANCESCHETTI, G. and R. LANARI (1999): Synthetic Aperture Radar Processing (CRC Press, Boca Raton, FL), 25-28. GOLDSTEIN, R.M., H.A. ZEBKER, C.L. WERNER, P.A. ROSEN and A. GABRIEL (1994): On the derivation of coseismic displacement field using differential radar interferome- try: the Landers earthquake, J. Geophys. Res., 99, 19618-19634. Table II (continued). GPS GPS DiffSAR Coherence GPS-SAR benchmarck displacement displacement difference number (mm) (mm) (mm) 36 −56 −38 0.2 −18 37 −37 −15 0.2 −22 38 −24 −56 0.2 32 39 −17 −9 0.2 −8 40 −25 −25 0.4 0 41 −15 −21 0.2 6 42 −9 −26 0.3 17