Next Article in Journal
A Neural Network Model for Estimation of Failure Stresses and Strains in Cohesive Soils
Previous Article in Journal
In Situ Skin Friction Capacity Modeling with Advanced Neuro-Fuzzy Optimized by Metaheuristic Algorithms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geological Exploration, Landslide Characterization and Susceptibility Mapping at the Boundary between Two Crystalline Bodies in Jajarkot, Nepal

1
Central Department of Geology, Tribhuvan University, Kirtipur, Kathmandu 44618, Nepal
2
NEA Engineering Company Limited, Trade Tower, Thapathali, Kathmandu 44600, Nepal
3
Department of Geology, Tri-Chandra Multiple Campus, Tribhuvan University, Kathmandu 44605, Nepal
*
Author to whom correspondence should be addressed.
Present address: Department of Civil and Environmental Engineering, University of Toledo, Toledo, OH 43606, USA.
Geotechnics 2022, 2(4), 1059-1083; https://doi.org/10.3390/geotechnics2040050
Submission received: 18 October 2022 / Revised: 16 November 2022 / Accepted: 25 November 2022 / Published: 1 December 2022
(This article belongs to the Topic Support Theory and Technology of Geotechnical Engineering)

Abstract

:
The geology of the Himalayas is intricated and intriguing. It features numerous tectonic bodies and structures too complex to interpret. Along with such mysteries it has too many common geohazards, such as landslides. In this study, a detailed geological investigation is carried out to overcome the discrepancies in structural interpretation, the nature of two crystalline bodies, and non-uniformity in geological mapping in the central Himalayan arc, in the Jajarkot district of Nepal. Along with the geological exploration and landslide characterization of the area, consequent landslide susceptibility mapping is performed considering 13 causative factors related to geology, topography, land use, hydrology, and the anthropogenic factor, using two bivariate statistical models. This study concludes that the two metamorphic crystalline bodies in the study area are most probably the klippen, due to the absence or erosion of the root zone. The field study revealed that haphazard road excavation without the consideration of geological and geotechnical features has caused shallow landslides. The results obtained from the susceptibility maps, with a varying range of susceptibility zones, are in good agreement with the spatial distribution of pre-historic landslides. The results of the susceptibility modeling are validated by calculating landslide density and plotting area under curves (AUC). The AUC value for the WOE, and the FR method, revealed an overall success rate of 79.42% and 77.62%, respectively.

1. Introduction

The Himalayan Orogeny is the perfect example of continent-continent collision that occurred obliquely at around 65 Ma [1], when the Indian plate moved towards the Eurasian plate. This collision was later followed by anticlockwise rotation of the Indian plate and collision around 50 Ma [2]. The conversion velocity of the Indian plate and Eurasian plate decreased significantly from 15 to 4 cm/year at about 50–35 Ma [3]. Still, the Indian plate is moving towards the north and is pushing the Eurasian plate with a movement speed of 2 cm/year, making the youngest mountain range tectonically still active [4,5], and making the area vulnerable to geohazards.
Geo-hazards such as landslides, debris flows, glacier lake outburst flood (GLOF), mudflows, earthquake and earthquake-triggered landslides, and debris flows are common in the fragile Himalayas [6]. The growing population and development of road networks, infrastructures, settlements, and lifelines have increased the frequency and impact of natural hazards, especially landslides and debris flows [7]. Worldwide, between the years 1994 and 2013, 218 million people were affected by natural disasters annually [8]. Among all the natural hazards, landslides and debris flows are the most frequent and are responsible for the destruction of lives, property, and natural settings. Landslides and debris flow account for 17% of all deaths caused by natural hazards [9]. Landslides are caused by factors such as incompetent lithological and geological conditions, rugged topography and relief, and are triggered by natural agents such as rainfall [10,11] and anthropogenic actions such as land use pattern [12]. Landslides are sometimes triggered by the movement of active faults and earthquakes [13]. Debris flows are often mobilized from shallow landslides [14,15] and later converted into extreme flows with the addition of water and the entrainment of solid mass along the downslope.
Landslides and debris flows are most frequent within the region of active tectonism with distinctive geological structures and differential lithostratigraphy. The lesser Himalayan region in the Himalayan arc has suffered from multiple events of mass movements because of fragile geology and high-intensity rainfall [16,17,18,19]. The topography around the Chheda Khola watershed in the Jajarkot district, Nepal, has been continuously hit by various events of mass movements. Such events have impacted the area with socio-economic losses; to mitigate the impact, potential landslide-prone areas should be identified [20]. A landslide susceptibility map indicates where future landslides are likely to occur by identifying the areas of previous landslide occurrences, and areas with similar or identical physical characteristics [21,22]. It divides a spatial area into zones of varying likelihoods of landslide occurrence [23]. The investigation of landslide contributing factors, and the selection of an effective approach, are critical for the preparation of a landslide susceptibility map. In the area of active tectonics, the identification of the detailed lithological units, precise geological structure placement, along with the preparation of the accurate landslide inventory, play a vital role in such mapping.
The Lesser Himalaya, also known as thrust-and-fold belt, has complex geological settings due to the folding and thrusting of crystalline sequences. The western and the central part of the Nepal Himalaya comprises numerous crystalline bodies, which are core of great Himalayan syncline or anticline. Among well-known crystalline in Nepal Himalaya, two metamorphic cryatalline bodies i.e., Jajarkot Crystalline Zone and Karnali Crystalline Zone, are located along the study area [24,25,26]. The Jajarkot Crystalline Zone consists of sheared phyllite, schists, and blasto-mylonitic augen gneiss [25]. The Karnali Crystalline Zone consists of kyanite sillimanite–bearing gneiss, calc-silicate gneiss, migmatitic gneiss, and augen gneiss [27]. Some of the researchers have intrepreted them as klippen, while some others have interpreted them as nappes [25,26,27]. Although many researchers have worked on the study area, this region still lacks uniform geological and structural interpretation. One of the major interests of this study is to investigate whether the crystalline bodies present in the study area are klippen or nappes. Further, there are no studies carried out on landslide characterization and landslide susceptibility assessment. To address these issues the following studies are carried out:
  • A detailed geological exploration is carried out to generate the precise geological map of the study area.
  • A field survey on landslide occurrence patterns is carried out.
  • The geological and geo-structural control on landslide occurrence are investigated.
  • Landslide susceptibility modeling is conducted using two bivariate statistical approaches, Frequency Ratio (FR) and Weight of Evidence (WOE), to explore and predict the probability distribution of future landslides in the area.

2. Study Area

The study area lies in the southern part of Jajarkot district, Karnali Province, Mid-Western Nepal, between latitude and longitude of 28.677557 N, 82.053348 E, and 28.838401 N, 82.110773 E, covering an area of approximately 190 km2 (Figure 1). It is about 600 km west of Kathmandu, Nepal’s capital, and is accessible by regular air services from Kathmandu to Birendranagar in Surkhet and Nepalgunj in Banke district. From both of these airstrips, there is a well-paved road leading up to Chheda, which is in the southern part of the study area, and a fairweather gravel road connects Chheda with Thalaha. Yet, most of the study areas are hardly accessible due to uneven topography, difficult terrain, and a lack of infrastructure.
The geomorphology of the study area is characterized by highly elevated mountains, extremely rugged topography, and river valleys with uneven relief. Peak, saddle, spur, ridge, valley, river valley, alluvial terrace, and frequent colluvial cones are the common geomorphic features. The topography of the area is highly dissected and rough with elevations ranging from 620 m to 2520 m. Geologically the study area lies in the Lesser Himalayan Crystalline zone consisting of sheared phyllite, schists, blasto-mylonitic, augen gneiss, kyanite-sillimanite–bearing gneiss, calc-silicate gneiss, migmatitic gneiss, and augen gneiss [28,29].

3. Methodology

Initially, a preliminary landslide inventory map was prepared using topographic maps and high-resolution Google Earth images. After that, fieldwork was conducted for detailed geological investigation via multiple routes. Systematic sample collection, in situ permeability testing, landslide characterization, and a validation of the preliminary landslide inventory map was carried out during the field investigation. The collected rock samples were further tested in the lab.
In this study, 13 causative factors were considered for landslide susceptibility mapping. Along with the field investigation of the geology and landslides, the required data for the investigation on causative factors of the landslide was collected from several sources. Furthermore, landslide susceptibility modeling was performed using two bivariate statistical methods: the Weight of Evidence (WOE) and the Frequency Ratio (FR). The results of these two methods were verified by calculating the landslide density and plotting the success rate curve.

3.1. Data Acquisition and Data Processing

The investigation on the causes of past landslides can be used as the source for predicting future landslides [30]. For this, it is necessary to prepare and analyze the landslide contributing factors. The Digital Elevation Model (DEM) of raster format 20 m × 20 m, prepared using the contour data provided by Department of Survey (DOS) [31] has been used to prepare some data base of the predisposing factors (slope, aspect, SPI, TWI). Based on the DEM, topography factors, slope gradient, slope aspect, slope shape, TWI, SPI, and hydrological factors including river and stream proximity data, were derived. Rainfall intensity in the study area was analyzed using the data obtained from the Department of Hydrology and Meteorology (DHM) [32]. Land use map and distances to road map were prepared using digital data from the DOS. Maps related to geology and tectonism were developed from field investigation. All these data layers were transferred into 20 m × 20 m pixel size raster format. The data processing was carried out using GIS.

3.2. Landslide Susceptibility Modeling

Landslide susceptibility modeling can be grouped into four categories: landslide inventories, heuristic, statistical and deterministic approaches, among which statistical methods have been widely applied. Likewise, this study applies two well-known bivariate statistical methods; the WOE and the FR. The WOE method is a quantitative method used to combine datasets used for various purposes [33,34,35]. In this method, weights for each landslide conditioning factor are calculated based on the presence or absence of landslides within the different classes of a causative factor [33,34]. The equations involved are given in Equations (1) and (2):
W + = ln P { F | L } P { F | L ¯ }    
W = ln P { F ¯ | L } P { F ¯ | L ¯ }
where W+ = weight for the presence of landslides, W = weight absence of landslides within a certain class of causative factor map. P (F/L) = conditional probability of F occurring given the presence of L, a bar above a symbol signifies the negation.
The weights can be calculated easily by cross tabulating the observed landslide map with the causative factor map using Equations (3) and (4).
W + = ln N { F L } / N { L } N { F L } ¯ / N { L ¯ }
W = ln N { F ¯ L } / N { L } N { F ¯ L } ¯ / N { L ¯ }
C = W+W is a contrast and gives the probability of occurrence of landslides for the particular area after combining all the conditioning factors. C is positive for the positive correlation and negative for non-correlation.
The frequency ratio is the ratio of the probability of an occurrence to the probability of a non-occurrence of landslides for any given set of attributes [36]. The lower the value of the frequency ratio, the weaker the correlation. The value of the frequency ratio can be estimated by using the approach followed by KC et al. [37]. The equation used is given in Equation (5).
F R = ln L c A c
where Lc is percentage of in particular class; is calculated as:
Lc = N (A|L)/N (L) × 100%,
where N (A|L) = the number of landslide pixels within class A. N (L) = the total number of landslide pixels in the map.
Ac is the percentage of the specific class A of the entire study area; it can be calculated as:
Ac = N (A)/Nt × 100%,
where N (A) = the number of pixels of a specific class. Nt = the total number of pixels of the entire study area on the map.

4. Results

Results can briefly be grouped into two categories: preparation of a new geological map accompanied by detailed geological investigation; and landslide characterization and landslide susceptibility mapping.

4.1. Geology

Based on lithological characteristics, stratigraphic setup, tectonic and structural features, and a consideration of the regional geological picture, three tectonic units have been recognized in the study area: Jajarkot Klippe, Lesser Himalaya, and Karnali Klippe. Although there is not uniformity in using the terms nappe and klippe for these two cryatalline bodies, the term klippe has been used here following the work of Dhital [35]. These three tectonic units consist of an unmetamorphosed to incipiently metamorphosed rock sequence that runs from south to north, with a layer of unmetamorphosed Lesser Himalayan rock sequence between two thrust sheets (Figure 2). Not only do the three units have distinct lithological characteristics, but they also have distinct metamorphic gradation and deformation.

4.1.1. Lithostratigraphy

The Karnali Klippe, assumed to be the eastern end of the Almora-Dadeldhura nappe in Western Nepal [38], is composed of high-grade metamorphic rocks resembling the rocks of the Higher Himalaya. It consists of tourmaline-bearing granitic augen gneisses, coarsely crystalline quartz gneisses, and kyanite-bearing banded gneisses, with garnetiferous schist in the south. The upper Karnali Klippe is composed of coarsely crystalline granitic augen gneiss intercalated with kyanite grade paragneiss with frequent pegmatitic intrusions, while the lower portion consists of kyanite-garnet gneiss, coarsely crystalline quartz-gneiss, and coarsely crystalline quartzites.
Despite the fact that Hagen [39] introduced the term “Jajarkot Nappe”, the current observation concludes that it is klippe. The four-rock formation of Jajarkot Klippe has well-established rock succession within the limits of their tectonic boundaries (Figure 2). They are Marka Formation, Dangri Formation, Chaurjhari Formation, and Thabang Formation, representing older to the younger sequence stratigraphically.
The Marka Formation [40] is comprised mainly of light grey to pale yellow micaceous quartzites with intercalations of blue-grey garnet schists. In the upper part of this formation, near the Thalaha Thrust (Figure 2), basic intrusion and graphitic schist intercalations are more common. The Marka Formation also contains a mappable meta-diamictite succession, the Chhera Diamictite, which has been mapped as a member due to constraint in lateral continuity of the succession. The lower contact of the Chhera Diamictite [40] is rather sharp, but the contact in the Bheri River section appears to be indistinct and truncated; this is inferred as disconformity [40], while the upper contact is transitional. This rock succession is made up of thin-to-medium bedded, gray, crenulated, calcareous meta-diamictite with angular to sub-rounded clasts of gray dolomite, gray gneiss, schists, and pale yellow and pink quartzite metamorphosed to garnet grade. The proportion of clast in meta-diamictite ranges from 1% to 12% and is strongly stretched parallel to foliation.
The Dangri Formation [41] consists of crenulated, grey, and grey-green garnetiferous schists and mica schists that alternate with very thick-bedded, medium to coarse-grained, laminated grey, yellow, and white quartzites. Both psammitic and pelitic schists predominate over quartzites (about 70%); the quartzite-schist transition is gradational. Intercalation of the gritty layer (psammitic schist) and pelitic schists gives rise to the gneissic appearance in this region.
The Chaurjhari Formation [42] is composed of medium-grained, coarsely crystalline muscovite-biotite-quartz schists, crenulated garnetiferous schists, and graphitic schists with light grey micaceous quartzite; the density and size of the garnet have increased in this formation. Garnet porphyroblasts in schists reach up to 4 mm in diameter.
The Thabang Formation [42] unit’s rocks are well represented by thinly to thickly bedded coarsely crystalline white, bluish grey to purple limestone marbles interbedded with garnetiferous mica schists, forming the core of the syncline (Figure 2).
The Lesser Himalayan meta-sediments between these two nappes consist of only one rock unit, the Ranimatta Formation, which is named after the village Ranimatta by Dhital [40] and has been retained for similar rock successions in the study area. The lower part of the Ranimatta Formation is composed of parallel laminated, infrequently ripple marked fine-to-medium grained pale grey, yellow, or white quartzite with medium to coarse-grained meta-basic intrusions. The upper part contains a gray to green-gray phyllitic meta-sandstone succession that is occasionally intercalated with quartzites.

4.1.2. Structures

Within the study area, there are four megascopic structures: the south-dipping Thalaha Thrust (TT), the gently north-dipping Main Central Thrust (MCT), and two large-scale gently plunging, non-cylindrical syncline and anticline fold, Suwa Gad Anticline (SA) and Bhoor Syncline (BS) (Figure 2). The Thalaha Thrust, named after the village Thalaha in the present study, is equivalent to the Mahabharat Thrust (MT) [43] of the Kathmandu Nappe and the Dubung Thrust of the Kahun Klippe [44]. The MCT overrides the Lesser Himalayan rocks and places it over the high-grade metamorphic rocks of the Higher Himalayan crystalline slab. The TT, on the other hand, juxtaposed the Lesser Himalayan meta-sedimentary sequence (footwall) with the Lesser Himalayan crystalline sequence as the Jajarkot Nappe rocks (hanging wall).
These structural discontinuities are delineated based on the presence of outcrop-scaled deformation structures such as parasitic folds trending N–S, crenulations trending E–W, and brittle shear zones, as well as microstructures and attitudes of the beds. Meso and microstructural observations, as well as measurements of asymmetrical boudins, projection of lineation data, and axis measurement of small parasitic folds, revealed top-to-the-south movement throughout the area, indicating southward propagation of the TT and the MCT. Figure 3 shows two structural features exhibiting top-to-the-south movement. Furthermore, multiple shear zones form the linear pattern along the eastern part of the area, from which a shear band is delineated (Figure 2).

4.2. Landslide Occurrence Pattern and Landslide Inventory

The study area is dominated by excavation-induced shallow landslides (Figure 4), with the volume of failure surface ranging from 20 cubic meters to 38,000 cubic meters covering the area from 25 square meters to 7600 square meters. Several deep-seated complex paleo-landslides with volumes exceeding 1000 cubic meters which were shoveling settlements, agricultural lands, and even human life (Figure 4) were discovered. Failures occurred primarily along colluvium, residual soil, and bedrock, with very few landslides occurring at soil-bedrock contacts. Approximately 50% of landslides are initiated in bedrock, while residual soils, colluvium, and soil-bedrock contact account for around 18%, 29%, and 3% of total landslides, respectively. The landslide deposition materials are dominated by coarse gravels (Figure 4), with coarse and fine deposits covering around 89.5% and 10.5% of the total landslide deposits, respectively. According to the USCS soil classification scheme, soil categories GP, GM, and GC account for approximately 63%, 4%, and 23% of total deposits, respectively, while fine soils ML and CL account for 5% of total deposits. In the study area, the translational type of failure appeared to be the most prevalent. Because of the entrainment of water on its way toward the downslope, most of the landslides have converted into debris flows.
A total of 789 landslides, rock falls, and debris flows were observed and recorded; the largest landslides reach about 300 m in length, 100 m in height, and an average of 15 m in depth. The smallest ones are a few meters in length, width, and depth. Out of the total area of 189.697 km2, landslides cover a total area of 1.1369 km2 i.e., 0.599% of the total area. The landslide inventory map is shown in Figure 5.

4.3. Landslide Conditioning Factors

Understanding the relationship between previous landslides and landslide contributing factors is essential for landslide susceptibility mapping. A spatial database of 13 landslide-related factors was created (Table 1, Figure 6. Later, all of the landslide conditioning factors were classified into various categories.
Slope gradient is one of the principal factors affecting landslide occurrence. The slope angle distribution of the topography in the study area ranges from 0 to 69.71°, which was categorized into seven different classes: 0–15°, 15°–25°, 25°–35°, 35°–45°, 45°–55°, 55°–65°, and >65° (Figure 6a). The study revealed that the landslides are most prone to topography with the slope angles ranging from 35° to 45° having a frequency ratio of 1.55 and weightage of 0.61.
The slope aspect controls moisture retention and it is also related to the attitude of bedding of the rock formation, which in turn affects the physical properties of slope material and its susceptibility to failure [45,46]. In the present study, the slope aspect has been classified into nine different classes (Figure 6b), as flat (0), northeast (22.5°–67.5°); east (67.5°–112.5°); southeast (112.5°–157.5°); south (157.5°–202.5°); southwest (202.5°–247.5°); west (247.5°–292.5°); northwest (292.5°–337.5°); and north (337.5°–22.5°) (Figure 7b). The study revealed that the south and southwest facing slopes are more prone to landslide susceptibility (Table 1). Such slopes of the basins are typically anti-dipping slopes, where rainfall and the continuous fall of weathered debris are typical, and the south-facing landscape is typically on the windward side of the monsoon rain.
Slope erosion processes can affect slope stability, mostly in steep terrain. This process concentrates or disperses surface and subsurface water in the landscape due to slope shape, which makes slope shape a landslide causative factor [47]. The spatial distribution of slope curvature has been classified into three types: linear (−0.001–0.001), concave (−0.001), and convex (>0.001) (Figure 6c). The slope with a concave shape, also known as a convergent landform, is usually the least stable because it develops relatively high pore-water pressure, increasing the driving force due to the concentration of water in a certain area, which eventually initiates a slope failure. The study revealed that the concave shape of the slope is most susceptible to landslides. It has a weightage value of 0.36 and a frequency ratio of 1.29.
The Topographic Wetness Index (TWI) is an index value that reflects the tendency of water to accumulate at any point within the catchment and the tendency of gravitational forces to move that water down the slope [48]. The TWI is commonly used to quantify topographic control on hydrological processes [49] and can be calculated as per the equation used by Youssef et al. [50] given in Equation (6):
TWI = ln ( A t a n β )
where A = specific catchment area and β = slope gradient.
TWI of the topography is classified into six classes: 2.05–5.4, 5.4–6.8, 6.8–8.77, 11.57–15.59, and 15.59–23.51 (Figure 6d, Table 1). Landslides are more likely to occur in topographic features with lower TWI values that favor the accumulation of water (Table 1).
The stream power index (SPI) describes the erosive strength of the stream on the slope. It is a measurement of the erosive power of flowing water based on the assumption that discharge is proportional to specific catchment area [50]. It can be calculated as:
SPI = (tan β × As)
where β = slope angle in degree.
As can be estimated by using the approach introduced by Hengl et al. [51] as:
As = (Am × P2Li)
where P = pixel size, Am = sum of drainage fraction from surrounding, ∑Li = total length of drainage pixels.
In the present study, three classes are considered: (1) low erosive; (2) high erosive; and (3) very high erosive (Figure 6e). The study showed that, as per the calculated weights, the very high erosive class comparatively contributes the most to landslides.
Stream proximity is the distance from the location of interest to the nearest stream. The closeness of stream or natural drainage has an important effect on slope stability [52]. For the present study, stream proximity is categorized into four classes: (1) 0–25 m; (2) 25–50 m; (3) 50–100 m; and (4) >100 m (Figure 7f). The landforms close to the stream i.e., within the proximity of 25 m are most susceptible to the events of mass movement (Table 1). Stream proximity appeared to be a major landslide contributing factor with significant weightage. Similarly, the distance from road alignments also appeared to be one of the contributing factors for the landslide occurrences.
The properties of the rock units and geological structure within the study area are important factors in landslide occurrence [52,53]. Depending on the characteristics and types, different rock units have varying attitudes toward landslides. The major rock types in the study area are schist, diamictite-schist, phyllite, and quartzite, with frequent basic intrusions (Figure 7). All the rock types in the study area are highly deformed, fractured, and weathered due to the presence of weathering-sensitive minerals such as micas and feldspar (Figure 7).
The majority of the area is covered by ductile incompetent rocks such as schist, meta-diamictite, and phyllites. They have undergone folding during various stages of Himalayan deformation. Furthermore, even competent rocks, such as quartzite and marble, are highly fractured and brecciated. Quartzite and marble are intercalated with incompetent rocks, such as phyllite and schists. This stratigraphic arrangement is prone to differential weathering. The Thabang Formation, which has interbedded marble and schist, and the Chaurjhari Formation, which has incompetent lithology, are the most prone to landslides (Figure 6g, Table 1).
Furthermore, the presence of a fault decreases the rockmass strength, thereby increasing the landslide susceptibility; faults also play an important role in abnormal groundwater conditions as there is high water seepage through the fault or even significant water flow can occur along fault line which can saturate the slopes and eventually trigger landslides [37,54,55]. Therefore, faults are one of the major landslides contributing factors. A digital thematic distance from faults map was prepared and classified into three classes: (1) 0–500 m; (2) 500–1000 m; and (3) >1000 m (Figure 6h). Table 1 shows that the areas closer to the faults are more susceptible to landslide occurrences.
Geological structures formed by both brittle and ductile deformation break the homogeneity of rock mass and reduce the shear strength. Likewise, the area closer to the fold’s axis is weaker and is also susceptible to landslides. The study area has one large-scale/regional syncline fold axis and one anticline fold axis. Furthermore, many small-scale anticline and syncline folds contribute to the instabilities. A digital thematic distance from syncline and anticline map was created and classified into three categories: close (0–100 m); nearby (100–500 m); and distant (>500 m) (Figure 6i,j).
Furthermore, land use can alter soil microbiological processes and have an impact on slope stability. It is classified into four categories in the current study based on data provided by Nepal’s Department of Survey (DOS): (1) barren lands; (2) cultivation and built-up area; (3) forest; and (4) water bodies and deposits. Figure 6l and Table 1 show that barren lands are most susceptible to landslides.
Lastly, rainfall is a major triggering agent of landslides in mountainous terrain [6,12]. When daily precipitation exceeds 144 mm, the likelihood of a landslide occurring in the Himalayas increases significantly [56]. The annual precipitation in the area ranges from 1600 mm to 2000 mm [32]. This factor is classified into three classes in the current study: 600 mm, 1800 mm, and 2000 mm (Figure 6m). According to the results of the in situ permeability tests, the coefficients of permeability (K), transmissivity (T), and storativity (S) for old colluvium materials are 1.88 × 10−2 m/s, 6.2 × 10−2 m2/sec, and 0.588, respectively. Likewise, the values of K, T, and S for residual soils are 2.73 × 10−4 m/s, 0.091 m2/sec, and 0.3841, respectively, and that for weathered bedrock are 5.62 × 10−5 m/s, 1.68 × 10−3 m2/sec, and 0.7248, respectively. With a significant amount of rainfall, the coefficients of permeability, storativity, and transmissivity of all these material categories are quite significant, and enough water could infiltrate to downslope which has incompetent rocks like schists and intercalation of incompetent and competent rocks like quartzite and schists. This process can lead to differential weathering in rocks and an increment in pore water pressure below the surface, which contributes to landslides.
Road construction can disturb the slopes, increasing the possibility of a landslide [57]. During the field survey, several landslides were observed along the roadside. A digital thematic distance from a road map was prepared for the present study and classified into four classes: (1) 0–100 m; (2) 100–300 m; (3) 300–500 m; and (4) greater than 500 m (Figure 6k). Calculation during landslide susceptibility modeling revealed that the area nearer to the road is more susceptible to landslides.

4.4. Landslide Susceptibility Modeling

The process of statistical approach combined with geo-mechanical understanding can be a powerful tool for landslide susceptibility modeling. The prior probability and the posterior probability are the major principles of the WOE [58]. Calculation of the ratio of the probability of occurrence to non-occurrence of landslides for a given set of an attribute is the major working principle for the FR. In the present research, the FR and WOE methods were used to determine the level of correlation between the location of landslides and landslide-conditioning factors. Prior probability was calculated first by the probability of an event of the same type occurring in the previous period, and then the same probability was used for an event’s future probability. The detailed results of the susceptibility assessment are shown in Table 1, including the spatial relationship between landslide and landslide contributing factor classes for both the FR and the WOE models.
The landslide susceptibility index (LSI) values for the WOE method varies from −1.45 to 13.22 and that for the FR ranges from 201 to 456. There is no general rule for susceptibility zone division. In the present study, the LSI index map from both approaches is categorized into very low, low, moderate, high, and very high susceptibility classes in such a way that each class has nearly equal area distribution (around 20%) based on LSI values. Table 2 shows the percentage distribution of overall area coverage by very low, low, medium, high, and very high susceptibility class and landslide distribution within each class. Figure 8 shows the landslide susceptibility models developed using the FR and the WOE methods.

4.5. Validation

Two basic assumptions are needed for the validation of landslide susceptibility maps: landslides are related to spatial information such as topography and geology, and future landslides will occur by a specific trigger, such as rainfall or an earthquake [59,60]. In the present study, these two necessary conditions are fulfilled as landslides are related to all the spatial information and they are also triggered by rainfall.
In the present case, approximately 80% of the recorded landslides are within high and very high susceptibility zones for both FR and WOE models. To evaluate the overall quality of the landslide susceptibility map, it is also necessary to compute the landslide density of each class [61]. The landslide density for the very high susceptible zone is 3.11 for the FR and 3.18 for the WOE, which is significantly higher than the landslide density for the other susceptible zones. There is a gradual decrease in landslide density from the very high to lower landslide susceptible zones, with significant differences in the values between the different susceptible zones (Table 2). These results are consistent with the field observations.
The validity of the landslide susceptibility map can also be estimated graphically and quantitatively using the success rate curve [59,62]. For the success rate curve, the cumulative percentage of observed landslide occurrence is plotted against the areal cumulative percentage in decreasing LSI values (Figure 9). The area under a curve for the WOE method is 0.7942, while that for the FR method is 0.7762, resulting in an overall success rate of 79.42% and 77.62%, respectively. This result also validates the landslide susceptibility map with the pre-existing slope instability conditions.

5. Discussion

For the present research, we considered the region between Chhedachaur to Saureni, the Southern part of Jajarkot district, Nepal. Despite an exploration of the area by many researchers, uniformity in the geological map of the area was lacking, due to which a detailed geological investigation was conducted as the first step. Moreover, there is still some confusion about the boundary of these two crystalline bodies. As we did not find any traces of the root zone within the study area, we support the concept of Klippe following the work of Dhital [35]. The detailed geodynamic studies give information on the placement of crystalline bodies [63,64,65,66], however there are no studies carried out related to geodynamics and research that incorporate geodetic data with crustal deformation in the study area.
The geological investigation included a detailed field survey, both along and across strike through many routes to delineate the lithological and tectonic boundary, which is later verified by petrographic analysis. The area revealed three tectonostratigraphic successions: the Karnali Klippe, the Lesser Himalaya, and the Jajarkot Klippe, most of which have their own well-defined litho-stratigraphic successions. The litho-stratigraphic subdivision of the rock sequences of the present work can be correlated with the work of Fuchs [67], Sharma et al. [68], and Stocklin [69], who have worked within, nearby, or in similar tectonic settings (Table 3). The outcome of this research is not in agreement with the recently published map of the Department of Mines and Geology, Nepal [70], in terms of structural boundaries and litho-stratigraphy.
In the second step, the spatial distribution pattern of the landslides was carried out. Most of the landslides were observed near the thrust and fold axis, close to roads and stream, within the formation comprising of inter-bedding of competent and incompetent lithology, and within the south and south-west dipping slopes having slope angles of 35°–45°. The most important observation was that fresh shallow landslides were observed along the road alignment, which is most probably due to the haphazard excavation of the road without geological and geotechnical consideration. In addition, the reduction of the shear strength by the development and/or propagation of the geological structures, and the differential weathering of the rock mass, contributed the most to the occurrence of landslides.
An important step to assess the landslide susceptibility is the selection of the predisposing factors, the understanding of their combined effects, and the relative importance of the selected factors on landslide distribution [71]. So, in the third step, landslide susceptibility modeling was carried out employing two different methods i.e., the WOE and the FR, by considering the 13 causative factors of slope angle, slope aspect, slope shape, SPI, TWI, lithology, distance from thrust, distance major syncline, distance from major anticline, stream proximity, distance from the road, land use, and mean annual precipitation. The statistical analysis results showed that the models performed well, as evidenced by the degree of fit and the percentage of occurrences of landslides in each susceptibility class. The high and very high susceptibility classes have a higher degree of fit, which decreases significantly towards the lower susceptibility class. Furthermore, the success rate obtained is good (79.42% for the WOE approach and 77.62% for the FR approach). Other studies also reported the success rate within this proximity, for example Dahal et al. [72] reported a success rating of 77.6% to 80.66% from the WOE method. Kayastha et al. [31] reported a success rate of 79% from the WOE method. KC et al. [37] found 75% and 71% success rates for the WOE and FR methods, respectively, from the research in Arun Tectonic Window, Nepal. Bijukchhen et al. [73] have reported a success rate of around 80% from the heuristic and bivariate statistical methods. Akgun [74] reported the success rate of logistic regression, frequency ratio, and analytical hierarchy process (AHP) methods at 81%, 76%, and 71%, respectively.
Some of the previous studies have concluded that the WOE method produced better results than the FR method [37,75], and the same is true for the current study area. The results show that the FR and WOE models are suitable for this study and can be successfully applied to landslide susceptibility zonation mapping in other mountainous areas with similar topographical and physical features.

6. Conclusions

In this study, detailed geologial investigation is carried out to develop the precise geological map of the southen part of Jajarkot district in between the two crystalline zones i.e., Karnali Crystalline Zone and Jajarkot Crystalline Zone, Nepal. Investigation on whether these crystalline bodies are klippen or nappes is performed. Furthermore, extensive survey on landslide distribution pattern in the area is conducted along with the landslide susceptibility mapping using two bivarate models i.e., Frequency Ratio (FR) and Weight of Evidence (WOE), considering 13 causative factors. The major findings of this study are listed below:
  • From detailed geological exploration, a new geological map of the study area is developed, which can be the reference for further research in this area. This map can be crucial for the planning of infrastructure development.
  • This study supports the concept that the crystalline bodies i.e., Karnali Crystalline Zone and Jajarkot Crystalline Zone, located around the study area are klippen.
  • The field survey shows that the haphazard road excavation without the consideration of geological and geotechnical features of the site has caused rampant shallow landslides.
  • Susceptibility modeling using both the FR and the WOE methods shows that the lithology is the most important conditioning factor for landslide occurrence, while, land use, distance from geological structures, road, SPI, slope, aspect and TWI, have played a moderate role.
  • The overall success rating of the FR and WOE method in the study area are 77.62% and 79.42%, respectively, with the 13 causative factors considered. The development of such susceptibility maps under similar geoenvironmental and topographical features, as in this study, and considering similar causative factors, can be useful for the mitigation of landslide hazards in other places of the world.

Author Contributions

Conceptualization, Y.B.S. and S.K.; methodology, D.K., H.D., Y.B.S., S.K. and G.A.; software, S.K., D.K. and H.D.; formal analysis, Y.B.S., S.K., D.K., H.D. and K.R.P.; investigation, Y.B.S., S.K. and G.A.; resources, K.R.P.; data curation, Y.B.S., S.K., G.A., D.K. and H.D.; writing—original draft preparation, Y.B.S., S.K., G.A., D.K. and H.D. writing—review and editing, D.K., S.K. and H.D.; visualization, Y.B.S., S.K., G.A. and H.D.; supervision, D.K. and K.R.P.; project administration, K.R.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data used in the present study is contained within the article.

Acknowledgments

The authors would like to thank Pawan Kumar Acharya and Aneeta Thapa for their cooperation and support during the field study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaeger, J.; Courtillot, V.; Tapponnier, P. Paleontological view of the ages of the Deccan Traps, the Cretaceous/Tertiary boundary, and the India-Asia collision. Geology 1989, 17, 316–319. [Google Scholar] [CrossRef]
  2. Besse, J.; Courtillot, V.; Pozzi, J.P.; Westphal, M.; Zhou, Y.X. Palaeomagnetic estimates of crustal shortening in the Himalayan thrusts and Zangbo suture. Nature 1984, 311, 621–626. [Google Scholar] [CrossRef]
  3. Copley, A.; Avouac, J.-P.; Royer, J.-Y. India–Asia collision and the Cenozoic slowdown of the Indian plate: Implications for the forces driving plate motions. J. Geophys. Res. 2010, 115, 1–14. [Google Scholar] [CrossRef] [Green Version]
  4. Avouac, J.-P. Mountain building, erosion and the seismic cycle in the Nepal Himalaya. Adv. Geophys. 2003, 46, 1–80. [Google Scholar]
  5. Searle, M.P. Structural evolution and sequence of thrusting in the High Himalayan, Tibetan-Tethys and Indus suture zones of Zanskar and Ladakh, Western Himalaya. J. Struct. Geol. 1986, 8, 923–936. [Google Scholar] [CrossRef]
  6. Pratt-Sitaula, B.; Burbank, D.W.; Heimsath, A.; Ojha, T. Landscape disequilibrium on 1000–10,000 year scales Marsyandi River, Nepal, central Himalaya. Geomorphology 2004, 58, 223–241. [Google Scholar] [CrossRef]
  7. Guzzetti, F. Landslide Hazard and Risk Assessment. Ph.D. Thesis, Rheinischen Friedrich Wilhelms-University Bonn, Bonn, Germany, 2005. [Google Scholar]
  8. CRED. The Human Cost of Natural Disasters 2015—A Global Perspective; CRED; UNISDR: Brussels, Belgium, 2015; p. 55. [Google Scholar]
  9. Herath, S.; Wang, Y. Case studies and national experiences. In Landslides–Disaster Risk Reduction; Springer: Berlin/Heidelberg, Germany, 2009; pp. 475–497. [Google Scholar]
  10. Cruden, D.M. A simple definition of a landslide. Bull. Int. Assoc. Eng. Geol. 1991, 41, 27–29. [Google Scholar] [CrossRef]
  11. Crosta, G.B.; Frattini, P. Rainfall-induced landslides and debris flows. Hydrol. Process. 2008, 22, 473–477. [Google Scholar] [CrossRef]
  12. Cannon, S.H.; Bigio, E.R.; Mine, E. A process for fire-related debris flow initiation, Cerro Grande fire, New Mexico. Hydrol. Process. 2001, 15, 3011–3023. [Google Scholar] [CrossRef]
  13. Fannin, R.J.; Rollerson, T.P. Debris flows: Some physical characteristics and behaviour. Can. Geotech. J. 1993, 30, 71–81. [Google Scholar] [CrossRef]
  14. Iverson, R.M.; Reid, M.E.; LaHusen, R.G. Debris flow mobilization from landslides. Annu. Rev. Earth Planet Sci. 1997, 25, 85–138. [Google Scholar] [CrossRef]
  15. Dangi, H.; Bhattarai, T.N.; Thapa, P.B. An approach of preparing earthquake induced landslide hazard map: A case study of Nuwakot District, central Nepal. J. Nepal Geol. Soc. 2019, 58, 153–162. [Google Scholar] [CrossRef]
  16. Dahal, B.K.; Dahal, R.K. Landslide hazard map: Tool for optimization of low-cost mitigation. Geoenviron. Disasters 2017, 4, 8. [Google Scholar] [CrossRef] [Green Version]
  17. KC, D.; Dangi, H.; Naqvi, M.W.; Hu, L. Landslide Mobilized Debris Flow at Kalli Village in Achham, Nepal: A Case Study. In Proceedings of the 2021 ASEE North Central Section Conference, Toledo, OH, USA, 19 March 2021. [Google Scholar]
  18. Joshi, I.; Chaudhary, R.K.; Diwakar, K.C. Engineering Perspective of Melamchi Disaster and Mitigation Methods. Zenodo 2021, 1–4. [Google Scholar] [CrossRef]
  19. KC, D.; Gautam, K.; Dangi, H.; Kadel, S.; Hu, L. Challenges in Tunneling in the Himalayas: A Survey of Several Prominent Excavation Projects in the Himalayan Mountain Range of South Asia. Geotechnics 2022, 2, 802–824. [Google Scholar] [CrossRef]
  20. Devkota, K.C.; Regmi, A.D.; Pourghasemi, H.R.; Yoshida, K.; Pradhan, B.; Ryu, I.C.; Althuwaynee, O.F. Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling–Narayanghat road section in Nepal Himalaya. Nat. Hazards 2013, 65, 135–165. [Google Scholar] [CrossRef]
  21. Van Westen, C.J.; Van Asch, T.W.; Soeters, R. Landslide hazard and risk zonation—Why is it still so difficult? Bull. Eng. Geol. Environ. 2006, 65, 167–184. [Google Scholar] [CrossRef]
  22. Kayastha, P.; Dhital, M.R.; De Smedt, F. Landslide susceptibility mapping using the weight of evidence method in the Tinau watershed, Nepal. Nat. Hazards 2012, 63, 479–498. [Google Scholar] [CrossRef]
  23. Anbalagan, R. Landslide hazard evaluation and zonation mapping in mountainous terrain. Eng. Geol. 1992, 32, 269–277. [Google Scholar] [CrossRef]
  24. Fuchs, G.; Frank, W. The geology of west Nepal between the rivers Kali Gandakiand Thulo Bheri. Geol. Bundesanst. 1970, 18, 1–103. [Google Scholar]
  25. Arita, K.; Shiraishi, K.; Hayashi, D. Geology of western Nepal and a comparison with Kumaun, India. Bull. Fac. Sci. Hokkaido Univ. 1984, 21, 1–20. [Google Scholar]
  26. Dhital, M.R. Juxtaposition of Greater and Lesser Himalayan nappes in west Nepal: Implications for delineating Main Central Thrust. Himal. Geol. 2022, 43, 231–240. [Google Scholar]
  27. La Roche, R.S.; Godin, L.; Crowley, J.L. Reappraisal of emplacement models for Himalayan external crystalline nappes: The Jajarkot klippe, western Nepal. Bull. Geol. Soc. Am. 2018, 130, 1041–1056. [Google Scholar] [CrossRef]
  28. Yoneshiro, T.; Kizaki, K. Metamorphism and thermal structure of the Karnali Klippe, western Nepal Himalaya. Bull. Dep. Geol. 1996, 5, 1–20. [Google Scholar]
  29. Hayashi, D.; Fujii, Y.; Yoneshiro, T.; Kizaki, K. Observations on the geology of the Karnali region, west Nepal. J. Nepal Geol. Soc. 1984, 4, 29–40. [Google Scholar]
  30. Wang, H.B.; Sassa, K. Comparative evaluation of landslide susceptibility in Minamata area, Japan. Environ. Geol. 2005, 47, 956–966. [Google Scholar] [CrossRef]
  31. DOS. Topographic Maps of Nepal; Department of Survey: Kathmandu, Nepal, 1996. [Google Scholar]
  32. DHM. Precipitation Map of Nepal; Department of Hydrology and Meterology: Kathmandu, Nepal, 2021. [Google Scholar]
  33. Bonham-Carter, G.F.; Agterberg, F.P.; Wright, D.F. Integration of geological datasets for gold exploration in Nova Scotia. Photogramm. Eng. Remote Sens. 1988, 54, 1585–1592. [Google Scholar]
  34. Bonham-Carter, G.F. Weights of evidence modelling: A new approach to mapping mineral potential. In Statistical Applications in the Earth Sciences. Geol. Surv. Can. 1989, 8–9, 171–183. [Google Scholar]
  35. Lee, S.; Ryu, J.H.; Won, J.S.; Park, H.J. Determination and application of the weights for landslide susceptibility mapping using an artificial neural network. Eng. Geol. 2004, 71, 289–302. [Google Scholar] [CrossRef]
  36. Bonham-Carter, G.F. Geographic Information Systems for Geoscientists: Modeling with GIS; Pergamon Press: New York, NY, USA, 1994; p. 398. [Google Scholar]
  37. KC, D.; Dangi, H.; Hu, L. Assessing Landslide Susceptibility in the Northern Stretch of Arun Tectonic Window, Nepal. Civ. Eng. 2022, 3, 525–540. [Google Scholar] [CrossRef]
  38. Kansakar, D.R.; Chitrakar, G.R. Report on Geology of the Bheri River Section. Dep. Mines Geol. Kathmandu 1984. Unpublished. [Google Scholar]
  39. Hagen, T. Report on the Geological survey of Nepal, Vol. 1: Preliminary reconnaissance. Denkschriften der Schweizerischen Naturforschenden Gesellschaft Memoires de la Societe Helvetique des Sciences Naturelles. 1969; Volume 84, p. 185. [Google Scholar]
  40. Dhital, M.R. Geology of the Nepal Himalaya: Regional Perspective of the Classic Collidedorogen; Springer: Berlin/Heidelberg, Germany, 2015; p. 487. [Google Scholar]
  41. Dhital, M.R.; Kizaki, K. Lithology and Stratigraphy of the Northern Dang, Lesser Himalaya. Bull. Coll. Sci. Univ. Ryukyus Okinawa 1987, 45, 183–244. [Google Scholar]
  42. Adhikary, T.P.; Sharma, T. Geological Report of Marma, Chaurjhari, and Jajarkot Area (Far Western Nepal). Dep. Mines Geol. Kathmandu 1983. Unpublished. [Google Scholar]
  43. Stöcklin, J.; Bhattarai, K.D. Geology of the Kathmandu Area and Central Mahabharat Range, Nepal Himalaya. Dep. Mines Geol. UNDP 1977. Unpublished. [Google Scholar]
  44. Paudyal, K.R.; Paudel, L.P. Geological study and root zone interpretation of the Kahun Klippe, Tanahun, and central Nepal. Himal. Geol. 2013, 34, 93–106. [Google Scholar]
  45. Dai, F.C.; Lee, C.F.; Li, J.X.Z.W.; Xu, Z.W. Assessment of landslide susceptibility on the natural terrain of Lantau Island, Hong Kong. Environ. Geol. 2001, 40, 381–391. [Google Scholar]
  46. Ghimire, M. Landslide occurrence and its relation with terrain factors in the Siwalik Hills, Nepal: Case study of susceptibility assessment in three basins. Nat. Hazards 2011, 56, 299–320. [Google Scholar] [CrossRef]
  47. Nefeslioglu, H.A.; Sezer, E.B.R.U.; Gokceoglu, C.; Bozkir, A.S.; Duman, T.Y. Assessment of landslide susceptibility by decision trees in the metropolitan area of Istanbul, Turkey. Math. Probl. Eng. 2010, 2010, 901095. [Google Scholar] [CrossRef] [Green Version]
  48. Moore, I.D.; Grayson, R.B. Terrain-based catchment partitioning and runoff prediction using vector elevation data. Water Resour. Res. 1991, 27, 1177–1191. [Google Scholar] [CrossRef]
  49. Sörensen, R.; Zinko, U.; Seibert, J. On the calculation of the topographic wetness index: Evaluation of different methods based on field observations. Hydrol. Earth Syst. Sci. 2006, 10, 101–112. [Google Scholar] [CrossRef] [Green Version]
  50. Youssef, A.M.; Pourghasemi, H.R.; El-Haddad, B.A.; Dhahry, B.K. Landslide susceptibility maps using different probabilistic and bivariate statistical models and comparison of their performance at Wadi Itwad Basin, Asir Region, Saudi Arabia. Bull. Eng. Geol. Environ. 2015, 75, 63–87. [Google Scholar] [CrossRef]
  51. Hengl, T.; Gruber, S.; Shrestha, D.P. Digital Terrain Analysis in ILWIS; International Institute for Geo-Information Science and Earth Observation Enschede: Enschede, The Netherlands, 2003; p. 56. [Google Scholar]
  52. Duman, T.Y.; Çan, T.; Emre, O.; Keçer, M.; Doğan, A.; Ateş, Ş.; Durmaz, S. Landslide inventory of northwestern Anatolia, Turkey. Eng. Geol. 2005, 77, 99–114. [Google Scholar] [CrossRef]
  53. Yalcin, A. An Investigation on Ardesen (Rize) Region on the Basis of Landslide Susceptibility. Ph.D. Thesis, Karadeniz Technical University, Trabzon, Turkey, 2005. [Google Scholar]
  54. Leir, M.; Michell, A.; Ramsay, S. Regional landslide hazard susceptibility mapping for pipelines in British Columbia. Geo-engineering for the society and its environment. In Proceedings of the 57th Canadian Geotechnical Conference and the 5th Joint CGS-IAH Conference, Waterloo, ON, Canada, 24–27 October 2004; pp. 1–9. [Google Scholar]
  55. Chen, Y.; Yu, J.; Shahbaz, K.; Xevi, E. A GIS-based sensitivity analysis of multi-criteria weights. In Proceedings of the 18th World IMACS/MODSIM Congress, Cairns, Australia, 13–17 July 2009. [Google Scholar]
  56. Dahal, R.K.; Hasegawa, S. Representative rainfall thresholds for landslides in the Nepal Himalaya. Geomorphology 2008, 100, 429–443. [Google Scholar] [CrossRef]
  57. Pradhan, B. Landslide susceptibility mapping of a catchment area using frequency ratio, fuzzy logic and multivariate logistic regression approaches. J. Indian Soc. Remote Sens. 2010, 38, 301–320. [Google Scholar] [CrossRef]
  58. Ding, N.; Melloni, L.; Tian, X.; Poeppel, D. Rule-based and word-level statistics-based processing of language: Insights from neuroscience. Lang. Cogn. Neurosci. 2017, 32, 570–575. [Google Scholar] [CrossRef] [Green Version]
  59. Chung, C.J.F.; Fabbri, A.G. Probabilistic prediction models for landslide hazard mapping. Photogramm. Eng. Remote Sens. 1999, 65, 1389–1399. [Google Scholar]
  60. Lee, S.; Pradhan, B. Landslide hazard mapping at Selangor, Malaysia using frequency ratio and logistic regression models. Landslides 2007, 4, 33–41. [Google Scholar] [CrossRef]
  61. Sarkar, S.; Kanungo, D.P. An integrated approach for landslide susceptibility mapping using remote sensing and GIS. Photogramm. Eng. Remote Sens. 2004, 70, 617625. [Google Scholar] [CrossRef]
  62. Van Westen, C.J.; Rengers, N.; Soeters, R. Use of geomorphological information in indirect landslide susceptibility assessment. Nat. Hazards 2003, 30, 399–419. [Google Scholar] [CrossRef]
  63. Hashimoto, M. Plate tectonics and crustal deformation around the Japanese Islands. J. Geophys. Res. 1993, 98, 16149–16166. [Google Scholar] [CrossRef]
  64. Müller, M.D.; Geiger, A.; Kahle, H.G.; Veis, G.; Billiris, H.; Paradissis, D.; Felekis, S. Velocity and deformation fields in the North Aegean domain, Greece, and implications for fault kinematics, derived from GPS data 1993–2009. Tectonophysics 2013, 597, 34–49. [Google Scholar] [CrossRef]
  65. Wu, Y.; Li, L.; Chen, C.; Liang, H.; Guo, N.; Li, Y. GNSS deformation characteristics of North China in the past two decades. Geod. Geodyn. 2021, 12, 392–398. [Google Scholar] [CrossRef]
  66. Lazos, I.; Papanikolaou, I.; Sboras, S.; Foumelis, M.; Pikridas, C. Geodetic Upper Crust Deformation Based on Primary GNSS and INSAR Data in the Strymon Basin, Northern Greece—Correlation with Active Faults. Appl. Sci. 2022, 12, 9391. [Google Scholar] [CrossRef]
  67. Fuchs, G. On the geology of the Karnali and Dolpo regions, west Nepal. Mitt. Geol. Ges. Wien 1974, 66, 21–35. [Google Scholar]
  68. Sharma, T.; Kansakar, D.R.; Kizaki, K. Geology and tectonics of the region between Kali Gandaki and Bheri Rivers in central west Nepal. Bull. Coll. Sci. Univ. Ryukus 1984, 38, 57–102. [Google Scholar]
  69. Stöcklin, J. Geology of Nepal and its regional frame. J. Geol. Soc. 1980, 137, 1–34. [Google Scholar] [CrossRef]
  70. DMG. Geological Map of Karnali Province; Department of mines and Geology: Kathmandu, Nepal, 2021. [Google Scholar]
  71. Conforti, M.; Ietto, F. Modeling shallow landslide susceptibility and assessment of the relative importance of predisposing factors, through a GIS-based statistical analysis. Geosciences 2021, 11, 333. [Google Scholar] [CrossRef]
  72. Dahal, R.K.; Hasegawa, S.; Nonomura, A.; Yamanaka, M.; Dhakal, S. DEM-based deterministic landslide hazard analysis in the Lesser Himalaya of Nepal. Georisk 2008, 2, 161–178. [Google Scholar] [CrossRef]
  73. Bijukchhen, S.M.; Kayastha, P.; Dhital, M.R. A comparative evaluation of heuristic and bivariate statistical modelling for landslide susceptibility mappings in Ghurmi-Dhad Khola, east Nepal. Arab. J. Geosci. 2013, 6, 2727–2743. [Google Scholar] [CrossRef]
  74. Akgun, A. A comparison of landslide susceptibility maps produced by logistic regression, multi-criteria decision, and likelihood ratio methods: A case study at Izmir, Turkey. Landslides 2012, 9, 93–106. [Google Scholar] [CrossRef]
  75. Pradhan, B.; Youssef, A.M. Manifestation of remote sensing data and GIS on landslide hazard analysis using spatial-based statistical models. Arab. J. Geosci. 2010, 3, 319–326. [Google Scholar] [CrossRef]
Figure 1. Location map of the study area: (a) Jajarkot district in Nepal around neighboring countries; (b) location of the study area within Jajarkot district and bordering districts of Nepal; (c) the study area with major streams.
Figure 1. Location map of the study area: (a) Jajarkot district in Nepal around neighboring countries; (b) location of the study area within Jajarkot district and bordering districts of Nepal; (c) the study area with major streams.
Geotechnics 02 00050 g001
Figure 2. Geological map of the study area representing the region between Chhedachaur to Saureni.
Figure 2. Geological map of the study area representing the region between Chhedachaur to Saureni.
Geotechnics 02 00050 g002
Figure 3. Meso-scale geological structures showing top-to-the-south movement: (a) S-type folded vein exposed in Dasera section; (b) Z-type drag fold exposed in Karkigaun section.
Figure 3. Meso-scale geological structures showing top-to-the-south movement: (a) S-type folded vein exposed in Dasera section; (b) Z-type drag fold exposed in Karkigaun section.
Geotechnics 02 00050 g003
Figure 4. Some photographic references of the landslides reported within the study area: (a) excavation-induced landslides with multiple scarps, (bd) are landslides that are converted into debris flows; (e,f) are Google Earth images of the Bhoor Landslide and Kaptola landslide, which are deep-seated complex paleo-landslides.
Figure 4. Some photographic references of the landslides reported within the study area: (a) excavation-induced landslides with multiple scarps, (bd) are landslides that are converted into debris flows; (e,f) are Google Earth images of the Bhoor Landslide and Kaptola landslide, which are deep-seated complex paleo-landslides.
Geotechnics 02 00050 g004aGeotechnics 02 00050 g004b
Figure 5. Landslide inventory map of the study area showing elevation ranges.
Figure 5. Landslide inventory map of the study area showing elevation ranges.
Geotechnics 02 00050 g005
Figure 6. Landslide conditioning factor maps of the study area: (a) slope angle; (b) slope aspect; (c) slope shape, (d) TWI; (e) SPI; (f) stream proximity; (g) lithology (h) distance from thrust; (i) distance from major anticline; (j) distance from major syncline; (k) distance from the road; (l) land use; and (m) precipitation.
Figure 6. Landslide conditioning factor maps of the study area: (a) slope angle; (b) slope aspect; (c) slope shape, (d) TWI; (e) SPI; (f) stream proximity; (g) lithology (h) distance from thrust; (i) distance from major anticline; (j) distance from major syncline; (k) distance from the road; (l) land use; and (m) precipitation.
Geotechnics 02 00050 g006aGeotechnics 02 00050 g006bGeotechnics 02 00050 g006cGeotechnics 02 00050 g006d
Figure 7. Photomicrographic reference of representative rock samples from the study area: (a) schist; (b) meta-diamictite; (c) garnetiferous schist; (d) micaceous quartzite; (e) amphibolite (meta-basic rocks); and (f) siliceous marble. plg = plagioclase feldspar.
Figure 7. Photomicrographic reference of representative rock samples from the study area: (a) schist; (b) meta-diamictite; (c) garnetiferous schist; (d) micaceous quartzite; (e) amphibolite (meta-basic rocks); and (f) siliceous marble. plg = plagioclase feldspar.
Geotechnics 02 00050 g007
Figure 8. Landslide susceptibility models prepared using: (a) WOE method, (b) FR method.
Figure 8. Landslide susceptibility models prepared using: (a) WOE method, (b) FR method.
Geotechnics 02 00050 g008
Figure 9. Success rate curve showing cumulative percentage of observed landslide occurrences versus cumulative percentage of decreasing landslide susceptibility index value.
Figure 9. Success rate curve showing cumulative percentage of observed landslide occurrences versus cumulative percentage of decreasing landslide susceptibility index value.
Geotechnics 02 00050 g009
Table 1. The spatial relationship between each landslide-related factor and landslide locations using WOE and FR.
Table 1. The spatial relationship between each landslide-related factor and landslide locations using WOE and FR.
FactorsWeight of EvidenceFrequency Ratio
Total Pixel CountLandslide Pixel CountW+WCS2 (W+)S2 (W)SCC/SC% Class Pixel% Landslide PixelFRRF
Slope
0–1537,767123−0.6080.039−0.650.00810.000370.092−7.007.964.350.550.09
15–25127,475489−0.4440.124−0.570.00200.000430.050−11.3926.8817.290.640.11
25–35172,7201002−0.0280.016−0.040.00100.000550.039−1.1136.4235.420.970.16
35–45103,7659590.441−0.1680.610.00100.000530.04015.2821.8833.901.550.26
45–5528,7952400.337−0.0260.360.00420.000380.0685.366.078.481.400.23
55–65358916−0.2930.002−0.290.06270.000350.251−1.170.760.570.750.12
>651300-0.000--0.00035--0.030.000.000.00
Aspect
Flat974829−0.6990.011−0.710.03450.000350.187−3.792.061.030.500.06
N 47,57391−1.1410.073−1.210.01100.000360.107−11.3910.033.220.320.04
NE 54,616144−0.8200.0055−0.890.00690.000370.086−10.3911.525.090.440.05
E 59,8914010.116−0.0180.130.00250.000410.0542.4812.6314.171.12 0.13
SE 63,4895100.300−0.0550.360.00190.000430.0497.2313.3918.031.350.16
South 62,6715490.387−0.0740.460.00180.000440.0489.6713.2219.411.460.16
SW 66,0225790.388−0.0800.470.00170.000440.04710.0013.9220.471.470.17
W 64,5273960.029−0.0050.0330.00250.000410.0540.6113.6114.001.030.12
NW 45,704130−0.7440.055−0.800.00770.000370.090−8.889.644.600.480.06
Curvature
<0118,9579210.263−0.1060.370.00100.000520.0409.1425.0832.561.300.43
0119,4617210.012−0.0040.020.00140.000470.0430.3625.1925.491.010.34
>0235,8231187−0.1710.145−0.320.00080.000610.038−8.2549.7341.960.840.28
SPI
Low erosive152,140583−0.439−0.226−0.210.00170.000440.047−4.5732.0820.610.640.21
high erosive157,321792−0.165−0.3240.160.00120.000490.0423.8033.1728.000.840.28
very high erosive164,78014540.400−0.7191.120.00060.000720.03829.6634.7551.401.480.49
TWI
2.05–5.4248,94112140.822−0.5571.380.00080.000620.03836.2252.4942.910.820.15
5.4–6.8139,7989991.207−0.4321.640.00100.000540.03941.5229.4835.311.200.22
6.8–8.7755,6914561.384−0.1711.550.00220.000420.05130.3011.7416.121.370.25
8.77–11.5721,4711281.005−0.0411.050.00780.000370.09111.534.534.521.000.18
11.57–15.596495280.726−0.0040.730.03580.000350.1903.831.370.990.720.13
15.59–23.51184540.3640.0050.360.25050.000350.500.720.390.140.360.07
Lithology
Badakada Formation72,664104−1.427−0.032−1.400.00960.000370.100−13.9515.383.690.240.02
Baskot Marble23,7437171.651−0.2891.940.00140.000470.04444.335.0225.415.060.50
Batule Formation114,9747700.121−0.3140.440.00130.000480.04210.2824.3327.291.120.11
Chhera Diamictite59,48568−1.652−0.019−1.630.01470.000360.123−13.3012.592.410.190.02
Karkigaun Schist19,3912600.822−0.0910.910.00390.000390.06613.954.109.212.250.22
Karnali Klippe41,79483−1.099−0.024−1.070.01200.000360.112−9.648.842.940.330.03
Marka Formation 140,523820−0.017−0.3390.320.00120.000500.0427.7429.7429.060.980.04
Land Use
Barren land52611141.312−0.0241.340.00890.000370.09713.831.111.040.940.38
Cultivation and Buildup136,610546−0.3960.132−0.530.00180.000440.048−11.0728.815.010.170.07
Forest323,20020040.045−0.0830.130.00050.001210.0413.0768.1518.380.270.11
Water Bodies91591651.123−0.0351.160.00610.000370.08114.311.931.510.780.32
Thrust (m)
0–50026,2312450.456−0.0340.490.00410.000380.0677.295.518.861.570.52
500–100026,825149−0.0720.004−0.080.00670.000370.084−0.905.665.270.930.31
>1000421,2852435−0.0320.222−0.250.00040.002550.054−4.6688.8386.070.970.32
Syncline
Close28,99957−1.106−1.014−1.090.01750.000360.1348.156.102.010.330.07
Nearby28,0094891.091−0.1851.280.00200.000420.050−25.475.9117.292.930.64
Distant417,2932283−0.081−1.6441.560.00040.001830.048−32.7787.9980.700.920.20
Anticline
Close97811370.868−0.0440.910.00740.000370.088−10.362.064.842.350.54
Nearby9729960.513−0.0290.540.01050.000360.104−5.192.053.391.650.38
Distant454,7312596−0.038−2.4962.460.00030.004290.068−35.9295.8991.760.960.22
Road (m)
0–100122,89113640.626−0.3600.990.00070.000680.03826.1026.9148.211.860.47
100–300132,958721−0.0960.035−0.130.00130.000470.043−3.30328.0425.490.910.23
300–50086,129407−0.2340.045−0.280.00240.000410.054−5.2118.1614.390.790.20
>500132,263337−0.8540.201−1.060.00290.000400.058−18.1627.8911.910.430.11
Stream (m)
0–2578,2406130.280−0.2400.520.00160.000450.04611.3516.5021.671.310.32
25–5062,0863740.016−0.1370.150.00260.000400.0562.7413.0913.221.010.24
50–100109,003586−0.099−0.2270.130.00170.000440.0472.7722.9820.710.900.22
>100224,9121256−0.060−0.5840.520.00080.000630.03813.8047.4344.400.940.23
Precipitation (mm/yr)
1600 52,34094−1.199−0.028−1.170.01060.000360.105−11.1511.033.320.300.13
1800328,43722590.149−1.6011.750.00040.001750.40737.2969.2679.851.150.50
200093,464476−1.153−1.0790.030.00210.000420.0500.5219.7016.830.850.37
Table 2. Observed landslide density in the different susceptibility zones of the landslide susceptibility zonation map.
Table 2. Observed landslide density in the different susceptibility zones of the landslide susceptibility zonation map.
MethodSusceptibility
Zones
Total Area
(sq. km)
Total
Area%
Landslide
Area in Zones
(sq. km)
Landslide
Area%
Landslide
Density
Degree of
Fit (%)
FRVery low38.5120.370.032.800.142.70
Low39.0320.650.087.510.367.14
Medium37.3019.730.1110.100.5110.05
High37.5019.840.2219.170.9718.98
Very High36.6819.410.6860.423.1161.13
WOEVery low36.1119.110.043.120.163.36
Low38.3220.280.065.070.255.14
Medium36.6019.370.119.390.489.97
High38.4820.360.1815.800.7815.96
Very High39.4720.890.7566.623.1965.58
Table 3. Litho-stratigraphic correlation of rock succession of the present study area with other existing literature.
Table 3. Litho-stratigraphic correlation of rock succession of the present study area with other existing literature.
Present ObservationHagen (1969) [39]Fuchs (1974) [67]Sharma et al. (1984) [68]Stocklin (1980) [69]Dhital (2015) [40]
Karnali KlippeHigher Himalayan Rocks Kathmandu NappeUpper Crystalline Nappe
MCT---
Lesser HimalayaRanimatta FormationHiunchuli ZoneChail Nappe Ranimatta Formation
--Thalaha Thrust (TT)------
Jajarkot KlippeThabang FormationJajarkot Nappe 1Lower Crystalline NappeChaurjhari GroupThabang FormationBhimphedi GroupBhainseddoban MarbleHigher Himalayan Crystallines
Chaurjhari FormationChaurjhari FormationRaduwa Formation
Dangri FormationHiunchuli Zone and Piuthan Zone (?)Chail Nappe-Marma Khola Thrust-
Kuncha Group
---Mahabharat Thrust-------MCT---
Sharda GroupDangri Formation
----Kapurkot Thrust----
Marka Formation with Chhera Diamictite member Dailekh GroupMarka Formation
Chhera Diamictite
Ranimatta Formation
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shahi, Y.B.; Kadel, S.; Dangi, H.; Adhikari, G.; KC, D.; Paudyal, K.R. Geological Exploration, Landslide Characterization and Susceptibility Mapping at the Boundary between Two Crystalline Bodies in Jajarkot, Nepal. Geotechnics 2022, 2, 1059-1083. https://doi.org/10.3390/geotechnics2040050

AMA Style

Shahi YB, Kadel S, Dangi H, Adhikari G, KC D, Paudyal KR. Geological Exploration, Landslide Characterization and Susceptibility Mapping at the Boundary between Two Crystalline Bodies in Jajarkot, Nepal. Geotechnics. 2022; 2(4):1059-1083. https://doi.org/10.3390/geotechnics2040050

Chicago/Turabian Style

Shahi, Yubraj Bikram, Sushma Kadel, Harish Dangi, Ganesh Adhikari, Diwakar KC, and Kabi Raj Paudyal. 2022. "Geological Exploration, Landslide Characterization and Susceptibility Mapping at the Boundary between Two Crystalline Bodies in Jajarkot, Nepal" Geotechnics 2, no. 4: 1059-1083. https://doi.org/10.3390/geotechnics2040050

APA Style

Shahi, Y. B., Kadel, S., Dangi, H., Adhikari, G., KC, D., & Paudyal, K. R. (2022). Geological Exploration, Landslide Characterization and Susceptibility Mapping at the Boundary between Two Crystalline Bodies in Jajarkot, Nepal. Geotechnics, 2(4), 1059-1083. https://doi.org/10.3390/geotechnics2040050

Article Metrics

Back to TopTop