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.     

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.     

EROSION PROCESSES AT THE HILLSLOPE SCALE

The results are presented of an intensive monitoring programme to determine eroded volumes due to observed rainfall events on a hillslope surface. The surface investigated has been reproduced through a digital elevation model and analysed in terms of drainage network, contributing area and slope by performing planimetric and altimetric analysis. Erosion maps were derived from a comparison of consecutive digital elevation models relative to time. These maps reveal the spatial and temporal evolution of the erosion process at the hillslope scale. The erosion was uniform across the surface, supporting the assumption of randomness in the erosion process commonly used in surface drainage development models. The observed value of erosion has been estimated at approximately 0.11 m/m² per year, with almost 500 mm of total annual rainfall.     

Subsurface flow responses of a small forested catchment in the Ouachita mountains

Despite considerable research performed on forested catchments in the Ouachita Mountains of Oklahoma and Arkansas, little information on hydrological processes in operation is available. Based on catchment physical characteristics, subsurface flow was thought to be an important hydrological process in the region. Therefore, this study was undertaken to determine the occurrence, rates, timing and volumes of subsurface flow, and to estimate the importance of subsurface flow as a streamflow generating process. Subsurface flow was collected from three hillslope sites on a 7.7 ha forested catchment. Hillslope sites drained through natural seepage faces located near stream channels. Subsurface flow was collected from three depths at each hillslope site, below the litter layer, below the a horizon, and within the B horizon (Bt21). Subsurface flow occurred and was measured during 11 of 31 rainfall events. Subsurface flow responded rapidly to the initiation of and to changes in intensity of rainfall at all depths. the rapid response was indicative of flow through soil macropores. B horizon subsurface flow commenced within 10 to 180 min of the initiation of rainfall. Multiple linear regression showed that the volume of subsurface flow generated during a given storm was directly related to rainfall depth and a 7-day antecedent precipitation index used to represent antecedent water content. About 67 per cent of the total subsurface flow collected during the study was produced in one large storm under wet antecedent conditions. the storm was equal to the 2-year, 24-hour storm for the region. Measured subsurface flow volumes were extended to the watershed scale to provide estimates of catchment-wide contributions to streamflow. It was estimated that subsurface flow contributed from 1 to 48 per cent of total quickflow measured at the catchment outlet. Based on the timing of subsurface flow, it was estimated that subsurface flow May, contribute up to 70 per cent of quickflow before and soon after peak flow.     

Simulation of Turbulent Flow in An Urban Forested Park Damaged by a Windstorm

The damage caused by windstorms to forest ecosystems is often very heterogeneous. In order to improve the stability of forested landscapes, it is of great importance to identify the factors responsible for this spatial variability. The structure of the landscape itself may play a role, through possible influences of canopy heterogeneities on the development of turbulence. For the purpose of investigating the role of landscape fragmentation on turbulence development, we used a numerical flow model with a k–ε turbulence scheme model, previously validated in simple cases with well-defined surface changes (roughness change and forest edge flow). A series of two- and three-dimensional simulations were performed over a heterogeneous urban forested park in Europe, which was severely damaged in various places by the Lothar windstorm in December 1999. The model shows the development of a region of strong turbulence, resulting from the generation of large wind shear at the top of the canopy. A sensitivity study shows how the location, extension and intensity of the region depend on canopy characteristics such as the leaf density, the nature of the edge or the presence of gaps and clearings. Simulations performed in conditions representative of the windstorm show that the location of the damaged areas corresponds very closely to the regions where the turbulent kinetic energy was above a certain threshold.     

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.     

Effects of hillslope geometry on spatial infiltration using the TOPMODEL and SCS-CN models

ABSTRACT

The estimation of infiltration is a main issue in runoff simulation. The geometry of hillslopes (plan shape and profile curvature) may affect the responses, as well as infiltration over the hillslopes. In this study, the equations of TOPMODEL (a topography-based model) were applied to complex hillslopes to develop the complex TOPMODEL. This model was coupled with the SCS-CN (Soil Conservation Service Curve Number) model to examine the effects of geometry on infiltration and derive a saturation excess-based curve number (CN). The effects of plan shape and profile curvature upon the spatial distribution of CN and infiltration were studied. The results show that convergent hillslopes have 15.4% less infiltration and divergent hillslopes have 7.8% more infiltration than parallel ones. The infiltration over concave hillslopes is 13.5% lower and infiltration over convex hillslopes 5.8% higher than for straight ones. The degree of convergence/divergence has a greater effect on the CN compared to that of profile curvature.