Topographical features of rainfall-triggered landslides in Mon State, Myanmar, August 2019: spatial distribution heterogeneity and uncommon large relative heights

Continuous 5-day (August 4–9, 2019) torrential rainfall in the monsoon season triggered more than 90 landslides on northwest-southeast extended mountain range of Mon State, Myanmar. In this study, remote sensing images, DEM, and limited fieldworks were used to create the landslide inventory. The topography features of these landslides are analyzed via ArcGIS. The largest one occurred on 9 August 2019 and caused 75 deaths and 27 buildings were damaged. This landslide occurred on gentle topography (slope angle, 23°) with long run-out, in which the angle of reach was relatively low (10°). The volume was 111,878 m3 was mainly composed of weathered granite and red soil and the sliding depth was approximately 7.5 m. Topographic characteristics including the relative slope height, angle of reach, and slope angle of source area of 35 landslides with areas > 4000 m2 were analyzed. The spatial distribution characteristics and topographic features of the 35 landslides below are distinguished: (1) the concentration of most of landslides on southwest-facing slopes showing the heterogeneous spatial distribution of landslide; (2) an uncommon landslide distribution in which more than half of landslide originates from upper slope; (3) the range of the angle of the source area (17°–38°) compatible with the internal friction angle of soils in tropical regions (17°–33°); and (4) the tangent of the angle of reach is generally smaller than 0.5 (angle of reach < 27°) shows a relative high mobility and the relation between landslide mobility and the slope angle of the landslide source area is similar to the one of earthquake-triggered landslides, even though the triggering mechanism, landslide type, and landslide volume are dramatically different.


Introduction
Landslides play a significant role to reshape the earth's surfaces and lead to hazards. Interactions between boundary conditions such as geology, geomorphology, hydrology, and triggering factors like weathering, rainfall, and human activity enhance their complexity. Myanmar is a mountainous country that contains two major ranges of mountains, the Western Ranges and the Eastern Highlands. The topography of mountain regions may influence rainfall patterns because mountains usually act as barriers to wind blows and enhance rainfall on the windward side (Kozo 1980;Parish 1982;O'Connor et al. 1994;Bougeault et al. 2001;Medina and Houze 2003). Infiltration during continuous torrential rainfall elevates pore pressure or destroys the capillary tension forces of soils and rocks and reduces the shear strength of slopes (Iverson 2000;Pailha and Pouliquen 2009;Rahardjo et al. 2010;Marino et al. 2020). Human activities in the mountainous regions of Myanmar related to modern infrastructure such as roads and buildings also have a significant role to increase landslide activities. Yearly, substantial life and socioeconomic losses in Myanmar caused by landslides are a major hazard affecting the country due to the geological, topographical, geomorphological, and meteorological conditions (ADPC 2009;Aung et al. 2017).
This study mainly focuses on the landslides triggered by heavy rainfall during the monsoon season of 2019 in Mon State, Myanmar. In total, 90 landslides were identified and an inventory was produced and documented. The influence of slope aspects that are relevant to metrological conditions on the landslide concentration was explored. Among them, the sloping relief of the 35 landslides having an area larger than 4000 m 2 was analyzed. The geological features (e.g., lithology and weathering conditions) relevant to the largest, catastrophic one, Thae Phyu Kone landslide which occurred on 9 August 2019, were illustrated and the causative factors related to this landslide were discussed in depth.

Site location
The study area lies in Paung County, Mon State, in the eastern region of Myanmar, and 270 km away from the southeast of Yangoon (Fig. 1a). The coastal region of the Gulf of Martaban (Fig. 1b) is approximately 5-10 km west of these ranges, which directly controls the climatic variability of the area. Hilly mountain ranges with flat terrain cover this area. The Mawlamyine Highway connecting Yangoon passes through the area along the foothill side (Fig. 1b, c). The area contains a relatively low thin mountain extending from the northwest to the southeast (Fig. 1c). Landslides are spatially distributed along the southern section of the low relief hill (Fig. 1c, d).

