Next Article in Journal
Testing Landscape, Climate and Lithology Impact on Carbon, Major and Trace Elements of the Lena River and Its Tributaries during a Spring Flood Period
Next Article in Special Issue
On Dam Failure Induced Seismic Signals Using Laboratory Tests and on Breach Morphology due to Overtopping by Modeling
Previous Article in Journal
MLE-Based Parameter Estimation for Four-Parameter Exponential Gamma Distribution and Asymptotic Variance of Its Quantiles
Previous Article in Special Issue
Tree-Ring Based Chronology of Landslides in the Shirakami Mountains, Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Hotspots and Decadal Evolution of Extreme Rainfall-Induced Landslides: Case Studies in Southern Taiwan

Department of Water Resources Engineering and Conservation, Feng Chia University, Taichung 40724, Taiwan
*
Author to whom correspondence should be addressed.
Water 2021, 13(15), 2090; https://doi.org/10.3390/w13152090
Submission received: 29 June 2021 / Revised: 24 July 2021 / Accepted: 27 July 2021 / Published: 30 July 2021
(This article belongs to the Special Issue Soil–Water Conservation, Erosion, and Landslide)

Abstract

:
The 2009 Typhoon Morakot triggered numerous landslides in southern Taiwan, and the landslide ratios in the Ailiao and Tamali river watershed were 7.6% and 10.7%, respectively. The sediment yields from the numerous landslides that were deposited in the gullies and narrow reaches upstream of Ailiao and Tamali river watersheds dominated the landslide recovery and evolution from 2010 to 2015. Rainfall records and annual landslide inventories from 2005 to 2015 were used to analyze the landslide evolution and identify the landslide hotspots. The landslide recovery time in the Ailiao and Tamali river watershed after 2009 Typhoon Morakot was estimated as 5 years after 2009 Typhoon Morakot. The landslide was easily induced, enlarged, or difficult to recover during the oscillating period, particularly in the sub-watersheds, with a landslide ratio > 4.4%. The return period threshold of rainfall-induced landslides during the landslide recovery period was <2 years, and the landslide types of the new or enlarged landslide were the bank-erosion landslide, headwater landslide, and the reoccurrence of old landslide. The landslide hotspot areas in the Ailiao and Tamali river watershed were 2.67–2.88 times larger after the 2009 Typhoon Morakot using the emerging hot spot analysis, and most of the new or enlarged landslide cases were identified into the oscillating or sporadic or consecutive landslide hotspots. The results can contribute to developing strategies of watershed management in watersheds with a dense landslide.

1. Introduction

Landslides induced by large earthquakes or extreme rainfall events have been the main reason for disasters in the past two decades in Taiwan. Typhoon Morakot in 2009 dumped around 2000 mm of rainfall in 3 days in southern Taiwan [1], resulting in severe landslide-related disasters, including the catastrophic deep-seated Xiaolin landslide [2] and the following dam failure [3]. Over a decade since the 2009 Typhoon Morakot, sediment-related disaster events still occurred in the Kaoping River watershed in southern Taiwan. Although most landslides in southern Taiwan had been gradually recovered, the hillslope was still under high landslide susceptibility.
The rate and location of landslide recovery after the large earthquake or extreme rainfall events play essential roles in developing the watershed management strategies for watersheds with a dense landslide. The landslide recovery in the watersheds with dense landslides after large earthquake events is related to the earthquake magnitude, geological settings, and fault distribution and characteristics [4,5,6], while recovery after extreme rainfall events were mostly related to the distribution of drainage network [7]. The sediment yield from landslides or debris flow in the watersheds with dense landslides is usually the dominant factor behind the geomorphologic evolution, particularly in the upstream watershed. The randomly deposited sediment in narrow upstream reaches usually results in rivers gradually becoming sinuous and inducing bank-erosion landslides. Sediment from bank-erosion landslides usually increases the sinuosity of narrow reaches and changes the geomorphology of the river in the upstream watershed.
Long-term geomorphologic landslide evolution in watersheds is strongly related to spatiotemporal landslide distribution [8], which can be observed using the spatiotemporal cluster analysis with the high-resolution digital elevation model (DEM) and multi-annual landslide inventories [4,6]. Several researchers have discussed the changes in the distribution and activeness of landslides after extreme rainfall-induced [7] or earthquake-induced [5,8,9,10,11] events and found that the spatiotemporal distribution and activeness of landslides were key factors behind the geomorphologic evolution of watersheds. Identifying landslide hotspots and cold spots using multi-annual landslide inventories can help researchers analyze landslide activeness and recovery after large earthquake-induced landslide disasters [12].
The space-time cluster analysis (abbreviated as spatiotemporal cluster analysis) in ArcGIS Pro software [13] is a useful analysis tool that can describe data’s spatial and temporal distribution patterns. This tool had been used to analyze the spread of the COVID-19 virus [14], road traffic accident occurrences [15], and the spread of air pollution [4,16] in recent years. Landslide disaster studies using the spatiotemporal cluster analysis have focused on discussing the long-term spatiotemporal distribution of disasters [5,8] and analyzing the relationship between disaster occurrence and related factors [6,9,17,18]. Spatiotemporal cluster analysis with multi-annual landslide inventories after extreme rainfall events can contribute to determining landslide hotspots and cold spots, identify locations where the landslide recovery was difficult, and analyze the reasons behind these factors. The use of spatiotemporal cluster analysis to observe landslide evolution trends and identify landslide clustering locations is more effective than only the spatial or temporal analysis of landslides.
The 2009 Typhoon Morakot (from 6–10 August 2009) caused the most severe rainfall-induced disaster event in the past two decades in Taiwan, and the return period accumulated 24 and 48 h of rainfall during the 2009 Typhoon Morakot in southern Taiwan exceeded 200 y [1]. The extreme rainfall event also caused numerous landslides and severe debris flow in southern Taiwan, and the landslide ratio (i.e., the ratio of the landslide area to watershed area) in the four sub-watersheds of the Kaoping River watershed after the typhoon exceeded 6.5% [1]. The geomorphologic evolution and developing trends of watersheds with dense landslides after 2009 Typhoon Morakot (abbreviated as after 2009) in southern Taiwan are worthy of discussion. The Ailiao river watershed (abbreviated as ARW) and Tamali river watershed (abbreviated as TRW) were the watersheds with the highest landslide ratio in southeastern and southwestern Taiwan after 2009. The ARW and TRW were selected to observe the landslide evolution from 2005 to 2015 and identify the landslide hotspots and cold spots using the spatiotemporal cluster analysis. The evolution characteristic of extreme rainfall-induced landslide events in Taiwan was also compared with that of large earthquake-induced landslide events in the world, and the cluster location and reason of new or enlarged landslides in the following years after 2009 were analyzed in the study.

