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.     

Contrasting metamorphic evolution of metasedimentary rocks from the Çine and Selimiye nappes in the Anatolide belt, western Turkey

Abstract P–T conditions, mineral isograds, the relation of the latter to foliation planes and kinematic indicators are used to elucidate the tectonic nature and evolution of a shear zone in an orogen exhumed from mid‐crustal depths in western Turkey. Furthermore, we discuss whether simple monometamorphic fabrics of rock units from different nappes result from one single orogeny or are related to different orogenies. Metasedimentary rocks from the Çine and Selimiye nappes at the southern rim of the Anatolide belt of western Turkey record different metamorphic evolutions. The Eocene Selimiye shear zone separates both nappes. Metasedimentary rocks from the Çine nappe underneath the Selimiye shear zone record maximum P–T conditions of about 7 kbar and >550 °C. Metasedimentary rocks from the overlying Selimiye nappe have maximum P–T conditions of 4 kbar and c. 525 °C near the base of the nappe. Kinematic indicators in both nappes are related to movement on the Selimiye shear zone and consistently show a top‐S shear sense. Metamorphic grade in the Selimiye nappe decreases structurally upwards as indicated by mineral isograds defining the garnet‐chlorite zone at the base, the chloritoid‐biotite zone and the biotite‐chlorite zone at the top of the nappe. The mineral isograds in the Selimiye nappe run parallel to the regional S_R foliation, parallel the Selimiye shear zone and indicate that the Selimiye shear zone formed during this prograde greenschist to lower amphibolite facies metamorphic event but remained active after the peak of metamorphism. ⁴⁰Ar/³⁹Ar mica ages and the tectonometamorphic relationship with the Eocene Cyclades–Menderes thrust, which occurs above the Selimiye nappe in the study area, suggests an Eocene age of metamorphism in the Selimiye nappe. Metasedimentary rocks of the Çine nappe 20–30 km north of the Selimiye shear zone record maximum P–T conditions of 8–11 kbar and 600–650 °C. An age of about 550 Ma is indicated for amphibolite facies metamorphism and associated top‐N shear in the orthogneiss of the Çine nappe. Our study shows that simple monophase tectonometamorphic fabrics do not always indicate a simple orogenic development of a nappe stack. Preservation in some areas and complete overprinting of those fabrics in other areas apparently occur very heterogeneously.     

Zircon fission‐track ages from Newfoundland—A proxy for high geothermal gradients and exhumation before opening of the Central Atlantic Ocean

Following Appalachian orogenesis, metamorphic rocks in central Newfoundland were exhumed and reburied under Tournaisian strata. New zircon fission‐track (ZFT) ages of metamorphic rocks below the Tournaisian unconformity yield post‐depositionally reset ages of 212–235 Ma indicating regional fluid‐absent reheating to at least ≥220°C. Post‐Tournaisian sedimentary thicknesses in surrounding basins show that burial alone cannot explain such temperatures, thus requiring that palaeo‐geothermal gradients increased to ≥30–40°C/km before final late Triassic accelerated cooling. We attribute these elevated palaeo‐geothermal gradients to localized thermal blanketing by insulating sediments overlying radiogenic high‐heat‐producing granitoids. Late Triassic rifting and magmatism before break up of Pangaea likely also contributed to elevated heat flow, as well as uplift, triggering late Triassic accelerated cooling and exhumation. Thermochronological ages of 240–200 Ma are seen throughout Atlantic Canada, and record rifting and basaltic magmatism on the conjugate margins of the Central Atlantic Ocean preceding the onset of oceanic spreading at ~190 Ma.     

Porphyroclast systems as kinematic indicators

Porphyroclasts of relatively strong minerals in mylonites commonly have an internal monoclinic shape symmetry defined by tails of dynamically recrystallized material. The geometry of a porphyroclast and its tails, called a ‘porhyroclast system’, can serve as a valuable indicator of the sense of vorticity. Porphyroclast systems have been divided into σ- and δ-types on the basis of the geometry of the tails. σ-Types have wedge-shaped recrystallized tails whose median lines lie on opposite sides of a reference plane parallel to the tails and containing the symmetry axis for the system. σ-Types are further subdivided into a σ_a-types, in which the porphyroclast is isolated in a relatively homogeneous matrix, and σ_b-types, in which the porphyroclast system is associated with a shear band foliation in the matrix. δ-Types typically have narrow recrystallized tails whose median lines cross the reference plane adjacent to the porphyroclast. Consequently, embayments of matrix material occur adjacent to the porphyroclasts and the tails display characteristic bends.A porphyroclast system in a mylonite develops when the relatively weak dynamically recrystallized grain aggregate in the porphyroclast mantle changes its shape due to non-coaxial flow in the adjacent matrix. This behaviour has been modelled in shear box experiments. Passive marker lines around rigid cylinders embedded in silicone putty were subjected to simple shear. The experiments were modified to simulate a change in recrystallization rate (R) with respect to rate of deformation (γ) by decreasing the diameter of the rigid cylinder during deformation at variable rates. The ratio R/γ appears to be one of the most important factors in determining which porphyroclast system will develop. At high R/γ values, flow of recrystallized material away from the porphyroclast is continuously appended by the production of new grains and wedge-shaped σ_a-type tails develop. At low R/γ values, relatively few new grains are added to the tails which become thinned and deflected by drag due to the spinning motion of the porphyroclast. In addition, most porphyroclast systems at low shear strains are of σ_a-type or lack monoclinic symmetry, whereas δ-types are only developed at high shear strain values. Complex porphyroclast systems, characterized by two generations of tails, are observed in many of the natural and model shear zones studied and may form due to fluctuating R/γ. Conditions that allow isolated σ_a- and δ-type porphyroclast systems to be used as sense of vorticity indicators are: the systems should have a monoclinic shape symmetry; matrix grain size should be small with respect to porphyroclast size; matrix fabric should be homogeneous; deformation history should be simple, and observations should be made on sections normal to the inferred bulk vorticity vector for the mylonite.     

δ objects as a gauge for stress sensitivity of strain rate in mylonites

Our understanding of the flow properties of deforming rocks in the Earth's lithosphere is mainly based on theoretical work and on the extrapolation of high-strain-rate experimental data to the low strain rates of rock deformation in nature. The geometry of structures in naturally deformed rocks can be an additional source of information on the rheology of the lithosphere. Flow experiments show that the geometry of a mantle of recrystallised material around a rigid object can be used to distinguish between a linear or power-law relation of differential stress and strain rate in strongly deformed rocks such as mylonites. This means that it is possible to use geometrical patterns in deformed rocks to constrain aspects of the rheology and to check the applicability of theoretical models and the extrapolation of laboratory data on creep in rocks to nature.     

Palaeostress analysis of the northern Nijar and southern Vera basins: constraints for the Neogene displacement history of major strike-slip faults in the Betic Cordilleras, SE Spain

This paper presents the results of a detailed structural analysis of the northern Nijar and southern Vera basins with special emphasis on the evolution of the regional stress field and the associated timing of movement of the Serrata, Gafarillos and Palomares strike-slip fault zones. These major fault zones control the Neogene deformation of the SE Internal Betic Cordilleras in Spain. Detailed stress analysis on Neogene sediments of the Vera and Nijar basins shows a strike-slip regime with NW–SE-oriented subhorizontal maximum principal stress (σ₁) during Tortonian and earliest Messinian times. Under the influence of this stress field, dextral displacement along the N090E-trending Gafarillos fault zone resulted in deformation of the sediments of the southern Sorbas and northeastern Nijar basins. During the early Messinian a clock-wise rotation of the stress field occurred. Stress analysis in rocks with late–early Messinian up to Quaternary ages in the Nijar and Vera basins indicates a strike-slip regime with N–S-oriented subhorizontal maximum principal stress (σ₁). Under the influence of this stress field the main activity along the N010E-striking Palomares strike-slip fault zone took place, resulting in deformation of the Neogene sediments of the southeastern Vera basin and culminating in a maximum sinistral displacement of more than 20 km. At the same time the stress field was not suitably oriented to exert a large shear component on the Gafarillos fault zone, which activity ended after the earliest Messinian. Fault and outcrop patterns of syntectonic Neogene sediments in the Vera basin show that displacement along the Palomares fault zone decreased at the end of the Middle Miocene although minor displacement phases may still have occurred during the Late Miocene and possibly even Pliocene. From the Middle Miocene onward, deformation in the Nijar basin was controlled by sinistral displacement along the N040E-trending Serrata strike-slip fault zone.     

Geology of the Bozdag area, central Menderes massif, SW Turkey: Pan-African basement and Alpine deformation

