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.     

Study of landslides at the archaeological site of Abakainon necropolis (NE Sicily) by geomorphological and geophysical investigations

This paper reports on the investigation of the old landslide body on which the necropolis of Abakainon, an Italian archaeological site, was built about 2400 years ago. This site is located in a geologically complex area of northeastern Sicily, where active tectonics has strongly modified the morphological landscape, as testified by the meandering path of watercourses, widespread presence of paleosurfaces, and deviations of crest line. Previous studies have hypothesized that the necropolis, suddenly destroyed by a catastrophic event in historical times, could have been built on an old landslide body. This represents the starting point of this paper, focused on geomorphological analysis and geophysical surveys, which allowed discovering peculiar elements in support of the landslide theory. In particular, the elaboration of a Digital Terrain Model allowed locating the main morphostructural elements, thus highlighting how the occurrence of landslides in the study area is strongly linked to recent tectonics. Passive seismic surveys, carried out through the acquisition of ambient noise, allowed detecting the presence of impedance contrasts, which were associated to peculiar features, such as the landslide body itself. Moreover, local strong directional effects testify the presence of faults, affecting even the necropolis area. Such results, further corroborated by the outcomes of an electrical resistivity survey, lay the foundations for interesting consideration on the possible cause of the destruction of this site, probably related to the occurrence of a historical earthquake in the first century AD.     

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.     

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.     

Relationships between natural terrain landslide magnitudes and triggering rainfall based on a large landslide inventory in Hong Kong

Rain-induced landslides are recognized as one of the most catastrophic hazards on hilly terrains. To develop strategies for landslide risk assessment and management, it is necessary to estimate not only the rainfall threshold for the initiation of landslides, but also the likely magnitudes of landslides triggered by a storm of a given intensity. In this study, the frequency distributions of both open hillside landslides and channelized debris flows in Hong Kong are established on the basis of the Enhanced Natural Terrain Landslide Inventory (ENTLI) with 19,763 records in Hong Kong up to 2013. The landslide magnitudes are measured in terms of the number, scar area, volume, or density of landslides. The mean values of the scar areas and volumes are 55.2 m² and 102.0 m³, respectively, for the open hillside landslides and 91.3 m² and 166.5 m³, respectively, for the channelized debris flows. Empirical correlations between the numbers, scar areas, and volumes of hillside landslides or channelized debris flows and the maximum rolling rainfall intensities of different periods have been derived. The maximum rolling 4- to 24-h rainfall amounts provide better predictions compared with those with the maximum rolling 1-h rainfall. Maximum rolling rainfall intensity-duration thresholds identifying the likely rainfall conditions that yield natural terrain landslides or debris flows of different magnitudes are also proposed. The initiation rainfall thresholds are identified as 75, 90, 100, 120, 150, 180, and 200 mm for the maximum rolling 1-, 2-, 4-, 6-, 8-, 12-, and 24-h rainfall, respectively.     

Assessing rockfall susceptibility in steep and overhanging slopes using three-dimensional analysis of failure mechanisms

Rockfalls strongly influence the evolution of steep rocky landscapes and represent a significant hazard in mountainous areas. Defining the most probable future rockfall source areas is of primary importance for both geomorphological investigations and hazard assessment. Thus, a need exists to understand which areas of a steep cliff are more likely to be affected by a rockfall. An important analytical gap exists between regional rockfall susceptibility studies and block-specific geomechanical calculations. Here we present methods for quantifying rockfall susceptibility at the cliff scale, which is suitable for sub-regional hazard assessment (hundreds to thousands of square meters). Our methods use three-dimensional point clouds acquired by terrestrial laser scanning to quantify the fracture patterns and compute failure mechanisms for planar, wedge, and toppling failures on vertical and overhanging rock walls. As a part of this work, we developed a rockfall susceptibility index for each type of failure mechanism according to the interaction between the discontinuities and the local cliff orientation. The susceptibility for slope parallel exfoliation-type failures, which are generally hard to identify, is partly captured by planar and toppling susceptibility indexes. We tested the methods for detecting the most susceptible rockfall source areas on two famously steep landscapes, Yosemite Valley (California, USA) and the Drus in the Mont-Blanc massif (France). Our rockfall susceptibility models show good correspondence with active rockfall sources. The methods offer new tools for investigating rockfall hazard and improving our understanding of rockfall processes.     

Location-based environmental factors contributing to rainfall-triggered debris flows in the Ba river catchment,northwest Viti Levu island,Fiji

Mass movements in tropical Pacific small island developing states (SIDS) can be devastating although studies are relatively few and contributing environmental factors are not often investigated in detail. On 25 January 2012, following 3 days of heavy monsoonal rainfall (c. 550 mm) during a La Niña episode, more than 150 debris flows were triggered in the western part of the Ba river catchment of northwest Viti Levu island, Fiji. Reconnaissance field survey and geographical information system (GIS) analyses using high-resolution satellite imagery were carried out to investigate factors that may have led to the occurrence of the debris flows in the catchment. We evaluated the correlation between the density of mass movements (number of mass movements/km²) and several continuous variables using data measured within the GIS. There was a weak but significant positive correlation between mass movement density and elevation (r = 0.38, p value < 0.01), cyclonic precipitation (r = 0.37, p value < 0.01) and stream density (r = 0.31, p value < 0.01). Ninety-three percent of the mass movements occur within a plantation of Pinus caribaea (Caribbean pine) on slopes oriented mainly to the northeast and east on (trade) windward slopes and may be significant factors for their development. Although forests generally have a stabilizing effect on slopes, the plantation at Ba was a mature stand on its second plantation cycle and is a species that has a shallow rooting system making it more susceptible to failure.     

Implementing landslide path dependency in landslide susceptibility modelling

Landslide susceptibility modelling—a crucial step towards the assessment of landslide hazard and risk—has hitherto not included the local, transient effects of previous landslides on susceptibility. In this contribution, we implement such transient effects, which we term “landslide path dependency”, for the first time. Two landslide path dependency variables are used to characterise transient effects: a variable reflecting how likely it is that an earlier landslide will have a follow-up landslide and a variable reflecting the decay of transient effects over time. These two landslide path dependency variables are considered in addition to a large set of conditioning attributes conventionally used in landslide susceptibility. Three logistic regression models were trained and tested fitted to landslide occurrence data from a multi-temporal landslide inventory: (1) a model with only conventional variables, (2) a model with conventional plus landslide path dependency variables, and (3) a model with only landslide path dependency variables. We compare the model performances, differences in the number, coefficient and significance of the selected variables, and the differences in the resulting susceptibility maps. Although the landslide path dependency variables are highly significant and have impacts on the importance of other variables, the performance of the models and the susceptibility maps do not substantially differ between conventional and conventional plus path dependent models. The path dependent landslide susceptibility model, with only two explanatory variables, has lower model performance, and differently patterned susceptibility map than the two other models. A simple landslide susceptibility model using only DEM-derived variables and landslide path dependency variables performs better than the path dependent landslide susceptibility model, and almost as well as the model with conventional plus landslide path dependency variables—while avoiding the need for hard-to-measure variables such as land use or lithology. Although the predictive power of landslide path dependency variables is lower than those of the most important conventional variables, our findings provide a clear incentive to further explore landslide path dependency effects and their potential role in landslide susceptibility modelling.     

Palynology of Triassic–Jurassic boundary sections in northern Switzerland

A first palynostratigraphic scheme of Upper Triassic deposits in northern Switzerland was established based on spore-pollen associations and dinoflagellate cyst records from the upper part of the Upper Triassic Klettgau Formation and the lower part of the Lower Jurassic Staffelegg Formation. Drill cores from the Adlerberg region (Basel Tabular Jura) and from Weiach (northern part of Canton Zurich) as well as from an outcrop at the Chilchzimmersattel (Basel Folded Jura) were studied and five informal palynological associations are distinguished. These palynological associations correlate with palynological association of the Central European Epicontinental Basin and the Tethyan realm and provide a stratigraphic framework for the uppermost Triassic sediments in northern Switzerland. Throughout the uppermost Triassic to Jurassic palynological succession a remarkable prominence of Classopollis spp. is observed. Besides Classopollis spp. the three Rhaetian palynological associations A to C from the Upper Triassic Belchen Member include typical Rhaetian spore-pollen and dinoflagellate taxa (e.g., Rhaetipollis germanicus, Geopollis zwolinskae, Rhaetogonyaulax rhaetica, and Dapcodinium priscum). Association B differs from association A in a higher relative abundance of the sporomorph taxa Perinopollenites spp. and the consistent occurrence of Granuloperculatipollis rudis and Ricciisporites tuberculatus. Spore diversity is highest in the late Rhaetian palynological association C and includes Polypodiisporites polymicroforatus. A Rhaetian age for the Belchen Member is confirmed by palynological associations A–C, but there is no record of the latest Rhaetian and the earliest Jurassic. In contrast to the Rhaetian palynological associations the Early Jurassic associations W and D include Pinuspollenites spp., Trachysporites fuscus (in association W), and Ischyosporites variegatus. In the view of the end-Triassic mass extinction and contemporaneous environmental changes the described palynofloral succession represents the pre-extinction phase (associations A and B) including a distinct transgression, the extinction phase (association C) associated with a regression, and the post-extinction phase (association W).     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.