Using wavelet tools to analyse seasonal variations from InSAR time-series data: a case study of the Huangtupo landslide

Synthetic aperture radar interferometry (InSAR) has proven to be a powerful tool for monitoring landslide movements with a wide spatial and temporal coverage. Interpreting landslide displacement time-series derived from InSAR techniques is a major challenge for understanding relationships between triggering factors and slope displacements. In this study, we propose the use of various wavelet tools, namely, continuous wavelet transform (CWT), cross wavelet transform (XWT) and wavelet coherence (WTC) for interpreting InSAR time-series information for a landslide. CWT enables time-series records to be analysed in time-frequency space, with the aim of identifying localized intermittent periodicities. Similarly, XWT and WTC help identify the common power and relative phase between two time-series records in time-frequency space, respectively. Statistically significant coherence and confidence levels against red noise (also known as brown noise or random walk noise) can be calculated. Taking the Huangtupo landslide (China) as an example, we demonstrate the capabilities of these tools for interpreting InSAR time-series information. The results show the Huangtupo slope is affected by an annual displacement periodicity controlled by rainfall and reservoir water level. Reservoir water level, which is completely regulated by the dam activity, is mainly in ‘anti-phase’ with natural rainfall, due to flood control in the Three Gorges Project. The seasonal displacements of the Huangtupo landslide is found to be ‘in-phase’ with respect to reservoir water level and the rainfall towards the front edge of the slope and to rainfall at the higher rear of the slope away from the reservoir.     

Accuracy assessment of InSAR derived input maps for landslide susceptibility analysis: a case study from the Swiss Alps

In recent years SAR interferometry has become a widely used technique for measuring altitude and displacement of the surface of the earth. Both these capabilities are highly relevant for landslide susceptibility studies. Although there are many problems that make the use of SAR interferometry less suitable for landslide inventory mapping, it’s use in landslide monitoring and in the generation of input maps for landslide susceptibility assessment looks very promising. The present work attempts to evaluate the usefulness and limitations of this technique based on a case study in the Swiss Alps. Input maps were generated from ERS repeat pass data using SAR interferometry. A land cover map has been generated by image classification of multi-temporal SAR intensity images. An InSAR DEM was generated and a number of maps were derived from it, such as slope-, aspect, altitude- and slope form classes. These maps were used to generate landslide and rockfall susceptibility maps, which give fairly well acceptable results. However, a comparison of the InSAR DEM with the conventional Swisstopo DEM, indicated significant errors in the absolute height and slope angles derived from InSAR, especially along the ridges and in the valleys. These errors are caused by low coherence mostly due to layover and shadow effects. Visual comparison of stereo images created from hillshading maps and corresponding DEMs demonstrate that a considerable amount of topographic details have been lost in the InSAR-derived DEM. It is concluded that InSAR derived input maps are not ideal for landslide susceptibility assessment, but could be used if more accurate data is lacking.     

Shallow landslide susceptibility assessment under future climate and land cover changes: A case study from southwest China

There is no doubt that land cover and climate changes have consequences on landslide activity, but it is still an open issue to assess and quantify their impacts. Wanzhou County in southwest China was selected as the test area to study rainfall-induced shallow landslide susceptibility under the future changes of land use and land cover (LULC) and climate. We used a high-resolution meteorological precipitation dataset and frequency distribution model to analyse the present extreme and antecedent rainfall conditions related to landslide activity. The future climate change factors were obtained from a 4-member multi-model ensemble that was derived from statistically downscaled regional climate simulations. The future LULC maps were simulated by the land change modeller (LCM) integrated into IDRISI Selva software. A total of six scenarios were defined by considering the rainfall (antecedent conditions and extreme events) and LULC changes towards two time periods (mid and late XXI century). A physically-based model was used to assess landslide susceptibility under these different scenarios. The results showed that the magnitude of both antecedent effective recharge and event rainfall in the region will evidently increase in the future. Under the scenario with a return period of 100 years, the antecedent rainfall in summer will increase by up to 63% whereas the event rainfall will increase by up to 54% for the late 21st century. The most considerable changes of LULC will be the increase of forest cover and the decrease of farming land. The magnitude of this change can reach + 22.1% (forest) and –9.2% (farmland) from 2010 until 2100, respectively. We found that the negative impact of climate change on landslide susceptibility is greater than the stabilizing effect of LULC change, leading to an over decrease in stability over the study area. This is one of the first studies across Asia to assess and quantify changes of regional landslide susceptibility under scenarios driven by LULC and climate change. Our results aim to guide land use planning and climate change mitigation considerations to reduce landslide risk.     

一种土的阻尼比近似计算方法

