Next Article in Journal
Recent Alpine Glacier Variability: Wind River Range, Wyoming, USA
Next Article in Special Issue
Architecture of Glaciotectonic Complexes
Previous Article in Journal
Health Effects Associated with Inhalation of Airborne Arsenic Arising from Mining Operations
Previous Article in Special Issue
Far-Field Deformation Resulting from Rheologic Differences Interacting with Tectonic Stresses: An Example from the Pacific/Australian Plate Boundary in Southern New Zealand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GIS-Based Landslide Susceptibility Mapping on the Peloponnese Peninsula, Greece

1
Department of Geography, Harokopio University, El. Venizelou 70, 17671 Athens, Greece
2
Department of Geological Sciences, University of KwaZulu-Natal, Westville, Private Bag X54001, South Africa
*
Author to whom correspondence should be addressed.
Submission received: 16 May 2014 / Revised: 12 August 2014 / Accepted: 13 August 2014 / Published: 20 August 2014
(This article belongs to the Special Issue Geological Mapping and Modeling of Earth Architectures)

Abstract

:
In this paper, bivariate statistical analysis modeling was applied and validated to derive a landslide susceptibility map of Peloponnese (Greece) at a regional scale. For this purpose, landslide-conditioning factors such as elevation, slope, aspect, lithology, land cover, mean annual precipitation (MAP) and peak ground acceleration (PGA), and a landslide inventory were analyzed within a GIS environment. A landslide dataset was realized using two main landslide inventories. The landslide statistical index method (LSI) produced a susceptibility map of the study area and the probability level of landslide occurrence was classified in five categories according to the best classification method from three different methods tested. Model performance was checked by an independent validation set of landslide events. The accuracy of the final result was evaluated by receiver operating characteristics (ROC) analysis. The prediction ability was found to be 75.2% indicating an acceptable susceptibility map obtained from the GIS-based bivariate statistical model.

Graphical Abstract

1. Introduction

Landslides are one of the major types of geo-hazards [1] as almost 9% of global natural disasters refer to landslides [2]. Despite advances in science and technology, these events continue to result in economic, human, and environmental losses worldwide. Landslide susceptibility (LS) is the propensity of soil or rock to produce various types of landslides [3,4]. LS is usually expressed through cartographic means. A LS map presents areas with the potential for landsliding in the future by combining some of the critical factors that contributed to the occurrence of past landslides [5]. Such a map is a valuable tool for assessing current and potential risks that can be used for developing early warning systems and mitigation plans, such as selecting the most suitable locations for construction of structures and roads.
Elevation, slope, aspect, lithology, land cover, mean annual precipitation (MAP) and peak ground acceleration (PGA) were selected as landslide-conditioning factors in our study. Among all parameters for LS zonation, elevation, slope and aspect have been recognized as the most important conditioning factors [6,7,8,9]. The elevation dataset is useful to classify the local relief and locate points of maximum and minimum heights within terrains. Generally, it is well justified through the literature [10,11] that slope gradients have a large impact on landsliding in Peloponnese. The aspect parameter is related to differential weathering, exposure to sunlight and drying winds and soil moisture [7]. Lithology also plays a key role in landslide activity since different lithologic units have different landslide susceptibility values [12]. Moreover, slope stability is also influenced by land cover. Finally, both seismicity and precipitation factors [13,14] have been used as conditioning factors in many LS zonation studies. It is, therefore, imperative to incorporate these parameters, while carrying out LS analysis in seismically active and frequently rainfall-influenced regions.
From the beginning of the 1970s, the interest of both geoscientists and engineering professionals in LS zonation and the increasing emphasis on the use of Geographic Information Systems (GIS) technology allowed for the development of many methods. According to general overviews, presented in the work of [15,16,17,18], these methods can be divided into two groups: qualitative and quantitative methods.
The qualitative methods depend on the knowledge and previous experience of the experts, and include geomorphologic analysis [19] and the use of index or parameter maps [20]. The quantitative methods depend on numerical expressions of the relationships between conditioning factors and landslide occurrence. They include geotechnical engineering approaches [21,22], statistical analysis, the artificial neural network (ANN) and neuro-fuzzy logic methods [23,24]. Some qualitative approaches, however, incorporate the idea of ranking and weighting the parameters involved, and therefore can be considered as semi-quantitative in nature [25].
In statistical analysis methods, the combination of parameters that have led to landslides in the past are determined statistically, and quantitative predictions are made for areas not affected by landslides but where similar background conditions exist [13]. Several researchers have applied these methods, which are either bivariate [26,27] or multivariate analyses [28,29]. In bivariate analysis, each individual landslide-influencing factor is combined with a landslide inventory map, and weight values based on landslide densities are calculated for its corresponding classes. A well-known and widely used bivariate analysis method for LS zonation is frequency ratio calculation, often referred as the landslide susceptibility index (LSI) [26,27,30,31,32,33,34]. Tien Bui et al. [35] evaluated and compared the results of applying the LSI and the (multivariate) logistic regression methods for estimating landslide susceptibility in the Hoa Binh province of Vietnam. Conforti et al. [36] applied the LSI method for drafting a landslide susceptibility map in the Vitravo River catchment (southern Italy). In addition, Polykretis et al. [37] compared the performance of a conventional statistical method like LSI and a soft computing method like ANNs in the Krathis and Krios drainage basins (northern Peloponnese). These three studies indicated a good prediction capability (from 73% up to 94%) for the LSI method.
The main aim of this paper was to produce a landslide susceptibility map at regional scale (1:500,000) using a quantitative method. The performance of this model was evaluated in Peloponnese peninsula, Greece. Furthermore, we implemented validation analysis to estimate the prediction ability of the proposed model.

2. Study Area

