A novel ensemble bivariate statistical evidential belief function with knowledge-based analytical hierarchy process and multivariate statistical logistic regression for landslide susceptibility mapping

This study compares the landslide susceptibility maps from four application models, namely, (1) the bivariate model of the Dempster–Shafer based evidential belief function (EBF); (2) integration of the EBF in the knowledge-based analytical hierarchy process (AHP) as a pairwise comparison model processed by using all available causative factors; (3) integration of the EBF in the knowledge-based AHP as a pairwise comparison model by using high nominated causative factor weights only; and (4) integrated EBF in the logistic regression (LR) as a multivariate model by using nominated causative factor weights only. These models were tested in Pohang and Gyeongju Cities (South Korea) by using the geographic information system GIS platform. In the first step, a landslide inventory map consisting of 296 landslide locations were prepared from various data sources. Then, a total of 15 landslide causative factors (slope angle, slope aspect, curvature, surface roughness, altitude, distance from drainages, stream power index, topographic wetness index, wood age, wood diameter, wood type, forest density, soil thickness, soil texture, and soil drainage) were extracted from the database and then converted into a raster. Final susceptibility maps exhibit close results from the two models. Models 1 and 3 predicted 82.3% and 80% of testing data during the analysis, respectively. Thus, Models 1 and 3 show better performance than LR. These resultant maps can be used to extend the capability of bivariate statistical based model, by finding the relationship between each single conditioning factor and landslide locations, moreover, the proposed ensemble model can be used to show the inter-relationships importance between each conditioning factors, without the need to refer to the multivariate statistic. The research outcome may provide powerful tools for natural hazard assessment and land use planning

Ionospheric TEC from the Turkish Permanent GNSS Network (TPGN) and comparison with ARMA and IRI models

Abstract: The present study investigates the ionospheric Total Electron Content (TEC) variations in the lower mid-latitude Turkish region from the Turkish permanent GNSS network (TPGN) and International GNSS Services (IGS) observations during the year 2016. The corresponding vertical TEC (VTEC) predicted by Auto Regressive Moving Average (ARMA) and International Reference Ionosphere 2016 (IRI-2016) models are evaluated to realize their effectiveness over the region. The spatial, diurnal and seasonal behavior of VTEC and the relative VTEC variations are modeled with Ordinary Least Square Estimator (OLSE). The spatial behavior of modeled result during March equinox and June solstice indicates an inverse relationship of VTEC with the longitude across the region. On the other hand, the VTEC variation during September equinox and December solstice including March equinox and June solstice are decreasing with increase in latitude. The GNSS observed and modeled diurnal variation of the VTEC show that the VTEC slowly increases with dawn, attains a broader duration of peak around 09.00 to 12.00 UT, and thereafter decreases gradually reaching minimum around 21.00 UT..

Semi-quantitative landslide risk assessment using GIS-based exposure analysis in Kuala Lumpur City

A semi-quantitative landslide-risk assessment method, which would provide a spatial estimate of future landslide risks in a densely populated area in Kuala Lumpur City, was presented in this study. This work focused on detail risk assessment by identifying the number of elements at risk. A medium-scale analysis was performed using geospatial based techniques. The estimation of rainfall threshold and the landslide hazard map used in the current work are obtained from the previous literature published by the same authors. Subsequently, the vulnerability value was generalized, and then a valid integration between elements at risk and the hazard map was conducted to determine the expected number of elements that would likely be under direct risk. Results showed that the approximate number of predicted affected elements per pixel, as a percentage of the settlement unit, is nearly 50% in residential areas, 35% in commercial buildings, 31% in industrial buildings, 31% in utility areas, and 18% in densely populated areas. Similarly, a significant percentage of predicted losses (27%) were found for the road network. The results showed the capability of the method to approximately predict the number of infrastructure elements and the population density under landslide risk in data-scarce environments

A novel ensemble decision tree-based CHi-squared Automatic Interaction Detection (CHAID) and multivariate logistic regression models in landslide susceptibility mapping

