基于要素分析和二元统计模型的区域滑坡危险等级制图——以国道212线陇南段为例

国道212线陇南段是我国地质灾害最发育的地区之一,绘制该区的滑坡危险等级地图对灾害管理和发展规划是极其必要的。基于滑坡的野外调查、机理研究和室内试验等工作,分析了滑坡与各种要素的相关性,选择控制滑坡的9个重要要素作为评价要素,利用GIS和二元统计的信息值模型和滑坡先验风险要素模型绘制了研究区的滑坡危险等级地图。最后,选用区内11个具有明显滑动位移的活动滑坡与滑坡危险等级地图比较,检验其可靠度。结果表明,活动的滑坡绝大部分都位于危险等级很高和高的范围内,说明两种模型的评价结果与研究区实际情况相吻合,同时也反映出信息值模型与实际情况更加相符。     

Rockfall hazard and risk analysis for Monte da Lua,Sintra, Portugal

The prediction of rockfall trajectories below a rock cliff is essential in susceptibility, hazard and risk maps, particularity close to populated areas. The Monte da Lua hill area in Portugal, a tourist destination close to the historic city of Sintra (UNESCO World Heritage), is a typical granite boulder chaos landscape where from time to time rockfalls occur, the last such event having occurred on 29 January 2002. This area is therefore suitable to develop a rockfall study in order to provide hazard and risk maps a basis for mitigation measures. A preliminary investigation of the area leads to the identification of 188 potentially dangerous boulders. Detailed locations and geotechnical characteristics in terms of geometry, strength and context were sampled for each boulder. Digital elevations at 1 × 1 m resolution, known rockfall trajectory and building locations are provided in a GIS project for the study together with the spatial database of boulder characteristics. The modelling approach was conducted in two steps: (1) discrimination of the boulders in terms of static and dynamic mobility behaviour with multivariate analysis; (2) stochastic simulation of rockfall trajectories. The rockfall trajectory algorithm proposed is straightforward and is only dependent on elevation data, initial location of boulders and a friction angle. Due to the slope of the area, it assumes that rockfall is always of the rolling or sliding type. The friction angle was calibrated on the basis of the rockfall travel distance recorded on 29 January 2002 and generates simulated “realistic” trajectories. A smaller friction angle increases all simulated trajectories, leading to more “pessimistic” scenarios. The combined analysis of trajectories and potential damage to buildings and discrimination in terms of static and dynamic behaviour provides a final table in which all 188 sampled boulders are classified in one of the five risk grades.     

Rockfall susceptibility zoning based on back analysis of rockfall deposit inventory in Gunung Kelir,Java

A rockfall susceptibility based on trajectory-energy/velocity approach needs release area or rockfall source. However, identification of rockfall source is not always possible for some areas in Indonesia. This paper presents a rockfall susceptibility zoning based on back analysis technique of rockfall deposit inventory in Gunung Kelir, Java. There were several steps in the rockfall susceptibility zoning: (1) rockfall deposit inventory, (2) rockfall simulation based on back analysis of rockfall deposit inventory, (3) sensitivity analysis, and (4) rockfall susceptibility zoning. The result suggests that the travel distance is affected by the spatial distribution of rockfall source, lithology or surface material, and topography (angle of slope and angle of aspect). Final trajectories were employed to generate landslide susceptibility map which may allow a policy maker to have an advanced consideration to achieve specified risk measures and evaluation of their cost efficiency to optimize budget and design. Application of rockfall susceptibility zoning based on back analysis of rockfall deposits is efficient where rockfall source information is unavailable.     

Rockfall detachment susceptibility map in El Hierro Island, Canary Islands, Spain

