1 Introduction

Landslides are a significant natural hazard that poses a threat to human lives, river health, and infrastructure in many regions globally [1]. The rise in human population and urbanization in hilly regions has led to a rise in the number of landslides, causing significant economic losses and loss of biodiversity globally [2].

Landslides are a significant natural hazard in India, causing extensive damage to property and biodiversity [3]. India’s unique geology and climatic conditions make it particularly vulnerable to landslides, with over 12.6% of the country’s landmass prone to landslides [4]. In recent years, the landslides' frequency and intensity in India have increased due to a lot of factors such as improper land use practices, deforestation, induced rainfall patterns, and climate change [5,6,7]. Landslide events in India have claimed 3971 lives over the period of 300 years in 248 fatal events. In addition to the loss of human lives, landslides in India have also caused socio-economic damage in 371 significant events [8]. Given the scale and severity of landslides in India, effective and comprehensive landslide risk management is essential. One of the essential steps in landslide risk management is the landslide hazard zonation mapping, which identifies areas susceptible to landslide occurrence. This can be accomplished by using advanced remote sensing and GIS techniques, which can provide valuable insights into the various factors contributing to landslides and help in identifying high-risk areas [9].

The application of geospatial technologies is a crucial strategy for determining landslide hazard zones. The incorporation of spatial data on geology, slope, elevation, soil, LULC, and other factors with remote sensing imagery provides an in-depth knowledge of landslide-prone areas [10,11,12]. This approach allows planners and decision-makers to mitigate the effects of landslides by identifying areas where preventive measures and early warning systems should be implemented.

Numerous research was done to map the landslide hazard zones utilizing remote sensing and GIS techniques [13,14,15,16,17,18,19,20,21,22]. For instance, Wang et al. in 2018 used remote sensing and GIS data to create a landslide susceptibility map of Yihuang County in China [23]. Similarly, Huang et al. in 2020 utilized remote sensing and GIS-based approaches to map landslide hazards in the Guangxi Province of China [24]. Another research by Taloor et al. [25] examined the landslide susceptibility of Doda Kishtwar ramban (DKR) region of Jammu and Kashmir based on the weighted overlay method where they found 46% of the total area over moderate to very high risk zone. Shit et al. in 2016 studied landslide risk in Maldevta, Dehradun, finding that around 33% of the area is categorized as having high to very high landslide risk using a weighted overlay model (WOM) [26].

In order to comprehend and manage landslide hazards, the current study focuses on the landslide hazard zonation mapping using remote sensing and GIS techniques. Planners and decision-makers can use a comprehensive understanding of landslide-prone areas provided by the integration of geospatial data and remote sensing imagery to mitigate the effects of landslides.

2 Study area

The Doon Valley is located in the middle of the Lesser Himalayas, and the Shiwalik mountain ranges in Uttarakhand [27]. The area is coordinated between 77° 38ʹ E and 78° 20ʹ E in longitude and 30° 01ʹ N to 30° 28ʹ N in latitude. The valley is comprised of two watersheds, Asan and Song [28]. The Asan watershed has an area of 701.15 km2, and the Song watershed has an area of 1040.49 km2 [29]. The outlet of the Asan watershed is in the Yamuna River, and the outlet of the Song watershed is in the Ganga River. The total area of the Doon Valley is 1741.64 km2, with elevations that range from 303 to 2764 m. The major rivers in the Doon Valley are Bindal, Bandal, Song, Asan, Tons Rispana, and Suswa. The region receives ~ 2000 mm of average annual rainfall [30]. In general, the region experiences hot summers and chilly winters. The yearly range of temperature is between − 3 °C in the winter and 45 °C in the summer. Dehradun, Mussoorie, Vikasnagar, Doiwala, and Rishikesh are the important cities in the study area. The Study Area map is indicated below in (Fig. 1).

Fig. 1
figure 1

The map of the study area

2.1 Methodology

3 Data and methods

