A Reliability Based Expert System for Assessment and Mitigation of Landslides Hazard Under Seismic Loading

Different models were developed for evaluating the probabilistic three-dimensional (3-D) stability analysis of earth slopes and embankments under earthquake loading using both the safety factor and the displacement criteria of slope failure.The probabilistic models evaluate the probability of failure under seismic loading considering the different sources of uncertainties involved in the problem. The models also take into consideration the spatial variabilities and correlations of soil properties. The developed models are incorporated in a computer program PTDDSSA.These analysis/design procedures are incorporated within a code named SARETL developed in this study for stability analysis and remediation of earthquake triggered landslides. In addition to the dynamic inertia forces, the system takes into consideration local site effects.The code is capable of assessing the landslide hazard affecting major transportation routes in the event of earthquakes and preparing earthquake induced landslide hazard maps (i.e., maps showing expected displacements and probability of slope/embankments failure) for different earthquake magnitudes and environmental conditions. It can also beused for proposing a mitigation strategy against landslides.     

Seismically induced landslide displacements: a predictive model

Newmark's model for predicting earthquake-induced landslide displacements provides a simple way to predict the coseismic displacements affecting a sliding mass subject to earthquake loading. In this model, seismic slope stability is measured in terms of critical acceleration, which depends on the mechanical soil properties, pore-pressure distribution, and slope geometry. The triggering seismic forces are investigated in terms of energy radiation from the source, propagation, and site effects, based on 190 accelerometric recordings from 17 Italian earthquakes with magnitudes between 4.5 and 6.8. The method is based on the calibration of relations having the general form of an attenuation law that relates the energy of the seismic forces to the dynamic shear resistances of the sliding mass to propagate the expected landslide displacements as an inverse function of the distance from the fault rupture; the amount of displacement computed through these relations provides a criterion to predict the occurrence of slope failures. Finally, maps showing, in a deterministic and a probabilistic way, the potential of seismically induced landslide displacements are displayed as a tool to provide seismic landslide scenarios and earthquake-induced landslide hazard maps, respectively.     

Seismic reliability analysis of earth slopes under short term stability conditions

Different models were developed for evaluating the probabilistic three-dimensional (3-D) stability analysis of earth slopes and embankments under earthquake loading. The 3-D slope stability model assumed is that of a simple cylindrical failure surface. The probabilistic models evaluate the probability of failure under seismic loading considering the randomness of earthquake occurrence, and earthquake induced acceleration and uncertainties stemming from the discrepancies between laboratory-measured and in-situ values of shear strength parameters. The models also takes into consideration the spatial variabilities and correlations of soil properties. The probabilistic analysis and design approach is capable of obtaining the 2-D and 3-D static and dynamic safety factors, the probability of slope failure, the earthquake induced acceleration coefficient, the yield acceleration coefficient, the earthquake induced displacement, and the probability of allowable displacement exceedance taking into account the local site effect. The approach is applied to a well known landslide case: Congress Street Landslide in Chicago. A sensitivity analysis was conducted on the different parameters involved in the models by applying those models to the Congress Street landslide considering different levels of seismic hazard. Also, a sensitivity analysis was carried out to study the sensitivity of computed results to input parameters of undrained shear strength, and corrective factors. A comparison was made between the different models of failure. The parametric study revealed that the hypocentral distance and earthquake magnitude have major influence on the earthquake induced displacement, probability of failure and dynamic 2-D and 3-D safety factors.     

地震作用下甘肃天水地区黄土-泥岩接触面滑坡机理

以甘肃天水地区渭河北岸大型黄土-泥岩滑坡为研究对象,通过野外地质调查、工程地质钻探、室内力学试验以及FLAC3D数值计算的方法,对该类滑坡的发育特征及地震作用下的稳定性进行分析评价。结果表明:该类滑坡多属于历史地震滑坡,其破坏类型为滑移拉裂型。通过对滑坡的稳定性进行计算分析可知,静力作用下滑坡处于稳定状态;地震作用下,稳定性明显降低,坡体可能再次滑动,预测滑面位于黄土-泥岩接触带位置,与野外调查观测结果一致。研究结果对天水地区该类滑坡的早期识别有重要参考意义,可为该类滑坡的防灾减灾提供科学依据。     

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.     

Hazard analysis of seismic submarine slope instability

To assess the risk associated with a submarine landslide, one must estimate the probability of slope failure and its consequences. This paper proposes a procedure to estimate the probability of earthquake-induced submarine slope failure (hazard) based on probabilistic seismic hazard analyses, ground response analyses and advanced laboratory tests. The outcomes from these analyses are treated in a probabilistic framework, with analytical simulations using mathematical techniques such as the first-order reliability method, Monte Carlo simulation and Bayesian updating. Fragility curves of slope failure during the earthquake (co-seismic) and after the earthquake (post-seismic) were developed in this study, and were shown to provide a clear and well-organized procedure to estimate the annual failure probability of a submarine slope under earthquake loading.     

GIS技术在区域地震滑坡危险性预测中的应用以龙陵地震滑坡为例

地震滑坡是一种有着严重危害的次生地震灾害形式,形成机制复杂,涉及因素众多。运用G IS丰富的空间分析功能,对地震滑坡的影响因素进行研究,并进行潜在地震滑坡区的预测,是地震滑坡研究领域的一种新的发展趋势。本文在对1976年龙陵地震引发的地震滑坡分布特征研究的基础上,结合前人有关中国西南地区地震滑坡特征的研究成果,应用G IS对该区潜在地震滑坡危险区进行了预测。     

Regional landslide susceptibility following the Mid NIIGATA prefecture earthquake in 2004 with NEWMARK’S sliding block analysis

This study proposes a calculation method for regional earthquake-induced landslide susceptibility that applies the permanent seismic displacement calculated using Newmark’s sliding block analysis with estimated vertical and horizontal seismic motions. The proposed method takes into account the direction of slope failure based on the specified slope azimuth. The study results reveal the importance of predominant slope failure direction using a simple infinite slope model subjected to earthquakes. The target area for the earthquake-induced landslide susceptibility analysis constituted a region of more than 2000 km² surrounding the epicenter of the Mid Niigata prefecture earthquake in 2004. An earthquake-induced landslide susceptibility map was created based on the proposed method with a specific combination of friction angle and cohesion, and the resulting data were compared to the landslide inventory map produced from aerial photographs following the Mid Niigata prefecture earthquake in 2004. To create the susceptibility map, geomaterial cohesion values for the slope were back-calculated to satisfy the minimum safety factor in the static state. This study also proposes a calculation method for the prediction rate and determines the back-calculated strength parameters of geomaterials. The proposed regional landslide susceptibility map will be useful for understanding potential slope failure locations and magnitude of damage, as well as for planning field investigation and preventing secondary disasters immediately after earthquakes.     

A subaqueous hazard map for earthquake-triggered landslides in Lake Zurich,Switzerland

The awareness of geohazards in the subaqueous environment has steadily increased in the past years and there is an increased need to assess these hazards in a quantitative sense. Prime examples are subaqueous landslides, which can be triggered by a number of processes including earthquakes or human activities, and which may impact offshore and onshore infrastructure and communities. In the literature, a plenitude of subaqueous landslide events are related to historical earthquakes, including cases from lakes in Switzerland. Here, we present an approach for a basin-wide earthquake-triggered subaquatic landslide hazard assessment for Lake Zurich, which is surrounded by a densely populated shoreline. Our analysis is based on high-resolution sediment-mechanical and geophysical input data. Slope stabilities are calculated with a grid-based limit equilibrium model on an infinite slope, which uses Monte Carlo sampled input data from a sediment-mechanical stratigraphy of the lateral slopes. Combined with probabilistic ground-shaking forecasts from a recent national seismic hazard analysis, subaquatic earthquake-triggered landslide hazard maps are constructed for different mean return periods, ranging from 475 to 9975 years. Our results provide a first quantitative landslide hazard estimation for the lateral slopes in Lake Zurich. Furthermore, a back-analysis of a case-study site indicates that pseudostatic accelerations in the range between 0.04 and 0.08 g were needed to trigger a well-investigated subaqueous landslide, dated to ~2210 cal. years B.P.     

Landslides induced by the 2017 Mw7.3 Sarpol Zahab earthquake (Iran)

Landslides are the main secondary effects of earthquakes in mountainous areas. The spatial distribution of these landslides is controlled by the local seismic ground motion and the local slope stability. While gravitational instabilities in arid and semi-arid environments are understudied, we document the landslides triggered by the Sarpol Zahab earthquake (November 12, 2017, Mw7.3, Iran/Iraq border), the largest event ever recorded in the semi-arid Zagros Mountains. An original earthquake-induced landslide inventory was derived, encompassing landslides of various sizes and velocities (from rapid disrupted rockfalls to slow-moving coherent landslides). This inventory confirms the low level of triggered landslides in semi-arid environments. It also displays clear differences in the spatial and volumetric distributions of earthquake-induced landslides, having 386 rockfalls of limited size triggered around the epicenter, and 9 giant (areas of ca. 10⁶ m²) active and ancient deep-seated landslides coseismically accelerated at locations up to 180 km from the epicenter. This unusual distant triggering is discussed and interpreted as an interaction between the earthquake source properties and the local geological conditions, emphasizing the key role of seismic ground motion variability at short spatial scales in triggering landslides. Finally, the study documents the kinematics of slow-moving ancient landslides accelerated by earthquakes, and opens up new perspectives for studying landslide triggering over short (~?1–10 years) and long-time (~?1000–10,000 years) periods.

     

环境地球科学之滑坡地震地质学

本文提出滑坡地震地质学学科概念,属于环境地球科学学科范畴。滑坡地震地质学是研究地震荷载下滑坡发生机理、分布与演化规律的学科;也是借助地震滑坡解决地震地质问题的学科。它是地质灾害学、地震地质学、工程地质学之间的一门交叉边缘学科。其也与固体地球物理学、构造地质学、构造地貌学、环境地质学、数值计算、地理信息系统、遥感、统计分析等学科密切相关。文章从科研层面与应用层面的分类大体构建了滑坡地震地质学的研究框架与体系。科研层面的目的是探索地震滑坡规律,应用层面的目的是地震滑坡防灾减灾。根据科研与应用层面的研究目的,将其研究内容大体分为两部分:（1）研究地震滑坡发生机制、分布样式、演化规律等;（2）建立地震滑坡与地震、发震构造、地震破裂过程等的关系式,进而利用地震滑坡解决地震地质问题,如地震问题（地震参数、地震动强度、地震烈度）、发震构造问题（性质、破裂过程）、地质问题（震区地貌演化）等。最后,从地震滑坡数据库,全球范围多震例的地震滑坡分析,地震滑坡与地震地质对象关联的机理研究,地震滑坡演化规律等方面对该学科进行了展望。总之,考虑到随着遥感与GIS等技术的发展成熟、地震滑坡防灾减灾的越来越迫切的需求、多学科的融合交叉的发展趋势等方面,滑坡地震地质学势必会受到越来越广泛的关注。     

CSEES: An expert system for analysis and design of cut slopes and embankments

 An expert system is developed for evaluating failure potential of cut slopes and embankments. The fuzzy sets theory is used with the modified Monte Carlo simulation technique to obtain the Slope Failure Potential Index (SFPI) incorporating factors affecting slope stability, such as geology, topography, geomorphology, precipitation, vegetation and drainage conditions. The developed Cut Slopes and Embankments Expert System (CSEES) includes a classification system for evaluating failure potential of cut slopes and embankments, and a data bank on landslides in Jordan. The proposed classification system and slope failure-potential method proved to be successful for the areas that experienced landslides in the past. The expert system can be used directly for areas with the same geological formations as those areas in which landslides occurred in the past and can be used for areas with other geological formations by modifying the rock type or foundation-material type factor incorporated in the expert system. Received: 18 February 1998 · Accepted: 19 October 1998     

Susceptibility assessments of the Wenchuan earthquake-triggered landslides in Longnan using logistic regression

