Landslides_7_4_2010
Landslides_7_4_2010
Landslides_7_4_2010
net/publication/225951175
CITATIONS READS
240 712
3 authors:
Daniele Spizzichino
Institute for Environmental Protection and Research (ISPRA)
191 PUBLICATIONS 1,495 CITATIONS
SEE PROFILE
All content following this page was uploaded by Carla Iadanza on 05 February 2018.
Abstract Landslides constitute one of the most important hazard models at a national level and for land use planning
natural hazards in Italy as they are widespread and result in and management.
considerable damage and fatalities every year. The Italian
Landslide Inventory (IFFI) Project was launched in 1999 with The Italian Landslide Inventory
the aim of identifying and mapping landslides over the entire The Italian Landslide Inventory (IFFI Project), with funding of
Italian territory. The inventory currently holds over 480,000 4.1 million euro from the Italian Government, is conducted by
landslides and has been available by means of Web services Institute for Environmental Protection and Research (ISPRA)-
since 2005. The aim of this study is to define quality indices Geological Survey of Italy and by the 21 Regions and Self-
for evaluation of the homogeneity and completeness of the Governing Provinces. ISPRA has the task of organizing and
IFFI database. In order to estimate the completeness of the coordinating the activities, producing guidelines, and disseminat-
landslide attribute information, a heuristic approach has been ing information via the Internet, while the regions and the self-
used to assign weighting values to significant parameters governing provinces collect historical documents, archive data,
selected from the landslide data sheet. The completeness and and map areas affected by landslides.
homogeneity of the landslide mapping has been evaluated by
means of three different methods: an area-frequency distribu-
tion analysis; the proximity of the landslides surveyed to urban Data collection and representation
areas; variation of the landslide index within the same The methodology applied to collect, harmonize, integrate, and
lithology. The quality indices have allowed identification of validate the landslide data is based on aerial photo interpretation,
areas with a high level of completeness and critical areas in historical documents, and field surveys. The classification of the
which the data collected have been underestimated or are not landslides is based on schemes by Varnes (1978), Cruden and
very accurate. The quality assessment of collected and stored Varnes (1996), and recommendations by the International
data is essential in order to use the IFFI database for definition Association of Engineering Geology (1990), the International
and implementation of landslide susceptibility models and for Geotechnical Societies UNESCO Working Party on Word Land-
land use planning and management. slide Inventory (1990, 1991, 1993a, b, 1994), and the International
Union of Geological Science Working Group on Landslides (1995).
Keywords Landslide . Inventory . Quality index . The IFFI integrates information from (a) national projects as
Natural hazards . GIS . Italy AVI-Inventory of information on sites historically affected by
landslides and floods in Italy for the period 1918–2000 (Guzzetti
Introduction et al. 1994), SCAI-Special project for the study of unstable towns,
The exposure of the Italian territory to natural hazards is CARG-Geological map of Italy, scale 1:50,000; (b) landslide
remarkably high. After earthquakes, landslides cause the greatest inventories by Regions, River Basin Authorities, research insti-
losses in terms of victims and damage to buildings, infrastructure, tutes, and universities; (c) river basin plans (PAI - Law 267/98);
and environmental and cultural heritage. Expansion of urban (d) Civil Defense Emergency Declarations (Law 225/92); (e)
areas and infrastructure over the last 50 years has led to a marked national and local public archives; and (f) scientific papers and
increase in landslide risk. In recent times, disastrous landslide technical reports.
events have occurred in Piemonte in 1994, in Versilia in 1996, in The information related to each landslide has been collected
Sarno and Quindici-Campania Region in 1998, in a large sector of using a landslide data sheet (Fig. 1) organized in three informa-
northern Italy in 2000, in Val Canale-Friuli Venezia Giulia Region tion levels: the first level contains basic data on landslide location,
in 2003, and in Messina in 2009 (ANPA-ARPAT 1998; ARPA type of movement and state of activity; the second level provides
Piemonte 2003; Calcaterra et al. 2000; Del Prete et al. 1998). data on geometry, geologic units, lithology, geotechnical proper-
The Italian Landslide Inventory (IFFI) Project was launched ties, land use, causes, and date of occurrence; the third level gives
in 1999 with the aim of identifying and mapping landslides detailed information on damage, investigation process, and
throughout Italy on the basis of standardized criteria. The remedial measures for risk reduction (Amanti et al. 2001).
inventory currently holds over 480,000 landslides covering an The IFFI geo-database contains geographic features of land-
area of approximately 20,500 km2, which is equivalent to 6.8% of slides (Fig. 2), an alphanumeric archive of attributes and a
the Italian territory. multimedia archive (photos, videos, documents in PDF format).
The aim of this paper is to apply methods for analyzing the Landslides are represented by a geo-referenced point (located at
quality of landslide databases and to specifically evaluate the highest point of the crown), by a polygon if the landslide may
the quality of the IFFI. This evaluation is essential for use of be mapped at the adopted survey scale or by a line when the
the database, both for the development of susceptibility and width is very narrow, as in the case of debris flows. The mapping
scale is 1:10,000 for most of Italy. However, a smaller scale area/total area × 100) calculated on a 1×1-km grid (Trigila
(1:25,000) has been used in mountainous and sparsely populated 2007).
areas. The landslide ID is a code which uniquely identifies each
landslide within the whole Italian territory. Web services for the dissemination of information via the internet
The Italian Landslide Inventory currently holds 482,272 The IFFI Inventory was published on the Internet in 2005
landslides. Figure 3 shows the Landslide Index (landslide (http://www.sinanet.apat.it/progettoiffi) with the aim of dissem-
inating landslide information to central and local public The WMS service of the IFFI has been developed according
administrations, research institutes, geologists, engineers, and to European Directive INSPIRE-Infrastructure for spatial informa-
citizens. Two services have been designed: the WebGIS tion in Europe-2007/2/CE (The European Parliament and the
application Cart@net-IFFI and the WMS service. Cart@net-IFFI Council of the European Union 2007) and the Open GIS
allows geographical data to be explored and detailed informa- Consortium specifications (Open Geospatial Consortium 2006).
tion to be obtained on the most important landslide param- This interoperable service enables the user to build a customized
eters. The application manages large amounts of spatial data map overlaying the IFFI landslides with other layers stored on the
with extremely rapid response times and integrates the vector user’s own computer or provided by different Internet map
datasets, served by ArcSDE and ArcIMS Server (ESRI) and servers. The user can access the WMS service through a WMS
Oracle DB, with the imagery served by the Image Web Server viewer, such as the one provided by the INSPIRE Geo-portal, or
(ER Mapper) in ECW (Enhanced Compressed Wavelet) format. through a Desktop GIS (e.g. ESRI ArcGIS®).
The user can view the landslides of the IFFI together with
other vector layers (e.g., urban areas—Corine Land Cover Quality assessment of the Italian landslide database
2000, roads and railways, administrative boundaries, and
drainage network) and raster layers (digital terrain model, Quality of landslide inventories
digital orthophoto TerraItaly it2000, Landsat satellite images, The quality of a landslide inventory depends mainly on the
and IGM topographic map; Fig. 4). method used to create the inventory. For example, the use
merely of aerial photo interpretation does not allow the photos (Carrara et al. 1995; Guzzetti 2002; Wills and Mc Crink
landslide date of occurrence to be determined and makes it 2002). In order to set up a more complete inventory, it is
difficult to identify dormant landslides which have been necessary to use a multi-temporal series of aerial photos
significantly remodeled by farming activities, covered by thick which enable reconstruction of the developmental dynamics of
vegetation, or are located in urban areas. The correct the landslides, in addition to allowing a better identification
identification and mapping of the landslides depends not and classification. With regard to the use of historical and
only on the degree of conservation of the morphological archive data, the greatest limitations are found in the
evidence (slopes, tensile cracks, etc.) but also on the imprecise description of the landslide location and the non-
experience of the photo interpreter and the scale of the aerial homogeneous spatial distribution of the landslides due to the
underestimation of events in sparsely populated areas and, With regard to the quality of the IFFI, while the data
therefore, with a lower number of manmade structures were being checked, all geometrical features and attribute
(Guzzetti et al. 1999). In addition, while there is a great deal information passed a series of controls before being stored in
of information on the damage, data on the type of landslide the geo-database. In particular, two sets of quality control
and its dynamics are often lacking (Carrara et al. 2003). The tasks were performed: the first set identified omissions,
integration of the historical data and the photo interpretation duplicate values in the ID field, digital mapping anomalies,
analyses enables the abovementioned limitations to be over- geo-positional errors, etc. The second set inspected the IFFI
come and to significantly increase the quality of the with other datasets: AVI (Inventory of information on sites
inventory. historically affected by landslides and floods in Italy for the
To date, there are no well-established scientific criteria in period 1918–2000), SCAI (Special project for the study of
the literature for evaluating the accuracy, reliability, and unstable towns), CARG (Geological map of Italy, scale
completeness of a landslide inventory. This evaluation is 1:50,000), and PAI (River basin plan).
carried out in some studies by comparing inventories This study has carried out a further evaluation of the
prepared for the same area (Ardizzone et al. 2002; Carrara completeness of the information contained in the alphanu-
et al. 1992; Galli et al. 2008), while in other studies, the meric database, and the homogeneity and level of detail of
frequency-area statistics are determined and analyzed. the landslide mapping, as described below.
Table 1 (continued)
Parameters Level of IFFI data sheet Landslides (No.) Landslides (%) Weight
Investigations and monitoring 3 1,523 0.32 3
Remedial measures 3 5,278 1.09 3
Photos – 5,841 1.21 3
Attached documents – 522 0.11 10
Fig. 6 Relationship between the rating of completeness Very high and the number of landslides for each Italian Region
order to evaluate the completeness of the alphanumeric previously existing inventories but also to the greater or lesser
archive of the IFFI, four classes of values have been defined: degree of use of the aerial photo interpretation, field surveys, or
Low for values <10; Medium 10–15; High 15–21; Very high ≥ 21 historical and archives data in the methodology adopted by the
(Table 2). individual Regions/Self-Governing Provinces. The data relative to
It has been seen from the analysis carried out at a three Regions of Southern Italy (Basilicata, Calabria and Sicily)
national scale that 75% of the landslides have a Medium value are underestimated compared with the actual instability situation
with respect to the completeness of the alphanumeric data- since, to date, the landslide survey activities been concentrated
base. At a regional scale, it has been seen that the quality mainly in urban areas or areas with main road and rail
depends on the number of landslides which populate each infrastructures.
regional database, such as, for example, in the case of The completeness and homogeneity of the landslide
Sardinia (Fig. 6). mapping has been evaluated in this study by means of three
different analysis methods: (a) area-frequency statistics; (b)
Completeness and homogeneity of the IFFI landslide mapping analysis of the proximity of the landslides surveyed to urban
There is a good level of detail in the landslide mapping for the areas; (c) analysis of the variation of the landslide index within
IFFI Project: over a total of 482,272 landslides surveyed, 377,016 the same lithology.
(78%) are represented by a polygon, 50,651 (11%) by a line and
only 53,236 (11%) by a point. If the distribution of the landslides Frequency-area statistics
throughout Italy is analyzed (Fig. 3), it is seen that there is a The use of this method is based on the assumption that landslides,
certain lack of homogeneity of the landslide data within the just like other natural hazards (e.g., earthquakes, forest fires), show
Alpine Arc due not only to the different levels of detail of the fractal or scale invariant size statistics, i.e., the cumulative
% Landslides
dNCL Region/ Quality
(distance > 750 m
¼ C0 Ab
L ð2Þ Self-Governing Province Index
dAL from urban areas)
Piemonte 44.5
with b ¼ a þ 1 (α exponent of the cumulative distribution). Valle d'Aosta 59.9
In order to obtain the noncumulative area-frequency distribu- Lombardia 66.4
tion starting from the cumulative distribution, use may be made of
Bolzano-Bozen 49.2
the derivatives method, which calculates the derivative of NCL(A) as
Trento 54.5
the gradient of the cumulative curve which interpolates a
determined number of points (generally 5) around a particular Veneto 32.6
value of A (Frattini et al. 2003). The area-frequency distribution Friuli-Venezia Giulia 63.7
may also be estimated by using the logarithmic binning method Liguria 32.3
(Guzzetti et al. 2003; Hovius et al. 1997, 2000). The data sample is Emilia Romagna 51.0
subdivided into intervals (bins) with constant width along a Toscana 52.9
logarithmic axis, that is, keeping the difference logsi–logsi−1 Umbria 59.6
constant, where si and si−1 are the ends of the ith interval). Since Marche 51.5
the bin widths grow with increasing landslide size, the logarithmic Lazio 39.1
binning method overcomes the problem that large landslides are
Abruzzo 43.2
generally very few.
Molise 71.6
The empirical frequency is calculated as the ratio between
the number of landslides which fall within each interval and Campania 47.4
the width of the interval itself. A diagram is prepared in the Puglia 62.2
bilogarithmic plane of si snii1 versus bsi , where ni is the number Basilicata 25.9
of landslides which fall within the ith interval and bsi is the Calabria 12.6
average value of the interval. In the case of the power law, the Sicilia 54.7
noncumulative distribution obtained by means of the loga- Sardegna 61.0
rithmic binning method is equivalent to a cumulative distri-
bution, so the gradient of the straight line which interpolates
the data in the bilogarithmic plane is equivalent to the
exponent α of the cumulative distribution. Quality Index
The area-frequency distribution is described by the High Landslide(%) > Mean + St. Dev.
power law of a certain interval. Below a certain area value, Medium Mean - St.Dev. < Landslide(%) < Mean + St.Dev.
there is a drop in frequency, known as rollover (deviation
Low Landslide(%) < Mean - St. Dev.
from power law scaling). The rollover which may be seen in
all the frequency curves shown in the literature is inter-
preted by some authors as real, that is, due to the physical Mean = 49; Standard Deviation = 14,5
The area-frequency distribution for the 377,016 landslides Self-Governing Provinces (min. area, max. area and most
of the IFFI represented by a polygon has been determined frequent area; Table 3) enabled the adopted work scale and
by means of the logarithmic binning method. The data the level of detail of the mapping resolution to be evaluated.
sample was classified in intervals with a constant width of The most complete inventories at a regional level correspond
log si log si1 ¼ 0:05. to values less than the most frequent landslide area.
The area-frequency distribution follows a power-law, with
the straight line interpolated in the bilogarithmic plane having Proximity to urban areas
a gradient of 2.07, which is equivalent to an exponent 3.07 of Landslide inventories are more accurate in areas close to
the noncumulative distribution (Fig. 7). The exponent urban centers or linear communication infrastructures since
obtained has an excellent correspondence with the values in the historical information linked to the damage caused by
the literature. Studies carried out on numerous inventories for landslides is more numerous, more detailed studies are
landslides throughout the world, in a multitude of climatic carried out as part of the territorial plans (e.g. River Basin
and physiographic regions, have shown a range of values for Plans - PAI) and access to the sites is easier for performance
the exponent b ¼ 2:5 0:5 (Guzzetti et al. 2002). of the ground surveys. The spatial distribution of the land-
The power-law scaling is valid for four area magnitudes, slides in the IFFI Inventory was, therefore, evaluated on the
from a maximum area of 40 km2 to a minimum area of 1× basis of the distance from the urban areas. For this purpose a
10–2 km2, which corresponds to the rollover. The most 750-m wide buffer was generated around the Italian urban
frequent area in the inventory is 2,990 m2. areas extracted from the Prior 10 k DB (Intesa GIS Project -
At a national scale, the inventory may be deemed to be scale 1:10,000) and the percentage of landslides which fall
complete for landslides with areas greater than 10,000 m2. The outside the buffer was calculated. Three quality classes were
analysis of the area-frequency distributions for the 21 Regions/ then defined on the basis of the mean and standard deviation
Fig. 8 Landslide mapping in Molise Region: the figure shows a full coverage of the territory: urban areas (in red); landslides far from urban areas (in blue); landslides
inside the 750 m buffer (in cyan)
of the 21 regional values (Table 4). It may be seen from the Landslide density and lithology
analyses that four regions have a Low quality index, which is Lastly, the spatial completeness of the IFFI was evaluated by
indicative of a non-uniform territorial coverage of the analyzing the variation of the landslide index within the
investigations, 14 Regions have a Medium index, and three same outcropping lithology. The Lithological Map 1:500,000
have a High index, as in the case of the Molise Region (Fig. 9), derived from the Geological Map of Italy by
(Fig. 8). Geological Survey of Italy (scale 1:500,000-1983), was used as
input by grouping the 128 geological units into 13 lithological include territories at an elevation >600 m (Trigila and
classes. The analysis was carried out within the mountainous- Iadanza 2008).
hilly territory (Fig. 10) which is the one in which landslides Figure 11 shows the landslide index calculated for the lithological
may potentially occur. This territory was obtained by using a class “Clays and clays & marls” in the mountainous-hilly territory.
20×20-m DEM. The alluvial plains include territories with The particularly low landslide index values for the “Clays and clays &
an elevation <300 m and slope <3°; the hills have a slope >3° marls” in Basilicata, Calabria and Sicily are indicative of an
or an elevation between 300 and 600 m; the mountains underestimate of landslides in these areas (Fig. 12).
Fig. 11 Landslide Index (%) related to clays, clays, and marls in mountainous-hilly territory
Fig. 12 Landslide Index % (calculated over a 250-m wide grid) related to clays, clays, and marls in Basilicata Region
exponent which has a good correspondence with the values susceptibility models. The analyses have also detected critical
in the literature. The frequency distribution differs from the areas in which the data collected have been underestimated
power-law (rollover) at areas equal to 10,000 m2. With regard or are not very accurate and for which it will be necessary to
to the analysis of the proximity of the landslide surveys to integrate the data in the future. However, the aims of the
the urban centers, it was seen that the territorial coverage of IFFI Project may be considered fully achieved bearing in
the investigations is almost uniform for the majority of the mind that in 1999, before the start of the Project, approx-
Regions. Comparison of the landslide indices for areas with imately 70,000 landslides had been recorded by the regions
the same outcropping lithology has enabled identification of and self-governing provinces.
critical areas in which the data collected has been under- The IFFI, in terms of the number of landslides surveyed, the
estimated, as in the case for the Basilicata Region. full coverage of the Italian territory, and the detailed landslide
The quality indices have allowed identification of repre- mapping, represents a capable tool for hazard assessment and
sentative areas in terms of homogeneity and reliability of the land use planning. The IFFI Project methodology could be
IFFI Project data, which may be used as training and test adopted at the European level for the realization of a unique
areas for the implementation of quantitative landslide and homogeneous inventory of landslides in Europe.