Next Article in Journal
Membrane Distillation for Wastewater Treatment: A Mini Review
Next Article in Special Issue
Performance Assessment of Event-Based Ensemble Landslide Susceptibility Models in Shihmen Watershed, Taiwan
Previous Article in Journal
Block-Greedy and CNN Based Underwater Image Dehazing for Novel Depth Estimation and Optimal Ambient Light
Previous Article in Special Issue
The Effect of Temporal Characteristics on Developing a Practical Rainfall-Induced Landslide Potential Evaluation Model Using Random Forest Method
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Water
Volume 13
Issue 24
10.3390/w13243479
water-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating the Landslide Stability and Vegetation Recovery: Case Studies in the Tsengwen Reservoir Watershed in Taiwan

by
Chun-Hung Wu
Department of Water Resources Engineering and Conservation, Feng Chia University, Taichung 40724, Taiwan
Water 2021, 13(24), 3479; https://doi.org/10.3390/w13243479
Submission received: 22 October 2021 / Revised: 2 December 2021 / Accepted: 2 December 2021 / Published: 7 December 2021
(This article belongs to the Special Issue Natural Disasters Occurrence, Reduction, and Restoration in Mountain Regions)
Browse Figures
Review Reports Versions Notes

Abstract

:
The sediment yield from numerous landslides triggered in Taiwan’s mountainous regions by 2009 Typhoon Morakot have had substantial long-term impacts on the evolution of rivers. This study evaluated the long-term evolution of landslides induced by 2001 Typhoon Nari and 2009 Typhoon Morakot in the Tsengwen Reservoir Watershed by using multiannual landslide inventories and rainfall records for the 2001–2017 period. The landslide activity, vegetation recovery time, and the landslide spatiotemporal hotspot analyses were used in the study. Severe landslides most commonly occurred on 35–45° slopes at elevations of 1400–2000 m located within 500 m of the rivers. The average vegetation recovery time was 2.29 years, and landslides with vegetation recovery times exceeding 10 years were most frequently retrogressive landslide, riverbank landslides in sinuous reaches, and the core area of large landslides. The annual landslide area decline ratios after 2009 Typhoon Morakot in Southern Taiwan was 4.75% to 7.45%, and the time of landslide recovery in the Tsengwen reservoir watershed was predicted to be 28.48 years. Oscillating hotspots and coldspots occupied 95.8% of spatiotemporal patterns in the watershed area. The results indicate that landslides moved from hillslopes to rivers in the 2001–2017 period because the enormous amount of sediment deposited in rivers resulted in the change of river geomorphology and the riverbank landslides.

1. Introduction

Several serious landslide disaster events were caused by heavy rainfall events between 2001 and 2010, including 2001 Typhoon Toraji [1], 2004 Typhoon Mindulle [2], 2004 Typhoon Aere, and 2009 Typhoon Morakot [3], in the mountainous areas of central and southern Taiwan. The landslide events are believed to be related to the 1999 Chichi earthquake [4]. The return period of heavy rainfall brought by Typhoon Morakot (5–10 August 2009), which caused approximately 45,000 landslides concentrated in central and southern Taiwan, was estimated to be over 200 years [3]. Moreover, loose deposits from the numerous landslides in mountainous areas continues to affect watershed evolution and landslide recovery. As of 2021, mountainous areas in southern Taiwan remain under high risk of landslides and debris flows.
The analysis of landslide recovery or landslide evolution had been widely used in the long-term observation of landslide distribution after large earthquake events, including the 1999 Chichi earthquake [5], 2005 Kashimir earthquake [6], 2008 Wenchuan earthquake [7,8], and 2015 Gorkha earthquake [9]. Landslide evolution in the years following a large earthquake or extreme rainfall event in a watershed with dense landslide cases is the key determinant of watershed management and the mitigation of secondary geohazards [10]. Some artificial factors, including land use [11,12] and road development [13], were the significant factors for the evolution and reoccurrence of landslide. Some studies have been conducted on the characteristics of landslide evolution, including landslide activity [6,7], the spatiotemporal distributions of landslides [8,14], landslide recovery characteristics [15], and landslide evolution trends [10,16], with regard to severe landslides induced by large earthquakes or extreme rainfall. However, few studies have explored landslide evolution by examining spatiotemporal hotspots. Lin et al. (2017) [17] explained the distribution of landslide hotspots and coldspots on a catchment scale by using multiannual records of landslides, rainfall, and earthquakes for the 2003–2012 period in Taiwan. Conventional hotspot analysis can only describe the clustering pattern and intensity of a particular spatial location for a specific time interval; it cannot explain underlying trends over time [18]. The spatiotemporal evolution of landslide distribution patterns is challenging to investigate but merits scholarly attention.
New analysis methods by using the multi landslide inventories, including spatiotemporal hotspot and landslide activity, had been used to explore the characteristic of landslide evolution after large earthquake or extreme rainfall events [6,7,19]. The emerging hot spot analysis in ArcGIS Pro software has been employed in tracing the spread of COVID-19 [20], locating the source area of pollutant emissions [16], and assessing spatiotemporal changes in fisheries [21]. Emerging hotspot analysis was used in an assessment of long-term landslide evolution in Taiwan after 2009 Typhoon Morakot [19]. The results explained the spatiotemporal pattern and distribution of landslide hotspots and coldspots. Landslide-concentrated areas in the upper reaches of watersheds were determined to have poor landslide recovery and be highly susceptible to further landslides [19]. Emerging hotspot analysis can both explain changes in the spatiotemporal pattern and location of landslide hot spots and evaluate the intensity and location of landslide clustering. Landslide activity data have been employed in assessing the rates of and differences in landslide recovery after large earthquake-induced landslides [6,7] but not after extreme rainfall-induced landslide events.
The present study analyzed characteristics of landslide evolution by assessing landslide activity, estimating vegetation recovery time, and detecting the pattern and spatiotemporal distribution of landslides by using multiannual landslide inventories and long-term rainfall records from 2001 to 2017. Moreover, the characteristics of landslide evolution after large earthquakes were compared with the corresponding characteristics after extreme rainfall events. The Tsengwen Reservoir Watershed (TRW) suffered the most serious landslide disasters of all reservoir watersheds in Taiwan following 2009 Typhoon Morakot. The landslide evolution characteristics in the watershed with dense rainfall-induced landslides are poorly understood, but these characteristics represent information essential for watershed management and disaster prevention.

2. Research Area