2. Research Areas

2.1. Ailiao River Watershed (ARW)

The Ailiao river watershed (abbreviated as ARW) is located in southwestern Taiwan (Figure 1 and Figure 2), and the area is 623.3 km2. The average elevation and average slope in the ARW are 1006 m and 30.5°. The average annual precipitation is 3716 mm based on the records of six rainfall stations from 2005 to 2015 in the neighborhood of ARW (Figure 2a). The average precipitation in the rainy seasons, i.e., from May to October, occupies > 90% of the average annual precipitation. The land use distribution in the ARW is dominated by forest, which occupies 80.8% of the total watershed area. The main geological settings in the ARW consist of the Chaochou Formation, the Pilushan Formation, the Alluvium, and the Kaoling Schist (Figure 2b). The total precipitation during the 2009 Typhoon Morakot in the ARW was 2977 mm, i.e., around 80% of the average annual precipitation. The 2995 landslide cases (Figure 2a) were induced by the 2009 Typhoon Morakot in the ARW, and the landslide ratio, i.e., the ratio of the landslide area to the watershed area, was estimated as 7.6%. The landslides after 2009 centralized in the northeast ARW, especially in the A01 (8.2 km2), A02 (6.7 km2), A03 (2.5 km2), A07 (3.2 km2), and A11 (9.5 km2) sub-watersheds (Figure 2b). The occupied percentage of the landslide cases with area > 100,000 m2, 1000–100,000 m2, and <1000 m2 to all landslide cases in 2009 in the ARW were 3.5%, 73.0%, and 23.0%, respectively. The relation between the landslide length to width ratio and the mean slope in the ARW is shown in Figure 3; 93.1% and 57.6% of the landslide cases in 2009 in the ARW were of the landslide length to width ratio > 1.0 and ranged from 1.0 to 5.0. The rainfall-triggered slides, including the rotational and translational slides and flows on the hillslope with the slope > 30 degree, were the main landslide types in the ARW.

2.2. Taimali River Watershed (TRW)

The Taimali River Watershed (abbreviated as TRW) is located in southeastern Taiwan (Figure 1 and Figure 4) and the area is 264.9 km2. The average elevation and slope in the TRW are 789.4 m and 30.4°, respectively. The average annual precipitation is 2185 mm based on the records of five rainfall stations from 2005 to 2015 in the neighborhood of TRW (Figure 4a). The average precipitation in the rainy seasons, i.e., from May to October, occupies 76% of the average annual precipitation. The land use distribution in the TRW consists of forest (81.59%), agricultural land (9.12%), water conservancy (4.32%), and others. The main geological settings in the TRW consist of three main strata, including the Chaochou Formation, the Pilushan Formation, and the Alluvium (Figure 4b). The total precipitation during the 2009 Typhoon Morakot in the TRW was 932.5 mm, i.e., around 42.7% of the average annual precipitation. The 1283 landslide cases (Figure 4a) were induced by 2009.
Typhoon Morakot in the TRW, and the landslide ratio was estimated as 10.7%. The landslide after 2009 centralized in the upstream TRW, especially in the T01 sub-watershed (121.6 km2). The occupied percentage of the landslide cases with area > 100,000 m2, 1000–100,000 m2, and <1000 m2 to all landslide cases in 2009 in the TRW were 4.2%, 71.2%, and 24.6%, respectively. The relation between the landslide length to width ratio and the mean slope in the TRW is shown in Figure 3; 98.1% and 64.3% of the landslide cases in 2009 in the TRW were of the landslide length to width ratio > 1.0 and ranged from 1.0 to 5.0. These data show that the majority landslide type of the landslide cases induced by the 2009 Typhoon Morakot in the TRW were rainfall-triggered slides on the steep slope.

3. Data and Methods

3.1. Annual Landslide Inventories

The annual landslide inventories from 2005 to 2015 produced by the Forestry Bureau in Taiwan were used in this study, and the minimum landslide area in the annual landslide inventories was 100 m2. Based on Varnes’ classification [19], the majority of the landslide cases induced by 2009 Typhoon Morakot in southern Taiwan were rotational slides, translational slides, and flows [20,21,22].

3.2. Effective Accumulated Rainfall Index (EAR)

The EAR index (unit: mm) was used to assess the landslide-induced strength of rainfall events. The EAR index, defined in Equation (1), is the summation of daily rainfall on the assessment day (Rt) and the 7-day antecedent rainfall before the assessment day. The K coefficient, representing the decay constant, was set to 0.7 based on Taiwanese landslide research [23]. Equation (1) is calculated as follows:
E A R t = i = 0 7 R t × K i
The rainfall records used to estimate the EAR index value were collected from the representative rainfall stations at watersheds. For inclusion, the representative rainfall stations had to be located within the watershed, and the rainfall records from 2005 to 2015 had to be available without any missing data. The representative rainfall stations in the two watersheds are Ali station in the ARW and Jinfong station in the TRW. The annual landslide inventories were used in this study. It is challenging to find data on the time and date of landslide occurrences and estimate the rainfall threshold to induce the landslide. Rather than determining the precise time and date of landslide occurrences, the EAR values used in this study serve as reference coefficients to understand the intensity of landslides induced by typhoons and other heavy rainfall events each year from 2005 to 2015.

3.3. Landslide Topographic Position Analysis

