Large landslides and their effect on sediment flux in South Westland,New Zealand

Landslides and runoff are dominant erosional agents in the tectonically active alpine South Westland area of New Zealand, characterized by high uplift rates and extreme orographic precipitation. Despite a high density of shallow debris slides and flows, the geomorphic imprints of deep‐seated bedrock failures are dominant and persistent. Over 50 large (>1 km²) landslides comprising rock slide[sol ]avalanches, complex rotational and rock‐block slides, wedge failures, and deep‐seated gravitational slope deformation were detected on air photos and shaded‐relief images. Major long‐term impacts on alpine rivers include (1) forced alluviation upstream of landslide dams, (2) occlusion of gorges and triggering of secondary riparian landslides, and (3) diversion of channels around deposits to form incised meandering gorges. Remnants of large prehistoric (i.e. pre‐1840) landslide deposits possibly represent the low‐frequency (in terms of total area affected yet dominant) end of the spectrum of mass wasting in the western Southern Alps. This is at odds with high erosion rates in an active erosional landscape. Large landslides appear to have dual roles of supplying and retaining sediment. The implications of these roles are that (1) previous models of (shallow) landslide‐derived sediment flux need to be recalibrated, and (2) geomorphic effects of earthquake‐induced landsliding may persist for at least 10² years. Copyright © 2005 John Wiley & Sons, Ltd.     

Large ice‐contact slope movements: glacial buttressing,deformation and erosion

Glaciers and slope movements may act simultaneously to erode and modify glaciated slopes. Undercutting by glaciers can destabilize slopes but the extent to which slope failure may progress prior to subsequent glacier withdrawal has not hitherto been considered. The traditional view has been that the buttressing effect of ice prevents slope movement. The problem with this view is that ice is one‐third the density of rock and flows under low applied stress. Consequently, failed slopes may move into the glacier if they exert a stress in excess of the resistance provided by the glacier. Slope movement rate depends on ice rheology and other factors influencing driving and resisting stresses. Simple viscous equations are used to investigate these variables. The equations predict that small (<125 000 m³) ice‐contact rockslides can deform ice at several mm/year, increasing to several m/year for very large (>10⁸ m³) rockslides. To test these estimates, field evidence is presented of slope movements in glaciated valleys of New Zealand; narrowing or squeezing of glaciers adjacent to unstable rock slopes is demonstrated and considered to be the result of slope movement. For one site, geomorphic mapping and slope movement monitoring data show that movement rates are of similar order of magnitude to those predicted by the viscous equations; closer agreement could be achieved with the application of modelling techniques that can more realistically model the complex slope geometries and stability factors encountered, or by obtaining additional empirical data to calibrate the models. This research implies that, while the concept of glacial debuttressing – the reduction of slope support from withdrawal of glaciers – is valid, complete debuttressing is not a prerequisite for the movement of ice‐contact rock slopes. These slope movements may contribute to the erosional processes of glaciers and the evolution of glaciated slopes in a previously unrecognized way. Copyright © 2012 John Wiley & Sons, Ltd.     

Rain‐induced deep‐seated catastrophic rockslides controlled by a thrust fault and river incision in an accretionary complex in the Shimanto Belt,Japan

To clarify the geological causes of rockslides induced by rainstorms in accretionary complexes, the geology and geomorphology of two large rockslides (volumes > 10⁶ m³) induced by the heavy rainfall of Typhoon Talas in the Shimanto Belt, Kii Mountains, Japan in 2011 are investigated. Our analysis reveals that thrusts with brittle crush zones controlled the occurrence of the rockslides. The properties and distribution of thrusts were poorly constrained before this study. Flooding during the rainstorm removed surface materials along rivers, allowing thorough geological mapping to be performed. Gravitationally deformed slopes were studied using GIS analysis of 1 m digital elevation models (DEMs) and fieldwork, and X‐ray diffraction (XRD) analysis, permeability, and direct shear tests were used to characterize the mineralogy and geotechnical properties of fault gouge. The Kawarabi thrust has a brittle crush zone up to 6 m thick and acts as the sliding surface for both landslides. The thrust dips 34° downslope and is cut by high‐angle faults and joints along one or both sides of each landslide body. Prior to failure, the upper part of the slope contained small scarps, suggesting that the slopes were already gravitationally deformed. The slope instability can be attributed to long‐term river erosion, which has undercut the slope and exposed the thrust at the base of the slope. The groundwater level, monitored in boreholes, suggests that the Kawarabi thrust is a barrier to groundwater flow. The weak and impermeable nature of the thrust played an essential role in the generation of gravitational slope deformation and catastrophic failure during periods of increased rainfall. Thrusts are a common feature of accretionary complexes, including in the Shimanto Belt, and the mechanism of slope failure stated above can be typical of rockslides in accretionary complexes and provide new insights into landslide disaster mitigation.     