In this present study, the integrated use of satellite and DEM data was used for the generation of spatial layers for the landslide hazard zonation mapping. The weighted overlay method was used to calculate a multi-criteria analysis between several raster layers. This method allows for the integration of multiple raster layers by reclassifying each layer to a common scale, assigning them relative weights based on their importance, and combining them to produce a single output raster layer. The process begins with defining the problem and identifying the relevant criteria. Next, the necessary raster layers are collected and preprocessed to ensure they have the same extent. Each raster layer is then reclassified to a standardized scale, typically from 1 to 9, where higher values represent greater risk or susceptibility. Weights are assigned to each criterion to reflect its relative importance, ensuring the total weights sum to 100%. Using the Weighted Overlay tool in ArcGIS, these reclassified layers are combined according to their assigned weights (1). The resulting output raster highlights areas of interest based on the combined influence of all criteria, facilitating informed decision-making and planning.

$$ {\text{Output}}\;{\text{Raster}} = \sum \left( {{\text{Reclassified}}\;{\text{Raster}} \times {\text{Weight}}} \right) $$
(1)

For the LULC classification, the supervised classification method was applied using the ERDAS Imagine software. The maximum likelihood classifier was used to identify each training samples. Based on the training samples, the LULC classes was defined. The ArcGIS software was used to create maps of the Doon Valley’s aspect, slope, elevation, geology, drainage density, vegetation, LULC, soil, and landslide hazard zone. The source of data and flow chart of the methodology explained in Table 1 and Fig. 2.

Table 1 Data source and data type
Fig. 2
figure 2

Flowchart

3.1 Flowchart of the methodology

See Fig. 2.

4 Results and discussion

4.1 Slope

The area’s steepness is indicated by the slope [31]. According to Fig. 3, the Doon Valley slope is categorized into five classes: (0° to ≤ 5°) is very gentle, having 50.01% of the total area, (> 5° to ≤ 15°) is gentle, having 21.21% of the total area, (> 15° to ≤ 30°) is Moderate having 16.77% of the total area, (> 30° to ≤ 45°) is steep having 10.84% of total area and (> 45° to 75°) is very steep having 1.17% of total area. The table and map make it clear that half of the area is having very gentle slope (Table 2). With little runoff and a gradual slope, the valley is favorable for groundwater penetration. Based on how frequently certain slope angles occur, slope categories are defined by slope maps. The geomorphological history of the landscape determines how the slope categories are distributed [32].

Fig. 3
figure 3

Slope map

Table 2 Slope

4.2 Elevation

According to Fig. 4 & Table 3, there are five different elevation classifications in the Doon Valley: (≤ 600 m) is very low, having 38.96% of the total area, (> 600 m to ≤ 900 m) is low, having 35.47% of the total area, (> 900 m to ≤ 1200 m) is moderate having 7.55% of the total area, > 1200 m to ≤ 1600 m) is high having 8.4% of the total area and > 1600 m to 2764 m) is very high having 9.62% of the total area. Most of the area has an elevation value that is low to very low, indicating that elevation and drainage density are related indirectly, whereas elevation and slope are related directly [33].

Fig. 4
figure 4

Elevation map

Table 3 Elevation

4.3 Drainage density

The total length of streams per unit area is what is referred to as the drainage density. The subdivision of the landscape and the catchment area’s potential for runoff are measured by drainage density. A greater degree of land subdivision indicates a greater probability of slope failure [34]. As indicated in Fig. 5, the drainage density is classified into five classes: (≤ 3 km/km2) is very low, having 11.94% of the total area, (> 3 km/km2 to ≤ 5 km/km2) is low having 35.68% of the total area, (> 5 km/km2 to ≤ 7 km/km2) is moderate having 37.79% of the total area, (> 7 km/km2 to ≤ 9 km/km2) is high having 11.35% of the total area and (> 9 km/km2 to 10.60 km/km2) is very high having 3.24% of the total area. The majority portion of the total area has a medium to low drainage density (Table 4). The red color on the map denotes the area with the highest drainage density. The subdivision of the landscape and the catchment area’s potential for runoff are measured by drainage density. A greater degree of land subdivision indicates a greater probability of slope failure.

Fig. 5
figure 5

Drainage density map

Table 4 Drainage density

4.4 Geology

