Landslide susceptibility mapping of vicinity of Yaka Landslide (Gelendost,Turkey) using conditional probability approach in GIS

On 19 February 2007, a landslide occurred on the Alaard?ç Slope, located 1.6 km south of the town of Yaka (Gelendost, Turkey.) Subsequently, the displaced materials transformed into a mud flow in E?lence Creek and continued 750 m downstream towards the town of Yaka. The mass poised for motion in the Yaka Landslide source area and its vicinity, which would be triggered to a kinetic state by trigger factors such as heavy or sustained rainfall and/or snowmelt, poise a danger in the form of loss of life and property to Yaka with its population of 3,000. This study was undertaken to construct a susceptibility mapping of the vicinity of the Yaka Landslide’s source area and to relate it to movement of the landslide mass with the goal of prevention or mitigation of loss of life and property. The landslide susceptibility map was formulated by designating the relationship of the effecting factors that cause landslides such as lithology, gradient, slope aspect, elevation, topographical moisture index, and stream power index to the landslide map, as determined by analysis of the terrain, through the implementation of the conditional probability method. It was determined that the surface area of the Goksogut formation, which has attained lithological characteristics of clayey limestone with a broken and separated base and where area landslides occur, possesses an elevation of 1,100–1,300 m, a slope gradient of 15 °–35 ° and a slope aspect between 0 °–67.5 ° and 157 °–247 °. Loss of life and property may be avoided by the construction of structures to check the debris mass in E?lence Creek, the cleaning of the canal which passes through Yaka, the broadening of the canal’s base area, elevating the protective edges along the canal and the establishment of a protective zone at least 10-m wide on each side of the canal to deter against damage from probable landslide occurrence and mud flow.     

The environmental impact of landslides on the population living on the eastern footslopes of the Aberdare ranges in Kenya: a case study of Maringa Village landslide

 Landslides are a recognized but poorly studied phenomena in the eastern footslopes of the Aberdare ranges in central Kenya. The Aberdare ranges traverse the Murang'a, Kiambu, Nyeri and Nyandarua districts of Kenya where fatal landslides have occurred in the past. A detailed study of a major landslide which occurred at Maringa village in Murang'a district on 30 April 1997, has been undertaken. The landslide buried three houses at the toe of a slope and killed all the 11 occupants in their sleep. It is shown that the geology and climate of the study area contribute to the prevalence of landslides not only in Murang'a but in all the andosol-covered areas throughout the eastern footslopes of Aberdare ranges. Rapid weathering of pyroclastic rocks in a warm and wet climate create a regolith which is generally weaker than the underlying better-cemented basalt and basaltic agglomerates. The landsliding occurred when the weathered pyroclastic regolith become saturated after a heavy rainstorm on high mechanically unstable slope which slid over the more stable basalts. Over-saturation of clay soils (andosol) which were also derived from weathering of pyroclastic rocks contributed to the slope failure. The results of the study show that besides fatalities and injuries, the landslide destroyed homes, plantations and roads. Received: 16 March 1998 · Accepted: 22 September 1998     

An approach for GIS-based statistical landslide susceptibility zonation—with a case study in the Himalayas

Landslide susceptibility zonation (LSZ) is necessary for disaster management and planning development activities in mountainous regions. A number of methods, viz. landslide distribution, qualitative, statistical and distribution-free analyses have been used for the LSZ studies and they are again briefly reviewed here. In this work, two methods, the Information Value (InfoVal) and the Landslide Nominal Susceptibility Factor (LNSF) methods that are based on bivariate statistical analysis have been applied for LSZ mapping in a part of the Himalayas. Relevant thematic maps representing various factors (e.g., slope, aspect, relative relief, lithology, buffer zones along thrusts, faults and lineaments, drainage density and landcover) that are related to landslide activity, have been generated using remote sensing and GIS techniques. The LSZ derived from the LNSF method, has been compared with that produced from the InfoVal method and the result shows a more realistic LSZ map from the LNSF method which appears to conform to the heterogeneity of the terrain.     

Landslide susceptibility zonation mapping and its validation in part of Garhwal Lesser Himalaya, India, using binary logistic regression analysis and receiver operating characteristic curve method

A landslide susceptibility zonation (LSZ) map helps to understand the spatial distribution of slope failure probability in an area and hence it is useful for effective landslide hazard mitigation measures. Such maps can be generated using qualitative or quantitative approaches. The present study is an attempt to utilise a multivariate statistical method called binary logistic regression (BLR) analysis for LSZ mapping in part of the Garhwal Lesser Himalaya, India, lying close to the Main Boundary Thrust (MBT). This method gives the freedom to use categorical and continuous predictor variables together in a regression analysis. Geographic Information System has been used for preparing the database on causal factors of slope instability and landslide locations as well as for carrying out the spatial modelling of landslide susceptibility. A forward stepwise logistic regression analysis using maximum likelihood estimation method has been used in the regression. The constant and the coefficients of the predictor variables retained by the regression model have been used to calculate the probability of slope failure for the entire study area. The predictive logistic regression model has been validated by receiver operating characteristic curve analysis, which has given 91.7% accuracy for the developed BLR model.     

Understanding rainfall-landslide relationships in man-modified environments: a case-history from Caramanico Terme, Italy

 The expansion of Caramanico Terme in this century has led to the urbanization of marginally stable valley slopes, and this has coincided with the apparent acceleration of landslide processes. Recent landslides on man-modified slopes were caused, but not necessarily triggered, by heavy precipitation (antecedent moisture was a more critical factor than the amount of storm rainfall). Because no important landslides on natural slopes in the same period were reported in the Caramanico area, a clear distinction must be made between natural settings and those modified by man when determining rainfall thresholds for predictive purposes. In recently urbanized mountainous environments, the thresholds used to assess landslide hazards should not be weighted too heavily on old historical records of precipitation and associated mass movements. Instead, more weight ought to be given to the period following the occurrence of any major anthropogenic and natural (e.g. high-magnitude earthquake) modification of slope setting. Received: 19 October 1996 · Accepted: 25 June 1997     

Validation of Landslide Susceptibility Maps; Examples and Applications from a Case Study in Northern Spain

A procedure for validating landslide susceptibility maps wasapplied in a study area in northern Spain and the results obtained compared. Validationwas used to carry out sensitivity analysis for individual variables and combinationsof variables. The validity of different map-making methods was tested, as well as theutility of different types of Favourability Functions. The results obtained show thatvalidation is essential to determine the predictive value of susceptibility maps. Italso helps to better select the most suitable function and significant variables, thus improving the efficiency of the mapping process. Validation based on a temporal strategy makes it possible to derive hazard maps from susceptibility maps.     

Landslide risk assessment using concepts of danger pixels and fuzzy set theory in Darjeeling Himalayas

Landslide risk assessment (LRA) is a key component of landslide studies. The landslide risk can be defined as the potential for adverse consequences or loss to human population and property due to the occurrence of landslides. The LRA can be regional or site-specific in nature and is an important information for planning various developmental activities in the area. LRA is considered as a function of landslide potential (LP) and resource damage potential (RDP). The LP and RDP are typically characterized by the landslide susceptibility zonation map and the resource map (i.e., land use land cover map) of the area, respectively. Development of approaches for LRA has always been a challenge. In the present study, two approaches for LRA, one based on the concept of danger pixels and the other based on fuzzy set theory, have been developed and implemented to generate LRA maps of Darjeeling Himalayas, India. The LRA map based on the first approach indicates that 1,015 pixels of habitation and 921 pixels of road section are under risk due to landslides. The LRA map derived from fuzzy set theory based approach shows that a part of habitat area (2,496 pixels) is under very high risk due to landslides. Also, another part of habitat area and a portion of road network (7,204 pixels) are under high risk due to landslides. Thus, LRA map based on the concept of danger pixels gives the pixels under different resource categories at risk due to landslides whereas the LRA map based on the concept of fuzzy set theory further refines this result by defining the degree of severity of risk to these categories by putting these into high and low risk zones. Hence, the landslide risk assessment study carried out using two approaches in this paper can be considered in cohesion for assessing the risks due to landslides in a region.     

Regional hydrological thresholds for landslides and floods in the Tiber River Basin (central Italy)

 The definition of landslide warning thresholds, based on the analysis of hydrological data, is proposed. In the Tiber River Basin of central Italy historical information on landslides and floods, for the period 1918–1990, was available from a nationwide bibliographical and archive inventory on geohydrological catastrophes. Hydrological data were obtained from mean daily discharge records at various gauging stations within the basin. Several hundred hydrological events, broadly defined as a series of consecutive days having mean daily discharge exceeding a predefined value, were identified. Hydrological parameters obtained from the discharge records were used to rank the events according to their probability to trigger mass movements or inundations and to define regional thresholds for the occurrence of landslides and floods. The proposed approach, not lacking limitations, has conceptual and operational advantages, among which is the possibility of using historical information on geohydrological catastrophes. Received: 20 November 1996 · Accepted: 25 June 1997     

Landslide susceptibility mapping: A comparison of logistic regression and neural networks methods in a medium scale study, Hendek region (Turkey)

Landslide susceptibility mapping is one of the most critical issues in Turkey. At present, geotechnical models appear to be useful only in areas of limited extent, because it is difficult to collect geotechnical data with appropriate resolution over larger regions. In addition, many of the physical variables that are necessary for running these models are not usually available, and their acquisition is often very costly. Conversely, statistical approaches are currently pursued to assess landslide hazard over large regions. However, these approaches cannot effectively model complicated landslide hazard problems, since there is a non-linear relationship between nature-based problems and their triggering factors. Most of the statistical methods are distribution-based and cannot handle multisource data that are commonly collected from nature. In this respect, logistic regression and neural networks provide the potential to overcome drawbacks and to satisfy more rigorous landslide susceptibility mapping requirements. In the Hendek region of Turkey, a segment of natural gas pipeline was damaged due to landslide. Re-routing of the pipeline is planned but it requires preparation of landslide susceptibility map. For this purpose, logistic regression analysis and neural networks are applied to prepare landslide susceptibility map of the problematic segment of the pipeline. At the end, comparative analysis is conducted on the strengths and weaknesses of both techniques. Based on the higher percentages of landslide bodies predicted in very high and high landslide susceptibility zones, and compatibility between field observations and the important factors obtained in the analyses, the result found by neural network is more realistic.     

Field variability of landslide model parameters

 A data set of parameters (slope, soil depth and soil shear strength) relevant to spatially distributed modelling of shallow landslides triggered by rain and snowmelt events was determined from field measurements in 250 grid elements of dimensions 25 m (downslope)×10 m (across slope) in an area of 250 m×250 m on a hillslope in Scotland. These data provide an unusually detailed basis for the evaluation of spatial variability and uncertainty in model parameterisation. The variations in slope and soil strength are represented adequately by normal distributions; a Weibull distribution is suggested for the soil depth data. The factor of safety calculated at each point in the grid was shown partially to identify observed landslides, with a number of false predictions of occurrence. Trend analysis and semivariogram analysis of the data set suggest that the use of kriging could improve upon this approach to landslide prediction by providing areal estimates of parameters at the grid element scale with associated error bounds. Received: 30 October 1996 · Accepted: 25 June 1997     

Hydrogeological and geotechnical assessments of the Kadifekale landslide area, İzmir, Turkey

 The study area is located in the south of İzmir city centre in an area built without planning permission or any overall city plan. In this area, a number of mass movements occurred in the past and the region is still an active landslide field at present. The real sliding factor is geological structure, which is made up of hard volcanic rocks overlying soft Neogene clayey soils, and forms a typical structure prone to sliding. Unplanned human activities change the hydrogeological and geotechnical stabilities of the geological formations in a negative way. Rain and water leaking either directly from the mains or septic holes infiltrate into the ground and act as one of the factors in causing landslides. It is clearly shown that the irregular urbanization in and around Kadifekale is one of the factors contributing to the landslides in the area. To prevent the occurrence of landslides in the study area an effective surface and under ground drainage should be established. Rain and wastewater should be removed from the area by separate systems. Slopes should be reduced, water-loving trees should be planted and construction of high rise buildings should be avoided. Received: 12 December 1997 · Accepted: 4 April 2000     

Evolution of the Academician Ridge Accommodation Zone in the central part of the Baikal Rift, from high-resolution reflection seismic profiling and geological field investigations

 New high-resolution seismic reflection data from the central part of Lake Baikal provide new insight into the structure and stratigraphy of Academician Ridge, a large intra-rift accommodation zone separating the Central and North Baikal basins. Four seismic packages are distinguished above the basement: a thin top-of-basement unit; seismic-stratigraphic unit X; seismic-stratigraphic unit A; and seismic-stratigraphic unit B. Units A and B were cored on selected key locations. The four packages are correlated with a series of deposits exposed on the nearby western shores: the Ularyar Sequence (Oligocene); the Tagay Sequence (Lower to Middle Miocene); the Sasa Sequence (Upper Miocene to Lower Pliocene); the Kharantsy Sequence (Upper Pliocene); and the Nyurga Sequence (Lower Pleistocene). Based on stratal relationships, sedimentary geometries, distribution patterns and principal morphostructural elements – both onshore and offshore – we propose a new palaeogeographic evolution model for the area. In this model progressive tectonic subsidence of the Baikal basins and successive pulses of uplift of various segments of the rift margins lead to: (a) formation of the ridge as a structural and morphological feature separating the Central and North Baikal basins during the Middle to Late Miocene; (b) gradual flooding of the main parts of the ridge and establishment of a lacustrine connection between the two rift basins during the Late Miocene; and (c) total submergence of the top parts of the crest of the ridge during the latest Pleistocene. This new model helps to better constrain numerous phases in the structural evolution of the Baikal Rift, in which the Academician Ridge as an accommodation zone plays a crucial role. Received: 26 November 1999 / Accepted: 12 March 2000     

Heavy metals in freshly deposited sediments of the Gomati River (a tributary of the Ganga River): effects of human activities

 The concentrations of various metals (Cr, Cu, Co, Fe, Mn, Ni, Pb, Zn, and Cd) were determined in recently deposited surface sediments of the Gomati River in the Lucknow urban area. Markedly elevated concentrations (milligrams per kilogram) of some of the metals, Cd (0.26–3.62), Cu (33–147), Ni (45–86), Pb (25–77), and Zn (90–389) were observed. Profiles of these metals across the Lucknow urban stretch show a progressive downstream increase due to additions from 4 major drainage networks discharging the urban effluents into the river. The degree of metal contamination is compared with the local background and global standards. The geoaccumulation index order for the river sediments is Cd>Zn>Cu>Cr>Pb. Significant correlations were observed between Cr and Zn, Cr and Cu, Cu and Zn and total sediment carbon with Cr and Zn. This study reveals that the urbanization process is associated with higher concentrations of heavy metals such as Cd, Cu, Cr, Pb, and Zn in the Gomati River sediments. To keep the river clean for the future, it is strongly recommended that urban effluents should not be overlooked before their discharge into the river. Received: 16 February 1996 · Accepted: 29 February 1996     

Landslide complex in the northern part of the Silesian Beskydy Mountains (Czech Republic)

The study brings new information concerning a large slope deformation in the northern part of the Silesian Beskydy Mountains in the territory of the Czech Republic. The studied area is a large-scale transformation of the forefront of a flysch nappe formed especially by rigid sandstones upon a plastic underlying nappe. As a consequence of the vertical contact of mechanically different formations and depositing conditions of rocks and tectonics on the nappe forefront, several generations of deep-seated slope deformations occurred accompanied by shallow landslides in the historic period. Structural measurement and geophysical sounding proved prevailing lateral spreading in the upper parts of the slopes as well as the occurrence of rotational landslides whose landslide scarps overlap with the line of overthrust of the Godula Nappe on the Těšín Nappe. Radiocarbon dating determined the period of the occurrence of rotational landslides as falling into the Subboreal up to Subatlantic.     

Landslide susceptibility mapping using GIS and digital photogrammetric techniques: a case study from Ardesen (NE-Turkey)

Ardesen is a settlement area which has been significantly damaged by frequent landslides which are caused by severe rainfalls and result in many casualties. In this study a landslide susceptibility map of Ardesen was prepared using the Analytical Hierarchy Process (AHP) with the help of Geographical Information Systems (GIS) and Digital Photogrametry Techniques (DPT). A landslide inventory, lithology–weathering, slope, aspect, land cover, shear strength, distance to the river, stream density and distance to the road thematics data layers were used to create the map. These layer maps are produced using field, laboratory and office studies, and by the use of GIS and DPT. The landslide inventory map is also required to determine the relationship between these maps and landslides using DPT. In the study field in the Hemsindere Formation there are units that have different weathering classes, and this significantly affects the shear strength of the soil. In this study, shear strength values are calculated in great detail with field and laboratory studies and an additional layer is evaluated with the help of the stability studies used to produce the landslide susceptibility map. Finally, an overlay analysis is carried out by evaluating the layers obtained according to their weight, and the landslide susceptibility map is produced. The study area was classified into five classes of relative landslide susceptibility, namely, very low, low, moderate, high, and very high. Based on this analysis, the area and percentage distribution of landslide susceptibility degrees were calculated and it was found that 28% of the region is under the threat of landslides. Furthermore, the landslide susceptibility map and the landslide inventory map were compared to determine whether the models produced are compatible with the real situation resulting in compatibility rate of 84%. The total numbers of dwellings in the study area were determined one by one using aerial photos and it was found that 30% of the houses, with a total occupancy of approximately 2,300 people, have a high or very high risk of being affected by landslides.     

Environmental isotopes in determining groundwater flow systems, northern part of Epirus, Greece

On the basis of the isotopic composition of water in the northern part of Epirus, Greece, from springs at different altitudes with well-defined recharge areas, the altitude effect on the δ¹⁸O value of groundwater is –0.142±0.003ö (100?m)^–1 and is uniform over the entire study area. Using the δ¹⁸O composition of surface water and groundwaters, the contribution of Ioannina Lake and the channel draining the lake water to the Kalamas River to the recharge of springs and boreholes was confirmed and quantitatively defined. In contrast, the Voidomatis and Vikos Rivers are not sources for recharge of the big springs along their banks. However, water from the Aoos River does replenish the aquifer in the unconsolidated deposits underlying the plain of Konitsa. In addition, limestones of Senonian–Late Eocene ages, dolomites, and limestones of the "Vigles" facies are hydraulically interconnected, and the limestones of the "Pantokrator" facies are hydraulically isolated from the other carbonate formations.     

A comparative study of conventional, ANN black box, fuzzy and combined neural and fuzzy weighting procedures for landslide susceptibility zonation in Darjeeling Himalayas

Landslides are one of the most destructive phenomena of nature that cause damage to both property and life every year, and therefore, landslide susceptibility zonation (LSZ) is necessary for planning future developmental activities. In this paper, apart from conventional weighting system, objective weight assignment procedures based on techniques such as artificial neural network (ANN), fuzzy set theory and combined neural and fuzzy set theory have been assessed for preparation of LSZ maps in a part of the Darjeeling Himalayas. Relevant thematic layers pertaining to the causative factors have been generated using remote sensing data, field surveys and Geographic Information System (GIS) tools. In conventional weighting system, weights and ratings to the causative factors and their categories are assigned based on the experience and knowledge of experts about the subject and the study area to prepare the LSZ map (designated here as Map I). In the context of objective weight assignments, initially the ANN as the black box approach has been used to directly produce an LSZ map (Map II). In this approach, however, the weights assigned are hidden to the analyst. Next, the fuzzy set theory has then been implemented to determine the membership values for each category of the thematic layer using the cosine amplitude method (similarity method). These memberships are used as ratings for each category of the thematic layer. Assuming weights of each thematic layer as one (or constant), these ratings of the categories are used for the generation of another LSZ map (Map III). Subsequently, a novel weight assignment procedure based on ANN is implemented to assign the weights to each thematic layer objectively. Finally, weights of each thematic layer are combined with fuzzy set derived ratings to produce another LSZ map (Map IV). The maps I–IV have been evaluated statistically based on field data of existing landslides. Amongst all the procedures, the LSZ map based on combined neural and fuzzy weighting (i.e., Map IV) has been found to be significantly better than others, as in this case only 2.3% of the total area is found to be categorized as very high susceptibility zone and contains 30.1% of the existing landslide area.     

Weathered bedrock as a source of sand and gravel aggregate in north-eastern New South Wales, Australia

 The most commonly recognized sources of sand and gravel aggregate are the deposits of coastal, fluvial and glacial sedimentary processes. It is not commonly recognized that weathered bedrock is also an important source of sand and gravel aggregate. In the case of weathered sedimentary bedrock the product may be indistinguishable from modern sedimentary materials and deposits may be misidentified. Batson's quarry is a major sand and gravel aggregate resource for the rapidly developing region of north-eastern New South Wales. The deposit has previously been assumed to be a young sedimentary deposit, but detailed field examination and comparison with surrounding rock types indicates that it has formed from weathered Mesozoic sandstone bedrock. The extraction of weathered bedrock aggregate deposits has distinctive environmental implications because, unlike sedimentary aggregates, they are not restricted by the geometry of sedimentary environments and do not necessarily interfere with surface drainage systems. This study shows that the recognition of the geological origin of a resource is essential in order to determine its potential geographical extent, and thereby maximize its utilization and minimize land use conflict. Received: 3 January 1996 / Accepted: 1 November 1996     

Rigid Unit Modes in disordered nepheline: a study of a displacive incommensurate phase transition

 Calculations of the Rigid Unit Modes (RUMs) allowed in the nepheline structure are used to explain the diffuse scattering previously seen in electron diffraction experiments. The RUM calculations also show that the modulation wavelength for incommensurate nephelines is essentially determined by the framework topology. X-ray diffraction is used to measure the intensity of the diffuse scattering as a function of temperature. The diffuse intensity increases sharply at 308 K. This effect is interpreted as being due to the softening of a phonon mode, indicating a phase transition. Measurements of this phase transition below the transition temperature are made using hard mode infrared spectroscopy. Received: 17 February 1999 / Revised, accepted. 15 October 1999     

Landslide susceptibility mapping of the slopes in the residual soils of the Mengen region (Turkey) by deterministic stability analyses and image processing techniques

The aim of the present study is to prepare a landslide susceptibility map of a region of about 120 km², between Gökcesu and Pazarköy (around Mengen, NW Turkey) at approximately 10 km north of the North Anatolian Fault Zone, where frequent landslides occur. For this purpose, mechanisms of the landslides were studied by two-dimensional stability analyses together with field observations, and the parameters controlling the development of such slides were identified. Field observations indicated that the failures generally developed within the unconsolidated and/or semiconsolidated soil units in forms of rotational, successive shallow landslides within the weathered zone in Mengen, Cukurca and Sazlar formations. Although consisting of residual soils, Capak and Gökdag formations do not exhibit landslides as the natural slopes formed on these, do not exceed the critical slope angles. Statistical evaluations and distribution of the landslides on the topographical map showed that such parameters as cohesion, angle of internal friction, slope, relative height, orientation of slopes, proximity to drainage pattern, vegetation cover and proximity to major faults were the common features on the landslides. Digital images were obtained to represent all these parameters on gray scale on the SPOT image and on the digital elevation model (DEM) of the area using image processing techniques. Soil mechanics tests were carried out on 36 representative samples collected from different units, and parameters were determined for two-dimensional stability analyses basing on “sensitivity approach” and for the preparation of digital shear strength map. In order to determine the critical slope angle values for the residual soils, a series of sensitivity analyses were realized by using two-dimensional deterministic slope stability analyses techniques for varying values of cohesion, angle of internal friction and slope height along with varying saturation conditions. According to the results of the sensitivity analyses, the Mengen formation was found to be most susceptible unit to landslides, covering about 33.5% of the region studied in terms of surface area. The distribution of the critical slopes were determined by superimposing the critical slope values from sensitivity analyses on slope map of the study area. On the other hand, iso-cohesion and iso-friction maps were produced by locating the values of cohesion and internal friction angles in a geographic coordinate system such that they coincide with sample locations on the DEM and by further interpolation of the values concerned. The pixel values were evaluated in gray scale from 0 to 255, 0 representing the lowest pixel value and 255 representing the highest. Sensitivity analyses on cohesion and angle of internal friction to investigate the effects of these parameters only on stability, revealed that cohesion was effective at a rate of 70% by itself, while angle of internal friction alone controlled the stability by a rate of 30%. The iso-cohesion and iso-friction maps previously obtained were digitally combined in these rates and a “shear strength map” was prepared. The geographic setting of the study area is such that northern slopes usually receive dense precipitation. In relation to this fact, about 42% of the landslides are due north. Thus, a slope orientation map was prepared using the DEM, and slopes facing north were evaluated as being more susceptible to sliding. Proximity to the drainage pattern was another important factor in the evaluation, as streams could adversely affect the stability by either eroding the toe or saturating the slope, or both. When considered together, in conjunction with the field observations, faults and landslides showed a close association. In the area, about 88% of the landslides were detected within an area closer than 250 m to major faults, therefore, a main discontinuity map was produced using the SPOT image of the region, and “proximity to major faults” was evaluated as a parameter as most of the landslides developed in areas where the vegetation was rather sparse. A vegetation cover map was therefore obtained from the SPOT image, and the areas with denser vegetation were considered to be less susceptible to sliding with respect to the areas with less or no vegetation. Having prepared the maps accounting for the distribution of critical slopes, shear strength properties, relative height, slope angle, orientation of the slopes, vegetation cover, proximity to the drainage pattern, geographic corrections were carried on each of these, and a potential failure map was obtained for the residual soils by superimposing all these maps. Next, a classification was performed on the final map and five relative zones of susceptibility were defined. When compared with this map, all of the landslides identified in the field were found to be located in the most susceptible zone. The performance of the method used in processing the images appears to be quite high, the zones determined on the map being the zones of relative susceptibility.