Next Article in Journal
Effects of Long-Term High Carbonate Alkalinity Stress on the Ovarian Development in Exopalaemon carinicauda
Next Article in Special Issue
Model Test and Numerical Simulation of Slope Instability Process Induced by Rainfall
Previous Article in Journal
A Modified GALDIT Method to Assess Groundwater Vulnerability to Salinization—Application to Rhodope Coastal Aquifer (North Greece)
Previous Article in Special Issue
Study on the Risk Assessment Method of Rainfall Landslide
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Analysis of the Potential of Landslides Induced by Combination between Rainfall and Earthquakes

1
Department of Land Management and Development, Chang Jung Christian University, Tainan 711301, Taiwan
2
Department of Business Administration, Chang Jung Christian University, Tainan 711301, Taiwan
3
Department of Aviation and Maritime Transportation Management, Chang Jung Christian University, Tainan 711301, Taiwan
4
Disaster Prevention Research Center, National Cheng Kung University, Tainan 701401, Taiwan
*
Author to whom correspondence should be addressed.
Water 2022, 14(22), 3691; https://doi.org/10.3390/w14223691
Submission received: 21 October 2022 / Revised: 10 November 2022 / Accepted: 11 November 2022 / Published: 15 November 2022
(This article belongs to the Special Issue Rainfall-Induced Geological Disasters)

Abstract

:
This study analyzed the potential of landslides induced by the interaction between rainfall and earthquakes. Dapu Township and Alishan Township in Chiayi County, southern Taiwan, were included as study areas. From satellite images and the literature, we collected data for multiple years and time series and then used the random forest data mining algorithm for satellite image interpretation. A hazard index for the interaction between earthquakes and rainfall (IHERI) was proposed, and an index for the degree of land disturbance (IDLD) was estimated to explore the characteristics of IHERI under specific natural environmental and slope land use conditions. The results revealed that among the investigated disaster-causing factors, the degree of slope land use disturbance, the slope of the natural environment, and rainfall exerted the strongest effect on landslide occurrence. When IHERI or IDLD was higher, the probability of a landslide also increased, and under conditions of a similar IDLD, the probability of landslides increased as the IHERI value increased, and vice versa. Thus, given the interaction between rainfall and earthquakes in the study area, the effect of the degree of slope land use disturbance on landslides should not be ignored. The results of a receiver operating characteristic (ROC) curve analysis indicated that the areas under the ROC curve for landslides induced by different trigger factors were all above 0.94. The results indicate that the area in which medium–high-level landslides are induced by an interaction between rainfall and earthquakes is large.

1. Introduction

Taiwan is an island located in an earthquake zone. Taiwan frequently experiences various natural disasters, such as typhoons and earthquakes. Furthermore, global climate change has adversely affected Taiwan in recent years by causing frequent extreme rainfall, short-term heavy rainfall, and disturbances in slope land use; these factors can result in landslides and substantial changes in mountainous areas. Therefore, studies should investigate the relationship between landslides caused by the interaction of factors that trigger rainfall and earthquakes, particularly those related to the target region’s natural environment and disturbances in slope land use, and the scale of landslides. The index of the interaction effect of the dual trigger factors and the degree of disturbance in slope land use should be determined; they can serve as preliminary standards that can be used to effectively regulate the degree of slope land use and development and thus prevent slope land disasters.
Some scholars have used a genetic adaptive neural network and performed texture analysis to interpret satellite images and investigate the relationship between land use and rainfall-induced landslides; these approaches can be beneficial for landslide monitoring [1,2,3]. Feng et al. [4] reported that the combination of a random forest (RF) with texture analysis exhibited a higher accuracy in mapping land covered by vegetation in cities than did the traditional maximum likelihood method. Stumpf et al. [5] indicated that the interpretation and classification of satellite images can be helpful for analyzing big data, obtaining cadastral information, and determining land cover and vegetation types and soil properties.
Zhong [6] used the 1999 Chi-Chi earthquake as a boundary and selected two typhoons before and after the earthquake to investigate the effect of the interaction between the earthquake and rainfall on slope stability. They determined that before and after the earthquake, the number and areas of landslides were positively correlated with the maximum rainfall intensity and the slope of a landslide caused by higher rainfall was lower than that of a landslide caused by lower rainfall. In addition, Lin et al. [7] evaluated the impact of the Chi-Chi earthquake on subsequent rainfall-induced landslides and determined that the rainfall-induced landslides were mostly distributed between 40° and 50° and between 30° and 40° after and before the earthquake, respectively. Tang et al. [8] reported that after the Wenchuan earthquake, subsequent torrential rain caused a 30% increase in landslide areas. Moreover, landslide points were more numerous after the earthquake than they were before the earthquake. Huang et al. [9] suggested that the use of seismic factors can increase the prediction accuracy for landslide occurrence. Chen et al. [10] determined that extreme rainfall events caused frequent landslide erosion, with extreme rainfall (i.e., a maximum 24-h rainfall level of >600 mm) accounting for 64%–79% of the average landslide erosion rate. Yang et al. [11] conducted a sensitivity analysis and employed the analytical hierarchy process and factor weighting to evaluate and compare the occurrence of landslides after and before earthquakes. They determined that landslides caused by strong earthquakes exhibited more spatial clustering. Jan et al. [12] demonstrated that the impact of typhoon rainfall depends not only on the amount of rainfall but also on its intensity. Li et al. [13] observed that new landslides tended to occur in low-altitude or low-slope areas under normal rainfall and in high-altitude areas under high rainfall. Tseng et al. [14] reported that the number and area of landslides induced by rainfall were positively correlated with the degree of land disturbance. Wistuba et al. [15] determined that earthquakes caused up to 50% of landslide events, either alone or in combination with rainfall, but earthquakes were rarely the sole trigger of landslides. Chen et al. [16] reported that the main trigger factors for landslides were cumulative rainfall in areas that were more prone to landslides and rainfall intensity in areas that were less prone to landslides. Extreme climate conditions have caused increases in the magnitude and frequency of landslides in both types of areas that are generally less and more prone to landslides. Valagussa et al. [17] analyzed three earthquake events and discovered that surface disturbance affects the intensity of landslides; the higher the intensity of a landslide is, the greater the surface disturbance is. Moreover, they proposed that in addition to surface disturbance, terrain type and lithology affect the intensity of landslides. Quesada-Román et al. [18] investigated the combined effects of earthquakes and rainfall on landslides and debris flow. They observed that the density of landslides was higher in areas closer to earthquake epicenters and that more rainfall and landslides occurred in areas with higher slopes and elevations. Ruggeri et al. [19] demonstrated the possibility to distinguish the difference between seismic and rainfall induced displacements of the slope. Bontemps et al. [20] reported that moderate-intensity earthquakes with a magnitude of >5.5 on the Richter scale increased the probability of landslide occurrence and rainfall events either before or after the earthquake and therefore led to more severe landslides. A combination of low-intensity earthquakes with a magnitude of 3.2–3.6 on the Richter scale and heavy rainfall reduced the stability of slopes. Liu et al. [21] analyzed the debris flow records of the Wenchuan, Lushan, and Jiuzhaigou earthquakes and determined that the critical values of postearthquake landslides were smaller than those of pre-earthquake rainfall.
Artificial intelligence and decision trees are commonly used for data mining when big data are analyzed. Breiman et al. [22] reported that an RF is a combination of decision tree classifiers; each tree is influenced by independently sampled random vector values, and the distribution is the same for all decision trees in a forest. Colditz et al. [23] reported that the RF algorithm yielded the best results for the area ratio assignment of the training samples for each class. Belgiu et al. [24] demonstrated that an RF is a classifier that aggregates several decision trees and can handle high data dimensionality and multiplicity. Lagomarsino et al. [25] Taalab et al. [26] reported that an RF can process a large amount of data. In earth science and landslide sensitivity research, an RF is considered an advanced and mature technology. Catani et al. [27] demonstrated that the dimension of parameter space, the mapping unit and the training process strongly influence the classification accuracy and the prediction results. Goetz et al. [28] presented a comparison of traditional statistical and novel machine learning models applied for landslide susceptibility modeling. Steger et al. [29] explored discrepancies between the predictive performance of a landslide susceptibility model and the geomorphic plausibility of subsequent landslide susceptibility maps. Chen et al. [30] used kernel logistic regression, an RF, and an alternating decision tree to build a potential map of groundwater sources. They observed the highest area under the curve (AUC) value and degree of accuracy for the RF.
Ercanoglu [31] used a neural network in combination with a geographic information system (GIS) to draw a potential map of slope land hazards on the basis of landslide locations identified on aerial photos and six parameters, namely slope, aspect, elevation, terrain, a humidity index, and a vegetation coverage index. Shahabi et al. [32] combined remote sensing and GISs to apply statistical models for the delineation of landslide-sensitive areas. They considered the following factors: slope, aspect, elevation, lithology, normalized difference vegetation index, vegetation, rainfall, distance from the fault, distance from a water system, and distance from a road. Tseng et al. [14] employed the time range before and after a typhoon that had recently caused a road slope landslide in the study area as a benchmark and interpreted images to determine surface changes before and after the occurrence of landslides to establish a landslide potential assessment model. Then, they used a GIS platform to plot a landslide potential map.
The present study used an RF algorithm to interpret satellite images, quantitatively analyzed landslide potential and hazard factors induced by the interaction between rainfall and earthquakes, and explored the characteristics of landslides induced by the interaction between rainfall and earthquakes in a specific natural environment under slope land utilization and development conditions in the region.