The TRW (Figure 1) is a watershed located in the upper reaches of the Tsengwen River in southwestern Taiwan. The average elevation in the TRW is 959 m, and the average slope is 29.2°; 49.5% of the total area has a slope of greater than 30°. According to 1/5000 basin geological maps of Taiwan [22], the main stratigraphical formations in the TRW (Figure 2 and Table 1) include the Miocene-era Changchihkeng formation, Pliocene-era Ailiaochiao formation, and Miocene- to Pliocene-era Tangenshan sandstone (occupying 34.2%, 16.6%, and 16.4% of the watershed, respectively) [17]. The TRW is also surrounded by nine faults (Figure 2). Regarding land use types, forest, agriculture, development, rivers, and bare land account for 80.6%, 10.4%, 1.4%, 4.3%, and 3.3% of the total area based on the land use maps produced in 2008 by National Land Surveying and Mapping Center in Taiwan.
The average (standard deviation) annual precipitation from 2001 to 2017, based on the records of Matoushan Rainfall Station (Figure 1), was 2998.8 (934) mm. Rainfall characteristics in the TRW are non-uniform in space and time. The 3 years with the most precipitation were 2005 (4892 mm), 2008 (4118 mm), and 2001 (3764 mm). Those with the least precipitation were 2003 (1533 mm), 2013 (1735 mm), and 2002 (1878 mm; Figure 3). The difference between the precipitation in 2003 and that in 2005 is 1.12 times the average annual precipitation. The total precipitation in the rainy season (May–October; 2749 mm) was 91.7% of the average annual precipitation (2998.2 mm). The precipitation in August 2009, the largest monthly total in the considered period, was 2425 mm—70.2% of the annual precipitation.
Landslide and debris flow were the two main geohazards in the TRW, and the 1999 Chichi earthquake and heavy rainfall events during the 2000s were the main landslide triggers. Typhoon Nari (9–10 September 2001) and Typhoon Morakot (5–10 August 2009) caused the two main rainfall-induced landslide events during the 2000s. In the TRW, the accumulated precipitation brought by 2001 Typhoon Nari and 2009 Typhoon Morakot was 990.5 and 2158.5 mm (62.5% of the annual precipitation in 2009), respectively. The landslide ratios (i.e., the ratio of landslide area to watershed area) after 2001 Typhoon Nari and 2009 Typhoon Morakot were estimated to be 1.17% and 2.88%, respectively. Events with a landslide ratio exceeding 1.0% are considered serious disaster events in Taiwan [3]; thus, these two typhoons are regarded as serious disaster events within the TRW. The sediment deposition volume in the Tsengwen Reservoir increased by approximately 1.28 × 108 m3 from 2001 to 2017 (Figure 3). According to the field survey results, the majority of landslides can be classified as shallow landslide [23].

3. Materials and Method

3.1. Materials

A digital elevation model (DEM) and multiannual landslide inventories were the main study materials. The DEM had a spatial resolution of 5 m, and the 5 m × 5 m grid was the basic analysis unit in the study. The landslide inventories, corresponding to the 2001–2017 period, were produced by Taiwan’s Forestry Bureau, and landslides were identified from images captured by the Formosat-2 satellite with spatial resolutions of 2–5 m. Multiannual inventories of landslides in the TRW (for the years 2001–2017) were examined, and the minimum landslide area in the inventories was 100 m2. Landslides with an area of greater than 10 hectares (ha) were considered large landslides. The main types of landslide [24] included in the annual landslide inventories from 2001 to 2017 in Taiwan were fall, slide, and flow [19]. The deep-seated landslide cases, similar to the deep-seated gravitational slope deformation in central Taiwan [25], were not listed in the annual landslide inventories. Data on daily rainfall from 2001 to 2017 were collected from Matoushan Rainfall Station (Figure 1). The TRW contains 14 subwatersheds (Figure 1). The subwatersheds in the upper, middle, and lower reaches of the river are labeled as T01–T03, T04–T08, and T09–T14, respectively.
A recovery characteristics comparison of landslides induced by large earthquake and extreme rainfall events was conducted to explore the characteristic of landslide evolution. The large earthquake-induced landslides considered were the landslide cases after the 2005 Kashmir Earthquake (Mw = 7.6) in Pakistan [9] and the 2008 Wenchuan Earthquake (Mw = 7.9) in China [7,8,16]. The extreme rainfall-induced landslides used in the study were landslides occurring in the Chishan River Watershed [10], Ailiao River Watershed [19], and Taimali River Watershed [19] after 2009 Typhoon Morakot, and the TRW before and after 2009 Typhoon Morakot.

3.2. Landslide Activity

Landslide activity, which is useful in analyzing the spatiotemporal evolution of landslides, has been assessed by using the presence or absence of landslides in specific years [7,8] or by examining the ratio of the active landslide area to the total landslide area [26]. Following the landslide activity assessment methods used in the evolution discussion of earthquake-induced landslide cases [6,14,27], the criteria for the five types of landslide activity considered are presented in Table 2. The meaning of landslide activity used in a single hillslope scale was the current moving condition of the hillslope [25], but that used in a watershed scale was the vegetation recovery and the stability of loose materials deposited on the hillslopes [6,14,27]. The meanings of extremely active, very active, and active landslide were that the loose material deposited on the hillslope was instability or the vegetation recovery was too poor to stabilize the hillslope. The inventories of landslides in the TRW (2001–2017) were appropriate materials for assessing and comparing the activity of landslides induced by 2001 Typhoon Nari and 2009 Typhoon Morakot. The landslide activity induced by the typhoons was evaluated using multiannual landslide inventories for the 2001–2008 and 2009–2017 periods.

3.3. Landslide Frequency and Vegetation Recovery Time

Landslide frequency in the study was defined as the total occurrence of landslide identified in each grid from 2001 to 2017. The area that had been identified as landslide at least once from 2001 to 2017 in the study had been named as the landslide-identified area. The total occurrence of landslide-identified from 2001 to 2017 can be considered as the maximum vegetation recovery time (in years). The landslide recovery time had been estimated by using different factors, including landslide density [28], landslide rate [29], number of landslides [30], and landslide area [6]. Vegetation recovery is a significant factor to assess the landslide stability in long-term landslide evolution analysis [31] because the reinforcement of vegetation root contributes to the hillslope stability, especially in the shallow landslide cases [32,33]. The total occurrence of identified landslide of a grid was three, which meant that the grid was identified as landslides three times over the 17-year period. If the three landslide-identified years were discontinuous (e.g., 2005, 2008, and 2012), the vegetation recovery time was considered as 1 year. However, if the three landslide-identified years were continuous (e.g., 2002–2004), the vegetation recovery time was considered as 3 years. Herein, the maximum recovery time was taken as the vegetation recovery time in the study. Vegetation recovery time can be employed as an index to assess the difficulty of vegetation recovery in a watershed and in individual landslide cases.

3.4. Emerging Hotspot Analysis

The analysis of landslide spatiotemporal hotspots was based on the landslide spatiotemporal cube (STC) model, which represented landslide clustering patterns in each location at various time intervals. The STC model was composed of numerous 5 m × 5 m landslides-identified grids in the TRW from 2001 to 2017 in chronological order. The time series of landslide clustering intensity, represented by the combination of basic units at the same location, described the temporal evolution of landslides.
The emerging spatiotemporal mining method in ArcGIS Pro software was used in the STC model to explore the mechanism of temporal landslide evolution in the TRW. The Getis-Ord Gi statistic [34] was used to estimate the clustering intensity and to classify hot spot patterns (Table 3). The hotspot classification was based on the landslide clustering intensity in the neighborhood of specific 5 m × 5 m grids in both time and space, and specific basic units were designated as hotspots or coldspots if they featured high and low values of landslide clustering, respectively. Under the emerging spatiotemporal mining method, the consistency and intensity of landslide clustering in each time step was calculated, as was the significance of the autocorrelation and dependence in the vicinity of specific 5 m × 5 m grids. The method also used a space-time implementation of the Mann-Kendall statistic [35] to measure the intensity of landslide clustering in the neighborhood of specific basic units. The time step was set as 1 year, and the neighborhood distance was set as 5, 25, 50, 100, and 200 m for further comparison in the study.