An ensemble algorithm of data mining decision tree (DT)-based CHi-squared Automatic Interaction Detection (CHAID) is widely used for prediction analysis in variety of applications. CHAID as a multivariate method has an automatic classification capacity to analyze large numbers of landslide conditioning factors. Moreover, it results two or more nodes for each independent variable, where every node contains numbers of presence or absence of landslides (dependent variable). Other DT methods such as Quick, Unbiased, Efficient Statistic Tree (QUEST) and Classification and Regression Trees (CRT) are not able to produce multi branches based tree. Thus, the main objective of this paper is to use CHAID method to perform the best classification fit for each conditioning factors, then, combined it with logistic regression (LR) to find the corresponding coefficients of best fitting function that assess the optimal terminal nodes. In the first step, a landslide inventory map with 296 landslide locations were extracted from various sources over the Pohang-Kyeong Joo catchment (South Korea). Then, the inventory was randomly split into two datasets, 70 % was used for training the models, and the remaining 30 % was used for validation purpose. Thirteen landslide conditioning factors were used for the susceptibility modeling. Then, CHAID was applied and revealed that some conditioning factors such as altitude, soil drain, soil texture and TWI, as terminal nodes and reflected the best classification fit. Then, a proposed ensemble technique was applied and the interpretations of the coefficients showed that the relationship between the decision tree branch nodes distance from drain, soil drain, and TWI, respectively, leads to better consequences assessment of landslides in the current study area. The validation results showed that both success and prediction rates, 75 and 79 %, respectively. This study proved the efficiency and reliability of ensemble DT and LR model in landslide susceptibility mapping

Estimation of rainfall threshold and its use in landslide hazard mapping of Kuala Lumpur metropolitan and surrounding areas

This paper presents rainfall-induced landslide thresholds and predicts landslide hazard in Kuala Lumpur metropolitan city and surrounding areas. Landslide events from 2000 to 2012 were collected. The long and short antecedent rainfall days were prepared for landslide and non-landslide days simultaneously. First, threshold analysis was conducted by using data obtained from rainfall stations located in highly urbanized areas of Kuala Lumpur metropolis. Six rainfall gauges were selected, and the study area was divided into six zones according to rainfall gauges: Taman Desa Station (TD-KL), Genting Klang (GK-KL), LDG Edinburgh Station (LDGE-KL), SG Raya Hulu Langat Station (SRHL-Slg), Puchong Drop Station (PD-Slg), and Bukit Antarabangsa (BTA-Slg). After the threshold analysis was conducted for different periods (10, 15, and 30 days) in each station, reliability index test was conducted to optimize the best threshold that limits the predicted events along the study period for each region. Second, the threshold analysis results were used as input in the Poisson probability model to estimate landslide temporal probability (P T). Third, the spatial probability (P S) analysis was prepared by using the evidential belief function multiplied by the P T to obtain the hazard maps for 1-, 3-, and 5-year scenarios. Finally, a validation process was conducted to test the prediction performance of the resultant hazard map for a 1- and 2-year prediction by using the landslide inventory of 2012 to early 2014, which was not included in the modeling of the hazard map. Results showed a valid correlation between the high and moderate hazardous areas for the six zones. The predicted hazard maps indicated a quantitative assessment of the prone areas and proved to be a valid disaster management tool. The produced hazard maps may play a vital role as input component in risk analysis

A novel integrated model for assessing landslide susceptibility mapping using CHAID and AHP pair-wise comparison

This article uses an integrated methodology based on a chi-squared automatic interaction detection (CHAID) model combined with analytic hierarchy process (AHP) for pair-wise comparison to assess medium-scale landslide susceptibility in a catchment in the Inje region of South Korea. An inventory of 3596 landslide locations was collected using remote sensing, and a random sample comprising 30% of these was used to validate the model. The remaining portion (70%) was processed by the nearest-neighbour index (NNI) technique and used for extracting the cluster patterns at each location. These data were used for model training purposes. Ten landslide-conditioning factors (independent variables) representing four main domains, namely (1) topology, (2) geology, (3) hydrology, and (4) land cover, were used to produce two landslide-susceptibility maps. The first landslide-susceptibility map (LSM1) was produced by overlaying the terminal nodes of the CHAID result tree. The second landslide-susceptibility map (LSM2) was produced using the overlay result of AHP pair-wise comparisons of CHAID terminal nodes. The prediction rate curve results were better with LSM2 (area under the prediction curve (AUC) = 0.80) than with LSM1 (AUC = 0.76). The results confirmed that the integrated hybrid model has superior prediction performance and reliability, and it is recommended for future use in medium-scale landslide-susceptibility mapping

