Effects of gradient,distance, curvature and aspect on steep burned and unburned hillslope soil erosion and deposition

Wildfire denudes vegetation and impacts chemical and physical soil properties, which can alter hillslope erosion rates. Post‐wildfire erosion can also contribute disproportionately to long‐term erosion rates and landscape evolution. Post‐fire hillslope erosion rates remain difficult to predict and document at the hillslope scale. Here we use ²¹⁰Pb_aex (lead‐210 mineral‐adsorbed excess) inventories to describe net sediment erosion on steep, convex hillslopes in three basins (unburned, moderately and severely burned) in mountainous central Idaho. We analyzed nearly 300 soil samples for ²¹⁰Pb_aex content with alpha spectrometry and related net sediment erosion to burn severity, aspect, gradient, curvature and distance from ridgetop. We also tested our data against models for advective, linear and non‐linear diffusive erosion. Statistically lower net soil losses on north‐ versus south‐facing unburned hillslopes suggest that greater vegetative cover and soil cohesion on north‐facing slopes decrease erosion. On burned hillslopes, erosion differences between aspects were less apparent and net erosion was more variable, indicating that vegetation influences erosion magnitude and fire drives erosion variability. We estimated net soil losses throughout the length of unburned hillslopes, including through a footslope transition to concave form. In contrast, on burned hillslopes, the subtle shift from convex to concave form was associated with deposition of a post‐fire erosion pulse. Such overall patterns of erosion and deposition are consistent with predictions from a non‐linear diffusion equation. This finding also suggests that concave sections of overall convex hillslopes affect post‐disturbance soil erosion and deposition. Despite these patterns, no strong relationships were evident between local net soil losses and gradient, curvature, distance from ridgetop, or erosion predicted with advection or diffusion equations. The observed relationship between gradient and erosion is therefore likely more complex or stochastic than often described theoretically, especially over relatively short timescales (60–100 years). Copyright © 2016 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.     

Hillslope response to knickpoint migration in the Southern Appalachians: implications for the evolution of post‐orogenic landscapes

The southern Appalachians represent a landscape characterized by locally high topographic relief, steep slopes, and frequent mass movement in the absence of significant tectonic forcing for at least the last 200 Ma. The fundamental processes responsible for landscape evolution in a post‐orogenic landscape remain enigmatic. The non‐glaciated Cullasaja River basin of south‐western North Carolina, with uniform lithology, frequent debris flows, and the availability of high‐resolution airborne lidar DEMs, is an ideal natural setting to study landscape evolution in a post‐orogenic landscape through the lens of hillslope–channel coupling. This investigation is limited to channels with upslope contributing areas >2.7 km², a conservative estimate of the transition from fluvial to debris‐flow dominated channel processes. Values of normalized hypsometry, hypsometric integral, and mean slope vs elevation are used for 14 tributary basins and the Cullasaja basin as a whole to characterize landscape evolution following upstream knickpoint migration. Results highlight the existence of a transient spatial relationship between knickpoints present along the fluvial network of the Cullasaja basin and adjacent hillslopes. Metrics of topography (relief, slope gradient) and hillslope activity (landslide frequency) exhibit significant downstream increases below the current position of major knickpoints. The transient effect of knickpoint‐driven channel incision on basin hillslopes is captured by measuring the relief, mean slope steepness, and mass movement frequency of tributary basins and comparing these results with the distance from major knickpoints along the Cullasaja River. A conceptual model of area–elevation and slope distributions is presented that may be representative of post‐orogenic landscape evolution in analogous geologic settings. Importantly, the model explains how knickpoint migration and channel–hillslope coupling is an important factor in tectonically‐inactive (i.e. post‐orogenic) orogens for the maintenance of significant relief, steep slopes, and weathering‐limited hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Asymmetric hillslope erosion following wildfire in Fourmile Canyon,Colorado