The Peloponnese peninsula is located in the southern part of Greece (Figure 1), and is connected with the mainland through the Isthmus of Corinth. The total area of Peloponnese is 21,439 km2, and its population stands at 1,086,935 inhabitants (Hellenic Statistical Agency, 2001). Agricultural, forest and semi-natural areas cover the main part of Peloponnese, whereas urban is the dominant land cover in the coastal zone of the peninsula. The climate is typical Mediterranean with a hot and relatively dry summer between June and August, and a wet season during autumn, winter and spring [38].
Peloponnese has a complex geomorphology, with mountainous inland, many coastal cliffs in the south, and basins, coastal beaches, lakes and inland basins at the western and southwestern coasts. The slopes vary from gentle to very steep, while the drainage network is well developed, and is highly controlled by fault tectonics. The study area belongs to an active tectonic zone manifested by faults, thrust zones and folds. The main rock units in the study area are (a) carbonate rocks (44%) and (b) Neogene sediments (22%). In these two types of rock units, the majority of landslides have occurred.
Figure 1. The location of the study area (Peloponnese peninsula) and the landslide inventory.
Figure 1. The location of the study area (Peloponnese peninsula) and the landslide inventory.
Geosciences 04 00176 g001
Peloponnese is a region highly damaged by the occurrence of severe natural disasters such as earthquakes, floods, landslides, and forest fires. Heavy rainfall and earthquakes have triggered several landslide events, mainly in northern and western areas of the peninsula [10,39,40]. Many serious events are related to major faults and to unstable zones located in steep slopes. In addition, the human interventions for the construction of roads have also played a key role in landslide activation.
Accordingly, the study area forms a complex physiographic region where all conditioning factors of landslides present a high spatial variability. Considering these conditions, and the fact that the past is the key to the future, it is evident that slope instability is one of the most severe hazards.

3. Data

In order to accomplish the LS analyses in the study area, a spatial database was designed and developed, and spatial analysis tools were implemented within GIS environment with the use of ArcGIS (ver. 9.3) software package. This database comprises two main parts: (a) the datasets with the background geographic conditions (slope, lithology, land cover, etc.) and (b) the landslide inventory dataset.
A spatial dataset that represents former landslides (presented as point features in the centroid of each pixel) is the most critical information layer in order to implement quantitative statistical analysis for LS assessment. In the present paper, two main landslide inventories were used: (a) an inventory maintained by Institute of Geology and Mineral Exploration (IGME) formed only from the recent historical records, covering a time period 1910 to 1995 [41]; and (b) an inventory developed on the basis of fieldwork and aerial photograph (provided from the Hellenic Military Geographical Service—HMGS, with scale 1:40,000) interpretation. This inventory contains landslide events that occurred from 1995 through 2003. The final landslide dataset consists of 282 landslides (217 and 65 landslides from the first and second inventory, respectively). This dataset (Figure 1) was randomly split into two separate groups: a training dataset (75% of the landslide inventory) and a validation dataset (25% of the landslide inventory). The training dataset was used for the implementation of statistical analysis, whereas the validation dataset was used during the verification of the results produced from the model.
The landslides under investigation have more or less the same characteristics: lateral based and downslope movement of soils or rocks. Seismic triggered landslides, occur in the vicinity of active faults, and usually related to other secondary seismic events, like soil liquefaction, subsidence of the coastal strip, and rock falls [42,43,44]. Subaqueous and liquefaction events are not included in this study. Rainfall triggered landslides usually are rapid short moving events, while slow short-moving type also occur including extensive instability zones [45]. They occur in gentle natural slopes where the translational type predominates. The occasional planar slip surfaces are located in the weathered zone of marls or flysch while ground water level reaches the surface of the slope during heavy rainfall. The most critical landslide—prone formations regarding lithology, and structure are flysch and neogene sediments, while schist and cherts significantly contribute in landslide phenomena [46]. Slides which usually take place in the gentle slope of flysch mantle are typically quite shallow and take form of a sheet of weathered zone sliding on a slip surface parallel to the ground [47]. In line with [48] and [49], in this paper the term landslide is used for translational and rotational earth slides, which were recorded in landslide inventories. These events vary consistently in volume, from some thousands of m3 to several million m3 [46], and depicting small to extremely large magnitude, according to [50] classification.
In this study, elevation, slope, aspect, PGA, MAP, lithology and land cover have been selected as the conditioning factors on landslide susceptibility. Although there are no standard guidelines for selecting these parameters [16], the nature of the study area, the scale of the analysis, and data availability were taken into account [18]. The seven factors used in current research were selected on the basis of the aforementioned criteria while literature outputs and general guidelines for GIS-based studies were also considered [7,13,18,51,52].
The Digital Elevation Model (DEM) was the key to create various topographic parameters related to the landslide activity of the study area. Here, we used the DEM from SRTM (Shuttle Radar Topography Mission) database, with cell size 90 m × 90 m. From this DEM, elevation (0–2367 m), slope angle (0°–54°) and slope aspect layers have been extracted.
The lithological layer (with cell size 250 m × 250 m) of the study area was created from the Geological Map of Greece with scale 1:500,000 (IGME, 1983). The land-cover layer (with cell size 250 m × 250 m) in the region of Peloponnese was based on CORINE program (Coordinate of Information on the Environment). The study area was classified by following the level-1 classification scheme of the CORINE data [53].
Previous studies emphasized the need of incorporating dynamic factors (PGA, MAP) with other “static” factors for LS zonation mapping in areas, whereas these factors are playing an important role not only in the reactivation of old landslides but in the development of new ones [54,55]. The precipitation data (339–1655 mm) used in this study refers to the mean annual precipitation during the period from 1950 through 1974 (source: Public Power Company—PPC, with about 30 active meteorological-rainfall stations in the study area, cell size: 250 m × 250 m). MAP is the average of the available long-term records [56]. The seismic factor (0.05–3.95 m/s2) layer was produced from the map of expected Peak Ground Acceleration (PGA) with 475-year return period (10% probability of exceedance in 50 years). PGA is the absolute maximum amplitude of recorded acceleration [57]. The source of this layer (with cell size 250 m × 250 m) was the Technical Chamber of Greece (TCG, 1992). Landslide distribution is strongly affected by seismicity and especially by ground acceleration, while magnitude—distance relations have been established for earthquake induced landslides [48].