Geology formation becomes important in the study of landslides. The 11 classes of geology formation were shown in Fig. 6: Alluvium, Blaini/Infrakrol, Chandpur, Damta, Krol Sandstone + Krol A, Low Grade, Lower Tal, Middle Shiwalik, Nagthat, Upper Shiwalik, and Upper Tal. The Alluvium Formation has 62.62% of the total area, followed by the Upper Shiwalik Formation (Table 5). The Krol, Nagthat, Blaini, and Chandpur formations are more prone to landslides as compared with the Shiwalik formations [35, 36].

Fig. 6
figure 6

Geology map

Table 5 Geology formation

4.5 Soil type

Soil Type is an important parameter for LHZ Mapping. The Soil Type was classified into five Classes (Fig. 7): Sandy-Skeletal, Loamy-Skeletal, Coarse-Loamy, Fine-Loamy, Clayey, and Waterbodies [37]. The Majority of the area has Loamy-Skeletal soil types followed by Coarse-Loamy (Table 6). Loamy-Skeletal and Clayey soil types are more susceptible to landslides due to the presence of colluvium as a parent material.

Fig. 7
figure 7

Soil map

Table 6 Soil types

4.6 Vegetation

Vegetation Types play an important role in landslide susceptibility mapping. The vegetation type map was created using NDVI, or Normalized Differential Vegetation Index. The formula for calculating the NDVI is (NIR-Red)/(NIR + Red) where (NIR) is (Band 8) and (Red) is (Band 4) of Sentinel-2B. The three classes of vegetation types are shown in Fig. 8: Dense Vegetation, Sparse Vegetation, and Barren Land. According to the vegetation types, Dense vegetation trees/plants with deep roots have contributed to an improvement in slope stability [38]. The Barren Land class, 23.1% of the total area, is more susceptible to landslide (Table 7).

Fig. 8
figure 8

Vegetation types map

Table 7 Vegetation types

4.7 Land use land cover LULC

Another essential factor for determining the Landslide Hazard Zone of any region is the LULC classification. It illustrates how anthropogenic activities, particularly urbanization, and agriculture, utilized land resources [39]. Due to climate change and land use practices, water resources are constantly under intense pressure. The LULC can be used to extrapolate the complicated physical processes occurring on the earth's surface [40, 41]. Sentinel-2B imagery was used to create the LULC map. For LULC classification, the ERDAS Imagine software uses the supervised classification method. There are 6 LULC classes, as illustrated in Fig. 9: built-up, Agriculture Land, Forest Area, Wasteland, Dry Riverbed, and Waterbodies. The map’s conclusion shows that Forest Area covered more than half of the total area (69.21%), followed by Agricultural Land (11.49%) (Table 8). Another dominant class is the Built-up, which takes up 9.09% of the space. The center of the valley has an anthropogenic load that directly affects the biodiversity and health of the rivers like Rispana, Bindal, Suswa, and Asan. The Wasteland covers the area 7.27% area of the Doon Valley. The two categories taking up the least amount of area are Waterbodies and Dry Riverbeds. Increasing tourist pressure in hilly region, changing climate are potentials risks due to landslide [42]. Tourism activities can lead to deforestation, soil erosion, and infrastructure stress, while climate change contributes to more intense and frequent rainfall. Together, these factors destabilize slopes, increasing the likelihood of landslides. Runoff from agricultural activities can significantly increase the risk of erosion and landslides in riverine areas. Excess water from irrigation and rainfall reduce soil stability. This, combined with the removal of vegetation for agriculture, can accelerate soil erosion and trigger landslides along riverbanks [43].

Fig. 9
figure 9

LULC map

Table 8 LULC classes

4.8 Aspect

The Aspect displays the slope’s direction. The Aspect at (0°–22.5°) is north, at (22.5°–67.5°) Northeast, and so on. For this study, the direction of the slope is on two sides, one west-facing and another southwest- south facing as shown in Fig. 10 & Table 9. In comparison to the east-facing slope, the west and south-facing slope has more vegetation cover and more moisture.

Fig. 10
figure 10

Aspect map

Table 9 Aspect classes

4.9 Landslide hazard zonation (LHZ) mapping