Infrequent, high‐magnitude events cause a disproportionate amount of sediment transport on steep hillslopes, but few quantitative data are available that capture these processes. Here we study the influence of wildfire and hillslope aspect on soil erosion in Fourmile Canyon, Colorado. This region experienced the Fourmile Fire of 2010, strong summer convective storms in 2011 and 2012, and extreme flooding in September 2013. We sampled soils shortly after these events and use fallout radionuclides to trace erosion on polar‐ and equatorial‐facing burned slopes and on a polar‐facing unburned slope. Because these radionuclides are concentrated in the upper decimeter of soil, soil inventories are sensitive to erosion by surface runoff. The polar‐facing burned slope had significantly lower cesium‐137 (¹³⁷Cs) and lead‐210 (²¹⁰Pb) inventories (p < 0.05) than either the polar‐facing unburned slope or equatorial‐facing burned slope. Local slope magnitude does not appear to control the erosional response to wildfire, as relatively gently sloping (~20%) polar‐facing positions were severely eroded in the most intensively burned area. Field evidence and soil profile analyses indicate up to 4 cm of local soil erosion on the polar‐facing burned slope, but radionuclide mass balance indicates that much of this was trapped nearby. Using a ¹³⁷Cs‐based erosion model, we find that the burned polar‐facing slope had a net mean sediment loss of 2 mm (~1 kg m^?2) over a one to three year period, which is one to two orders of magnitude higher than longer‐term erosion rates reported for this region. In this part of the Colorado Front Range, strong hillslope asymmetry controls soil moisture and vegetation; polar‐facing slopes support significantly denser pine and fir stands, which fuels more intense wildfires. We conclude that polar‐facing slopes experience the most severe surface erosion following wildfires in this region, indicating that landscape‐scale aridity can control the geomorphic response of hillslopes to wildfires. Copyright © 2018 John Wiley & Sons, Ltd.     

The effects of slope shape and polyacrylamide application on runoff,erosion and nutrient loss from hillslopes under simulated rainfall

Natural hillslopes are mostly composed of complex slope shapes, which significantly affect soil erosion. However, existing studies have mainly focused on uniform slopes to simplify complex hillslopes, and the mechanisms responsible for the influence of slope shape on soil and nutrient losses are still not well understood, especially in the application of soil improvers to reduce soil loss. To investigate the effects of slope shape and polyacrylamide (PAM) application on runoff, soil erosion and nutrient loss, this study conducted artificial field rainfall experiments involving two PAM application rates and nine slope shapes. The results indicate that the average amount of soil loss from convex slopes was 1.5 and 1.3 times greater than that from concave and uniform slopes, respectively, and the average amount of ammonia nitrogen loss and phosphate loss increased by 24.0%–58.6%. Soil and nutrient losses increased as the convexity of the convex slopes increased. For runoff, there was little difference between concave and convex slopes, but the runoff amount for both slopes was greater than that for uniform slopes. After PAM application, the soil loss decreased by more than 90%, and the nutrient loss decreased by 28.2%–68.1%. The application of PAM was most effective in reducing soil erosion and nutrient loss from convex slopes, and it is recommended to appropriately increase the PAM application rate for convex slopes. A strong linear relationship between ammonia nitrogen and phosphate concentrations and sediment concentrations was found in the runoff on slopes with no PAM application. However, this linear relationship weakened for slopes with PAM application. The findings of this study may be valuable for optimizing nonpoint source pollution management in basins.     

Erosion rates in the southern oregon coast range: Evidence for an equilibrium between hillslope erosion and sediment yield

The relationship of hillslope erosion rates and sediment yield is often poorly defined because of short periods of measurement and inherent spatial and temporal variability in erosion processes. In landscapes containing hillslopes crenulated by alternating topographic noses and hollows, estimates of local hillslope erosion rates averaged over long time periods can be obtained by analysing colluvial deposits in the hollows. Hollows act as local traps for a portion of the colluvium transported down hillslopes, and erosion rates can be calculated using the age and size of the deposits and the size of the contributing source area. Analysis of colluvial deposits in nine Oregon Coast Range hollows has yielded average colluvial transport rates into the hollows of about 35cm³cm^?1yr^?1 and average bedrock lowering rates of about 0.07 mm yr^?1 for the last 4000 to 15000 yr. These rates are consistent with maximum bedrock exfoliation rates of about 0.09 mm yr^?1 calculated from six of the hollows, supporting the interpretation that exfoliation rates limit erosion rates on these slopes. Sediment yield measurements from nine Coast Range streams provide similar basin-wide denudation rates of between 0.05 and 0.08mm yr^?1, suggesting an approximate steady-state between sediment production on hillslopes and sediment yield. In addition, modern sediment yields are similar in basins varying in size from 1 to 1500 km², suggesting that erosion rates are spatially uniform and providing additional evidence for an approximate equilibrium in the landscape.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