4. Statistical Method: Landslide Susceptibility Index (LSI)

One of the most common methods used in the statistical analysis of landslides is the frequency ratio method. This method is based on the relationship between the spatial distribution of landslides and each conditioning parameter [7]. In the present paper, we estimate this relationship with the calculation of landslide susceptibility index (LSI). This method calculates the LSI for each category of all conditioning factors (e.g., land cover, lithology, slope, aspect, elevation, etc.), which are selected for the case study. Thus, supposing that j is a category within the factor i, then the LSI for this category (LSIi,j) is defined as follows [13]:
Geosciences 04 00176 i001
with Ni,j the number of landslides in category j of the factor i, Ai,j the area of this category, NT the total number of landslides and AT the total area under investigation. Thus, the LSI presents the relative susceptibility to landslide occurrence. If a category is highly correlated to landslides, the area associated with this category will have a high positive LSI value. A negative LSI value for a specific category is an indicator of low landslide density in this class. Thus, for a conditioning factor to be useful for landslide susceptibility mapping, its categories should provide a range of LSI values. The overall susceptibility S for each pixel is defined as:
Geosciences 04 00176 i002
where LSIi is the susceptibility for the factor i, and n is the total number of the factors. This analysis involves three main steps:
(a)
Categorization of all landslide- conditioning factors. In this step, the natural breaks (jenks) categorization in five discrete classes was implemented for the factors with continuous values (elevation, PGA, MAP), except for slope factor whose categorization was executed in a manual way based on its presented values. For the categorical factors (landcover, lithology, and aspect) we preserved all the classes of the nominal scale.
(b)
Calculation of the landslide density in each class and of the area for each category using GIS-based overlay functions. In this stage, the calculation of LSI for all factor categories based on Equation (1) was implemented (Table 1).
(c)
Finally, the integrated landslide susceptibility map (Figure 2) was created by combining the LSI values of multiple factors by means of GIS overlay analysis, and the Equations (1) and (2). We classified this map into five discrete categories: “Very Low”, “Low”, “Moderate”, “High” and “Very High” landslide susceptibility according to the following classification methods: (a) equal interval; (b) standard deviation and (c) natural breaks (jenks).
The standard deviation (σn), the range (rangen) and the maximum value (maxn) of LSI for all factors were also calculated (Table 2) in order to interpret the importance of each factor [13].
Table 1. Categories, landslide densities and susceptibility values (LSI) for landslide related factors.
Table 1. Categories, landslide densities and susceptibility values (LSI) for landslide related factors.
Layers (Factors)Categories (Classes)Landslide DensityLSI
Land coverArtificial surfaces0.092.04
Agricultural areas0.660.42
Forest and semi-natural land0.25−0.79
LithologyPhyllites/Gneiss (metamorphic)0.070.16
Limestones—Marbles0.20−0.77
Schists (metamorphic)0.01−0.51
Neogene0.440.69
Tertiary0.12−0.06
Flysch0.090.06
Cherts—Schists0.060.36
MAP<663 mm0.03−1.69
663–884 mm0.300.14
885–1079 mm0.430.29
1080–1295 mm0.19−0.03
>1295 mm0.05−0.22
PGA<1.16 m/s20.02−0.14
1.16–1.88 m/s20.01−0.98
1.89–2.53 m/s20.17−0.81
2.54–3.11 m/s20.670.53
>3.11 m/s20.13−0.28
Elevation<234 m0.360.17
234–524 m0.270.11
525–851 m0.22−0.03
852–1244 m0.14−0.13
>1244 m0−2.60
Slope<5°0.31−0.16
5°–10°0.300.19
11°–15°0.220.23
16°–20°0.120.08
>20°0.05−0.77
AspectNorth0.310.31
East0.21−0.15
South0.16−0.47
West0.320.16
Table 2. Landslide susceptibility index (LSI) properties of the conditioning factors.
Table 2. Landslide susceptibility index (LSI) properties of the conditioning factors.
FactorLSIminLSImaxLSIrangeLSIst. dev
Land cover−0.792.04 *2.831.42
Lithology−0.770.691.470.50
MAP−1.690.291.980.80
PGA−0.980.531.510.60
Elevation−2.600.172.771.18
Slope−0.770.230.990.41
Aspect−0.470.310.780.30
LSImin: Minimum value of LSI, LSImax: Maximum value of LSI, LSIrange: Range of LSI values, LSIst. dev: Standard deviation of LSI values. Asterisk indicates very limited area cover of this class.

5. Results—Validation