4. Results

4.1. Multiannual Rainfall and Landslide Data

The Matoushan Rainfall Station (Figure 1) was used as the representative rainfall station in the TRW. The average annual rainfall recorded by the station from 2001 to 2017 (2998 mm) is 1.06 times that of the corresponding level from 1969 to 2020 (2835 mm) [36]. Table 4 presents data on daily and annual rainfall and on annual landslides in the TRW from 2001 to 2017. The annual rainfall distribution in the TRW was clearly non-uniform in time. The difference between the lowest annual rainfall (1533 mm in 2003) and the highest annual rainfall (4892 mm in 2005) was 3359 mm, which is 1.12 times the average annual rainfall. The daily rainfall corresponding to return periods of 1.11, 2, 5, 10, 20, 50, 100, and 200 years, based on station records from 1969 to 2020, were estimated to be 144.5, 300.6, 447.9, 545.4, 639.0, 760.1, 850.9, and 941.3 mm, respectively [33]. The maximum daily rainfall from 2001 to 2008 of 545.5 mm was observed twice, once during 2001 Typhoon Nari and once during 2005 Typhoon Haitang, and the return period was estimated to be 10–20 years. The maximum daily rainfall during Typhoon Morakot of 1206.0 mm was the highest amount of rainfall recorded at Matoushan Rainfall Station. The return period of daily rainfall from 2010 to 2017 was <10 years, and the maximum daily rainfall of 502.5 mm was recorded during 2015 Typhoon Soudelor.
The landslide area in 2001 (565.2 ha) was the highest for 2001–2008, and that in 2009 (1392.9 ha) marked the historical high recorded in the TRW. The average annual landslide area corresponding to the 2009–2017 period was 2.25 times that of the 2001–2008 period. The annual average number of landslides from 2009 to 2017 was 2.29 times larger than that from 2001 to 2008.
Southern Taiwan suffered serious landslide disaster during 2009 Typhoon Morakot, and the TRW was located in this region. The number of landslides in 2009 was 5.3 times larger than that in 2008. The number of the landslide cases with area greater than 100,000 m2, 1000–100,000 m2, and <1000 m2 in 2009 were 29, 1114, and 7, respectively. The average landslide length/width ratio in all landslide cases induced by 2009 Typhoon Morakot was 7.2, and this information showed that the majority type of landslide was rainfall-triggered slide.

4.2. Landslide Activity

Landslide activity data provide insights into the stability of rainfall-induced landslides for the analysis of the spatial distribution of landslide evolution. Figure 4 presents the spatial distribution of landslide activity after 2001 Typhoon Nari and 2009 Typhoon Morakot in the TRW. Areas of extremely active, very active, active, dormant, and inactive landslides after 2001 Typhoon Nari constituted 0.04, 40.8, 277.9, 939.5, and 561.4 ha, respectively. The corresponding areas after 2009 Typhoon Morakot were 54.0, 82.0, 424.7, 1577.0, and 762.2 ha, respectively. The accumulated area of extremely active, very active, and active landslides after 2009 Typhoon Morakot was 1.76 times that after 2001 Typhoon Nari, indicating that the landslides induced by Typhoon Morakot were more difficult to recover. The extremely active landslide area in subwatersheds in the upper reaches of the TRW increased substantially from 0 ha after 2001 Typhoon Nari to 42 ha after 2009 Typhoon Morakot, reflecting the difficulty of landslide recovery in the upstream subwatersheds.
A comparison of landslide activity after the 2005 Kashmir Earthquake [6] and 2009 Typhoon Morakot revealed that the proportions of unstable (extremely active, very active, or active) landslides and stable landslides (dormant or inactive landslides) were 28% and 72% in 2018 (13 years after the 2005 Kashmir earthquake). However, the corresponding proportions 8 to 9 years after typhoons in the TRW were 17.5–19.3% and 80.7–82.5%, respectively. These results indicate that vegetation recovery of landslides induced by 2009 Typhoon Morakot in the TRW was easier than vegetation recovery of landslides induced by 2005 Kashimir earthquake in Pakistan [6].

4.3. Relationship between Landslide Ratio and Geomorphological Factors

The relationship between earthquake-induced landslide area and geomorphological control factors, including elevation, slope, and distance to the river, had been discussed in analyses of the evolution of earthquake-induced landslides [6,7,8,18]. The temporal distribution of landslide ratios with respect to elevation, slope, and distance to the river from 2001 to 2017 in the TRW is shown in Figure 5. The landslide ratio distribution for the 2001–2007 period was used to observe the evolution of landslides induced by 2001 Typhoon Nari, and that for the 2009–2017 period was used to examine the evolution of landslides induced by 2009 Typhoon Morakot.
The landslide ratios were concentrated at 1800 to 2200 m elevation intervals after 2001 Typhoon Nari. However, from 2003 to 2007, they were concentrated at 800–1400 m elevation intervals. The landslide ratios at elevations below 1000 m and above 1600 m decreased notably during this period. After 2009 Typhoon Morakot, the landslide ratios were concentrated at 600 m to 1400 m elevation intervals. From 2010 to 2017, they were concentrated at 1200 m to 1800 m elevation intervals. The landslide ratios at elevations below 1200 m decreased notably, and those at elevations above 1800 m exhibited small changes from 2010 to 2017.
Only the distribution of extremely active landslides after 2009 Typhoon Morakot is plotted herein (Figure 5) because the area of extremely active landslides after 2001 Typhoon Nari was quite small (0.04 ha). Extremely active landslides were clustered at 1400 to 2000 m elevation intervals after Typhoon Morakot. The extremely active or very active landslides were clustered at moderately high elevations after large earthquakes. For example, they were clustered from 1000 to 1600 m after the 2005 Kashmir Earthquake [6] and below 2000 m after the 2008 Wenchuan Earthquake [7]. Overall, the clustered area and locations of extremely active or very active landslides after large earthquakes or extreme rainfall events was very similar.
The landslide ratios were concentrated at two slope intervals, 35° to 45° and 55° to 60°, after 2001 Typhoon Nari and 2009 Typhoon Morakot. The landslide ratios were clustered at the same intervals from 2003 to 2007 and from 2010 to 2017. The extremely active landslides were concentrated at 35–45° slope intervals after 2009 Typhoon Morakot. Extremely active or very active landslides were also clustered at 30–40° after the 2005 Kashmir Earthquake [6] and at 30° to 50° after the 2008 Wenchuan Earthquake [7].
In Taiwan, loose debris from landslides after heavy rainfall events was generally deposited in rivers and gullies [10,15]. The TRW was divided into 10 buffer areas with 100 m distance intervals from the rivers. The landslide ratios in all river buffer intervals after 2001 Typhoon Nari ranged from 0.88% to 1.63%, but those after 2009 Typhoon Morakot were all greater than 1.86%. The highest landslide ratio in all river buffer intervals after 2001 Typhoon Nari was 1.63% within 1000 m of the river, and that after 2009 Typhoon Morakot was 6.3% within 100 m of the river.
Extremely active landslides were clustered along the fault line after the 2005 Kashmir Earthquake [6] but were concentrated along rivers in the TRW after 2009 Typhoon Morakot. This indicates that large pieces of loose material from landslides were transported by excess surface runoff and flood discharge and deposited in the rivers, thereby inducing riverbank landslides and change of river geomorphology.