Geological settings
Myanmar lies continuously as Indo Myanmar Ranges in Eastern Himalayan Syntaxis of a continental collision of the Indian and Eurasian plate extended southwards to the north from the Indonesian Arc System to the Eastern Himalayan orogenic belt (Khin et al. 2017;Sloan et al. 2017). Northern Myanmar consists of exposed Indian Plate with various slices of faults also overlain by the Naga hills. Nevertheless, the active subduction of the Indian plate develops transcurrent faults, whereas the pre-collision margin of the Indian Eurasian plate can be evidenced by the most prominent Magok metamorphic belt (Searle et al. 2007;Khin et al. 2017;Sloan et al. 2017), which mainly composed of the oldest rock units in Myanmar, and lies in the major tectonic structural belt. The Magok metamorphic belt (Fig. 2a) mainly composed of gneiss, schist, phyllite, graywacke, dolomite, limestone, and metasedimentary rocks (Mitchell et al. 2012). Intruded igneous rocks are common within this belt. This region frequently experiences landslides and mass wasting due to steep slopes with soils or highly weathered rock masses and torrential monsoons (ADPC 2009).
The study area (red square of Fig. 2b) lies in the Magok metamorphic belt where the carbonate sedimentary rocks and intruded granitic bodies outcrop. The carbonate rocks of phlogopite and diopside marbles are Carboniferous to Devonian time. Between these rocks, intruded granitic/pegmatitic body can be observed. The age of the slightly foliated granite is 70.0 ± 0.9 Ma (Searle et al. 2007). Slightly elongated pegmatitic nodules of quartz-tourmaline-garnet rock with an abundant leucocratic outer zone indicate foliation forming after the granite.

Precipitation
Myanmar is a country with a highly complex spatial hydroclimatic phenomenon with variation in multiple topographical influences of the environment. Mon State, Myanmar, is situated in the tropical monsoon climate region. The temperature is not varied significantly. July and August are the peak of monsoon season and heavy rainfall hit the study area during this period. Since the largest landslide in the study area occurred on 9 August 2019, 8:00 a.m., we illustrate 5-day (4-9 Aug.) daily rainfall records of three rainfall stations (Thaton, Hpa-an, and Mawlayine stations, Fig. 1c) in Mon State (Fig. 3). According to the Department of Metrology and Hydrology (Myanmar), these daily rainfall events were collected during the previous 24 h up to 9:30 a.m. each day. The Mawlamyine station is the closest station to the study area and shows a high rainfall amount (356 mm) from 9:30 a.m., 8 August, to 9:30 a.m., 9 August, which was immediately before the occurrence of the largest landslide event. Thaton and Hpa-an stations in the northern part of the study area recorded relatively low rainfall compared to the rainfall recorded at Mawlamyine station, which shows strong daily rainfall recorded heterogeneity.

Spatial distribution of the identified 90 landslides
A landslide inventory is a detailed description of the spatial distribution, location, and type of landslides (Ayalew et al. 2004;Alessio et al. 2005;Korup 2005;Devoli et al. 2007). We mapped the areal extension of 90 landslides (Fig. 4) that were identified from Google Earth Pro images of 2019/9/26. Landslides were mapped as polygons in kmz format (Online Resource 1). Most of the identified landslides were concentrated on the slopes along the uphill side (northeastern side) of the Mawlamyine Highway. We considered that all of these landslides occurred along with heavy rainfall from the 4th to 9th of August due to lack of landslide activity in past Google Earth images. Most of the landslides in this study area were shallow, debris slides with long transportation paths (Fig. 5). Landslide mapping is not difficult because this disrupted landslide type can be identified easily from the vegetation cover condition. However, human distortions for agriculture and other activities should be carefully excluded. The landslides in this study were mapped at a fixed eye elevation about 500 m. The mapped landslide polygons were further analyzed using ArcGIS software. It should be noted that separating the source, transportation, and deposition areas precisely is difficult because the landslides are small. The size distribution of the 90 landslides is illustrated in Fig. 6. Among them, the topography features of 35 landslides that have areas larger than 4000 m 2 were characterized and will be introduced in Section "Topography features of the 35 landslides with areas > 4000 m2." The largest Thae Phyu Kone landslide was studied via field investigation and will be introduced in more detail in Section "Thae Phye Kone landslide".

