Hillslope Angle Research Articles

A Mechanistic Model and Experiments on Bedrock Incision and Channelization by Rockfall

AbstractRockfall and rock avalanches are common in steep terrain on Earth and potentially on other planetary bodies such as the Moon and Mars. Since impacting rocks can damage exposed bedrock as they roll and bounce downhill, rockfall might be an important erosive agent in steep landscapes, even in the absence of water. We developed a new theory for rockfall‐driven bedrock abrasion using the ballistic trajectories of rocks transported under gravity. We calibrated this theory using laboratory experiments of rockfall over an inclined bedrock simulant. Both the experiments and the model demonstrate that bedrock hillslopes can be abraded by dry rockfall, even at gradients below the angle of repose, depending on the bedrock roughness. Feedback between abrasion and topographic steering of rockfall can produce channel‐like forms, such as bedrock chutes, in the absence of water. Particle size has a dominant influence on abrasion rates and runout distances, while the hillslope angle has a comparatively minor influence. Rockfall transport is sensitive to bedrock roughness; terrain with high friction angles can trap rocks creating patches of rock cover that affect subsequent rockfall pathways. Our results suggest that dry rockfall can play an important role in eroding and channelizing steep, rocky terrain on Earth and other planets, such as crater degradation on the Moon and Mars.

Periglacial resurfacing of hillslopes and channels with large boulders in the Virginia Appalachians

AbstractLarge, resistant, quartz‐rich boulders deposited on hillslopes and in channels armour the landscape, trap sediment and influence hillslope angle and erodibility. In the Virginia Appalachians, such boulders are a significant component of hillslopes and channels. Establishing the timing of and processes responsible for bedrock fracture and boulder deposition is a critical piece of understanding the landscape as a system. In this study, we use cosmogenic 10Be exposure dating to resolve the timing of boulder deposition at three sites in the Virginia Valley and Ridge province: Gap Mountain, Brush Mountain and Little Stony Creek, and at one site in the Virginia Blue Ridge: Devil's Marbleyard. The correlation between measured boulder exposure ages (101.7 ± 6.9 ka to 10.8 ± 0.8 ka; n = 23) and the Wisconsin Glacial Stage and subsequent Laurentide Ice Sheet (LIS) deglaciation (~115–11.7 ka) suggests a periglacial origin for deposition of large hillslope and channel boulders in the Virginia Appalachians. The lack of boulder exposure ages corresponding to the Last Interglacial Stage or following Wisconsin LIS retreat suggests interglacial non‐deposition and stability. The absence of exposure ages from the penultimate Illinoian or older Quaternary Glacial Stages suggests that periglacial hillslope processes allow the landscape to be resurfaced with large boulders during each return to cold climate conditions. This cyclic resurfacing of hillslopes and channels is an example of how climatic oscillations insert disequilibrium into the landscape cycle and contributes to our appreciation of the timescales over which contemporary climate change may impact boulder dominated landscapes in rapidly warming alpine and arctic environments.

Transient Response and Adjustment Timescales of Channel Width and Angle of Valley‐Side Slopes to Accelerated Incision

AbstractStudying bedrock rivers during their transient states helps understand the response of a fluvial system to the changing boundary conditions. Although theoretical studies predict how river form adjusts to changes in incision or rock uplift rates, field constraints on the timescale of this adjustment are limited. We investigated the transient behavior of channels and hillslopes and estimated the adjustment times of channel width and angle of valley‐side slopes to accelerated incision based on knickpoint travel time. We documented channel slopes, channel widths, and hillslope angles along six rivers around an active normal fault in Iwaki, Japan, and identified river sections in a transient state. Channel slopes and basin‐averaged erosion rates determined from 10Be concentrations are distinct between rivers near and distant from the fault, suggesting that past increases in fault throw rates triggered the knickpoint formation and the observed transient response. Adjustment time for width is at least 105 years and can take 106 years after the knickpoint passage. Adjustment time for hillslope angles is generally shorter than for channel width. However, the hillslope adjustment may take longer than previously reported if the effect of width adjustment on hillslope angles is significant or the complete adjustment of hillslope angles is considered. The fact that channel slope, channel width, and hillslope angle have distinct adjustment times underlines the importance of correctly identifying river sections that are fully adjusted to the new boundary conditions when inferring erosion or relative uplift rates for bedrock rivers.

