Permanent Scatterers for landslide investigations: outcomes from the ESA-SLAM project

Within the SLAM project (Service for Landslide Monitoring), launched in 2003 by the European Space Agency (ESA) the Permanent Scatterers (PS) technique, a multi-image interferometric approach, coupled with the interpretation of aerial-photos and optical satellite images, was carried out for landslide investigations. The PS analysis was applied at a regional scale as support for landslide inventory mapping and at local scale for the monitoring of single well-known slope movements. For the integration of the PS measurements within a landslide inventory the Arno river basin (Italy) was chosen as test site for the presence of a high number of mass movements (to date about 300 areas at high landslide risk and more than 27,000 individual landslides mapped by the institutional authorities). About 350 SAR images have been interferometrically processed by means of the PS technique, with the detection of about 600,000 PS. The use of optical images contributed spatial meaning to the point-wise information provided by the PS, making it easier to identify terrain features related to slope instability and the landslide boundaries. Here we describe the employed methodology and its impact in the updating of a preexisting landslide inventory. 6.8% of the total number of landslides were characterized by ground displacement measurements from the PS: 6.1% of already mapped landslides and 0.8% of new unstable areas detected through the PS analysis. Moreover, most of the PS are located in urban areas, showing that the proposed methodology is suitable for landslide mapping in areas with a quite high density of urbanization, but that over vegetated areas it still suffers from the limitations induced by the current space-borne SAR missions (e.g. temporal de-correlation). On the other hand, the use of InSAR for the monitoring of single slow landslides threatening built-up areas has provided satisfactory results, allowing the measurement of superficial deformations with high accuracy on the landslide sectors characterized by a good radar reflectivity and coherence.     

Ground deformation monitoring by using the Permanent Scatterers Technique: The example of the Oltrepo Pavese (Lombardia, Italy)

The applicability of the Permanent Scatterers Synthetic Aperture Radar Interferometry (PSInSAR) technique for detecting and monitoring ground displacements was tested in the Oltrepo Pavese territory (Northern Italy, southern Lombardia), which could be representative of similar geological contexts in the Italian Apennines. The study area, which extends for almost 1100 km², is characterized by a complex geological and structural setting and the presence of clay-rich sedimentary formations. These characteristics make the Oltrepo Pavese particularly prone to several geological hazards: shallow and deep landslides, subsidence and swelling/shrinkage of the clayey soils. The PSInSAR technique used in this study overcomes most of the limitations of conventional interferometric approaches by identifying, within the area of interest, a set of “radar benchmarks” (PS), where very precise displacement measurements can be carried out. More than 90,000 PS were identified by processing Synthetic Aperture Radar (SAR) images acquired from 1992 to 2001 by the European Remote Sensing satellites (ERS). The PSInSAR application at a sub-regional scale detected slow ground deformations ranging from + 5 to − 16 mm/year, and resulting from various processes (landslides, swelling/shrinkage of clay soils and water pumping). The PS displacements were analysed by collecting data obtained through geological, geomorphologic field surveys, geotechnical analysis of the soils and the information was integrated within a landslide inventory and the damaged building inventory. Despite the limited number of landslide bodies with PS (7% of the inventoried landslides), the PS data helped to revise the state of activity of several landslides. Furthermore, some previously unknown unstable slopes were detected. Two areas of uplift and two areas of subsidence were identified.     

Statistical modelling of Europe-wide landslide susceptibility using limited landslide inventory data

In many regions, the absence of a landslide inventory hampers the production of susceptibility or hazard maps. Therefore, a method combining a procedure for sampling of landslide-affected and landslide-free grid cells from a limited landslide inventory and logistic regression modelling was tested for susceptibility mapping of slide- and flow-type landslides on a European scale. Landslide inventories were available for Norway, Campania (Italy), and the Barcelonnette Basin (France), and from each inventory, a random subsample was extracted. In addition, a landslide dataset was produced from the analysis of Google Earth images in combination with the extraction of landslide locations reported in scientific publications. Attention was paid to have a representative distribution of landslides over Europe. In total, the landslide-affected sample contained 1,340 landslides. Then a procedure to select landslide-free grid cells was designed taking account of the incompleteness of the landslide inventory and the high proportion of flat areas in Europe. Using stepwise logistic regression, a model including slope gradient, standard deviation of slope gradient, lithology, soil, and land cover type was calibrated. The classified susceptibility map produced from the model was then validated by visual comparison with national landslide inventory or susceptibility maps available from literature. A quantitative validation was only possible for Norway, Spain, and two regions in Italy. The first results are promising and suggest that, with regard to preparedness for and response to landslide disasters, the method can be used for urgently required landslide susceptibility mapping in regions where currently only sparse landslide inventory data are available.     