4.4. Landslide Frequency and Vegetation Recovery Time

The distribution of landslide frequency from 2001 to 2017 in the TRW is shown in Figure 6. The total landslide-identified area was 26.4 km2 (5.5% of the total area), and the average landslide frequency was 2.29. The area with landslide frequency = 1 constituted 16.28 km2, which was 61.8% of the landslide-identified area. The area with landslide frequency greater than 10 constituted 0.394 km2. Those landslides were most commonly retrogressive landslide or riverbank landslides in the sinuous reaches and the core area of the large landslide cases.
The landslide frequency in the subwatersheds in the upper reaches of the TRW (2.28–3.73) differed substantially from the distributions in the subwatersheds in the middle and lower reaches (1.82–1.97 and 1.72–2.35, respectively). The average vegetation recovery time in the subwatersheds in the upper reaches was 1.41 to 1.59 times longer than that in the subwatersheds in the middle and lower reaches. The subwatersheds in the upper reaches had the highest landslide ratio (Figure 6b,c) and smallest catchment area, indicating that the large amount of sediment yield from numerous landslides was deposited in the narrow upper reaches, and that resulted in the change of river flowing path and the riverbank landslides.
Large landslides and riverbank landslides in the sinuous reaches had the longest vegetation recovery time of all landslides in the TRW. The 2009 landslide inventory contained 22 large landslide cases, the mean vegetation recovery time of which was 3.29 years. Two large landslide cases with mean vegetation recovery times of greater than 6 years occurred in the sinuous upper reaches of the subwatersheds. The large landslide occurring in a gully source area in the T01 subwatershed (Figure 7) can be a typical example to explain the vegetation recovery condition. The area and relief of the large landslide were 35.7 ha and 503.2 m, and the average slope was 42.1°. The stratigraphical formation and lithology of the large landslide were the Nanchuang formation and the sandstone and shale. The activity type of the large landslide was very active. The mean vegetation recovery time of the large landslide was 10.25 years. The vegetation recovery time in the boundary area of the large landslide was less than 3 years, indicating that landslides were easily re-induced in this boundary area. However, the vegetation recovery time in the core area of the large landslide case was over 7 years, and vegetation recovery in the core area of the large landslide was difficult.
The river buffer area was defined as the area within 300 m of the river. The landslide-identified area occupied 7.5% of the total area in the river buffer area in the TRW, and the mean vegetation recovery time in the river buffer area was 2.65 years. The proportion of the landslide-identified area and the landslide frequency in the river buffer area were slighter larger than those in the TRW. The enormous amount of loose material from landslides was deposited in the confluence of river waters and sinuous reaches, resulting in riverbank landslides, especially in the upper reaches of landslide-concentrated subwatersheds. The riverbank landslide occurring in the confluence downstream of the T01 and T02 subwatersheds (Figure 8) is a typical example demonstrating the impacts of large sediment deposits in sinuous reaches on landslide recovery. The area and relief of the large landslide were 0.92 ha and 269.7 m, and the average slope was 32.5°. The stratigraphical formation and lithology of the large landslide were the Changchihkeng formation and the sandstone and shale. The activity type of the large landslide was active and very active. The riverbank landslide was located in a sinuous reach with the sinuosity index = 1.62, and the sediment deposition depth from 2009 to 2011 in this sinuous reach was estimated as 4.32 m. The mean vegetation recovery time of the riverbank landslide (Figure 8) was 6.79 years (2.56 times longer than that in the TRW). The considerable amount of sediment deposition in the river in the vicinity of the riverbank landslide resulted in the formation of river meanders and several new riverbank landslides.

4.5. Landslide Spatiotemporal Hotspot Distribution and Trend

The total number of 5 m × 5 m grids in the landslide STC model, collected from the aggregated spatiotemporal data corresponding to all landslide grids from 2001 to 2017 in the TRW, was 17,861,645. The spatiotemporal hot spot distribution in the TRW is displayed in Figure 9. Neighborhood distance settings in analyses of spatiotemporal landslide hotspots have not yet been suggested. Therefore, these settings were based on a performance comparison of results obtained through analyses performed using various neighborhood distances. Specifically, the distances were set to 5, 25, 50, 100, and 200 m.
The no pattern area in the analysis with a 5 m neighborhood distance occupied 85.6% of the TRW, but that with a neighborhood distance of ≥25 m dropped to 2.4–3.4% (Figure 9d). This result indicates that the 5 m neighborhood distance setting in the analysis was excessively short for detecting landslide hotspots. Furthermore, the 5 m setting yielded only five spatiotemporal patterns for the entire TRW. By contrast, when neighborhood distances of ≥25 m were used, 10–12 spatiotemporal patterns were generated (Figure 9a–c). The main hotspot patterns in the results obtained using a 5 m neighborhood distance were sporadic hotspot (SHS) and consecutive hotspot (CHS). The corresponding patterns in the results obtained using a neighborhood distance of ≥25 m were oscillating hotspot (OHS) and oscillating coldspot (OCS). The spatiotemporal distribution of hotspots in the results obtained using 25, 50, 100, and 200 m neighborhood distances were similar yet distinct from that in the results obtained using the 5 m neighborhood distance. Therefore, the neighborhood distance setting in spatiotemporal analyses of landslide hotspots was suggested to be ≥25 m from the result in the study.
The numbers of 5 m × 5 m grids corresponding to coldspot trends in the analyses performed using neighborhood distances of 200 and 25 m were 2.46 and 3.80 times the number of bins corresponding to hotspot trends. These results indicate that the landslide clustering intensity in the TRW was low. The summed proportion of OCS and OHS in analyses with neighborhood distances of ≥25 m in the TRW was greater than 95.8%. In other words, the dominant spatiotemporal landslide pattern in the TRW was oscillating. The spatiotemporal hot spot pattern of the large landslide in the T01 subwatershed and the riverbank landslide in the T03 subwatershed are presented in Figure 7d and Figure 8d, respectively. Oscillating hotspots and coldspots characterized these two landslides.

5. Discussion

