1. Introduction
Landslides are disasters which often occur in hilly or mountainous places in the rainy season. They are classified as a category of disasters triggered by hazards related to geological processes (earthquakes, volcanic eruptions, floods in mountainous areas) and meteorology (heavy rains, storms, or typhoons) [
1,
2]. According to a statistical report on disasters in the last decade (2006–2015) produced by the International Federation of Red Cross (IFRC), landslides account for 4.7% (177 events) of the total number of reported natural hazards [
1]. The highest density of landslide occurrences is in Asia, accounting for 118 of the 177 reported landslides occurring worldwide [
1]. A report by the The World Bank, “Natural Disaster Hotpots: A Global Risk Analysis” (2005), noted that Taiwan is a small country with an area of approximately 36,000 km
2 and a population of about 23 million people, but is classified as a high-risk country for such disasters [
3]. Taiwan ranked first of the top 15 countries (based on land area) with regard to being the most exposed to multiple hazards (affected by three or more disasters) [
3]. In other statistics examining countries which are at relatively high mortality risk from multiple hazards, Taiwan also ranked first of the top 35, with 95.1% of its population living in areas at risk [
3]. Besides earthquakes and typhoons, landslides are the most common natural disasters in Taiwan. Landslide inventory reports of The GEODAC (Global Earth Observation and Data Analysis Center) show that the total landslide area in Taiwan increased from 169.93 km
2 in 2004 to 392.44 km
2 in 2015, and the area was highest in 2010 (reaching 506.09 km
2), corresponding to the post-Typhoon Morakot period in August 2009. Most landslides here are triggered by typhoons and heavy rains, and sometimes by large earthquakes.
M. J. Crozier [
4] proposed that the factors influencing the landslide occurrence could be classified into two main categories, preparatory and triggering. The triggering group is external factors that cause immediate changes in the strength/stress of the slopes, resulting in movement, such as earthquakes, volcanic eruptions, typhoons, and heavy rainfalls, whereas the preparatory group is potential factors, such as slope failures, material composition, and coverage, which cause sliding due to gravity [
5,
6,
7,
8]. Landslides are one of the three major geological hazards (earthquakes, landslides, subsidence) occurring in Taiwan, and typhoons and earthquakes are the two main triggering factors. However, more than 10 typhoons hit Taiwan annually, and most landslides are caused by typhoons. Understanding the relationship between the rainfall thresholds that can trigger landslides and landslide susceptibility is very significant to enable better warnings of landslides during the typhoon and rainy season. Therefore, the research project is divided into two phases as follows: phase one is to produce the landslide susceptibility map using potential or preparatory factors, phase two is to find the correlation between the precipitation thresholds that can trigger landslide occurrences and the landslide susceptibility index (also called the relationship between triggering and preparatory factors). This article describes phase one of the research project.
There are many studies of landslide susceptibility modeling that have been developed in Taiwan and around the world using different approaches, such as statistical (multivariate statistics, bivariate statistics) [
7,
8,
9,
10,
11,
12,
13], probabilistic [
14], deterministic [
15], and even heuristic approaches [
16,
17,
18,
19]. These studies focus on developing models for calculating the landslide susceptibility index using various parameters of preparatory factors (slopes, lithology, land-use, and so on). Besides, some other studies also use a rainfall parameter in their models [
20,
21]. Generating a landslide susceptibility map by using the AHP method has been carried out by many researchers worldwide [
16,
17,
18,
19,
22,
23]. Most of these studies applied expert opinions for dividing the classes of each parameter (factor) and for setting the values of pair-wise comparisons (the pair of classes or parameters). The AHP helps decision-makers find the best choice and thus better understand problems. Based on mathematics and psychology, AHP was developed by Saaty in the 1970s and has been expanded and supplemented ever since [
24]. AHP provides an accurate framework for structuring a problem that needs to be resolved. It has a special application in group decision-making, and is used around the world in many critical contexts, such as government, business, industry, healthcare, education, and even disaster planning [
19,
25,
26,
27,
28,
29]. Although commonly used, the AHP often has limitations because of the inability to incorporate uncertainties and inherent inaccuracies in the data, which are related to the mapping between the perceptions and judgment of the decision maker and the precise numbers used in the AHP model. After decomposing an unstructured situation into small parts and building a hierarchy, it is important to establish pairwise comparative analysis of the importance of various indicators in decision making. In most previous AHP studies, the authors made these pairwise comparisons based on a synthesis of expert opinions or experiences [
19,
23,
25,
26,
27,
28]. However, in this study, we assign numerical values for subjective comparisons of the importance of pairwise comparisons by comparing the correlation between landslide area and distribution area of elements, using bivariate analysis and correlation analysis statistics. On the other hand, the different areas have different characteristics of landslide factors, and hence the probability of landslides is not the same. In some cases, the landslide occurrences did not follow the theoretical rules, such as the highest slope angle did not have the highest probability of landslides, or an area with lower drainage density had lower landslide probability. Beside the slope angle, the probability of landslides is also influenced by the ground material’s hardness or coverage, and so on. Therefore, it is not objective just to rely on expert opinions to set the values of pair-wise comparisons. To solve this problem, historical landslide data was used to calculate the values of pair-wise comparisons (classes or parameters) using bivariate statistics and correlation statistics. This is a new approach to apply the AHP in the calculating weights of landslide factors.
M. T. Suzen et al. carried out a survey of 145 studies on the modeling of landslide susceptibility (from 1986 to 2007), and found that four major groups of parameters were used to analyze the landslide susceptibility index, including geological parameters (geology/lithology, fault/lineaments, strata-slope interactions), topographical parameters (drainage, surface roughness, elevation, slope angle, aspect, curvature, slope length) geotechnical parameters (soil thickness, soil texture), and environmental parameters (land-use/land-cover, anthropogenic parameters, rainfall, positions within catchment) [
30]. Of these, the five most commonly used parameters were slope angle (97.9%), geology/lithology (92.31%), land-use/land-cover (75.52%), drainage (72.73%), and aspect (59.44%) [
30]. This suggests that the above five categories are important factors for landslides. In addition, N. Casagli et al. also suggested that the three key factors causing landslides were topography (slope failure), material composition (geology/lithology), and land-use [
2]. Based on the above analysis, five parameters, including lithology, slope angle, aspect, drainage density and land-use, were used to analyze the landslide susceptibility index in this study.
Study Area
The Chenyulan river watershed is located in the central part of Taiwan, in the island’s major fault zone (slate belt) (
Figure 1). This river is a tributary of the Zhoushui river. There are important fault systems in this area, such as Chenyulan, Dili, Kuling Chiao, Shenmu, and so on. The average elevation is 1680 m above sea level, and the average slope is 34.5°. The erosion activities here are complex, and occur at high densities. The area of landslides is 16,598,445 m
2, or 4.2% of the total of 392,444,312 m
2 in Taiwan. This area is also affected by the tropical monsoon climate of Taiwan, which is warm and humid all year. The rainy season is from May to October, and typhoons often occur in the summer and autumn, from June to October, more frequently in August and September, with on average four direct hits per year [
31]. The annual precipitation is more than 2500 mm. The rainfall in Taiwan often concentrates in mountainous areas, and therefore the Chen-Yu-Lan river watershed is also one of the areas with the highest rainfall in Taiwan. Summer rainfall accounts for around 80% of the annual precipitation.
The study area’s lithology can be divided into the following: metamorphic rocks are distributed in the east, sedimentary rocks are located in the west, terrace deposits (gravel, sand, clay) are scattered from the central to the south, and alluvium is concentrated along the main river in the north. Geological structures here are very complicated, with many faults and fold systems which are distributed from the north to the south of the area. The fault systems are in two main directions, including the north-south (Chenyulanchi fault, Shanshilchia fault, Dili fault, and Kulingchiao fault) and the east-west (Shenmu fault, Shihpachekeng fault, and Erhyu fault). The fold systems are arranged in the directions of the north–south (Pinglin syncline) and the north northeast–south southwest (Hoshe anticline and Tungfushan syncline). The Chen-Yu-Lan River is located on a large fault (Chenyulanchi), where a recent significant earthquake (Chi-Chi) occurred in 1999. For all these reasons, this area is one of the most vulnerable to landslides in Taiwan.
Based on a statistical report of the IFRC, the highest number of landslides occurred in the period from 2008 to 2010 (accounting for 72 of the 177 events occurring worldwide) [
1]. This period also corresponds to some very strong typhoons in Taiwan (Ex: Typhoon Kalmaegi 2008, Typhoon Morakot in 2009). As typhoons that landed on the main island are often accompanied by enormous amounts of rainfall, they often trigger shallow landslides, such as debris flows or slides, and these are two common categories of landslides in the Chen-Yu-Lan river watershed. Typhoon Kalmaegi, which formed on 13 July, brought heavy rainfall of approximately 1200 mm over a three-day period (16 July to 18 July 2009), caused a loss of lives (25 fatalities and one missing) and destruction of about 332.3 million USD of property [
31]. Typhoon Morakot, which formed on 2 August and dissipated on 12 August 2009, hitting Taiwan from the 8 August to 9 August, brought much heavier rainfall of 3000 mm caused a greater loss of lives (677 deaths, four severely injured, 22 missing, and 25 unidentified bodies) and damage of about 274 million USD. Based on sustained winds ranging from 118 to 156 km/h (Typhoon Kalmaegi: 120 km/h; Typhoon Morakot: 140 km/h), the Regional Specialized Meteorological Center Tokyo (RSMC Tokyo) ranked Typhoon Kalmaegi and Typhoon Morakot in the typhoon group of the Tropical Cyclone Intensity Scale. However, the impacts (damages and fatalities) of these two typhoons in Taiwan are very different (as mentioned above) because of the huge difference in their rainfall (Typhoon Kalmaegi: 1200 mm; Typhoon Morakot: approximate 3000 mm). Due to this difference, in this study, we considered Typhoon Kalmaegi as a normal event and Typhoon Morakot is an extreme event which triggered landslides. Landslide data related to these events are thus used to generate a landslide susceptibility map and evaluate the accuracy of the predictions of our combined model.
4. Conclusions
Producing a disaster prediction map (landslide susceptibility map) is one of the important tasks that disaster researchers are interested in and need to carry out. They always try to develop statistical analysis models so that they can have the most accurate predictions of landslides in the future. Landslide researchers have constructed many models for landslide prediction by various methods, such as multivariate statistics, binary statistics, the AHP, neutral networks, etc. Especially, in using the AHP method for landslide interpretation, researchers often use the statistics of experts’ opinion and experience to establish weights or compare the importance for pairwise indicators (classes of each landslide factor).
The purpose of this paper is to present a new approach and via AHP modeling to generate landslide prediction maps (LSMs). In this study, a combination of the AHP, bivariate, and correlation statistics model was developed. Instead of using experts’ opinions to build a weighting matrix for the classes of factors, the bivariate analysis method (considering the landslide density of each class) is used to analyze this weight. In addition, correlation statistics were also used to evaluate the correlation coefficients of the landslide factors (C) by the AHP model. Finally, a multivariate analysis model incorporating the correlation coefficients of the landslide factors was applied to produce the LSM map.
To assess the accuracy of the above integration model, the LM data (the post-event landslide inventory maps, the during-event landslide inventory maps) associated with the events was used to evaluate the binary classification method, Kappa index, and success rate curves. The results showed that for normal events, such as Kalmaegi, LSM’s success rate is high (over 77% for post-event, and over 70% for during-event), whereas for extreme events, such as Morakot, it is less accurate (over 65% for post-event, and approximately 60% for during-event). This assessment proposed that further consideration of the rainfall data in landslide predictive models is needed in future work. Therefore, in the forthcoming work, we will continue to use this approach in developing a model for calculating landslide hazard indices, where the correlation coefficient between LSI (built on preparatory factors) and rainfall (a triggering factor) will be considered in the calculation model.