Topography features of the 35 landslides with areas > 4000 m 2
Topographical factors account for the instability of slopes (Montgomery and Dietrich 1994;Talebi et al. 2008;Hess et al. 2017;Cogan and Gratchev 2019). Understanding the topographical factors of the studied shallow landslides is significantly important. The 35 landslides, whose areas were larger than 4000 m 2 , are illustrated in Fig. 7. The topography features including slope aspect, relative slope height, angle of reach, and slope angle of the landslide source area were analyzed. The 30*30-m pixel SRTM (Shuttle Radar Topography Mission) was used for the topographical analysis. The topographical parameters used in this study are presented in Fig. 8. The relative slope height (h/H) is the slope height (h) divided by the elevation difference between the ridge, top, and the valley floor (H) deals with the topographical position (Lee 2013). Angle represents the slope angle of landslide source area. The angle connecting from the top of the source area to the bottom of the deposition area is termed the angle of reach ( ) (Heim 1932), which evaluates the mobility of landslides (Hsu 1975;Iverson 1997;Okura et al. 2003).

Relative slope height (h/H)
The main common feature of rainfall-induced landslides is a favorable location on the lower slope, which can be characterized by a small value of relative slope height (Lee 2013). Herein, the landslides are analyzed via relative slope height (Fig. 9). Less than half of the 35 landslides were concentrated at the lower slope (h/H < 0.5). There were still many landslides located at the upslope even near the ridges of slopes, which was not expected.
Distribution of slope angle of landslide source area ( ) and angle of reach ( ) The most important topographic factor that dominates the initiation of shallow landslides is the slope gradient (Mulder and Alexander 2001). The distribution of the slope angle of the landslide source area ( ) of the 35 landslides is illustrated in Fig. 10a. The determined ranged from 17° to 38° with a mean value of 27° (Online Resource 2).
The angle of reach of landslides is influenced by the debris properties, geometry of the flowing path, and presence of obstacles (Corominas 1996). Figure 10b shows the distribution of the angle of reach β of the 35 landslides. The smallest angle of reach β of the studied landslide (the largest the Thae Phye Kone Landslide) is 10°, which will be introduced later. The angles of reaches are 31°, 31°, 34°, and 34° for landslides IDs 22, 25, 28, and 33, respectively (Fig. 7). The transportation paths of these landslides were roughly perpendicular to the main creek and the runouts were largely constrained, accounting for the high angle of reach. The angle of the source area was relatively higher than the angle of reach (Fig. 9a, b).

Fig. 3
Five-day (a) accumulated rainfall; (b) daily rainfall amounts in Mon State close to the deadly Thae Phye Kone landslide, which occurred on 8:00 a.m. on 9 August 2019, rainfall data of each day recorded by 9:30 a.m. of last 24 h. The Mawlamyine station recorded the highest rainfall from 9:30 a.m., 8 August, to 9:30 a.m., 9 August, which is significantly close to the study area (DMH 2019) Thae Phye Kone landslide On 9 August 2019, at approximately 8:00 a.m., after long-lasting heavy rainfall (Fig. 3), a massive landslide hit Thae Phye Kone village (16°33′40.95″N latitude and 97°34′46.22″E longitude), north of Mottama, Mon State. Thirty buildings were damaged, and a total of 75 people were killed due to this slope failure, together with debris slides of red soil and rock fragments (Fig. 11a, b). The huge mass of debris on the highway buried several vehicles and passengers. Twenty-seven buildings were damaged by sliding debris (BBC 2019).
Based on the satellite images, pre-event SRTM, and site investigations, the features of the Thae Phyue Kone landslide were characterized (Table 1). The Thae Phyue Kone landslide is a weathered  This landslide was investigated via field investigation which will be described in Section "Thae Phye Kone landslide" rock/soil slide. The weathered rock fragments and residual soils slid down from the ridge of the hill to the toe of the slope. Soils and weathered rock debris spread widely on flat areas where the vehicles on the road and buildings were buried suddenly. A pagoda located on the lower part of the slope was badly damaged by the big boulders (Fig. 12). Based on the pre-and post-satellite images, the landslide covered 48,558 m 2 . The area was calculated by pixel counting via ArcGIS of the polygon-based projected area. The maximum projected length of this landslide was approximately 400 m, and the maximum projected width is approximately 170 m.
The large boulders in the deposition area are mainly composed of granite with different weathering degrees. The soils are red in color and clay-rich, which came from the thick residual soils on the top of a ridge (Fig. 12).
The whole landslide area is subdivided into three different parts based on slope geometry and material properties. The crown (point c), source area (points c to b), and deposition area (points b to a) of the Thae Phyue Kone landslide were identified and marked in Fig. 13a. The deposition area with large rock fragments was also marked. The profile of the Thae Phyue Kone landslide is illustrated in Fig. 13b. The exposed source area is close to the slope top, and the relative slope height is 0.7. The exposed source area was 5922 m 2 , and the slope angle of the landslide source area was 23°. Since part of the source could still be buried by the rock fragments, the true source area could be larger than the exposed source area. We define a transitional/depositional area, together with the rock fragment area (Fig. 13a), which means that the area includes erosion but is covered by moving debris. The transitional/depositional area (5912 m 2 ) and rock fragment area (3132 m 2 ) cover 9044 m 2 , and the total source area could be approximately 14,966 m 2 .
A weathered granitic rock mass with several vertical and inclined joints can be observed on the landslide crown. Spheroidal weathering could dominate the mass wasting process. The large boulders of rock fragments near the source area are mainly the rock cores of the spheroidal weathering process. The angle of repose of the rock fragments is approximately 20°, and the rock to debris slide is presumed. The deposition area lies on the toe part of the landslide and the angle of reach is only 10°. The deposition area is approximately 33,619 m 2 , and the average thickness is estimated approximately 3 m from the height of buried buildings. The details of the depth estimation and volume calculation are presented in Table 2. If the deposition depth of the transitional/deposition area (5912 m 2 ) is 3 m too and the deposition depth of the rock fragment area (3132 m 2 ) is 5 m, the deposition volume (after expansion) of this landslide is approximately 134,253 m 3 . Since the estimated source area is 14,966 m 2 , the average sliding depth is approximately 7.5 m if an expansion ratio of 1.2 (Guzzetti et al. 2009) is assumed (i.e., the landslide volume before expansion is 111,878 m 3 ).