2. Research Methods

2.1. RF

An RF is a multifunctional machine learning algorithm. The algorithm contains numerous decision tree classifiers without any relationships. An RF can handle a large amount of input data, smooth out errors in class determination, and maintain accuracy by classifying data when data are insufficient or missing [33]. An RF selects a subset to predict the growth of a decision tree and determine the growth of each tree in accordance with the bootstrap training set. The bootstrapping method involves repeatedly sampling a limited sample with replacement to establish a new sample that can adequately represent the distribution of the population [34,35]. The randomly sampled bootstrap samples and randomly selected input factors are used to construct multiple CARTs (Classification and Regression Tree) in accordance with the maximum depth, and generate forests with the lowest generalization error. Each CART in the random forest grows independently and uses different random samples of bootstraps. Therefore, each CART has no relationship. After the establishment of the random forest, when a new sample is input, each decision tree in the forest judges separately and the result is dependent on majority voting. When creating a random forest, there is no over-fitting problem [22,35,36]. The RF method is widely used in different fields and has provided satisfactory results in studies analyzing the potential of landslide occurrence [37,38]. In addition, the RF method is suitable for data with many variables or insufficient resolution, and it does not lead to the problem of overfitting [22].
In the present study, we used the Train Random Trees Classifier (ESRI, 2021) in ArcGIS (10.5.1) to interpret and classify satellite images, and the RF module [39] developed in the R language was used to establish a landslide potential assessment model.

2.2. Texture Analysis

Texture is an image characteristic, and it can be analyzed to distinguish groups of images. Mathematically, texture is defined as the correlation of the grayscale value or color space of adjacent pixels in an image or the visual representation of image grayscale values and color changes with spatial positions, including edges, shapes, stripes, and color blocks [40]. This study used the image processing software ERDAS IMAGINE to analyze image texture. Four aspects of textures can be analyzed using ERDAS IMAGINE, namely variance, skewness, kurtosis, and mean Euclidean distance [41].

2.3. Accuracy Assessment

To determine whether the accuracy of the satellite image interpretation results met the requirement for land use classification data, we adopted the error matrix method and evaluated the accuracy of the image classification. After the error matrix analysis and calculation are completed, several commonly used accuracy evaluation indicators can be calculated, including producer’s accuracy (PA), user’s accuracy, overall accuracy (OA), and coefficient of agreement (Kappa index) [42,43]. Kappa values are between 0 and 1; the larger the Kappa value is, the higher the classification accuracy is [44].

2.4. Receiver Operating Characteristic Curve

A receiver operating characteristic (ROC) curve can be used to describe all the possible combinations of correct and incorrect [45]. The advantage of using an ROC curve is that its value is not affected by the number of classified samples and can be used to evaluate the recognition of patterns. For ROC curves, true positive (TP) rates (TPR) are used as ordinates and false positive (FP) rates (FPR) are used as abscissas. To plot an ROC curve, the TPR, true negative (TN) rate, FPR, and false negative (FN) rate should be calculated. The FPR represents the ratio of all negative samples that are wrongly judged as positive, as presented in Formula (1). The TPR represents the ratio of all positive samples that are correctly judged to be positive, as presented in Formula (2).
F P R = F P / ( F P + T N )
T P R = T P / ( T P + F N )
This study used the AUC to determine the advantages and disadvantages of the evaluation model. If the area is assumed to be obtained in the space of 1 × 1, the AUC must be between 0 and 1. The higher the AUC value is, the better the discriminative ability is, and vice versa [45].

3. Study Areas

This study was conducted in Dapu Township and Alishan Township in Chiayi County, southern Taiwan (Figure 1). Typhoon Morakot occurred in Taiwan in August 2009 and caused record-breaking disasters. The accumulated rainfall in Alishan Township reached 3060 mm and resulted in 619 deaths and 76 missing individuals. At 1:43:3.2 pm on 2 June 2013, an earthquake with a magnitude of 6.5 on the Richter scale and a depth of 14.5 km occurred 29.3 km east of the Nantou County Government headquarters; it killed 4 people and injured 19 people. This earthquake also caused rail bending, rockfalls, and landslides.
Alishan Township is located in the eastern part of Chiayi County, Taiwan. This township has the largest area and lowest population density in Chiayi County. Dapu Township is located in the southeast corner of Chiayi County. This township has the second lowest population density in Chiayi County. Alishan Township and Dapu Township in Chiayi County are located to the north and south of the Tropic of Cancer, respectively. Alishan Township is mainly located in the middle of the Central Mountains. Its altitude varies considerably. The terrain is high in the east and low in the west. The altitude is approximately 360–3997 m. The temperature mostly ranges from 6.3 °C to 15 °C. The average minimum and maximum temperatures are 7.2 °C and 15.5 °C, respectively. Dapu Township is located in the Alishan Mountains, with an altitude of approximately 230–1100 m. An alpine terrain, hills and plains, and a basin are present in the east, west, and center, respectively. The annual average maximum temperature is 28 °C. The annual precipitation in Alishan Township and Dapu Township is approximately 4186 and 2236 mm, respectively [46]. The main river flowing through the study area is Zengwen River, which is the fourth longest stream in Taiwan. Its source is 2609 m above sea level, with a total length of 138.5 km and a drainage area of 1176.7 km2. Its main tributaries are the Tanaku River and Puyanu, Caolan, Houbori, Cailiao, and Guantian streams. Zengwen River is rich in water resources and flows through the Zengwen, Nanhua, and Wushantou reservoirs. This river is also used as a water supply and for power generation and attracts tourists to the local area [47].

4. Potential of the Induction of Landslides by the Interaction between Rainfall and Earthquakes

4.1. Interpretation and Classification of Satellite Images before and after Rainfall or Earthquakes in the Study Area and Extraction of Landslide Data

In this study, we collected Formosa Satellite 2, SPOT-5, SPOT-6, and SPOT-7 images with image resolutions of 6 × 6 M, 8 × 8 M, and 10 × 10 M, respectively. We used ERDAS IMAGINE to fuse various images into a multispectral image and then performed image positioning. Subsequently, the editor function of ArcGIS was used to delineate the parts covered by clouds or shadows in the study area and then crop and remove the clouds. To ensure the integrity of the satellite images before and after each rainfall or typhoon in the study area and reduce the impact of cloud cover, we first distinguished between the time periods before and after each occurrence of rainfall or each earthquake. For each time period, if a preprocessed image had a blank block because of cloud removal, the Mosaic to New Raster function in ArcGIS (10.5.1) was used to map and merge preprocessed satellite images belonging to the same time period and with similar image dates to obtain a satellite image with low or no cloud coverage. According to the method reported by Chen et al. [10], we used ERDAS IMAGES [41] to analyze the texture of the aforementioned multispectral images to improve the accuracy of interpretation and classification.
With reference to Chen et al. [1], we used seven classifications, namely water, roads, buildings, streamways, bare ground, green cover, and cash crops (fruit trees and tea plants), as interpretive classifications. The ArcGIS image classification module (segmentation and classification toolset) in the random tree classifier (train random trees) was used for the classification and interpretation of the satellite images. To confirm the accuracy of the collected satellite images after interpretation, we used the aforementioned coefficient of agreement (Kappa) and OA as the bases for accuracy evaluation and conducted random sampling of the selected images. Seven interpretation classifications were determined, and 25 points (interpretation grids) were randomly selected as checkpoints. High-resolution aerial photos and on-site survey data were used for cross comparison. The Kappa values and OA of the satellite image interpretation and classification results before and after different rainfall and earthquake events in the study area are summarized in Table 1. The average Kappa value and OA of all 37 satellite image interpretation results were 0.66 and 70.1%, respectively, and the overall image classification accuracy was medium to medium high.
To accurately determine the locations of landslide areas in the study area, the images of the aforementioned rainfall and earthquake events were used to manually identify the locations of exposed areas. The bare ground in the images before and after an earthquake event was subtracted to determine the location of the landslide area. Among the types of slope failure, debris slides were the easiest and most reliable type to be identified on satellite images since vegetation was effectively stripped off from the slopes. Therefore, debris slides are the major landslides mapped in our study.

4.2. Selection of Landslide Hazard Factors