提出了一种用多边形逼近滞回曲线计算土的阻尼比的方法,该方法可直接利用试验数据进行计算,算法可重复利用,使得阻尼比的计算工作得以简化。将其与传统的椭圆拟合方法计算结果进行了对比,由于实测的滞回圈不是标准的椭圆,两者的计算结果存在一定的差异。     

利用InSAR技术研究黄土地区滑坡分布

InSAR技术能够获取大面积、连续、高精度的地表垂直形变信息,可用来监测地震、火山、滑坡等自然灾害造成的地表形变。文章介绍了InSAR技术在监测陕北黄土地区滑坡中的应用,首先进行野外地质勘察和TM光学遥感影像解译,接着通过EnviSat SAR数据差分干涉处理,获取研究区干涉形变场,提取出滑坡位移量,最后详细分析黄草湾至董家寺沿线一带的滑坡变形范围,并划定出了4个有一定变形的重点监视区。     

A statistical approach to landslide risk modelling at basin scale: from landslide susceptibility to quantitative risk assessment

A procedure for landslide risk assessment is presented. The underlying hypothesis is that statistical relationships between past landslide occurrences and conditioning variables can be used to develop landslide susceptibility, hazard and risk models. The latter require also data on past damages. Landslides occurred during the last 50 years and subsequent damages were analysed. Landslide susceptibility models were obtained by means of Spatial Data Analysis techniques and independently validated. Scenarios defined on the basis of past landslide frequency and magnitude were used to transform susceptibility into quantitative hazard models. To assess vulnerability, a detailed inventory of exposed elements (infrastructures, buildings, land resources) was carried out. Vulnerability values (0–1) were obtained by comparing damages experienced in the past by each type of element with its actual value. Quantitative risk models, with a monetary meaning, were obtained for each element by integrating landslide hazard and vulnerability models. Landslide risk models showing the expected losses for the next 50 years were thus obtained for the different scenarios. Risk values obtained are not precise predictions of future losses but rather a means to identify areas where damages are likely to be greater and require priority for mitigation actions.     

基于XGBoost模型的三峡库区燕山乡滑坡易发性评价与区划

滑坡易发性评价是精细化滑坡灾害风险评价的基础。为了提升滑坡易发性评价模型的精度和稳健性,以三峡库区万州区燕山乡为例,选取工程地质岩组、堆积层厚度等九个影响因子构建滑坡易发性评价指标体系,应用信息量模型定量分析滑坡发育与指标之间的关系。在此基础上,随机选取70%/30%的滑坡样本作为训练/验证数据集,应用极致梯度提升模型(extreme gradient boosting, XGBoost)开展易发性评价。随后从模型预测精度和模型稳定性两方面将其与决策树模型(decision tree, DT)和梯度提升树模型(gradient boosting decision tree, GBDT)进行对比。结果表明：研究区堆积层滑坡主要受长江水系、堆积层厚度和工程地质岩组影响。XGBoost模型具有最高的准确率(94.3%)和预测精度(97.3%)。在模型稳定性验证中,平均预测精度最高(97.3%),优于DT(91.3%)和GBDT(95.7%),模型标准差和变异系数均为0.01,低于其余两种模型。XGBoost在区域滑坡易发性评价与制图中得到了可靠的结果,为滑坡灾害空间预测提供了新的技术支撑。     

Evaluation of the consistency of landslide susceptibility mapping: a case study from the Kankai watershed in east Nepal

GIS-based landslide susceptibility maps for the Kankai watershed in east Nepal are developed using the frequency ratio method and the multiple linear regression technique. The maps are derived from comparing observed landslides with possible causative factors: slope angle, slope aspect, slope curvature, relative relief, distance from drainage, land use, geology, distance from faults and mean annual rainfall. The consistency of the maps is evaluated using landslide density analysis, success rate analysis and spatially agreed area approach. The first two analyses produce almost identical quantitative results, whereas the last approach is able to reveal spatial differences between the maps and also to improve predictions in the agreed high landslide-susceptible area.     

Geo-diversity and its hydrological response in relation to landslide susceptibility in the Himalaya: a GIS-based case study