Application of an evidential belief function model in landslide susceptibility mapping

The objective of this paper is to exploit the potential application of an evidential belief function model to landslide susceptibility mapping at Kuala Lumpur city and surrounding areas using geographic information system (GIS). At first, a landslide inventory map was prepared using aerial photographs, high resolution satellite images and field survey. A total 220 landslides were mapped and an inventory map was prepared. Then the landslide inventory was randomly split into a testing dataset 70% (153 landslides) and remaining 30% (67 landslides) data was used for validation purpose. Fourteen landslide conditioning factors such as slope, aspect, curvature, altitude, surface roughness, lithology, distance from faults, ndvi (normalized difference vegetation index), land cover, distance from drainage, distance from road, spi (stream power index), soil type, precipitation, were used as thematic layers in the analysis. The Dempster–Shafer theory of evidence model was applied to prepare the landslide susceptibility maps. The validation of the resultant susceptibility maps were performed using receiver operating characteristics (ROC) and area under the curve (AUC). The validation results show that the area under the curve for the evidential belief function (the belief map) model is 0.82 (82%) with prediction accuracy 0.75 (75%). The results of this study indicated that the EBF model can be effectively used in preparation of landslide susceptibility maps

Extraction of soil moisture from RADARSAT-1 and its role in the formation of the 6 December 2008 landslide at Bukit Antarabangsa, Kuala Lumpur

Active microwave has a huge potential in the estimation of soil moisture especially over large areas where the meteorological observations are seldom. The large contrast in dielectric constant between different types of soil is considered as the main factor for measuring the moisture content. This study is aimed at the extraction of soil moisture over the areas of Bukit Antarabangsa, Malaysia using active microwave remote sensing technique in order to examine the impact of moisture content dynamically on landslides occurrence, which have been a basic challenge that threaten Bukit Antarabangsa area, particularly in falling of monsoon seasons. This study addressed a specific event that took place in 6 December 2008 due to a very high level of precipitation that resulted in a raise in ground water table causing the occurrence of landslide. One Radarsat-1 image acquired in July 2008 before the landslide was used for generating the moisture content map. The resultant moisture content map showed a reasonable distribution of the moisture concentrated over the forest areas which has previous records landslides. Moreover, it was found that the previous landslide events were within the high moisture zone indicating the presence of high moisture content. Subsequently, three moisture maps were extracted from Landsat-7 ETM+, which were then used for validation process. A statistically based validation technique was used by calculating area under the curve that correlates the high moisture values of three images. In order to validate the Landsat-7 ETM+ moisture content, monthly rainfall data was plotted against the high moisture values derived from three Landsat-7 images. The validation result indicated an acceptable compatibility. The spatial relation between high moisture areas in Landsat-7 ETM+ images along the year resulted in a good fitting in the high–low moisture distribution areas with sensitivity ranged of 60–70 %. Finally, the moisture content map generated by Radarsat-1 was validated using a landslide inventory map. The resultant validation produced an area under curve of 0.704 (70 %)

Spatio-temporal analysis of oil spill impact and recovery pattern of coastal vegetation and wetland using multispectral satellite Landsat 8-OLI imagery and machine learning models