This study analyzed three environmental hazards that can induce landslides, namely the natural environment, slope land use disturbance, and trigger factors. Elevation, slope, slope aspect, geology, distance from the water system, and distance from the fault were initially selected as natural environmental hazards. Road density, building density, crop density, green coverage, and bare land density were selected as hazard factors for slope land use disturbance. Rainfall and earthquake events were included as trigger factors. The product of effective accumulated rainfall (EAR) and the maximum 3-h rolling rainfall intensity (I3Rmax) were used as the rainfall indicator. Peak ground acceleration (PGA) was used as the earthquake indicator. The methods for estimating and grading the various hazard factors were in accordance with those reported by Tseng et al. [14] and are explained in the following.

4.2.1. Natural Environmental Factors

  • Elevation
We used ArcGIS Spatial Analyst to determine the digital elevation mode (DEM) in the study area and ArcGIS Zonal Statistics to estimate the average elevation value in each basic grid of the study area. The elevation was divided into seven grades, with an interval of 250 m, and then graded and coded. The smaller the code was, the higher the elevation value was.
B.
Slope
We used the DEM of the study area and the slope analysis tool (Slope) in ArcGIS to estimate the average slope value of each basic grid and divided the slope into six grades, with an interval of 10%. The sixth grade slope was coded as 1, and the first grade slope was coded as 6. The smaller the code was, the greater the slope was.
C.
Aspect
According to the method reported by Tseng et al. [14], we used the ArcGIS aspect analysis tool (Aspect) to estimate the average aspect value for each basic grid. Because the slope aspect of the land with exposed mudstone in southwestern Taiwan mainly faces south, south was coded as 1 and southeast and southwest were coded as 2. The aspects were divided into five grades, and the flat land was coded as 0.
D.
Distance from the river
According to Moayedi et al. [48], an area more than 700 m away from a river is less likely to experience a landslide. Hong et al. [49] classified the distance of an area from a river into the levels of 100 m or less, 100–300 m, 300–500 m, 500–700 m, and above 700 m. Therefore, to match the distribution range of the study area, we used the buffer analysis tool in ArcGIS to divide the study area into six levels. The closer an area is to a water system, the higher the possibility of a landslide is. Thus, an area ≤100 m away from a river was coded as 1. The greater the distance is, the lower the probability of a landslide is. Thus, an area >1000 m away from a river was coded as 6.
E.
Geology
Geological map data [50] provided by the Central Geological Survey of the Ministry of Economic Affairs were used to obtain geological names and rock properties corresponding to each grid. The strength grades of various geologies in the area were classified in accordance with the relationship between compressive strength and rock strength grades proposed by the International Society of Rock Mechanics [51]. A higher strength grade was coded as 1, whereas a lower strength grade was coded as 6. The strength grade was divided into six levels [1].
F.
Distance from the fault
Torkashvand et al. [52] reported that the areas within 500 m of a fault have a total landslide area of 39.18% and that the areas further than 3500 m from the fault have high slope stability. In addition, Guo et al. [53] indicated that when the distance from a fault exceeds 2500 m, it does not significantly affect the occurrence of landslides. Therefore, the present study used the buffer analysis tool in ArcGIS to divide the distance from the fault into three levels: less than 500 m, 500–3000 m, and more than 3000 m.

4.2.2. Disturbance Factor of Slope Land Utilization

Slope land use varies with time or space. Therefore, this study interpreted the aforementioned satellite images to determine the slope land use disturbance and obtain information on landslides in the study area for different periods. In accordance with the method reported by Chen et al. [1,43], we selected five items, namely road density, building density, crop density, green coverage rate, and bare land density, as slope land use disturbance factors affecting the occurrence of landslides in the study area. We used the image classification module (segmentation and classification tools) in ArcGIS to calculate the proportions of slope land use disturbance factors in each basic grid (40 m × 40 m) of the study area after image interpretation. The size of the interpretation grid was equal to the resolution of each satellite image.
The method used to estimate the degree of slope land use disturbance under specific natural environmental conditions in the study area was based on the method reported by Chen et al. [43] and is presented in Equation (3).
I D L D = G D C G E C = ( W D C R D C ) ( W E C R E C ) ,
where the grade of the disturbance condition (GDC) is the score of the aforementioned slope land use disturbance in each grid and the grade of the environmental conditions (GEC) is the score of the natural environment of the sloping land in each grid. GDC is the ratio of the area of each slope land use disturbance factor in the basic grid, and REC represents the coding of the natural environmental factors of each slope land in the basic grid.

4.2.3. Rain Trigger Factors

Rainfall data before and after a typhoon or rainstorm occurred in the study area in the past few years were obtained from the rain gauge station [46] of the Central Weather Bureau. The inverse distance weighted (IDW) method of the ArcGIS spatial interpolation calculation tool was used to estimate the distribution of the rainfall trigger indicators in the study area as a whole and in the basic grid.

4.2.4. Earthquake Trigger Factors

The earthquake factor index of landslides that occurred in the study area was determined as the PGA from the free-field strong seismic observation network of the Seismological Center of the Central Weather Bureau [54]. We used the IDW function of ArcGIS Spatial Analyst to estimate the distribution of seismic factors in the study area, and the results were evaluated for the basic grid.

4.3. Weight Analysis of Hazard-Causing Factors

This study used the RF method to analyze the explanatory power of the hazard factors and used the mean decrease accuracy of the RF analytical results to investigate the effects of the characteristic variables and calculate the explanatory power of each hazard factor. To determine the relationship between the slope land use disturbance factors and natural environmental factors, we used the aforementioned set of natural environmental factor classification codes. We used the aforementioned set of indicators to investigate the effect of the correlation between slope land use disturbance and natural environmental factors on landslides [55]. Building density, road density, and bare land density were all positively correlated with landslides, whereas crop density and green coverage were negatively correlated with landslides. Elevation, slope, aspect, and distance coding from a water system were negatively correlated with the occurrence of landslides, whereas distance coding from the fault and geological coding were positively correlated with the occurrence of landslides.
The average reduction in the precision of slope land use disturbance factors and natural environmental factors was divided by the total average reduction in precision to estimate the explanatory power of each hazard factor. After obtaining the square root and adding the correlation coefficient, we obtained positive and negative correlation values. The effects of the characteristic variables of slope land use disturbance factors and natural environmental factors are summarized in Table 2 (a) and Table 2 (b), respectively.
As shown in Table 2 (a), the larger the correlation value was, the stronger the effect of the factor on the occurrence of landslides was. The impact score of bare land density on landslide occurrence was 5. By contrast, the impact score of the green coverage rate with the smallest correlation value was 1. The impact scores of the slope land use disturbance factors on landslide occurrence are listed in Table 3.
The aforementioned findings related to the hierarchical coding of natural environmental factors indicated that the higher the probability of a landslide was, the smaller the coding value was. According to the positive and negative correlation values of the natural environmental factors listed in Table 2 (b), the correlation value of the slope coding was the smallest. Therefore, the score of the impact of the slope coding on landslide occurrence was set to 6. By contrast, the correlation value of the distance from the fault coding was the highest, and the score of its impact on landslide occurrence was set to 1. The impact scores of the related natural factors on landslide occurrence are listed in Table 4.
In this study, the EAR × I3Rmax-normalized data of each single rainfall occurrence and the PGA-normalized data of each single earthquake were integrated with the EAR × I3Rmax- and PGA-normalized data of postrainfall earthquakes and postearthquake rainfall. Then, an RF was used to analyze the degree of interpretation. Input data included the normalized values of EAR × I3Rmax and the PGA of each basic grid under each rainfall occurrence or each earthquake in the study area and the landslide condition of the corresponding basic grid after each rainfall occurrence or each earthquake. The mean reduction in precision was used to determine the explanatory degree of characteristic variables and estimate the explanatory power of different trigger factors. The degree of influence of EAR × I3Rmax and PGA characteristic variables under different trigger factors is summarized in Table 5. The explanatory power of each trigger factor item was set as the score value, and the postrainfall earthquake or postearthquake rainfall score value was obtained by adding up their individual explanatory powers. The effect of rainfall trigger factors on landslide occurrences in the study area was slightly higher than that of earthquake trigger factors.

4.4. Establishment and Discussion of Hazard Indicators of the Interactive Correlation between Rainfall- and Earthquake-Induced Landslides

4.4.1. Establishment of Hazard Indicators of Interaction between Rainfall and Earthquakes

This study explored the effect of the interaction between rainfall and earthquake trigger factors on landslide occurrences in slope land in the study area. This study referred to the degree of land development and utilization proposed by Chen et al. [1] and Chen et al. [56] and adjusted the hazard potential index to enable consideration of slope land disturbance. The IHERI of a basic grid number i after each rainfall occurrence or earthquake in the study area is defined as follows:
[ I H E R I ] i = [ ( C E + C R E ) F ( T E ) + ( C R + C E R ) F ( T R ) ] i
  • F(TR): Standardized values of rainfall factors in each basic grid
  • F(TE): Standardized values of seismic factors in each basic grid
  • CR: Estimated score value for a single rainfall-induced landslide
  • CE: Estimated score value for a single earthquake-induced landslide
  • CRE: Score value calculated for a postrainfall earthquake-induced landslide
  • CER: Score value calculated for a postearthquake rainfall-induced landslide