A semi-quantitative heuristic methodology is developed to map a rockfall detachment susceptibility zonation of El Hierro Island (Canary Archipelago). The rationalized procedure, which we called non-weighted bounded indicators, is based on overlapping thematic maps of conditioning factors to mass movement, which are appropriately and individually rescaled and then composed by addition to obtain a susceptibility numerical index through a GIS. As the consistency of the geomorphological analysis depends on the expert subjective criteria and the appropriate interpretation of the landscape, the use of this methodology reduces subjectivity and quantifies the degree of susceptibility. The main factors affecting the mass movement phenomena (rockfalls events), also recognized in the field and, therefore, considered in the presented GIS arrangement, are slope, profile curvature, lithology, vegetation cover and dykes density. To calculate the slope threshold or minimum angle characteristic of rockfall source areas, mixed Gaussian slope frequency decomposition is used. The curvature index reveals stepwise areas. Qualitative geomechanical characteristics are linked to a quantitative index according to a volcanic lithological-complexes classification. Both destabilization (root-wedging) and stabilization effects are considered into the vegetation cover index. The dyke density index incorporates the bearing rock capacity decrease produced in the halo around a dyke network intrusion. Slope, curvature and vegetation indexes thresholds have been fitted following field observations. A rockfall detachment susceptibility map is obtained and classified based on the histogram maxima. The rockfall inventory, based on rockfall events reported within the island, was used for the model validation. A 12?% of the whole island shows medium to very high susceptibility.     

Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs

Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.     

编目图数据形式对区域崩塌易发性评价结果的影响

本文以北京市怀柔区为例,通过现场调查,对688处崩塌灾害分别以面数据和点数据的形式获取了两套编目图。根据现场调查和资料分析,选取岩性、地形、断裂和道路建设作为该区崩塌灾害的主控因素,采用频率比(FR)模型对崩塌灾害的易发性进行了评价。为了对评价结果的预测性进行检验,采用随机分割法,选取了415处崩塌用于频率比模型的计算,剩余的273处崩塌用于评价结果预测性的验证。预测曲线表明,基于崩塌面数据的评价结果比基于点数据的评价结果具有明显的优越性。根据基于面数据的频率比模型评价结果,可以将研究区的崩塌灾害易发性划分为5个等级:较低易发(占全区14%)、低易发(占全区20%)、中等易发(占全区27%)、高易发(占全区22%)和极高易发(占全区17%)。相关工作和结论可以为区域地质灾害易发性评价中编目图的编制提供参考,并为怀柔区区域国土利用和防灾减灾提供指导。     

基于GIS的白龙江流域甘肃段滑坡易发性评价

在甘肃省白龙江流域地质灾害资料收集及现场调查的基础上, 统计分析了该区滑坡发育与地层岩性、坡度、坡向、高程、断裂、植被等因素之间的关系, 建立了白龙江流域滑坡易发性评价指标体系。采用基于GIS的层次分析法评价模型, 完成了滑坡易发性分区评价, 将研究区滑坡按易发程度划分为高易发区、中易发区、低易发区和极低易发区, 其中, 高易发区占研究区总面积的13.59%, 主要分布在断裂带、白龙江两侧以及软弱岩土体分布的区域; 中易发区占27.85%;主要分布在白龙江支流以及主要道路两侧的一定范围内; 低易发区占33.09%, 主要分布在海拔相对较高、植被覆盖度较高、基本上无断裂带通过的区域; 其余区域为极低易发区, 占25.46%。对比分析显示评价结果与实际滑坡发育情况吻合, 可以较好地反映区内滑坡灾害发育的总体特征。      

Landslide susceptibility assessment using the bivariate statistical analysis and the index of entropy in the Sibiciu Basin (Romania)

The Sibiciu Basin is located in Romania between the Buzău Mountains and the Buzau Subcarpathians (Curvature Carpathians and Subcarpathians). The geology of the basin consists of Paleogene flysch deposits represented by an alternation of sandstones, marls, clays and schists and Neogene deposits represented by marls, clays and sands. The area is affected by different types of landslides (shallow, medium-deep and deep-seated failures). In Romania, in the last decades, direct and indirect methods have been applied for landslide susceptibility assessment. The most utilized before 2000 were based on qualitative approaches. This study evaluates the landslide susceptibility in the Sibiciu Basin using a bivariate statistical analysis and an index of entropy. A landslide inventory map was prepared, and a susceptibility estimate was assessed based on the following parameters which influence the landslide occurrence: slope angle, slope aspect, curvature, lithology and land use. The landslide susceptibility map was divided into five classes showing very low to very high landslide susceptibility areas.     