© 2020, by the authors. Oil spills are a global phenomenon with impacts that cut across socio-economic, health, and environmental dimensions of the coastal ecosystem. However, comprehensive assessment of oil spill impacts and selection of appropriate remediation approaches have been restricted due to reliance on laboratory experiments which offer limited area coverage and classification accuracy. Thus, this study utilizes multispectral Landsat 8-OLI remote sensing imagery and machine learning models to assess the impacts of oil spills on coastal vegetation and wetland and monitor the recovery pattern of polluted vegetation and wetland in a coastal city. The spatial extent of polluted areas was also precisely quantified for effective management of the coastal ecosystem. Using Johor, a coastal city in Malaysia as a case study, a total of 49 oil spill (ground truth) locations, 54 non-oil-spill locations and Landsat 8-OLI data were utilized for the study. The ground truth points were divided into 70% training and 30% validation parts for the classification of polluted vegetation and wetland. Sixteen different indices that have been used to monitor vegetation and wetland stress in literature were adopted for impact and recovery analysis. To eliminate similarities in spectral appearance of oil-spill-affected vegetation, wetland and other elements like burnt and dead vegetation, Support Vector Machine (SVM) and Random Forest (RF) machine learning models were used for the classification of polluted and nonpolluted vegetation and wetlands. Model optimization was performed using a random search method to improve the models' performance, and accuracy assessments confirmed the effectiveness of the two machine learning models to identify, classify and quantify the area extent of oil pollution on coastal vegetation and wetland. Considering the harmonic mean (F1), overall accuracy (OA), User's accuracy (UA), and producers' accuracy (PA), both models have high accuracies. However, the RF outperformed the SVM with F1, OA, PA and UA values of 95.32%, 96.80%, 98.82% and 95.11%, respectively, while the SVM recorded accuracy values of F1 (80.83%), OA (92.87%), PA (95.18%) and UA (93.81%), respectively, highlighting 1205.98 hectares of polluted vegetation and 1205.98 hectares of polluted wetland. Analysis of the vegetation indices revealed that spilled oil had a significant impact on the vegetation and wetland, although steady recovery was observed between 2015-2018. This study concludes that Chlorophyll Vegetation Index, Modified DifferenceWater Index, Normalized Difference Vegetation Index and Green Chlorophyll Index vegetation indices are more sensitive for impact and recovery assessment of both vegetation and wetland, in addition to Modified Normalized Difference Vegetation Index for wetlands. Thus, remote sensing and Machine Learning models are essential tools capable of providing accurate information for coastal oil spill impact assessment and recovery analysis for appropriate remediation initiatives

Landslide susceptibility mapping using certainty factor, index of entropy and logistic regression models in GIS and their comparison at Mugling-Narayanghat road section in Nepal Himalaya.

Landslide susceptibility maps are vital for disaster management and for planning development activities in the mountainous country like Nepal. In the present study, landslide susceptibility assessment of Mugling–Narayanghat road and its surrounding area is made using bivariate (certainty factor and index of entropy) and multivariate (logistic regression) models. At first, a landslide inventory map was prepared using earlier reports and aerial photographs as well as by carrying out field survey. As a result, 321 landslides were mapped and out of which 241 (75 %) were randomly selected for building landslide susceptibility models, while the remaining 80 (25 %) were used for validating the models. The effectiveness of landslide susceptibility assessment using GIS and statistics is based on appropriate selection of the factors which play a dominant role in slope stability. In this case study, the following landslide conditioning factors were evaluated: slope gradient; slope aspect; altitude; plan curvature; lithology; land use; distance from faults, rivers and roads; topographic wetness index; stream power index; and sediment transport index. These factors were prepared from topographic map, drainage map, road map, and the geological map. Finally, the validation of landslide susceptibility map was carried out using receiver operating characteristic (ROC) curves. The ROC plot estimation results showed that the susceptibility map using index of entropy model with AUC value of 0.9016 has highest prediction accuracy of 90.16 %. Similarly, the susceptibility maps produced using logistic regression model and certainty factor model showed 86.29 and 83.57 % of prediction accuracy, respectively. Furthermore, the ROC plot showed that the success rate of all the three models performed more than 80 % accuracy (i.e. 89.15 % for IOE model, 89.10 % for LR model and 87.21 % for CF model). Hence, it is concluded that all the models employed in this study showed reasonably good accuracy in predicting the landslide susceptibility of Mugling–Narayanghat road section. These landslide susceptibility maps can be used for preliminary land use planning and hazard mitigation purpose