where F(TR) or F(TE) can be obtained from the spatial distribution data of the trigger factors (EAR × I3Rmax or PGA) of each rainfall occurrence or each earthquake in the study area. The values of CR, CE, CRE, and CER can be calculated separately using the aforementioned RF algorithm through a weight analysis. The trigger factors of single rainfall, single earthquake, earthquake after a previous rainfall event, and rainfall after a previous earthquake event induced landslides in the slope land in the study area. The larger the IHERI was, the higher the possibility of the interaction between rainfall and earthquakes inducing a landslide was.
To investigate the relationship between the grid ratio of the landslides (the number of basic grids with landslides divided by the number of basic grids without landslides) induced by the aforementioned rainfall or earthquake events (13 in total) in the study area and IHERI, we first used SPSS Cluster Analysis [55] to categorize the IHERI values into five levels and then plotted the graph for the IHERI against the landslide grid ratio (Figure 2). The landslide grid ratio in the study area tended to increase with an increase in the degree of interaction between rainfall and earthquake events.

4.4.2. Classification of Slope Land Use Disturbance Degree and Its Influence on Landslides in the Study Area

The data collected after 13 rainfall- or earthquake-induced landslides in the study area were merged, and data for 2,339,322 basic grids were obtained. Equation (3) was used to estimate the IDLD of each basic grid. By using SPSS Cluster Analysis, we categorized the IDLD values into five grades. We plotted the IDLD values against the 13 rainfall- or earthquake-induced landslide grid ratios in the study area (the number of basic grids with landslides divided by the number of basic grids where landslides occur) to determine the degree of slope land use disturbance (Figure 3). The landslide grid ratio of the slope land with a low degree of disturbance was 0.006, whereas that of the slope land with a high degree of disturbance was as high as 0.93. The results indicated that the landslide grid ratio of the slope land in the study area increased with the degree of slope land use disturbance.

4.4.3. Interaction between the Hazard Indexes of Rainfall- and Earthquake-Induced Landslides and the Index of the Degree of Slope Land Use Disturbance

This study investigated the relationships among IHERI, IDLD, and landslide occurrence in each basic grid in the study area. A schematic of the interval settings is presented in Figure 4 and the schematic obtained after dividing the study area into 9 equidistant interval grids (from letter A to I) is presented in Figure 5. The red ○’s and gray ⨯’s in the figure represent basic grids with and without landslides, respectively. The basic grid with landslides represents the total number of landslides induced by all 13 rainfall or earthquake events in the study area, which yielded a total of 37,197 records. The grid points in the figure where a landslide did not occur were determined using a 1:1 ratio, which was obtained through random sampling from all basic grids where a landslide did not occur.
The number of basic grids and the ratio of collapsed grids to noncollapsed grids for each interval are presented in Figure 5 and summarized in Table 6. In each interval, with the exception of I, which did not have data, only the A and B grids without landslides were larger than the grids with landslides. The ratio of landslide grids (the ratio of the number of landslide grids to the number of nonlandslide grids or the ratio of the number of landslide grids to the total number of grids in an interval) sequentially increased from A to H. When the IHERI or IDLD was larger, the induced landslide ratio increased. When the degree of slope land use disturbance was similar in the study area, the landslide ratio, which indicated the probability of a landslide, increased with the IHERI value. Similarly, when the IHERI was similar in the study area, the probability of a landslide increased with the degree of slope land use disturbance. The interaction between rainfall and earthquake events in the study area indicates that the effect of the degree of slope land use disturbance on landslide occurrence is notable.

4.5. Establishment and Verification of Rainfall- and Earthquake-Induced Landslide Potential Assessment Models and Potential Map Drawing

4.5.1. Establishment of a Landslide Potential Assessment Model

On the basis of the data collected on land surface changes and disasters after 13 rainfall or earthquake events in the study area between 2004 and 2018 and with reference to Wang [39], we used the RF algorithm developed in R language to establish a landslide potential assessment model with consideration of four trigger factors. After the collection and analysis of 839,047 records of single rainfall events (Typhoon Mindulle, 0726 heavy rain, Nanmadu typhoon, 0520 heavy rain, and Megi typhoon), 1,033,319 records of single earthquake events (0226 earthquake, 0602 earthquake, 0206 earthquake, 1122 earthquake, and 0320 earthquake), 163,189 records of postrainfall earthquake events (1105 earthquake after Morakot and 1108 earthquake after Fanapi), and 303,767 records of postearthquake rainfall events (Kalmaegi after 0305 earthquake), with reference to Chen et al. [56], we randomly sampled the records without landslides on the basis of the number of landslides and then constructed training and testing data sets at a ratio of 7:3. By performing an error matrix analysis, we determined the PA and OA and used them as evaluation factors for establishing precision criteria. The results are summarized in Table 7.
The PA of both the training and testing sets was above 83%, and the average OA was 95.7. The finding indicated that the evaluation model could achieve a high degree of accuracy for landslide classification. In addition, both the PA and OA of potential landslides induced by a single rainfall event, a single earthquake, a postrainfall earthquake, and postearthquake rainfall were above 83%, indicating excellent accuracy (Table 8, Table 9, Table 10 and Table 11).

4.5.2. Verification of the Landslide Potential Assessment Model and Drawing of a Potential Map

We used the ROC to determine the effectiveness of applying the established landslide potential assessment model in the study area. We employed the AUC as a reference standard. The AUC values of the results for landslides induced by a single rainfall event, a single earthquake, a postrainfall earthquake, and postearthquake rainfall in the study area were 0.95, 0.99, 0.94, and 0.97, respectively. The prediction for induced landslides was highly accurate. We used the equal interval classification method of ArcGIS to classify the values for potential landslides induced by rainfall or earthquakes into four categories: 0–0.25, 0.25–0.5, 0.5–0.75, and 0.75–1. These categories correspond to low potential, medium-low potential, medium-high potential, and high potential, respectively. The higher the potential value is, the greater the probability of a landslide occurrence is. Figure 6, Figure 7, Figure 8 and Figure 9 present the potential maps of the study area after the occurrence of landslides induced by a single rainfall event, a single earthquake, a postrainfall earthquake, and postearthquake rainfall.
This study referred to information provided by the Central Geological Survey [50] of the Ministry of Economic Affairs and the Soil and Water Conservation Bureau [57] of the Executive Yuan to obtain a distribution map of the historical disaster locations in the study area and an overlay map of potential landslides induced by the aforementioned rainfall and earthquake events. The landslide potential maps of Alishan Township and Dapu Township in the study area overlapped with historical landslide range results, as illustrated in Figure 10 and Figure 11, respectively.
In this study, we calculated the landslide area ratios of medium–high-level landslides induced by rainfall and earthquakes in Alishan Township and Dapu Township, respectively, to determine the degree of damage caused by different trigger factors. The area ratio of medium–high-level landslides was calculated as the area of historical disaster regions with landslide potential values of >0.5 divided by the total historical disaster area. To prevent cloud cover or shadows from affecting the results, we combined image grids affected by cloud cover or shadows before and after the 13 rainfall or earthquake events in the study area. Thus, the total grid number of satellite images before and after the rainfall and earthquake events was identical. The area ratios of medium–high-level landslides induced by rainfall and earthquakes in Alishan Township are shown in Table 12 and Figure 12, respectively, and those of Dapu Township are presented in Table 13 and Figure 13, respectively. With the exception of the area ratios of landsides induced by Typhoon Mindulle in 2004 in Alishan Township and 0520 heavy rainfall in Dapu Township, the landslide areas affected by the interaction between rainfall and earthquakes were large. A larger area ratio of medium–high-level landslides was observed in Alishan Township after typhoon Minduli in 2004 and Dapu Township after 0520 heavy rainfall. This may be because the average amount of rainfall was higher than that during trigger events (Table 14), which caused the probability of landslides during rainfall events to be higher than that during the dual trigger events. In addition, on the basis of the results presented in Table 12, Table 13 and Table 14, we plotted the trigger index values (EAR × I3Rmax or PGA) against the landslide area ratios of medium–high-level landslides induced by a single rainfall or a single earthquake event in the study area (Figure 14 and Figure 15, respectively). Irrespective of whether the average value of the EAR × I3Rmax of a single rainfall event was larger or the maximum value of the PGA of a single earthquake was larger, the landslide area ratio of the medium–high-level landslides exhibited an increasing trend. This finding indicates that the landslide potential estimated in this study is reasonable.

5. Conclusions