Large scale landslides triggered by Quaternary tectonics in the Acambay graben,Mexico

Detailed analysis was conducted on large‐scale gravitational‐tectonic deformations and landslides in the Acambay graben, an intra‐arc basin in the trans‐Mexican volcanic belt (TMVB). Field mapping and remote sensing revealed the slope instability of the northern graben boundary induced by the Acambay‐Tixmadejé fault. Two major landslides of 0·1 km³ and 0·05 km³ in volume were identified and their characteristics were analyzed according to the role of tectonics, mechanism of slope failure, and possible triggering factors. Quaternary faulting played a major role in increasing the local relief, and the activity of the Acambay‐Tixmadejé fault represents the main geomorphic factor conditioning the gravitational movements. Moreover, displacements along this fault generated sliding surfaces and reduced the strength of the rock mass. The two landslides are classified as large‐scale rotational slides involving volcanic rocks of late Miocene‐Pleistocene age. Since the Acambay graben is a seismogenic area with a known maximum horizontal ground acceleration of 0·5 g, a strong earthquake could be ascribed as the possible triggering mechanism of the landslides. Our work represents the first analysis of large gravitational slope movements in tectonically active regions in Mexico, a process that can be common in the intra‐arc basins of the TMVB, where active tectonic, seismicity, weak altered volcanic rocks, and heavy rains affect the slope stability. Copyright © 2010 John Wiley & Sons, Ltd.     

An integrated hydrological model for rain‐induced landslide prediction

This paper describes an extension to the Combined Hydrology And Stability Model (CHASM) to fully include the effects of vegetation and slope plan topography on slope stability. The resultant physically based numerical model is designed to be applied to site‐specific slopes in which a detailed assessment of unsaturated and saturated hydrology is required in relation to vegetation, topography and slope stability. Applications are made to the Hawke's Bay region in New Zealand where shallow‐seated instability is strongly associated with spatial and temporal trends in vegetation cover types, and the Mid‐Levels region in Hong Kong, an area subject to a variety of landslide mechanisms, some of which may be subject to strong topographic control. An improved understanding of process mechanism, afforded by the model, is critical for reliable and appropriate design of slope stabilization and remedial measures. Copyright © 2002 John Wiley & Sons, Ltd.     

Modeling past and future alpine permafrost distribution in the Colorado Front Range

Rock glaciers, a feature associated with at least discontinuous permafrost, provide important topoclimatic information. Active and inactive rock glaciers can be used to model current permafrost distribution. Relict rock glacier locations provide paleoclimatic information to infer past conditions. Future warmer climates could cause permafrost zones to shrink and initiate slope instability hazards such as debris flows or rockslides, thus modeling change remains imperative. This research examines potential past and future permafrost distribution in the Colorado Front Range by calibrating an existing permafrost model using a standard adiabatic rate for mountains (0·5 °C per 100 m) for a 4 °C range of cooler and warmer temperatures. According to the model, permafrost currently covers about 12 per cent (326·1 km²) of the entire study area (2721·5 km²). In a 4 °C cooler climate 73·7 per cent (2004·4 km²) of the study area could be covered by permafrost, whereas in a 4°C warmer climate almost no permafrost would be found. Permafrost would be reduced severely by 93·9 per cent (a loss of 306·2 km²) in a 2·0 °C warmer climate; however, permafrost will likely respond slowly to change. Relict rock glacier distribution indicates that mean annual air temperature (MAAT) was once at least some 3·0 to 4·0 °C cooler during the Pleistocene, with permafrost extending some 600–700 m lower than today. The model is effective at identifying temperature sensitive areas for future monitoring; however, other feedback mechanisms such as precipitation are neglected. Copyright © 2005 John Wiley & Sons, Ltd.     