Multi-decadal erosion rates from glacierized watersheds on Mount Baker, Washington, USA, reveal topographic, climatic, and lithologic controls on sediment yields

Understanding land surface change in and sediment export out of proglacial landscapes is critical for understanding geohazard and flood risks over engineering timescales and characterizing landscape evolution over geomorphic timescales. We used automated Structure from Motion software to process historical aerial photographs and, with modern lidar data, generated a high-resolution DEM time series with coverage over 10 glacierized watersheds on Mount Baker, Washington, USA for the time period between 1947 and 2015. We measured basin-wide sediment yields and sediment redistribution on hillslopes and in stream channels. Slopes within most measured erosion sites are above theoretical and observed debris-flow thresholds. We observed significant erosion of hillslopes and limited deposition on hillslopes and in stream channels. Sediment delivery ratios during time periods with net erosion averaged 0.73. We determined, consistent with previous field observations, that debris flows originating from moraines are a primary erosion mechanism in proglacial zones on Mount Baker. Time series measurements indicate that temporal variability in erosion rates is associated with climate oscillations, with higher erosion rates during cooler-wetter periods. Basin-wide sediment yield is positively correlated with lithology (r2 = 0.54), hillslope angle (r2 = 0.52), drainage area (r2 = 0.82), and negatively correlated with stream channel slope (r2 = 0.67). Topographic differences between high and low yielding basins indicate that spatial variability in erosion on Mount Baker is sensitive to Pleistocene and Holocene glacial and volcanic activity. Specific sediment yields in six basins averaged 4600 ton/km2/yr, consistent with global measurements in glacierized catchments. Specific sediment yield decreased with increasing basin area, with total loads in the downstream main stem Nooksack River estimated between 480 and 820 ton/km2/yr. Proglacial sediment yields account for between 18 and 32 % of total sediment load in the main stem Nooksack River and exceed contributions by bluff and terrace erosion, which account for between 8 and 13 % of total load. Our findings indicate that erosion in glacierized basins is sensitive to decadal climate oscillations and that high proglacial sediment yields provide an important contribution to river systems downstream, particularly in catchments where upland topography and lithology is favorable.

Experimental investigation on near‐field edge wave run‐ups generated by landslides in narrow reservoirs

Edge waves generated by subaerial landslide‐tsunamis (SLTs) often have potentially huge amplitudes that can endanger human lives, offshore structures, coastal cities and port facilities around reservoirs. Maximum edge wave run‐up (RuM) is the most important parameter for SLT hazard mitigation in an engineering context. This study aimed to investigate the effects of the hill slope angle α on RuM, the first edge wave run‐up Ru1 and the first edge wave decay in the near field Ru1(r/h) through 84 laboratory experiments, which were conducted in a three‐dimensional river channel with α values varying from 15° to 40°. For smaller α values, we found that RuM increased with the increase in α while for large α values exceeding 30°, RuM decreased with the increase in α. An empirical equation considering the effect of α was proposed for the prediction of RuM in the near field, which was well verified by a real case. Additionally, the first edge wave run‐up decay along the shore is also discussed in this paper. This work could have some favourable implications for the preliminary hazard assessment of SLTs.

Response of modern fluvial sediments to regional tectonic activity along the upper Min River, eastern Tibet

Abstract. The deposition of fluvial sediments in tectonically active areas is mainly controlled by tectonics, climate, and associated Earth surface processes; consequently, fluvial sediments can provide a valuable record of changes in regional climate and tectonic activity. In this study, we conducted a detailed analysis of the grain-size distribution in modern fluvial sediments from the upper Min River, eastern Tibet. These data, combined with information on regional climate, vegetation, hydrology, geomorphology, lithology, and fault slip rate, indicate that modern regional tectonic activity along upper Min River can be divided into three segments. Specifically, fluvial sediments in the Minjiangyuan–Diexi segment are dominated by silts (< 63 µm, 70.2 %), agreeing with low runoff, low rainfall, and high vegetation cover and revealing a windblown origin influenced by the arid and windy climate. These observations are consistent with the low hillslope angle and low relief, all indicating weak activity along the Minjiang Fault. The coarse-grained fraction (> 250 µm) of fluvial sediments in the Diexi–Wenchuan and Wenchuan–Dujiangyan segments increases stepwise downstream, although runoff and rainfall do not change significantly. These patterns correlate well with increases in both regional relief and hillslope angles. Together, these observations imply that regional tectonic activity along the Maoxian–Wenchuan Fault becomes more pervasive downstream along the Min River. The occurrence of well-sorted and well-rounded pebbles of fluvial sediments downstream of Dujiangyan must be related to the long-time scouring and sorting by rivers. This study marks the first development of a new research approach that can characterize regional tectonic activity by analysis of grain-size distribution of fluvial sediments collected from tectonically active regions.

Lithological and structural controls on river profiles and networks in the northern Sierra Nevada (California, USA)

Abstract In this study, the strong lithological heterogeneity of the northern Sierra Nevada (California, USA) is exploited to elucidate the role of lithology on river profiles and patterns at the mountain-range scale. The analyses indicate that plutonic, metavolcanic, and quartzite bedrock generally host the steepest river reaches, whereas gentle reaches flow across non-quartzite metasedimentary rocks and fault zones. In addition, the largest immobile boulders are often in the steepest reaches, suggesting that wide joint spacing plays a role in creating steep channels, and a positive relationship between boulder size and hillslope angle highlights the coupling of the hillslope and fluvial systems. With respect to river network configurations, dendritic patterns dominate in the plutonic bedrock, with channels aligned down the slope of the range; in contrast, river reaches in the metamorphic belts are mainly longitudinal and parallel to the structural grain. River profiles and patterns in the northern Sierra Nevada, therefore, bear a strong lithological imprint related to differential erosion. These observations indicate that attempts to infer uplift and tilting of the range based on the gradients and orientations of paleochannel remnants should first account for the effect of bedrock erodibility. Indeed, the differences in gradients of Tertiary paleochannel remnants used to argue for late Cenozoic uplift of the range can be wholly explained by differences in lithology.

Cosmogenic 10Be in river sediment: where grain size matters and why

Abstract. Concentrations of in-situ-produced cosmogenic 10Be in river sediment are widely used to estimate catchment-average denudation rates. Typically, the 10Be concentrations are measured in the sand fraction of river sediment. However, the grain size of bedload sediment in most bedrock rivers covers a much wider range. Where 10Be concentrations depend on grain size, denudation rate estimates based on the sand fraction alone are potentially biased. To date, knowledge about catchment attributes that may induce grain-size-dependent 10Be concentrations is incomplete or has only been investigated in modelling studies. Here we present an empirical study on the occurrence of grain-size-dependent 10Be concentrations and the potential controls of hillslope angle, precipitation, lithology, and abrasion. We first conducted a study focusing on the sole effect of precipitation in four granitic catchments located on a climate gradient in the Chilean Coastal Cordillera. We found that observed grain size dependencies of 10Be concentrations in the most-arid and most-humid catchments could be explained by the effect of precipitation on both the scouring depth of erosion processes and the depth of the mixed soil layer. Analysis of a global dataset of published 10Be concentrations in different grain sizes (n=73 catchments) – comprising catchments with contrasting hillslope angles, climate, lithology, and catchment size – revealed a similar pattern. Lower 10Be concentrations in coarse grains (defined as “negative grain size dependency”) emerge frequently in catchments which likely have thin soil and where deep-seated erosion processes (e.g. landslides) excavate grains over a larger depth interval. These catchments include steep (> 25∘) and humid catchments (> 2000 mm yr−1). Furthermore, we found that an additional cause of negative grain size dependencies may emerge in large catchments with weak lithologies and long sediment travel distances (> 2300–7000 m, depending on lithology) where abrasion may lead to a grain size distribution that is not representative for the entire catchment. The results of this study can be used to evaluate whether catchment-average denudation rates are likely to be biased in particular catchments.

Short communication: Increasing vertical attenuation length of cosmogenic nuclide production on steep slopes negates topographic shielding corrections for catchment erosion rates

Abstract. Interpreting catchment-mean erosion rates from in situ produced cosmogenic 10Be concentrations in stream sediments requires calculating the catchment-mean 10Be surface production rate and effective mass attenuation length, both of which can vary locally due to topographic shielding and slope effects. The most common method for calculating topographic shielding accounts only for the reduction of nuclide production rates due to shielding at the surface, leading to catchment-mean corrections of up to 20 % in steep landscapes, and makes the simplifying assumption that the effective mass attenuation length for a given nuclide production mechanism is spatially uniform. Here I evaluate the validity of this assumption using a simplified catchment geometry with mean slopes ranging from 0 to 80∘ to calculate the spatial variation in surface skyline shielding, effective mass attenuation length, and the total effective shielding factor, defined as the ratio of the shielded surface nuclide concentration to that of an unshielded horizontal surface. For flat catchments (i.e., uniform elevation of bounding ridgelines), the effect of increasing vertical attenuation length as a function of hillslope angle and skyline shielding exactly offsets the effect of decreasing surface production rate, indicating that no topographic shielding correction is needed when calculating catchment-mean vertical erosion rates. For dipping catchments (as characterized by a plane fit to the bounding ridgelines), the catchment-mean surface nuclide concentrations are also equal to that of an unshielded horizontal surface, except for cases of extremely steep range-front catchments, where the surface nuclide concentrations are counterintuitively higher than the unshielded case due to added production from oblique cosmic ray paths at depth. These results indicate that in most cases topographic shielding corrections are inappropriate for calculating catchment-mean erosion rates, and are only needed for steep catchments with nonuniform distributions of quartz and/or erosion rate. By only accounting for shielding of surface production, existing shielding approaches introduce a slope-dependent systematic error that could lead to spurious interpretations of relationships between topography and erosion rate.

Impulse waves generated by subaerial landslides of combined block mass and granular material

Experimental results, based on two-dimensional (2-D) laboratory tests, for impulse waves generated by subaerial landslides of combined solid block and granular materials are presented in this study. The results are compared with those of individual models of pure solid block and granular landslides. By considering an identical slide mass and release height, the effect of the mass ratio m* on wave generation was investigated. The mass ratio is defined as m* = MB/M, where MB is the mass of the solid block, M = MB + MG is the total mass and MG is the mass of the granular portion. The experimental results show that the combined landslides with m* = 0.6 and 0.8 in this study generally produce much larger impulse waves in the impact zone compared with those triggered by pure solid block landslides (m* = 1) and pure granular landslides (m* = 0). This suggests that the primary wave amplitudes of impulse waves might have been underestimated in previous laboratory tests with solely solid or granular assemblies when using the same slide mass and release height. The larger primary wave amplitudes of the combined landslides compared with those of the pure solid block landslides are mainly attributed to the relatively large thickness of the combined landslides and the continuous motion of the granular portion, as the inside solid block stops at the hill slope toe. Compared with pure granular landslides, combined landslides generally have a larger Froude number and slide thickness, which account for the larger primary wave amplitudes. The latter plays a quantifiably more important role – nearly two times that of the former. The effects of the hill slope angle α, and the grain size of the granular materials D, were also studied in this work. Four different hill slope angles (α = 22.5°, 30°, 37.5° and 45°) and four different mean grain diameters (D = 5 mm, 10 mm, 20 mm and 30 mm) were tested. By comparing the measured maxima of the wavelet spectra obtained via a wavelet transform method with those reconstructed by means of the phase celerity of the solitary wave cs and the group celerity of linear wave cg(f), it is found that the measured maxima travel at either cs or cg(f), depending on the landslide type and the hill slope angle α. However, the celerity of the measured maxima for the high frequency is lower than cg(f) and cs for the impulse waves generated by combined landslides when α ≤ 37.5° in this study. Further investigations show that the wave amplitude decreases rapidly during wave propagation, following an exponential function. The decrease in the wave amplitude during propagation is mainly attributed to dispersion.

Topographic Loads Modified by Fluvial Incision Impact Fault Activity in the Longmenshan Thrust Belt, Eastern Margin of the Tibetan Plateau

AbstractWhether external or internal forces of the Earth control the behaviors of upper‐crustal faults in a fold‐and‐thrust belt has been debated for decades. The Longmenshan thrust belt (LTB) along the eastern margin of the Tibetan Plateau may provide insights into such a debate, as the central segment of the LTB has relatively uniform shortening strains yet various fluvial incision capability along the strike. This tectonic setting enables a better assessment of the effects of external forces on fault activity. We analyzed the variations of topography, fluvial incision intensity, coseismic slips, and coseismic landslides along the central LTB. The Longmen subsegment in the northern half has higher elevation and lower fluvial incision intensity than the Hongkou subsegment in the south. We calculated the topographic stresses on the faults ruptured during the 2008 Wenchuan earthquake and found topographic introduced normal stress increase may explain the coseismic slip partitioning onto two subfaults along the Longmen subsegment. Our results indicate that fluvial incision may have produced the along‐strike variations of the topography, which may further produce the different rupture behavior. In addition, the mean hillslope angle along the central LTB appeared to be at the critical condition prior to the 2008 Wenchuan earthquake and reduced significantly by the earthquake, indicating that geomorphic indices may vary with different stages in an earthquake cycle. Therefore, scrutinizing the mean hillslope angle and other geomorphic indices may help identify potential seismic hazards in an active fault system.

Tectonic controls of Holocene erosion in a glaciated orogen

Abstract. Recent work has highlighted a strong, worldwide, alpine glacial impact on orogen erosion rates over the last 2 Ma. While it may be assumed that glaciers increased erosion rates when active, the degree to which past glaciations influence Holocene erosion rates through the adjustment of topography is not known. In this study, we investigate the influence of long-term tectonic and post-glacial topographic controls on erosion in a glaciated orogen: the Olympic Mountains, USA. We present 14 new 10Be and 26Al analyses which constrain Holocene erosion rates across the Olympic Mountains. Basin-averaged erosion rates scale with basin-averaged values of 5 km local relief, channel steepness, and hillslope angle throughout the range, similar to observations from non-glaciated orogens. These erosion rates are not related to mean annual precipitation or the marked change in Pleistocene alpine glacier size across the range, implying that glacier modification of topography and modern precipitation parameters do not exert strong controls on these rates. Rather, we find that despite spatial variations in glacial modification of topography, patterns of recent erosion are similar to those from estimates of long-term tectonic rock uplift. This is consistent with a tectonic model where erosion and rock uplift patterns are controlled by the deformation of the Cascadia subduction zone.

The influence of vegetation on debris-flow initiation during extreme rainfall in the northern Colorado Front Range

Abstract We explored regional influences on debris-flow initiation throughout the Colorado Front Range (Colorado, USA) by exploiting a unique data set of more than 1100 debris flows that initiated during a 5 day rainstorm in 2013. Using geospatial data, we examined the influence of rain, hillslope angle, hillslope aspect, and vegetation density on debris-flow initiation. In particular we used a greenness index to differentiate areas of high tree density from grass and bare soil. The data demonstrated an overwhelming propensity for debris-flow initiation on south-facing hillslopes. However, when the debris-flow density was analyzed with respect to total rainfall and greenness we found that most debris flows occurred in areas of high rainfall and low tree density, regardless of hillslope aspect. These results indicate that present-day tree density exerts a stronger influence on debris-flow initiation locations than aspect-driven variations in soil and bedrock properties that developed over longer time scales.

Assessing uncertainty and sensor biases in passive microwave data across High Mountain Asia

Abstract Snowfall comprises a significant percentage of the annual water budget in High Mountain Asia (HMA), but snow-water equivalent (SWE) is poorly constrained due to lack of in-situ measurements and complex terrain that limits the efficacy of modeling and observations. Over the past few decades, SWE has been estimated with passive microwave (PM) sensors with generally good results in wide, flat, terrain, and lower reliability in densely forested, complex, or high-elevation areas. In this study, we use raw swath data from five satellite sensors — the Special Sensor Microwave/Imager (SSMI) and Special Sensor Microwave Imager/Sounder (SSMIS) (1987–2015, F08, F11, F13, F17), Advanced Microwave Scanning Radiometer — Earth Observing System (AMSR-E, 2002–2011), AMSR2 (2012–2015), and the Global Precipitation Measurement (GPM, 2014–2015) — in order to understand the spatial and temporal structure of native sensor, topographic, and land cover biases in SWE estimates in HMA. We develop a thorough understanding of the uncertainties in our SWE estimates by examining the impacts of topographic parameters (aspect, relief, hillslope angle, and elevation), land cover, native sensor biases, and climate parameters (precipitation, temperature, and wind speed). HMA, with its high seasonality, large topographic gradients and low relief at high elevations provides an excellent context to examine a wide range of climatic, land-cover, and topographic settings to better constrain SWE uncertainties and potential sensor bias. Using a multi-parameter regression, we compare long-term SWE variability to forest fraction, maximal multi-year snow depth, topographic parameters, and long-term average wind speed across both individual sensor time series and a merged multi-sensor dataset. In regions where forest cover is extensive, it is the strongest control on SWE variability. In those regions where forest density is low ( While forest cover and elevation have been integrated into many SWE algorithms, wind speed and long-term maximal snow depth have not. Our results show that wind redistribution of snow can have impacts on SWE, especially over large, flat, areas. Using our regression results, we have developed an understanding of sensor-specific SWE uncertainties and their spatial patterns. The uncertainty maps developed in this study provide a first-order approximation of SWE-estimate reliability for much of HMA, and imply that high-fidelity SWE estimates can be produced for many high-elevation areas.

Particle transport over rough hillslope surfaces by dry ravel: Experiments and simulations with implications for nonlocal sediment flux

Past studies of hillslope evolution have typically assumed that soil creep processes are governed by a linear relationship between local hillslope angle and transport distance. The assumption of “linear diffusion” has fallen out of favor because, when coupled with an expression of mass continuity, it yields unrealistic hillslope profiles. As a consequence, a better understanding of the mechanics of sediment transport is needed. Here we report results from a series of flume experiments performed to investigate sediment transport by dry ravel, a common soil creep process in arid and semiarid environments. We find that, at gentle slopes, transport distances follow distributions characteristic of local transport. As gradients steepen, a fraction of the particles begins to exhibit nonlocal transport, and that fraction increases rapidly with slope. A stochastic discrete element model that couples an effective friction term with a shock term reproduces the results from the flume experiments, suggesting that it can be used to explore the nature of particle transport on rough surfaces. The model predicts that exponential distributions of transport distances on gentle slopes evolve into quasi‐uniform distributions on steep slopes, and the transition occurs as slopes approach the angle of repose. Our results support previous findings that the angle of repose represents a threshold between friction and inertial regimes. In addition, we propose that the angle of repose represents a fuzzy boundary between local and nonlocal transport.

Hillslope response to tectonic forcing in threshold landscapes

ABSTRACTHillslopes 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 (10Be)‐based dataset of catchment‐averaged erosion rates combined with a 500 km2 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.

Landslide causes: Human impacts on a Himalayan landslide swarm

Geographers are commonly called upon to diagnose the causes of natural disasters and provide guidance for policy makers. This case study concerns one of the landslide and flood disasters that afflicted the Himalaya in the later Monsoon of 2010, the Almora District (Uttarakhand), landslide swarm, which affected thousands of kilometres of roadway and caused damages >$ 125,000,000. Causes were sought on a 7.4km reach of roadway that had been part of a previous study. Individual landslides were surveyed along with a suite of potential environmental indicators. These results show that the numbers of landslides on the roadbed were not unusual, so research switched to exploring factors that may have enhanced their volumes. Local opinion considers that the major causes were human impacts, especially recent urban development, or geological weaknesses. However, results show that human impacts such as building construction, deforestation, and the collapse of agricultural terraces were associated with smaller landslides, while geological factors showed no significant correlations with landslide volume. Road engineering factors such as the location of another road upslope and the collapse of roadcut retaining wall also proved unimportant. However, the height of the roadcut, channel incision undermining the roadcut and, especially, the steepness of the hillside upslope proved significant. Despite the absence of any direct link between human activity and landslide volumes in Almora, a comparison of the 2010 results from this survey and that of another reach of hill road, which was also part of the 1985 survey, showed important differences in landslide numbers.In 1985, the Kilbury Road, which cuts through reserved forest, carried 153 landslides to the Almora Lower Mall’s 88. In 2010, there were only 9 landslides on the Kilbury Road but 108 on the suburban Almora Lower Mall. The conclusion is that human impacts in Almora have combined to prevent the landslides healing as they have in the undisturbed forest along the Kilbury Road. However, past surveys confirm that landslide numbers on the Almora Lower Mall have remained steady since 1985. Since, the impacts of this rainfall event were manifest only as increased landslide outfall, not increased numbers, the chief influences of outfall volume remain hillslope angle, height of the initial roadcut, undermining by channel incision and the fact that these hillsides exist close to criticality, their limits of stability.

A 1-D mechanistic model for the evolution of earthflow-prone hillslopes

[1] In mountainous terrain, deep-seated landslides transport large volumes of material on hillslopes, exerting a dominant control on erosion rates and landscape form. Here, we develop a mathematical landscape evolution model to explore interactions between deep-seated earthflows, soil creep, and gully processes at the drainage basin scale over geomorphically relevant (>103 year) timescales. In the model, sediment flux or incision laws for these three geomorphic processes combine to determine the morphology of actively uplifting and eroding steady state topographic profiles. We apply the model to three sites, one in the Gabilan Mesa, California, with no earthflow activity, and two along the Eel River, California, with different lithologies and varying levels of historic earthflow activity. Representative topographic profiles from these sites are consistent with model predictions in which the magnitude of a dimensionless earthflow number, based on a non-Newtonian flow rheology, reflects the magnitude of recent earthflow activity on the different hillslopes. The model accurately predicts the behavior of earthflow collection and transport zones observed in the field and estimates long-term average sediment fluxes that are due to earthflows, in agreement with historical rates at our field sites. Finally, our model predicts that steady state hillslope relief in earthflow-prone terrain increases nonlinearly with the tectonic uplift rate, suggesting that the mean hillslope angle may record uplift rate in earthflow-prone landscapes even at high uplift rates, where threshold slope processes normally limit further topographic development.

Characteristics of concentrated flow hydraulics for rangeland ecosystems: implications for hydrologic modeling

ABSTRACTConcentrated flow is often the dominant source of water erosion following disturbance on rangelands. Because of the lack of studies that explain the hydraulics of concentrated flow on rangelands, cropland‐based equations have typically been used for rangeland hydrology and erosion modeling, leading to less accurate predictions due to different soil and vegetation cover characteristics. This study investigates the hydraulics of concentrated flow using unconfined field experimental data over diverse rangeland landscapes within the Great Basin Region, United States. The results imply that the overall hydraulics of concentrated flow on rangelands differ significantly from those of cropland rills. Concentrated flow hydraulics on rangelands are largely controlled by the amount of cover or bare soil and hillslope angle. New predictive equations for concentrated flow velocity (R2 = 0·47), hydraulic friction (R2 = 0·52), and width (R2 = 0·4) representing a diverse set of rangeland environments were developed. The resulting equations are applicable across a wide span of ecological sites, soils, slopes, and vegetation and ground cover conditions and can be used by physically‐based rangeland hydrology and erosion models to estimate rangeland concentrated flow hydraulic parameters. Published in 2011. This article is a US Government work and is in the public domain in the USA.

Bioturbation on a south-east Australian hillslope: estimating contributions to soil flux

The contribution of bioturbation to downslope soil transport is significant in many situations, particularly in the context of soil formation, erosion and creep. This study explored the direct flux of soil caused by Aphaenogaster ant mounding, vertebrate scraping and tree-throw on a wildfire-affected hillslope in south-east Australia. This included the development of methods previously applied to Californian gopher bioturbation, and an evaluation of methods for estimating the volume of soil displaced by tree-throw events. All three bioturbation types resulted in a net downslope flux, but any influence of hillslope angle on flux rates appeared to be overshadowed by environmental controls over the spatial extent of bioturbation. As a result, the highest flux rates occurred on the footslope and lower slope. The overall contribution of vertebrate scraping (57.0 ± 89.4 g m−1 yr−1) exceeded that of ant mounding (36.4 ± 66.0 g m−1 yr−1), although mean rates were subject to considerable uncertainty. Tree-throw events, which individually cause major disturbance, were limited in their importance by their scarcity relative to faunalturbation. However, tree-throw might be the dominant mechanism of biotic soil flux on the mid-slope provided that it occurs at a frequency of at least 2–3 events ha−1 yr−1. Although direct biotic soil flux appears to be geomorphologically significant on this hillslope, such transport processes are probably subordinate to other impacts of bioturbation at this site such as the enhancement of infiltration following wildfire. Copyright © 2011 John Wiley & Sons, Ltd.