The LHZ mapping is a crucial tool for mitigating the risk associated with landslides. Different spatial data sources, including topography, geology, LULC, vegetation, and soil characteristics, must be integrated during the process [44]. The LHZ map gives a thorough picture of how the landslide risk is distributed spatially within an area and can be used to guide decision-making for disaster management, land use planning, and infrastructure development. The LHZ map was created using the Weighted Overlay method in ArcGIS Software. Weighted Overlay is a multi-criteria analysis between several raster layers. For this study, the eight raster layers were used and given the rank (Ri) from 8 to 1, where 8 means highest and 1 means lowest. Each raster layer is further divided into different classes and given the weightage (Wi) from 1 (lowest) to 8 (highest). Also, the percent influence (%) of each raster layer from the overall 100% is mentioned in Table 10.

Table 10 Rank, influence and weightage of factors for landslide

As depicted in Fig. 11, there are five classes that make up the Landslide Hazard Zones: very high hazard zone, high hazard zone, medium hazard zone, low hazard zone, and very low hazard zone. According to Table 11, 16.64% of the total area is located near and above the Main Boundary Thrust/Fault (MBT), which is prone to landslides, and is classified as being in the high to very high hard zone. The upper stretches of the Doon Valley rivers are more significant to landslides than the lower stretches because steep slope, high elevation and presence of Krol, Nagthat, Blaini, and Chandpur geological formations. River ecosystems can be significantly impacted by landslides, which can change the water's physical and chemical properties as well as the river's biodiversity. Landslides can result in significant amounts of sediment and debris entering rivers, which lowers water quality and raises turbidity [45]. The loss of crucial fish breeding habitat and the deposition of fine sediment in river beds are both potential effects of the increased sediment load [46]. Recent event of cloud burst on 19th and 20th August, 2022 have triggered the landslide in many regions of Doon Valley and washes away the banks of river Bandal and Song [47]. Increasing population and uncontrolled tourism in the valley is contributing to the occurrence of landslides by placing additional pressure on vulnerable landscapes and altering natural terrain [48]. Not limited to this, landslides also cause changing in land uses of a particular area. The removal of vegetation cover due to landslides may cause more soil erosion, which in turn degrades the soil’s quality and limits its capacity to support plant development [49]. In general, landslides can have a major impact on an area’s LULC pattern, leading to a loss of vegetation cover, an increase in rocky areas and barren land, as well as a decrease in agricultural output. These effects may have long-term consequences for an area's environment and biodiversity, emphasizing the necessity for efficient landslide risk management and mitigation measures.

Fig. 11
figure 11

Landslide hazard zonation map with photographs of landslides in Doon Valley

Table 11 LHZ table

5 Conclusion

This study demonstrates a clear link between the factors: LULC, Slope, Geology, Soil Type, Elevation, Drainage Density, Vegetation, and Aspect, all are significant for the LHZ mapping. In-depth knowledge of the numerous landslide-causing factors and identification of the high-risk areas can be achieved through the use of advanced remote sensing and GIS techniques. The approach for the LHZ mapping includes the creation of thematic data layers, the incorporation of spatial data, and the validation of outcomes. The findings of this study clearly demonstrate that the area near and above the Main Boundary Thrust/Fault (MBT) is more susceptible to landslide. Also, the upper stretches of the Doon Valley rivers are much more vulnerable to landslides for habitat and biodiversity loss than the lower stretches. Anthropogenic activities such as construction of new houses/buildings and development of new roads/highways near the steep slope and very steep slope regions of the Doon Valley should be managed or planned sustainably. Additionally, understanding geological formations with changing land uses is vital for determining the region's susceptibility to landslides.

Overall, the LHZ mapping is a crucial part of risk assessment and management plans. It offers a reliable tool for decision-makers to comprehend possible dangerous locations and put in place the necessary safeguards to reduce the effects of landslides. Further research in this area is necessary to improve the accuracy of landslide hazard zonation maps and develop better mitigation strategies. For effective future use of the LHZ mapping, consider the following key recommendations: integrate multi-temporal satellite imagery to analyze land use and vegetation changes over time, enhancing trend identification and prediction. Apply machine learning algorithms for advanced pattern recognition and predictive modeling. Combine geological, hydrological, meteorological, and land use data for a comprehensive landslide risk model. Engage local communities in data collection and validation, and conduct regular field surveys to ensure map accuracy.