Investigating landslides with space-borne Synthetic Aperture Radar (SAR) interferometry

This paper is addressed to readers without advanced knowledge of remote sensing. It illustrates some current and potential uses of satellite Synthetic Aperture Radar interferometry (InSAR) for landslide assessment. Data acquired by SAR systems can provide 3D terrain models and be used to assist in regional scale investigations, e.g. aimed at evaluation of susceptibility of slopes to failure. Under favourable environmental conditions, the innovative Permanent Scatterers (PS) technique, which overcomes several limitations of conventional SAR differential interferometry (DInSAR) applications in landslide studies, is suitable for monitoring slope deformations with millimetric precision. The PS technique combines the wide-area coverage typical of satellite imagery with the capability of providing displacement data relative to individual image pixels. With the currently available radar satellites, however, only very slow ground surface displacements can be reliably detected and measured. The presented case study of a landslide from the Liechtenstein Alps indicates that the most attractive and reliable contribution provided by this remote sensing technique lies in the possibility of (i.) wide-area qualitative distinction between stable and unstable areas and (ii.) qualitative (relative) hazard zonation of large, slow landslides based on the identification of segments characterised by different movement rates. Since only the radar line of sight projection of the displacements can be detected, a quantitative exploitation of the PS data is possible only where sufficient ground truth is available. In site specific or single landslide investigations the PS data can represent a very useful complementary data source with respect to the information acquired through ground based observations and in situ surveying. However, the difficulties associated with the feasibility assessments of the applicability of SAR data to local scale problems, as well as with the interpretation of PS results, require a close collaboration between landslide experts and specialists in advanced processing of radar satellite data. The interpretation of the exact geotechnical significance of small, radar sensed ground surface deformations is challenging, especially where ground truth is lacking. Although any ground deformation is potentially of interest to an engineering geologist, detection of movements in both vertical and horizontal directions is needed in the case of landslides to evaluate slope failure mechanisms. With their high radar viewing angles, however, the current space-borne systems can detect only a fraction of the horizontal component of movement. It is expected that the upcoming SAR dedicated missions with new sensors and different acquisition geometries, combined with the rapid developments in the field of advanced radar data processing, will allow a full 3D reconstruction of deformation data and help to further reduce the current limitations of the PS and similar DInSAR approaches.     

Activity and kinematic behaviour of deep-seated landslides from PS-InSAR displacement rate measurements

Large landslides and deep-seated gravitational slope deformations (DSGSD) represent an important geo-hazard in relation to the deformation of large structures and infrastructures and to the associated secondary landslides. DSGSD movements, although slow (from a few millimetres to several centimetres per year), can continue for very long periods, producing large cumulative displacements and undergoing partial or complete reactivation. Therefore, it is important to map the activity of such phenomena at a regional scale. Ground surface displacements at DSGSD typically range close to the detection limit of monitoring equipment but are suitable for synthetic aperture radar (SAR) interferometry. In this paper, permanent scatterers (PSInSAR?) and SqueeSAR? techniques are used to analyse the activity of 133 DSGSD, in the Central Italian Alps. Statistical indicators for assigning a degree of activity to slope movements from displacement rates are discussed together with methods for analysing the movement and activity distribution within each landslide. In order to assess if a landslide is active or not, with a certain degree of reliability, three indicators are considered as optimal: the mean displacement rate, the activity index (ratio of active PS, displacement rate larger than standard deviation, overall PS) and the nearest neighbor ratio, which allows to describe the degree of clustering of the PS data. According to these criteria, 66% of the phenomena are classified as active in the monitored period 1992–2009. Finally, a new methodology for the use of SAR interferometry data to attain a classification of landslide kinematic behaviour is presented. This methodology is based on the interpretation of longitudinal ground surface displacement rate profiles in the light of numerical simulations of simplified failure geometries. The most common kinematic behaviour is rotational, amounting to 41 DSGSDs, corresponding to the 62.1% of the active phenomena.     

