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.     

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.     

Probabilistic tsunami hazard assessment for Messina Strait Area (Sicily, Italy)

The general modular Bayesian procedure is applied to provide a probabilistic tsunami hazard assessment (PTHA) for the Messina Strait Area (MSA), Italy. This is the first study in an Italian area where the potential tsunamigenic events caused by both submarine seismic sources (SSSs) and submarine mass failures (SMFs) are examined in a probabilistic assessment. The SSSs are localized on active faults in MSA as indicated by the instrumental data of the catalogue of the Italian seismicity; the SMFs are spatially identified using their propensity to failure in the Ionian and Tyrrhenian Seas on the basis of mean slope and mean depth, and using marine geology background knowledge. In both cases the associated probability of occurrence is provided. The run-ups were calculated at key sites that are main cities and/or important sites along the Eastern Sicily and the Southern Calabria coasts where tsunami events were recorded in the past. The posterior probability distribution combines the prior probability and the likelihood calculated in the MSA. The prior probability is based on the physical model of the tsunami process, and the likelihood is based on the historical data collected by the historical catalogues, background knowledge, and marine geological information. The posterior SSSs and SMFs tsunami probabilities are comparable and are combined to produce a final probability for a full PTHA in MSA.     

SARETL: an expert system for probabilistic displacement-based dynamic 3-D slope stability analysis and remediation of earthquake triggered landslides

 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. In the 3-D analysis, the critical and total slope widths become two new and important parameters. The probabilistic models evaluate the probability of failure under seismic loading and consider the different sources of uncertainties involved in the problem, i.e. uncertainties stemming from the discrepancies between laboratory-measured and in situ values of shear strength parameters, randomness of earthquake occurrence, and earthquake-induced acceleration. The models also take into consideration the spatial variabilities and correlations of soil properties. The developed models are incorporated in a computer program, PTDDSSA (probabilistic three-dimensional dynamic slope stability analysis). These developed analysis/design procedures are incorporated within a code named SARETL (stability analysis and remediation of earthquake-triggered landslides) that was developed in this study for stability analysis and remediation of earthquake-triggered landslides. In addition to the dynamic inertia forces; the developed system takes into consideration the local site effects. The code is capable of: 1. Prediction of permanent deformations that result from landslides under seismic loading using both probabilistic and deterministic approaches. 2. The assessment of landslide hazard that affects major transportation routes in the event of earthquakes, and the preparation of earthquake-induced landslide hazard maps (i.e. maps that show expected displacements and probability of slope/embankments failure) for different earthquakes magnitudes and environmental conditions. 3. Proposing a mitigation strategy against landslides and suggesting guidelines for remedial measures. The developed expert system is applied to a major highway case study. Design maps are developed for the highway under seismic loading. Received: 18 March 1999 · Accepted: 11 October 1999     

Estimation of the Variability of Tailings Dams Properties in Order to Perform Probabilistic Assessment

Dams are often used for storing mine tailings. These structures present failure mechanisms that can lead to major risks for surrounding populations. In view to managing these risks, it is necessary to take a probabilistic approach in predicting their behaviour. The application of such approaches is limited by the difficulty of obtaining experimental data to estimate the variability of the parameters and conditioned by the relevance of the probability models chosen to represent this variability. This article proposes (1) a probabilistic modelling of the index properties of the mine tailings constituting these dams based on statistical analyses and (2) a method using dynamic penetration tests to estimate on site the mine tailings friction angle and its variability. This method, applied to chilean tailings dams proposes a single model for all tailings dams in order to associate a probability law to the effective friction angle (ϕ′). The procedure is illustrated on the probabilistic study of slope stability carried out at the global scale of a dam and also at the local scale of each of its constituent layers.     

小坡度海底土层地震液化诱发滑移分析方法

地震可使海底砂质、粉质土层液化并导致上部土层的滑移。基于有效应力有限元动力分析方法和Newmark刚性滑块理论,提出了一种计算海底小坡度（≤5o）土层地震液化引起侧向滑移的简化方法。该方法将波浪荷载简化为海底恒定的上覆压力和初始孔压,忽略了海水粘性对海底土层地震反应的影响,利用改进的Seed孔压模型进行动力分析和液化判别,用Newmark滑块理论计算了土层侧向滑移。通过算例和对比分析,研究了海水深度和土层坡度对侧向滑移的影响,表明该方法的有效性,可为近海工程场地地震地质灾害评价提供参考数据。     

A method for producing digital probabilistic seismic landslide hazard maps

The 1994 Northridge, California, earthquake is the first earthquake for which we have all of the data sets needed to conduct a rigorous regional analysis of seismic slope instability. These data sets include: (1) a comprehensive inventory of triggered landslides, (2) about 200 strong-motion records of the mainshock, (3) 1:24 000-scale geologic mapping of the region, (4) extensive data on engineering properties of geologic units, and (5) high-resolution digital elevation models of the topography. All of these data sets have been digitized and rasterized at 10 m grid spacing using ARC/INFO GIS software on a UNIX computer. Combining these data sets in a dynamic model based on Newmark's permanent-deformation (sliding-block) analysis yields estimates of coseismic landslide displacement in each grid cell from the Northridge earthquake. The modeled displacements are then compared with the digital inventory of landslides triggered by the Northridge earthquake to construct a probability curve relating predicted displacement to probability of failure. This probability function can be applied to predict and map the spatial variability in failure probability in any ground-shaking conditions of interest. We anticipate that this mapping procedure will be used to construct seismic landslide hazard maps that will assist in emergency preparedness planning and in making rational decisions regarding development and construction in areas susceptible to seismic slope failure.     

Stochastic seismic slope stability assessment using polynomial chaos expansions combined with relevance vector machine

This paper presents probabilistic assessment of seismically-induced slope displacements considering uncertainties of seismic ground motions and soil properties.A stochastic ground motion model representing both the temporal and spectral non-stationarity of earthquake shakings and a three-dimensional rotational failure mechanism are integrated to assess Newmark-type slope displacements.A new probabilistic approach that incorporates machine learning in metamodeling technique is proposed,by combining relevance vector machine with polynomial chaos expansions(RVM-PCE).Compared with other PCE methods,the proposed RVM-PCE is shown to be more effective in estimating failure probabilities.The sensitivity and relative influence of each random input parameter to the slope displacements are discussed.Finally,the fragility curves for slope displacements are established for sitespecific soil conditions and earthquake hazard levels.The results indicate that the slope displacement is more sensitive to the intensities and strong shaking durations of seismic ground motions than the frequency contents,and a critical Arias intensity that leads to the maximum annual failure probabilities can be identified by the proposed approach.     

Tsunami hazard assessment of Canada

We present a preliminary probabilistic tsunami hazard assessment of Canadian coastlines from local and far-field, earthquake, and large submarine landslide sources. Analyses involve published historical, palaeotsunami and palaeoseismic data, modelling, and empirical relations between fault area, earthquake magnitude, and tsunami run-up. The cumulative estimated tsunami hazard for potentially damaging run-up (≥1.5 m) of the outer Pacific coastline is ~40–80 % in 50 years, respectively one and two orders of magnitude greater than the outer Atlantic (~1–15 %) and the Arctic (<1 %). For larger run-up with significant damage potential (≥3 m), Pacific hazard is ~10–30 % in 50 years, again much larger than both the Atlantic (~1–5 %) and Arctic (<1 %). For outer Pacific coastlines, the ≥1.5 m run-up hazard is dominated by far-field subduction zones, but the probability of run-up ≥3 m is highest for local megathrust sources, particularly the Cascadia subduction zone; thrust sources further north are also significant, as illustrated by the 2012 Haida Gwaii event. For Juan de Fuca and Georgia Straits, the Cascadia megathrust dominates the hazard at both levels. Tsunami hazard on the Atlantic coastline is dominated by poorly constrained far-field subduction sources; a lesser hazard is posed by near-field continental slope failures similar to the 1929 Grand Banks event. Tsunami hazard on the Arctic coastline is poorly constrained, but is likely dominated by continental slope failures; a hypothetical earthquake source beneath the Mackenzie delta requires further study. We highlight areas susceptible to locally damaging landslide-generated tsunamis, but do not quantify the hazard.     

Probabilistic Assessment of Earthquake Insurance Rates for Turkey

A probabilistic model is presented to obtain a realistic estimate of earthquake insurance rates for reinforced concrete buildings in Turkey. The model integrates information on seismic hazard and information on expected earthquake damage on engineering facilities in a systematic way, yielding to estimates of earthquake insurance premiums. In order to demonstrate the application of the proposed probabilistic method, earthquake insurance rates are computed for reinforced concrete buildings constructed in five cities located in different seismic zones of Turkey. The resulting rates are compared with the rates currently charged by the insurance companies. The earthquake insurance rates are observed to be sensitive to the assumptions on seismic hazard and damage probability matrices and to increase significantly with increasing violation of the code requirements.     