Hillslope response to tectonic forcing in threshold landscapes

Hillslopes are thought to poorly record tectonic signals in threshold landscapes. Numerous previous studies of steep landscapes suggest that large changes in long‐term erosion rate lead to little change in mean hillslope angle, measured at coarse resolution. New LiDAR‐derived topography data enables a finer examination of threshold hillslopes. Here we quantify hillslope response to tectonic forcing in a threshold landscape. To do so, we use an extensive cosmogenic beryllium‐10 (¹⁰Be)‐based dataset of catchment‐averaged erosion rates combined with a 500 km² LiDAR‐derived 1 m digital elevation model to exploit a gradient of tectonic forcing and topographic relief in the San Gabriel Mountains, California. We also calibrate a new method of quantifying rock exposure from LiDAR‐derived slope measurements using high‐resolution panoramic photographs. Two distinct trends in hillslope behavior emerge: below catchment‐mean slopes of 30°, modal slopes increase with mean slopes, slope distribution skewness decreases with increasing mean slope, and bedrock exposure is limited; above mean slopes of 30°, our rock exposure index increases strongly with mean slope, and the prevalence of angle‐of‐repose debris wedges keeps modal slopes near 37°, resulting in a positive relationship between slope distribution skewness and mean slope. We find that both mean slopes and rock exposure increase with erosion rate up to 1 mm/a, in contrast to previous work based on coarser topographic data. We also find that as erosion rates increase, the extent of the fluvial network decreases, while colluvial channels extend downstream, keeping the total drainage density similar across the range. Our results reveal important textural details lost in 10 or 30 m resolution digital elevation models of steep landscapes, and highlight the need for process‐based studies of threshold hillslopes and colluvial channels. Copyright © 2012 John Wiley & Sons, Ltd.     

A chronology of post‐glacial landslides suggests that slight increases in precipitation could trigger a disproportionate geomorphic response

Large, deep‐seated landslides are common throughout the south‐eastern San Juan Mountains of Colorado and New Mexico, but their timing and initiation are not well understood. Determining when the landslides occurred would aid in clarifying the mechanisms for initiating landslides in the region and would help us to understand post‐glacial landscape evolution. We studied seven pre‐historic landslides located within the Tertiary volcanic rocks of the San Juan Volcanic Field. The landslides range in area from ~0.8 km² to ~11.3 km² and most are located in areas that were previously mapped as having been ice‐covered during the last glaciation. Landslide deposits were dated using a variety of methods including surface exposure dating (chlorine‐36, ³⁶Cl), radiocarbon dating of basal bog sediments and organic material buried in soils, and relative soil development. The resulting limiting ages range from approximately 14 ka to 2 ka and show that deep‐seated landsliding has occurred throughout the post‐glacial period. This broad range in ages is inconsistent with our initial hypothesis, which proposed that landslides were likely the result of debuttressing of glacial walls during glacial retreat. Furthermore, the timing of landslides does not seem to correlate with documented post‐glacial climatic shifts. Therefore, we conclude that landsliding in the region was the result of wetter than normal periods lasting months to years acting on weak bedrock preconditioned to failure and prepared by glacial debuttressing. Our findings suggest that the study area is likely still susceptible to deep‐seated landsliding and may become even more prone to large‐scale slope failure if future climate change increases precipitation in the San Juan Mountains. Copyright © 2017 John Wiley & Sons, Ltd.     

Slope control on the frequency distribution of shallow landslides and associated soil properties,North Island,New Zealand

The extrapolation of results from field trials to larger areas of land for purposes of regional impact assessment is an important issue in geomorphology, particularly for landform properties that show high stochastic variability in space and time, such as shallow landslide erosion. It is shown in this study, that by identifying the main driver for spatial variability in shallow landslide erosion at field scales, namely slope angle, it is possible to develop a set of generic functions for assessing the impact of landslides on selected soil properties at larger spatial scales and over longer time periods. Research was conducted within an area of pastoral soft‐rock Tertiary hill country in the North Island of New Zealand that is subject to infrequent high intensity rainfall events, producing numerous landslides, most of which are smaller than several hundred square metres in size and remove soil to shallow depths. All landslides were mapped within a 0·6 km² area and registered to a high resolution (2 m) slope map to show that few landslides occur on slopes < 20° and 95% were on slopes > 24°. The areal density of landslides from all historical events showed an approximately linear increase with slope above 24°. Integrating landslide densities with soil recovery data demonstrates that the average value of a soil property fluctuates in a ‘saw‐tooth’ fashion through time with the overall shape of the curve controlled by the frequency of landslide inducing storm events and recovery rate of the soil property between events. Despite such fluctuations, there are gradual declines of 7·5% in average total carbon content of topsoil and 9·5% in average soil depth to bedrock, since the time of forest clearance. Results have application to large‐scale sediment budget and water quality models and to the New Zealand Soil Carbon Monitoring System (CMS). Copyright © 2012 John Wiley & Sons, Ltd.     