The new landslide inventory of Tuscany (Italy) updated with PS-InSAR: geomorphological features and landslide distribution

In this paper, the updating of the landslide inventory of Tuscany region is presented. To achieve this goal, satellite SAR data processed with persistent scatter interferometry (PSI) technique have been used. The updating leads to a consistent reduction of unclassified landslides and to an increasing of active landslides. After the updating, we explored the characteristics of the new inventory, analysing landslide distribution and geomorphological features. Several maps have been elaborated, as sliding index or landslide density map; we also propose a density-area map to highlight areas with different landslide densities and sizes. A frequency-area analysis has been performed, highlighting a classical negative power-law distribution. We also explored landslide frequency for lithology, soil use and several morphological attributes (elevation, slope gradient, slope curvature), considering both all landslides and classified landslide types (flows, falls and slides).     

Landslide hazard and risk mapping at catchment scale in the Arno River basin

We present the methodologies adopted and the outcomes obtained in the analysis of landslide risk in the basin of the Arno River (Central Italy) in the framework of a project sponsored by the Basin Authority of the Arno River, started in the year 2002 and completed at the beginning of 2005. In particular, a complete set of methods and applications for the assessment of landslide susceptibility and risk are described and discussed. A new landslide inventory of the whole area was realized, using conventional (aerial-photo interpretation and field surveys) and non-conventional methods (e.g. remote sensing techniques such as DInSAR and PS-InSAR). The great majority of the mapped mass movements are rotational slides (75%), solifluctions and other shallow slow movements (17%) and flows (5%), while soil slips, and other rapid landslides, seem less frequent everywhere within the basin. The relationships between landslide characteristics and environmental factors have been assessed through statistical analysis. As expected, the results show a strong control of land cover, lithology and morphology on landslide occurrence. The landslide frequency-size distribution shows a typical scaling behaviour already underlined in other landslide inventories worldwide. The assessment of landslide hazard in terms of probability of occurrence in a given time, based for mapped landslides on direct and indirect observations of the state of activity and recurrence time, has been extended to landslide-free areas through the application of statistical methods implemented in an artificial neural network (ANN). Unique conditions units (UCU) were defined by the map overlay of landslide preparatory factors (lithology, land cover, slope gradient, slope curvature and upslope contributing area) and afterwards used to construct a series of model vectors for the training and test of the ANN. Various different ANNs were selected throughout the basin, until each UCU was assigned a degree of membership to a susceptibility and a hazard class. Model validation confirms that prediction results are very good, with an average percentage of correctly recognized mass movements of about 85%. The analysis also revealed the existence of a large number of unmapped mass movements, thus contributing to the completeness of the final inventory. Temporal hazard was estimated via the translation of state of activity in recurrence time and hence probability of occurrence. The following intersection of hazard values with vulnerability and exposure figures, obtained by reclassification of digital vector mapping at 1:10,000 scale, lead to the definition of risk values for each terrain unit for different periods of time into the future. The final results of the research are now undergoing a process of integration and implementation within land planning and risk prevention policies and practices at local and national level.     

时序InSAR技术三峡库区藕塘滑坡稳定性监测与状态更新

坡体表面形变是表征坡体稳定性的重要信息,因此,非常有必要对滑坡多发区域进行时序常规变形监测.近年来,星载合成孔径雷达数据由于其覆盖范围大、形变监测精度高的特点,被越来越多的用于山区滑坡识别与探测.首先介绍了联合分布式目标与点目标的时序InSAR方法,并将该方法应用于分析覆盖三峡藕塘滑坡的2007年至2011年的19景ALOS PALSAR数据和2015年至2018年的47景Sentinel-1数据,提取了数据覆盖时间段内的藕塘地区的变形速率.发现相比于2007年至2011年,2015年至2018年新增三处不稳定斜坡.进一步对滑坡的时序变形分析表明,降雨和水位变化是坡体稳定性最大的两个影响因素.实验证明时序InSAR方法可以作为常规形变手段来识别与监测三峡库区等地区潜在的滑坡,为防灾减灾提供支持与依据.      