The results of the statistical (LSI) analysis are summarized in Table 1 and Table 2. The parameters with the strongest relationship with landslide occurrence were land cover (LSImax = 2.04, LSIrange = 2.83, LSIst. dev = 1.42), elevation (LSImax = 0.17, LSIrange = 2.77, LSIst. dev = 1.18) and MAP (LSImax = 0.29, LSIrange = 1.98, LSIst. dev = 0.80). For these factors the subclasses “agricultural areas”, “<234 m” and “885–1079 mm” are the most prone to landslide occurrences performing LSI values 0.42, 0.17 and 0.29 respectively. It is noted that the LSI value for the “artificial surfaces” land cover was very high (LSI = 2.04). However, this category has very limited cover in the study area. In contrast, the subclasses “forest and semi-natural land”, “>1244 m”, and “<663 mm” exhibited the lowest LSI values (−0.79, −2.60, and −1.69 respectively). The aspect was the parameter with the weakest relationship with landside occurrence (LSImax = 0.31, LSIrange = 0.78, LSIst. dev = 0.30). For the remaining factors, significant high LSI values were calculated for the following categories: (a) “lithology unit: Neogene” (LSI = 0.69); (b) “peak ground acceleration: 2.54–3.11 m/s2” (LSI = 0.53); and (c) “slope angle: 11°–15°” (LSI = 0.23). Furthermore, low LSI values were calculated for the following categories: (a) “lithology unit: limestones” (LSI = −0.77); (b) “peak ground acceleration: 1.16–1.88 m/s2” (LSI = −0.98); and (c) “slope angle: >20°” (LSI = −0.77).
The output LS map (Figure 2) from the LSI model shows that 25% (5319 km2) and 16% (3391 km2) of the study area were classified as “High” and “Very High” susceptibility zones, respectively. The same map also shows that the zones of “High” and “Very High” susceptibility are located around the northern (coastal) and western part of the study area. Finally, the overlay of the final LS map with the landslide training dataset indicated that 48%, 32% and 15% (total: 95%) of the landslide events fall within “Very High”, “High” and “Moderate” landslide susceptibility zones (in total 65% of the study area), respectively. It is notable that according to the used model only 5% of the landslide test set falls in “Low” susceptibility zone, with no landslides into the “Very Low” susceptibly zone.
We have also proceeded with a standard validation analysis using data from landslide events, not included in the initial spatial database (validation data set), in order to estimate the overall performance of the LS model in the study area. For the validation of the output from our analysis, the receiver operating characteristics (ROC) curve was drawn, and the area under curve (AUC) value was calculated for the proposed model. In practice, the AUC performs very well and is often used when a general measure of predictiveness is desired [58]. ROC analysis is considered as a powerful method for the validation of landslide susceptibility models [59,60]. The AUC value ranges from 0.5 to 1.0. The ideal model yields an AUC value close to 1.0 (perfect fit), whereas a value close to 0.5 indicates an inaccurate model (random fit).
As was aforementioned, the landslide susceptibility map was classified into five discrete categories according to equal interval, standard deviation and natural breaks (jenks) classification methods. In order to evaluate these classification methods, the results from ROC analysis for the three methods were computed and compared (Table 3). From this process, the natural breaks classification was shown to be the most efficient classification method, as it presented good overall results for landslide and non-landslide cases (high values of sensitivity and specificity, and a satisfactory AUC value). So the final map was classified based on this classification method (Figure 2).
Figure 2 shows the ROC curve of LSI model for the validation dataset. The AUC value of 0.752 indicates a reasonable prediction ability of the model.
Table 3. ROC analysis results for the classification methods.
Table 3. ROC analysis results for the classification methods.
ROC Analysis ResultsClassification Method
Equal IntervalStandard DeviationNatural Breaks (Jenks)
Accuracy64.3%62.9%65.7%
Sensitivity54.3%72.9%81.4%
Specificity74.3%52.9%50%
Positive Cases Missed321913
Negative Cases Missed183335
AUC value0.7090.7420.752
Figure 2. (a) The landslide susceptibility map produced by the LSI model and (b) receiver operating characteristics (ROC) curve for the LSI model.
Figure 2. (a) The landslide susceptibility map produced by the LSI model and (b) receiver operating characteristics (ROC) curve for the LSI model.
Geosciences 04 00176 g002

6. Discussion and Conclusions