Linking weathering and rock slope instability: non‐linear perspectives

Weathering is linked complexly to the erosion and evolution of rock slopes. Weathering influences both the strength of rock slopes and the stresses that act upon them. While weathering has often been portrayed in an over‐simplified way by those studying rock slope instability, in reality it consists of multiple processes, acting over different spatial and temporal scales, with many complex inter‐linkages. Through a demonstration of the sources of non‐linearities in rock slope weathering systems and their implications for rock slope instability, this paper proposes five key linkages worthy of further study. Copyright © 2012 John Wiley & Sons, Ltd.     

Erosional power in the Swiss Alps: characterization of slope failure in the Illgraben

Landslides and rockfalls are key geomorphic processes in mountain basins. Their quantification and characterization are critical for understanding the processes of slope failure and their contributions to erosion and landscape evolution. We used digital photogrammetry to produce a multi‐temporal record of erosion (1963–2005) of a rock slope at the head of the Illgraben, a very active catchment prone to debris flows in Switzerland. Slope failures affect 70% of the study slope and erode the slope at an average rate of 0.39 ± 0.03 m yr¯¹. The analysis of individual slope failures yielded an inventory of ~2500 failures ranging over 6 orders of magnitude in volume, despite the small slope area and short study period. The slope failures form a characteristic magnitude–frequency distribution with a rollover and a power‐law tail between ~200 m³ and 1.6 × 10⁶ m³ with an exponent of 1.65. Slope failure volume scales with area as a power law with an exponent of 1.1. Both values are low for studies of bedrock landslides and rockfall and result from the highly fractured and weathered state of the quartzitic bedrock. Our data suggest that the magnitude–frequency distribution is the result of two separate slope failure processes. Type (1) failures are frequent, small slides and slumps within the weathered layer of highly fractured rock and loose sediment, and make up the rollover. Type (2) failures are less frequent and larger rockslides and rockfalls within the internal bedded and fractured slope along pre‐determined potential failure surfaces, and make up the power‐law tail. Rockslides and rockfalls of high magnitude and relatively low frequency make up 99% of the total failure volume and are thus responsible for the high erosion rate. They are also significant in the context of landscape evolution as they occur on slopes above 45° and limit the relief of the slope. Copyright © 2012 John Wiley & Sons, Ltd.     

Topographic correlations with soil and regolith thickness from shallow‐seismic refraction constraints across upland hillslopes in the Valles Caldera,New Mexico

How rock is weathered physically and chemically into transportable material is a fundamental question in critical‐zone science. In addition, the distribution of weathered material (soil and intact regolith) across upland landscapes exerts a first‐order control on the hydrology of watersheds. In this paper we present the results of six shallow seismic‐refraction surveys in the Redondo Mountain region of the Valles Caldera, New Mexico. The P‐wave velocities corresponding to soil (≤ 0.6 km s^?1) were inferred from a seventh seismic survey where soil‐thickness data were determined by pit excavation. Using multivariable regression, we quantified the relationships among slope gradient, aspect, and topographic wetness index (TWI) on soil and regolith (soil plus intact regolith) thicknesses. Our results show that both soil and regolith thicknesses vary inversely with TWI in all six survey areas while varying directly with slope aspect (i.e. thicker beneath north‐facing slopes) and inversely with slope gradient (i.e. thinner beneath steep slopes) in the majority of the survey areas. An empirical model based on power‐law relationships between regolith thickness and its correlative variables can fit our inferred thicknesses with R² ‐values up to 0.880 for soil and 0.831 for regolith in areas with significant topographic variations. These results further demonstrate the efficacy of shallow seismic refraction for mapping and determining how soil and regolith variations correlate with topography across upland landscapes. Copyright © 2016 John Wiley & Sons, Ltd.     