Heterogeneity of rainfall and landslide spatial distribution
In the study area, a thin mountain extending from northwest to southeast ( Fig. 14; green line is a ridge separating the southwestand northeast-facing slopes) could induce rainfall heterogeneity and dominate the spatial distribution of landslides. In total, 77 out of 90 landslides were located on the southwest side of the ridge (Fig. 14). We suggest that this spatial distribution of rainfallinduced landslides can be correlated with rainfall pattern heterogeneity in the study area. High temporal and spatial heterogeneity in rainfall can be observed in Fig. 3, wherein stations Thaton and Hpa-an recorded lower daily rainfall (61 and 69 mm, respectively) than the Mawlamyine station (356 mm) during the last 24 h before 9:30 a.m. on August 9, 2019. This temporal and spatial rainfall heterogeneity is observed clearly from radar data (Fig. 15a-d). Figure 15a shows that rainfall only accumulated on the northern and southern sides of the study area at 18:48 p.m. on August 7, 2019. At 15:51 p.m. on August 8, 2019, the rainfall area concentrated more to the study area (Fig. 15b). The rainfall started to increase near the study area at 18:47 p.m. on August 8, 2019 (Fig. 15c), and the rainfall intensity reached approximately 100 mm/h on the southwestern side of the study area at 7:03 a.m. on August 9, 2019, which was approximately 2 h before the deadly Thae Phye Kone landslide events. It can be observed that the rainfall intensity is lower on the northeastern side.
Rainfall patterns globally and locally are spatially heterogeneous due to topography effects (Smith 1979). The topographical conditions in the study area include sequences of mountain ranges, and valleys lie very close to the Gulf of Martaban and Bay of Bengal. In southern Myanmar, the direction of monsoon wind in this region is from southwest to northeast (Aung et al. 2017). On the windward side, higher rainfall due to the lifting of air masses releases an increase in rainfall but decreases the precipitation on the other side of mountain ranges (Smith 1979;Hill 1983). Mountains act as a barrier to wind flow and automatically facilitate the rainfall process. This type of barrier mountain range is found in many other regions of the entire world and is termed the orographic effect (Pierrehumbert and Wyman 1985;Smith 1989). The enhanced rainfall on the windward side can percolate more into the subsurface of slopes and reduce the stability of slopes (Parish 1982;McCauley and Sturman 1999;Houze 2012). However, it is difficult to conclude that the spatial distribution features of the studied landslides are due to orographic effects  due to the lack of dense precipitation station arrays near the study area.