This study used the random tree image classification module of a GIS platform to interpret the satellite images of the study area before and after the occurrence of landslides induced by four types of events, namely a single rainfall event, a single earthquake, a postrainfall earthquake, and postearthquake rainfall. We used the texture analysis module of ERDAS IMAGINE to perform image interpretation and obtain information on surface changes and landslide areas in the study area. The major mapped landslides are debris slides. The average Kappa value (0.66) and average OA (70.1%) of the 37 satellite images before and after the 13 rainfall and earthquake events in the study area, which were determined using seven interpretive classifications, namely water, roads, buildings, cash crops, green coverage, river channels, and bare land, reveal that the interpretation accuracy was medium to medium high. In addition, we established four landslide potential assessment models for landslides induced by different trigger factors by using an RF. We established a hazard index (IHERI) for landslides induced by postrainfall earthquakes and postearthquake rainfall and estimated the IDLD. This study then explored the characteristics of IHERI that affect landslides under specific natural environmental factor and slope land use disturbance conditions. The weight analysis results indicate that of the slope land use disturbance factors, bare land density exerted the strongest effect. Of the natural environmental factors, slope coding exerted the strongest effect. Compared with earthquakes, rainfall exerted a stronger effect. Irrespective of whether landslides were induced by a single rainfall event, a single earthquake, or both rainfall and earthquakes, the grid ratio for landslides in the study area tended to increase with an increase in the degree of slope land use disturbance. When IHERI or IDLD was larger, the induced landslide ratio increased. Thus, when the interaction between rainfall and earthquakes is considered, the impact of slope land use disturbance on landslide hazards should be noted. Moreover, regardless of whether single rainfall, single earthquake, or both rainfall and earthquake events were considered, the accuracy of each landslide potential model constructed using the RF method was above 83%, which is considered excellent. The AUC values in the assessment results for landslides induced by single rainfall, single earthquake, postrainfall earthquake, and postearthquake rainfall events were all above 0.94, indicating a higher risk of landslide induction by single rainfall and earthquake events. The forecasts had a high accuracy. With the exception of when the average rainfall in the study area was higher than that during the dual trigger events, the area ratio of medium–high-level landslides induced by the interaction between rainfall and earthquakes was large. According to the estimated landslide potential in this study, irrespective of whether the average EAR × I3Rmax value of a single rainfall event is larger or the maximum PGA value of a single earthquake is larger, the area of medium–high-level landslides tends to increase.

Author Contributions