Divergent evolution in fluviokarst landscapes of central Kentucky

Central Kentucky is characterized by a mixture of karst and ?uvial features, typically manifested as mosaic of karst‐rich/channel‐poor (KRCP) and channel‐rich/karst‐poor (CRKP) environments. At the regional scale the location and distribution of KRCP and CRKP areas are not always systematically related to structural, lithological, topographic, or other controls. This study examines the relationship of KRCP and CRKP zones along the Kentucky River gorge area, where rapid incision in the last 1·5 million years has lowered local base levels and modi?ed slopes on the edge of the inner bluegrass plateau. At the scale of detailed ?eld mapping on foot within a 4 km² area, the development of karst and ?uvial features is controlled by highly localized structural and topographic constraints, and can be related to slope changes associated with retreat of the Kentucky River gorge escarpment. A conceptual model of karst/?uvial transitions is presented, which suggests that minor, localized variations are suf?cient to trigger a karst–?uvial or ?uvial–karst switch when critical slope thresholds are crossed. Copyright © 2004 John Wiley & Sons, Ltd.     

岩浆岩和变质岩边坡岩体结构类型的划分

在研究岩浆岩和变质岩特点的基础上,根据岩浆岩岩体结构体的大小、有无缓倾的结构面及岩浆岩散粒状边坡岩体的特征,将岩浆岩边坡岩体结构划分为4大类7个亚类;根据部分正变质岩边坡岩体结构与岩浆岩边坡岩体结构类似,部分副变质岩中发育的板理、千枚理、片理类似沉积岩中的层理,以及边坡岩体具有似层状岩体结构的特征,把变质岩边坡岩体的结构划分为5大类11个亚类。岩体结构的划分为边坡稳定性研究、边坡病害治理及其工程效果评价提供了依据。     

Influence of bedrock structures on the spatial pattern of erosional landforms in small alpine catchments

Structural settings and lithological characteristics are traditionally assumed to influence the development of erosional landforms, such as gully networks and rock couloirs, in steep mountain rock basins. The structural control of erosion of two small alpine catchments of distinctive rock types is evaluated by comparing the correspondences between the orientations of their gullies and rock couloirs with (1) the sliding orientations of potential slope failures mechanisms, and (2) the orientation of the maximum joint frequency, this latter being considered as the direction exploited primarily by erosion and mass wasting processes. These characteristic orientations can be interpreted as structural weaknesses contributing to the initiation and propagation of erosion. The morphostructural analysis was performed using digital elevation models and field observations. The catchment comprised of magmatic intrusive rocks shows a clear structural control, mostly expressed through potential wedges failure. Such joint configurations have a particular geometry that encourages the development of gullies in hard rock, e.g. through enhanced gravitational and hydrological erosional processes. In the catchment underlain by sedimentary rocks, penetrative joints that act as structural weaknesses seem to be exploited by gullies and rock couloirs. However, the lithological setting and bedding configuration prominently control the development of erosional landforms, and influence not only the local pattern of geomorphic features, but the general morphology of the catchment. The orientations of the maximum joint frequency are clearly associated with the gully network, suggesting that its development is governed by anisotropy in rock strength. These two catchments are typical of bedrock‐dominated basins prone to intense processes of debris supply. This study suggests a quantitative approach for describing the relationship between bedrock jointing and geomorphic features geometry. Incorporation of bedrock structure can be relevant when studying processes governing the transfer of clastic material, for the assessment of sediment yields and in landforms evolution models. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of regional material yield from coastal landslides based on historical digital terrain modelling