This study was aimed at assessing the landslide susceptibility at a regional scale (1:500,000) on the Peloponnese peninsula, Greece. By applying a bivariate statistical analysis, implemented in a GIS environment, the relationships between landslide events and geo-environmental factors were assessed and shown on the susceptibility map. Despite the model, the choice of these factors plays a major role in the relative accuracy of the outcomes. Limited emphasis has been directed towards selection. The literature indicates that the most common conditioning factors are; lithological units, tectonic features, slope angle, proximity to (road or drainage) networks, land cover and rainfall distribution [3,61]. All of these factors are considered to be related to landslide occurrence. However, there are additional factors that may be arguably as influential.
In this study, seven geo-environmental factors (elevation, slope, aspect, PGA, MAP, lithology and land cover) were selected. All the factor layers were harmonized to cell size 250 m × 250 m, apart from elevation layer which was harmonized to cell size 90 m × 90 m.
Statistical methods are generally considered the most appropriate method for LS mapping at regional scales because they are objective, reproducible and easily updatable [13]. The LSI model is a data-driven bivariate statistical approach in which each parameter is analyzed individually and the calculation and application are easy and fast [62]. The bivariate statistical technique is one of the most preferable models at this scale [63]. However, many researchers argued that multivariate models—although more complex—are superior to bivariate methods as they predict the landslide susceptibility in areas with limited training data available better [64]. In data-driven LS modeling, all the problems of landslide inventory maps will automatically be inherited to the final susceptibility maps [16]. In general, a reliable, accurate and complete landslide inventory map will provide high-quality statistical models [35]. In many countries—including Greece—this kind of extended landslide inventory is not available.
The validation with the use of the empiric ROC area for the LSI model was estimated to be 0.752 for the validation dataset (Figure 3). Thus, there is 75.2% agreement between prepared LS map and landslide locations of the validation dataset, which is a reasonable result, taking into consideration the scale of analysis. Recently, LS analyses in the international literature have used ROC analysis, not only to validate the landslide susceptibility mapping models, but also to compare their prediction capabilities. Bivariate statistical analysis is a common method at this scale [64,65,66], with fair to satisfactory results (AUC values from 0.59 up to 0.76). Several researchers have also compared the LSI model with expert-based qualitative models using different data sets, and finding the LSI model to be superior [26,27,67].
In addition, a previous study [25] in the same area with the use of a semi-quantitative (expert-based fuzzy weighting—EFW) model produced rather worse results (AUC value: 0.70). Comparing the final LS map of this study with this one of [25], it is derived that in both maps 9% (1917 km2), 7% (1566 km2) and 1% (286 km2) of the study area belong to “Moderate”, “High” and “Very High”, respectively, susceptibility zones. On the other hand, the differences between two LS maps for these three susceptibility zones are presented to cover 26% (5518 km2) of the study area. The differences for the remaining susceptibility zones (“Very Low” and “Low”) cover a restricted part of the study area (6% or 1335 km2).
The implementation of the bivariate statistical model in the study area revealed that the high susceptibility areas are mainly concentrated in the western Peloponnese (Figure 2). Additionally, a linear pocket of high susceptibility values is located in the northern coastal zone of the study area. According to the final LS map, the “Very High” susceptibility zone covers a significant part of the study area (16% of the total area).
The most important factors for the LS zonation in the study area are land cover, elevation and MAP. It seems that the incorporation of dynamic factors (precipitation, seismicity) in this regional analysis was more or less beneficial to the assessment of landslide susceptibility. Most of the landslide events are detected in “agricultural areas”, in areas composed of “neogene”, in areas with elevation lower than 234 m, slope angle lower than 10°, north or west facing, moderate levels of annual precipitation (885–1079) and high seismic acceleration (2.54–3.11 m/s2).
The findings of our analysis are acceptable for this scale and better than other semi-quantitative approaches [25]. Thus, in a future work, we intend to combine this method with expert-based modeling in a “hybrid approach”.
Some limitations and assumptions of the applied method have to be pointed out. First of all, as the final result is given in a medium-scale map, the determination of the exact extent of the slope instability areas demands further site-specific research. At large scales, more exhaustive datasets and detailed geotechnical information are required. A second limitation is related to the landslide inventory dataset. At a regional scale, this dataset does not include the total amount of the landslide events within the study area. Furthermore, another limitation of this study is that the developed method does not take into account the mutual relationships between conditioning factors. Finally, it should be mentioned that, according to our analysis, the output LS map presents only the predicted spatial distribution of landsliding and not its temporal probability. Therefore, the result of this study should be used in the preliminary hazard mapping and guide quantitative risk analysis at a detailed scale.
Despite these limitations, the produced LS map could be very useful to community and local officials for choosing suitable locations for future land-use planning and implementation of developments. Additionally, planners and developers could use this map to identify roads and settlements subject to damage by future landslides, and take drastic measures for preventing the landslide events.

Acknowledgments

Authors thank the anonymous reviewers for their comments and constructive review to this paper.

Author Contributions