Conceptualization, C.-M.T., Y.-R.C. (Yie-Ruey Chen), C.-M.C., Y.-L.Y., Y.-R.C. (Yu-Ru Chen) and S.-C.H.; methodology, C.-M.T., Y.-R.C. (Yie-Ruey Chen) and C.-M.C.; investigation Y.-R.C. (Yu-Ru Chen) and S.-C.H.; writing—original draft preparation, C.-M.T. and Y.-R.C.; writing—review, C.-M.T., Y.-R.C. (Yie-Ruey Chen), C.-M.C., Y.-L.Y., Y.-R.C. (Yu-Ru Chen) and S.-C.H.; editing, C.-M.T. and Y.-R.C. (Yie-Ruey Chen); supervision, C.-M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by grants from the Taiwan National Science and Technology Council (NSTC 110-2625-M-309-001 and NSTC 111-2625-M-309-001). The authors wish to express their appreciation to Disaster Prevention Research Center, NCKU, for providing the technical support of ENVI.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, Y.R.; Chen, J.W.; Hsieh, S.C.; Ni, P.N. The Application of Remote Sensing Technology to the Interpretation of Land Use for Rainfall-Induced Landslides Based on Genetic Algorithms and Artificial Neural Networks. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2009, 2, 87–95. [Google Scholar] [CrossRef]
  2. Chue, Y.S.; Chen, J.W.; Chen, Y.R. Rainfall-Induced Slope Landslide Potential and Landslide Distribution Characteristics Assessment. J. Mar. Sci. Technol. 2015, 23, 705–716. [Google Scholar] [CrossRef]
  3. Chen, Y.C.; Chang, K.T.; Lee, H.Y.; Chiang, S.H. Average landslide erosion rate at the watershed scale in southern Taiwan estimated from magnitude and frequency of rainfall. Geomorphology 2015, 228, 756–764. [Google Scholar] [CrossRef]
  4. Feng, Q.; Liu, J.; Gong, J. UAV Remote Sensing for Urban Vegetation Mapping Using Random Forest and Texture Analysis. Remote Sens. 2015, 7, 1074–1094. [Google Scholar] [CrossRef] [Green Version]
  5. Stumpf, A.; Malet, J.P.; Delacourt, C. Correlation of satellite image time-series for the detection and monitoring of slow-moving landslides. Remote Sens. Environ. 2017, 189, 40–55. [Google Scholar] [CrossRef]
  6. Zhong, Y.Y. Comparison on the Characteristics of Rainfall-Induced Landslides before and after the Chi-Chi Earthquake. Master’s Thesis, National Taiwan University, Taipei, Taiwan, 2005. (In Chinese). [Google Scholar] [CrossRef]
  7. Lin, C.W.; Liu, S.H.; Lee, S.Y.; Liu, C.C. Impacts of the Chi-Chi earthquake on subsequent rainfall-induced landslides in central Taiwan. Eng. Geol. 2006, 86, 87–101. [Google Scholar] [CrossRef]
  8. Tang, C.; Zhu, J.; Qi, X.; Ding, J. Landslides induced by the Wenchuan earthquake and the subsequent strong rainfall event: A case study in the Beichuan area of China. Eng. Geol. 2011, 122, 22–23. [Google Scholar] [CrossRef]
  9. Huang, H.P.; Yang, K.C.; Lin, B.W. Statistical evaluation of the effect of earthquake with other related factors on landslide susceptibility: Using the watershed area of Shihmen reservoir in Taiwan as a case study. Environ. Earth Sci. 2013, 69, 2151–2166. [Google Scholar] [CrossRef]
  10. Chen, Y.R.; Tsai, K.J.; Hsieh, S.C.; Ho, Y.L. Evaluation of Landslide Potential due to Land Use in the Slope. Electron. J. Geotech. Eng. 2015, 20, 4277–4292. [Google Scholar]
  11. Yang, Z.H.; Lan, H.X.; Geo, X.; Li, L.P.; Meng, Y.S.; Wu, Y.M. Urgent landslide susceptibility assessment in the 2013 Lushan earthquake-impacted area, Sichuan Province, China. Nat. Hazards 2015, 75, 2467–2487. [Google Scholar] [CrossRef]
  12. Jan, C.D.; Yang, S.Y.; Su, Y.W.; Haung, W.S. Investigation about rainfall-induced shallow landslides in CYL and TWR watersheds, Taiwan. Environ. Earth Sci. 2016, 75, 898. [Google Scholar] [CrossRef]
  13. Li, C.; Wang, M.; Liu, K. A decadal evolution of landslides and debris flows after the Wenchuan earthquake. Geomorphology 2018, 323, 1–12. [Google Scholar] [CrossRef]
  14. Tseng, C.M.; Chen, Y.R.; Wu, S.M. Scale and spatial distribution assessment of rainfall-induced landslides in a catchment with mountain roads. Nat. Hazards Earth Syst. Sci. 2018, 18, 687–708. [Google Scholar] [CrossRef] [Green Version]
  15. Wistuba, M.; Malik, I.; Krzemień, K.; Gorczyca, E.; Sobucki, M.; Wrońska-Wałach, D.; Gawior, D. Can low-magnitude earthquakes act as a triggering factor for landslide activity? Examples from the Western Carpathian Mts, Poland. CATENA 2018, 171, 359–375. [Google Scholar] [CrossRef]
  16. Chen, C.W.; Tung, Y.S.; Liou, J.J.; Li, H.C.; Cheng, C.T.; Chen, Y.M.; Oguchi, T. Assessing landslide characteristics in a changing climate in northern Taiwan. CATENA 2019, 175, 263–277. [Google Scholar] [CrossRef]
  17. Valagussa, A.; Marc, O.; Frattini, P.; Crosta, G.B. Seismic and geological controls on earthquake-induced landslide size. Earth Planet. Sci. Lett. 2019, 506, 268–281. [Google Scholar] [CrossRef]
  18. Quesada-Román, A.; Fallas-López, B.; Hernández-Espinoza, K.; Stoffel, M.; Ballesteros-Cánovas, J.A. Relationships between earthquakes, hurricanes, and landslides in Costa Rica. Landslides 2019, 16, 1539–1550. [Google Scholar] [CrossRef]
  19. Ruggeri, P.; Fruzzetti, V.M.E.; Ferretti, A.; Scarpelli, G. Seismic and Rainfall Induced Displacements of an Existing Landslide: Findings from the Continuous Monitoring. Geosciences 2020, 10, 90. [Google Scholar] [CrossRef] [Green Version]
  20. Bontemps, N.; Lacroix, P.; Larose, E.; Jara, J.; Taipe, E. Rain and small earthquakes maintain a slow-moving landslide in a persistent critical state. Nat. Commun. 2020, 11, 780. [Google Scholar] [CrossRef] [Green Version]
  21. Liu, S.; Hu, K.; Zhang, Q.; Zhang, S.; Hu, X.; Tang, D. Quantitative Analysis of the Effects of an Earthquake on Rainfall Thresholds for Triggering Debris-Flow Events. Front. Earth Sci. 2021, 9, 417. [Google Scholar] [CrossRef]
  22. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  23. Colditz, R.R. An Evaluation of Different Training Sample Allocation Schemes for Discrete and Continuous Land Cover Classification Using Decision Tree-Based Algorithms. Remote Sens. 2015, 7, 9655–9681. [Google Scholar] [CrossRef] [Green Version]
  24. Belgiu, M.; Drăguţ, L. Random forest in remote sensing: A review of applications and future directions. ISPRS J. Photogramm. Remote Sens. 2016, 114, 24–31. [Google Scholar] [CrossRef]
  25. Lagomarsino, D.; Tofani, V.; Segoni, S.; Catani, F.; Casagli. A Tool for Classification and Regression Using Random Forest Methodology: Applications to Landslide Susceptibility Mapping and Soil Thickness Modeling. Environ. Model. Assess. 2017, 22, 201–214. [Google Scholar] [CrossRef]
  26. Taalab, K.; Cheng, T.; Zhang, Y. Mapping landslide susceptibility and types using Random Forest. Big Earth Data 2018, 2, 159–178. [Google Scholar] [CrossRef]
  27. Catani, F.; Lagomarsino, D.; Segoni, S.; Tofani, V. Landslide susceptibility estimation by random forests technique: Sensitivity and scaling issues. Nat. Hazards Earth Syst. Sci. 2013, 13, 2815–2831. [Google Scholar] [CrossRef] [Green Version]
  28. Goetz, J.N.; Brenning, A.; Petschko, H.; Leopold, P. Evaluating machine learning and statistical prediction techniques for landslide susceptibility modelling. Comput. Geosci. 2015, 81, 1–11. [Google Scholar] [CrossRef]
  29. Steger, S.; Brenning, A.; Bell, R.; Petschko, H.; Glade, T. Exploring discrepancies between quantitative validation results and the geomorphic plausibility of statistical landslide susceptibility maps. Geomorphology 2016, 262, 8–23. [Google Scholar] [CrossRef]
  30. Chen, W.; Li, Y.; Tsangaratos, P.; Shahabi, H.; Ilia, I.; Xue, W. Groundwater spring potential mapping using artificial intelligence approach based on kernel logistic regression, random forest, and alternating decision tree models. Appl. Sci. 2020, 10, 425. [Google Scholar] [CrossRef] [Green Version]
  31. Ercanoglu, M. Landslide susceptibility assessment of SE Bartin (West Black Searegion, Turkey) by artificial neural networks. Nat. Hazards Earth Syst. Sci. 2005, 5, 979–992. [Google Scholar] [CrossRef]
  32. Shahabi, H.; Khezri, S.; Ahmad, B.B.; Hashim, M. Landslide susceptibility mapping at central Zab basin, Iran: A comparison between analytical hierarchy process frequency ratio and logistic regression models. Catena 2014, 115, 55–70. [Google Scholar] [CrossRef]
  33. Ho, T.K. The random subspace method for constructing decision forests. IEEE Trans. Pattern Anal. Mach. Intell. 1998, 20, 832–844. [Google Scholar] [CrossRef] [Green Version]
  34. Efron, B. Computers and the Theory of Statistics: Thinking the Unthinkable. SIAM Rev. 1979, 21, 460–480. [Google Scholar] [CrossRef]
  35. Cutler, A.; Cutler, D.R.; Stevens, J.R. Random Forests. In Ensemble Machine Learning; Zhang, C., Ma, Y., Eds.; Springer: Boston, MA, USA, 2012; pp. 157–175. [Google Scholar]
  36. Zhang, H.; Wu, P.B.; Yin, A.J.; Yang, X.H.; Zhang, M.; Gao, C. Prediction of soil organic carbon in an intensively managed reclamation zone of eastern China: A comparison of multiple linear regressions and the random forest model. Sci. Total Environ. 2017, 592, 704–713. [Google Scholar] [CrossRef]
  37. Youssef, A.M.; Pourghasemi, H.R.; Pourtaghi, Z.S.; Al-Katheeri, M.M. Landslide susceptibility mapping using random forest, boosted regression tree, classification and regression tree, and general linear models and comparison of their performance at Wadi Tayyah Basin, Asir Region, Saudi Arabia. Landslides 2016, 13, 839–856. [Google Scholar] [CrossRef]
  38. Chen, W.; Peng, J.B.; Hong, H.Y.; Shahabi, H.; Pradhan, B.; Liu, J.Z.; Zhu, A.; Pie, X.; Duan, Z. Landslide susceptibility modelling using GIS-based machine learning techniques for Chongren County, Jiangxi Province, China. Sci. Total Environ. 2018, 626, 1121–1135. [Google Scholar] [CrossRef]
  39. Wang, Y.L. The Assessment of Rainfall-induced Landslides Using Data Mining Techniques. Master’s Thesis, Chang Jung Christian University, Tainan, Taiwan, 2018. [Google Scholar]
  40. Sklansky, J. Image Segmentation and Feature Extraction. IEEE Trans. Syst. Man Cybern. 1978, 8, 238–247. [Google Scholar] [CrossRef]
  41. ERDAS. ERDAS IMAGINE® Tour Guide; ERDAS World Headquarter: Atlanta, GA, USA, 2011. [Google Scholar]
  42. Verbyla, D.L. A gentic adaptive neural network approach to pricing options: A simulation analysis. J. Comput. Intell. Financ. 1995, 6, 13–23. [Google Scholar]
  43. Chen, Y.R.; Ni, P.N.; Tsai, K.J. Construction of a sediment disaster risk assessment model. Environ. Earth Sci. 2013, 70, 115–129. [Google Scholar] [CrossRef]
  44. Landis, J.R.; Koch, G.G. The measurement of observer agreement for categorical data. Biometrics 1977, 33, 159–174. [Google Scholar] [CrossRef] [Green Version]
  45. Fawcett, T. An introduction to ROC analysis. Pattern Recognit. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  46. Central Weather Bureau. 2021. Available online: https://www.cwb.gov.tw (accessed on 1 January 2021).
  47. Water Resources Agency, Ministry of Economic Affairs. 2021. Available online: https://www.wra.gov.tw/ (accessed on 1 January 2021).
  48. Moayedi, H.; Khari, M.; Bahiraei, M.; Foong, L.K.; Bui, D.T. Spatial assessment of landslide risk using two novel integrations of neuro-fuzzy system and metaheuristic approaches; Ardabil Province, Iran. Geomat. Nat. Hazards Risk 2020, 11, 230–258. [Google Scholar] [CrossRef] [Green Version]
  49. Hong, H.; Pradhan, B.; Jebur, M.N.; Bui, D.T.; Xu, C.; Akgun, A. Spatial prediction of landslide hazard at the luxi area (China) using support vector machines. Environ. Earth Sci. 2016, 75, 40. [Google Scholar] [CrossRef]
  50. Central Geological Survey, Ministry of Economic Affairs. 2021. Available online: https://www.moeacgs.gov.tw/ (accessed on 1 January 2021).
  51. ISRM. Rock Characterization, Testing and Monitoring; Pergamon Press: Oxford, UK, 1981. [Google Scholar]
  52. Torkashvand, A.M.; Irani, A.; Sorur, J. The preparation of landslide map by landslide numerical risk factor (LNRF) model and geographic information system (GIS). Egypt. J. Remote Sens. Space Sci. 2014, 17, 159–170. [Google Scholar] [CrossRef] [Green Version]
  53. Guo, C.; Montgomery, D.R.; Zhang, Y.; Wang, K.; Yang, Z. Quantitative assessment of landslide susceptibility along the xianshuihe fault zone, tibetan plateau, China. Geomorphology 2015, 248, 93–110. [Google Scholar] [CrossRef]
  54. Seismological Center, Central Weather Bureau. 2021. Available online: https://scweb.cwb.gov.tw/zh-tw/earthquake/data (accessed on 1 January 2021).
  55. SPSS. SPSS 14.0 Brief Guide; SPSS Inc.: Chicago, IL, USA, 2005. [Google Scholar]
  56. Chen, Y.R.; Chang, C.M.; Tseng, C.M.; Hung, M.Z.; Hsieh, S.C. New model of hazard potential index based on slope land disturbance and regional carrying capacity: A case study of Namaxia District in Kaohsiung City, Taiwan. Environ. Earth Sci. 2021, 80, 440. [Google Scholar] [CrossRef]
  57. Soil and Water Conservation Bureau, Council of Agriculture, Executive Yuan. 2021. Available online: https://246.swcb.gov.tw/ (accessed on 1 January 2021).