无人机影像在高陡边坡危岩体调查中的应用

在高陡边坡危岩体的调查中,复杂的地形条件经常限制工作的正常开展,如何快速准确地获取地质灾害信息一直是地质灾害调查研究中的难点之一。以往的研究中对无人机遥感技术在黄土、高原等地区应用有所报道,但对西南地区高陡边坡危岩体灾害调查的研究尚无报道。文章以锦屏二级水电站出线场边坡落石灾害所在区域为例,将无人机摄影测量技术应用于高陡边坡危岩体调查中,通过无人机倾斜摄影获取高分辨率遥感影像,开展遥感影像三维建模,进行地质灾害遥感解译,总结了无人机遥感系统在高陡边坡危岩体调查的技术流程。通过三维实景模型,精确地分析了落石灾害的空间分布、失稳模式及演化过程,查明了区域内危岩隐患点的分布特征;基于三维点云模型,提取出地质灾害体的属性信息,测得落石方量为11.7 m3,采用最小二乘法进行平面拟合,得到落石两组主控结构面产状为275.4°∠31.2°、103.5°∠63.3°。实践表明,无人机遥感技术在高陡边坡地区落石灾害调查中具有明显的可行性和优越性,可以较好地应用于高陡边坡危岩体调查中。     

Rockfall susceptibility zoning at a large scale: From geomorphological inventory to preliminary land use planning

Rockfall, up to several hundreds of cubic meters, is a frequent and rapid landslide which menaces extensive areas in mountainous territories. Rockfall susceptibility zoning map at a large scale (1:5000–1:25 000) can be the first tool for land use planning in order to manage rockfall risk. A methodology allowing to analyze susceptibility in extensive areas with optimum cost/benefit relationship is needed. This work analyzes rockfall susceptibility in an extensive rocky mountain of the Principality of Andorra (Pyrenees Mountains), first on the rock slope and then on the exposed area located below. The rockfall record, obtained by means of geomorphological analysis, supplies the main data to analyze the susceptibility on the rock slope. An additional historical inventory verifies the accuracy of rockfall sizes recorded by means of the geomorphological analysis. According to the classification recommended by the Guidelines of Joint Technical Committee, the density of rockfall features on the rock slope assesses susceptibility in four levels. Subsequently, susceptibility on exposed areas has been analyzed by means of reach probability of rock blocks analysis using empirical models. Data acquired from thirteen recent events, from 1999 to 2004, have been used to verify the accuracy of the two empirical models mainly used (reach angle and shadow angle). Five reach probability limits (1, 0.5, 0.25, 0.01, and 0) establish boundaries between susceptibility levels. The resulting rockfall susceptibility zoning map allows: (a) to identify land areas and human elements exposed to rockfalls and, (b) to establish several exposition levels. This map can be a useful and cost-effective tool for administrations responsible to manage natural risk in order to guide urban grow in extensive areas or decide upon work programs based on in-depth analysis (hazard and risk).     

Assessment of landslide susceptibility on the natural terrain of Lantau Island, Hong Kong

 Steep terrain and the high frequency of tropical rainstorms make landslide occurrence on natural terrain a common phenomenon in Hong Kong. For example, more than 800 slope failures were triggered by a rainstorm in November 1993 on Lantau Island, Hong Kong. Maps of recent landslides interpreted from aerial photographs, in combination with a geographical information system, were used to evaluate the frequency and distribution of landslides, with particular reference to such physical parameters as lithology, slope gradient, slope aspect, elevation, vegetation cover, and proximity to drainage line, all of which are considered to be influential in the occurrence of landslides. A stepwise logistic regression model was obtained between landslide susceptibility and the above mentioned physical parameters. The study area has been classified into five classes of relative landslide susceptibility, namely, very low, low, moderate, high, and very high, based on this methodology. Received: 17 December 1999 · Accepted: 21 March 2000     

Assessing rockfall susceptibility in steep and overhanging slopes using three-dimensional analysis of failure mechanisms

Rockfalls strongly influence the evolution of steep rocky landscapes and represent a significant hazard in mountainous areas. Defining the most probable future rockfall source areas is of primary importance for both geomorphological investigations and hazard assessment. Thus, a need exists to understand which areas of a steep cliff are more likely to be affected by a rockfall. An important analytical gap exists between regional rockfall susceptibility studies and block-specific geomechanical calculations. Here we present methods for quantifying rockfall susceptibility at the cliff scale, which is suitable for sub-regional hazard assessment (hundreds to thousands of square meters). Our methods use three-dimensional point clouds acquired by terrestrial laser scanning to quantify the fracture patterns and compute failure mechanisms for planar, wedge, and toppling failures on vertical and overhanging rock walls. As a part of this work, we developed a rockfall susceptibility index for each type of failure mechanism according to the interaction between the discontinuities and the local cliff orientation. The susceptibility for slope parallel exfoliation-type failures, which are generally hard to identify, is partly captured by planar and toppling susceptibility indexes. We tested the methods for detecting the most susceptible rockfall source areas on two famously steep landscapes, Yosemite Valley (California, USA) and the Drus in the Mont-Blanc massif (France). Our rockfall susceptibility models show good correspondence with active rockfall sources. The methods offer new tools for investigating rockfall hazard and improving our understanding of rockfall processes.     

Multi-scenario Rockfall Hazard Assessment Using LiDAR Data and GIS

Transportation corridors that pass through mountainous or hilly areas are prone to rockfall hazard. Rockfall incidents in such areas can cause human fatalities and damage to properties in addition to transportation interruptions. In Malaysia, the North–South Expressway is the most significant expressway that operates as the backbone of the peninsula. A portion of this expressway in Jelapang was chosen as the site of rockfall hazard assessment in multiple scenarios. Light detection and ranging techniques are indispensable in capturing high-resolution digital elevation models related to geohazard studies. An airborne laser scanner was used to create a high-density point cloud of the study area. The use of 3D rockfall process modeling in combination with geographic information system (GIS) is a beneficial tool in rockfall hazard studies. In this study, a 3D rockfall model integrated into GIS was used to derive rockfall trajectories and velocity associated with them in multiple scenarios based on a range of mechanical parameter values (coefficients of restitution and friction angle). Rockfall characteristics in terms of frequency, height, and energy were determined through raster modeling. Analytic hierarchy process (AHP) was used to compute the weight of each rockfall characteristic raster that affects rockfall hazard. A spatial model that considers rockfall characteristics was conducted to produce a rockfall hazard map. Moreover, a barrier location was proposed to eliminate rockfall hazard. As a result, rockfall trajectories and their characteristics were derived. The result of AHP shows that rockfall hazard was significantly influenced by rockfall energy and then by frequency and height. The areas at risk were delineated and the hazard percentage along the expressway was observed and demonstrated. The result also shows that with increasing mechanical parameter values, the rockfall trajectories and their characteristics, and consequently rockfall hazard, were increased. In addition, the suggested barrier effectively restrained most of the rockfall trajectories and eliminated the hazard along the expressway. This study can serve not only as a guide for a comprehensive investigation of rockfall hazard but also as a reference that decision makers can use in designing a risk mitigation method. Furthermore, this study is applicable in any rockfall study, especially in situations where mechanical parameters have no specific values.     

Landslide susceptibility mapping based on GIS and information value model for the Chencang District of Baoji,China