High relative slope height, gentle slope angle of the landslide source area, and high mobility
As shown in Section "Topography features of the 35 landslides with areas > 4000 m 2 ," the distinct topographic features (Figs. 9 and 10) of the 35 landslides included (1) high relative slope height (h/H), (2) gentle slope angle of the landslide source area ( ), and (3) high mobility (small angle of reach; ): (1) High relative slope height (h/H): Rainfall-induced landslides are favorably located on the lower slope (Lee 2013) because the groundwater table upslope is generally deeper than that of the downslope. It is interesting to find that almost half of the studied landslides out of 35 (area having > 4000 m 2 ) are located near the top of the slope (h/H > 0.5). Among them, 4 started from the top of the slope (Fig. 9), which is not common for rainfall-triggered landslides.
(2) Gentle slope angle of the landslide source area ( ): The slope angle of the landslide source area δ ranged from 17 to 38°, with a mean value of 26° (Online Resource 2). This relatively gentle slope angle of the landslide source area ( ) is correlated to the high relative slope height (h/H). Generally, the slope angle near the ridge is relatively gentle. Therefore, soil erosion was minimized near the ridge because the slope-steepness factor and slope-length factor (factors in USELE) are low (Wischmeier and Smith 1978;Renard et al. 1997). Thick weathered regolith can therefore accumulate easily at gentle topography near the ridge, although the soil production rate will decrease when the soil depth increases (Heimsath et al. 1999).
The slope angle of the landslide source area positively correlates to the friction angle of the sliding surface of landslides. Rigo et al. (2011) collected the residual shear strength of tropical soil throughout the world, including Southeast Asia, Africa, and South America, from different studies. The internal friction angle of the migmatitic, micaceous saprolitic soil of South America ranges from 17° to 33°. The main lithology are migmatitic granites of the study area, and the documented internal friction angles are compatible with the slope angle of the landslide source area in the study area (ranging from 17° to 38°).
(3) High mobility (small angle of reach; ): A tan roughly equals 0.5 for saturated, poorly sorted debris flows along unconfined paths (Iverson 1997). Among the 35 studied landslides, 28 landslides had tan (angle of reach) less than 0.5 ( < 27 • ) (Online Resource 2). The angle of reach of the studied landslides was generally smaller than the angle of the landslide source area ( ) (Fig. 9a,  Based on the three features observed, possible explanations are discussed. First, the angle of the source area is relatively higher than the angle of reach (Fig. 9a, b), implying that the role of fluid (runoff mixed with debris) is significant . The low angle of reach (10°) of the Thae Phye Kone landslide could be related to the low strength of weathered clay-rich soils mixed with runoff or groundwater. This promotes the rock (debris) to slide into debris (soil) flow with long run-out. The mechanism of undrained loading lowering the angle of reach due to an increase in pore water pressure (Hutchinson and Bandhari 1971;Sassa et al. 1991) can be another explanation of the observed long run-out.
Second, for dry slopes without cohesion, the critical slope angle equals the drain friction angle of the sliding debris for infinite slope instability. When the groundwater table is identical to the slope surface and the groundwater flow is parallel to the slope surface, the critical slope angle will be half of the drain friction angle of the sliding debris (Lambe and Whitman 1969). Fig. 13 The scarp, source, and deposition area of Thae Phyue Kone landslide. (a) Plan view; (b) profile. The areas are divided using Google Earth images (Fig. 11b) and a digital elevation model. Two more areas were identified: transitional/deposition area and rock fragment area, where the source area could be covered by debris or rock fragments The slope angle of the landslide source area compatible with the documented friction angle of saprolitic soil suggests that the failure of the studied landslides cannot be caused by fully saturated slopes with strong seepage force. Finally, rainfall-triggered landslides usually originate from an increase in the groundwater table (pore pressure elevation). The groundwater table near the top of the slope is generally low, and the pore pressure rise due to the elevated groundwater table is not a good explanation for the occurrence of landslides upslope, as observed in this study. Alternatively, during infiltration, wetting front migration on unsaturated upslopes reduces matrix suction, and upslopes that become unstable are frequently observed in tropical areas (Fredlund et al. 1978;Fredlund and Rahardjo 1993;Rahardjo et al. 1995), even though the sliding surface is well above the groundwater table. The thick weathered regolith accumulated on gentle topography near the ridge can enhance infiltration during rainfall. We speculate that the saturation of thick weathered regolith from the slope surface due to rainfall infiltration in the study area could episodically promote shallow landslides on the top of slopes. However, this argument has not been well constrained, and convincing evidence could be obtained via further detailed site investigation and monitoring projects, which is out of the scope of this study.
Relation between the slope angle of the landslide source area and angle of reach (landslide mobility) The mobility of landslides is critical and complex phenomena from the viewpoint of landslide hazard mitigation. Landslide mobility is mainly controlled by topographical and morphological factors such as volume of the sliding mass, properties of debris, motion mechanism, pore water pressures, geometry or regularity of the traveling path, and presence of obstacles such as vegetation and boulders (Heim 1932;Hsu 1975;Coraminas 1996;Hunter and Fell 2003). Zhan et al. (2017) presented a positive correlation between the tangent of the angle of reach and the slope angle of the landslide source area ( tan ) of 38 channelized rock avalanches induced by the 2008 Wenchuan earthquake. The data points and regression line are shown in Fig. 16 (red circles and straight line). Similarly, rainfall induced 35 Myanmar landslides in Fig. 16 (blue circles) which suggests that the tangent of the angle of reach ( tan ) correlated positively with the tangent of the slope angle of the landslide source area. Landslide mobility is decreased (the angle of reach increases) with an increasing slope angle in the landslide source area reflecting the fact that the slope angle of the landslide source area correlates the friction angle in sliding debris positively (Okura et al. 2003;Zhan et al. 2017). This illustrated that the kinetic energy of the sliding mass Fig. 16 Relationship between the tangent of the angle of reach (tan β or H/B) and tangent of the slope angle of landslide source area ( tan ) where red circles are data of earthquake-induced landslides (Zhan et al. 2017), and blue circles are data on rainfall-induced landslides (this study) consumed by internal friction is positively correlated to the slope angle of the landslide source area (Iverson 1997). The data scattering in Fig. 16 implies that the angle of reach could be controlled by other factors, such as strength heterogeneity and local topography aside from the valley rather than slope gradient, which could be an obstacle for the mobility of landslides and increase the angle of reach. Another possible reason are obstructions by dense vegetation and boulders for shallow landslides (Hunter and Fell 2003). It is interesting to find that the data points of earthquake-induced and rainfall-induced landslides are quite similar in Fig. 16. Notably, the landslide volumes of these two datasets are dramatically different. This issue deserved more study in the future.

Conclusions
We created an inventory of a total of 90 landslides triggered by heavy rainfall during the monsoon season of 2019 in Mon State, Myanmar. Among them, 35 landslides with areas > 4000 m 2 were carefully analyzed using Google Earth imagery and SRTM with ArcGIS software. The topographic characteristics of these 35 landslides, including the relative slope height, angle of reach, and angle of the source area, were documented. The biggest and most catastrophic Thae Phye Kone landslide was studied via field investigations, and the characteristics were documented in more detail. The critical findings can be summarized in the following headings: 1. On 9 August 2019, the biggest Thae Phye Kone landslide, triggered by long-lasting heavy rainfall (more than 3 days with approximately 100 mm rainfall per day) killed 75 people and 27 buildings were damaged. The landslide was mainly composed of weathered granite and residual soils in red, located on a gentle slope (slope angle 23°) near the slope top (relative height = 0.7). The volume was approximately 111,878 m 3 , and the average sliding depth is approximately 7.5 m. The angle of reach was quite low (~ 10°), indicating the high mobility of this catastrophic landslide. 2. A highly heterogeneous rainfall pattern from August 4 to 9 was observed. The mountain range extending from northwest to southeast could contribute to rainfall heterogeneity. In total, 77 out of 90 landslides located on the southwest side of the ridge showed strong spatial distribution heterogeneity of landslide occurrence, which could be related to the heterogeneous rainfall pattern.
3. More than half of 35 landslides with area > 4000 m 2 were located on the upper slopes (relative height > 0.5) including 4 landslides on the slope top, which is uncommon for rainfall-triggered landslides. The slope angle of the source area ranges from 17° to 38°, which is compatible with the internal friction angles of saprotlitic soils in tropical regions (17°-33°). The groundwater table rise or seepage force caused by groundwater flow in saturated soils fails to account for the observations. 4. The tangent of the angle of reach is generally smaller than 0.5 (angle of reach < 27°), which indicates a relatively high mobility. The relation between tan (angle of reach) and tan (slope angle of source area) of the studied rainfall-triggered landslides surprisingly resembles the relation of earthquake-induced landslides, although the triggering mechanism, landslide type, and landslide volume are dramatically different.