Late Cenozoic geomorphologic signal of Andean forearc deformation and tilting associated with the uplift and climate changes of the Southern Atacama Desert (26°S–28°S)

We analyze remarkable examples of the large ( 10,000 km²) and local-scale ( 100 km²) landscape forms related to Late Cenozoic geomorphologic evolution of the Andean forearc region in the Southern Atacama Desert. We also consider the continental sedimentary deposits, so-called “Atacama Gravels”, which are related to the degradation of the landscape during the Neogene. Our analysis integrates 1:50,000 field cartography, Landsat TM images observations,  1:1000 sedimentary logging data, and 50 m horizontal resolution topographic data to reconstruct the Late Cenozoic geomorphologic evolution of this region and discuss the factors that control it, i.e., Miocene aridification of the climate and Neogene Central Andean uplift. We determine that the Precordillera was already formed in the Oligocene and most of the present-day altitude of the Precordillera was reached before that time. Afterward, five episodes of geomorphologic evolution can be differentiated: (1) the development of an Oligocene deep incised drainage system cutting the uplifted Precordillera (up to 2000 m of vertical incision) and connecting it to the Ocean; followed by (2) the infilling of deep incised valleys by up to 400 m of Atacama Gravels. This infill started in the Early Miocene with the development of fluvial deposition and finished in the Middle Miocene with playa and playa lake depositions. We propose that playa-related deposition occurs in an endorheic context related to tectonic activity of the Atacama Fault System and Coastal Cordillera uplift. However, the upward sedimentologic variation in the Atacama Gravels evidences a progressive aridification of the climate. Subsequently, we have identified the effects of the Middle–Upper Miocene slow tectonic deformation: the Neogene Andean uplift is accommodated by a tilting or flexuring of the inner-forearc (Central Depression and Precordillera) related to some hundreds of meters of uplift in the Precordillera. This tilting or flexuring results in (3) the Middle Miocene re-opening of the valley network to the Pacific Ocean. Upper Miocene aridification, from arid to hyperarid, induces alluvial fans backfilling in the Central Depression (4) resulting in up to 50 m of Atacama Gravel deposition. Finally, in response to an increase in the rate of tilting, a new phase of vertical incision (up to 800 m in the Precordillera) allows the development of the canyon that crosses the forearc (5).     

Sedimentary impacts from landslides in the Tachia River Basin, Taiwan

A case study of coseismic landslides and post-seismic sedimentary impacts of landslides due to rainfall events was conducted in the Tachia River basin, Taichung County, central Taiwan. About 3000 coseismic landslides occurred in the basin during the M_L 7.3 Chi-Chi earthquake in 1999. The deposits from these landslides provided material for numerous debris flows induced by subsequent rainfall events. The estimated 4.1 × 10⁷ m³ of landslide debris produced in the upland area caused sediment deposition in riverbeds, and flash floods inundated downstream areas with sediment during torrential rains. The landslide frequency-size distributions for the coseismic landslides and the subsequent rainfall-induced landslides were analyzed to determine the sediment budgets of the post-seismic geomorphic response in the landslide-dominated basin. Both the coseismic and the rainfall-induced landslides show a power–law frequency-size distribution with a rollover. It was found that the rainfall-induced landslide magnitude was smaller than the coseismic one, and that both have comparable negative scaling exponents in cumulative form, of about − 2.0 for larger landslides (> 10^− 2 km²). This may be attributed to ongoing movement or reactivation of old landslides, and a natural stabilisation of small landslides between 10^− 4 and 10^− 2 km². It is proposed that the characteristics of geological formations and rainfall as well as changes in landslide area are reflected in the power–law distribution.     

Gigantic low-gradient landslides in the northern periphery of the Crimean Mountains (Ukraine)

