Were merely storm-landslides driven by the 2015-2016 Niño in the Mendoza River valley?

Since Holocene time, above-mean precipitations recorded during the El Niño warm ENSO phase have been linked to the occurrence of severe debris flows in the arid Central Andes. The 2015–2016 El Niño, for its unusual strength, began driving huge and dangerous landslides in the Central Andes (32°) in the recent South Hemisphere summer. The resulting damages negatively impacted the regional economy. Despite this, causes of these dangerous events were ambiguously reported. For this reason, a multidisciplinary study was carried out in the Mendoza River valley. Firstly, a geomorphological analysis of affected basins was conducted, estimating morphometric parameters of recorded events such as velocity, stream flow, and volume. Atmospheric conditions during such events were analyzed, considering precipitations, snow cover, temperature range, and the elevation of the zero isotherm. Based on our findings, the role of El Niño on the slope instability in the Central Andes is more complex in the climate change scenario. Even though some events were effectively triggered by intense summer rainstorm following expectations, the most dangerous events were caused by the progressive uplifting of the zero isotherm in smaller basins where headwaters are occupied by debris rock glaciers. Our research findings give light to the dynamic coupled system ENSO–climate change–landslides (ECCL) at least in this particular case study of the Mendoza River valley. Landslide activity in this Andean region is driven by wetter conditions linked to the ENSO warm phase, but also to progressive warming since the twentieth century in the region. This fact emphasizes the future impact of the natural hazards on Andean mountain communities.     

Precursory slope distress prior to the 2010 Mount Meager landslide,British Columbia

In 2010, the south flank of Mount Meager failed catastrophically, generating the largest (53 ± 3.8 × 10⁶ m³) landslide in Canadian history. We document the slow deformation of the edifice prior to failure using archival historic aerial photographs spanning the period 1948–2006. All photos were processed using Structure from Motion (SfM) photogrammetry. We used the SfM products to produce pre-and post-failure geomorphic maps that document changes in the volcanic edifice and Capricorn Glacier at its base. The photographic dataset shows that the Capricorn Glacier re-advanced from a retracted position in the 1980s then rapidly retreated in the lead-up to the 2010 failure. The dataset also documents 60 years of progressive development of faults, toe bulging, and precursory failures in 1998 and 2009. The 2010 collapse was conditioned by glacial retreat and triggered by hot summer weather that caused ice and snow to melt. Meltwater increased pore water pressures in colluvium and fractured rocks at the base of the slope, causing those materials to mobilize, which in turn triggered several secondary failures structurally controlled by lithology and faults. The landslide retrogressed from the base of the slope to near the peak of Mount Meager involving basement rock and the overlying volcanic sequence. Elsewhere on the flanks of Mount Meager, large fractures have developed in recently deglaciated areas, conditioning these slopes for future collapse. Potential failures in these areas have larger volumes than the 2010 landslide. Anticipated atmospheric warming over the next several decades will cause further loss of snow and glacier ice, likely producing additional slope instability. Satellite- and ground-based monitoring of these slopes can provide advanced warning of future landslides to help reduce risk in populated regions downstream.     

Sensitivity analysis of automatic landslide mapping: numerical experiments towards the best solution

The automatic detection of landslides after major events is a crucial issue for public agencies to support disaster response. Pixel-based approaches (PBAs) are widely used in the literature for various applications. However, the accuracy of PBAs in the case of automatic landslide mapping (ALM) is affected by several issues. In this study, we investigated the sensitivity of ALM using PBA through digital terrain models (DTMs). The analysis, carried out in a study area of Poland, consisted of the following steps: (1) testing the influence of selected DTM resolutions for ALM, (2) assessing the relevance of diverse landslide morphological indicators for ALM, and (3) assessing the sensitivity to landslide features for a selected size of moving window (kernel) calculations for ALM. Ultimately, we assessed the performance of three classification methods: maximum likelihood (ML), feed-forward neural network (FFNN), and support vector machine (SVM). This broad analysis, as combination of grid cell resolution, surface derivatives calculation, and performance classification methods, is the challenging aspect of the research. The results of almost 500 experimental tests provide valuable guidelines for experts performing ALM. The most important findings indicate that feature sensitivity in the case of kernel size increases with coarser DTM resolution; however, the peak of the optimal feature performance for the selected study area and landslide type was demonstrated for a resolution of 20 m. Another finding indicated that in combining a set of topographic variables, the optimal performance was acquired for a DTM resolution of 30 m and the support vector machine classification. Moreover, the best performance of the identification is represented for SVM classification.     

Residual-state creep of clastic soil in a reactivated slow-moving landslide in the Three Gorges Reservoir Region,China

We study the creep properties of clastic soil in residual state. The intact samples are taken from a reactivated slow-moving landslide in the Three Gorges Reservoir Region, China. Firstly, the patterns of the landslide movement are analysed based on recent monitoring data, which indicate that the soil within the shear zone is undergoing two deformation processes: a creep phase, characterised by different creep rates, and a dormant phase. We then study the creep behaviour of the soil samples through a series of ring shear creep tests under various shear stress conditions. The creep response depends strongly on the ratio of the shear stress to the residual strength, and the normal effective stress, whereas the creep rate decreases due to strength regain. The long-term strength of the clastic soil is close to the residual strength. Therefore, the residual strength obtained from conventional shear test, which is less time consuming than creep test, can be used in long-term stability analyses of creeping landslides.     

Study on the deformation and failure modes of filling slope in loess filling engineering: a case study at a loess mountain airport

Lvliang airport is a typical loess filling engineering located in 20.5 km north of Lvliang City in Shanxi Province, China. By the end of March 2012, 14 fissures extending more than 7.5 m were observed in a loess-filled slope, of which the longest fissure is up to 82 m. Field monitoring and laboratory tests have been performed to investigate the slope failure modes. The test program includes wetting tests on unsaturated compacted samples and stress path tests on saturated samples. Field monitoring and observations show that differential settlement caused by non-homogeneity in compacted loess density might lead to the formation of fissures in the loess-filled slope. It was founded that the wetting deformation contributed to the development of differential settlement. Fissures are the essential factor for the loess-filled slope failure. Four deformation stages exhibit in the loess-filled slope prior final failure including development of the fissures, softening of the compacted loess, creeping of the slope leading edge and fissuring of the trailing edge and formation of the through-sliding surface. Development of the sliding surface mainly includes upward and downward expansion of the fissures. Upward expansion is a wetting failure process in loading condition, while downward expansion is a load-off wetting process. In addition, development of the sliding surface is accelerated by softening of the compacted soils as a result of water infiltration. Therefore, the key for taking countermeasures against filling landslides is to monitor and control the development of differential settlement and fissures in the slope shoulders. Digging out and extra-filling the fissures are an effective way for preventing these landslides.     

Application of incomplete similarity theory to the estimation of the mean velocity of debris flows

The mean velocity of debris flow is one of the most important parameters in the design of mitigation structures and in quantitative risk analysis. This study develops a model to predict the mean debris flow velocity observed in the field by applying the incomplete similarity argument. An equation for estimating the Darcy-Weisbach resistance coefficient for debris flows with a volumetric sediment concentration larger than 0.19 is accordingly derived using 128 sets of observation data from nine Chinese gullies, in which both the effect of the volumetric sediment concentration and channel slope on resistance are considered. The derived equation is then verified and compared against five previously published equations by using 61 sets of published observation data from six gullies located in four countries. The applications of the proposed equation are discussed, and the improvements made using the proposed equation are clearly very significant when compared with the previously published equations.     

Identification of shear strength and seismic coefficient by back analyzing surficial slides in the 2004 Mid-Niigata prefecture earthquake

Sliding of natural and artificial slopes generally occurs during or following strong earthquakes. Such sliding is greatly affected by a combination of geological conditions and earthquake loading. Earthquake-induced landslides often cause more damage to infrastructure and human lives than the earthquake itself. Pseudo-static analysis is widely implemented as one of several design methods used in engineering practice to assess the seismic stability of natural and artificial slopes. However, the most important issue of pseudo-static analysis is to select the most appropriate method for measuring seismic coefficient. In order to investigate this, back analysis was conducted for surficial slides subjected to strong ground motion during the 2004 Mid-Niigata prefecture earthquake in Japan. This paper surveyed the stochastic properties of earthquake-induced surficial slides and clearly showed that the obtained results were applicable to back analysis of shear strength and seismic coefficient. In back analysis, soil properties such as soil strength and density and sliding depth were assumed as random variables owing to their uncertainties. Seismic coefficient is also assumed to be a random variable and varies with distance from the epicenter fault line. The analysis of 4504 recorded surficial slides clearly shows a unique relationship of landslide occurrence ratio with slope angle and distance from the epicenter fault line. This study’s results enhance the calculation of the shear strength of weathered soil covering slopes and the horizontal seismic coefficient through back analysis procedure. By considering possible stochastic properties of variables, some case studies were implemented in the back analysis.     

Comparison of single- and dual-permeability models in simulating the unsaturated hydro-mechanical behavior in a rainfall-triggered landslide

Landslide-prone slopes in earthquake-affected areas commonly feature heterogeneity and high permeability due to the presence of cracks and fissures that were caused by ground shaking. Landslide reactivation in heterogeneous slope may be affected by preferential flow that was commonly occurred under heavy rainfall. Current hydro-mechanical models that are based on a single-permeability model consider soil as a homogeneous continuum, which, however, cannot explicitly represent the hydraulic properties of heterogeneous soil. The present study adopted a dual-permeability model, using two Darcy-Richards equations to simulate the infiltration processes in both matrix and preferential flow domains. The hydrological results were integrated with an infinite slope stability approach, attempting to investigate the hydro-mechanical behavior. A coarse-textured unstable slope in an earthquake-affected area was chosen for conducting artificial rainfall experiment, and in the experiment slope, failure was triggered several times under heavy rainfall. The simulated hydro-mechanical results of both single- and dual-permeability model were compared with the measurements, including soil moisture content, pore water pressure, and slope stability conditions. Under high-intensity rainfall, the measured soil moisture and pore water pressure at 1-m depth showed faster hydrological response than its simulations, which can be regarded as a typical evidence of preferential flow. We found the dual-permeability model substantially improved the quantification of hydro-mechanical processes. Such improvement could assist in obtaining more reliable landslide-triggering predication. In the light of the implementation of a dual-permeability model for slope stability analysis, a more flexible and robust early warning system for shallow landslides hazard in coarse-textured slopes could be provided.