The topographic position analysis method can be used to explain the main inducing factors of landslides [24]. Three parameters of the landslide on the hillslope, including the distance between the ridge and the crown of the landslide (DP), the distance between the stream and the toe of the landslide (DB), and the distance between the ridge of the hillslope and river (DH), are used to explain the relative location of the landslide in the hillslope. The bubble plot is frequently used to draw the result of the topographic position analysis using the normalized distance from a landslide to the ridge (DP/DH) as the X-axis, the normalized distance from a landslide to the stream (DB/DH) as Y-axis, and the size of the bubble as the landslide area. If the bubbles are located in the upper-left portion of the bubble plot (DP/DH < 0.5 and DB/DH > 0.5), the landslide cases are located near the ridge and possibly induced by earthquake events [24]. If the bubbles are located in the lower-left portion of the bubble plot (DP/DH > 0.5 and DB/DH < 0.5), the landslide cases are located near the stream and possibly induced by rainfall or flooding events [1,21].

3.4. Spatiotemporal Cluster Analysis Method

We used the emerging hot spot analysis in the space-time cluster analysis tool in the ArcGIS Pro software to analyze the landslide evolution and identify the landslide hotspots and cold spots from 2005 to 2015. The emerging hot spot analysis tool can detect eight hotspot or cold spot trends, and the definition of the eight hot spot or cold spot trends had been described in Table 1 (revised from [14]). The emerging hot spot analysis was widely used in observing the evolution of the natural or artificial phenomenon but has still rarely been used to analyze the landslide evolution. The analysis unit in the study is a 5 m × 5 m grid, and the time step is a year. The clustering intensity of landslide in each analysis unit was estimated using the Getis-Ord Gi statistic [25], which considered the clustering intensity value for each analysis unit within the context of the values for the neighboring analysis unit. In the study, the neighborhood distance of the analysis unit was set as 25 m.

4. Decadal Analyses Results

4.1. Rainfall Distribution and Landslide Ratio

The EAR distributions from 2005 to 2015 in the two watersheds are shown in Figure 5 and Table 2. The average EAR value from 2005 to 2015 was 39.7 in the ARW and 29.7 in the TRW. The highest EAR values (EARh) in the ARW and TRW were 1926.9 and 1123.5 on 8 August 2009. The EARh and EARa (the average of the three highest EAR values in each year) from 2005 to 2008 in the two watersheds were larger than those from 2010 to 2015. The return periods of the top ten daily rainfall events from 2005 to 2008 in the two watersheds were estimated to be 10–50 years, and those from 2010 to 2015 were estimated to only be <2 year. The data demonstrated that the EAR value and the return periods of rainfall events from 2005 to 2008 in the two watersheds were larger than those from 2010 to 2015.
The landslide ratios in 2009 in the ARW and TRW (Figure 5) were their historical peaks. The average landslide ratios in the ARW and TRW from 2005 to 2008 were 1.6% and 2.2%, respectively, and those from 2010 to 2015 were 4.3% and 5.9%, respectively. The trends of the landslide ratios in the two watersheds after 2009 were oscillating instead of stably decaying. The EARh and EARa in 2011 and 2013 in the ARW were smaller than those from 2005 to 2007, but the landslide ratio increased in 2011 and 2013. Other similar examples are shown in comparing the EARh, EARa, and landslide ratios in 2011 and 2013 in the TRW. This data implied that landslides were more easily induced after the 2009 Typhoon Morakot. The rainfall factor was possibly not the only landslide-inducing factor in the two watersheds after 2009.

4.2. Landslide Statistical Data

The research period was divided into three periods (i.e., 2005–2008, 2009, and 2010–2015) to analyze the changes in landslide distribution before and after 2009. The landslides’ statistical data from 2005 to 2015 in the two watersheds are shown in Figure 6 and Table 3. The area and number of landslides from 2010 to 2015 in the two watersheds were larger than those from 2005 to 2008. From 2005 to 2015, the year with the most landslides was 2009, but the year with the most landslide numbers was 2013. In the ARW, for example, the landslide area in 2013 was 42% smaller than that in 2009, but the landslide number in 2013 was 31% higher than that in 2009. The same trend was observed in 2013 in the TRW. This data implies that most of the landslides induced by 2009 Typhoon Morakot gradually recovered, but some new landslides occurred in the two watersheds in 2013.
This study analyzed the landslide distribution at the sub-watershed scale to find the sub-watersheds in which landslides were induced in the years following the 2009 Typhoon Morakot. The landslide evolution trend index (abbreviated as LET) in this study was defined as the average change ratio of the landslide area from 2010 to 2015, and the LET was estimated in each sub-watershed of the two watersheds (Figure 7). A negative LET value indicates that the total landslide area in this sub-watershed gradually decreases, while a positive LET value indicates that the total landslide area gradually increases. The average LET value in the sub-watersheds was −0.022 and −0.072 km2/year in the ARW and TRW.
The sub-watersheds with positive LET values were located upstream of ARW and TRW. There were 13 and 2 sub-watersheds with the positive LET values in the ARW and TRW, respectively, and the landslide ratio of the 15 sub-watersheds after 2009 was larger than 4.4%. There were six sub-watersheds with the LET values > 0.05, including A01, A02, A03, A07, and A11 in the ARW and T01 in the TRW, and the landslide ratio of the six sub-watersheds after 2009 was greater than 12.1%. The watershed areas in the A01, A02, A03, A07, A11, and T01 sub-watersheds were 8.2, 6.7, 24.7, 31.7, 9.5, and 121.6 km2, respectively, and the landslide ratios after the 2009 Typhoon Morakot were 27.8%, 21.2%, 26.2%, 21.5%, 12.1%, and 20.7%, respectively. These results imply that the landslides in the sub-watersheds with a landslide ratio of >4.4% after 2009 in the ARW and TRW were difficult to recover and were easily induced or re-induced from 2010 to 2015.

4.3. Landslide Topographic Position Analysis

