The effects of logging on frequency and distribution of landslides in three watersheds on Vancouver Island, British Columbia

Three hundred and sixty three landslides in three watersheds that totaled 382 km² were identified from air photographs, beginning at a date that preceded logging to the present. The three watersheds all lie on Vancouver Island; however, they have different precipitation regimes, topography, and amounts logged. Landslide areas in the watersheds varied in size from 200 m² to more than 1 ha. Nearly 80% of the landslides were debris slides; 15% were debris flows, and the remainder primarily rock falls. Following logging, the number of landslides increased substantially in all watersheds although the amount of increase was variable: approximately 11, 3, and 16 times in Macktush Creek, Artlish River, and Nahwitti River, respectively. Other analyses of changes in landslide density also produced highly variable results, with the number of landslides increasing between 2.4 and 24 times. Further, 2–12 times more landslides reached streams following logging activities. Densities for landslides impacting streams increased for the period of record from 1.5 to 10 times following logging activities. The densities were substantially greater where only landslides that reached streams since development began in a watershed were considered. Roads had the greatest spatial impact in the watersheds (compared to their total area), with frequencies determined to have increased by 27, 12, and 94 times for Macktush, Artlish, and Nahwitti, respectively. The results highlight the relative impact of roads and their role in slope stability.     

The frequency and distribution of recent landslides in three montane tropical regions of Puerto Rico

Landslides are common in steep mountainous areas of Puerto Rico where mean annual rainfall and the frequency of intense storms are high. Each year, landslides cause extensive damage to property and occasionally result in loss of life. Average population density is high, 422 people/km², and is increasing. This increase in population density is accompanied by growing stress on the natural environment and physical infrastructure. As a result, human populations are more vulnerable to landslide hazards. The Blanco, Cibuco, and Coamo study areas range in surface area from 276 to 350 km² and represent the climatologic, geographic, and geologic conditions that typify Puerto Rico. Maps of recent landslides developed from 1:20,000-scale aerial photographs, in combination with a computerized geographic information system, were used to evaluate the frequency and distribution of shallow landslides in these areas. Several types of landslides were documented—rainfall-triggered debris flows, shallow soil slips, and slumps were most abundant. Hillslopes in the study area that have been anthropogenically modified, exceed 12° in gradient, are greater than 300 m in elevation, and face the east-northeast, are most prone to landsliding. A set of simplified matrices representing geographic conditions in the three study areas was developed and provides a basis for the estimation of the spatial controls on the frequency of landslides in Puerto Rico. This approach is an example of an analysis of the frequency of landslides that is computationally simple, and therefore, may be easily transferable to other settings.     

Size and spatial distribution of loess slides on the Chinese Loess Plateau

ABSTRACT

The size and spatial distribution of loess slides are important for estimating the yield of eroded materials and determining the landslide risk. While previous studies have investigated landslide size distributions, the spatial distribution pattern of landslides at different spatial scales is poorly understood. The results indicate that the loess slide distribution exhibits a power-law scaling across a range of the size distribution. The mean landslide size and size distribution in the different geomorphic types are different. The double Pareto and inverse gamma functions can coincide well with the empirical probability distribution of the loess slide areas and can quantitatively reveal the rollover location, maximum probability, and scaling exponents. The frequency of loess slides increases with mean monthly precipitation. Moreover, point distance analysis showed that > 80% of landslides are located < 3 km from other loess slides. We found that the loess slides at the two study sites (Zhidan and Luochuan County) in northern Shaanxi Province, China show a significant clustered distribution. Furthermore, analysis results of the correlated fractal dimension show that the landslides exhibit a dispersed distribution at smaller spatial scales and a clustered distribution at larger spatial scales.     

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.     

Magnitude–frequency relationships of debris flows and debris avalanches in relation to slope relief

Probability of occurrence, hazard intensity and encounter probability are key parameters in the quantitative risk analysis (QRA) of landslides. All are strongly dependent on magnitude of the landslides. As a result, magnitude–frequency analysis should be a part of QRA. Deriving representative magnitude–frequency relationships for debris avalanches and debris flows, however, is difficult. One key problem is illustrated with the example of a unique database from the coastal region of British Columbia, Canada, which was compiled entirely from detailed ground investigations. The magnitude of debris avalanches and debris flows is not an independent statistical quantity, but a function of the scale of a given slope, as characterized by the slope length. Thus, attempting to derive probability and magnitude for a given location or sub-region from a regionally-derived magnitude–frequency curve may lead to incorrect predictions. The same problem is pertinent to the application of the same approach to any type of landslide in which the largest combined dimension of the source volume (including entrainment) is of the same order as the length of the slope. It is recommended that greater emphasis be placed on site-specific geological observations, at the expense of generalized statistics.     