In this study a Wenchuan earthquake-induced landslide susceptibility assessment was carried out in the Longnan area in northwestern China using a GIS-based logistic regression model. This region has frequently been affected by landslides in the past, and was intensively affected by the 5.12 Wenchuan earthquake which received considerable international attention. The data used for this study consist of the landslides triggered by the Wenchuan earthquake and a landslide pre-disposing factor database. Information regarding the landslide causative factors came from additional data sources, such as a digital elevation model (DEM) with a 30 × 30 m² resolution, orthophotos, geological and land-use maps, precipitation records, and information on peak ground acceleration data from the 2008 earthquake. The statistical analysis of the relationship between the Wenchuan earthquake-triggered landslides and pre-disposing factors showed the great influence of lithological and topographical conditions on slope failures. The quality of susceptibility mapping was validated by splitting the study area into training and validation sections. The prediction capability analysis demonstrated that the landslide susceptibility map could be used for land planning as well as emergency planning by local authorities.     

Analysis of seismic hazard in landslide-prone regions: criteria and example for an area of Daunia (southern Italy)

In relation to the assessment of earthquake-induced landslide hazard, this paper discusses general principles and describes implementation criteria for seismic hazard estimates in landslide-prone regions. These criteria were worked out during the preparation of a hazard map belonging to the official Italian geological cartography and they are proposed as guidelines for future compilation of similar maps. In the presented case study, we used a procedure for the assessment of seismic hazard impact on slope stability adopting Arias intensity Ia as seismic shaking parameter and critical acceleration a _c as parameter representing slope strength to failures induced by seismic shaking. According to this procedure, after a preliminary comparison of estimated historical maximum values of Ia with values proposed in literature as landslide-triggering thresholds, a probabilistic approach, based on the Newmark’s model, is adopted: it allows to estimate the minimum critical acceleration a _c required for a slope to keep under a prefixed value, the probability of failures induced by seismic shakings expected in a given time interval. In this way, one can prepare seismic hazard maps where seismic shaking is expressed in an indirect way through a parameter (the critical acceleration) representing the “strength” that seismic shakings mobilise in slope materials (strength demand) with a prefixed exceedance probability. This approach was applied to an area of Daunia (Apulia—southern Italy) affected by frequent landslide phenomena. The obtained results indicate that shakings with a significant slope destabilisation potential can be expected particularly in the north-western part of the area, which is exposed to the seismic activity of Apennine tectonic structures.     

强震作用下顺层岩质斜坡动力失稳机制及启动速度研究

岩质斜坡地震稳定性评价和地震滑坡启动速度计算是边坡地震失稳防治中的难点。结合5•12四川汶川地震诱发的大型高速滑坡,重点以唐家山高速滑坡为例,提出强震作用下中陡倾顺层岩质斜坡动力失稳机制为拉裂－楔劈－滑移－剪断,并着重阐述了拉裂面形成机制、“楔劈”岩块的杠杆作用和碎屑岩块的滚动摩擦效应。对应于该失稳模式,采用弹性力学理论,推导出地震作用下顺层岩体斜坡锁固段的破坏判据和突发剪断时锁固段岩体的形变能公式。考虑锁固段岩体形变能释放的优势方向,运用能量转化原理,给出了更为合理和精确的顺层岩质滑坡突发启动速度公式,进而确定滑坡启程速度。其研究结果为滑坡灾害范围和冲击损害程度提供有效定量计算方法。     

地震作用下黄土斜坡失稳及运动过程的离散元模拟

在山区,地震诱发滑坡造成的人员伤亡往往占到地震总伤亡人数的一半以上,尤其是在黄土地区,由于地震造成的黄土滑坡具有运动距离远等特征,其灾难性更严重。针对地震诱发黄土斜坡的失稳峰值加速度、临界位移和运动距离等问题,本文利用离散元（PFC）方法,通过对室内三轴实验的应力-应变曲线进行标定和高精度航拍三维地形数据进行转化,以宁夏海口村黄土斜坡为研究对象,开展三维地震作用下斜坡失稳破坏和运动过程的数值模拟研究。通过监测不同位置在地震作用下的应力分量,计算监测点的p、q值并与室内三轴实验所得到p-q破坏线进行比较,获得斜坡破坏过程中的应力路径。结合颗粒的监测位移,获得斜坡失稳破坏的临界位移,得出该黄土斜坡失稳峰值加速度为0.135 g,临界位移为50 cm。同时根据不同地面环境条件,预测了在不同摩擦系数下斜坡失稳后的危害范围,为黄土地区边坡工程的抗震设计及防震减灾工作提供一种新的可视化方法。     

