Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

Spatial data for landslide susceptibility, hazard, and vulnerability assessment: An overview

The aim of this paper is to discuss a number of issues related to the use of spatial information for landslide susceptibility, hazard, and vulnerability assessment. The paper centers around the types of spatial data needed for each of these components, and the methods for obtaining them. A number of concepts are illustrated using an extensive spatial data set for the city of Tegucigalpa in Honduras. The paper intends to supplement the information given in the “Guidelines for Landslide Susceptibility, Hazard and Risk Zoning for Land Use Planning” by the Joint ISSMGE, ISRM and IAEG Technical Committee on Landslides and Engineered Slopes (JTC-1). The last few decades have shown a very fast development in the application of digital tools such as Geographic Information Systems, Digital Image Processing, Digital Photogrammetry and Global Positioning Systems. Landslide inventory databases are becoming available to more countries and several are now also available through the internet. A comprehensive landslide inventory is a must in order to be able to quantify both landslide hazard and risk. With respect to the environmental factors used in landslide hazard assessment, there is a tendency to utilize those data layers that are easily obtainable from Digital Elevation Models and satellite imagery, whereas less emphasis is on those data layers that require detailed field investigations. A review is given of the trends in collecting spatial information on environmental factors with a focus on Digital Elevation Models, geology and soils, geomorphology, land use and elements at risk.     

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

An operational multiscale system for hazards prediction, mapping, and response

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Methodology for preliminary assessment of Natech risk in urban areas

Concern for natural hazard-triggered technological disasters (Natech disasters) in densely populated and industrialized areas is growing. Residents living in urban areas subject to high natural hazard risk are often unaware of the potential for secondary disasters such as hazardous materials releases from neighboring industrial facilities, chemical storage warehouses or other establishments housing hazardous materials. Lessons from previous disasters, such as the Natech disaster during the Kocaeli earthquake in Turkey in 1999 call for the need to manage low frequency/high consequence events, particularly in today’s densely populated areas. However, there is little guidance available on how local governments and communities can assess Natech risk. To add to the problem, local governments often do not have the human or economic resources or expertise to carry out detailed risk assessments. In this article, we propose a methodology for preliminary assessment of Natech risk in urban areas. The proposed methodology is intended for use by local government officials in consultation with the public. The methodology considers possible interactions between the various systems in the urban environment: the physical infrastructure (e.g., industrial plants, lifeline systems, critical facilities), the community (e.g., population exposed), the natural environment (e.g., delicate ecosystems, river basins), and the risk and emergency management systems (e.g., structural and nonstructural measures). Factors related to vulnerability and hazard are analyzed and qualitative measures are recommended. Data from hazardous materials releases during the Kocaeli, Turkey earthquake of August 17, 1999 are used as a case study to demonstrate the applicability of the methodology. Limitations of the proposed methodology are discussed as well as future research needs.

Implementing low-cost landslide risk reduction: a pilot study in unplanned housing areas of the Caribbean

Landslides pose a serious physical and environmental threat to vulnerable communities living in areas of unplanned housing on steep slopes in the Caribbean. Some of these communities have, in the past, had to be relocated, at costs of millions of dollars, because of major slides triggered by tropical storm rainfall. Even so, evidence shows that: (1) risk reduction is a marginal activity; (2) there has been minimal uptake of hazard maps and vulnerability assessments and (3) there is little on-the-ground delivery of construction for risk reduction. This article directly addresses these issues by developing a low-cost approach to the identification of the potential pore pressure changes that trigger such slides we seek to address these three commentaries directly. A complex 45–60° slope site in St Lucia, West Indies was selected as a pilot for a modelling approach that uses numerical models (FLAC and CHASM) to verify the need for surface water management to effectively reduce landslide risk. Following the model confirmation, a series of drains were designed and constructed at the site. Post-construction evidence indicates the methodology to be sound, in that the site was stable in subsequent 1-in-1 to 1-in-4 year rainfall events. A critical feature of the approach is that it is community-based from data acquisition through to community members participating in construction.     

三峡库区生基包滑坡监测及成果分析

生基包滑坡监测属于三峡库区奉节县三期地质灾害监测预警项目之一,该滑坡位于长江左岸,临近人口稠密的安坪乡集镇,地理位置重要。三峡水库175 m蓄水后,其变形破坏特征有何表现？对航道安全运营有无潜在威胁？是否会对滑坡体上的重要建筑及村民生产生活构成危害？针对这些问题,首先分析了滑坡的工程地质特征及主要的影响因素;其次,确立以4种监测手段为主、人工巡查为辅的监测方案;通过对大地变形GPS、深部位移、滑坡推力等几种监测方法的运用及对其成果进行分析研究,以实例说明其在滑坡监测中的应用;再次,结合宏观人工巡查进行对照分析,以充分说明大地变形、深部位移和滑坡推力监测在实际运用中的可行性;最后,根据监测结论提出对生基包滑坡防治的建议。     

乌鲁木齐河流域地下水水质监测网设计

文章运用地下水易污性编图及污染源分布图法进行了乌鲁木齐河流域地下水水质监测网设计。共设计了130监测孔,现有46个监测孔,另需要84个新的监测孔。按监测类型分为面源监测点22个,点源监测点87个,重点水源地与泉水监测点21个。按监测运行分长期监测点55个,流域普查监测点75个。普查监测点监测频率为1次/5年,长期监测点监测频率为1次/年。首期有针对性地在污染严重的柴窝堡新化厂排污区、乌鲁木齐河谷老排污区、米泉污灌区、米泉工业污染区、老龙河污染区取了25个污染水样测试分析,结果显示地下水已经受到严重污染。