High‐resolution historical (1942) and recent (1994) digital terrain models were derived from aerial photographs along the Big Sur coastline in central California to measure the long‐term volume of material that enters the nearshore environment. During the 52‐year measurement time period, an average of 21 000 ± 3100 m³ km^?1 a^?1 of material was eroded from nine study sections distributed along the coast, with a low yield of 1000 ± 240 m³ km^?1 a^?1 and a high of 46 700 ± 7300 m³ km^?1 a^?1. The results compare well with known volumes from several deep‐seated landslides in the area and suggest that the processes by which material is delivered to the coast are episodic in nature. In addition, a number of parameters are investigated to determine what influences the substantial variation in yield along the coast. It is found that the magnitude of regional coastal landslide sediment yield is primarily related to the physical strength of the slope‐forming material. Coastal Highway 1 runs along the lower portion of the slope along this stretch of coastline, and winter storms frequently damage the highway. The California Department of Transportation is responsible for maintaining this scenic highway while minimizing the impacts to the coastal ecosystems that are part of the Monterey Bay National Marine Sanctuary. This study provides environmental managers with critical background data on the volumes of material that historically enter the nearshore from landslides, as well as demonstrating the application of deriving historical digital terrain data to model landscape evolution. Published in 2005 by John Wiley & Sons, Ltd.     

Parameterizing a physically based shallow landslide model in a data poor region

Shallow landslides and consequent debris flows are an increasing concern in the Western Ghats of Kerala, India. Their increased frequency has been associated with deforestation and unfavourable land‐use practices in cultivated areas. In order to evaluate the influence of vegetation on shallow slope failures a physically based, dynamic and distributed hydrological model (STARWARS) coupled with a probabilistic slope stability model (PROBSTAB) was applied to the upper Tikovil River basin (55·6 km²). It was tuned with the limited evidence of groundwater conditions during the monsoon season of 2005 and validated against observed landslide activity in the hydrological year 2001–2002. Given the data poor conditions in the region some modifications to the original model were in order, including the estimation of parameters on the basis of generalized information from secondary sources, pedo‐transfer functions, empirical equations and satellite remote sensing data. Despite the poor input, the model captured the general temporal and spatial pattern of instability in the area. Sensitivity analysis proved root cohesion, soil depth and angle of internal friction as the most dominant parameters influencing slope stability. The results indicate the importance of root cohesion in maintaining stability and the critical role of the management of rubber plantations in this. Interception and evapotranspiration showed little influence on the development of failure conditions. The study also highlights the importance of high resolution digital terrain models for the accurate mechanistic prediction of shallow landslide initiation. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatial patterns of suspended sediment yields in a humid tropical watershed in Costa Rica

An Erratum has been published for this article in Hydrological Processes 16(5) 2002, 1130–1131. Humid tropical regions are often characterized by extreme variability of fluvial processes. The Rio Terraba drains the largest river basin, covering 4767 km², in Costa Rica. Mean annual rainfall is 3139±419_sd mm and mean annual discharge is 2168±492_sd mm (1971–88). Loss of forest cover, high rainfall erosivity and geomorphologic instability all have led to considerable degradation of soil and water resources at local to basin scales. Parametric and non‐parametric statistical methods were used to estimate sediment yields. In the Terraba basin, sediment yields per unit area increase from the headwaters to the basin mouth, and the trend is generally robust towards choice of methods (parametric and LOESS) used. This is in contrast to a general view that deposition typically exceeds sediment delivery with increase in basin size. The specific sediment yield increases from 112±11·4_sd t km^?2 year^?1 (at 317·9 km² on a major headwater tributary) to 404±141·7_sd t km^?2 year^?1 (at 4766·7 km²) at the basin mouth (1971–92). The analyses of relationships between sediment yields and basin parameters for the Terraba sub‐basins and for a total of 29 basins all over Costa Rica indicate a strong land use effect related to intensive agriculture besides hydro‐climatology. The best explanation for the observed pattern in the Terraba basin is a combined spatial pattern of land use and rainfall erosivity. These were integrated in a soil erosion index that is related to the observed patterns of sediment yield. Estimated sediment delivery ratios increase with basin area. Intensive agriculture in lower‐lying alluvial fans exposed to highly erosive rainfall contributes a large part of the sediment load. The higher elevation regions, although steep in slope, largely remain under forest, pasture, or tree‐crops. High rainfall erosivity (>7400 MJ mm ha^?1 h^?1 year ^?1) is associated with land uses that provide inadequate soil protection. It is also associated with steep, unstable slopes near the basin mouth. Improvements in land use and soil management in the lower‐lying regions exposed to highly erosive rainfall are recommended, and are especially important to basins in which sediment delivery ratio increases downstream with increasing basin area. Copyright © 2001 John Wiley & Sons, Ltd.