Figure 1. Study areas.
Figure 1. Study areas.
Water 14 03691 g001
Figure 2. The hazard index IHERI against landslide grid ratio.
Figure 2. The hazard index IHERI against landslide grid ratio.
Water 14 03691 g002
Figure 3. The degree of land disturbance IDLD against landslide grid ratio.
Figure 3. The degree of land disturbance IDLD against landslide grid ratio.
Water 14 03691 g003
Figure 4. Schematic of the interval setting method of the relationships among the IHERI, IDLD, and landslide occurrence of each basic grid.
Figure 4. Schematic of the interval setting method of the relationships among the IHERI, IDLD, and landslide occurrence of each basic grid.
Water 14 03691 g004
Figure 5. Relationships among the IHERI, IDLD, and landslide occurrence after various rainfall and earthquake events in the study area.
Figure 5. Relationships among the IHERI, IDLD, and landslide occurrence after various rainfall and earthquake events in the study area.
Water 14 03691 g005
Figure 6. Single rainfall event-induced landslide potential.
Figure 6. Single rainfall event-induced landslide potential.
Water 14 03691 g006
Figure 7. Single earthquake-induced landslide potential.
Figure 7. Single earthquake-induced landslide potential.
Water 14 03691 g007
Figure 8. Postrainfall earthquake-induced landslide potential.
Figure 8. Postrainfall earthquake-induced landslide potential.
Water 14 03691 g008
Figure 9. Postearthquake rainfall-induced landslide potential.
Figure 9. Postearthquake rainfall-induced landslide potential.
Water 14 03691 g009
Figure 10. Landslide potential map of Alishan Township after rainfall and earthquakes overlapped by historical landslides.
Figure 10. Landslide potential map of Alishan Township after rainfall and earthquakes overlapped by historical landslides.
Water 14 03691 g010
Figure 11. Landslide potential map of Dapu Township after rainfall or earthquakes overlapped by historical landslides.
Figure 11. Landslide potential map of Dapu Township after rainfall or earthquakes overlapped by historical landslides.
Water 14 03691 g011aWater 14 03691 g011b
Figure 12. Landslide area ratio of medium–high-level landslides induced by rainfall and earthquakes in Alishan Township.
Figure 12. Landslide area ratio of medium–high-level landslides induced by rainfall and earthquakes in Alishan Township.
Water 14 03691 g012
Figure 13. Landslide area ratio of medium–high-level landslides induced by rainfall and earthquakes in Dapu Township.
Figure 13. Landslide area ratio of medium–high-level landslides induced by rainfall and earthquakes in Dapu Township.
Water 14 03691 g013
Figure 14. Landslide area ratios of medium–high-level landslides against average EAR × I3Rmax and maximum PGA induced by a single rainfall or a single earthquake event in Alishan Township.
Figure 14. Landslide area ratios of medium–high-level landslides against average EAR × I3Rmax and maximum PGA induced by a single rainfall or a single earthquake event in Alishan Township.
Water 14 03691 g014
Figure 15. Landslide area ratios of medium–high-level landslides against average EAR × I3Rmax and maximum PGA induced by a single rainfall or a single earthquake event in Dapu Township.
Figure 15. Landslide area ratios of medium–high-level landslides against average EAR × I3Rmax and maximum PGA induced by a single rainfall or a single earthquake event in Dapu Township.
Water 14 03691 g015
Table 1. Coefficient of agreement and overall accuracy of satellite image interpretation results.
Table 1. Coefficient of agreement and overall accuracy of satellite image interpretation results.
Satellite ImageKappaOA (%)Kappa
(Mosaic)
OA (%)
(Mosaic)
NoYearDateBefore/After EventTownship
1200410 Feb.Before Typhoon MindulleAlishan0.6065.7
2200410 JulyAfter Typhoon MindulleAlishan0.7577.1
320085 Jan.Before the 0305 Earthquake and Typhoon KalmaegiAlishan0.6568.00.6567.5
4200810 Jan.Before the 0305 Earthquake and Typhoon KalmaegiAlishan0.6466.9
5200821 JulyAfter the 0305 Earthquake and Typhoon KalmaegiAlishan0.7577.1
620085 Jan.Before the 0305 Earthquake and Typhoon KalmaegiDapu0.6368.0
7200821 JulyAfter the 0305 Earthquake and Typhoon KalmaegiDapu0.7073.1
8200912 Apr.Before Typhoon Morakot and the 1105 earthquakeDapu0.6670.3
920096 Nov.After Typhoon Morakot and the 1105 earthquakeDapu0.6266.9
10201011 Apr.Before the 0726 heavy rainDapu0.6770.9
1120104 Aug.After the 0726 heavy rain Dapu0.6771.4
1220104 Aug.Before Typhoon Fanapi and the 1108 EarthquakeDapu0.6569.1
13201027 Dec.After Typhoon Fanapi and the 1108 EarthquakeDapu0.6569.1
14201127 JulyBefore the Typhoon NanmadolDapu0.7174.90.6872.0
15201116 Aug.Before the Typhoon NanmadolDapu0.6469.1
16201127 Sep.After the Typhoon NanmadolDapu0.6166.9
17201210 Feb.Before the 0226 earthquakeDapu0.7880.0
1820127 Mar.After the 0226 earthquakeDapu0.6670.2
1920132 JuneBefore the 0602 earthquakeAlishan0.6871.4
20201329 JuneAfter the 0602 earthquakeAlishan0.6265.70.6366.3
2120134 JulyAfter the 0602 earthquakeAlishan0.6366.9
22201528 Feb.Before the 0520 heavy rainAlishan0.6668.6
23201510 JuneAfter the 0520 heavy rainAlishan0.6870.9
2420155 Mar.Before the 0520 heavy rainDapu0.7374.9
25201510 JuneAfter the 0520 heavy rainDapu0.7073.1
26201528 Nov.Before the 0206 earthquakeAlishan0.6165.7
27201614 Feb.After the 0206 earthquakeAlishan0.6466.90.6870.3
28201630 Mar.After the 0206 earthquakeAlishan0.7273.7
29201528 Nov.Before the 0206 earthquakeDapu0.6368.0
30201630 Mar.After the 0206 earthquakeDapu0.6670.3
31201630 Mar.Before Typhoon MegiDapu0.6669.7
32201619 Nov.After Typhoon MegiDapu0.6468.6
33201718 Oct.Before the 1122 earthquakeAlishan0.6064.60.6266.3
34201717 Nov.Before the 1122 earthquakeAlishan0.6468.0
35201816 Jan.After the 1122 earthquakeAlishan0.6669.7
36201717 Nov.Before the 0320 earthquakeDapu0.6568.6
3720189 Apr.After the 0320 earthquakeDapu0.7274.3
Total AverageKappa = 0.66    OA = 70.1%
Table 2. Effects of the characteristic variables of slope land use disturbance factors and natural environmental factors.
Table 2. Effects of the characteristic variables of slope land use disturbance factors and natural environmental factors.
(a) Slope Land Use Disturbance Factor
ItemDegree of Influence of Characteristic Variables
Mean Decrease AccuracyExplanatory PowerCorrelation Value
Road Density39.880.150.39
Building Density37.590.140.37
Bare Density48.000.180.42
Crop Density61.420.22−0.48
Green Coverage84.350.31−0.56
Total271.241
(b) Natural Environment Factors
ItemDegree of Influence of Characteristic Variables
Mean Decrease AccuracyExplanatory PowerCorrelation Value
Elevation Code103.090.13−0.36
Slope Code235.930.29−0.54
Aspect Code90.900.11−0.33
Geology Code117.760.140.37
Distance Code from Water system150.720.18−0.42
Distance Code from Fault120.810.150.39
Total819.211
Table 3. Impact scores of slope land use disturbance factors on landslides.
Table 3. Impact scores of slope land use disturbance factors on landslides.
Slope Use Disturbance FactorGreen CoverageCrop DensityBuilding DensityRoad DensityBare Density
Score12345
Table 4. Impact scores of natural environmental factors on landslide occurrence.
Table 4. Impact scores of natural environmental factors on landslide occurrence.
Natural Environment FactorsDistance Code from FaultGeology CodeAspect CodeElevation CodeDistance Code from Water SystemSlope Code
Score123456
Table 5. Degree of influence of EAR × I3Rmax and PGA under different trigger factors.
Table 5. Degree of influence of EAR × I3Rmax and PGA under different trigger factors.
Trigger FactorIndexInfluence Degree of Characteristic Variables
Mean Decrease AccuracyExplanatory PowerScore
Single RainEAR × I3Rmax90.970.270.27
Single EarthquakePGA54.280.160.16
Post-Earthquake RainfallEAR × I3Rmax48.340.140.32
PGA60.880.18
post-rainfall EarthquakeEAR × I3Rmax45.500.130.25
PGA39.340.12
Total339.311
Table 6. The IHERI, IDLD, and number of grids with and without landslides and the ratio of the landslide grids in each interval corresponding to the occurrence of landslides.
Table 6. The IHERI, IDLD, and number of grids with and without landslides and the ratio of the landslide grids in each interval corresponding to the occurrence of landslides.
Interval No.Number of Grids with LandslideNumber of Grids without LandslideLandslide Grid Ratio
Number of Grids with Landslide/
Number of Grids without Landslide
Number of Grids with Landslide/
Total Number of Grids in the Interval
A16,37333,3260.490.33
B165622970.720.42
C6444891.320.57
D14,795100114.800.94
E11324724.100.96
F216636.000.97
G22053073.500.99
H1761176.000.99
I00
Table 7. Producer accuracy and overall accuracy for training and testing data sets.
Table 7. Producer accuracy and overall accuracy for training and testing data sets.
Trigger Factor(a) Single Rainfall(b) Single Earthquake
Accuracy TrainingTestingTrainingTesting
PA of Grids with Landslide91.61%83.14%98.66%95.96%
PA of Grids without Landslide99.40%92.17%99.90%95.63%
Overall Accuracy100%89.84%99.92%95.70%
Trigger Factor(c) Post-rainfall Earthquake(d) Post-earthquake Rainfall
Accuracy TrainingTestingTrainingTesting
PA of Grids with Landslide98.23%87.02%96.26%89.99%
PA of Grids without Landslide99.68%88.01%98.21%92.22%
Overall Accuracy100%88.27%99.84%91.67%
Table 8. Evaluation accuracy of single rainfall-induced landslide potential.
Table 8. Evaluation accuracy of single rainfall-induced landslide potential.
Single Typhoon or RainfallTyphoon Mindulle0726 Heavy RainTyphoon Nanmadol
CoverageAlishan TownshipDapu TownshipDapu Township
Evaluation Result
LandslideNon-landslideLandslideNon-landslideLandslideNon-landslide
Actual SituationLandslide608391969725221
Non-landslide20,948225,031399377,505227976,647
PA of Grids with Landslide86.88%99.71%99.81%
PA of Grids without Landslide91.48%95.10%97.11%
Overall Accuracy91.36%95.14%97.13%
Single Rainfall or Earthquake0520 Heavy Rain0520 Heavy RainTyphoon Megi
CoverageAlishan TownshipDapu TownshipDapu Township
Evaluation Result
LandslideNon-landslideLandslideNon-landslideLandslideNon-landslide
ActualLandslide3958467101217757738
Non-landslide18,868242,051534170,094329678,541
PA of Grids with Landslide89.45%85.11%93.82%
PA of Grids without Landslide92.77%92.92%95.96%
Overall Accuracy92.71%92.80%95.96%
Table 9. Evaluation accuracy of single earthquake-induced landslide potential.
Table 9. Evaluation accuracy of single earthquake-induced landslide potential.
Single Earthquake0226 Earthquake0602 Earthquake0206 Earthquake
CoverageDapu TownshipAlishan TownshipAlishan Township
Evaluation Result
LandslideNon-landslideLandslideNon-landslideLandslideNon-landslide
Actual SituationLandslide6872266021573839122
Non-landslide315576,99010,587236,64511,263249,499
PA of Grids with Landslide96.90%97.68%96.92%
PA of Grids without Landslide96.06%95.72%95.68%
Overall Accuracy96.07%95.77%95.70%
Single Earthquake0206 Earthquake1122 Earthquake0320 Earthquake
CoverageDapu TownshipAlishan TownshipDapu Township
Evaluation Result
LandslideNon-landslideLandslideNon-landslideLandslideNon-landslide
Actual SituationLandslide47527313964845
Non-landslide219779,73114,714253,350321976,404
PA of Grids with Landslide94.62%99.81%98.98%
PA of Grids without Landslide97.32%94.51%95.96%
Overall Accuracy97.30%94.57%95.98%
Table 10. Evaluation accuracy of postrainfall earthquake-induced landslide potential.
Table 10. Evaluation accuracy of postrainfall earthquake-induced landslide potential.
Post-Rainfall Earthquake1105 Earthquake after Typhoon Morakot1108 Earthquake after Typhoon Fanapi
CoverageDapu TownshipDapu Township
Evaluation Result
LandslideNon-landslideLandslideNon-landslide
Actual SituationLandslide156811696224
Non-landslide12,97567,610510274,832
PA of Grids with Landslide93.11%97.57%
PA of Grids without Landslide83.90%93.62%
Overall Accuracy84.09%93.67%
Table 11. Evaluation accuracy of postearthquake rainfall-induced landslide potential.
Table 11. Evaluation accuracy of postearthquake rainfall-induced landslide potential.
Post-Earthquake RainfallTyphoon Kalmaegi after 0305 EarthquakeTyphoon Kalmaegi after 0305 Earthquake
CoverageAlishan TownshipDapu Township
Evaluation Result
LandslideNon-landslideLandslideNon-landslide
Actual SituationLandslide334818491067
Non-landslide15,377221,58810,27552,018
PA of Grids with Landslide94.79%93.14%
PA of Grids without Landslide93.51%83.51%
Overall Accuracy93.53%83.65%
Table 12. Landslide area ratios of medium–high-level landslides induced by rainfall and earthquakes in Alishan Township.
Table 12. Landslide area ratios of medium–high-level landslides induced by rainfall and earthquakes in Alishan Township.
TriggerSingle RainfallSingle EarthquakePost-Earthquake Rainfall
Landslide Area or
Area Ratio
Typhoon Mindulle0520 Heavy Rainfall0602 Earthquake0206 Earthquake1122 EarthquakeTyphoon Kalmaegi after 0305 Earthquake
Total area of historical disaster areas (hectares)1394.61394.61394.61394.61394.61394.6
Area (hectares) in historical disaster area with medium-high landslide potential627.9480.1441.4397.3414.7509.1
Landslide area ratio of medium-high landslide potential0.450.3440.3170.2850.2970.365
Table 13. Landslide area ratios of medium–high-level landslides induced by rainfall and earthquakes in Dapu Township.
Table 13. Landslide area ratios of medium–high-level landslides induced by rainfall and earthquakes in Dapu Township.
TriggerSingle RainfallSingle EarthquakePost-Earthquake RainfallPost-Rainfall Earthquake
Landslide
Area or Area
Ratio
0726
Heavy Rain
Typhoon Nanmadol0520
Heavy Rain
Typhoon Megi0226
Earthquake
0206
Earthquake
0320
Earthquake
Typhoon Kalmaegi after 0305 Earthquake1105 Earthquake after Typhoon Morakot1108 Earthquake after Typhoon Fanapi
Total area of historical disaster areas (hectares)482.6482.6482.6482.6482.6482.6482.6482.6482.6482.6
Area (hectares) in historical disaster area with medium-high landslide potential7868.2100.891.372.478.384.7152.9132.999.4
Landslide area ratio of medium-high landslide potential0.1620.1410.2090.1890.150.1620.1760.3170.2750.206
Table 14. EAR × I3Rmax and PGA after rainfall and earthquakes.
Table 14. EAR × I3Rmax and PGA after rainfall and earthquakes.
(a) Alishan Township
Rainfall or (and) EarthquakeEAR × I3Rmax (mm2/3-h)PGA (gal)
MaximumMinimumMeanMaximumMinimumMean
Typhoon Mindulle138,44939,99386,407
0520 Heavy Rainfall85,60819,76745,671
0602 Earthquake 44534116
0206 Earthquake 21252126
1122 Earthquake 3183294
Typhoon Kalmaegi after 0305 Earthquake 131,57218,51882,08188735
(b) Dapu Township
Rainfall or (and) EarthquakeEAR × I3Rmax (mm2/3-h)PGA (gal)
MaximumMinimumMeanMaximumMinimumMean
0726 Heavy Rainfall23,269471311,161
Typhoon Nanmadol10,29948567259
0520 Heavy Rainfall67,37540,66758,007
Typhoon Megi78,59320,65939,822
0226 Earthquake 682555
0206 Earthquake 18572138
0320 Earthquake 24377146
Typhoon Kalmaegi after 0305 Earthquake193,711115,644164,8771691355
1105 Earthquake After Typhoon Morakot437,136158,541295,497681635
1108 Earthquake after Typhoon Fanapi35,863926820,9622082369
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tseng, C.-M.; Chen, Y.-R.; Chang, C.-M.; Yang, Y.-L.; Chen, Y.-R.; Hsieh, S.-C. Statistical Analysis of the Potential of Landslides Induced by Combination between Rainfall and Earthquakes. Water 2022, 14, 3691. https://doi.org/10.3390/w14223691

AMA Style

Tseng C-M, Chen Y-R, Chang C-M, Yang Y-L, Chen Y-R, Hsieh S-C. Statistical Analysis of the Potential of Landslides Induced by Combination between Rainfall and Earthquakes. Water. 2022; 14(22):3691. https://doi.org/10.3390/w14223691

Chicago/Turabian Style

Tseng, Chih-Ming, Yie-Ruey Chen, Chwen-Ming Chang, Ya-Ling Yang, Yu-Ru Chen, and Shun-Chieh Hsieh. 2022. "Statistical Analysis of the Potential of Landslides Induced by Combination between Rainfall and Earthquakes" Water 14, no. 22: 3691. https://doi.org/10.3390/w14223691

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop