Next Article in Journal
Photogrammetry in the Study of the Antique and Medieval Archaeological Site of Markul (Village Markula, Ochamchira Region, Republic of Abkhazia) in Northwestern Colchis Area (Black Sea Coast of the Caucasus)
Next Article in Special Issue
Debris Flows Occurrence in the Semiarid Central Andes under Climate Change Scenario
Previous Article in Journal
Tectono-Sedimentary Cenozoic Evolution of the El Habt and Ouezzane Tectonic Units (External Rif, Morocco)
Previous Article in Special Issue
Linking Soil Hydrology and Creep: A Northern Andes Case
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Impact of Quality of Digital Elevation Models on the Result of Landslide Susceptibility Modeling Using the Method of Weights of Evidence

by
Mirosław Kamiński
Polish Geological Institute–National Research Institute, Rakowiecka 4, 02-519 Warszawa, Poland
Geosciences 2020, 10(12), 488; https://doi.org/10.3390/geosciences10120488
Submission received: 29 September 2020 / Revised: 25 November 2020 / Accepted: 30 November 2020 / Published: 3 December 2020
(This article belongs to the Special Issue Landslides and Granular Flows on Earth)

Abstract

:
The paper discusses the impact that the quality of the digital elevation model (DEM) has on the final result of landslide susceptibility modeling (LSM). The landslide map was developed on the basis of the analysis of archival geological maps and the Light Detection and Ranging (LiDAR) digital elevation model. In addition, complementary field studies were conducted. In total, 92 landslides were inventoried and their degree of activity was assessed. An inventory of the landslides was prepared using a 1-m-LiDAR DEM and field research. Two digital photogrammetric elevation models with an elevation pixel resolution of 20 m were used for landslide susceptibility modeling. The first digital elevation model was obtained from a LiDAR point cloud (DEM–airborne laser scanning (ALS)), while the second model was developed based on archival digital stereo-pair aerial images (DEM–Land Parcel Identification System (LPIS)). Both models were subjected to filtration using a Gaussian low-pass filter to reduce errors in their elevation relief. Then, using ArcGIS software, a differential model was generated to illustrate the differences in morphology between the models. The maximum differences in topographic elevations between the DEM–ALS and DEM–LPIS models were calculated. The Weights-of-Evidence model is a geostatistical method used for the landslide susceptibility modeling. Six passive factors were employed in the process of susceptibility generation: elevation, slope gradient, exposure, topographic roughness index (TRI), distance from tectonic lines, and distance from streams. As a result, two landslide susceptibility maps (LSM) were obtained. The accuracy of the landslide susceptibility models was assessed based on the Receiver Operating Characteristic (ROC) curve index. The area under curve (AUC) values obtained from the ROC curve indicate that the accuracy of classification for the LSM–DEM–ALS model was 78%, and for the LSM–LPIS–DEM model was 73%.

1. Introduction

Landslides are a significant geodynamic threat in many areas of the world. They cause major economic losses as well as threaten human life [1]. The Flysch Carpathians area in Poland, especially the Dynów Foothills, is predisposed to the formation of various types of mass movements, mainly landslides [2,3,4]. The purpose of this research was to develop an optimal model of landslide susceptibility, taking into account digital photogrammetric elevation models of different qualities, which were used for the modeling process. For this purpose, two digital elevation models were used. The first was developed based on filtration and classification of the cloud of elevation points originating from aerial laser scanning (ALS), which is hereafter referred to as DEM–ALS [5]. The other model was developed based on digital stereo-pairs of aerial images, which is hereafter referred to as DEM–LPIS (Land Parcel Identification System) [6,7,8,9,10]. The general availability of DEM–ALS and airborne imagery derived from DEM in geodetic centers in Poland was the reason for using these for the comparison process. A 20 m resolution was chosen for this study because global research into the impact of the resolution on the landslide susceptibility model shows this to be optimal [11]. In addition, the quality of the airborne imagery was low, which made it difficult to obtain greater detail.
When looking through the literature, it can be seen that many researchers have focused on the impact of the digital elevation model resolution on the final correctness of landslide susceptibility classifications [12,13]. It was found that DEM resolution is an important criterion and that it can significantly affect the accuracy of landslide susceptibility maps. It was noticed that the higher the DEM spatial resolution, the more strongly the results depended on topographic details [14,15,16]. It is worth noting that the digital elevation model (DEM) is also a key part of the basic spatial data used to analyze landslide susceptibility. The digital elevation model is the basis for generating factors determining the formation of landslides. These include, among others, the slope angle, the spatial orientation of slopes, the elevation, the curvature, and the topographic roughness index (TRI) [17]. DEMs are often generated using data obtained from different remote sensors, including optical imaging sensors, light detection and ranging (LiDAR), and synthetic aperture radar (SAR) [18]. The qualities of DEM-derived factors often depend on the spatial resolution of DEMs. This has been widely discussed in the literature [19]. Therefore, the choice of DEM is important for the assessment of landslide susceptibility.
The quality of the DEMs is essential for assessing their suitability and determines the quality of the geomorphometric analysis [20,21,22]. As a result of the variety of available DEMs, it is necessary to investigate their use in different scenarios. It is necessary to remember that working with digital data requires paying particular attention to their quality. Small errors in DEMs can produce large errors in the derived elevation attributes [23], especially second-order derivatives such as curvature [24]. DEM accuracy depends on the type of topography and ruggedness of the elevation, as well as the type of vegetation [25], the methods for collecting elevation data, the method for DEM generation, the type of DEM grid, and DEM resolution [26,27]. The issue of error analysis in DEMs is relevant and discussed in the literature [28,29]. This study focuses on the usefulness of different DEMs with the same resolution for landslide susceptibility studies. These analyses consisted of checking their vertical accuracy.
The impact that the quality of digital elevation models based on various sources of data has on landslide susceptibility models has not yet been widely discussed in the professional literature. In my opinion, landslide susceptibility modeling is more influenced by the quality of data (especially in digital elevation models) than the number of factors used.
The Weights-of-Evidence (WoE) method is a geostatistical method used for landslide susceptibility modeling [30,31,32,33,34,35,36,37,38]. This is a logarithmic form of the Bayesian probability model, often used in geomorphological risk assessment [33,39,40,41,42,43]. Two-dimensional statistical analysis methods, such as WoE, are amongst many of those used to assess landslide susceptibility. GIS-based geostatistical methods have become very popular in evaluating landslide susceptibility [44,45,46,47,48,49,50,51,52,53,54,55,56]. The Receiver Operating Characteristic (ROC) curve was used to assess the quality of landslide susceptibility classification maps [57,58]. Therefore, this research hopes to contribute to the scientific discussion on the problem presented in the article.

2. Study Area

The study area, the Dynów Foothills, is located in southeastern Poland in the San River valley (Figure 1). The topography of the area follows its geological structure. The elevation relief features wide ridges rising to elevations exceeding 400 m above sea level. The slopes of the ridges are cut by stream valleys. The upstream areas are usually trough-like and pass downstream into V-shaped gully valleys. Their lower reaches have narrow and flat bottoms. The slopes of the ridges are usually convex or straight and gently inclined at an angle of up to 20°, except in the headwaters areas, where they are trough-shaped. The relief is enriched with numerous landslides. In terms of geological structure, the Dynow Foothills are built of overthrust and folded flysch formations of the Skole Nappe [59].
Thick covers of Quaternary deposits, which favor the formation of shallow landslides, occur in the study area. These are mainly tills of various origins, as well as loesses. Their thicknesses are variable and can be up to 10 m. In addition, the flysch rock lithology accompanied by tectonic features promotes the formation of structural landslides. In the discussed area, the youngest flysch rocks are Oligocene shales and sandstones of the Krosno Beds. They consist of carbonate-cemented sandstones (Figure 1). They are interbedded with grey marly shales. These layers are underlain by Eocene variegated shales represented by clayey red, green, grey and spotted rocks, i.e., the so-called variegated shales [59]. They are predisposed to the formation of landslides due to their physical properties (cleavage, plasticity, impermeability). The variegated shales are underlain by the Cretaceous sandstone–shale flysch. The sandstones are usually thin-bedded and interlayered with ashy grey mudshales [59]. As a result of its physical properties, the flysch is highly susceptible to landsliding. In terms of tectonic structure, the research area is located within the overthrust zone of the Jawornik–Dubiecko overturned fold onto the Chyrzynka–Huta syncline. This zone makes this region of the Dynow Foothills more susceptible to the formation of landslides by weakening the structure of the rocks and allowing rainwater migration within them.

3. Methods and Materials

3.1. Archival Digital Stereo-Pair Aerial Images

In order to precisely determine the extent of a landslide and to present its surface relief, an accurate photogrammetric digital elevation model was used. It was purchased at the National Surveying and Cartographic Documentation Centre. Ten digital stereo-pairs of black-and-white aerial images, at a scale of 1:10,000, taken in 2009 as part of the Land Parcel Identification System (LPIS) project, were used to develop the model. The digital elevation model was generated on the Dephos photogrammetric station using an automatically generated 20 × 20 m grid of elevation points. The grid was subjected to stereoscopic editing consisting in determining the values of the grid points (the so-called pickets). The measured values of the height of the grid pickets were enriched with structural relief lines, such as scarp edges, lines of ridges and watercourses, and areas excluded due to dense vegetation. Using the Inverse Distance Weight (IDW) algorithm, a digital elevation model was generated in a 20 × 20 m grid structure. It is very common to use an algorithm to interpolate the spatial data [60]. The average fitting error for the photogrammetric model in the grid is presented in Table 1.

3.2. Airborne Laser Scanning (ALS)—LiDAR Data

The elevation relief was examined based on elevation data from the IT system of the country’s protection against extraordinary threats (ISOK, http://www.gugik.gov.pl/projekty/isok). The products that LiDAR (Light Detection and Ranging) obtained from that project were made available by the National Surveying and Cartographic Documentation Centre in Warsaw. The data acquired during laser scanning are saved as files in the ASPRS Accuracy Standards for Digital Geospatial Data’s Laser File Format (LAS). The scanning density was four points per square meter. The basic procedures for processing the point cloud obtained as a result of airborne laser scanning (ALS) are their classification and filtration [61]. The elevation data were acquired in September 2014. The process of generating DEM was based on the point cloud classification. The point cloud was divided into the following layers: ground, low vegetation, medium vegetation, high vegetation, and buildings. Then, the points representing the ground were filtered from the point cloud, using the tools available in the QCoherent software. The next step was the interpolation of elevation points using the algorithms. Two algorithms were tested: Inverse Distance Weighting (IDW) and the Adaptive Triangulated Irregular Network (TIN) model algorithm. Both algorithms provide a similar classification for point clouds describing land use for agriculture, e.g., areas on which a single building, shrub, or tree is located. The Adaptive TIN algorithm works better than the IDW algorithm in terms of the points recorded by the laser beam being reflected from ground, vehicles, and bridges. The Adaptive TIN Ground Point Cloud Task is an automated algorithm designed to separate points that have a high probability of being ground points from other points. The algorithm divides the task area into cells whose X and Y dimensions are defined by a seed parameter. For each cell, a “best” candidate ground point is selected. These “seed” points are then used to construct a Triangulated Irregular Network (TIN). The algorithm then iterates, attempting to add additional points to the TIN, based on several inclusion criteria. The iterations conclude when either the user-specified stopping criteria is met or no additional points were added during the previous iteration. As a result, two digital elevation models were developed in the grid structure at mesh sizes of 1 mm and 20 m. The first was used to identify landslides, the second for landslide susceptibility modeling.

3.3. Gaussian Low-Pass Filter

On the basis of the literature [62], the Gaussian low-pass filter was selected for further generalization of the surface of both digital elevation models. Low-pass filters are the averaging filters with certain weights. By adjusting the filter matrix row correctly, you can remove disturbances and errors of varying magnitudes from a digital elevation model of grid structure. To generalize the surface of the digital elevation models, four Gaussian low-pass filters were tested: STd_dev_0.391.ker; STd_dev_0.625.ker; STd_dev_1.0.ker; and STd_dev_1.6.ker, which are available in the ER Mapper software. The best results of the digital elevation model generalization were obtained after applying the STd_dev_1.6.ker 9 filter [63]. The obtained results of elevation generalization were compared in the ArcGis software with the hypsometry of the topographic map. It was assumed that the topographic map at the scale of 1:10,000 is a sufficiently accurate reflection of the topography+. The operation of the Gaussian low-pass filter relies on removing details from the image (e.g., differences between adjacent pixels), thereby smoothing the image. This filter was used to generalize the relief in both the DEM–ALS and DEM–LPIS digital elevation models.

3.4. Weights-of-Evidence Method

The WoE method is one of the multi-criteria decision-making methods. It is a log-linear version of Bayesian General Theory. This method has been used in many fields of science, e.g., in medical diagnostics and geology (determining potential sites of mineral deposits) [64,65,66,67], as well as for determining slopes susceptible to landsliding [68]. In probability theory and statistics, Bayes’ theorem describes the probability of an event, based on prior knowledge of conditions that might be related to the event.
This is expressed by the formula
P ( L / F ) = P ( L F ) P ( F ) = P ( L ) P ( F / L ) P ( F )
where
P (L) is an a priori probability of a given phenomenon (landslide occurrence probability);
F is a factor determining the occurrence of a given phenomenon;
P (F/L) is a conditional probability;
P (L/F) is a probability that a given phenomenon will occur due to factor F (a posteriori).
The a priori probability can be defined as the ratio of the area occupied by the landslide to the entire study area.
For practical reasons, the odds expression is introduced, which is defined as the ratio of the probability of occurrence of a given phenomenon to the probability of its not occurring.
This is expressed by the formula
O d d s = O = P ( L ) 1 P ( L ) .
As a result of further transformations, conditional odds are determined, and then the so-called association indicators (weights of “premises”) W+ and W−, i.e., Weights-of-Evidence (which gave the name to the method), are calculated based on their logarithmic quantities. The difference between W + and W− is referred to as contrast (C). This is one of the elements, which determines the acceptance or rejection of a class of variables adopted for modeling.
W + = ln P { B / D } P { B / D }
W = ln P { B / D } P { B / D }
where
P is a landslide occurrence probability;
B is a factor determining the occurrence of a given phenomenon;
D is a factor indicating the existence or absence of landslide.
The W+ positive weight shows more evidence than chance alone and is, therefore, directly proportional to the influence that the respective variable has on the likelihood of landslide occurrence. Contrast (C) is the difference between two weights (C = W+ − W−). This is one of the elements that determines the acceptance or rejection of a class of variables adopted for modeling.

4. Results

4.1. Landslide Inventory

The landslide inventory was made using archival digital aerial images and a DEM–ALS [69]. Supplementary field studies were also carried out [70,71]. Ninety-two landslides with a total area of 225,182 m2 were surveyed and their activity model developed (Figure 2). In order to develop a model of landslide activity, three field trips were made. Active slopes of landslides were located using a Real Time Kinematic (RTK GPS) receiver. The new fractures observed both on the landslide surface and on buildings and roads were the main criteria for landslide activity. Many landslides had new fractures in the main and minor scarps. This mainly indicated the activity of landslides.
For comparison purposes, the landslide ratio/index and the landslide density index were calculated [2].
O p = P o P t ×   100 %
where
Op—Landslide ratio/index [%];
PO—Landslide area [km2];
Pt—Surface area [km2].
G = n P t
where
G—Landslide density ratio/index [number/km2];
n—Number of landslides [number];
Pt—Surface area [km2].
The landslide ratio/index was 24.7%, while the landslide density ratio/index was 1.316/km2 [2]. For comparison, the density ratio/index given in the literature for the Dynow Foothills area was 0.193/km2 [2]. The reason for such a large difference was the greater number of inventory landslides.
The largest landslide clusters occur in the western and northeastern parts of the study area. In the western part, the largest landslide, number 36, registered in the western part of the study area, covers an area of over 2.5 ha (The actual size of this landslide is much greater because only its northeastern part extends into the area under study. It is developed totally within the Inoceramian Beds of Cretaceous flysch.
Most of the landslides occur in the northwestern part of the study area. These are obsequent landslides developed mainly within a dense network of stream valleys overgrown with dense forests (Figure 3). Many of these landslides threaten local roads and buildings (Figure 4). There may be a phenomenon of stream capture between landslides number 10 and number 11 because their main slopes are only about 3 m apart. Landslides numbers 1 and 3 slide downslope to a stream whose flow direction in the channel is strongly disturbed by the advancing landslide toes. As a result of their increased activity, the stream might be blocked by the toes of both landslides and a dam lake might form. In relation to the strata dip, landslide number 3 is a consequent landslide. Its length is over 240 m and it covers an area of 2.7 ha (Figure 5). The main scarp is about 20 m high. It developed within flysch deposits of the Krosno Beds. The colluvial material consists of tills with rock debris. It threatens two residential buildings located near the landslide. Over the last 10 years, the landslide has not shown significant activity. Landslide number 1 is located on the southwestern slope. It is a typical rotational landslide, i.e., the sliding rock masses have rotated. As a result of the movement direction of rock material in relation to the position of bedrock layers, it is an obsequent landslide. Its length is over 350 m and its area is about 3.5 ha. The main scarp exceeds 1.5 m in height. The landslide developed between the Eocene variegated shales and the Oligocene flysch of the Krosno Beds. It was activated in spring 2001. It activated as a result of the springtime melting of wet snow on the slope for a long period. As a result, it damaged a local road and destroyed one outbuilding. A number of other residential buildings found themselves in a zone of immediate high risk of damage. Most of the landslides in this part of the study area are not large in terms of surface area. However, they are mostly active and show clearly developed main scarps and numerous, recently formed fractures.
For these reasons, they pose a great threat to local infrastructure. Landslides located in the northeastern part of the research area are larger than those found in its western part. The largest of them is landslide number 62, which occupies an area of over 2 ha. It threatens local roads and numerous residential and utility buildings. It is a subsequent landslide that formed on a boundary between four lithological complexes: Cretaceous shale–sandstone flysch of the Inoceramian Beds, Eocene variegated shales and shale flysch of the Hieroglyphic Beds, and shale–sandstone flysch of the Oligocene Krosno Beds. The other landslides are slightly smaller; in most cases, they are classified as obsequent and subsequent landslides with distinct main scarps. Some of these landslides were formed in forests, e.g., number 86 and number 82, while others developed within arable fields and meadows, which is why they are a big problem for agriculture and forestry.

4.2. Landslide Conditioning Factors and Their Analysis

Landslide susceptibility modeling comprised four main stages. The first stage consisted of collecting data and constructing the GIS database. These were the layers with landslides and nine predictive factors: lithology, faults, streams, elevation, slope gradient, spatial orientation of slopes, curvature, roughness index (TRI), and terrain position index (TPI) (Figure 6). On the basis of field research and an analysis of archival materials, a model of landslide activity was developed. One hundred and ninety-three active landslide scarps were identified. Then, the dataset was divided randomly into two subsets. The first subset of data (97 slopes) was used to develop susceptibility models, while the second subset of data (96) was used to verify these models.
The modeling was carried out using the ArcSDM (Spatial Data Modeler) module, which is an extension of ArcGIS, developed by the Geological Survey of Canada [72]. This module automatically calculates test chi-square positive and negative weight values, variances, contrast, and posteriori probability (Table 2 and Table 3). The first three factors were created in a vector form in the geobase, while the remaining six were developed as a result of DEM–ALS and DEM–LPIS processing. All layers were converted to a raster format in a grid structure, with a 20 m pixel size. In the second stage, a chi-squared test of independence was performed for the nine predictive factors. The test chi-square is based on a comparison of the observed values (obtained in the study) and theoretical values (calculated assuming that there is no relationship between the variables). Large differences indicate the existence of dependencies between the variables [32].
Layers such as lithology, curvature, and topographic position index (TPI) showed a conditional relationship. Therefore, in a later stage of modeling, these were eliminated. Ultimately, only six predictive factors were used: faults, streams, elevation, slope gradient, spatial orientation of slopes, and roughness index (TRI). The third stage covered the use of the Weights-of-Evidence method to weigh the individual predictors. The calculation of weights for predictive quantities includes crossing operations between raster maps and a point map of the occurrence of active landslide scarps. As a result, two maps of landslide susceptibility were generated. The fourth stage involved estimating the accuracy and validation of the obtained landslide susceptibility models.

4.3. Comparison of the Elevation Differences between the Digital Elevation Models (ALS and LPIS)

A comparison of the quality of the photogrammetric digital elevation models was made in ArcGis.10.6 software. Two digital elevation models, LPIS and ALS, with an elevation pixel resolution of 20 m, were used for the analysis. Both models were subjected to a low-pass Gaussian filtration process. Using the Arc Map tools, a differential elevation model between the DEM–ALS and DEM–LPIS models was generated (Figure 7). In the western part of the differential model, two areas (number 1 and number 2) were distinguished, where there are the largest discrepancies between the models. This research showed a spatial distribution of elevation changes between the DEM–ALS and DEM–LPIS models. Their maximum values ranged from −21 m to +15 m. Such considerable differences in elevations between the models are mainly caused by the masking effect of vegetation on the surface (Figure 8).
Area number 1 is marked on the orthophotomap. It is covered with dense vegetation. Area number 2 is more diverse in terms of vegetation cover. Meadows and arable land are distinguished here, in addition to high and medium vegetation. In other areas, covered mostly with forests, a significant increase in differences in elevations between both models was noted. In addition, a morphological cross-section was constructed for area number 1 (Figure 9). In the upper part of the cross-section, running through the meadows, slight differences in elevation of about 2 m were found. In turn, in the middle and lower parts of the cross-section, crossing the valleys of streams overgrown with trees, there was a significant increase in elevation differences. The maximum differences in the cross-section ranged from −8 m to + 7 m.

4.4. Landslide Susceptibility Maps

The susceptibility maps were developed by means of the Weights-of-Evidence geostatistical method using two different 20-m-resolution digital elevation models. The model allowed us to divide the study area into five morphodynamic categories: very weakly susceptible, weakly susceptible, moderately susceptible, highly susceptible, and very highly susceptible (Figure 9).
Stable areas, i.e., very poorly susceptible and poorly susceptible areas, are represented mainly by the stream valleys consisting of Holocene and Quaternary terraces. The same category refers to the flattening of hills, composed largely of Cretaceous sandstone–shale rocks and Oligocene shale–sandstone rocks. The remaining area was included in the category “unstable” (moderately, highly, and very highly susceptible). The most susceptible areas occur mainly on the hill slopes and predominantly follow the course of geological strata, mainly of the variegated shales and the Cretaceous shale–sandstone flysch of the Inoceramian Beds. There are several key differences between the LSM–ALS and LSM–LPIS landslide susceptibility models. These differences relate to the distribution and size of the areas occupied by individual susceptibility classes. The greatest number of the most weakly susceptible areas occurs in the LSM–LPIS model (Figure 10). Low and medium susceptibility areas occupy larger areas in the LSM–ALS model than in the LSM–LPIS model, whereas highly susceptible areas cover larger areas in the LSM–LPIS model than in the LSM–ALS model. Very highly susceptible areas occupy similar surface areas in both models (Figure 11).
The modeling results showed that the greatest number of landslides occurs in the classes “high” (12.9%) and “very high” (29.03%) in the LSM–ALS model (Figure 11).

5. Validation

The characteristics of the Receiver Operating Characteristic (ROC) curve were used to validate the two landslide susceptibility models obtained using the WoE geostatistical method. The ROC approach shows values between the true positive rate (TPR) and false positive rate (FPR) in determining the correctness of landsliding susceptibility. The y-axis indicates TPR, while the x-axis indicates FPR. The TPR is the pixels to which landslide areas were correctly assigned, while the FPR is the pixels to which landslide areas were incorrectly assigned. The Area Under Curve (AUC) is a measure indicating the accuracy of the landslide susceptibility models. The AUC values obtained from the ROC curve indicate that the classification accuracy was 78% for the LSM–DEM–ALS model and 73% for the LSM–LPIS–DEM model. The resulting ROC curve is shown in Figure 12.

6. Discussion

It is clear that the DEM–ALS and DEM–LPIS digital photogrammetric elevation models differ in terms of the accuracy of relief details. On the basis of the research, it should be stated that the DEM–ALS model is more accurate than the DEM–LPIS model. There are two main reasons for the quality deterioration of the DEM–LPIS model. The first is related to the methods of generating digital elevation models based on the digitization of stereo-pairs of aerial images. This process is carried out by the operator on a digital photogrammetric station. It is always a possibility that the operator will incorrectly measure both the elevation point values of the measuring grid and the structural lines of the relief, e.g., scarp edges and stream lines. The second reason for errors in the model is related to the dense vegetation covering the area, e.g., forests and shrubs, and compact settlements. The differential model comparing the DEM–ALS and DEM–LPIS models showed that the greatest differences occurred in areas covered with vegetation, which ranged from +15 m to −21 m and are considered large errors. The author is of the opinion that the elevation models generated digitally from aerial images are significantly burdened by the effect of land cover; therefore, they should be used to a limited extent in spatial analyses [73]. It should also be added that the DEM–ALS model is not without flaws. There are at least two factors that affect the accuracy of surface imaging. The first is the correct filtration of the points representing the elevation surface. In most cases, this is done automatically. However, in this way, points of low vegetation, for example, may be classified incorrectly as the surface [74]. The second factor is the selection of a method of elevation point interpolation. This is a very important stage in generating digital elevation models, which requires a separate discussion. The present study used the Adaptive TIN algorithm, which is found in QCoherent software. The research presented in this paper refers to a wide research stream focused on the impact that the resolution and quality of digital elevation models has on the result of landslide susceptibility [75,76,77,78].
This study also demonstrated the usefulness of the Weights-of-Evidence geostatistical method in evaluating landslide susceptibility in mountainous areas [79,80]. The main advantage of this method is the ability to integrate many predictive factors and check, by means of the chi-square test, which are the most reliable for susceptibility modeling. The susceptibility maps differed in the distribution and size of areas occupied by individual landslide susceptibility classes. These differences were affected by the different digital elevation models used. The LSM–LPIS susceptibility model is characterized, for instance, by an uneven distribution of pixels representing high and very high susceptibility. This is an obvious error, which results from the poor quality of the DEM–LPIS digital model. Regarding the LSM–LPIS model, very poorly susceptible areas also covered over 60% of the study area. This percentage value seems to be too high as compared with 36% from the LSM–ALS model.

7. Conclusions

The Weights-of-Evidence method used in the study proved to be useful in landslide susceptibility modeling. Two digital elevation models, each with a 20 m pixel resolution, were employed for the modeling. The first model was developed on the basis of the classification, filtration, and interpolation of the elevation point cloud acquired from DEM–ALS aerial laser scanning, while the other was developed on a digital photogrammetric station based on digital stereo-pairs of aerial images at the scale of 1:10,000. Both models were filtered in the raster format with a Gaussian low-pass filter. The filter generalized the relief and reduced the impact of various artifacts on the accuracy of the models. A comparative analysis between the digital elevation models showed that the DEM–ALS model is more accurate than the DEM–LPIS model. The precision of photogrammetric measurements is mainly dependent on the quality of the aerial photographs. The photographic resolution is a function of the optical quality of an image and is influenced by the resolving power of the film and camera lens, image motion during exposure, atmospheric conditions, and the film processing conditions. The errors that occurred in the DEM–LPIS model were mainly due to the masking effect of vegetation, which effectively limits the ability to see details in the area examined. The second factor affecting the accuracy of the model was the subjective assessment of the elevations of the measuring grid points and the structural lines appearing in the aerial images, e.g., slope ridge lines and stream lines, by the operator of the photogrammetric station. It is worth noting that the density and type of vegetation may be a source of errors in both digital elevation models. On the other hand, errors in the DEM–ALS model may be generated both in the process of point cloud classification and its interpolation. Using the AUC curve, the classification accuracy was assessed to be 78% for the LSM–DEM–ALS model and 73% for the LSM–LPIS–DEM model. Roads, residential buildings, and forests fell within a range of very high- and high-risk classes. This research helps us to better understand the process by which various errors form and their impact on the final result of landslide susceptibility modeling. The presented methodology for developing landslide susceptibility maps has practical applications in terms of spatial planning for local self-government authorities and for banks, developers, and insurance companies. In order to improve the methodology, we plan to test other geostatistical and machine-learning methods in terms of their impact on the accuracy of landslide susceptibility evaluation. In addition, the effect of different resolutions of LiDAR data and methods of interpolation and filtration will be examined.

Funding

The final editorial work was supported by the PGI-NRI statutory funds (Project No. 62.9012.1942.00.0).

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Aleotti, P.; Chowdhury, R. Landslide hazard assessment: Summary review and new perspectives. Bull. Eng. Geol. Environ. 1999, 58, 21–44. [Google Scholar] [CrossRef]
  2. Rączkowski, W.; Mrozek, T. Activating of landsliding in the Polish Flysch Carpathians by the end of the 20th century. Studia Geomorphol. Carpatho Balc. 2002, 36, 91–111. [Google Scholar]
  3. Bober, L. Landslide areas in the Polish part of the Carpathians flysh and their relationship with geological structure. Biult. Inst. Geolog. 1984, 340, 115–162. [Google Scholar]
  4. Kamiński, M. Landslide GIS analysis for selected area of the Dynów foreland. Arch. Photogramm. Cartogr. Remote Sens. 2006, 16, 279–287. [Google Scholar]
  5. Derron, M.H.; Jaboyedoff, M. LIDAR and DEM technique for landslides monitoring and characterization. Nat. Hazard Earth Syst. Sci. 2010, 10, 1877–1879. [Google Scholar] [CrossRef]
  6. Greve, C. Digital Photogrammetry an Addendum to the Manual of Photogrammetry. In American Society for Photogrammetry and Remote Sensing; Asprs Publications: Bethesda, MD, USA, 1996; pp. 121–123. [Google Scholar]
  7. Genelli, G.; Giusti, E.; Pizzaferri, G. Photogrammetric Technique for the Investigation of the Corniglio Landslide. In Applied Geomorphology, Theory and Practice; John Wiley & Sons: Hoboken, NJ, USA, 2002; pp. 23–28. [Google Scholar]
  8. Mora, P.; Baldi, P.; Casula, G.; Fabris, M.; Ghirotti, M.; Mazzioni, E.; Pesci, A. Global positioning systems and digital photogrammetry for monitoring of mass movements: Application to the Ca’ di Malta landslide (northern Apennines, Italy). Eng. Geol. 2003, 68, 103–121. [Google Scholar] [CrossRef]
  9. Dewitte, O.; Jasselette, J.-C.; Cornet, Y.; Eeckhaut, M.V.; Collignon, A.; Poesen, J.; Demoulin, A. Tracking landslide displacements by multi-temporal DTMs: A combined aerial sterophotogrammetric and LiDAR approach in western Belgium. Eng. Geol. 2010, 99, 11–22. [Google Scholar] [CrossRef]
  10. Kamiński, M. Application od photogrammetric methods to assess the dynamics of mass movements selected exmaples from Poland. Arch. Photogramm. Cartogr. Remote Sens. 2012, 24, 11–122. [Google Scholar]
  11. Meen, R.S.; Nachappa, G.T. Impact of Spatial Resolution of Digital Elevation Model on Landslide Susceptibility Mapping: A case Study in Kullu Valley, Himalayas. Geosciences 2019, 9, 67–74. [Google Scholar]
  12. Mohd, W.M.; Abdullah, M.A.; Hashim, S. Evaluation of Vertical Accuracy of Digital Elevation Models Generated from Different Sources: Case Study of Ampang and Hulu Langat, Malaysia. FIG Congress XXV.1-17. Geosciencies 2014, 9, 2–18. [Google Scholar]
  13. Miner, A.S.; Flentje, P.; Mazengarb, C.; Windle, D.J. Landslide Recognition using LiDAR derived Digital Elevation Models-Lessons learnt from selected Australian examples. In Proceedings of the Geologically Active Proceedings 11th IAEG Congregalia, Auckland, New Zealand, 5–10 September 2010; Volume 2, p. 352. [Google Scholar]
  14. Wang, B.; Shi, W.; Liu, E. Robust methods for assessing the accuracy of linear interpolated DEM. Int. J. Appl. Earth Obs. Geoinf. 2015, 34, 198–206. [Google Scholar] [CrossRef]
  15. Wu, B.; Tang, S.; Zhu, Q.; Tong, K.; Hu, H.; Li, G. Geometric integration of high-resolution satellite imagery and airborne LiDAR data for improved geopositioning accuracy in metropolitan areas. ISPRS J. Photogramm. Remote Sens. 2015, 109, 139–151. [Google Scholar] [CrossRef]
  16. Chang, K.T.; Doub, J.; Changc, Y.; Kuoa, C.P.; Xua, K.M.; Liud, J.K. Spatial resolution effects of digital terrain models on landslide susceptibility analysis. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, XLI-B8, 45–49. [Google Scholar]
  17. Komac, M. A landslide susceptibility model using the analytical hierarchy process method and multivariate statistics in perialpine Slovenia. Geomorphology 2006, 74, 17–28. [Google Scholar] [CrossRef]
  18. Wang, Q.; Guo, Y.; Li, W.; He, J.; Wu, Z. Predictive modeling of landslide hazards in Wen County, northwestern China based on information value, weights-of-evidence, and certainty factor. Geomat. Nat. Hazards Risk 2019, 10, 820–835. [Google Scholar] [CrossRef] [Green Version]
  19. Schlögel, R.; Marchesini, I.; Alvioli, M.; Reichenbach, P.; Rossi, M.; Malet, J.P. Optimizing landslide susceptibility zonation: Effects of DEM spatial resolution and slope unit delineation on logistic regression models. Geomorphology 2018, 301, 10–20. [Google Scholar] [CrossRef]
  20. Desmet, P.J.J. Effects of interpolation errors on the analysis of DEMs. Earth Surf. Process. Landf. 1997, 22, 563–580. [Google Scholar] [CrossRef]
  21. Florinsky, I.V. Accuracy of local topographic variables derived from digital elevation models. Int. J. Geogr. Inf. Sci. 1998, 12, 47–61. [Google Scholar] [CrossRef]
  22. Chaplot, V.; Darboux, F.; Bourennane, H.; Leguédois, S.; Silvera, N.; Phachomphon, K. Accuracy of interpolation techniques for the derivation of digital elevation models in relation to landform types and data density. Geomorphology 2006, 77, 126–141. [Google Scholar] [CrossRef]
  23. Fisher, P.F.; Tate, N.J. Causes and consequences of error in digital elevation models. Prog. Phys. Geogr. 2006, 30, 467–489. [Google Scholar] [CrossRef]
  24. Wise, S.M. The effect of GIS interpolation errors on the use of DEMs in geomorphology. In Landform Monitoring, Modeling and Analysis; Lane, S.N., Richards, K.S., Chandler, J.H., Eds.; Wiley: Chichester, UK, 1998; pp. 139–164. [Google Scholar]
  25. Thomas, J.; Prasannakumar, V.; Vineetha, P. Suitability of space-borne digital elevation models of different scales in topographic analysis: An example from Kerala, India. Environ. Earth Sci. 2015, 73, 1245–1263. [Google Scholar] [CrossRef]
  26. Zhang, W.; Montgomery, D. Digital elevation model grid size, landscape representation, and hydrologic simulations. Water Resour. Res. 1994, 30, 1019–1028. [Google Scholar] [CrossRef]
  27. Li, Z.L. A comparative study of the accuracy of digital terrain models (DTMs) based on various data models. ISPRS J. Photogramm. Remote Sens. 1994, 49, 2–11. [Google Scholar] [CrossRef]
  28. Bolstad, P.V.; Stowe, T. An evaluation of DEM accuracy: Elevation, slope, and aspect. Photogramm. Eng. Remote Sens. 1994, 60, 1327–1332. [Google Scholar]
  29. Edson, C.; Wing, M.G. LiDAR Elevation and DEM Errors in Forested Settings. Mod. Appl. Sci. 2015, 9, 139–157. [Google Scholar] [CrossRef] [Green Version]
  30. Bonham-Carter, G.F.; Agtberg, F.P.; Wright, D.F. Weight of evidence modeling: A new approach to mapping mineral potential. Statistical Applications in the Earth Sciences. In Geological Survey of Canada; Agterberg, F.P., Bonham-Carter, G.F., Eds.; Paper; F.A. Acland: Ottawa, ON, Canada, 1989; Volume 89, pp. 171–183. [Google Scholar]
  31. Bonham-Carter, G. Geographic Information System for Geoscientists: Modeling with GIS. In Computer Methods in the Geoscience; Pergamon Press: Oxford, UK, 1994; Volume 13, pp. 398–399. [Google Scholar]
  32. Agterberg, F.P.; Cheng, Q. Conditional independence test for weights of evidence modeling. Nat. Resour. Res. 2002, 11, 249–255. [Google Scholar] [CrossRef]
  33. Lee, S.; Choi, J.; Min, K. Landslide susceptibility analysis and verification using the Bayesian probability model. Environ. Geol. 2002, 43, 120–131. [Google Scholar] [CrossRef]
  34. Chung, J.O.F.; Fabbri, A.G. Validation of spatial prediction models for landslide hazard mapping. Nat. Hazard 2003, 30, 451–472. [Google Scholar] [CrossRef]
  35. Mrozek, T.; Poli, S.; Sterlascchini, S.; Zabuski, L. Landslide susceptibility assessment. A case study from the Beskid Niski Mts., Carpathains, Poland. In Proceedings of the Conference “Risks Caused by the Geodynamic Phenomena in Europe“, Wysowa, Poland, 20–22 May 2004; Volume 15, pp. 10–16. [Google Scholar]
  36. Neuhauser, B.; Terhorst, B. Landslide susceptibility assessment using ‘‘weights-of-evidence’’ applied to a study area at the Jurassic escarpment (SW-Germany). Geomorphology 2007, 86, 12–24. [Google Scholar] [CrossRef]
  37. Armas, I. Weights of Evidence method for landslide susceptibility mapping. Prahova Subcarpathians, Romania. Nat. Hazards 2011, 6, 31–38. [Google Scholar] [CrossRef]
  38. Kamiński, M. Application of LiDAR date to assess tha landslide susceptibility map using Weights of Evidence method—An example from Podhale region (Southern Poland). In Proceedings of the International Society for Photogrammetry and Remote Sensing ISPRS Congress, Prague, Czech Republic, 12–19 July 2016; Volume XLI-B1, pp. 25–29. [Google Scholar]
  39. Dai, F.C.; Lee, C.F.; Ngai, Y.Y. Landslide risk assessment and management: An overview. Eng. Geol. 2002, 64, 65–87. [Google Scholar] [CrossRef]
  40. Lee, S.; Choi, J.; Min, K. Probabilistic landslide hazard mapping using GIS and remote sensing data at Boun, Korea. Int. J. Remote Sens. 2004, 25, 2037–2052. [Google Scholar] [CrossRef]
  41. Nguyen, Q.P.; Bui, H.B. Landslide hazard mapping using Bayesian approach in GIS: Case study in YangSan area, Korea. In Proceedings of the International Symposium on GeoInformatics for Spatial-Infrastructure Development in Earth and Allied Sciences (GIS-IDEAS), Hanoi, Vietnam, 16–18 September 2004; Available online: http://gisws.media.osaka-cu.ac.jp/gisideas04/viewpaper.php?id=2 (accessed on 2 December 2020).
  42. Kamiński, M. Landslide susceptibility map in regionak scale an example from of San Valley in the Dynów Foothill. Biul. Państw. Inst. Geol. 2012, 452, 109–118. [Google Scholar]
  43. Ilia, I.; Tsangaratos, P. Applying weight of evidence method and sensitivity analysis to produce a landslide susceptibility map. Landslides 2016, 13, 379–397. [Google Scholar] [CrossRef]
  44. Varnes, D.J. Landslide hazard zonation: A review of principles and practice of the United Nations Educational, Scientific and Cultural Organization. Nat. Hazards 1984, 3, 127–138. [Google Scholar]
  45. Carrara, A.; Cardinali, M.; Detti, R.; Guzzetti, F.; Pasqui, V.; Reichenbach, P. GIS techniques and statistical models in evaluating landslide hazard. Earth Surf. Proc. Land 1991, 16, 427–445. [Google Scholar] [CrossRef]
  46. Van Westen, C.J. Application of Geographic Information Systems to Landslide Hazard Zonation. Ph.D. Thesis, Technical University Delft, Enschede, The Netherlands, 1993. [Google Scholar]
  47. Van Westen, C.J. Statistical Landslide Hazard Analysis, ILWIS 2.1 for Windows Application Guide; ITC Publication: Enschede, The Netherlands, 1997; pp. 73–84. [Google Scholar]
  48. Remondo, J.; Gonzales, A.; De Teran, J.R.D.; Fabrii, A.; Chung, C.J. Validation of landslide susceptibility maps; examples and applications from a case study in northern Spain. Nat. Hazards 2003, 1, 1–13. [Google Scholar] [CrossRef]
  49. Chung, C.J.; Fabbri, A.G. Probabilistic prediction models for landslide hazard mapping. Photogramm. Eng. Remote Sens. 1999, 65, 1389–1399. [Google Scholar]
  50. Glade, T. Vulenrability assessmnet in landslide risk analysis. Die Erde 2003, 134, 123–146. [Google Scholar]
  51. Glade, T.; Anderson, M.G.; Crozier, M.J. (Eds.) Landslide Hazard and Risk; Wiley: Chichester, UK, 2005; pp. 125–134. [Google Scholar]
  52. Van Westen, C.; Rengers, N.; Soeters, R. Use of geomorphological information in indirect landslide susceptibility assessment. Nat. Hazards 2003, 30, 399–419. [Google Scholar] [CrossRef]
  53. Van Westen, C.J.; Van Asch, T.W.J.; Soeters, R. Landslide hazard and risk zonation: Why is it still so difficult. Bull Eng. Geol. Environ. 2006, 65, 167–184. [Google Scholar] [CrossRef]
  54. Guzzetti, F. Landslide Hazard and Risk Assessment. Ph.D. Thesis, Rheinische Friedrich-Wilhelm-Universitaet Bonn, Bonn, Germany, 2006. Available online: http://hss.ulb.uni-bonn.de/dissonline/mathnatfak/2006/guzzettifausto/index.htm (accessed on 15 June 2009).
  55. Kamiński, M. Landslide susceptibility map: A case study from the Jodłówka region (Dynów Foothill)). Geol. Rev. 2007, 55, 779–784. [Google Scholar]
  56. Kouli, M.; Loupasakis, C.; Soupios, P.; Vallianatos, F. Landslide hazard zonation in high risk areas of Rethymno Prefecture, Crete Island, Greece. Nat. Hazards 2010, 52, 599–621. [Google Scholar] [CrossRef]
  57. Swets, J.A. ROC analysis applied to the evaluation of medical imaging techniques. Invest. Radiol. 1979, 14, 109–112. [Google Scholar] [CrossRef] [PubMed]
  58. Hanley, J.A. Receiver operating characteristic (ROC) methodology: The state of the art. Crit Rev. Diagn. Imaging 1989, 243. [Google Scholar]
  59. Kamiński, M.; Piotrowska, K. Detailed Geolgical Map of Poland at the Scale of 1:50,000 Kańczuga Sheet with Explanations. 2014. Available online: http:/bazadata.pgi.gov.pl/data/smgp/arkusze_skany/smgp 1006.jpg (accessed on 14 April 2014).
  60. Bartier, M.P.; Keller, P.C. Multivariate interpolation to incorporate thematic surface data using inverse distanece weighting. Comput. Geosci. 1996, 22, 795–799. [Google Scholar] [CrossRef]
  61. Jaboyedoff, M.; Oppikofer, T.; Abellán, A.; Derron, M.H.; Loye, A.; Metzger, R.; Pedrazzini, A. Use of LIDAR in landslide investigations: A review. Nat. Hazards 2012, 61, 5–28. [Google Scholar] [CrossRef] [Green Version]
  62. Milledge, G.D.; Stuart, L.; Warburton, J. Digital filtering of generic topographic data in geomorphological research. Earth Surf. Process. Landf. 2009, 34, 63–74. [Google Scholar] [CrossRef]
  63. Good, I.J. Statistical Evidence. In Encyclopedia of Statistics; Kotz, S., Johnson, N.L., Eds.; Wiley: New York, NY, USA, 1988; Volume 8, pp. 651–655. [Google Scholar]
  64. Bonham-Carter, G.F.; Agtberg, 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]
  65. Agterberg, F.P.; Bonham-Carter, G.F.; Wright, D.F. Statistical Pattern Integration for Mineral Exploration. In Computer Applications in Resource Estimation Prediction and Assessment for Metals and Petroleum; Gaal, G., Merriam, D.F., Eds.; Pergamon Press: Oxford, UK, 1990. [Google Scholar]
  66. Weed, D.L. Weight of evidence: A review of concept and methods. Risk Anal. 2005, 25, 1545–1557. [Google Scholar] [CrossRef]
  67. Fan, D.; Cui, X.; Yuan, D.; Wang, J.; Yang, J.; Wang, S. Weight of Evidence Method and Its Applications and Development. Procedia Environ. Sci. 2011, 11, 1412–1418. [Google Scholar] [CrossRef] [Green Version]
  68. Dahal, R.K.; Hasegawa, S.; Nonomura, A.; Yamanaka, M.; Masuda, T.; Nishino, K. GIS-Based Weights-of-Evidence Modelling of Rainfall-Induced Landslides in Small Catchments for Landslide Susceptibility Mapping. Environ. Geol. 2008, 54, 311–324. [Google Scholar] [CrossRef]
  69. Schulz, W.H. Landslides Mapped Using LIDAR Imagery, Seattle, Washington. US Geol. Surv. Open File Rep. 2004, 11, 1–13. [Google Scholar]
  70. Ardizzone, F.; Cardinali, M.; Galli, M.; Guzzetti, F.; Reichenbach, P. Identification and mapping of recent rainfall-induced landslides using elevation data collected by airborne LiDAR. Nat. Hazards Earth Syst. Sci. 2007, 7, 637–650. [Google Scholar] [CrossRef] [Green Version]
  71. Bell, R.; Petschko, H.; Röhrs, M.; Dix, A. Assessment of landslide age, landslide persistence and human impact using airborne laser scanning digital terrain models. Geogr. Ann. 2012, 94, 135–156. [Google Scholar] [CrossRef]
  72. Kemp, L.D.; Bonham-Carter, G.F.; Raines, G.L. Arc-WofE: ArcView Extension for Weights of Evidence Mapping. Ottawa, ON, Canada, 1999. Available online: http://ntserv.gis.nrcan.gc.ca/sdm/ (accessed on 25 February 2002).
  73. Fabeis, M.; Pesci, A. Automated extraction in digital aerial photogrammetry: Precisions and validation for mass movement monitoring. Ann. Geophys. 2005, 48, 973–988. [Google Scholar]
  74. Stephen, E.; Reutebuch, R.; McGaughey, J.; Hans-Erik, A.; Ward, W.C. Accuracy of a high-resolution lidar terrain model under a conifer forest canopy. Can. J. Remote Sens. 2003, 29, 527–535. [Google Scholar] [CrossRef]
  75. Cama, M.; Conoscenti, C.; Lombardo, L.; Rotigliano, E. Exploring relationships between grid cell size and accuracy for debris-flow susceptibility models: A test in the Giampilieri catchment Sicily, Italy. Environ. Earth Sci. 2016, 753, 238–254. [Google Scholar] [CrossRef]
  76. Chang, K.T.; Merghali, A.; Yunus, A.P.; Pham, B.T.; Dou, J. Evaluating scale effects of topographic variables in landslide susceptibility models using GIS-based machine learning techniques. Sci. Rep. 2019, 9, 142–157. [Google Scholar] [CrossRef] [Green Version]
  77. Arnone, E.; Francipane, A.; Scarbaci, A.; Puglisi, C.; Noto, L.V. Effect of raster resolution and polygon-conversion algorithm on landslide susceptibility mapping. Environ. Model Softw. 2016, 84, 467–481. [Google Scholar] [CrossRef]
  78. Achour, Y.; Pourghasemi, H.R. How do machine learning techniques help in increasing accuracy of landslide susceptibility maps? Geosci. Front. 2019. [Google Scholar] [CrossRef]
  79. Barbieri, G.; Cambuli, P. The weight of evidence statistical method in landslide susceptibility mapping of the Rio Pardu Valley (Sardinia, Italy). In Proceedings of the 18th World IMACS/MODSIM Congress, Cairns, Australia, 13–17 July 2009; pp. 13–17. [Google Scholar]
  80. Neuhauser, B.; Damm, B.; Terhorst, B. GIS-based assessment of landslide susceptibility on the base of the Weights-of Evidence model. Landslides 2012, 9, 511–528. [Google Scholar] [CrossRef]
Figure 1. Study area location and geological setting.
Figure 1. Study area location and geological setting.
Geosciences 10 00488 g001
Figure 2. Landslide distribution map and active scarp.
Figure 2. Landslide distribution map and active scarp.
Geosciences 10 00488 g002
Figure 3. Landslides in a Light Detection and Ranging (LiDAR) image before and after filtration.
Figure 3. Landslides in a Light Detection and Ranging (LiDAR) image before and after filtration.
Geosciences 10 00488 g003
Figure 4. The locations of buildings and roads against the background of landslide settlement.
Figure 4. The locations of buildings and roads against the background of landslide settlement.
Geosciences 10 00488 g004
Figure 5. 3D visualization of landslide number 1. Buildings that are endangered by the landslide are visible. Visualizations were developed in ENVI LiDAR software using point cloud filtration products.
Figure 5. 3D visualization of landslide number 1. Buildings that are endangered by the landslide are visible. Visualizations were developed in ENVI LiDAR software using point cloud filtration products.
Geosciences 10 00488 g005
Figure 6. Diagram of a landslide susceptibility map development.
Figure 6. Diagram of a landslide susceptibility map development.
Geosciences 10 00488 g006
Figure 7. (A) Orthophotomap. (B) Differential elevation model comparing the DEM–ALS and DEM –LPIS models. The two areas highlighted are the areas with the largest differences in elevation between the DEM–ALS and DEM–LPIS models.
Figure 7. (A) Orthophotomap. (B) Differential elevation model comparing the DEM–ALS and DEM –LPIS models. The two areas highlighted are the areas with the largest differences in elevation between the DEM–ALS and DEM–LPIS models.
Geosciences 10 00488 g007
Figure 8. Morphological cross-section.
Figure 8. Morphological cross-section.
Geosciences 10 00488 g008
Figure 9. Landslide susceptibility maps against the road network. (A) LSM–ALS susceptibility, (B) LSM–LPIS susceptibility.
Figure 9. Landslide susceptibility maps against the road network. (A) LSM–ALS susceptibility, (B) LSM–LPIS susceptibility.
Geosciences 10 00488 g009
Figure 10. Spatial visualization of maps of landslide susceptibility (A,B). (C) Percentage list comparing landslide susceptibility for the LSM–ALS and LSM–LPIS models.
Figure 10. Spatial visualization of maps of landslide susceptibility (A,B). (C) Percentage list comparing landslide susceptibility for the LSM–ALS and LSM–LPIS models.
Geosciences 10 00488 g010
Figure 11. Histogram representing the distribution of observed landslides falling into various susceptibility classes of different landslide susceptibility zonation maps for the LSM–ALS and LSM–LPIS models.
Figure 11. Histogram representing the distribution of observed landslides falling into various susceptibility classes of different landslide susceptibility zonation maps for the LSM–ALS and LSM–LPIS models.
Geosciences 10 00488 g011
Figure 12. Validation of susceptibility maps.
Figure 12. Validation of susceptibility maps.
Geosciences 10 00488 g012
Table 1. Summary of the results of aerial triangulation.
Table 1. Summary of the results of aerial triangulation.
RMS x (m)RMS y (m)RMS Total (m)The Number of Aerial ImagesYear
0.440.360.32122009
Table 2. Spatial relationship between each landslide conditioning factor and landslide using the LSM–ALS model.
Table 2. Spatial relationship between each landslide conditioning factor and landslide using the LSM–ALS model.
ClassArea (km2)PointsW+s (W+)W−s (W−)Cs (C)Stud (C)
Distance to fault (m)
0–5062662−0.21080.70980.00490.1025−0.21570.7171−0.3008
50–10067601−1.17671.00140.02340.1020−1.20011.0066−1.1922
100–150718720.20910.7111−0.00390.10250.21300.71850.2964
150–20074180
>200757830.61110.5823−0.01430.10310.62540.59131.0577
Elevation (m.a.l.s)
230–23511841−2.39631.00040.10860.1021−2.50481.0056−2.4909
235–25523687−1.14140.37850.17900.1054−1.32040.3929−3.3605
255–312355212−1.00740.28920.27780.1085−1.28530.3088−4.1615
312–363473628−0.44530.18950.25410.1202−0.69930.2245−3.1155
363–415592046−0.17010.14800.17920.1394−0.34920.2034−1.7176
Aspect
N1154140.25370.2689−0.03670.10960.29040.29041.0001
N–E2308280.25390.1901−0.08590.12000.33980.22491.5111
S–E3461420.25400.1552−0.15550.13420.40940.20521.9953
S–W4615500.13930.1422−0.12710.14490.26640.20301.3121
N–W5768650.17900.1247−0.27980.17470.45880.21472.1374
Slope (degrees)
0–28220
2–521266−1.21400.40880.16770.1048−1.38160.4221−3.2735
5–13376818−0.68560.23630.25050.1125−0.93610.2617−3.5773
13–20549246−0.12050.14810.12030.1394−0.24080.2034−1.1841
20–547107680.01380.1219−0.03050.18340.04430.22020.2012
Topographic Ruggedness Index
(−160)–(−2)11523−1.29600.57810.08760.1031−1.38360.5872−2.3561
(−2)–(−1)23055−1.47880.44770.20100.1043−1.67980.4597−3.6543
(−1)–0345714−0.85270.26780.25400.1098−1.10680.2894−3.8240
0–1461028−0.44520.18960.25400.1203−0.69930.2245−3.1150
1–12576250−0.08590.14200.09810.1451−0.18410.2030−0.9066
Distance to streams (m)
0–556841−1.84691.00070.05660.1020−1.90361.0059−1.8924
55–69912100.17800.3180−0.01840.10710.19630.33550.5851
69–9035650.42930.4504−0.01860.10420.44790.46230.9690
90–12049170.44300.3807−0.02720.10530.47010.39501.1902
Table 3. Spatial relationship between each landslide conditioning factor and landslide using the LSM–LPIS model.
Table 3. Spatial relationship between each landslide conditioning factor and landslide using the LSM–LPIS model.
ClassArea (km2)PointsW+s (W+)W−s (W−)Cs (C)stud (C)
Distance to fault (m)
0–50626661−0.00230.12870.00930.2595−0.01160.2896−0.0399
50–1006760660.00060.1237−0.00390.31780.00450.34100.0131
100–150718770−0.00190.12010.02210.4103−0.02400.4275−0.0561
150–200741872−0.00540.11840.10250.5027−0.10790.5165−0.2089
>2007578740.00080.1168−0.02830.71050.02910.72000.0404
Elevation (m.a.l.s)
230–2359031−2.17751.00060.08410.1021−2.26161.0057−2.2487
235–25522679−0.89780.33400.15910.1066−1.05680.3506−3.0144
255–312378618−0.71680.23630.26890.1125−0.98570.2617−3.7667
312–363533341−0.23310.15680.20930.1332−0.44250.2058−2.1504
363–415684563−0.05170.12660.10030.1699−0.15200.2119−0.7173
Aspect
N8377−0.16700.37960.01420.1059−0.18120.3941−0.4599
N–E2236290.27660.1869−0.09810.12180.37470.22311.6794
S–E3412430.24740.1535−0.16090.13670.40830.20551.9869
S–W4578530.16140.1382−0.16510.15140.32650.20501.5931
N–W6019660.10660.1238−0.19500.18030.30160.21871.3787
Slope (degrees)
0–28070
2–520627−1.07350.37860.16220.1060−1.23570.3932−3.1429
5–13368118−0.70700.23630.26690.1132−0.97390.2620−3.7170
13–20531144−0.17650.15140.17470.1382−0.35120.2050−1.7134
20–546804700.04220.1201−0.10180.19330.14400.22760.6325
Topographic Ruggedness Index (TRI)
(−160)–(−2)7822−1.28600.87800.05760.2031−1.68370.5630−2.1321
(−2)–(−1)207810−0.72460.31700.13040.1078−0.85500.3348−2.5536
(−1)–0373420−0.61680.22420.25050.1147−0.86730.2518−3.4437
0–1533244−0.18160.15140.18080.1382−0.36240.2050−1.7682
1–126778670.00050.1228−0.00120.18350.00180.22080.0080
Distance to streams (m)
0–556841−1.89481.00070.06000.1021−1.95481.0059−1.9433
55–6915828−0.64940.35450.08480.1060−0.73410.3700−1.9844
69–901931120.44240.28960.08080.1084−0.52330.3092−1.6924
90–12010168980.06000.10210.89480.00071.95481.0059−1.9433
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kamiński, M. The Impact of Quality of Digital Elevation Models on the Result of Landslide Susceptibility Modeling Using the Method of Weights of Evidence. Geosciences 2020, 10, 488. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10120488

AMA Style

Kamiński M. The Impact of Quality of Digital Elevation Models on the Result of Landslide Susceptibility Modeling Using the Method of Weights of Evidence. Geosciences. 2020; 10(12):488. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10120488

Chicago/Turabian Style

Kamiński, Mirosław. 2020. "The Impact of Quality of Digital Elevation Models on the Result of Landslide Susceptibility Modeling Using the Method of Weights of Evidence" Geosciences 10, no. 12: 488. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10120488

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop