Consistently over the years, particularly during monsoon seasons, landslides and related geohazards in Bhutan are causing enormous damage to human lives, property, and road networks. The determination of thresholds for rainfall triggered landslides is one of the most effective methods to develop an early warning system. Such thresholds are determined using a variety of rainfall parameters and have been successfully calculated for various regions of the world at different scales. Such thresholds can be used to forecast landslide events which could help in issuing an alert to civic authorities. A comprehensive study on the determination of rainfall thresholds characterizing landslide events for Bhutan is lacking. This paper focuses on defining event rainfall–duration thresholds for Chukha Dzongkhag, situated in south-west Bhutan. The study area is chosen due to the increase in frequency of landslides during monsoon along Phuentsholing-Thimphu highway, which passes through it and this highway is a major trade route of the country with the rest of the world. The present threshold method revolves around the use of a power law equation to determine event rainfall–duration thresholds. The thresholds have been established using available rainfall and landslide data for 2004–2014. The calculated threshold relationship is fitted to the lower boundary of the rainfall conditions leading to landslides and plotted in logarithmic coordinates. The results show that a rainfall event of 24 h with a cumulated rainfall of 53 mm can cause landslides. Later on, the outcome of antecedent rainfall varying from 3–30 days was also analysed to understand its effect on landslide incidences based on cumulative event rainfall. It is also observed that a minimum 10-day antecedent rainfall of 88 mm and a 20-day antecedent rainfall of 142 mm is required for landslide occurrence in the area. The thresholds presented can be improved with the availability of hourly rainfall data and the addition of more landslide data. These can also be used as an early warning system especially along the Phuentsholing–Thimphu Highway to prevent any disruptions of trade.

Landslides are one of the most destructive and most recurring natural calamities in the Himalayan region. Their occurrence leads to immense damage to infrastructure and loss of land, human lives, and livestock. One of the most affected regions is the Bhutan Himalayas, where the majority of the landslides are rainfall-induced. The present study aims to determine the hazard and risk associated with rainfall-induced landslides for the Phuentsholing region located in the southwestern part of the Bhutan Himalayas. The work involves developing a landslide risk map using hazard and vulnerability maps utilizing landslide records from 2004 to 2014. The landslide hazard map was generated by determining spatial and temporal probabilities for the study region. The spatial probability was computed by analyzing the landslide contributing factors like geology, slope, elevation, rainfall, and vegetation based on comprehensive field study and expertise about the area. The contributing factors were divided into various classes and the percentage of landslide occurrence under each class was calculated to understand its contributing significance. Thereafter, a weighted linear combination approach was used in a GIS environment to develop the spatial probability map which was multiplied with temporal probabilities based on regional rainfall thresholds already determined for the region. Consequently, vulnerability assessment was conducted using key elements at risk (population, land use/land cover, proximity to road, proximity to stream) and the weights were provided based on expert judgment and comprehensive field study. Finally, risk was determined and the various regions in the study area were categorized as high, medium, and low risk. Such a study is necessary for low-economic countries like Bhutan which suffers from unavailability of extensive data and research. The study is conducted for a specific region but can be extended to other areas around the investigated area. The tool can serve as an indicator for the civil authorities to analyze the risk posed by landslides due to the rapid infrastructure development in the region.