GIS analysis to assess landslide susceptibility in a fluvial basin of NW Sicily (Italy)

Landslide hazard assessment, effected by means of geostatistical methods, is based on the analysis of the relationships between landslides and the spatial distributions of some instability factors. Frequently such analyses are based on landslide inventories in which each record represents the entire unstable area and is managed as a single instability landform. In this research, landslide susceptibility is evaluated through the study of a variety of instability landforms: landslides, scarps and areas uphill from crown. The instability factors selected were: bedrock lithology, steepness, topographic wetness index and stream power index. The instability landform densities computed for all the factors, which were arranged in Unique Condition Unit, allowed us to derive a total of three prediction images for each landslide typology. The role of the instability factors and the effects generated by the use of different landforms were analyzed by means of: a) bivariate analysis of the relationships between factors and landslide density; b) predictive power validations of the prediction images, based on a random partition strategy.The test area was the Iato River Basin (North-Western Sicily), whose slopes are moderately involved in flow and rotational slide landslides (219 and 28, respectively). The area is mainly made up of the following complexes: Numidian Flysch clays (19%, 1%), Terravecchia sandy clays (5%, 1%), Terravecchia clayey sands (3%, 0.3%) and San Cipirello marly clays (9%, 0%). The steepness parameter shows the highest landslide density in the [11–19°] class for both the typologies (8%, 1%), even if the density distributions for rotational slides are right-asymmetric and right-shifted. We obtained significant differences in shape when we used different instability landforms. Unlike scarps and areas uphill from crowns, landslide areas produce left-asymmetric and left-shifted density distributions for both the typologies. As far as the topographic wetness index is concerned, much more pronounced differences were detected among the instability landforms of rotational slides. In contrast, the flow landslides produce normal-like density distributions. The latter and the rotational slide landslide areas produce the highest density values in the class [5.5–6.7], despite an abrupt decreasing trend starting from the first class [3.2–4.4], which is generated by the density values of the rotational slide scarps and areas uphill from crowns. The stream power index at the foot of the slopes, which was automatically derived using a GIS-procedure, shows a positive correlation with the landslide densities marked by the maximum classes: [4.8–6.0] for flows, and [6.0–7.2] for rotational slides. The validation procedure results confirmed that the choice of instability landform influences the results of the susceptibility analysis. Furthermore, the validation procedure indicates that: a) the predictive models are generally satisfactory; b) scarps and zones uphill from crown areas are the most diagnostically unstable landforms, for flow and rotational slide landslides respectively.     

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.     

Landslide inventory in a rugged forested watershed: a comparison between air-photo and field survey data

Landslide inventories are routinely compiled by means of aerial photo interpretation (API). When examining photo pairs, the forest canopy (notably in old-growth forest) hides a population of “not visible” landslides. In the present study, we attempt to estimate how important is the contribution of landslides not detectable from aerial photographs to the global mass of sediment production from mass failures on forested terrain of the Capilano basin, coastal British Columbia. API was coupled with intensive fieldwork for identification and measurement of all landslides. A 30-year framework was adopted. We show that “not visible” landslides can represent up to 85% of the total number of failures and account for 30% of the volume of debris mobilised. Such percentages display high sub-basin variability with rates of sediment production varying by one order of magnitude between two sub-basins of the study area. This is explained qualitatively by GIS-based analysis of slope frequency distributions, drainage density, and spatial distribution of surficial materials. Such observations find further support in the definitions of transport-limited and supply-limited basins. As a practical consideration to land managers, we envisage that supplementary fieldwork for landslide identification is mandatory in transport-limited systems only. Fieldwork has demonstrated that gully-related failures have a greater importance than one could expect from API.     

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.     

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.     

A diatom-based reconstruction of drought intensity, duration, and frequency from Moon Lake, North Dakota: a sub-decadal record of the last 2300 years