Mapping of slow landslides on the Palos Verdes Peninsula using the California landslide inventory and persistent scatterer interferometry

Extremely slow landslides, those with a displacement rate <16 mm/year, may be imperceptible without proper instrumentation. These landslides can cause infrastructure damage on a long-term timescale. The objective is to identify these landslides through the combination of information from the California landslide inventory (CLI) and ground displacement rates using results from persistent scatterer interferometry (PSI), an interferometric synthetic aperture radar (InSAR) stacking technique, across the Palos Verdes Peninsula in California. A total of 34 ENVISAT radar images (acquired between 2005 and 2010) and 40 COSMO-SkyMed radar images (acquired between 2012 and 2014) were processed. An InSAR landslide inventory (ILI) is created using four criteria: minimum PS count, average measured ground velocity, slope angle, and slope aspect. The ILI is divided into four categories: long-term slides (LTSs), potentially active slides (PASs), relatively stable slopes (RSSs), and unmapped extremely slow slides (UESSs). These categories are based on whether landslides were previously mapped on that slope (in the CLI), if persistent scatterers (PSs) are present, and whether PSs are unstable or stable. The final inventory includes 263 mapped landslides across the peninsula, of them 67 landslides were identified as UESS. Although UESS exhibit low velocity and are relatively small (average area of 8865 m² per slide), their presence in a highly populated area such as the Palos Verdes Peninsula could lead to destruction of infrastructure and property over the long term.     

Effectiveness of visually analyzing LiDAR DTM derivatives for earth and debris slide inventory mapping for statistical susceptibility modeling

Landslide inventories are the most important data source for landslide process, susceptibility, hazard, and risk analyses. The objective of this study was to identify an effective method for mapping a landslide inventory for a large study area (19,186 km²) from Light Detection and Ranging (LiDAR) digital terrain model (DTM) derivatives. This inventory should in particular be optimized for statistical susceptibility modeling of earth and debris slides. We compared the mapping of a representative set of landslide bodies with polygons (earth and debris slides, earth flows, complex landslides, and areas with slides) and a substantially complete set of earth and debris slide main scarps with points by visual interpretation of LiDAR DTM derivatives. The effectiveness of the two mapping methods was estimated by evaluating the requirements on an inventory used for statistical susceptibility modeling and their fulfillment by our mapped inventories. The resulting landslide inventories improved the knowledge on landslide events in the study area and outlined the heterogeneity of the study area with respect to landslide susceptibility. The obtained effectiveness estimate demonstrated that none of our mapped inventories are perfect for statistical landslide susceptibility modeling. However, opposed to mapping polygons, mapping earth and debris slides with a point in the main scarp were most effective for statistical susceptibility modeling within large study areas. Therefore, earth and debris slides were mapped with points in the main scarp in entire Lower Austria. The advantages, drawbacks, and effectiveness of landslide mapping on the basis of LiDAR DTM derivatives compared to other imagery and techniques were discussed.     

Monitoring of post-failure landslide deformation by the PS-InSAR technique at Lubietova in Central Slovakia

Interferometric synthetic aperture radar data from ERS and ENVISAT sensors were utilized in the analysis of the post-failure deformations in the area of Lubietova town in Central Slovakia. The catastrophic landslide of 1977 together with surrounding landslides in the Lubietova area were analysed with the help of persistent scatterers (PS) technique in order to evaluate recent and past deformations of the unstable slopes. Although long-term precise geodetic monitoring of the 1977 landslide revealed differential deformations inside the sliding mass, due to the lack of the PS located inside the landside caused by temporal decorrelation, unfortunately, these records could not be directly compared. The adjacent landslides with sufficient number of PS were analysed by transformation of the line of sight displacements recorded by the sensors to the slope vector direction. This procedure allowed identification of the precise boundaries of the actively moving landslide parts and the updating of the landslide inventory for the Lubietova area.     