Assessment and inventory of landslide susceptibility are essential for the formulation of successful disaster mitigation plans. The objective of this study was to assess landslide susceptibility in relation to geo-diversity and its hydrological response in the Lesser Himalaya with a case study using Geographic Information System (GIS) technology. The Dabka watershed, which constitutes a part of the Kosi Basin in the Lesser Himalaya, India, in the district of Nainital, has been selected for the case illustration. The study constitutes three GIS modules: geo-diversity informatics, hydro informatics and landslide informatics. Through the integration and superimposing of spatial data and attribute data of all three GIS modules, Landslide Susceptibility Index (LSI) has been prepared to identify the level of susceptibility for landslide hazards. This resonance study, carried out over a period of five years (2007–2011), found that areas of most stressed geo-diversity (comprising very steep slopes above 30°, geology of Lower Krol and Lariakanta formation, geomorphology of moist areas and debris sites, land use of barren land with a very high drainage frequency and spring density) have a high landslide susceptibility because of high rate of average runoff (33 l/s/km²), flood magnitude (307.28 l/s/km²), erosion (398 tons/km²) and landslide density (5–10 landslides/km²). The areas of least stressed geo-diversity (comprising gentle slopes below 10°, geology of Kailakhan and Siwalik formation, geomorphology of depositional terraces, land use of dense forest with low drainage frequency and spring density) have the lowest landslide susceptibility because of the low rate of average runoff (6.27 l/s/km²), flood magnitude (20.49 l/s/km²), erosion (65.80 tons/km²) and landslide density (1–2 landslides/km²).     

PFR model and GiT for landslide susceptibility mapping: a case study from Central Alborz, Iran

In northern parts of Iran such as the Alborz Mountain belt, frequent landslides occur due to a combination of climate and geologic conditions with high tectonic activities. This results in millions of dollars of financial damages annually excluding casualties and unrecoverable resources. This paper evaluates the landslide susceptible areas in Central Alborz using the probabilistic frequency ratio (PFR) model and Geo-information Technology (GiT). The landslide location map in this study has been generated based on image elements interpreted from IRS satellite data and field observations. The display, manipulation and analysis have been carried out to evaluate layers such as geology, geomorphology, soil, slope, aspect, land use, distance from faults, lineaments, roads and drainages. The validation group of actual landslides and relative operation curve method has been used to increase the accuracy of the final landslide susceptibility map. The area under the curve evaluates how well the method predicts landslides. The results showed a satisfactory agreement of 91% between prepared susceptibility map and existing data on landslide locations.     

Neural networks and landslide susceptibility: a case study of the urban area of Potenza

For those working in the field of landslide prevention, the estimation of hazard levels and the consequent production of thematic maps are principal objectives. They are achieved through careful analytical studies of the characteristics of landslide prone areas, thus, providing useful information regarding possible future phenomena. Such maps represent a fundamental step in the drawing up of adequate measures of landslide hazard mitigation. However, for a complete estimation of landslide hazard, meant as the degree of probability that a landslide occurs in a given area, within a given space of time, detailed and uniformly distributed data regarding their incidence and causes are required. This information, while obtainable through laborious historical research, is usually partial, incomplete and uneven, and hence, unsatisfactory for zoning on a regional scale. In order to carry this out effectively, the utilization of spatial estimation of the relative levels of landslide hazard in the various areas was considered opportune. These areas were classified according to their levels of proneness to landslide activity without taking recurrence periods into account. Various techniques were developed in order to obtain upheaval numerical estimates. The method used in this study, which was applied in the area of Potenza, is based on techniques derived from artificial intelligence (Artificial Neural Network—ANN). This method requires the definition of appropriate thematic layers, which parameterize the area under study. These are recognized by means of specific analyses in a functional relationship to the event itself. The parameters adopted are: slope gradient, slope aspect, topographical index, topographical shape, elevation, land use and lithology.     

An inverse method to determine the optimal stress from imperfect fault data

An inverse method is developed to determine the optimal stress inversion from imperfect fault data. It is specifically designed to process fault data at variable division—either with observed slip directions or with observed fault/slip senses or with both. The method has the flexibility of processing data with a few fault/slip senses incorrectly determined—a case we occasionally meet in the field. A vast number of artificial 10-datum examples are processed to illustrate the reliability of the method and to examine the effects of the data set at variable division on stress estimation. The results show that the precision of solution decreases with the increase in either of two factors, the percentage of the data limited to slip sense and the range of measurement error, both resulting in a broader dispersion of estimated stress vectors relative to the assigned stress vector in the sigma space. Thus, in order to keep the dispersion within a range of 5° requires the percentage of data for which only slip sense is known to be less than about 20%. This in turn allows for a measurement error range in the data of 10° in real space, much larger than real measurement errors. We believe that the proposed method would be of practical value for stress inversion of fault/slip data.     

GIS based landslide susceptibility mapping using a fuzzy logic approach: A case study from Ghurmi-Dhad Khola area, Eastern Nepal