The determination of statistical and deterministic hydrological landslide-triggering thresholds

 Hydrological landslide-triggering thresholds separate combinations of daily and antecedent rainfall or of rainfall intensity and duration that triggered landslides from those that failed to trigger landslides. They are required for the development of landslide early warning systems. When a large data set on rainfall and landslide occurrence is available, hydrological triggering thresholds are determined in a statistical way. When the data on landslide occurrence is limited, deterministic models have to be used. For shallow landslides directly triggered by percolating rainfall, triggering thresholds can be established by means of one-dimensional hydrological models linked to the infinite slope model. In the case of relatively deep landslides located in topographic hollows and triggered by a slow accumulation of water at the soil-bedrock contact, simple correlations between landslide occurrence and rainfall can no longer be established. Therefore real-time failure probabilities have to be determined using hydrological catchment models in combination with the infinite slope model. Received: 15 October 1997 · Accepted: 25 June 1997     

天水市麦积区税湾地震黄土滑坡特征及其形成机制

以天水市税湾地震黄土滑坡为例, 依据野外调查和室内测试结果, 总结天水地区历史地震黄土滑坡特点, 剖析地震黄土滑坡发生的力学机制, 初步提出历史地震黄土滑坡的识别标志。税湾滑坡及柳沟右岸滑坡群属典型的地震黄土滑坡, 具有规模大、滑动面切割深、滑坡坡度小、成群成带分布和高位下滑等特点, 可作为识别历史地震黄土滑坡的重要标志。税湾滑坡及柳沟右岸滑坡群坡体具有明显的黄土/泥岩二元斜坡结构, 极易沿黄土/泥岩接触面滑动。当坡体受到地震力作用时, 地震产生的循环动荷载一方面降低滑坡岩土体的抗剪强度, 另一方面改变滑坡体的力学状态, 坡体应力平衡遭到破坏, 地震力增加坡体下滑力、减小坡体抗滑力, 导致坡体失稳发生滑坡。目前, 税湾滑坡处于欠稳定状态, 遇地震或强降雨有可能再次失稳下滑, 因而有必要进一步开展地震黄土滑坡的成灾模式研究, 为潜在强震区防灾减灾提供科学依据。      

同震滑坡发生概率研究——新一代地震滑坡危险性模型

地震滑坡发生真实概率研究基本空白。本研究创新性的利用贝叶斯概率方法与机器模型开展了中国地震滑坡危险性真实概率研究,制作了第一代中国地震滑坡危险性概率图。基于9个地震案例开展研究,包括1999年台湾集集、2005年克什米尔、2008年汶川、2010年玉树、2013年芦山、2013岷县、2014鲁甸、2015尼泊尔、2017九寨沟地震,这9次地震中7次发生在中国,2005年克什米尔与2015尼泊尔地震均发生在中国邻区,可以更好的控制模型预测精度。这些地震事件均有详细完整的,利用面要素标识的地震滑坡数据,包括306 435处真实的地震滑坡记录。考虑到真实的地震滑坡发生区域,滑坡面积规模的差别,滑坡与不滑样本的比例等因素,共选取了5 117 000个模型训练样本。选择绝对高程、相对高差、坡度、坡向、斜坡曲率、坡位、地形湿度指数、土地覆盖类型、植被覆盖度、与断层距离、地层、年均降水量、地震动峰值加速度共13个地震滑坡影响因子。采用贝叶斯概率方法与机器学习模型相结合,建立地震滑坡发生的多因素影响模型,得到各个连续因子的权重与分类因子的各个分类的权重。再将模型应用到整个中国研究区,地震动峰值加速度因子为触发因子。分别考虑研究区在经历不同地震动峰值加速度（0.1~1 g,每0.1 g一个结果,共10个结果）下的地震滑坡发生真实概率。此外,还结合中国地震动峰值加速度分布图,得到了中国地震动峰值加速度背景下的地震滑坡发生真实概率分布。