The main objective of this study was to apply a statistical (information value) model using geographic information system (GIS) to the Chencang District of Baoji, China. Landslide locations within the study area were identified using reports and aerial photographs, and a field survey. A total of 120 landslides were mapped, of which 84 (70 %) were randomly selected for building the landslide susceptibility model. The remaining 36 (30 %) were used for model validation. We considered a total of 10 potential factors that predispose an area to a landslide for the landslide susceptibility mapping. These included slope degree, altitude, slope aspect, plan curvature, geomorphology, distance from faults, lithology, land use, mean annual rainfall, and peak ground acceleration. Following an analysis of these factors, a landslide susceptibility map was produced using the information value model with GIS. The resulting landslide susceptibility index was divided into five classes (very high, high, moderate, low, and very low) using the natural breaks method. The corresponding distribution area percentages were 29.22, 25.14, 15.66, 15.60, and 14.38 %, respectively. Finally, landslide locations were used to validate the results of the landslide susceptibility map using areas under the curve (AUC). The AUC plot showed that the susceptibility map had a success rate of 81.79 % and a prediction accuracy of 82.95 %. Based on the results of the AUC evaluation, the landslide susceptibility map produced using the information value model exhibited good performance.     

边坡滚石运动轨迹分段循环算法

边坡滚石是在基本建设中常遇到的一种工程灾害。其防护结构的设计依据主要是边坡滚石的运动轨迹,因而要对这种灾害进行防治,首先要确定滚石的运动轨迹。根据落石运动全过程的三个阶段,即滚动(滑动)运动阶段、抛物运动阶段和与坡面碰撞阶段,利用分段循环算法分别得到了三个运动阶段的运动速度计算公式。以某岩质边坡构建相似模型,将该理论公式应用于落石运动轨迹预测,验证了公式的可行性和有效性。与现有方法相比,分段循环算法具有计算原理简单明确、易于使用等优点。该方法可用于斜坡区滚石运动特征的预测,并可作为斜坡区滚石灾害防治的依据。     

Landslide susceptibility mapping by frequency ratio and logistic regression methods: an example from Niksar–Resadiye (Tokat,Turkey)

Turkey confronts loss of life and large economic losses due to natural disasters caused by its morphologic structure, geographical placement, and climate characteristics. The Kuzulu (Koyulhisar) landslide, which caused loss of life and property on 17th March 2005, occurred in an area near the country’s most important active fault, the North Anatolian Fault Zone. To mitigate and prevent landslide damages, prediction of landslide susceptibility areas based on probabilistic methods has a great importance. The purpose of this study was to produce a landslide susceptibility map by the logistic regression and frequency ratio methodologies for a 733-km² area near the North Anatolian Fault Zone from the southeast of Niksar to Resadiye in Tokat province. Conditioning parameters, such as elevation, slope gradient, slope aspect, distance to streams, roads, and faults, drainage density, and fault density, were used in the analysis. Before susceptibility analysis, the landslides observed in the area were separated into two groups for use in analysis and verification, respectively. The susceptibility maps produced had five different susceptibility classes such as very low, low, moderate, high, and very high. To test the performance of the susceptibility maps, area under curve (AUC) approach was used. For the logistic regression method, the AUC value was 0.708; while for the frequency rate method, this value was 0.744. According to these AUC values, it could be concluded that the two landslide susceptibility maps obtained were successful.     

Landslide susceptibility analysis with a bivariate approach and GIS in Northern Iran

Globally, landslides cause hundreds of billions of dollars in damage and hundreds of thousands of deaths and injuries each year. A landslide susceptibility map describes areas where landslides are likely to occur in the future by correlating some of the principal factors that contribute to landslides with the past distribution of landslides. A case study is conducted in the mountainous northern Iran. In this study, a landslide susceptibility map of the study area was prepared using bivariate method with the help of the geographic information system. Area density (bivariate) method was used to weight landslide-influencing data layers. 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 hazard classes: very low, low, moderate, high, and very high. The percentage distribution of landslide susceptibility degrees was calculated. It was found that about 26% of the study area is classified as very high and high hazard classes.     

A susceptibility mapping model of earthquake-triggered slope geohazards based on geo-spatial data in mountainous regions