Large-scale, low-gradient ancient landslides estimated at 5.4–18.9 km² in area and  0.2–1.2 km³ in volume have been studied in the northern hilly periphery of the Crimean Mountains (Ukraine). They originated on slopes along wide water gaps of rivers (Belbek, Kacha, Alma and Biyuk–Karasu) crossing the cuestas of the northern foothills. The slopes generally consist of slightly northward tilting Miocene (mainly Sarmatian) limestones overlying weak, clay-rich Lower Neogene–Palaeogene substratum with a significant content of smectite. Although the region is characterised by the least active contemporary morphodynamics within the Crimean Mountains, the landslides which were studied are of the same size or even larger than various types of landslides occupying active geomorphic domains of the highest mountain range in the southernmost part of the peninsula. The landslides are generally a spreading type, but the sliding mechanics were probably very complex, involving toppling, rotational slides, gravitational folding and translational block slides. All the landslides which were studied are located in the vicinity of regional faults and three of them have headscarps aligned along faults. A common feature is also a location close (within several km) to the Mesozoic suture zone which is the most important tectonic feature in the northern periphery of the Crimean Orogene. This suture was formerly classified as aseismic; however, evidence of strong, low-frequency palaeoearthquakes was collected during the last decade within both the Mesozoic suture and the low-lying northern part of the Crimean Peninsula. Radiocarbon dating of deposits associated with the landslides has revealed at least two phases of increased landslide-activity during the Late Glacial chronozone and Holocene epoch. The main landslide phase presumably took place at some time between the Late Glacial and Atlantic chronozones. Minor reactivation of landslide toes occurred during the Subatlantic chronozone and some of them have been active up to recent times. The first major landslide phase was possibly triggered by an earthquake, whereas late Holocene activity can be attributed both to seismic and hydroclimatic factors.     

Geomorphological study of long-term erosion on a tropical volcanic ocean island: Tahiti-Nui (French Polynesia)

A quantitative geomorphological study has been made on 27 river basins in Tahiti-Nui volcanic island (French Polynesia) to reconstruct the erosional evolution of a young oceanic island subjected to heavy tropical rainfall. Tahiti-Nui is composed of a main shield volcano cut by two huge landslides on each side of a main E–W rift zone. The northern landslide depression was rapidly buried by the construction of a second shield, the late activity of which overflowed the crest and then filled the southern landslide depression. The island is now volcanically inactive and is deeply dissected by erosion. The present geometries of the river basins are first compared using dimensionless parameters derived from a digital elevation model. The original volcanic surfaces are then reconstructed to estimate the volumes removed by erosion and determine the average rates of long-term erosion. The basins developed on the flanks of the main shield are wider, shallower, and gentler than the basins incising the post-landslide second shield, indicating a higher degree of evolution. Rainfall concentration on the windward (eastern) side of the island also contributed to increase the vertical lowering of the volcanic relief and the enlargement of the valleys. The magnitude of erosion, however, is neither directly linked with the age of the units incised nor with the differential amounts of rainfall. Erosion rates determined over the last 1 Myr range between 10^− 3 km³ kyr^− 1 and 0.25 km³ kyr^− 1. The highest values occur in the basins incising the main E–W rift zone and/or the lateral rims of the northern and southern landslide depressions. Long-term dissection has thus been enhanced along the geological discontinuities of the eruptive system. Deep erosion was first constrained along the axis of the main E–W rift zone, where numerous dykes compartmentalize the volcanic structure into large unstable blocks. Dykes most probably acted as mechanical discontinuities along which shallow gravitational landslides recurrently occurred. Such mass-wasting episodes produced significant amounts of debris, partly preserved as highly indurated sedimentary breccias of various ages exposed at various locations. Subsequent dissection of Tahiti-Nui was enhanced to the north and to the south, leading to the rapid evolution of the Papenoo and Taharuu drainage systems over the last 500 kyr. Long-term dissection on Tahiti-Nui has been responsible for the removal of at least 350 km³ of volcanic material from the surface, and for the partial exhumation of a shallow intrusive complex partly composed of coarse-grained plutonic rocks (gabbros and syenites) in the central part of the eruptive system. Structurally controlled erosion is thus a key component of landscape evolution on such high-relief oceanic tropical islands.     

Spatial patterns of old, deep-seated landslides: A case-study in the northern Ethiopian highlands

During the last decade, slope failures were reported in a 500 km² study area in the Geba–Werei catchment, northern Ethiopia, a region where landslides were not considered an important hazard before. Field observations, however, revealed that many of the failures were actually reactivations of old deep-seated landslides after land use changes. Therefore, this study was conducted (1) to explore the importance of environmental factors controlling landslide occurrence and (2) to estimate future landslide susceptibility. A landslide inventory map of the study area derived from aerial photograph interpretation and field checks shows the location of 57 landslides and six zones with multiple landslides, mainly complex slides and debris flows. In total 14.8% of the area is affected by an old landslide. For the landslide susceptibility modelling, weights of evidence (WofE), was applied and five different models were produced. After comparison of the models and spatial validation using Receiver Operating Characteristic curves and Kappa values, a model combining data on elevation, hillslope gradient, aspect, geology and distance to faults was selected. This model confirmed our hypothesis that deep-seated landslides are located on hillslopes with a moderate slope gradient (i.e. 5°–13°). The depletion areas are expected on and along the border of plateaus where weathered basalts rich in smectite clays are found, and the landslide debris is expected to accumulate on the Amba Aradam sandstone and upper Antalo limestone. As future landslides are believed to occur on inherently unstable hillslopes similar to those where deep-seated landslides occurred, the classified landslide susceptibility map allows delineating zones where human interventions decreasing slope stability might cause slope failures. The results obtained demonstrate that the applied methodology could be used in similar areas where information on the location of landslides is essential for present-day hazard analysis.     

Episodic wood loading in a mountainous neotropical watershed

The Upper Rio Chagres drains 414 km² of steep, mountainous terrain in central Panama. A tropical air mass thunderstorm on 10 July 2007 produced a flood across the basin that peaked at 720 m³ s^− 1 at a headwaters gage draining 17.5 km² and 1710 m³ s^− 1 at a downstream gage draining 414 km². The storm also triggered numerous landslides in the upper basin, which facilitated the formation of large logjams along portions of the channel where transport capacity of wood was reduced by a change in channel geometry such as a bend or channel expansion. During field work in February 2008, we characterized three jams with surface areas of 400–2450 m²; two of these jams resulted in storage of substantial (1100–8200 m³) sediment wedges upstream. We returned to these sites in March 2009 to document changes in the logjams and sediment storage. Drawing on observations made in the basin since 2002, and site visits during 2008 and 2009, we suggest that jams such as these last two years or less. We propose that wood dynamics in the Upper Chagres alternate between brief periods of moderate wood load in the form of large logjams and much longer periods of essentially no wood load, a situation that contrasts with the more consistent wood loads in catchments of similar size in temperate environments and with limited studies of more consistent wood load in tropical catchments with no landslides.     

Representative rainfall thresholds for landslides in the Nepal Himalaya

Measuring some 2400 km in length, the Himalaya accommodate millions of people in northern India and Pakistan, Nepal, Bhutan, and parts of other Asian nations. Every year, especially during monsoon rains, landslides and related natural events in these mountains cause tremendous damage to lives, property, infrastructure, and environment. In the context of the Himalaya, however, the rainfall thresholds for landslide initiation are not well understood. This paper describes regional aspects of rainfall thresholds for landslides in the Himalaya. Some 677 landslides occurring from 1951 to 2006 were studied to analyze rainfall thresholds. Out of the 677 landslides, however, only 193 associated with rainfall data were analyzed to yield a threshold relationship between rainfall intensity, rainfall duration, and landslide initiation. The threshold relationship fitted to the lower boundary of the field defined by landslide-triggering rainfall events is I = 73.90D^− 0.79 (I = rainfall intensity in mm h^− 1 and D = duration in hours), revealing that when the daily precipitation exceeds 144 mm, the risk of landslides on Himalayan mountain slopes is high. Normalized rainfall intensity–duration relationships and landslide initiation thresholds were established from the data after normalizing rainfall-intensity data with respect to mean annual precipitation (MAP) as an index in which N_I = 1.10D^− 0.59 (N_I = normalized intensity in h^− 1). Finally, the role of antecedent rainfall in causing landslides was also investigated by considering daily rainfall during failure and the cumulative rainfall to discover at what point antecedent rainfall plays an important role in Himalayan landslide processes. Rainfall thresholds presented in this paper are generalized so they can be used in landslide warning systems in the Nepal Himalaya.     