How long is a hillslope?

Hillslope length is a fundamental attribute of landscapes, intrinsically linked to drainage density, landslide hazard, biogeochemical cycling and hillslope sediment transport. Existing methods to estimate catchment average hillslope lengths include inversion of drainage density or identification of a break in slope–area scaling, where the hillslope domain transitions into the fluvial domain. Here we implement a technique which models flow from point sources on hilltops across pixels in a digital elevation model (DEM), based on flow directions calculated using pixel aspect, until reaching the channel network, defined using recently developed channel extraction algorithms. Through comparisons between these measurement techniques, we show that estimating hillslope length from plots of topographic slope versus drainage area, or by inverting measures of drainage density, systematically underestimates hillslope length. In addition, hillslope lengths estimated by slope–area scaling breaks show large variations between catchments of similar morphology and area. We then use hillslope length–relief structure of landscapes to explore nature of sediment flux operating on a landscape. Distinct topographic forms are predicted for end‐member sediment flux laws which constrain sediment transport on hillslopes as being linearly or nonlinearly dependent on hillslope gradient. Because our method extracts hillslope profiles originating from every ridgetop pixel in a DEM, we show that the resulting population of hillslope length–relief measurements can be used to differentiate between linear and nonlinear sediment transport laws in soil mantled landscapes. We find that across a broad range of sites across the continental United States, topography is consistent with a sediment flux law in which transport is nonlinearly proportional to topographic gradient. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Physically based modelling of sheet erosion (detachment and deposition processes) in complex hillslopes

A one‐dimensional uncoupled model governed by this research is a physics‐based modelling of the rainfall‐runoff induced erosion process. The presented model is composed of three parts of a three‐dimensional (3D) hillslope geometry, a nonlinear storage (kinematic wave) model for hillslope hydrological response, and an unsteady physically based surface erosion model. The 3D hillslope geometry model allows describing of the hillslope morphology by defining their plan shape and profile curvature. By changing these two topographic parameters, nine basic hillslope types are derived. The modelling of hillslope hydrological response is based on a flow continuity equation as the relation of discharge and flow depth is passed on kinematic wave approximation. The erosion model is based on a mass conservation equation for unsteady flow. The model assumes that suspended sediment does not affect flow dynamics. The model also accounts for the effect of flow depth plus loose soil depth on soil detachment. The presented model was run for two different precipitations, slope content, and length, and results were plotted for sediment detachment/deposition rate. Based on the obtained results, in hillslopes with convex and straight profile curvatures, sediment detachment only occurred in the whole length of the hillslope. However, in concave ones, sediment detachment and deposition only occurred together in hillslope. The hillslopes with straight profiles and convergent plans have the highest rate of detachment. Also, results show that most detachment rates occur in convex profile curvatures, which are about 15 times more than in straight profiles. Copyright © 2015 John Wiley & Sons, Ltd.     

Sediment delivery from hillslopes and the Universal Soil Loss Equation: some perceptions and misconceptions

The Universal Soil Loss Equation (USLE) or the revised USLE (RUSLE) are often used together with sediment delivery ratios in order to predict sediment delivery from hillslopes. In using sediment delivery ratios for this purpose, it is assumed that the sediment delivery ratio for a given hillslope does not vary with the amount of erosion occurring in the upslope area. This assumption is false. There is a perception that hillslope erosion is calculated on the basis that hillslopes are, in effect, simply divided into 22·1 m long segments. This perception fails to recognize the fact the inclusion of the 22·1 m length in the calculation has no physical significance but simply produces a value of 1·0 for the slope length factor when slopes have a length equal to that of the unit plot. There is a perception that the slope length factor is inappropriate because not all the dislodged sediment is discharged. This perception fails to recognize that the USLE and the RUSLE actually predict sediment yield from planar surfaces, not the total amount of soil material dislocated and removed some distance by erosion within an area. The application of the USLE/RUSLE to hillslopes also needs to take into account the fact that runoff may not be generated uniformly over that hillslope. This can be achieved by an equation for the slope length factor that takes account of spatial variations in upslope runoff on soil loss from a segment or grid cell. Several alternatives to the USLE event erosivity index have been proposed in order to predict event erosion better than can be achieved using the EI₃₀ index. Most ignore the consequences of changing the event erosivity index on the values for the soil, crop and soil conservation protection factors because there is a misconception that these factors are independent of one another. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

The effect of landform variation on vegetation patterning and related sediment dynamics

Semi‐arid ecosystems are often spatially self‐organized in typical patterns of vegetation bands with high plant cover interspersed with bare soil areas, also known as ‘tiger bush’. In modelling studies, most often, straight planar slopes were used to analyse vegetation patterning. The effect of slope steepness has been investigated widely, and some studies investigated the effects of microtopography and hillslope orientation. However, at the larger catchment scale, the overall form of the landscape may affect vegetation patterning and these more complex landscapes are much more prevalent than straight slopes. Hence, our objective was to determine the effect of landform variation on vegetation patterning and sediment dynamics. We linked two well‐established models that simulate (a) plant growth, death and dispersal of vegetation, and (b) erosion and sedimentation dynamics. The model was tested on a straight planar hillslope and then applied to (i) a set of simple synthetic topographies with varying curvature and (ii) three more complex, real‐world landscapes of distinct morphology. Results show banded vegetation patterning on all synthetic topographies, always perpendicular to the slope gradient. Interestingly, we also found that movement of bands – a debated phenomenon – seems to be dependent on curvature. Vegetation banding was simulated on the slopes of the alluvial fan and along the valley slopes of the dissected and rolling landscapes. In all landscapes, local valleys developed a full vegetation cover induced by water concentration, which is consistent with observations worldwide. Finally, banded vegetation patterns were found to reduce erosion significantly as compared to other vegetation configurations. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

Lithological controls on hillslope sediment supply: insights from landslide activity and grain size distributions

The volumes, rates and grain size distributions of sediment supplied from hillslopes represent the initial input of sediment delivered from upland areas and propagated through sediment routing systems. Moreover, hillslope sediment supply has a significant impact on landscape response time to tectonic and climatic perturbations. However, there are very few detailed field studies characterizing hillslope sediment supply as a function of lithology and delivery process. Here, we present new empirical data from tectonically‐active areas in southern Italy that quantifies how lithology and rock strength control the landslide fluxes and grain size distributions supplied from hillslopes. Landslides are the major source of hillslope sediment supply in this area, and our inventory of ~2800 landslides reveals that landslide sediment flux is dominated by small, shallow landslides. We find that lithology and rock strength modulate the abundance of steep slopes and landslides, and the distribution of landslide sizes. Outcrop‐scale rock strength also controls the grain sizes supplied by bedrock weathering, and influences the degree of coarsening of landslide supply with respect to weathering supply. Finally, we show that hillslope sediment supply largely determines the grain sizes of fluvial export, from catchments and that catchments with greater long‐term landslide rates deliver coarser material. Therefore, our results demonstrate a dual control of lithology on hillslope sediment supply, by modulating both the sediment fluxes from landslides and the grain sizes supplied by hillslopes to the fluvial system. Copyright © 2017 John Wiley & Sons, Ltd.     

Post‐fire thin debris flows: sediment transport and numerical modelling

The creation of a hydrophobic layer in the soil during ?res in semi‐arid environments inhibits the in?ltration of rainfall. This leads to increased rates of runoff and associated sediment transport. When the hydrophobic layer is deposited beneath the soil surface, a perched water table develops which may cause thin (1–2 cm) hillslope failures that are distinguishable from features caused by rilling and sheet?ow. Evidence for these failures was observed after a ?re near Santa Barbara, California. The amount of sediment eroded from some hillslopes was substantial, with 290 kg of sediment per metre width of hillslope delivered to the valley ?oor. The mechanics of these failures are examined with a numerical model that incorporates a stability analysis with subsurface ?ow routing along a typical hillslope pro?le. The model correctly predicts the location of the failures as well as the rainfall amount necessary to trigger them. Copyright © 2003 John Wiley & Sons, Ltd.     

The effect of accelerated soil erosion on hillslope morphology

Intensive agricultural land use can have detrimental effects on landscape properties, greatly accelerating soil erosion, with consequent fertility loss and reduced agricultural potential. To quantify the effects of such erosional processes on hillslope morphology and gain insight into the underlying dynamics, we use a twofold approach. First, a statistical analysis of topographical features is conducted, with a focus on slope and gradient distributions. The accelerated soil erosion is shown to be fingerprinted in the distribution tails, which provide a clear statistical signature of this human-induced land modification. Theoretical solutions are then derived for the hillslope morphology and the associated creep and runoff erosion fluxes, allowing us to distinguish between the main erosional mechanisms operating in disturbed and undisturbed areas. We focus our application on the landscape at the Calhoun Critical Zone Observatory in the US Southern Piedmont, where severe soil erosion followed intensive cotton cultivation, resulting in highly eroded and gullied hillslopes. The observed differences in hillslope morphologies in disturbed and undisturbed areas are shown to be related to the disruption of the natural balance between soil creep and runoff erosion. The relaxation time required for the disturbed hillslopes to reach a quasi-equilibrium condition is also investigated. © 2019 John Wiley & Sons, Ltd.     

Gully catchments as a sediment sink,not just a source: Results from a long‐term (~12 500 year) sediment budget

Sediment delivery from hillslopes to trunk streams represents a significant pathway of mass transfer in the landscape, with a large fraction facilitated by gully systems. The internal gully geomorphic dynamics represent a considerable gap in many landscape and empirical erosion models, therefore a better understanding of these processes over longer timescales (10–10⁴ years) is needed. This study analyses the sediment mass balance and storage dynamics within a headwater gully catchment in central Europe over the last ~12 500 years. Human induced erosion resulted in hillslope erosion rates ~2.3 times higher than under naturally de‐vegetated conditions (during the Younger Dryas), however the total sediment inputs to the gully system (and therefore gully aggradation), were similar. Net gully storage has consistently increased to become the second largest term in the sediment budget after hillslope erosion (storage is ~45% and ~73% of inputs during two separate erosion and aggradation cycles). In terms of the depletion of gully sediment storage, the sediment mass balance shows that export beyond the gully fan was not significant until the last ~500 years, due to reduced gully fan accommodation space. The significance of storage effects on the gully sediment mass balance, particularly the export terms, means that it would be difficult to determine the influences of human impact and/or climatic changes from floodplain or lake sedimentary archives alone and that the sediment budgets of the headwater catchments from which they drain are more likely to provide these mechanistic links. Copyright © 2015 John Wiley & Sons, Ltd.     

Linking landscape morphological complexity and sediment connectivity

Connectivity relates to the coupling of landforms (e.g. hillslopes and channels) and the transfer of water and sediment between them. The degree to which parts of a catchment are connected depends largely on the morphological complexity of the catchment's landscape. Landscapes can have very different and distinct morphologies, such as terraces, V‐shaped valleys or broad floodplains. The objective of this study is to better understand and quantify the relation between landscape complexity and catchment connectivity. We hypothesize that connectivity decreases with increasing landscape morphological complexity. To quantify the connectivity–complexity relationship virtual digital elevation models (DEMs) with distinct morphologies were used as inputs into the landscape evolution model LAPSUS to simulate the sediment connectivity of each landscape. Additionally, the hypothesis was tested on six common real DEMs with widely different morphologies. Finally, the effects of different rainfall time series on catchment response were explored. Simulation results confirm the hypothesis and quantify the non‐linear relation. Results from the exploration of sediment connectivity in response to sequences of rainfall events indicate that feedback between erosion and deposition are more important for certain landscape morphologies than for others: for a given rainfall input, a more effective sediment connectivity and erosion response may be expected from rolling or V‐shaped catchments than from dissected or stepped landscapes. Awareness of the differences in the behaviour and response of different morphologies to catchment processes provides valuable information for the effective management of landscapes and ecosystems through efficiently designed soil and water conservation measures. Copyright © 2013 John Wiley & Sons, Ltd.     

Fire effects on rangeland hydrology and erosion in a steep sagebrush‐dominated landscape

Post‐fire runoff and erosion from wildlands has been well researched, but few studies have researched the degree of control exerted by fire on rangeland hydrology and erosion processes. Furthermore, the spatial continuity and temporal persistence of wildfire impacts on rangeland hydrology and erosion are not well understood. Small‐plot rainfall and concentrated flow simulations were applied to unburned and severely burned hillslopes to determine the spatial continuity and persistence of fire‐induced impacts on runoff and erosion by interrill and rill processes on steep sagebrush‐dominated sites. Runoff and erosion were measured immediately following and each of 3 years post‐wildfire. Spatial and temporal variability in post‐fire hydrologic and erosional responses were compared with runoff and erosion measured under unburned conditions. Results from interrill simulations indicate fire‐induced impacts were predominantly on coppice microsites and that fire influenced interrill sediment yield more than runoff. Interrill runoff was nearly unchanged by burning, but 3‐year cumulative interrill sediment yield on burned hillslopes (50 g m^?2) was twice that of unburned hillslopes (25 g m^?2). The greatest impact of fire was on the dynamics of runoff once overland flow began. Reduced ground cover on burned hillslopes allowed overland flow to concentrate into rills. The 3‐year cumulative runoff from concentrated flow simulations on burned hillslopes (298 l) was nearly 20 times that measured on unburned hillslopes (16 l). The 3‐year cumulative sediment yield from concentrated flow on burned and unburned hillslopes was 20 400 g m^?2 and 6 g m^?2 respectively. Fire effects on runoff generation and sediment were greatly reduced, but remained, 3 years post‐fire. The results indicate that the impacts of fire on runoff and erosion from severely burned steep sagebrush landscapes vary significantly by microsite and process, exhibiting seasonal fluctuation in degree, and that fire‐induced increases in runoff and erosion may require more than 3 years to return to background levels. Published in 2008 by John Wiley & Sons, Ltd.     

Spatial discretization of large watersheds and its influence on the estimation of hillslope sediment yield

The combined use of water erosion models and geographic information systems has facilitated soil loss estimation at the watershed scale. Tools such as the Geo‐spatial interface for the Water Erosion Prediction Project (GeoWEPP) model provide a convenient spatially distributed soil loss estimate but require discretization to identify hillslopes and channels. In GeoWEPP, the TOpographic PArameteriZation (TOPAZ) model is used as an automated procedure to extract a watershed boundary, hillslopes and channels from a digital elevation model (DEM). Previous studies in small watersheds have shown that the size of the hillslopes and the channel distribution affect the model estimates, but in large watersheds, the effects on the soil loss estimates have yet to be tested. Therefore, the objective of this study was to evaluate the effect of discretization on the hillslope sediment yield estimates using GeoWEPP in two large watersheds (>10 km²). The watersheds were selected and discretized varying the TOPAZ parameters [critical source area (CSA) and minimum source channel length (MSCL)] in a 30‐m resolution digital elevation model. The drainage networks built with TOPAZ were compared with each other using the drainage density index. The results showed that the discretization affected hillslope sediment yield estimates and their spatial distribution more than the total runoff. The drainage density index and the hillslope sediment yield were proportional but inversely related; thus, soil loss estimates were highly affected by the spatial discretization. As a result of this analysis, a method to choose the CSA and MSCL values that generates the greatest fraction of hillslopes having profile lengths less than 200 m was developed. This slope length condition is particularly crucial when using the WEPP and GeoWEPP models, in order for them to produce realistic estimates of sheet and rill erosion. Finally, and as a result of this analysis, a more reliable method was developed for selecting the TOPAZ channel network parameters (CSA and MSCL). Copyright © 2015 John Wiley & Sons, Ltd.