The main purpose of this paper is to develop a susceptibility mapping model of earthquake-trigged slope geohazards based on geo-spatial data in mountainous region. A rapid-response mapping model of slope geohazards triggered by earthquakes that can be widely applied in practice can be constructed by a simple method of weighted superposition of related accessible factors. Based on site investigations into geological disasters in Wenchuan County caused by the Wenchuan earthquake on 12 May 2008, the main influencing factors of the earthquake triggering slope geohazards were studied, including peak ground acceleration (PGA) and 6 geo-spatial data factors: distance from streams, distance from highway, slope gradient, slope position, normalised difference vegetation index (NDVI) and micro-landforms. The order of the relative importance of each factor was determined by regression analysis, then the weight of each factor contributing to slope geohazards was calculated by using the analytic hierarchy process (AHP). The susceptibility mapping model of slope geohazards was derived from the map overlaying principle in ArcGIS platform, and finally, the results were re-classified to obtain the map concerning susceptibility areas of slope geohazards. The results, by comparing the geohazards points with the field survey, proved that over 63% geohazards points were located in high and very high susceptibility areas, while only about 14% geohazards points were located in very low and low susceptibility areas. The model showed a relative agreement between the points on the susceptibility mapping and onsite points of slope geohazards. Above all, predictions about the slope geohazards at meizoseismal areas of earthquakes occurring in mountainous regions can be based on geo-spatial data; so the model could be applied to study other areas which may be prone to the slope geohazards during the earthquakes.     

Landslide susceptibility zonation using GIS and evidential belief function model

The main goal of this paper is to generate a landslide susceptibility map through evidential belief function (EBF) model by using Geographic Information System (GIS) for Qianyang County, Shaanxi Province, China. At first, a detailed landslide inventory map was prepared, and the following ten landslide-conditioning factors were collected: slope angle, slope aspect, curvature, plan curvature, profile curvature, altitude, distance to rivers, geomorphology, lithology, and rainfall. The landslides were detected from the interpretation of aerial photographs and supported by field surveys. A total of 81 landslides were randomly split into the following two parts: the training dataset 70 % (56 landslides) were used for establishing the model and the remaining 30 % (25 landslides) were used for the model validation. The ArcGIS was used to analyze landslide-conditioning factors and evaluate landslide susceptibility; as a result, a landslide susceptibility map was generated by using EBF and ArcGIS 10.0, thus divided into the following five susceptibility classes: very low, low, moderate, high, and very high. Finally, when we validated the accuracy of the landslide susceptibility map, both the success-rate and prediction-rate curve methods were applied. The results reveal that a final susceptibility map has the success rate of 83.31 % and the prediction rate of 79.41 %.     

Rockfall hazard assessment along a road in the Sorrento Peninsula,Campania, southern Italy

Rockfalls are common in the steep and vertical slopes of the Campania carbonate massifs and ridges, and frequently represent the main threat to the anthropogenic environment, potentially damaging urban areas, scattered houses, roads, etc. Despite the generally limited volumes involved, the high velocity of movement (from few to tens of metres per second) poses rockfalls among the most dangerous natural hazards to man. Evaluating the rockfall hazard is not an easy task, due to the high number of involved factors, and particularly to the difficulty in determining the properties of the rock mass. In this paper, we illustrate the assessment of the rockfall hazard along a small area of the Sorrento Peninsula (Campania region, southern Italy). Choice of the site was determined by the presence of a road heavily frequented by vehicles. In the area, we have carried out detailed field surveys and software simulations that allow generating simple rockfall hazard maps. Over twenty measurement stations for geo-mechanical characterization of the rock mass have been distributed along a 400-m-long slope of Mount Vico Alvano. Following the internationally established standards for the acquisition of rock mass parameters, the main kinematics have been recognized, and the discontinuity families leading to the different failures identified. After carrying out field experiments by artificially releasing a number of unstable blocks on the rock cliff, the rockfall trajectories along the slope were modelled using 2-D and 3-D programs for rockfall analysis. The results were exploited to evaluate the rockfall hazard along the threatened element at risk.