The study used the landslide topographic position analysis to examine the landslide evolution before and after 2009 in the ARW and TRW. The A03 (LET = 0.32 km2/y), A31 (LET = −0.31 km2/y), and T01 (LET = 0.43 km2/y) sub-watersheds were selected for comparison of landslide evolution before and after 2009 (Table 4 and Figure 8 and Figure 9). The area in the A31 sub-watershed was 33.9 km2, and the landslide area and landslide ratio in 2009 in the A31 sub-watershed were 2.8 km2 and 8.3%. The ratio of landslide area from 2009 to 2015 in the upslope, mid-slope, and downslope were 19.4%, 25.5%, and 38.2%, respectively, in the ARW and 27.6%, 29.8%, and 31.1% in the TRW, respectively. The landslide located in the downslope was the most difficult to recover from 2009 to 2015 in the slope.
A similar trend was also found in the A03, A31, and T01 sub-watersheds. The ratio of landslide area from 2009 to 2015 in the downslope was 79.7%, 23.0%, and 68.5% in the A03, A31, and T01 sub-watersheds, respectively. Figure 8 and Figure 9 show that a reduction was observed in the number of upslope, mid-slope, and downslope landslides in the sub-watersheds, but the landslides in 2015 were concentrated in the downslope area. From 2009 to 2015, a large cluster of small-area landslides occurred downslope in the sub-watersheds, with poor recovery. Most of the landslides in the A03 and T01 sub-watersheds in 2015 were centered in areas with a normalized distance to a ridge of >0.7, meaning that the inducing factors should be related to the bank-erosion landslide, which was possibly induced by the sinuous rivers with huge amounts of deposited sediment.

4.4. Spatiotemporal Landslide Hotspot Analysis

The landslide ratios in the ARW and TRW after 2009 were 7.6% and 10.7%, and those were the top two highest landslide ratios in the watershed scale in Taiwan. It is interesting to understand the evolution of numerous landslides and compare the characteristic of landslide distribution before and after 2009 in the two watersheds. The evolutions of the landslide from 2005 to 2015 in the ARW and TRW were observed from the spatiotemporal landslide hotspot analyses (Table 5 and Figure 10) in the study.
The total areas from 2010 to 2015 in the two watersheds are 1.15–2.23 times larger than those from 2005 to 2008, and the increases in the landslide hot spot areas from before to after 2009 in the two watersheds were evident. The landslide hot spot areas from 2010 to 2015 in the two watersheds are 2.67–2.88 times larger than those from 2005 to 2008, and the landslide cold spot area is 1.73–1.93 times larger. This result means that the total time of areas identified as landslides from 2010 to 2015 is longer than that from 2005 to 2008. The landslide recovery was more difficult, and the landslide was easier to be clustered after than before 2009 Typhoon Morakot.
The no pattern detected area means that the time of area identified as a landslide is shorter than 90% of the research period (Table 1). The occupied percentages of the no pattern detected areas from 2005 to 2015 in the two watersheds are 53.0–56.5%, but those from 2005 to 2008 and from 2010 to 2015 are only 15.9% to 20.2%. This data means that 36.3–37.1% of landslide areas in the two watersheds recovered to the non-landslide areas in 4 to 9 years.
The landslide hot spots are centralized in OHS, SHS, and CHS in each research period, while the landslide cold spots are centralized in OCS, CCS, and SCS. The landslide hot spots and cold spots were reclassified into the main hot spots, the main cold spots, no pattern detected, and others. The main hot spots included OHS, SHS, and CHS, the main cold spots included OCS, CCS, and SCS, and the others included all the other hot spots and cold spots except the main hot spots and cold spots. After 2009, the main hot spots constituted 34.0–41.9% of all hot spots, whereas the main cold spots accounted for 31.6–37.8% of all cold spots.
Figure 11 and Figure 12 present the main hot spots and cold spots from 2005 to 2015 in the two watersheds. The upstream sub-watersheds with dense landslide distributions were the main hot spot cluster areas in the two watersheds.
The main hot spots from 2005 to 2008 were discretely distributed in the upstream sub-watersheds of ARW and TRW, and those from 2010 to 2015 were densely clustered in the upstream of ARW and TRW, particularly in the A01 and T01 sub-watersheds.
Obvious increases in the average landslide ratios from after to before the 2009 Typhoon Morakot in the two watersheds were noted. The CHS were the hot spots that exhibited the largest area expansion from after to before the 2009 Typhoon Morakot, and the OCS were the cold spots that exhibited the largest area reduction. The CHS percentage increased by 7.5% to 16.3% from after to before the 2009 Typhoon Morakot, and the OCS percentage decreased from 11.4% to 21.8%. This means that the recovery of landslides induced by 2009 Typhoon Morakot was slower than that before 2009.
The A01 sub-watershed was selected as the representative sub-watershed to explain the distribution of the main hot spots and cold spots in the study. The strata in A01 comprise the Pilushan and Chaochou formations (62.6% and 37.4%, respectively) from the Eocene epoch and Middle Miocene sub-epoch, respectively. The lithology of the Pilushan formation comprises slate with metasandstone and igneous rock, whereas that of the Chaochou formation is argillite and slate with an alternation of metasandstone or argillite. The main hot spots and main cold spots in the A01 sub-watershed increased substantially after Typhoon Morakot. From 2005 to 2008, 2010 to 2015, and 2005 to 2015, the main hot spots in the A01 sub-watershed constituted 3.0%, 17.0%, and 12.5%, respectively, and the main cold spots constituted 3.8%, 5.9%, and 4.2%, respectively. The main hotspots from 2010 to 2015 in the A01 sub-watershed were concentrated in the headwater landslides, bank-erosion landslides in sinuous reaches, and reoccurrence of older (from 2005 to 2008) landslides.
Mechanisms and triggering factors of landslide events, landslide areas with poor recovery, and geomorphological evolution trends can be explained, located, and predicted using the distributions of landslide hot spots and cold spots that were constructed through spatiotemporal analysis. The results of the spatiotemporal analysis over the various periods have different implications. Specifically, the hot spot and cold spot distributions from 2005 to 2015, from 2005 to 2008, and from 2010 to 2015 in the ARW and TRW represent the long-term landslide evolution.

5. Discussion