Landslides in the Central Coalfield (Cantabrian Mountains, NW Spain): Geomorphological features, conditioning factors and methodological implications in susceptibility assessment

A geomorphological study focussing on slope instability and landslide susceptibility modelling was performed on a 278 km² area in the Nalón River Basin (Central Coalfield, NW Spain). The methodology of the study includes: 1) geomorphological mapping at both 1:5000 and 1:25,000 scales based on air-photo interpretation and field work; 2) Digital Terrain Model (DTM) creation and overlay of geomorphological and DTM layers in a Geographical Information System (GIS); and 3) statistical treatment of variables using SPSS and development of a logistic regression model. A total of 603 mass movements including earth flow and debris flow were inventoried and were classified into two groups according to their size. This study focuses on the first group with small mass movements (10⁰ to 10¹ m in size), which often cause damage to infrastructures and even victims. The detected conditioning factors of these landslides are lithology (soils and colluviums), vegetation (pasture) and topography. DTM analyses show that high instabilities are linked to slopes with NE and SW orientations, curvature values between − 6 and − 0.7, and slope values from 16° to 30°. Bedrock lithology (Carboniferous sandstone and siltstone), presence of Quaternary soils and sediments, vegetation, and the topographical factors were used to develop a landslide susceptibility model using the logistic regression method. Application of “zoom method” allows us to accurately detect small mass movements using a 5-m grid cell data even if geomorphological mapping is done at a 1:25,000 scale.     

Quantitative landslide risk assessment and mapping on the basis of recent occurrences

A quantitative procedure for mapping landslide risk is developed from considerations of hazard, vulnerability and valuation of exposed elements. The approach based on former work by the authors, is applied in the Bajo Deba area (northern Spain) where a detailed study of landslide occurrence and damage in the recent past (last 50 years) was carried out. Analyses and mapping are implemented in a Geographic Information System (GIS).The method is based on a susceptibility model developed previously from statistical relationships between past landslides and terrain parameters related to instability. Extrapolations based on past landslide behaviour were used to calculate failure frequency for the next 50 years. A detailed inventory of direct damage due to landslides during the study period was carried out and the main elements at risk in the area identified and mapped. Past direct (monetary) losses per type of element were estimated and expressed as an average ‘specific loss’ for events of a given magnitude (corresponding to a specified scenario). Vulnerability was assessed by comparing losses with the actual value of the elements affected and expressed as a fraction of that value (0–1).From hazard, vulnerability and monetary value, risk was computed for each element considered. Direct risk maps (€/pixel/year) were obtained and indirect losses from the disruption of economic activities due to landslides assessed. The final result is a risk map and table combining all losses per pixel for a 50-year period. Total monetary value at risk for the Bajo Deba area in the next 50 years is about 2.4 × 10⁶ Euros.     

Clay slides in the Målselv valley, northern Norway: Characteristics, occurrence, and triggering mechanisms

The distribution of a large number of clay slides in the Målselv valley, northern Norway, is analysed and put into context with the stratigraphic organization of the valley-fill sediments. About 32% of the landslides larger than 10⁴ m³ occur close to the valley margins, where mud is either exposed or at shallow depth. About 57% of the landslides occur mid-valley, where relatively thin (< 15 m) coarse-grained deltaic sediments overlay fine-grained marine and glaciomarine sediments, and about 11% of the landslides occur in front of ice-contact deposits. The slide-prone areas are all characterized by the occurrence of heterogeneous sediments (interbedded clay, silt and sand), in addition to the presence of steep slopes eroded by rivers. The heterogeneous nature of the sediments probably enhanced groundwater drainage and leaching of salts from the clay, increasing sensitivity. Thus, the distribution and organization of the valley-fill sediments and groundwater drainage probably controlled the position of the slide scars and sliding planes. Since deglaciation of the valley (11,500 BP–present), isostatic rebound has enhanced fluvial incision and the creation of steep slopes due to a fall in relative sea level of up to 80 m.Arcuate-shaped, ‘bottleneck’ landslide scars ranging from c. 10⁴ to 10⁷ m³ in size is the dominant morphological signature of the slides, typical for quick clay slides or earth flows involving fluid mud. Their most common triggering mechanism is probably erosion at the toe slopes by the river Målselv or its tributaries. River erosion close to the valley margin, where glaciomarine and marine sediments are present, seems to give the most severe slides. The overall valley-fill organization controls the distribution of clay slides, which may apply to other fjord valleys having similar sediment distribution.     

Prolonged impact of earthquake-induced landslides on sediment yield in a mountain watershed: The Tanzawa region, Japan

To determine for how long a landslide affects sediment discharge, the sediment yields of 15 check-dam basins were compared with the time series of landslide distributions in a mountain basin in the Tanzawa region, central Japan. The distribution of sediment yield was quantitatively estimated from deposition in the sediment pools of check dams. The relationship between the landslide history and sediment discharge in the Nakagawa River basin was examined for an approximately 80-year period. Two major landslide events occurred during this period: the 1923 Kanto Earthquake and the 1972 disaster caused by heavy rainfall. The resulting trend in sediment discharge of the whole basin, estimated using reservoir sedimentation in the Miho Dam at its base, was nearly constant, with high sediment discharge (2897 m³ km^− 2 yr^− 1) in the intervening quarter-century, despite the recovery of vegetation on landslide areas in this period. Comparisons of the landslide distributions resulting from the two disasters, the sediment yields of check-dam basins, and the sediment discharge of the whole basin indicate that recent sediment discharge contains landslide debris that was originated by the Kanto Earthquake that occurred over 80 years ago. Thus, to understand high sediment discharge, it is essential to investigate not only the current basin condition and recent events, but also the landslide history of the basin for at least the previous 100 years.     

Dust emission by off-road driving: Experiments on 17 arid soil types, Nevada, USA

Field experiments were conducted in Nellis Dunes Recreational Area (Clark County, Nevada, USA) to investigate emission of dust produced by off-road driving. Experiments were carried out with three types of vehicles: 4-wheelers (quads), dirt bikes (motorcycles) and dune buggies, on 17 soil types characteristic for a desert environment. Tests were done at various driving speeds, and emissions were measured for a large number of grain size fractions. This paper reports the results for two size fractions of emissions: PM10 (particles < 10 μm) and PM60 (particles < 60 μm). The latter was considered in this study to be sufficiently representative of the total suspendable fraction (TSP). Off-road driving was found to be a significant source of dust. However, the amounts varied greatly with the type of soil and the characteristics of the top layer. Models predicting emission of dust by off-road driving should thus consider a number of soil parameters and not just one key parameter. Vehicle type and driving speed are additional parameters that affect emission. In general, 4-wheelers produce more dust than dune buggies, and dune buggies, more than dirt bikes. Higher speeds also result in higher emissions. Dust emitted by off-road driving is less coarse than the parent sediment on the road surface. Off-road driving thus results in a progressive coarsening of the top layer. Exceptions to this are silty surfaces with no, or almost no, vegetation. For such surfaces no substantial differences were observed between the grain size distribution of road dust and emitted dust. Typical emission values for off-road driving on dry desert soils are: for sandy areas, 30–40 g km^− 1 (PM10) and 150–250 g km^− 1 (TSP); for silty areas, 100–200 g km^− 1 (PM10) and 600–2000 g km^− 1 (TSP); for drainages, 30–40 g km^− 1 (PM10) and 100–400 g km^− 1 (TSP); and for mixed terrain, 60–100 g km^− 1 (PM10) and 300–800 g km^− 1 (TSP). These values are for the types of vehicles tested in this study and do not refer to cars or trucks, which produce significantly more dust.     

Time series assessment on landslide occurrences in an area undergoing development

This study assesses the influence of development on landslide occurrences in a rapidly developing area, Kota Kinabalu in Sabah Malaysia, across three assessment years (1978, 1994 and 2010). Two development indicators, land use and road density, were used to measure the influence of development on landslide occurrence. Land use was classified into four categories (barren, forest, developed and other), and road density was classified into low (< 50 m/40 000 m²), moderate (50–150 m/40 000 m²) and high (> 150 m/40 000 m²). Landslide density analysis was used to calculate the concentration of landslide for the different land use and road density categories. The number of landslides in developed areas increased from 19 landslides/100 km² in 1978, to 29 landslides/100 km² in 1994 and to 50 landslides/100 km² in 2010, mirroring an increase in land use for development purposes from 8 per cent in 1978 to 27 per cent in 2010. Landslide density also gradually increased in the high road density class from 10 landslides/100 km² in 1978, to 30 landslides/100 km² in 1994 and 62 landslides/100 km² in 2010. These results show that road construction activities influence landslide occurrences.     

Wave train model for knickpoint migration

Rivers respond to a drop in their base level by incising the topography. The upstream propagation of an incision, as usually depicted by a knickpoint migration, is thought to depend on several parameters such as the drainage area, lithology, and the amplitude of the base level drop. We first investigate the case of the Messinian Salinity Crisis that was characterized by the extreme base level fall (1500 m) of the Mediterranean Sea at the end of the Miocene. The response of drainage areas of three orders of magnitude (10³ to 10⁶ km²) highlights the dominant role of the drainage area (with a square root relationship) in controlling the knickpoint migration after a base level fall. A compilation of mean rates of knickpoint propagation for time durations ranging from 10² to 10⁷ years displays a similar relationship indicating that successive wave trains of knickpoint can migrate in a river: first, wave trains linked to the release of the alluvial cover and then, wave trains related to the bedrock incision, which correspond to the real time response of rivers. Wave trains with very low retreat rates (long lived knickpoints > 1 My) rather correspond to the response time of regional landscape.     

Susceptibility assessment of shallow landslides on Oahu, Hawaii, under extreme-rainfall events

The deterministic Stability INdex MAPping (SINMAP) model, which integrates a mechanistic infinite-slope stability model and a hydrological model, was applied to assess susceptibility of slopes in 32 shallow-landslide-prone watersheds of the eastern to southern areas of Oahu, Hawaii, USA. Input to the model includes a 10-m Digital Elevation Model (DEM), an inventory of storm-induced landslides that occurred from 1949 to 2006, and listings of soil-strength and hydrological parameters including transmissivity and steady-state recharge. The study area of ca. 384 km² was divided into four calibration regions with different geotechnical and hydrological characteristics. All parameter values were separately calibrated using observed landslides as references. The study used a quasi-dynamic scenario of soil wetness resulting from extreme daily rainfall events with a return period of 50 years. The return period was based on almost-90-year-long (1919–2007) daily rainfall records from 26 raingauge stations in the study area. Output of the SINMAP model includes slope-stability-index-distribution maps, slope-versus-specific-catchment-area charts, and statistical summaries for each region.The SINMAP model assessed susceptibility at the locations of all 226 observed shallow landslides and classified these susceptible areas as unstable. About 55% of the study area was predicted as highly unstable, highlighting a critical island problem. The SINMAP predictions were compared to an existing debris-flow-hazard map. Areas classified as unstable in the current study were classified as low-to-moderate and moderate-to-high debris-flow hazard risks by the prior mapping. The slope-stability maps provided by this study will aid in explaining the causes of known landslides, making emergency decisions, and, ultimately mitigating future landslide risks. The maps may be further improved by incorporating heterogeneous and anisotropic soil properties and spatial and temporal variation of rainfalls as well as by improving the accuracy of the DEM and the locations of shallow landslide initiation.     

Probabilistic landslide risk analysis considering direct costs in the area north of Lisbon (Portugal)

The purpose of the present study is the analysis of landslide risk for roads and buildings in a small test site (20 km²) in the area north of Lisbon (Portugal). For this purpose, an evaluation is performed integrating into a GIS information obtained from multiple sources: (i) landslide hazard; (ii) elements at risk; and (iii) vulnerability. Landslide hazard is assessed on a probabilistic basis for three different types of slope movement (shallow translational slides, translational slides and rotational slides), based on some assumptions such as: (i) the likelihood of future landslide occurrence can be measured through statistical relationships between past landslide distribution and specified spatial data sets considered as landslide predisposing factors; and (ii) the rainfall combination (amount–duration) responsible for past slope instability within the test site will produce the same effects (i.e. same type of landslides and similar total affected area), each time they occur in the future. When the return period of rainfall triggering events is known, different scenarios can be modelled, each one ascribed to a specific return period. Therefore, landslide hazard is quantitatively assessed on a raster basis, and is expressed as the probability for each pixel (25 m²) to be affected by a future landslide, considering a rainfall triggering scenario with a specific return period. Elements at risk within the test site include 2561 buildings and roads amounting to 169 km. Values attributed to elements at risk were defined considering reconstruction costs, following the guidelines of the Portuguese Insurance Institute. Vulnerability is considered as the degree of loss to a given element resulting from the occurrence of a landslide of a given magnitude. Vulnerability depends not only on structural properties of exposed elements, but also on the type of process, and its magnitude; i.e., vulnerability cannot be defined in absolute terms, but only with respect to a specific process (e.g. vulnerability to shallow translational slides). Therefore, vulnerability was classified for the three landslide groups considered on hazard assessment, taking into account: (i) landslide magnitude (mean depth, volume, velocity); (ii) damage levels produced by past landslide events in the study area; and (iii) literature. Finally, a landslide risk analysis considering direct costs was made in an automatic way crossing the following three layers: (i) Probabilistic hazard map for a landslide type Z, considering a particular rainfall triggering scenario whose return period is known; (ii) Vulnerability map (values from 0 to 1) of the exposed elements to landslide type Z; and (iii) Value map of the exposed elements, considering reconstruction costs.     

LIDAR monitoring of mass wasting processes: The Radicofani landslide, Province of Siena, Central Italy

The Radicofani Basin, stretching about 30 km NW–SE, is an intra-Central Apennine basin connected to Pliocene–Pleistocene extensional tectonics. It consists of an Early to Middle Pliocene succession including essentially shelf pelites. In the Radicofani area, province of Siena (Tuscany region), morphodynamic processes are very frequent with widespread badlands and rapidly evolving mudflows. In order to evaluate the general instability of the Radicofani area, geological and geomorphological surveys were carried out. The 1954, 1990 and 2003 aerial surveys allowed a comparison of the changes in the various morphological aspects of the study area, which suggested an increase in slope instability with time. A new complex translational landslide evolving into mudflows, activated during the winter of 2003, was monitored using an experimental system based on terrestrial LIDAR (Light Detection and Ranging) and GPS (Global Positioning System) technologies. This system allowed the monitoring of the morphologic and volumetric evolution of the landslide. A comparison of the monitoring data of October 2004, June 2005, May 2006 and May 2007 points out that the evolution is characterised by the sliding of displaced materials. A volume of about 1300 m³ of materials was removed during the period 2004–2005, 300 m³ for 2005–2006, and 400 m³ for 2006–2007. The greater initial mass movement probably reflects a greater static imbalance during the early period of landslide movement and increased rainfall. Therefore, the proposed monitoring system methodology allows the numerical evaluation of the landslide morphological evolution and to validate the landslide evolution model based on geological and geomorphological field surveys.     

Morphology and internal structure of a dormant landslide in a hilly area: The Collinabos landslide (Belgium)

This study attempts to reconstruct the history of the Collinabos landslide, a landslide with a fresh morphology that is representative for more than 150 dormant, deep-seated (> 3 m) landslides in the Flemish Ardennes (Belgium). A geomorphological map was created based on LIDAR (Light Detection and Ranging)-derived maps and detailed field surveys. The map showed that the landslide consisted of three zones with significant differences in surface topography. The northern landslide zone 1 is characterised by at least five reverse slopes, whereas zones 2 and 3, the southern landslide zones, have only two reverse slopes and a convex foot. Electric resistivity profiles measured in zones 1 and 2 revealed that the differences in surface topography were not related to differences in internal structure as both parts of the landslide were initiated as a rotational earth slide with a surface of rupture at 15 m deep, where the displaced material broke apart in two blocks. However, two shear surfaces of reactivations within landslide debris were only distinguished in the accumulation area of zone 1. The observed differences in surface morphology can be caused by a temporary conversion of a forest into cropland in zone 2. It is suggested that reverse slopes of smaller reactivations within landslide debris were obliterated during the agricultural activities. AMS radiocarbon dating of organic material found in ponds located in reverse slopes generally resulted in relatively recent dates (i.e. 1400–1950 Cal AD) suggesting that several of the small local reactivations occurred in that period. One dating at 8700–8440 Cal BP of organic matter collected in a reverse slope in zone 1 suggests that an initiation under periglacial conditions cannot be excluded for the Collinabos landslide. By combining different technologies, this study provides valuable information for a better understanding of dormant landslides.     

Kinematics of a landslide derived from archival photogrammetry and GPS data

The evolution of the Patigno landslide, a deep-seated gravitational slope deformation in the Northern Apennine range (Italy), was investigated using archival photogrammetry, and GPS observations from a permanent station located inside the landslide. Historical aerial photographs of the area taken in 1975 (scale 1:13,000), 1987 (1:13,000) and 2004 (1:30,000) were co-registered into the same reference frame using an unconventional method based on the detection of homologous points in multi-temporal models. Three DTMs were produced using a Digital Photogrammetric Workstation and compared. The displacement vectors of 293 points in the landslide were determined. The average movement velocity of the GPS station since 2004 (about 3.5 cm yr^− 1) agrees with the mean displacement rate obtained from photogrammetry. Furthermore, the estimated velocity varies seasonally between 3 and 6 cm yr^− 1. This variability correlated with rainfall suggests that the creep of the landslide is influenced by hydrological factors.     

LiDAR-derived DEM evaluation of deep-seated landslides in a steep and rocky region of Japan

In steep and rocky terrains, their rough surfaces make it difficult to create landslide inventories even with detailed maps/images produced from airborne LiDAR data. To provide objective clues in locating deep-seated landslides, the surface textures of a 5 km² steepland area in Japan was investigated using the eigenvalue ratio and slope filters calculated from a very high resolution LiDAR-derived DEM. The range of filter values was determined for each of a number of surface features mapped in the field and these included: cracked bedrock outcrops, coarse colluvial deposits, gently undulating surfaces, and smooth surfaces. Recently active slides commonly contained patches of ground in which deposition and erosion occurred together near the erosion front, or where cracked bedrock outcrops and coarse colluvial deposits coexisted under a gently undulating surface. The characteristic eigenvalue and slope filter values representing this sliding process were applied to maps of the DEM derived filter values to extract potential sites of recent landslide activity. In addition, the relationships between the filter values of deep-seated landslides at various stages of evolution within the field mapped area were extended to the entire study area, to assess the contribution that landslide evolution makes to change in the landscape as a whole. While landslide components made up the steepest as well as the gentlest parts of the landscape depending on their evolutionary stage, landslides were constantly coarsened and steepened by progressive erosion, probably initiated by river bank erosion at the foot of slopes.