Regional-scale landslide activity and landslide susceptibility zonation in the Nepal Himalaya

Landslide susceptibility zonation mapping is a fundamental procedure for geo-disaster management in tropical and sub-tropical regions. Recently, various landslide susceptibility zonation models have been introduced in Nepal with diverse approaches of assessment. However, validation is still a problem. Additionally, the role of various predisposing causative parameters for landslide activity is still not well understood in the Nepal Himalaya. To address these issues of susceptibility zonation and landslide activity, about 4,000 km² area of central Nepal was selected for regional-scale assessment of landslide activity and susceptibility zonation mapping. In total, 655 new landslides and 9,229 old landslides were identified with the study area with the help of satellite images, aerial photographs, field data and available reports. The old landslide inventory was “blind landslide database” and could not explain the particular rainfall event responsible for the particular landslide. But considering size of the landslide, blind landslide inventory was reclassified into two databases: short-duration high-intensity rainfall-induced landslide inventory and long-duration low-intensity rainfall-induced landslide inventory. These landslide inventory maps were considered as proxy maps of multiple rainfall event-based landslide inventories. Similarly, all 9,884 landslides were considered for the activity assessment of predisposing causative parameters. For the Nepal Himalaya, slope, slope aspect, geology and road construction activity (anthropogenic cause) were identified as most affective predisposing causative parameters for landslide activity. For susceptibility zonation, multivariate approach was considered and two proxy rainfall event-based landslide databases were used for the logistic regression modelling, while a relatively recent landslide database was used in validation. Two event-based susceptibility zonation maps were merged and rectified to prepare the final susceptibility zonation map and its prediction rate was found to be more than 82 %. From this work, it is concluded that rectification of susceptibility zonation map is very appropriate and reliable. The results of this research contribute to a significant improvement in landslide inventory preparation procedure, susceptibility zonation mapping approaches as well as role of various predisposing causative parameters for the landslide activity.     

Landslide susceptibility mapping using rough sets and back-propagation neural networks in the Three Gorges, China

In the Three Gorges of China, there are frequent landslides, and the potential risk of landslides is tremendous. An efficient and accurate method of generating landslide susceptibility maps is very important to mitigate the loss of lives and properties caused by these landslides. This paper presents landslide susceptibility mapping on the Zigui-Badong of the Three Gorges, using rough sets and back-propagation neural networks (BPNNs). Landslide locations were obtained from a landslide inventory map, supported by field surveys. Twenty-two landslide-related factors were extracted from the 1:10,000-scale topographic maps, 1:50,000-scale geological maps, Landsat ETM + satellite images with a spatial resolution of 28.5 m, and HJ-A satellite images with a spatial resolution of 30 m. Twelve key environmental factors were selected as independent variables using the rough set and correlation coefficient analysis, including elevation, slope, profile curvature, catchment aspect, catchment height, distance from drainage, engineering rock group, distance from faults, slope structure, land cover, topographic wetness index, and normalized difference vegetation index. The initial, three-layered, and four-layered BPNN were trained and then used to map landslide susceptibility, respectively. To evaluate the models, the susceptibility maps were validated by comparing with the existing landslide locations according to the area under the curve. The four-layered BPNN outperforms the other two models with the best accuracy of 91.53 %. Approximately 91.37 % of landslides were classified as high and very high landslide-prone areas. The validation results show sufficient agreement between the obtained susceptibility maps and the existing landslide locations.     

Landslides in the North of Lisbon Region (Portugal): Conditioning and triggering factors