Diatom assemblages preserved in sediment cores from closed-basin lakes can provide high-resolution records of past hydrologic and climatic conditions, including long-term patterns in the intensity, duration, and frequency of droughts. At Moon Lake, a closed-basin lake in eastern North Dakota, a comparison of diatom-inferred salinity and the precipitation-based Bhalme-Mooley Drought Index (BMDI) over the last 100 years was highly significant, suggesting that the diatom record contains a sensitive archive of past climatic conditions. A sub-decadal record of inferred salinity for the past 2300 years indicates that extreme droughts of greater intensity than those during the 1930s 'Dust Bowl' were more frequent prior to A.D. 1200. This high frequency of extreme droughts persisted for centuries and was most pronounced from A.D. 200–370, A.D. 700–850 and A.D. 1000–1200. A pronounced shift to generally wetter conditions with less severe droughts of shorter duration occured at A.D. 1200. This abrupt change coincided with the end of the 'Medieval Warm Period' (A.D. 1000–1200) and the onset of the 'Little Ice Age' (A.D. 1300–1850).     

Short-term changes in the magnitude, frequency and temporal distribution of floods in the Eastern Mediterranean region during the last 45 years — Nahal Oren, Mt. Carmel, Israel

Short-term changes in Eastern Mediterranean precipitation affecting flow regime were documented in Nahal Oren, a 35 km² ephemeral stream in Mt. Carmel, a 500 m high mountain ridge located at the NW coast of Israel. The rainy winter of the Mediterranean type climate (Csa) in Mt. Carmel is characterized by average annual rainfall of 550 mm at the coastal plain to 750 mm at the highest elevation while the summer is hot and dry. Stream flow generates after accumulated rainfall of 120–150 mm while “large floods”, defined as flows with peak discharge of > 5 m³ s^− 1 and/or > 150,000 m³ in volume, are generated in response to rainfall of over 100 mm. Hence, large floods in Nahal Oren stream occur not earlier than December. Precipitation and flow data were divided into two sub-periods: 1957–1969 and 1991–2003 and compared to each other. The results indicate a clear increase in the frequency of large floods, their magnitudes and volumes during the second period with no parallel change in the annual precipitation. Similarly, an increase in storm rainfall–runoff ratio from < 5% to > 15% and a decrease in the threshold rainfall for channel flow by 16–25% were documented. These short-term variations in flooding behavior are explained by the clear decrease in the length of the rainy season and by the resulting significant shortening in the duration of the dry-spells. The increase in the number of large rainfall events and the large floods in each hydrological year together with the increasing number of years with no floods indicate strengthening of their uncertainty of behavior.     

On the spatial relationship between landslides and causative factors on Lantau Island, Hong Kong

This paper presents a statistical approach to study the spatial relationship between landslides and their causative factors at the regional level. The approach is based on digital databases, and incorporates such methods as statistics, spatial pattern analysis, and interactive mapping. Firstly, the authors propose an object-oriented conceptual model for describing a landslide event, and a combined database of landslides and environmental factors is constructed by integrating the various databases within such a conceptual framework. The statistical histogram, spatial overlay, and dynamic mapping methods are linked together to interactively evaluate the spatial pattern of the relationship between landslides and their causative factors. A case study of an extreme event in 1993 on Lantau Island indicates that rainfall intensity and the migration of the center of the rainstorm greatly influence the occurrence of landslides on Lantau Island. A regional difference in the relationship between landslides and topography is identified. Most of the landslides in the middle and western parts of the island occurred on slopes with slope angles of 25–35°, while in the eastern part, the corresponding range is 30–35°. Overlaying landslide data with land cover reveals that a large number of landslides occurred in the bareland and shrub-covered area, and in the transition zones between different vegetation types. The proposed approach can be used not only to analyze the general characteristics of such a relationship, but also to depict its spatial distribution and variation, thereby providing a sound basis for regional landslide prediction.     

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.     

The statistical relationship between unconfined compressive strengths and the frequency distributions of slope gradients — A case study in northern Hungary