The landslide recovery rate and landslide activity are the keys to understanding the evolution of watersheds with frequent and dense landslides. Studies have addressed the rate of landslide recovery induced by large earthquake events [6,7,8], but few have examined the rate of recovery from landslides induced by extreme rainfall events [19]. The present study discussed and compared differences in the distribution rate of recovery from active landslides induced by large earthquakes and extreme rainfall events.
The majority of landslide types after the 2005 Kashimir earthquake were rock falls and rock slides [6], and those after the 2008 Wenchuan earthquake were debris slides [7]. The majority of landslide types after 2009 Typhoon Morakot in Taiwan were debris slides. Herein, the landslide area decline ratio was defined as the landslide area in the specific year to the landslide area in the year a large earthquake or extreme rainfall event occurred. For example, the landslide areas in the TRW in 2009 and 2010 were 1392.9 and 775.6 ha, and the landslide area decline ratio in 2010 was 55.7%. The landslide area decline ratios for the large earthquake- and extreme rainfall-induced landslides are shown in Figure 10. For each landslide event, the linear fitting equation was estimated. The slope coefficient of this equation can be regarded as the average annual landslide area decline ratio, which was used in the comparison of landslide recovery.
The annual landslide area decline ratios corresponding to the large earthquake events (5.02–9.56%) were slightly larger than those corresponding to the extreme rainfall events (4.75–7.45%). The R2 values of the linear fitting equations for large earthquake-induced landslides exceeded those of the corresponding equations for extreme rainfall-induced landslides. This indicates that recovery from large earthquake-induced landslides was more stable than recovery was from extreme rainfall-induced landslides. If the declining landslide area in the TRW after Typhoon Morakot follows the average annual ratio of landslide area decline, 28.48 years would be required for a return to the pretyphoon landslide area.
The characteristic and evolution comparison of landslide induced by large earthquake events in Northern Pakistan and Sichuan, China and 2009 Typhoon Morakot in TRW are shown in Table 5. The main reasons for the difference of landslide evolution between the large earthquake-induced or extreme rainfall-induced landslide cases were the landslide-inducing factors and the deposition location of loose materials yield from landslide. The inducing factor for the earthquake-induced slope instability was the ground shaking force to reduce the shear strength of slope [37], while that for the rainfall-induced slope instability was the infiltration water [38]. The extreme active, very active, and active landslide after large earthquake-induced landslide events were located in the neighboring area of faults, including the Muzaffarbad fault in North Pakistan [6] and the Longmenshan fault in China [7,16], but those after extreme rainfall-induced landslide events in the TRW were located in the neighboring area of rivers and gullies (Figure 4). The percentage of the extreme active, very active, and active landslide areas located in the area within 500 m of rivers or the source area of large landslide cases after 2001 Typhoon Nari in the TRW was 74.9%, and those after 2009 Typhoon Morakot was 85.3%. The pattern and distribution of large earthquake-induced [6,7] or extreme rainfall-induced landslide spatiotemporal hotspot demonstrated that the landslide hotspot had been moving from mountain regions to the neighborhood area of rivers, but the time needed was the main difference between the large earthquake-induced or extreme rainfall-induced landslide events. The landslide ratio in the area within 500 m to rivers in the first year after 2009 Typhoon Morakot in the TRW ranged from 2.4% to 6.3% (Figure 5c), which was much larger than that in the first year after the 2008 Wenchuan earthquake-induced landslide events in China [16]. The evolution of disaster types after extreme rainfall-induced landslide events in the TRW was from numerous landslide cases in the hillslope to the river geomorphological changes and riverbank landslide due to much sediment deposition in the rivers, especially in the upstream watersheds. This was also the main reasons for the difficulty of landslide recovery in the T01 and T02 subwatersheds. The key to reduce the disaster in the following years after extreme rainfall-induced landslide events based on the results in the study was how to control the huge amount of loose material yield from numerous landslides and that deposited in the rivers.

6. Conclusions

This study evaluated landslide activity and vegetation recovery time and detected the spatiotemporal hotspots of extreme rainfall-induced landslides in the TRW by using annual landslide inventories and long-term rainfall records from 2001 to 2017. Areas of extreme rainfall-induced landslides and landslides in general decreased consistently following 2001 Typhoon Nari and 2009 Typhoon Morakot. The return period of rainfall events caused a slight increase in the area of number of landslides; during 2002–2008, it was greater than 10 years. In the 2010–2017 period, it was greater than 2 years. The area of extremely active landslides in the TRW after 2009 Typhoon Morakot was notably larger than that after 2001 Typhoon Nari, and extremely active landslides were clustered in the subwatersheds in the upper reaches of the TRW. The study also discussed the relationship between temporal landslide distribution and geomorphological factors, including elevation, slope, and distance to the river. Landslides in the years following typhoon events were concentrated at 1400 to 2000 m elevations on 35° to 45° slopes within 500 m of the river. The average vegetation recovery time in the TRW was 2.29 years, and landslides with vegetation recovery times of greater than 10 years were commonly retrogressive landslide, riverbank landslides in the sinuous reaches, and landslides in the core area of large landslides. The time required to recover from landslides in the subwatersheds in the upper reaches was 1.41 to 1.59 times longer than the time required to recover from landslides in subwatersheds in the middle and lower reaches. The main spatiotemporal pattern of landslides in the TRW, including in the subwatersheds in the upper reaches, was characterized by oscillating hotspots and coldspots. The annual landslide area decline ratio in the TRW after Typhoon Morakot was estimated to be 94.08%, and approximately 28.5 years is required for a return to the pretyphoon landslide area. Findings on the characteristics of landslide recovery, the distribution of landslide activity, and the spatiotemporal patterns of landslides are useful for watershed management and disaster prevention in areas where rainfall-induced landslides commonly occur.

Funding

This research was funded by the Ministry of Science and Technology of Taiwan (R.O.C.), grant number MOST 109-2625-M-035-005 and the Project Research in Feng Chia University (grant number: 20H00710 and 21H00710).

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. Cheng, J.D.; Haung, Y.C.; Wu, H.L.; Yen, J.L.; Chang, C.H. Hydrometeorological and landuse attributes of debris flows and debris floods during typhoon Toraji, July 29–30, 2001 in central Taiwan. J. Hydrol. 2005, 306, 161–173. [Google Scholar] [CrossRef]
  2. Milliman, J.D.; Lee, T.Y.; Huang, J.C.; Kao, S.J. Impact of catastrophic events on small mountainous rivers: Temporal and spatial variations in suspended- and dissolved-solid fluxes along the Choshui River, central western Taiwan, during typhoon Mindulle, July 2–6, 2004. Geochim. Cosmochim. Acta 2017, 205, 272–294. [Google Scholar] [CrossRef]
  3. Wu, C.H.; Chen, S.C.; Chou, H.T. Geomorphologic characteristics of catastrophic landslides during typhoon Morakot in the Kaoping Watershed, Taiwan. Eng. Geol. 2011, 123, 13–21. [Google Scholar] [CrossRef]
  4. 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]
  5. Lin, C.Y.; Chuang, C.W.; Lin, W.T.; Chou, W.C. Vegetation recovery and landscape change assessment at Chiufenershan landslide area caused by Chichi earthquake in central Taiwan. Nat. Hazards 2010, 53, 175–194. [Google Scholar] [CrossRef]
  6. Shafique, M. Spatial and temporal evolution of coseismic landslides after the 2005 Kashmir earthquake. Geomorphology 2020, 362, 107228. [Google Scholar] [CrossRef]
  7. Yang, W.; Qi, W.; Wang, M.; Zhang, J.; Zhang, Y. Spatial and temporal analyses of post-seismic landslide changes near the epicentre of the Wenchuan earthquake. Geomorphology 2017, 276, 8–15. [Google Scholar] [CrossRef]
  8. Li, C.R.; Wang, M.; Liu, K. A decadal evolution of landslides and debris flows after the Wenchuan earthquake. Geomorphology 2018, 323, 1–12. [Google Scholar] [CrossRef]
  9. Kincey, M.E.; Rosser, N.J.; Robinson, T.R.; Densmore, A.L.; Shrestha, R.; Pujara, D.S. Evolution of coseismic and post-seismic landsliding after the 2015 Mw 7.8 Gorkha earthquake, Nepal. J. Geophys. Res. Earth Surf. 2021, 126, e2020JF005803. [Google Scholar] [CrossRef]
  10. Wu, C.H. Landslide susceptibility based on extreme rainfall-induced landslide inventories and the following landslide evolution. Water 2019, 11, 2609. [Google Scholar] [CrossRef] [Green Version]
  11. Glade, T. Landslide occurrence as a response to land use change: A review of evidence from New Zealand. Catena 2003, 51, 297–314. [Google Scholar] [CrossRef] [Green Version]
  12. Van Beek, L.P.H.; Van Asch, T.W. Regional assessment of the effects of land-use change on landslide hazard by means of physically based modelling. Nat. Hazards 2004, 31, 289–304. [Google Scholar] [CrossRef]
  13. Phillips, C.; Hales, T.; Smith, H.; Basher, L. Shallow landslides and vegetation at the catchment scale: A perspective. Ecol. Eng. 2021, 173, 106436. [Google Scholar] [CrossRef]
  14. Valenzuela, P.; Domínguez-Cuesta, M.J.; García, M.A.M.; Jiménez-Sánchez, M. A spatio-temporal landslide inventory for the NW of Spain: BAPA database. Geomorphology 2017, 293, 11–23. [Google Scholar] [CrossRef] [Green Version]
  15. Wu, C.H. Comparison and evolution of extreme rainfall-induced landslides in Taiwan. ISPRS Int. J. Geo-Inf. 2017, 6, 367. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, M.; Tang, C.; Xiong, J.; Shi, Q.Y.; Gong, L.F.; Wang, X.D.; Tie, Y. The long-term evolution of landslide activity near the epicentral area of the 2008 Wenchuan earthquake in China. Geomorphology 2020, 367, 107317. [Google Scholar] [CrossRef]
  17. Lin, S.C.; Ke, M.C.; Lo, C.M. Evolution of landslide hotspots in Taiwan. Landslides 2017, 14, 1491–1501. [Google Scholar] [CrossRef]
  18. Cheng, S.; Zhang, B.; Peng, P.; Yang, Z.; Lu, F. Spatiotemporal evolution pattern detection for heavy-duty diesel truck emissions using trajectory mining: A case study of Tianjin, China. J. Clean. Prod. 2020, 244, 118654. [Google Scholar] [CrossRef]
  19. Wu, C.H.; Lin, C.Y. Spatiotemporal Hotspots and Decadal Evolution of Extreme Rainfall-Induced Landslides: Case Studies in Southern Taiwan. Water 2021, 13, 2090. [Google Scholar] [CrossRef]
  20. Purwanto, P.; Utaya, S.; Handoyo, B.; Bachri, S.; Sari, I.; Sastro, K.; Utomo, B.; Aldianto, Y.E. Spatiotemporal analysis of COVID-19 spread with emerging hotspot analysis and space–time cube models in East Java, Indonesia. ISPRS Int. J. Geo-Inf. 2021, 10, 133. [Google Scholar] [CrossRef]
  21. Everett, B.I.; Fennessy, S.T.; van den Heever, N. Using hotspot analysis to track changes in the crustacean fishery off KwaZulu-Natal, South Africa. Reg. Stud. 2021, 41, 101553. [Google Scholar] [CrossRef]
  22. Sinotech Consultants. Geological Investigation and Database Construction for the UpstreamWatershed of Flood-Prone; Central Geological Survey, Ministry of Economic Affairs, R.O.C: Taipei, Taiwan, 2007. [Google Scholar]
  23. Tsai, K.J.; Lee, M.H.; Chen, Y.R.; Chen, X.Y. A Study on the Sediment Quanity Variation of Large Scale Landslides Induced by Extremely Heavy Rainfall Occurred at Zengwen Watershed in Southern Taiwan. Min. Metall. 2017, 61, 8–24. (In Chinese) [Google Scholar]
  24. Varnes, D.J. Slope Movement Types and Processes. In Special Report 176: Landslides: Analysis and Control; Schuster, R.L., Krizek, R.J., Eds.; Transportation and Road Research Board, National Academy of Science: Washington, DC, USA, 1978; pp. 11–33. [Google Scholar]
  25. Lin, H.H.; Lin, M.L.; Lu, J.H.; Chi, C.C.; Fei, L.Y. Deep-seated gravitational slope deformation in Lushan, Taiwan: Transformation from cleavage-controlled to weakened rockmass-controlled deformation. Eng. Geol. 2020, 264, 105387. [Google Scholar] [CrossRef]
  26. Tang, C.X.; VanWesten, C.J.; Tanyas, H.; Jetten, V.G. Analysing post-earthquake landslide activity using multi-temporal landslide inventories near the epicentral area of the 2008 Wenchuan earthquake. Nat. Hazards Earth Syst. Sci. 2016, 16, 2641–2655. [Google Scholar] [CrossRef] [Green Version]
  27. Hervás, J.; Barredo, J.I.; Rosin, P.L.; Pasuto, A.; Mantovani, F.; Silvano, S. Monitoring landslides from optical remotely sensed imagery: The case history of Tessina landslide, Italy. Geomorphology 2003, 54, 63–75. [Google Scholar] [CrossRef]
  28. Marc, O.; Behling, R.; Andermann, C.; Turowski, J.M.; Illien, L.; Roessner, S. Long-term erosion of the Nepal Himalayas by bedrock landsliding: The role of monsoons, earthquakes and giant landslides. Earth Surf. Dyn. 2019, 7, 107–128. [Google Scholar] [CrossRef] [Green Version]
  29. Barth, S.; Geertsema, M.; Bevington, A.R.; Bird, A.L.; Clague, J.J.; Millard, T. Landslide response to the 27 October 2012 earthquake (MW 7.8), southern Haida Gwaii, British Columbia, Canada. Landslides 2020, 17, 517–526. [Google Scholar] [CrossRef] [Green Version]
  30. Saba, S.B.; van der Meijde, M.; van derWerff, H. Spatiotemporal landslide detection for the 2005 Kashmir earthquake region. Geomorphology 2010, 124, 17–25. [Google Scholar] [CrossRef]
  31. Yang, W.; Qi, W.; Zhou, J. Decreased post-seismic landslides linked to vegetation recovery after the 2008 Wenchuan earthquake. Ecol. Indic. 2018, 89, 438–444. [Google Scholar] [CrossRef]
  32. Masi, E.B.; Segoni, S.; Tofani, V. Root Reinforcement in Slope Stability Models: A Review. Geosciences 2021, 11, 212. [Google Scholar] [CrossRef]
  33. Saito, H.; Uchiyama, S.; Teshirogic, K. Rapid vegetation recovery at landslide scars detected by multitemporal high-resolution satellite imagery at Aso volcano, Japan. Geomorphology 2021, 398, 107989. [Google Scholar] [CrossRef]
  34. Getis, A.; Ord, J.K. The analysis of spatial association by use of distance statistics. Geogr. Anal. 1992, 24, 189–206. [Google Scholar] [CrossRef]
  35. Mann, H.B. Nonparametric tests against trend. Econometrica 1945, 13, 245. [Google Scholar] [CrossRef]
  36. Sheen, S.W. Frequency Analysis of Annual 1-day Maximum Rainfall in the Upstream Drainage Area of Zeng-Wen River. Crop. Environ. Bioinform. 2016, 13, 1–12. (In Chinese) [Google Scholar]
  37. Kamp, U.; Growley, B.J.; Khattak, G.A.; Owen, L.A. GIS-based landslide susceptibility mapping for the 2005 Kashmir earthquake region. Geomorphology 2008, 101, 631–642. [Google Scholar] [CrossRef]
  38. Chiu, Y.J.; Lee, H.Y.; Wang, T.L.; Yu, J.; Lin, Y.-T.; Yuan, Y. Modeling Sediment Yields and Stream Stability Due to Sediment-Related Disaster in Shihmen Reservoir Watershed in Taiwan. Water 2019, 11, 332. [Google Scholar] [CrossRef] [Green Version]
  39. Korup, O.; Wang, G.H. Chapter 8-Multiple Landslide-Damming Episodes. In Landslide Hazards, Risks, and Disasters; Elsevier: Amsterdam, The Netherlands, 2015; pp. 241–261. [Google Scholar]
  40. Chan, H.C.; Chang, C.C.; Chen, S.C.; Wei, Y.S.; Wang, Z.B.; Lee, T.S. Investigation and analysis of the characteristics of shallow landslides in mountainous areas of Taiwan. J. Chin. Soil Water Conserv. 2015, 46, 19–28. (In Chinese) [Google Scholar]
Water 13 03479 g001 550
Figure 1. The distribution of elevation, rivers, and landslide after 2001 Typhoon Nari and 2009 Typhoon Morakot in the Tsengwen Reservoir watershed (Abbreviated as TRW). The down-left figure shows the location of TRW in Taiwan, and the up-left figure shows the names of sub-watersheds in the TRW.
Figure 1. The distribution of elevation, rivers, and landslide after 2001 Typhoon Nari and 2009 Typhoon Morakot in the Tsengwen Reservoir watershed (Abbreviated as TRW). The down-left figure shows the location of TRW in Taiwan, and the up-left figure shows the names of sub-watersheds in the TRW.
Water 13 03479 g001
Water 13 03479 g002 550
Figure 2. The distribution of geological settings and faults in the TRW.
Figure 2. The distribution of geological settings and faults in the TRW.
Water 13 03479 g002
Water 13 03479 g003 550
Figure 3. The distribution of monthly precipitation, landslide ratio, and reservoir deposition ratio from 2001 to 2017 in the TRW. The up, middle, and down figures are the distribution of reservoir deposition ratio, annual landslide ratio, and monthly precipitation from 2001 to 2017 in the TRW.
Figure 3. The distribution of monthly precipitation, landslide ratio, and reservoir deposition ratio from 2001 to 2017 in the TRW. The up, middle, and down figures are the distribution of reservoir deposition ratio, annual landslide ratio, and monthly precipitation from 2001 to 2017 in the TRW.
Water 13 03479 g003
Water 13 03479 g004 550
Figure 4. The distribution of landslide activity in the TRW. (a) is the distribution of landslide activity after 2001 Typhoon Nari and (b) is after 2009 Typhoon Morakot.
Figure 4. The distribution of landslide activity in the TRW. (a) is the distribution of landslide activity after 2001 Typhoon Nari and (b) is after 2009 Typhoon Morakot.
Water 13 03479 g004
Water 13 03479 g005 550
Figure 5. The relationship between landslide ratios and geomorphological factors, including elevation in (a), slope in (b), and distance to the rivers in (c), from 2001 to 2017 in the TRW. The EA line (green) means the extremely active landslide after 2009 Typhoon Morakot.
Figure 5. The relationship between landslide ratios and geomorphological factors, including elevation in (a), slope in (b), and distance to the rivers in (c), from 2001 to 2017 in the TRW. The EA line (green) means the extremely active landslide after 2009 Typhoon Morakot.
Water 13 03479 g005
Water 13 03479 g006 550
Figure 6. The landslide frequency from 2001 to 2017 (a) and the landslide ratio in each subwatershed from 2009 to 2017 (b,c) in the TRW.
Figure 6. The landslide frequency from 2001 to 2017 (a) and the landslide ratio in each subwatershed from 2009 to 2017 (b,c) in the TRW.
Water 13 03479 g006
Water 13 03479 g007 550
Figure 7. The aerial images in 2009 (a) and 2017 (b), the distribution of vegetation recovery time (c), and the spatiotemporal hotspots (d) in the large landslide in the T01 subwatersheds (Figure 6).
Figure 7. The aerial images in 2009 (a) and 2017 (b), the distribution of vegetation recovery time (c), and the spatiotemporal hotspots (d) in the large landslide in the T01 subwatersheds (Figure 6).
Water 13 03479 g007
Water 13 03479 g008 550
Figure 8. The aerial images in 2009 (a) and 2017 (b), the distribution of vegetation recovery time (c), and the spatiotemporal hotspots (d) in the riverbank landslide in the T03 subwatershed (Figure 6).
Figure 8. The aerial images in 2009 (a) and 2017 (b), the distribution of vegetation recovery time (c), and the spatiotemporal hotspots (d) in the riverbank landslide in the T03 subwatershed (Figure 6).
Water 13 03479 g008
Water 13 03479 g009 550
Figure 9. The landslide spatiotemporal hotspot distributions (a–c) and the occupied percentage (d) in the TRW.
Figure 9. The landslide spatiotemporal hotspot distributions (a–c) and the occupied percentage (d) in the TRW.
Water 13 03479 g009
Water 13 03479 g010 550
Figure 10. The area declining ratio fitting equation of large earthquake-induced and extreme rainfall-induced landslide events.
Figure 10. The area declining ratio fitting equation of large earthquake-induced and extreme rainfall-induced landslide events.
Water 13 03479 g010
Table
Table 1. Statistical data of landslide induced by 2001 Typhoon Nari and 2009 Typhoon Morakot in each stratigraphical division in the TRW.
Table 1. Statistical data of landslide induced by 2001 Typhoon Nari and 2009 Typhoon Morakot in each stratigraphical division in the TRW.
Stratigraphical DivisionLithologyOcc. Per. (%)Landslide Ratio (%)
In 2001In 2009
Water Area-3.070.000.27
AlluviumGravel, sand, silt, and clay0.700.801.96
Ailiaochiao FormationThin alternation of siltstone and shale16.650.802.52
Changchihkeng FormationSandstone, sandstone interbedded with shale34.161.373.48
Chutouchi FormationMuddy sandstone, interbedded muddy sandstone and shale2.760.365.43
Hunghuatzu FormationThick-bedded fine sandstone or siltstone, thick alternation of fine sandstone and siltstone2.480.744.13
Maopu ShaleShale with thin-bedded sandstone0.990.175.69
Nanchuang FormationAlternation of sandstone and shale12.441.932.91
Peiliao ShaleShale and sandy shale5.580.341.00
Sanming ShaleShale, shale intercalated with thin-bedded siltstone0.494.725.30
Terrace DepositsSand, silt, mud and gravel1.190.020.00
Tangenshan SandstoneThick-bedded massive sandstone and muddy sandstone16.391.452.65
Yenshuikeng ShaleMassive shale, occasionally intercalated with thin-bedded siltstone3.110.442.13
Note: The Occ. Per. means the occupied percentage of each stratigraphical division area to the watershed area.
Table
Table 2. The criteria of landslide activity in the study.
Table 2. The criteria of landslide activity in the study.
Activity TypeCriteria
Extremely ActiveLandslides present in the annual inventories from 2001 to 2008 or from 2009 to 2017.
Very ActiveLandslides present in the annual inventories from 2005 to 2008 or from 2013 to 2017.
ActiveLandslides present in the annual inventories of 2007 to 2008, or only in the annual inventory of 2008 for 2001 Typhoon Nari. Landslides present in the annual inventories of 2016 to 2017, or only in the annual inventory of 2017 for 2009 Typhoon Morakot.
DormantLandslides present in one or more annual inventories of 2001 to 2007 or 2009 to 2016 but absent in the annual inventory of 2008 or 2017.
InactiveLandslides present in the annual inventory of 2001 or 2009 but absent in the annual inventories from 2002 to 2008 or 2010 to 2017.
Table
Table 3. The classifications and definition of emerging landslide hotspots and coldspots in the study.
Table 3. The classifications and definition of emerging landslide hotspots and coldspots in the study.
PatternDefinition
Consecutive
(CHS or CCS) *		A landslide grid with a single uninterrupted run of statistically significant hotspot or coldspot grids in the final year during the research time period. The landslide grid has never been a statistically significant hotspot or coldspot before the final hotspot or coldspot run.
Diminishing
(DHS or DCS)		A landslide grid that has been a statistically significant hotspot or coldspot for 90% of the research time period, including the final year. In addition, the clustering intensity of landslides in each year is decreasing (increasing) overall, and that decrease (increase) is statistically significant.
Historical
(HHS or HCS)		The most recent year is not hotspot or coldspot, but at least 90% of the research time period has been a statistically significant hotspot or coldspot.
Intensifying
(IHS or ICS)		A landslide grid that has been a statistically significant hotspot or coldspot for 90% of the research time period. In addition, the clustering intensity of landslide for each year increased (decreased) overall, and that increase (decrease) was statistically significant.
New
(NHS or NCS)		A landslide grid identified as a statistically significant hotspot or coldspot since the first year of the research time period but was not previously identified as a statistically significant hotspot or coldspot.
Oscillating
(OHS or OCS)		A statistically significant hotspot or coldspot for the final year that has a history of also being a statistically significant coldspot or hotspot during a prior year. Less than 90% of the research time period has been a statistically significant hotspot or coldspot.
Persistent
(PHS or PCS)		A landslide grid that has been a statistically significant hotspot or coldspot for 90% of the research time period with no discernible trend indicating an increase or decrease in the clustering intensity of landslide over time.
Sporadic
(SHS or SCS)		A landslide grid that is an on-again then off-again hotspot or coldspot. Less than 90% of the research time period has been a statistically significant hotspot or coldspot, and none of the time-step intervals have been a statistically significant coldspot or hotspot.
No pattern detected (No)The analysis grid does not fit any definition of hotspot or coldspot classifications
*: CHS and CCS are the abbreviations of consecutive hotspot and consecutive coldspot. The regulation of abbreviation was applied to every hotspot and coldspot pattern in the study.
Table
Table 4. The statistical data of annual rainfall and landslide from 2001 to 2017 in the TRW.
Table 4. The statistical data of annual rainfall and landslide from 2001 to 2017 in the TRW.
YearAnnual Rainfall (mm)Daily Rainfall (mm)Return Period (Years)Landslide Area (ha)Landslide Number
20013764545.510 to 20565.2324
2002187896.5<1.11375.3241
20031533233.01.11 to 2117.4194
20042593502.55 to 10154.4294
20054892545.510 to 20252.1375
20063774436.55 to 10108.0169
20073917267.01.11 to 2216.0255
20084112501.55 to 10153.5217
200934541206.0>2001392.91150
20102492379.52 to 5775.61091
20112225241.51.11 to 2422.4434
20123569433.02 to 5375.5513
20131736381.02 to 5388.3565
20142077172.51.11 to 2352.0461
20153050502.55 to 10240.8345
20163478491.55 to 10270.4407
20172428334.02 to 5312.6425
Note: The daily rainfall meant the maximum daily rainfall in the specific year, and the return period was based on the daily rainfall.
Table
Table 5. The characteristic and evolution comparison of large earthquake-induced and extreme rainfall-induced landslide events.
Table 5. The characteristic and evolution comparison of large earthquake-induced and extreme rainfall-induced landslide events.
Events2005 Kashimir Earthquake2008 Wenchuan Earthquake2009 Typhoon Morakot
Research areaNorthern PakistanSichuan, ChinaTRW in the study
ScaleML = 7.6ML = 7.9Heavy rainfall with return period greater than 200 years
Landslide triggersground shaking force to reduce the shear strength of slopeinfiltration water to reduce the shear strength of slope
Landslide typesdebris fall, debris flow, debris slide, rock falls (majority), rock topple and rock slides (majority) [6]Rock falls, rock slides, rock flows, debris falls, debris slides (majority), and debris flows [7,16,26]debris falls, debris slide (majority), and debris flows
Sediment deposition locationconcentrated in the down-hillslope to rivers near the fault linesconcentrated in the down-hillslope to rivers
Landslide volume (m3)2.16 × 109 [6]5–15 × 109 [39]6.65 × 106 [40] *
Distribution of extremely active and very active landslide
Elevation1000m to 1600 mlower than 2000 m1200 m to 1800 m
Slope30° to 40°30° to 50°30° to 40°
Fault or riversclose to the Muzaffarbad faultsclose to the Longmenshan faults-
Rivers-100 to 400 m to the rivers<300 m to the rivers
Long-term landslide evolution
Annual landslide area decline ratios5.02%9.56%7.45%
*: The landslide volume induced by 2009 Typhoon Morakot in the TRW was estimated by using the empirical equations [40].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, C.-H. Evaluating the Landslide Stability and Vegetation Recovery: Case Studies in the Tsengwen Reservoir Watershed in Taiwan. Water 2021, 13, 3479. https://doi.org/10.3390/w13243479

AMA Style

Wu C-H. Evaluating the Landslide Stability and Vegetation Recovery: Case Studies in the Tsengwen Reservoir Watershed in Taiwan. Water. 2021; 13(24):3479. https://doi.org/10.3390/w13243479

Chicago/Turabian Style

Wu, Chun-Hung. 2021. "Evaluating the Landslide Stability and Vegetation Recovery: Case Studies in the Tsengwen Reservoir Watershed in Taiwan" Water 13, no. 24: 3479. https://doi.org/10.3390/w13243479

APA Style

Wu, C. -H. (2021). Evaluating the Landslide Stability and Vegetation Recovery: Case Studies in the Tsengwen Reservoir Watershed in Taiwan. Water, 13(24), 3479. https://doi.org/10.3390/w13243479

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