Development of a probabilistic approach for rock wedge failure

For rock slope engineering, uncertainty and variability are inherent in data collected on orientation and strength of discontinuities, yielding a range of results. Unfortunately, conventional deterministic analysis based on the factor of safety concept, requires a fixed representative value for each parameter without regard to the degree of uncertainty involved. Therefore, the deterministic analysis fails to properly represent uncertainty and variability, so common in engineering geology studies. To overcome this shortcoming, the probabilistic analysis method was proposed and used for more than a decade in rock slope stability analysis. However, most probabilistic analyses included a deterministic model as part of the analysis procedure causing subsequent problems, which went uncorrected. The objectives of this paper are to develop a solution for these difficulties in probabilistic analyses and to propose an appropriate simulation procedure for the probabilistic analysis of rock wedge failures. As part of the solution, probability of kinematic instability and probability of kinetic instability are evaluated separately to provide a proper, combined evaluation for failure probability. To evaluate the feasibility of this new probabilistic approach, the procedure is applied to a practical example, a major, highway rock cut in North Carolina, USA. Results of the probabilistic approach are compared to those of the deterministic analysis; findings are significantly different, indicating that the deterministic analysis does not depict rock slope variations, particularly where significant scatter in parameter data occurs.     

Seismic landslide analysis: Gurpinar (Istanbul) as a case history

Slope failures triggered by earthquakes are among the most important soil mechanics problems. In this study, static and pseudostatic analyses of slope stability for earthquake conditions were carried out in the Gurpinar area. In situ testing (SPT) was carried out and laboratory samples were obtained from six boreholes (maximum depth 50.0 m) to determine soil classification and strength characteristics. Geophysical studies (seismic refraction and MASW) were also carried out in the area to estimate the structure and stiffness strength characteristics of the slope to 50.0 m depth. All field and laboratory data were used to determine the mechanical and structural (geometrical) behavior of the slope. In order to solve the slope stability problem, three soil slope models were considered. Pseudostatic analysis was carried out to estimate the earthquake acceleration seismic hazard in the region. These analyses showed that, while there was no potential slope instability under static load conditions, some problems would appear with increasing earthquake acceleration. A geotechnical slope improvement project is proposed for the study area.     

基于Newmark模型的概率地震滑坡危险性模型参数优化与应用:以鲁甸地震区为例

应用概率地震危险性评价模型进行地震滑坡危险性区划,是解决潜在地震诱发滑坡危险性评价中震源不确定性与诱发滑坡时空不确定性的有效方法 .通过理论分析,结合鲁甸地震区的实际情况,对基于力学原理的Newmark滑块位移模型与概率地震滑坡危险性分析方法中的参数的不确定性问题进行了分析,将斜坡岩土体地震作用下的强度衰减效应、地震加速度地形放大效应、断层破碎带效应融合到了斜坡累积位移计算模型中,进行了模型计算参数的优化.改进后的分析模型,更好地反映了高陡斜坡地形与断层破碎带对地震滑坡灾害发育的控制作用,在鲁甸地震区域滑坡应用中,优化模型中的滑坡失稳极高风险区与实际地震滑坡分布表现出了较好的一致性,在超越概率2%的滑坡失稳概率分布中,鲁甸地区包谷垴-小河断裂、鲁甸-昭通断裂带及牛栏江河谷地带地震滑坡高-极高风险区分布面积增幅十分显著.因此,在Newmark滑块位移模型中考虑地震动参数与岩土参数动态响应规律与变量间的定量关系,对于提高区域斜坡稳定性分析的可靠性具有重要意义.     

Influence of the spatial variability of shear strength parameters on rainfall induced landslides: a case study of sandstone slope in Japan

Rainfall-induced landslides frequently occur in humid temperate regions worldwide. Research activity in understanding the mechanism of rainfall-induced landslides has recently focused on the probability of slope failure involving non-homogeneous soil profiles. This paper presents probabilistic analyses to assess the stability of unsaturated soil slope under rainfall. The influence of the spatial variability of shear strength parameters on the probability of rainfall-induced slope failure is conducted by means of a series of seepage and stability analyses of an infinite slope based on random fields. A case study of shallow failure located on sandstone slopes in Japan is used to verify the analysis framework. The results confirm that a probabilistic analysis can be efficiently used to qualify various locations of failure surface caused by spatial variability of soil shear strength for a shallow infinite slope failure due to rainfall.     

Probabilistic seismic and tsunami hazard analysis for design criteria development of a coastal area in the city of Banda Aceh

Both seismic and tsunami hazards design criteria are essential input to the rehabilitation and long-term development of city of Banda Aceh Post Sumatra 2004 (M _w=9.3) disaster. A case study to develop design criteria for future disaster mitigation of the area is presented. The pilot study consists of probabilistic seismic and tsunami hazard analysis. Results of the probabilistic seismic hazard analysis indicates that peak ground acceleration at baserock for 10 and 2% probability of exceedance in 50 years is 0.3 and 0.55 g, respectively. The analysis also provides spectral values at short (T=0.2 s) and long period (T=1.0 s) motions. Some non-linear time-domain earthquake response analyses for soft, medium, and hard site-class were conducted to recommend design response spectra for each site-class. In addition, tsunami inundation maps generated from probabilistic tsunami hazard analysis were developed through tsunami wave propagation analysis and run-up numerical modeling associated with its probability of tsunamigenic earthquake source potential. Both the seismic and tsunami hazard curve and design criteria are recommended as contribution of this study for design criteria, as part of the disaster mitigation effort in the development process of the city. The methodology developed herein could be applied to other seismic and tsunami disaster potential areas.     

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.     

基于极限分析上限方法的海底斜坡地震稳定性

极限分析上限方法在海底斜坡稳定性评价中受到了广泛关注,但已有成果未考虑地震荷载以及多土层海底斜坡可能出现的局部破坏机制。基于上限定理,推导了拟静力水平地震条件下多土层海底斜坡外力功率与内能耗散率平衡方程;结合强度折减技术和最优化方法,求解了海底斜坡整体和局部地震稳定性安全系数,并实现了多土层海底斜坡的局部滑动面搜索;通过典型算例分析,验证了本文方法的有效性。在此基础上,探讨了不同水平地震条件下两种组合土层海底斜坡的整体和局部稳定性,通过与数值解对比,其结果可以较准确地评价海底斜坡稳定性并有效预测滑移面位置。最后,将极限分析上限方法应用于一海底斜坡工程实例。     

Probabilistic assessment of soil liquefaction considering spatial variability of CPT measurements

The conventional liquefaction potential assessment methods (also known as simplified methods) profoundly rely on empirical correlations based on observations from case histories. A probabilistic framework is developed to incorporate uncertainties in the earthquake ground motion prediction, the cyclic resistance prediction, and the cyclic demand prediction. The results of a probabilistic seismic hazard assessment, site response analyses, and liquefaction potential analyses are convolved to derive a relationship for the annual probability and return period of liquefaction. The random field spatial model is employed to quantify the spatial uncertainty associated with the in-situ measurements of geotechnical material.     

Probabilistic analysis of landslide potential of an inclined uniform soil layer of infinite length: application

An extended probabilistic model that is a modification of the Chen et al. (2007,) model for evaluating the failure probability of an inclined soil layer with an infinite length was developed in the present paper, and then applied to evaluate the occurrence probability of landslide-related debris flow in Tungmen gully located in the eastern Taiwan, which occurred a devastating debris flow in 1990. The statistical properties of hydrogeological parameters were collected and summarized, and then used to evaluate the landslide-related debris-flow probabilities at various relative water depths for Tungmen gully by using the probabilistic model. Under the assumption that the soil is saturated, the soil’s cohesion is negligible and the specific gravity of the solid particles of soils is a constant, a simplified probabilistic critical slope equation for the stability of an infinite slope of soils was also developed, and used to estimate the occurrence probability of debris flow. The result shows that probabilistic landslide analysis for an infinite slope could provide a suitable approximation for the risk analysis of debris flow mobilization at a given gully.     

CST-based first order second moment method for probabilistic slope stability analysis

Correlated random variables are commonly involved in probabilistic slope stability analysis, such as reliability analysis of slopes with spatially variable soil properties. This paper proposes a simple Correlated Sampling Technique (CST) for generating samples of correlated random variables. The CST firstly produces correlated standard-normally distributed samples through linear combinations of independent standard-normally distributed samples. Correlated arbitrarily distributed samples can then be obtained by the Nataf transformation. The CST was combined with FOSM (named CST-based FOSM) for probabilistic slope stability analysis. The slope reliabilities of a single-layered cohesive soil slope and a high earth and rockfill dam were analyzed to illustrate the CST-based FOSM. These illustrative examples indicated that the CST-based FOSM can accurately estimate the slope reliability indices with considerably fewer simulations (especially in the case of low failure probability) compared with direct MCS, and the slope reliability was sensitive to the correlation of the strength parameters.