Detailed geomorphological mapping carried out in 5 sample areas in the North of Lisbon Region allowed us to collect a set of geological and geomorphological data and to correlate them with the spatial occurrence of landslide. A total of 597 slope movements were identified in a total area of 61.7 km², which represents about 10 landslides per km².The main landslide conditioning factors are: lithology and geological structure, slope angle and slope morphology, land use, presence of old landslides, and human activity.The highest landslide density occurs in Cretaceous marls and marly limestones, but the largest movements are in Jurassic clays, marls and limestones.The landslide density is higher on slopes with gradients above 20 °, but the largest unstable area is found on slopes of 10 ° to 15 °, thus reflecting the presence of the biggest slope movements. There is a correlation between landslides and topographical concavities, a fact that can be interpreted as reflecting the significance of the hydrological regime in slope instability.Concerning land use, the highest density of landslides is found on slopes covered with shrub and undergrowth vegetation.About 26% of the total number of landslides are reactivation events. The presence of old landslides is particularly important in the occurrence of translational slides and complex and composite slope movements.20% of the landslide events were conditioned by anthropomorphic activity. Human's intervention manifests itself in ill-consolidated fills, cuts in potentially unstable slopes and, in a few cases, in the changing of river channels.Most slope movements in the study area exhibit a clear climatic signal. The analysis of rainfall distribution in periods of recognised slope instability allows the distinction of three situations: 1) moderate intensity rainfall episodes, responsible for minor slope movements on the bank of rivers and shallow translational slides, particularly in artificial trenches; 2) high intensity rainfall episodes, originating flash floods and most landslides triggered by bank erosion; 3) long-lasting rainfall periods, responsible for the rise of the groundwater table and triggering of landslides with deeper slip surfaces.     

Landslide susceptibility mapping by using fuzzy logic: a case study of Cham-gardalan catchment,Ilam, Iran

Landslides are introduced as regional movements, which influence different engineering structures such as roads, railways, and dams and cause the person’s death. Identification of landslide zones may decrease the financial losses and human injuries or deaths. This study tries to achieve a landslide susceptibility mapping in Cham-gardalan catchment by weighting the main criteria and the membership functions of fuzzy logic. For this, we applied the best relationship function between the presence and absence of landslides as well as a collection of the elements. At first, the landslide points were identified by the means of some components those of satellite images, topographical (1:50,000) and geographical (1:100,000) maps, field visits, and Google Earth software followed by the preparation of landslide distribution maps. Then, all effective landslide factors such as percentage of slope, slope aspect, height, geology, land uses, distance from roads, distance from drainages, distance from breakage, and precipitation map have been utilized in order to conduct the fuzzy analyses. Landslide susceptibility map was performed by fuzzy operators (Gamma, Product, Sum, Or, And) in the study area. After fuzzificating and weighting, the effective criteria of landslides were determined through fuzzy Gamma operators with the landaus of 0.2, 0.5, 0.8, and 0.9 and by comparing final maps for making an appropriate model of landslide susceptibility mapping. The regional susceptibility map represents the landslide-prone areas in five categories those of very low, low, moderate, high, and very high. Our results indicated that among the applied operators, Gamma with landau of 0.9 can be used as an appropriate method for mapping the landslide susceptibility due to the suitable fuzzification of given criteria based on landslide distribution maps. In addition, the elements of road, percentage of slope, distance from drainage, and geology were recognized as the most important factors for occurring the landslides.     

Application of an incomplete landslide inventory, logistic regression model and its validation for landslide susceptibility mapping related to the May 12, 2008 Wenchuan earthquake of China

The main purpose of this paper is to present the use of multi-resource remote sensing data, an incomplete landslide inventory, GIS technique and logistic regression model for landslide susceptibility mapping related to the May 12, 2008 Wenchuan earthquake of China. Landslide location polygons were delineated from visual interpretation of aerial photographs, satellite images in high resolutions, and verified by selecting field investigations. Eight factors, including slope angle, slope aspect, elevation, distance from drainages, distance from roads, distance from main faults, seismic intensity and lithology were selected as controlling factors for earthquake-triggered landslide susceptibility mapping. Qualitative susceptibility analyses were carried out using the map overlaying techniques in GIS platform. The validation result showed a success rate of 82.751 % between the susceptibility probability index map and the location of the initial landslide inventory. The predictive rate of 86.930 % was obtained by comparing the additional landslide polygons and the landslide susceptibility probability index map. Both the success rate and the predictive rate show sufficient agreement between the landslide susceptibility map and the existing landslide data, and good predictive power for spatial prediction of the earthquake-triggered landslides.     