This paper focuses on the question of whether there is a deterministic connection between the slope gradient and unconfined compression strength (UCS) as a lithological factor and on describing the nature of this connection. Moreover, the authors determined the sensitive statistical parameters in the statistical surface analyses. The surface analyses were carried out in an area extending over 1500 km² in northern Hungary and containing parts of the uplifted Palaeozoic and Mesozoic basement and the semi-consolidated material of the Palaeogene and Neogene molasse sediments.The 67 geological formations of the area were grouped into 10 petrophysical categories characterised by unconfined compressive strength as a petrophysical parameter. The geological database was the digitalized geological map of North Hungarian Region (1 : 100 000). The digital topographic database was based on 10 m contour lines of 1 : 50 000 maps, the digital elevation model was generated by kriging interpolation. Three topographic models were created with resolutions of 25 × 25, 50 × 50 and 100 × 100 m per pixels.Evident correlation can be shown between the UCS and the relative frequency of the slope gradient. The adequate regression procedure is the power regression for low slope category values while logarithmic regression is applicable at high slope angles. Based on the characteristic of the relationship and the value of r², slope category intervals can be identified the relative frequency of which is proven to be determined by the UCS. These intervals are found to be between 4–10%, 10–16%, 16–22%, 22–44% and over 44%.Using the determination equations of slope gradient between 4–10% and over 44%, the UCS of the bedrock can be calculated approximately as the average value of the two calculated results. So the quotient of the frequency of these two category intervals can be regarded as an important morphometric index for a given bedrock.From the aspect of petrophysical characteristics, the rock with UCS between 6 and 86 MPa were proven to be deterministic for the slope development i.e. in the determination of the slope category frequencies. Applying the relative frequency of slope gradients a relative erosion resistance of the petrophysical categories can be calculated.Considering the determination coefficients, among the statistical parameters of the distribution of slope category values the standard deviation, mode, mean and median proved to be determined by the UCS of the bedrock.     

An analytical approximation of the distribution of total hydrocarbon discoveries as a function of discovery number

Recently, Manly's method has been successfully applied to hydrocarbon exploration modeling in order to approximate the expected value and the standard deviation of the total amount of hydrocarbons discovered. This method is much faster than running prolonged simulations normally required by the probabilistic model of the hydrocarbon discovery process, and the results are very accurate. This paper extends the usefulness of the approximation method by developing an approximate analytical model of the whole probability distribution of the total volume of hydrocarbons discovered. The mean and the standard deviation from Manly's approximation are used to help set the parameters of a family of beta distributions, to represent the distributions of the total amount of hydrocarbons discovered from the beginning to the end of the exploration process in an area. Three real datasets—the Nova Scotian Shelf from offshore northeastern Canada, the Bistcho Play, and the Zama Play from northwestern Canada—are chosen to verify the methodology developed. Confidence intervals of the forecast for each number of discovered fields are constructed from the analytical approximation and compared with confidence intervals generated by the simulation. Sensitivity analyses are performed to show that the idea of using a family of beta distributions is a robust approximation.     

Slope stability and landslides in proglacial varved clays of western Estonia

During the last decade the frequency of landslides at river valley slopes eroding into the glaciolacustrine plain in western Estonia has grown considerably. We studied in detail nine recent landslides out of 25 known and recorded sliding events in the area. All landslides occurred at the river banks in otherwise almost entirely flat areas of proglacial deposits capped with marine sands. Glaciolacustrine varved clay is the weakest soil type in the area and holds the largest landslides. Slope stability modelling shows that critical slope gradient for the clay is ≥ 10° and for the marine sand ≥ 20°. Fluvial erosion is the main process in decreasing slope stability at the outer bends of the river meanders. An extra shear stress generated by groundwater flow following the high stand of the groundwater level or rapid water level drawdown in the river channels are responsible for triggering the landslides. Consecutive occurrence of small-scale slides has a direct effect in triggering the large, retrogressive complexes of slides in the glaciolacustrine clay. A landslide hazard zonation map was composed based on digital elevation model and the data on spatial distribution of glaciolacustrine clays and marine sands, and on existing and critical slope angles of these deposits.     

Hillslope-to-valley transition morphology: New opportunities from high resolution DTMs