Christos Chalkias conceived of the study, designed the analysis, carried out its implementation and drafted the manuscript. Maria Ferentinou participated in the design and implementation of the study and in the preparation of the validation dataset. Christos Polykretis participated in the spatial database creation, and the implementation of statistical analysis. Christos Chalkias, Maria Ferentinou and Christos Polykretis participated in the revision of the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zillman, J. Report by the Director of Meteorology on the Bureau of Meteorology’s Forecasting and Warning Performance for the Sydney Hailstorm of 14 April 1999. Bureau of Meteorology: Melbourne, Australia, 1999. Available online: http://www.bom.gov.au/inside/services_policy/ storms/sydney_hail/hail_report.shtml (accessed on 28 April 2014).
  2. Gokceoglu, C.; Sonmez, H.; Nefeslioglu, H.A.; Duman, T.Y.; Can, T. The 17 March 2005 Kuzulu landslide (Sivas, Turkey) and landslide-susceptibility map of its near vicinity. Eng. Geol. 2005, 81, 65–83. [Google Scholar] [CrossRef]
  3. Guzzetti, F.; Carrara, A.; Cardinali, M.; Reichenbach, P. Landslide hazard evaluation: A review of current techniques and their application in a multi-scale study central Italy. Geomorphology 1999, 31, 181–216. [Google Scholar] [CrossRef]
  4. Crozier, M.J.; Glade, T. Landslide Hazard and Risk: Issues, Concepts and Approach. In Landslide Risk Assessment; Glade, T., Anderson, M.G., Crozier, M.J., Eds.; John Wiley: New York, NY, USA, 2005; pp. 1–40. [Google Scholar]
  5. Santacana, N.; Baeza, B.; Corominas, J.; de Paz, A.; Marturia, J. A GIS-Based Multivariate Statistical Analysis for Shallow Landslide Susceptibility Mapping in La Pobla de Lillet Area (Eastern Pyrenees Spain). Nat. Hazards 2003, 30, 281–295. [Google Scholar]
  6. Ayalew, L.; Yamagishi, H. The application of GIS-based logistic regression for landslide susceptibility mapping in the Kakuda-Yahiko Mountains Central Japan. Geomorphology 2005, 65, 15–31. [Google Scholar] [CrossRef]
  7. Magliulo, P.; Lisio, A.; Russo, F.; Zelano, A. Geomorphology and landslide susceptibility assessment using GIS and bivariate statistics: A case study in southern Italy. Nat. Hazards 2008, 47, 411–435. [Google Scholar] [CrossRef]
  8. Rahmani, M.; Hosieni, S.A.; Molla Aghajanzadeh, A.; Mashayekhan, A.; Moslemi Mahalleh, S.M. Landslide hazard zonation using the MCE method (Case study: Agh Mashhad, Mazandaran Province, Iran). Int. J. AgriSci. 2013, 3, 240–244. [Google Scholar]
  9. Ahmed, M.F.; Rogers, J.D.; Ismail, E.H. A regional level preliminary landslide susceptibility study of the upper Indus river basin. Eur. J. Remote Sens. 2014, 47, 343–373. [Google Scholar] [CrossRef]
  10. Rozos, D.; Bathrellos, G.D.; Skillodimou, H.D. Comparison of the implementation of rock engineering system and analytic hierarchy process methods, upon landslide susceptibility mapping, using GIS: A case study from the Eastern Achaia County of Peloponnesus, Greece. Environ. Earth Sci. 2011, 63, 49–63. [Google Scholar]
  11. Papathanassiou, G.; Valkaniotis, I.; Ganas, A.; Pavlides, S. GIS-based statistical analysis of the spatial distribution, of earthquake-induced landslides in the island, of Lefkada, Ionian Islands, Greece. Landslides 2013, 10, 771–783. [Google Scholar] [CrossRef]
  12. Carrara, A.; Cardinali, M.; Detti, R.; Guzzetti, F.; Pasqui, V.; Reichenbach, P. GIS techniques and statistical models in evaluating landslide hazard. Earth Surf. Process. Landf. 1991, 16, 427–445. [Google Scholar] [CrossRef]
  13. He, Y.; Beighley, R.E. GIS-based regional landslide susceptibility mapping: A case study in southern California. Earth Surf. Process. Landf. 2008, 33, 380–393. [Google Scholar] [CrossRef]
  14. Pareek, N.; Sharma, M.L.; Arora, M.K. Impact of seismic factors on landslide susceptibility zonation: a case study in part of Indian Himalayas. Landslides 2010, 7, 191–201. [Google Scholar] [CrossRef]
  15. Van Westen, C.J.; Rengers, N.; Soeters, R. Use of geomorphological information in indirect landslide susceptibility assessment. Nat. Hazards 2003, 30, 399–419. [Google Scholar] [CrossRef]
  16. Ayalew, L.; Yamagishi, H.; Maruib, H.; Takami, K. Landslides in Sado Island of Japan: Part II. GIS-based susceptibility mapping with comparisons of results from two methods and verifications. Eng. Geol. 2005, 81, 432–445. [Google Scholar]
  17. Thiery, Y.; Malet, J.P.; Sterlacchini, S.; Puissant, A.; Maquaire, O. A Landslide susceptibility assessment by bivariate methods at large scales: Application to a complex mountainous environment. Geomorphology 2007, 92, 38–59. [Google Scholar] [CrossRef] [Green Version]
  18. Yalcin, A. GIS-based landslide susceptibility mapping using analytical process and bivariate statistics in Ardesen (Turkey): Comparisons of results and confirmations. Catena 2008, 72, 1–12. [Google Scholar] [CrossRef]
  19. Listo, F.D.L.R.; Carvalho Vieira, B. Mapping of risk and susceptibility of shallow-landslide in the city of Sao Paulo, Brazil. Geomorphology 2012, 169–170, 30–44. [Google Scholar] [CrossRef]
  20. Ruff, M.; Czurda, K. Landslide susceptibility analysis with a heuristic approach in the Eastern Alps (Vorarlberg, Austria). Geomorphology 2008, 94, 314–324. [Google Scholar] [CrossRef]
  21. Pathak, S.; Poudel, R.K.; Kansakar, B.R. Application of Probabilistic Approach in Rock Slope Stability Analysis—An Experience from Nepal Disaster Mitigation of Debris Flows, Slope Failures and Landslides. In Proceedings of the INTERPRAEVENT International Symposium Disaster Mitigation of Debris Flows, Slope Failures and Landslides, Niigata, Japan, 25–27 September 2006; pp. 797–802.
  22. Thiebes, B.; Bell, R.; Glade, T. Deterministic Landslide Susceptibility Analysis Using SINMAP—Case Study in the Swabian Alb Germany. In Proceedings of the Conference—Geomorphology for the Future, Obergurgl, Austria, 2–7 September 2007; pp. 177–184.
  23. Ferentinou, M.; Chalkias, C.; Sakellariou, M. Landslide Susceptibility Mapping in National Scale and Preliminary Risk Analysis Applying Computational Methods in a GIS Environment. In Proceedings of the 9th Panhellenic Geographical Conference, Athens, Greece, 5–7 November 2010; pp. 682–689.
  24. Pourghasemi, H.R.; Pradhan, B.; Gokceoglu, C. Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz watershed, Iran. Nat. Hazards 2012, 63, 965–996. [Google Scholar] [CrossRef]
  25. Chalkias, C.; Ferentinou, M.; Polykretis, C. GIS supported landslide susceptibility modeling at regional scale: An expert-based fuzzy weighting method. ISPRS Int. J. Geo-Inf. 2014, 3, 523–539. [Google Scholar] [CrossRef]
  26. Akgun, A. A comparison of landslide susceptibility maps produced by logistic regression, multi-criteria decision, and likelihood ratio methods: A case study at İzmir, Turkey. Landslides 2012, 9, 93–106. [Google Scholar]
  27. Demir, G.; Aytekin, M.; Akgün, A.; Ikizler, S.B.; Tatar, O. A comparison of landslide susceptibility mapping of the eastern part of the North Anatolian Fault Zone (Turkey) by likelihood-frequency ratio and analytic hierarchy process methods. Nat. Hazards 2013, 65, 1481–1506. [Google Scholar] [CrossRef]
  28. Mousavi, S.Z.; Kavian, A.; Soleimani, K.; Mousavi, S.R.; Shirzadi, A. GIS-based spatial prediction of landslide susceptibility using logistic regression model. Geomat. Nat. Hazards Risk 2011, 2, 33–50. [Google Scholar]
  29. He, S.; Pan, P.; Dai, L.; Wang, H.; Liu, J. Application of kernelbased Fisher discriminant analysis to map landslide susceptibility in the Qinggan River delta, Three Gorges, China. Geomorphology 2002, 171–172, 30–41. [Google Scholar]
  30. Mezughi, T.H.; Akhir, J.M.; Rafek, A.G.; Abdullah, I. Landslide Susceptibility Assessment using Frequency Ratio Model Applied to an Area along the E-W Highway (Gerik-Jeli). Am. J. Environ. Sci. 2011, 7, 43–50. [Google Scholar]
  31. Yalcin, A.; Reis, S.; Aydinoglu, A.C.; Yomralioglu, T. A GIS-based comparative study of frequency ratio, analytical hierarchy process, bivariate statistics and logistics regression methods for landslide susceptibility mapping in Trabzon, NE Turkey. Catena 2011, 85, 274–287. [Google Scholar] [CrossRef]
  32. Abay, A.; Barbieri, G. Landslide Susceptibility and Causative Factors Evaluation of the Landslide Area of Debresina, in the Southwestern Afar Escarpment, Ethiopia. J. Earth Sci. Eng. 2012, 2, 133–144. [Google Scholar]
  33. Mondal, S.; Maiti, R. Integrating the Analytical Hierarchy Process (AHP) and the frequency ratio (FR) model in landslide susceptibility mapping of Shiv-khola watershed, Darjeeling Himalaya. Int. J. Disaster Risk Sci. 2013, 4, 200–212. [Google Scholar] [CrossRef]
  34. Paulín, L.; Bursik, G.; Hubp, M.; Mejía, J.L.; Aceves Quesada, F.L.M.P. A GIS method for landslide inventory and susceptibility mapping in the Río El Estado watershed, Pico de Orizaba volcano, México. Nat. Hazards 2014, 71, 229–241. [Google Scholar] [CrossRef]
  35. Tien Bui, D.; Lofman, O.; Revhaug, I.; Dick, O. Landslide susceptibility analysis in the Hoa Binh province of Vietnam using statistical index and logistic regression. Nat. Hazards 2011, 59, 1413–1444. [Google Scholar] [CrossRef]
  36. Conforti, M.; Robustelli, G.; Muto, F.; Critelli, S. Application and validation of bivariate GIS-based landslide susceptibility assessment for the Vitravo river catchment (Calabria, south Italy). Nat. Hazards 2012, 61, 127–141. [Google Scholar] [CrossRef]
  37. Polykretis, C.; Ferentinou, M.; Chalkias, C. A comparative study of landslide susceptibility mapping using landslide susceptibility index and artificial neural networks in the Krios River and Krathis River catchments (northern Peloponnesus, Greece). Bull. Eng. Geol. Environ. 2014. [Google Scholar] [CrossRef]
  38. Pnevmatikos, J.D.; Katsoulis, B.D. The changing rainfall regime in Greece and its impact on climatological means. Meteorol. Appl. 2006, 13, 331–345. [Google Scholar] [CrossRef]
  39. Stournaras, G.; Tsimpidis, G.; Tsoumanis, P.; Yannatos, G.; Guillande, R. Instability Phenomena in the Neogenne Deposits of Northern Peloponnese. Bull. Eng. Geol. Environ. 1998, 57, 65–68. [Google Scholar]
  40. Papathanassiou, G. Estimating slope failure potential in an earthquake prone area: A case study at Skolis Mountain, NW Peloponnesus, Greece. Bull. Eng. Geol. Environ. 2012, 71, 187–194. [Google Scholar]
  41. Koukis, G.; Tsiambaos, G.; Sabatakakis, N. Landslides Movements in Greece: Engineering Geological Characteristics and Environmental Consequences. In Proceedings of the International Symposium of Engineering geology and the Environment, Athens, Greece, 23–27 June 1997; Balkema: Athens, Greece, 1997; pp. 789–793. [Google Scholar]
  42. Papadopoulos, G.A.; Plessa, A. Magnitude-distance relations for earthquake-induced landslides in Greece. Eng. Geol. 2000, 58, 377–386. [Google Scholar] [CrossRef]
  43. Koukouvelas, I.K.; Doutsos, T. The Effects of Active Faults on the Generation of Landslides in NW Peloponnese, Greece. In Engineering Geology and the Environment; Marinos, P.G., Koukis, G.C., Tsiambaos, G.C., Stournaras, G.C., Eds.; Balkema: Rotterdam, The Netherlands, 1997; pp. 799–804. [Google Scholar]
  44. Koukouvelas, I. The Egion fault, earthquake-related and long-term deformation, Gulf of Corinth, Greece. J. Geodyn. 1998, 6, 501–513. [Google Scholar] [CrossRef]
  45. Dounias, G.; Fikiris, I.F.; Kollios, A.; Kabbadas, M.; Marinos, P. The Nemea Landslide. Trigger Mechanism—Remedial Measures. In Proceedings of the 5th Panellenic Conference of Geotechnics and Geoenvironmental Engineering, Xanthi, Greece, 31 May–2 June 2006.
  46. Sabatakakis, N.; Koukis, G.; Vassiliades, E.; Lainas, S. Landslide susceptibility zonation in Greece. Nat. Hazard. 2013, 65, 523–543. [Google Scholar] [CrossRef]
  47. Christoulas, S.; Kalteziotis, N.; Gassios, E.; Sabatakakis, N.; Tsiambaos, G. Instability Phenomena in Weathered Flysch in Greece. In Proceedings of the 5th International Symposium on Landslides, Lausanne, Switzerland, 10–15 July 1988; Volume 1, pp. 103–108.
  48. Cruden, D.M.; Varnes, D.J. Landslide types and processes. In Landslides: Investigation and Mitigation, Special Report; Turner, A.K., Schuster, R.L., Eds.; National Academy Press: Washington, DC, USA, 1996; pp. 36–75. [Google Scholar]
  49. Dikau, R.; Brunsden, D.; Schrott, L.; Ibsen, M. Landslide Recognition: Identification, Movement and Causes; Wiley: Chichester, UK, 1996. [Google Scholar]
  50. Fell, R. Landslide risk assessment and acceptable risk. Can. Geotech. J. 1994, 31, 261–272. [Google Scholar] [CrossRef]
  51. Hasekioǧullari, G.D.; Ercanoglu, M. A new approach to use AHP in landslide susceptibility mapping: A case study at Yenice (Karabuk, NW Turkey). Nat. Hazards 2012, 63, 1157–1179. [Google Scholar] [CrossRef]
  52. Das, O.; Sonmez, H.; Gokceoglu, C.; Nefeslioglu, H. A study on the effect of the ground acceleration on landslide susceptibility maps. Geophys. Res. Abstr. 2010, 12. EGU2010-5885-1. [Google Scholar]
  53. Heymann, Y.; Steenmans, C.; Croisille, G.; Bossard, M. CORINE Land Cover: Technical Guide; Office for Official Publications of the European Communities: Luxembourg City, Luxembourg, 1994. [Google Scholar]
  54. Sato, H.P.; Hasegawa, H.; Fujiwara, S.; Tobita, M.; Koarai, M.; Une, H.; Iwahashi, J. Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery. Landslides 2007, 4, 113–122. [Google Scholar] [CrossRef]
  55. Owen, L.A.; Kamp, U.; Khattak, G.A.; Harp, E.L.; Keefer, D.K.; Bauer, M.A. Landslides triggered by the 8 October 2005 Kashmir earthquake. Geomorphology 2008, 94, 1–9. [Google Scholar]
  56. Huade, G.; Wilson, J.L.; Makhnin, O. Geostatistical mapping of mountain precipitation incorporating autosearched effects of terrain and climatic characteristics. J. Hydrometeorol. 2005, 6, 1018–1031. [Google Scholar]
  57. Arjun, C.R.; Kumar, A. Artificial neural network-based estimation of peak ground acceleration. ISET J. Earthq. Technol. 2009, 46, 19–28. [Google Scholar]
  58. Fawcett, T. An introduction to ROC analysis. Pattern Recognit. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  59. Nandi, A.; Shakoor, A. Preparation of a Landslide Susceptibility Map of Summit County, Ohio, USA, Using Numerical Models. In Proceedings of the 10th IAEG Congress, Nottingham, UK, 6–10 September 2006. Paper No. 660.
  60. Frattini, P.; Crosta, G.; Carrara, A. Techniques for evaluating the performance of landslide susceptibility models. Eng. Geol. 2010, 111, 62–72. [Google Scholar] [CrossRef]
  61. Donati, L.; Turrini, M.C. An objective method to rank the importance of the factors predisposing to landslides with the GIS methodology: Application to an area of the Apennines (Valnerina; Perugia, Italy). Eng. Geol. 2002, 63, 277–289. [Google Scholar] [CrossRef]
  62. Süzen, M.L.; Doyuran, V. Data driven bivariate landslide susceptibility assessment using geographical information systems: A method and application to Asarsuyu catchment, Turkey. Eng. Geol. 2004, 71, 303–321. [Google Scholar] [CrossRef]
  63. Van Westen, C.J.; Rengers, N.; Terlien, M.T.J.; Soeters, R. Prediction of the occurrence of slope instability phenomena through GIS-based hazard zonation. Geologische. Rundschau. 1997, 86, 404–414. [Google Scholar] [CrossRef]
  64. Nandi, A.; Shakoor, A. A GIS-based landslide susceptibility evaluation using bivariate and multivariate statistical analyses. Eng. Geol. 2009, 110, 11–20. [Google Scholar] [CrossRef]
  65. Piacentini, D.; Troiani, F.; Soldati, M.; Notarnicola, C.; Savelli, D.; Schneiderbauer, S.; Strada, C. Statistical analysis for assessing shallow-landslide susceptibility in South Tyrol (south-eastern Alps, Italy). Geomorphology 2012, 151–152, 196–206. [Google Scholar]
  66. Schicker, R.; Moon, V. Comparison of bivariate and multivariate statistical approaches in landslide susceptibility mapping at a regional scale. Geomorphology 2012, 161–162, 40–57. [Google Scholar] [CrossRef]
  67. Reis, S.; Yalcin, A.; Atasoy, M.; Nisanci, R.; Bayrak, T.; Erduran, M.; Sancar, C.; Ekercin, S. Remote sensing and GIS-based landslide susceptibility mapping using frequency ratio and analytical hierarchy methods in Rize province (NE Turkey). Environ. Earth Sci. 2012, 66, 2063–2073. [Google Scholar]

Share and Cite

MDPI and ACS Style

Chalkias, C.; Ferentinou, M.; Polykretis, C. GIS-Based Landslide Susceptibility Mapping on the Peloponnese Peninsula, Greece. Geosciences 2014, 4, 176-190. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences4030176

AMA Style

Chalkias C, Ferentinou M, Polykretis C. GIS-Based Landslide Susceptibility Mapping on the Peloponnese Peninsula, Greece. Geosciences. 2014; 4(3):176-190. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences4030176

Chicago/Turabian Style

Chalkias, Christos, Maria Ferentinou, and Christos Polykretis. 2014. "GIS-Based Landslide Susceptibility Mapping on the Peloponnese Peninsula, Greece" Geosciences 4, no. 3: 176-190. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences4030176

Article Metrics

Back to TopTop