基于SBAS-InSAR技术的深切割高山峡谷区滑坡灾害早期识别

近年来,高山峡谷区滑坡灾害频频发生,给人民生命和财产安全带来严重威胁。针对多数学者利用SAR单轨道数据对高山峡谷区滑坡进行早期识别,存在SAR成像几何畸变造成部分滑坡不能识别、识别结果不全面等问题。为全面准确的对高山峡谷区滑坡隐患进行早期识别,文章采用SBAS-InSAR技术,以东川小江沿线两侧深切割高山峡谷区为研究区,通过升降轨SAR数据结合互补的方式进行滑坡灾害隐患识别,引入高分辨率光学影像等作为辅助识别,最终共识别出18处滑坡灾害体,其中5处为高风险潜在滑坡,并对三类典型潜在滑坡进行分析。分析结果表明:利用升降轨SAR数据结合互补的方式,能有效避免SAR单轨道数据在高山峡谷地区产生的几何畸变问题,同时,该方法能更为准确全面地对高山峡谷区滑坡隐患进行早期识别,为防灾减灾事业及政府部门决策提供一种有效的手段。     

Long-term measurements using an integrated monitoring network to identify homogeneous landslide sectors in a complex geo-environmental context (Lago,Calabria, Italy)

Mass movements of different thicknesses and types affect the slope where the Greci and Piscopie villages lie (Lago, Calabria, Italy). The study concerns a moderately urbanized area where the lithological, structural, and hydrogeological conditions have produced a large and complex unstable area. In the study area, a multi-temporal field survey has allowed for detailed landslide inventory mapping and definition of the geometric and kinematic characteristics of mass movements with long-term monitoring. In particular, an integrated monitoring network was progressively put in place, since 1996 to measure superficial and deep displacements (GPS stations and inclinometers). The network dataset shows that the mass movements in the study area are characterized by an extremely slow or very slow rate of displacement. The elements acquired by a long-term monitoring of the deep displacements point out that the landslide bodies are prone to develop new failure surfaces progressively shallower, with a consequent increase of the hazard. Moreover, medium-deep and deep-seated landslides may rapidly accelerate in relation to rainfall conditions and significantly affect structures and infrastructures. The integration of the geological and geomorphological knowledges with the monitoring data allows for distinguishing six different homogeneous sectors on the slope, as well as the relative geometric and kinematic characteristics, and the type of mass movements. In the study area, which is representative of several aspects of other areas, the method proposed and the knowledge acquired by long-term monitoring could be useful to define mass movement mechanisms, geotechnical models, and risk mitigation strategies.     

Use of GIS-based fuzzy logic relations and its cross application to produce landslide susceptibility maps in three test areas in Malaysia

Landslides are one of the most frequent and common natural hazards in Malaysia. Preparation of landslide susceptibility maps is one of the first and most important steps in the landslide hazard mitigation. However, due to complex nature of landslides, producing a reliable susceptibility map is not easy. For this reason, a number of different approaches have been used, including direct and indirect heuristic approaches, deterministic, probabilistic, statistical, and data mining approaches. Moreover, these landslides can be systematically assessed and mapped through a traditional mapping framework using geoinformation technologies. Since the early 1990s, several mathematical models have been developed and applied to landslide hazard mapping using geographic information system (GIS). Among various approaches, fuzzy logic relation for mapping landslide susceptibility is one of the techniques that allows to describe the role of each predisposing factor (landslide-conditioning parameters) and their optimal combination. This paper presents a new attempt at landslide susceptibility mapping using fuzzy logic relations and their cross application of membership values to three study areas in Malaysia using a GIS. The possibility of capturing the judgment and the modeling of conditioning factors are the main advantages of using fuzzy logic. These models are capable to capture the conditioning factors directly affecting the landslides and also the inter-relationship among them. In the first stage of the study, a landslide inventory was complied for each of the three study areas using both field surveys and airphoto studies. Using total 12 topographic and lithological variables, landslide susceptibility models were developed using the fuzzy logic approach. Then the landslide inventory and the parameter maps were analyzed together using the fuzzy relations and the landslide susceptibility maps produced. Finally, the prediction performance of the susceptibility maps was checked by considering field-verified landslide locations in the studied areas. Further, the susceptibility maps were validated using the receiver-operating characteristics (ROC) success rate curves. The ROC curve technique is based on plotting model sensitivity—true positive fraction values calculated for different threshold values versus model specificity—true negative fraction values on a graph. The ROC curves were calculated for the landslide susceptibility maps obtained from the application and cross application of fuzzy logic relations. Qualitatively, the produced landslide susceptibility maps showed greater than 82% landslide susceptibility in all nine cases. The results indicated that, when compared with the landslide susceptibility maps, the landslides identified in the study areas were found to be located in the very high and high susceptibility zones. This shows that as far as the performance of the fuzzy logic relation approach is concerned, the results appeared to be quite satisfactory, the zones determined on the map being zones of relative susceptibility.     