The search for the optimal spatial scale for observing landforms to understand physical processes is a fundamental issue in geomorphology. Topographic attributes derived from Digital Terrain Models (DTMs) such as slope, curvature and drainage area provide a basis for topographic analyses. The slope–area relationship has been used to distinguish diffusive (hillslope) from linear (valley) processes, and to infer dominant sediment transport processes. In addition, curvature is also useful in distinguishing the dominant landform process. Recent topographic survey techniques such as LiDAR have permitted detailed topographic analysis by providing high-quality DTMs. This study uses LiDAR-derived DTMs with a spatial scale between 1 and 30 m in order to find the optimal scale for observation of dominant landform processes in a headwater basin in the eastern Italian Alps where shallow landsliding and debris flows are dominant. The analysis considered the scaling regimes of local slope versus drainage area, the spatial distribution of curvature, and field observations of channel head locations. The results indicate that: i) hillslope-to-valley transitions in slope–area diagrams become clearer as the DTM grid size decreases due to the better representation of hillslope morphology, and the topographic signature of valley incision by debris flows and landslides is also best displayed with finer DTMs; ii) regarding the channel head distribution in the slope–area diagrams, the scaling regimes of local slope versus drainage area obtained with grid sizes of 1, 3, and 5 m are more consistent with field data; and iii) the use of thresholds of standard deviation of curvature, particularly at the finest grid size, were proven as a useful and objective methodology for recognizing hollows and related channel heads.     

GIS-based landslide susceptibility mapping for the 2005 Kashmir earthquake region

The M_w 7.6 October 8, 2005 Kashmir earthquake triggered several thousand landslides throughout the Himalaya of northern Pakistan and India. These were concentrated in six different geomorphic–geologic–anthropogenic settings. A spatial database, which included 2252 landslides, was developed and analyzed using ASTER satellite imagery and geographical information system (GIS) technology. A multi-criterion evaluation was applied to determine the significance of event-controlling parameters in triggering the landslides. The parameters included lithology, faults, slope gradient, slope aspect, elevation, land cover, rivers and roads. The results showed four classes of landslide susceptibility. Furthermore, they indicated that lithology had the strongest influence on landsliding, particularly when the rock is highly fractured, such as in shale, slate, clastic sediments, and limestone and dolomite. Moreover, the proximity of the landslides to faults, rivers, and roads was also an important factor in helping to initiate failures. In addition, landslides occurred particularly in moderate elevations on south facing slopes. Shrub land, grassland, and also agricultural land were highly susceptible to failures, while forested slopes had few landslides. One-third of the study area was highly or very highly susceptible to future landsliding and requires immediate mitigation action. The rest of the region had a low or moderate susceptibility to landsliding and remains relatively stable. This study supports the view that (1) earthquake-triggered landslides are concentrated in specific zones associated with event-controlling parameters; and (2) in the western Himalaya deforestation and road construction contributed significantly to landsliding during and shortly after earthquakes.     

Magnitude and frequency of landsliding in a large New Zealand catchment

Knowledge of long-term average rates of erosion is necessary if factors affecting sediment yields from catchments are to be understood. Without such information, it is not possible to assess the potential influence of extreme storms, and, therefore, to evaluate the relative importance of various components of a sediment budget. A study of the sediment budget for the Waipaoa catchment, North Island, New Zealand, included evaluation of long-term rates of landsliding for six landslide-prone land systems in the catchment. The number of landslides per unit area generated by each of several storms was counted on sequential aerial photographs and correlated with the magnitude of the corresponding storm. The resulting relationships were combined with magnitude–frequency relationships derived for storms from 70- to 100-year rainfall records in the area to estimate a long-term magnitude–frequency relationship for landsliding for each land system. The long-term average values of the areal landslide frequency (number of slides per unit area per unit time) were then calculated from these relationships. The volumes of a sample of landslide scars were measured in the field, and the proportion of slides that deliver sediment to channels was determined from aerial photographs. These measurements then allowed calculation of the long-term average rate of sediment production to streams from landslides for different land systems and types of vegetation. Results suggest that shallow landslides currently contribute about 15±5% of the suspended sediment load in the Waipaoa River above the Kanakanaia gauging station, and that 75% of the sediment production from the landslides occurs during storms with recurrence intervals of less than 27 years. Reforestation of 6.3% (93 km²) of the slide-prone lands in the catchment between 1990 and 1995 resulted in a calculated decrease in slide-derived sediment of 10%. Calculations suggest that reforestation of an additional 3% (66 km²) of the catchment in areas with the most sensitive combinations of land system and storm regime could decrease the total sediment inputs from landsliding by about 20%.