Landslides cause extensive loss of life and property in the Nepal Himalaya. Since the late 1980s, different mathematical models have been developed and applied for landslide susceptibility mapping and hazard assessment in Nepal. The main goal of this paper is to apply fuzzy logic to landslide susceptibility mapping in the Ghurmi-Dhad Khola area, Eastern Nepal. Seven causative factors are considered: slope angle, slope aspect, distance from drainage, land use, geology, distance from faults and folds, soil and rock type. Likelihood ratios are obtained for each class of causative factors by comparison with past landslide occurrences. The ratios are normalized between zero and one to obtain fuzzy membership values. Further, different fuzzy operators are applied to generate landslide susceptibility maps. Comparison with the landslide inventory map reveals that the fuzzy gamma operator with a γ-value of 0.60 yields the best prediction accuracy. Consequently, this operator is used to produce the final landslide susceptibility zonation map.     

Influence of seismic acceleration on landslide susceptibility maps: a case study from NE Turkey (the Kelkit Valley)

Particularly in the last decade, landslide susceptibility and hazard maps have been used for urban planning and site selection of infrastructures. Most of the procedures for preparing of landslide susceptibility maps need high-quality landslide inventory map. Although the rainfall and seismic activities are accepted as triggering factor for landslides, designation of the triggering factor for each landslide in the inventory is almost impossible when well-documented records are unavailable. Therefore, during preparation of landslide susceptibility map, whole landslide records in the inventory map are used together without classifying based on the triggering factors. Although seismic activity is accepted as a triggering factor, possible effect of the use of seismic activity on production of landslide susceptibility map was investigated in this study, and the subject is open to discussion. For this purpose, a series of stability analyses based on circular failure and infinite slope model were performed considering different pseudostatic conditions. The results of analyses show that gentle slopes have higher susceptibility to failure than steeper ones, even if their stability conditions (susceptibilities) are similar for static condition. The seismic forces acting on failure surfaces may not be sufficiently taken into consideration in the conventionally prepared landslide susceptibility maps. Employing the general decreasing trend in stability condition based on slope face angle and the seismic acceleration, a new procedure was introduced for preparing of the landslide susceptibility map for a scenario earthquake. The prediction performance of occurring landslides increased after the procedure was applied to the conventionally prepared landslide susceptibility map. According to the threshold independent spatial performance analyses of the proposed methodology and the produced landslide susceptibility maps, the area under ROC curve values were calculated as 0.801, 0.933, and 0.947 for the maps prepared by considering conventional method and scenario earthquakes having M _w values of 5.5 and 7.5, respectively.     

A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

A heuristic approach to global landslide susceptibility mapping

Landslides can have significant and pervasive impacts to life and property around the world. Several attempts have been made to predict the geographic distribution of landslide activity at continental and global scales. These efforts shared common traits such as resolution, modeling approach, and explanatory variables. The lessons learned from prior research have been applied to build a new global susceptibility map from existing and previously unavailable data. Data on slope, faults, geology, forest loss, and road networks were combined using a heuristic fuzzy approach. The map was evaluated with a Global Landslide Catalog developed at the National Aeronautics and Space Administration, as well as several local landslide inventories. Comparisons to similar susceptibility maps suggest that the subjective methods commonly used at this scale are, for the most part, reproducible. However, comparisons of landslide susceptibility across spatial scales must take into account the susceptibility of the local subset relative to the larger study area. The new global landslide susceptibility map is intended for use in disaster planning, situational awareness, and for incorporation into global decision support systems.     

一种确定导管架平台群桩p-y数据的方法

导管架海洋平台群桩基础的桩头约束是一种弹性约束。针对具有弹性约束桩头的群桩,提出了一种利用非线性地基梁群桩模型,通过迭代计算确定导管架平台群桩p-y数据的方法。该方法首先依据单桩p-y曲线,利用具有弹性约束桩头的群桩模型,计算桩头在总荷载作用下的水平位移;再利用Poulos相互作用系数确定由于群桩相互作用引起的桩头附加水平位移,将总荷载作用下的桩头水平位移与桩头附加水平位移叠加后的结果作为迭代计算的初始桩头水平位移。然后,依据桩头荷载与初始桩头水平位移,通过对单桩p-y数据的标定,确定其修正系数Ym,进而得到与第1次计算对应的群桩p-y数据。在接下来的迭代计算中,利用每次更新后的Ym确定该次计算使用的群桩p-y数据,并据此由群桩计算模型计算桩头位移,通过对单桩p-y数据标定确定相应的Ym,直到第i次与第i-1次计算出的Ym(i)和Ym(i-1)之间相对误差小于允许误差为止。由于该方法考虑群桩效应的p-y数据,且借助具有弹性约束桩头的群桩模型进行分析,从而使计算结果能客观反映具有弹性约束桩头群桩之间的相互作用与变形特性。