The prediction of landslide recovery in watersheds with dense landslides could be the key factor for watershed management. The characteristic of landslide recovery in the following years after the large earthquake or extreme rainfall events are worth comparing and discussing. We explained the recovery characteristic of extreme rainfall-induced landslides by comparing the landslide recovery conditions after the 2005 Kashmir earthquake [5], the 2008 Wenchuan earthquake [4,9], and the 2009 Typhoon Morakot in this study. The time, location, and rate of landslide recovery after the large earthquake or extreme rainfall events are the key discussion topics in this study.
The oscillating period was observed after the large earthquake or extreme rainfall events based on the annual landslide area data. The oscillating period can be defined as that the annual landslide area and landslide number in this period is an oscillating trend instead of a stable decline trend. The oscillating period for the serious earthquake-induced landslide events ranged from 3 to 5 years. The extreme rainfall-induced landslide events in the study were estimated as 5 years (Figure 4 and Table 6, from 2010 to 2014). The landslide in the watersheds in the oscillating period was active and easily induced, re-induced, or enlarged. The average annual landslide area decline rates (abbreviated as LAD) after 2014 were larger than that during the oscillating period (from 2010 to 2014), and the average LAD during or after the oscillating period in this study was also larger than those from the large earthquake-induced landslide events. This means that the recovery rate of the extreme rainfall-induced landslide was faster than that of large earthquake-induced landslide.
The location and reason of new or enlarged landslides after the large earthquake or extreme rainfall events are worth discussing and comparing. The new or enlarged landslides in the following years after the 2005 Kashimir earthquake (including the active, very active, and extremely active landslides in [5]) were mostly located along the Muzaffarabad fault or in the high fractured and jointed rocks areas, or along with the drainage network, or in the source of the river and large landslide. Moreover, the new or enlarged landslides in the following years after the 2008 Wenchuan earthquake (the active landslides in [4]) were located in deep gullies, the source of debris flow and large landslides. Three factors, including the geological setting, the drainage network, and the landslide area, dominate the rate of landslide recovery after 2009 in the ARW and TRW in the study.
The statistical data and distribution of landslide evolution in the ARW and TRW are shown in Table 7 and Figure 13. The new or enlarged landslide in the following years after 2009 centralized in the northeast ARW and upstream TRW. The strata in the northeast ARW comprise 62.6% Pilushan formation (metasandstone and igneous rock) and 37.4% Chaochou formations (argillite and slate with an alternation of metasandstone or argillite), and three faults and anticlines also pass through the northeast ARW. The strata in the upstream TRW comprise the Chaochou formations (sandstone), kaolinite schist, and Pilushan formations (metasandstone and igneous rock), and three faults and anticlines also pass through the northeast TRW. Fractured slate, sandstone, or argillite are the main geological composition in the northeast ARW and upstream TRW, and also explain the reasons for the centralization of new or enlarged landslides in this area.
The landslide evolution results from the comparison of landslide inventories in two different years can be classified into three types, including recovered landslides, not recovered landslides, and new or enlarged landslides (Figure 13). The recovered landslide area from the comparison between 2009 and 2010 (Table 7) was the area identified as landslide in 2009 but not in 2010, and the not recovered landslide area were the areas identified as landslide in 2009 and 2010. The new or enlarged landslide area was the area identified as landslide in 2010 but not in 2009. The new or enlarged landslide in the ARW and TRW also centralized along with the drainage network, particularly in the upstream watersheds. Hugh sediment yield from the landslide in the upstream watershed with dense landslide should be the main reason. The landslide volume was estimated the empirical equations from Taiwan [26] for the landslide area < 106 m2 and Italy [27] for the landslide area ≧ 106 m2]. The landslide volume induced by 2009 Typhoon Morakot was estimated as 65.0 × 106 m3 and 224.5 × 106 m3 in the ARW and TRW. There were 848 landslide cases after 2009 in the northeast upstream ARW, including A01, A02, A03, A07, and A11 sub-watersheds, and 1138 landslide cases in the upstream TRW, i.e., the T01 sub-watershed. The landslide volume was estimated as 5.2 × 106 m3 in the northeast upstream of the ARW and 223.9 × 106 m3 in the T01 sub-watershed. Huge sediment was yielded, deposited in the narrow reaches, and dominated the evolution of landslide and river geomorphology in the northeast ARW and upstream TRW.
Huge sediment in the upstream watershed was continuously transported into the gullies and rivers and also resulted in the frequent occurrence of new or enlarged landslides in the neighborhood of gullies and rivers from 2010 to 2015 in the ARW. Moreover, 51.3%, 54.0%, and 58.2% of the landslide areas after 2009 in the ARW had been recovered in 2010, 2013, and 2015, respectively. The new or enlarged landslide area from 2010 to 2015 in the ARW showed a continuously increasing trend. The occupied percentage of a new or enlarged landslides located in the neighborhood of gullies from 2010 to 2015 was 53.9–56.1%, and the area of a new or enlarged landslide located in the neighborhood of the river from 2010 to 2015 also showed an increasing trend.
The centralization of new or enlarged landslides in the neighborhood of large landslides was mentioned [4,5], and it was also observed in the study. The occupied percentage of a new or enlarged landslide located in the neighborhood of large landslide cases from 2010 to 2015 in the ARW was 31.2–35.3%, particularly in the northeast ARW.
The dentification of landslide hot spots using the emerging hot spot analysis in this study can show the clustering strength of old, new, and enlarged landslides in space and time and provide a potential landslide location. The advantage of the identification of landslide hot spots using the emerging hot spot analysis is that we can estimate the maintenance time of landslides from the classification of landslide hot spots, and this information also contributes to making the priority of watershed management measures in the watersheds with dense landslides. The management strategy for the watersheds with huge sediment yield should be implemented considering the landslide evolution trend. The landslide evolution cases in the ARW and TRW in Taiwan demonstrated that controlling the sediment in the drainage network and the landslide boundary should be the priority after the extreme rainfall-induced landslide events.

6. Conclusions

This study used the rainfall analysis, spatiotemporal landslide hotspot analyses, and comparison analysis of large earthquake- and extreme rainfall-induced landslide evolution to understand the characteristic of rainfall-induced landslide evolution, which was useful in assessing the landslide activeness after an extreme rainfall event. We used the EAR to assess the landslide-induced strength of rainfall events from 2005–2015, and the EAR values in the ARW and TRW were larger than before after the 2009 Typhoon Morakot. The landslide evolution trend index (LET) was used to assess the recovery ratio of landslide area after 2009, and the LET value in most of the sub-watersheds in the ARW and TRW were ranged 0.022–0.072 km2/year. However, some sub-watersheds in the ARW and TRW, particularly in the upstream watershed with the landslide ratio > 4.4%, were still of LET value > 0.05 km2/year after 2009. The landslides downslope of sub-watersheds with positive LET values in the ARW and TRW after 2009 were easily induced, re-induced, or enlarged and difficult to recover based on the landslide topographic position analysis. Most of the new or enlarged landslides in the ARW and TRW after 2009 were classified into oscillating or sporadic or consecutive landslide hotspots and centralized along with the drainage network or large landslide boundary. The watersheds with dense landslides needed to spend 3–5 years, i.e., the oscillating period in the study, to achieve the stable landslide recovery based on the comparison of landslide recovery after the large earthquake or extreme rainfall events. The landslide area decline rates in the ARW and TRW after 2009 were 1.6–2.5 times larger after than during the oscillating period. The new or enlarged landslides after 2009 in the ARW and TRW was centralized in the huge sediment-deposited, narrow, and sinuous reaches or the boundary of a large landslide in the upstream watersheds with a geological composition of fractured slate, sandstone, or argillite. The findings from the study point out that the watershed management strategies in the watershed with dense landslides after the extreme rainfall-induced landslide events should be emphasized to control the huge sediment yield from the numerous landslides, particularly in the upstream watersheds.

Author Contributions

Conceptualization and methodology, C.W. and C.L.; software, C.L.; writing, C.W. and C.L.; supervision and funding acquisition, C.W. Both authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Ministry of Science and Technology, Taiwan; Grant number: MOST 108-2625-M-035-003-(Taiwan).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Financial supports from the Ministry of Science and Technology of Taiwan (R.O.C.) under contract MOST 108-2625-M-035-003- are appreciated.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. 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]
  2. Wu, C.-H.; Chen, S.-C.; Feng, Z.-Y. Formation, failure, and consequences of the Xiaolin landslide dam, triggered by extreme rainfall from Typhoon Morakot, Taiwan. Landslides 2013, 11, 357–367. [Google Scholar] [CrossRef] [Green Version]
  3. Feng, Z.-Y. The seismic signatures of the surge wave from the 2009 Xiaolin landslide-dam breach in Taiwan. Hydrol. Process. 2011, 26, 1342–1351. [Google Scholar] [CrossRef]
  4. Li, J.; Wang, N.; Wang, J.; Li, H. Spatiotemporal evolution of the remotely sensed global continental PM2.5 concentration from 2000–2014 based on Bayesian statistics. Environ. Pollut. 2018, 238, 471–481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Shafique, M. Spatial and temporal evolution of co-seismic landslides after the 2005 Kashmir earthquake. Geomorphology 2020, 362, 107228. [Google Scholar] [CrossRef]
  6. Yunus, A.P.; Fan, X.; Tang, X.; Jie, D.; Xu, Q.; Huang, R. Decadal Vegetation Succession from MODIS Reveals the Spatio-Temporal Evolution of Post-Seismic Landsliding after the 2008 Wenchuan Earthquake. Remote Sens. Environ. 2020, 236, 111476. [Google Scholar] [CrossRef]
  7. Wu, C. Landslide Susceptibility Based on Extreme Rainfall-Induced Landslide Inventories and the Following Landslide Evolution. Water 2019, 11, 2609. [Google Scholar] [CrossRef] [Green Version]
  8. Valenzuela, P.; Domínguez-Cuesta, M.J.; Mora García, M.A.; Jiménez-Sánchez, M. A Spatial-Temporal Landslide Inventory for the NW of Spain: BAPA Database. Geomorphology 2017, 293, 11–23. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, S.H.; Lin, C.W.; Tseng, C.M. A Statisticalmodel for the Impact of the 1999 Chi-Chi Earthquake on the Subsequent Rain-fall-Induced Landslides. Eng. Geol. 2013, 156, 11–19. [Google Scholar] [CrossRef]
  10. Tang, C.; Jiang, Z.; Li, W. Seismic Landslide Evolution and Debris Flow Development: A Case Study in the Hongchun Catchment, Wenchuan Area of China. Eng. Geol. Soc. Territ. 2015, 2, 445–449. [Google Scholar] [CrossRef]
  11. Zhang, S.; Zhang, L. Impact of the 2008 Wenchuan earthquake in China on subsequent long-term debris flow activities in the epicentral area. Geomorphology 2017, 276, 86–103. [Google Scholar] [CrossRef]
  12. 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]
  13. ESRI. [WWW Document]. Doc. ArcGIS. 2018. Available online: https://pro.arcgis.com/en/pro-app/tool-reference/analysis/enrich-layer.htm (accessed on 15 May 2021).
  14. Barboza, G.E.; Schiamberg, L.B.; Pachl, L. A spatiotemporal analysis of the impact of COVID-19 on child abuse and neglect in the city of Los Angeles, California. Child. Abus. Negl. 2021, 116, 104740. [Google Scholar] [CrossRef]
  15. Cheng, Z.; Zu, Z.; Lu, J. Traffic Crash Evolution Characteristic Analysis and Spatiotemporal Hotspot Identification of Urban Road Intersections. Sustainability 2018, 11, 160. [Google Scholar] [CrossRef] [Green Version]
  16. Nielsen, C.; Amrhein, C.G.; Shah, P.S.; Stieb, D.M.; Osornio-Vargas, A.R. Space-time hot spots of critically ill small for gestational age newborns and industrial air pollutants in major metropolitan areas of Canada. Environ. Res. 2020, 186, 109472. [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. Shou, K.-J.; Lin, J.-F. Evaluation of the extreme rainfall predictions and their impact on landslide susceptibility in a sub-catchment scale. Eng. Geol. 2020, 265, 105434. [Google Scholar] [CrossRef]
  19. Varnes, D.J. Introduction to Landslides: Analysis and Control; Transportation and Road Research Board; National Academy of Science: Washington, DC, USA, 1978; pp. 11–33. [Google Scholar]
  20. Lu, S.Y.; Lin, C.Y.; Hwang, L.S. Spatial Relationships between Landslides and Topographical Factors at the Liukuei Experi-mental Forest, Southwestern Taiwan after Typhoon Morakot. Taiwan J. For. Sci. 2011, 26, 399–408. [Google Scholar]
  21. Wu, C.Y.; Tsai, C.W.; Chen, S.C. Topographic Characteristic Analysis of Landslides in Kaoping River Watershed. J. Chin. Soil Water Conserv. 2016, 47, 156–164. (In Chinese) [Google Scholar]
  22. Lin, Y.T.; Chang, K.C.; Yang, C.J. Object-Based Classification for Detecting Landslides and Vegetation Recovery—A Case at Baolai, Kaohsiung. J. Chin. Soil Water Conserv. 2018, 49, 98–109. (In Chinese) [Google Scholar]
  23. Chen, C.-W.; Saito, H.; Oguchi, T. Analyzing rainfall-induced mass movements in Taiwan using the soil water index. Landslides 2016, 14, 1031–1041. [Google Scholar] [CrossRef]
  24. Meunier, P.; Hovius, N.; Haines, J.A. Topographic site effects and the location of earthquake induced landslides. Earth Planet. Sci. Lett. 2008, 275, 221–232. [Google Scholar] [CrossRef]
  25. Getis, A.; Ord, J.K. The Analysis of Spatial Association by Use of Distance Statistics. Geogr. Anal. 2010, 24, 189–206. [Google Scholar] [CrossRef]
  26. 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]
  27. Guzzetti, F.; Ardizzone, F.; Cardinali, M.; Galli, M.; Reichenbach, P.; Rossi, M. Distribution of landslides in the Upper Tiber River basin, central Italy. Geomorphology 2008, 96, 105–122. [Google Scholar] [CrossRef]
Figure 1. Location of Taiwan, Ailiao river watershed (abbreviated as ARW), and Taimali river watershed (abbreviated as TRW).
Figure 1. Location of Taiwan, Ailiao river watershed (abbreviated as ARW), and Taimali river watershed (abbreviated as TRW).
Water 13 02090 g001
Figure 2. The distribution of elevation, landslide after 2009 Typhoon Morakot, and sub-watersheds (a), geological settings (b), in the ARW.
Figure 2. The distribution of elevation, landslide after 2009 Typhoon Morakot, and sub-watersheds (a), geological settings (b), in the ARW.
Water 13 02090 g002
Figure 3. Relationship between the ratio of landslide length to width and mean slope of the landslide cases induced by 2009 Typhoon Morakot in the ARW and TRW.
Figure 3. Relationship between the ratio of landslide length to width and mean slope of the landslide cases induced by 2009 Typhoon Morakot in the ARW and TRW.
Water 13 02090 g003
Figure 4. The distribution of elevation, landslide after 2009 Typhoon Morakot, and sub-watersheds (a), geological settings (b) in the TRW.
Figure 4. The distribution of elevation, landslide after 2009 Typhoon Morakot, and sub-watersheds (a), geological settings (b) in the TRW.
Water 13 02090 g004
Figure 5. The distribution of effective accumulated rainfall index (abbreviated as EAR) value (black bar) and landslide ratio (dash line) from 2005 to 2015 in the ARW (up figure) and TRW (down figure).
Figure 5. The distribution of effective accumulated rainfall index (abbreviated as EAR) value (black bar) and landslide ratio (dash line) from 2005 to 2015 in the ARW (up figure) and TRW (down figure).
Water 13 02090 g005
Figure 6. The area (solid line) and number (dash line) of landslide in the ARW (black) and TRW (red) from 2005 to 2015.
Figure 6. The area (solid line) and number (dash line) of landslide in the ARW (black) and TRW (red) from 2005 to 2015.
Water 13 02090 g006
Figure 7. The average landslide density and landslide evolution trend index value from 2005 to 2015 in the ARW (A) and TRW (B).
Figure 7. The average landslide density and landslide evolution trend index value from 2005 to 2015 in the ARW (A) and TRW (B).
Water 13 02090 g007
Figure 8. The topographic position analysis of landslide in 2008 (ad), 2009 (be), and 2015 (cf) in the A03 (up figures) and A31 (down figures) subwatersheds in the ARW.
Figure 8. The topographic position analysis of landslide in 2008 (ad), 2009 (be), and 2015 (cf) in the A03 (up figures) and A31 (down figures) subwatersheds in the ARW.
Water 13 02090 g008
Figure 9. The topographic position analysis of landslide in 2008 (a), 2009 (b), and 2015 (c) in the T01 sub-watershed in the TRW.
Figure 9. The topographic position analysis of landslide in 2008 (a), 2009 (b), and 2015 (c) in the T01 sub-watershed in the TRW.
Water 13 02090 g009
Figure 10. The occupied percentage of landslide hot spots and cold spots from 2005 to 2015 (black line), from 2005 to 2008 (blue line), from 2010 to 2015 (red line), and in the ARW (a) and TRW (b).
Figure 10. The occupied percentage of landslide hot spots and cold spots from 2005 to 2015 (black line), from 2005 to 2008 (blue line), from 2010 to 2015 (red line), and in the ARW (a) and TRW (b).
Water 13 02090 g010
Figure 11. The main landslide hot spot and cold spot from 2005 to 2015 (a), from 2005 to 2008 (b), and from 2010 to 2015 (c) in the ARW.
Figure 11. The main landslide hot spot and cold spot from 2005 to 2015 (a), from 2005 to 2008 (b), and from 2010 to 2015 (c) in the ARW.
Water 13 02090 g011
Figure 12. The main landslide hot spot and cold spot from 2005 to 2015 (a), from 2005 to 2008 (b), and from 2010 to 2015 (c) in the TRW.
Figure 12. The main landslide hot spot and cold spot from 2005 to 2015 (a), from 2005 to 2008 (b), and from 2010 to 2015 (c) in the TRW.
Water 13 02090 g012
Figure 13. The landslide evolution from 2009 to 2010 (a), to 2013 (b), and to 2015 (c) in the ARW. The left down plot in each figure is the landslide evolution in the northeast ARW, including the A01, A02, A03, A07, and A11 sub-watersheds.
Figure 13. The landslide evolution from 2009 to 2010 (a), to 2013 (b), and to 2015 (c) in the ARW. The left down plot in each figure is the landslide evolution in the northeast ARW, including the A01, A02, A03, A07, and A11 sub-watersheds.
Water 13 02090 g013
Table 1. The classifications and definition of emerging landslide hot spot and cold spot in the study.
Table 1. The classifications and definition of emerging landslide hot spot and cold spot in the study.
ClassificationDefinition
Consecutive
(CHS or CCS)
A landslide location with a single uninterrupted run of statistically significant hot spot or cold spot areas in the final year during the research time period. The landslide location has never been a statistically significant hot spot or cold spot before the final hot spot or cold spot run.
Diminishing
(DHS or DCS)
A landslide location that has been a statistically significant hot spot or cold spot for 90% of the research time period, including the final year. In addition, the clustering intensity of landslide in each year is decreasing (increasing) overall and that decrease (increase) is statistically significant.
Historical
(HHS or HCS)
The most recent year is not hot spot or cold spot, but at least 90% of the research time period has been a statistically significant hot spot or cold spot.
Intensifying
(IHS or ICS)
A landslide location that has been a statistically significant hot spot or cold spot 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 location identified as a statistically significant hot spot or cold spot since the first year of the research time period but was not previously identified as a statistically significant hot spot or cold spot.
Oscillating
(OHS or OCS)
A statistically significant hot spot or cold spot for the final year that has a history of also being a statistically significant cold spot or hot spot during a prior year. Less than 90% of the research time period have been statistically significant hot spot or cold spot.
Persistent
(PHS or PCS)
A landslide location that has been a statistically significant hot spot or cold spot 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 location that is an on-again then off-again hotspot or cold spot. Less than 90% of the research time period have been statistically significant hot spot or cold spot, and none of the time-step intervals have been statistically significant colds pot or hot spot.
No pattern detected (No)The analysis area does not fit any definition of hot spot or cold spot classifications
Note: The CHS and CCS are the abbreviations of consecutive hot spot and consecutive cold spot. The regulation of abbreviation is applied to each hot spot and cold spot in the study.
Table 2. The statistical data of the EAR values from 2005 to 2015 in the two watersheds.
Table 2. The statistical data of the EAR values from 2005 to 2015 in the two watersheds.
Year20052006200720082009201020112012201320142015
in the ARW
EARh1481.7784.8707.2631.91926.9894.5163.51104.4658.5438.5609.9
EARa1353.7700.6680.9616.51755.8697.8150.1943.6587.4340.7551.1
in the TRW
EARh997.0663.4723.1485.21123.5759.3374.6766.2576.2304.8367.3
EARa856.3620.6593.6449.6944.0588.9332.0670.3494.2250.8351.0
Note: The EARh means the highest EAR value, and the EARa means the average of the three highest EAR values in each year.
Table 3. The average area and number of landslides in the two watersheds.
Table 3. The average area and number of landslides in the two watersheds.
Year2005 to 20152005 to 200820092010 to 2015
ARWAverage landslide area (km2)22.810.148.427.1
Average landslide number1902.81132.323552341.2
TRWAverage landslide area (km2)13.35.828.415.8
Average landslide number766.3482.01100900.2
Table 4. The topographic position analysis results in the ARW and TRW.
Table 4. The topographic position analysis results in the ARW and TRW.
TimeB2009A2015B2009A2015B2009A2015
ARWA03 sub-watershedA31 sub-watershed
UA4.6413.563.692.630.130.290.280.210.010.130.040.02
MA0.812.392.200.610.765.272.881.920.151.310.480.21
DA2.4912.165.374.650.360.740.640.590.241.220.50.28
TRWT01 sub-watershed
UA2.054.783.371.320.090.230.20.15
MA0.383.191.090.950.454.760.030.01
DA1.595.922.561.840.371.110.90.76
Note: “B” and “A” mean that the average landslide area before 2009, i.e., from 2005 to 2008 and after 2009, i.e., from 2010 to 2015. UA, MA, and DA mean the upslope, mid-slope, and downslope landslide area (km2).
Table 5. The statistical data of spatiotemporal landslide hot spots and cold spots in the ARW and TRW.
Table 5. The statistical data of spatiotemporal landslide hot spots and cold spots in the ARW and TRW.
WatershedARW TRW
year05–1505–0810–1505–1505–0810–15
HS (km2)17.67.822.58.84.812.8
CS (km2)13.510.320.36.85.59.5
No (km2)35.14.68.120.32.34.9
Total (km2)66.222.850.935.912.627.1
Note: The 05–15, 05–08, and 10–15 mean from 2005 to 2015, from 2005 to 2008, and from 2010 to 2015. The HS and CS mean the total area of all hot spots and cold spots, and the NO means the no pattern detected area.
Table 6. Comparison of the average annual landslide area decline rate from the large earthquake events and the extreme rainfall events.
Table 6. Comparison of the average annual landslide area decline rate from the large earthquake events and the extreme rainfall events.
EventsOscillating PeriodLADO (km2/Year)LADA (km2/Year)
2005 Kashimir Earthquake2005–20100.320.99
2008 Wenchuan Earthauke2008–20110.410.96
2009 Typhoon Morakot (ARW)2009–20144.246.78
2009 Typhoon Morakot (TRW)2009–20142.376.04
Note: LAD means the average annual landslide area decline rate (km2/year), and LADO and LADA mean the LAD during the oscillating period and after the oscillating period.
Table 7. Statistical data of landslide evolution from 2009 to 2010, 2013, and 2015 in the ARW.
Table 7. Statistical data of landslide evolution from 2009 to 2010, 2013, and 2015 in the ARW.
Landslide TypesR AreaNR AreaNE
AreaGully-RelatedRiver-RelatedLarge-Related
2009–201024.823.54.32.410.021.34
2009–201326.122.25.73.070.051.84
2009–201528.120.27.24.040.082.54
Note: The unit of area in this table is km2. R, NR, and NE mean the recovered, not recovered, and new and enlarged landslide, and the gully-related, river-related, and large-related mean the NE landslide located in the neighborhood of gully, river, and large landslide.
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.; Lin, C. Spatiotemporal Hotspots and Decadal Evolution of Extreme Rainfall-Induced Landslides: Case Studies in Southern Taiwan. Water 2021, 13, 2090. https://doi.org/10.3390/w13152090

AMA Style

Wu C, Lin C. Spatiotemporal Hotspots and Decadal Evolution of Extreme Rainfall-Induced Landslides: Case Studies in Southern Taiwan. Water. 2021; 13(15):2090. https://doi.org/10.3390/w13152090

Chicago/Turabian Style

Wu, Chunhung, and Chengyi Lin. 2021. "Spatiotemporal Hotspots and Decadal Evolution of Extreme Rainfall-Induced Landslides: Case Studies in Southern Taiwan" Water 13, no. 15: 2090. https://doi.org/10.3390/w13152090

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