2008年汶川地震滑坡详细编目及其空间分布规律分析

最新研究成果表明, 2008年5月12日汶川MS 8.0级地震触发了超过197000处滑坡。首先,基于GIS与遥感技术构建了汶川地震滑坡的3类编目图,分别为单体滑坡面分布数据、滑坡中心点位置和滑坡后壁点位置。构建方法为基于地震前后高分辨率遥感影像的目视解译方法,区分单体滑坡并圈定其边界,对滑坡后壁进行识别与定点,并开展了部分滑坡的野外验证工作。这些滑坡分布在一个面积大约为110000km2的区域内,滑坡总面积约为1160km2。选择一个面积约为44031km2的区域作为研究区,区内滑坡数量为196007个,滑坡面积为1150.622km2,这是最详细完整的汶川地震滑坡编录成果,也是单次地震事件触发滑坡最多的记录。其次,开展研究区内的地震滑坡空间分布规律的研究。基于滑坡面与滑坡中心点分别构建滑坡空间分布面积密度图与点密度图,结果表明:滑坡多沿着映秀北川断裂分布,多发生在断裂的上盘。滑坡的高密度区位于映秀北川同震地表破裂的南西段(映秀镇与北川县之间)的上盘区域,这一区域恰对应着逆冲分量为主的断裂上盘,表明逆冲断裂对上盘区域发生滑坡的极强烈的控制作用,而该区域正是形变最大的区域,因此说明是地震滑坡发生的强烈控制作用。基于滑坡面密度(LAP)、滑坡中心点密度(LCND)与滑坡后壁点密度(LTND)这3个衡量指标,使用统计分析方法,评价了汶川地震滑坡与地震参数、地质参数、地形参数的关系。结果表明:LAP、LCND与LTND这3个衡量指标与坡度、地震烈度与PGA存在明显的正相关关系; 与距离震中、距离映秀北川同震地表破裂存在负相关关系; 斜坡曲率越接近0,滑坡越不易发生; LAP、LCND与LTND的高值高程区间为1200～3000m; 滑坡发生的优势坡向为E、SE、S方向; 滑坡发育的易发岩性为砂岩与粉砂岩(Z)、花岗岩; 滑坡与坡位的相关关系不太明显。统计结果还表明LCND与LTND两个衡量指标的差异对地震与地质因子不敏感,而对地形因子较敏感。最后将本文的统计结果与以往的汶川地震滑坡空间分布规律统计成果进行了一些对比,对比结果表明,对于某些因子,如高程、岩性、距离震中、距离映秀北川断裂的统计分析结果,采用不完整的滑坡分布数据或点数据,与采用较完整的滑坡分布面数据会有一定的差异,这种差异并未出现在针对坡度与坡向等因子的统计对比结果中。总之,作者认为一个完备、详细的地震滑坡分布面要素编目图是地震滑坡空间分布规律定量分析、危险性定量分析与滑坡控制的地震区地貌演化研究的重要基础,否则,与实际情况相比,得到统计结果会有一定的偏差,本文的研究成果与以往成果的对比结果证明了这一点。