The Menderes massif consists of a Precambrian Core Series that preserves evidence for a polymetamorphic history and a Paleozoic/Mesozoic Cover Series that experienced only the Alpine tectonometamorphic evolution. Structural, petrographic, and geochronologic investigations in the central Menderes massif demonstrate that (a) part of the metamorphic and structural evolution of the Precambrian basement is older than the undeformed 551±1.4-Ma-old Birgi metagranite, and (b) inferred Alpine fabrics overprinting the Cover Series largely have the same attitudes as the old structures in the much older Core Series. The inferred Alpine fabrics include both contractional and extensional structures. Contraction under greenschist to amphibolite facies conditions resulted in the imbrication of the Core and Cover Series and generated an inverted metamorphic sequence by north-directed thrusting. During Alpine extension, most of the south-dipping thrust faults were reactivated as extensional shear zones under decreasing greenschist facies conditions.     

Computer experiments to investigate complex fibre patterns in natural antitaxial strain fringes

Antitaxial non‐deforming strain fringes from Lourdes, France, show complex quartz, calcite and chlorite fibre patterns that grew around pyrite in a slate during non‐coaxial progressive deformation. Development of these fringes was modelled using a computer program ‘Fringe Growth 2.0’ which can simulate incremental growth of crystal fibres around core‐objects of variable shape. It uses object‐centre paths as input, which are obtained from fibre patterns in thin section. The numerical experiments produced fibre patterns that show complex intergrowth of displacement‐controlled, face‐controlled and intermediate fibres similar to those in the natural examples. The direction of displacement‐controlled growth is only dependent on the relative movement between core‐object and fringe, so that core‐object rotation with respect to the fringe influences the fibre patterns and produces characteristic asymmetric fibre curvature. Object‐centre paths should be used for kinematic analysis of strain fringes instead of single fibres since these paths represent the fringe as a whole. The length along the path can be interpreted in terms of finite strain and path curvature in terms of rigid body rotation of fringes with respect to an external reference frame.     

Empirical prediction of coseismic landslide dam formation

In this study we develop an empirical method to estimate the volume threshold for predicting coseismic landslide dam formation using landscape parameters obtained from digital elevation models (DEMs). We hypothesize that the potential runout and volume of landslides, together with river features, determine the likelihood of the formation of a landslide dam. To develop this method, a database was created by randomly selecting 140 damming and 200 non‐damming landslides from 501 landslide dams and > 60 000 landslides induced by the M_w 7.9 2008 Wenchuan earthquake in China. We used this database to parameterize empirical runout models by stepwise multivariate regression. We find that factors controlling landslide runout are landslide initiation volume, landslide type, internal relief (H) and the H/L ratio (between H and landslide horizontal distance to river, L). In order to obtain a first volume threshold for a landslide to reach a river, the runout regression equations were converted into inverse volume equations by taking the runout to be the distance to river. A second volume threshold above which a landslide is predicted to block a river was determined by the correlation between river width and landslide volume of the known damming landslides. The larger of these two thresholds was taken as the final damming threshold. This method was applied to several landslide types over a fine geographic grid of assumed initiation points in a selected catchment. The overall prediction accuracy was 97.4% and 86.0% for non‐damming and damming landslides, respectively. The model was further tested by predicting the damming landslides over the whole region, with promising results. We conclude that our method is robust and reliable for the Wenchuan event. In combination with pre‐event landslide susceptibility and frequency–size assessments, it can be used to predict likely damming locations of future coseismic landslides, thereby helping to plan emergency response. Copyright © 2014 John Wiley & Sons, Ltd.     

Rock slope instability assessment using spatially distributed structural orientation data in Darjeeling Himalaya (India)

We discuss a geographic information system (GIS)‐based methodology for rock slope instability assessment based on geometrical relationships between topographic slopes and structural discontinuities in rocks. The methodology involves (a) regionalization of point observations of orientations (azimuth and dip) of structural discontinuities in rocks in order to generate a digital structural model (DStM), (b) testing the kinematical possibility of specific modes of rock slope failures by integrating DStMs and digital elevation model (DEM)‐derived slope and aspect data and (c) computation of stability scenarios with respect to identified rock slope failure modes. We tested the methodology in an area of 90 km² in Darjeeling Himalaya (India) and in a small portion (9 km²) within this area with higher density of field structural orientation data. The results of the study show better classification of rock slope instability in the smaller area with respect to known occurrences of deep‐seated rockslides than with respect to shallow translational rockslides, implying that structural control is more important for deep‐seated rockslides than for shallow translational rockslides. Results of scenario‐based analysis show that, in rock slopes classified to be unstable, stress‐induced rock slope instability tends to increase with increasing level of water saturation. The study demonstrates the usefulness of spatially distributed data of orientations of structural discontinuities in rocks for medium‐ to small‐scale classification of rock slope instability in mountainous terrains. Copyright © 2010 John Wiley & Sons, Ltd.