Landscape evolution and glaciation of the Rwenzori Mountains, Uganda: Insights from numerical modeling

The Rwenzori Mountains are a high alpine mountain chain, about 40 × 80 km in size, just north of the equator in the western branch of the East African Rift System in Africa. The central part of the mountain chain is located in Uganda, and the highest peak, the Margherita Peak with 5119 m, lies on the border to the Democratic Republic of Congo. Topography is very pronounced, with steeply incised valleys and clear glacial landforms in the upper part of the mountain chain. The Rwenzori Mountains are an unusually high mountain chain located in the extensional setting of the East African Rift System, and the large elevation poses a challenging problem for geodynamists to explain.We have used the landscape evolution model ULTIMA THULE, which combines hillslope diffusion, fluvial erosion, and glacial abrasion and is driven by a climate driver, simulating the variations in temperature, precipitation, and relief over several glacial cycles. With a simulation time of 800 ka, we test the hypothesis of climate-tectonic interactions on the uplift of the Rwenzori Mountains.Our results show that a moderate cooling of around 6° causes widespread glaciation of the high mountain regions as observed during the peak glacial phases, and that morphological processes degrading the landscape allow for a tectonic uplift rate of around 0.5 mm a^− 1.     

Comparison of conceptual landscape metrics to define hillslope-scale sediment delivery ratio

The aim of this study was to evaluate four metrics to define the spatially variable (regionalised) hillslope sediment delivery ratio (HSDR). A catchment model that accounted for gully and streambank erosion and floodplain deposition was used to isolate the effects of hillslope gross erosion and hillslope delivery from other landscape processes. The analysis was carried out at the subcatchment (~ 40 km²) and the cell scale (400 m²) in the Avon-Richardson catchment (3300 km²), south-east Australia. The four landscape metrics selected for the study were based on sediment travel time, sediment transport capacity, flux connectivity, and residence time. Model configurations with spatially-constant or regionalised HSDR were calibrated against sediment yield measured at five gauging stations. The impact of using regionalised HSDR was evaluated in terms of improved model performance against measured sediment yields in a nested monitoring network, the complexity and data requirements of the metric, and the resulting spatial relationship between hillslope erosion and landscape factors in the catchment and along hillslope transects. The introduction of a regionalised HSDR generally improved model predictions of specific sediment yields at the subcatchment scale, increasing model efficiency from 0.48 to > 0.6 in the best cases. However, the introduction of regionalised HSDR metrics at the cell scale did not improve model performance. The flux connectivity was the most promising metric because it showed the largest improvement in predicting specific sediment yields, was easy to implement, was scale-independent and its formulation was consistent with sedimentological connectivity concepts. These properties make the flux connectivity metric preferable for applications to catchments where climatic conditions can be considered homogeneous, i.e. in small-medium sized basins (up to approximately 3000 km² for Australian conditions, with the Avon-Richardson catchment being at the upper boundary). The residence time metric improved model assessment of sediment yields and enabled accounting for climatic variability on sediment delivery, but at the cost of greater complexity and data requirements; this metric might be more suitable for application in catchments with important climatic gradients, i.e. large basins and at the regional scale. The application of a regionalised HSDR metric did not increase data or computational requirements substantially, and is recommended to improve assessment of hillslope erosion in empirical, semi-lumped erosion modelling applications. However, more research is needed to assess the quality of spatial patterns of erosion depicted by the different landscape metrics.     

RUNOFF PRODUCTION AND EROSION PROCESSES ON A DEHESA IN WESTERN SPAIN*

ABSTRACT. Runoff generation and soil erosion were investigated at the Guadalperalón experimental watershed (western Spain), within the land‐use system known as dehesa, or open, managed evergreen forests. Season and type of surface were found to control runoff and soil‐loss rates. Five soil units were selected as representative of surface types found in the study area: hillslope grass, bottom grass, tree cover, sheep trails, and shrub cover. Measurements were made in various conditions with simulated rainfall to gain an idea of the annual variation in runoff and soil loss. Important seasonal differences were noted due to surface cover and moisture content of soil, but erosion rates were determined primarily by runoff. Surfaces covered with grass and shrubs always showed less erosion; surfaces covered with holm oaks showed higher runoff rates, due to the hydrophobic character of the soils. Concentrations of runoff sediment during the simulations confirmed that erosion rates at the study site depended directly on the sediment available on the soil surface.     

Modelling the impact of land-use change and farm dam construction on hillslope sediment delivery to rivers at the regional scale

Landscapes in southeastern Australia have changed dramatically since the spread of European colonisation in the 19th century. Due to widespread forest clearance for cultivation and grazing, erosion and sediment yields have increased by a factor of more than 150. In the 20th century, erosion and sediment yield were reduced again due to an increasing vegetative cover. Furthermore, during the last decades, thousands of small farm dams were constructed to provide drinking water for cattle. These dams trap a lot of sediment, thereby further reducing sediment delivery from hillslopes to river channels. Changes in sediment delivery since European colonisation are documented in sediment archives. Within this study, these changing rates in hillslope erosion and sediment delivery were modelled using a spatially distributed erosion and sediment delivery model (WATEM/SEDEM) that was calibrated for Australian ecosystems using sediment yield data derived from sedimentation rates in 26 small farm dams. The model was applied to the Murrumbidgee river basin (30,000 km²) under different land-use scenarios. First, the erosion and sediment yield under pre-European land-use was modelled. Secondly, recent land-use patterns were used in the model. Finally, recent land-use including the impact of farm dams and large reservoirs was simulated. The results show that the WATEM/SEDEM model is capable of predicting the intensity of the geomorphic response to changes in land-use through time. Changes in hillslope erosion and hillslope sediment delivery rates are not equal, illustrating the non-linear response of the catchment. Current hillslope sediment supply to the river channel network is predicted to be 370% higher compared to the pre-European settlement period, yet farm dams have reduced this back to 2.5 times the pre-19th century values. The role of larger reservoirs is even more important as they have reduced the current sediment supply downstream to their pre-European values, thus completely masking the increased hillslope erosion rates from land-use change. However, the model does so far not include valley widening and sediment storage in river systems. Therefore, modelled rates of sediment delivery are lower than observed values.     

Soil-stratigraphic techniques in the study of soil and landform evolution across the Southern Alps, New Zealand

The Southern Alps lie along the convergent Pacific-Indian plate boundary. Geomorphically distinct eastern, axial and western regions reflect the east-west gradient in tectonic uplift (1 to 10 mm a⁻¹) and precipitation (600 to 10,000 mm a⁻¹). The eastern region is divided into front-ange and basin-and-range subregions. Soil-sequence studies on terraces established temporal contrasts in pedogenesis within and between eastern and western regions encompassing Entisols, Inceptisols and Spodosols. On Late Pleistocene and early Holocene terraces Dystrochrepts are persistent soils in the eastern region and Aquods in the western region. These soil sequences are used in the interpretation of relative soil age, stratigraphy and erosion history in hill and mountain drainage basins of the eastern and western regions. In the subhumid to humid eastern front-range subregion, simple soil forms occur as catenary sequences, and there is little evidence of erosion following the destruction of forests in the last millenium. Mollisols are dominant in the subhumid, and Dystrochrepts in humid areas, respectively. Soil-debris mantle regoliths date from the early Holocene and are still developing on slopes. The soil pattern on mountain slopes in the humid, eastern basin-and-range subregion is a complex array of simple, eroded, composite and compound soils. This pattern has resulted from erosion following forest destruction within the last millenium. The oldest surface or buried forest soils are Dystrochrepts dating from the Late Pleistocene to early Holocene. Wind erosion of these low-fertility soils contributes to the loessial sediments in which younger soils have formed. In the western region, soil patterns and soil stratigraphy indicate continous instability with a complex pattern of highly leached, shallow Orthents and bedrock outcrops on slopes. The soils are eroded from slopes within 2 ka. These contrasts in soil development and erosion periodicity in the eastern and western regions of the Southern Alps parallel the east-west contrasts in erosion rates of ca. 1–10 mm a⁻¹.     

A frost “buzzsaw” mechanism for erosion of the eastern Southern Alps, New Zealand

In the Southern Alps, New Zealand, large gradients in precipitation (< 1 to 12 m year^− 1) and rock uplift (< 1 to 10 mm year^− 1) produce distinct post-glacial geomorphic domains in which landslide-driven sediment production dominates in the wet, rapid-uplift western region, and rockfall controls erosion in the drier, low-uplift eastern region. Because the western region accounts for < 25% of the active orogen, the dynamics of erosion in the extensive eastern region are of equal importance in estimating the relative balance of uplift and erosion across the Southern Alps. Here, we assess the efficacy of frost cracking as the primary rockfall mechanism in the eastern Southern Alps using air photo and topographic analysis of scree slopes, cosmogenic radionuclide dating of headwalls, paleo-climate data, and a numerical model of headwall temperature. Currently, active scree slopes occur at a relatively uniform mean elevation ( 1450 m) and their distribution is independent of hillslope aspect and rock type, consistent with the notion that frost cracking (which is maximized between − 3 and − 8 °C) may control rockfall erosion. Headwall erosion rates of 0.3 to 0.9 mm year^− 1, measured using in-situ ¹⁰Be and ²⁶Al in the Cragieburn Range, confirm that rockfall erosion is active in the late Holocene at rates that roughly balance rock uplift. Models of the predicted depth of frost activity are consistent with the scale of fractures and scree blocks in our field sites. Also, vegetated, paleo-scree slopes are ubiquitous at elevations lower than active scree slopes, consistent with the notion that lower temperatures during the last glacial advance induced pervasive rockfall erosion due to frost cracking. Our modeling suggests temporally-averaged peak frost cracking intensity occurs at 2300 m a.s.l., the approximate elevation of the highest peaks in the central Southern Alps, suggesting that the height of these peaks may be limited by a “frost buzzsaw.”     

Infiltration on mountain slopes: a comparison of three environments

Water is well established as a major driver of the geomorphic change that eventually reduces mountains to lower relief landscapes. Nonetheless, within the altitudinal limits of continuous vegetation in humid climates, water is also an essential factor in slope stability. In this paper, we present results from field experiments to determine infiltration rates at forested sites in the Andes Mountains (Ecuador), the southern Appalachian Mountains (USA), and the Luquillo Mountains (Puerto Rico). Using a portable rainfall simulator–infiltrometer (all three areas), and a single ring infiltrometer (Andes), we determined infiltration rates, even on steep slopes. Based on these results, we examine the spatial variability of infiltration, the relationship of rainfall runoff and infiltration to landscape position, the influence of vegetation on infiltration rates on slopes, and the implications of this research for better understanding erosional processes and landscape change.Infiltration rates ranged from 6 to 206 mm/h on lower slopes of the Andes, 16 to 117 mm/h in the southern Appalachians, and 0 to 106 mm/h in the Luquillo Mountains. These rates exceed those of most natural rain events, confirming that surface runoff is rare in montane forests with deep soil/regolith mantles. On well-drained forested slopes and ridges, apparent steady-state infiltration may be controlled by the near-surface downslope movement of infiltrated water rather than by characteristics of the full vertical soil profile. With only two exceptions, the local variability of infiltration rates at the scale of 10° m overpowered other expected spatial relationships between infiltration, vegetation type, slope position, and soil factors. One exception was the significant difference between infiltration rates on alluvial versus upland soils in the Andean study area. The other exception was the significant difference between infiltration rates in topographic coves compared to other slope positions in the tabonuco forest of one watershed in the Luquillo Mountains. Our research provides additional evidence of the ability of forests and forest soils to preserve geomorphic features from denudation by surface erosion, documents the importance of subsurface flow in mountain forests, and supports the need for caution in extrapolating infiltration rates.     

Post-wildfire erosion response in two geologic terrains in the western USA

Volumes of eroded sediment after wildfires vary substantially throughout different geologic terrains across the western United States. These volumes are difficult to compare because they represent the response to rainstorms and runoff with different characteristics. However, by measuring the erosion response as the erodibility efficiency of water to detach and transport sediment on hillslopes and in channels, the erosion response from different geologic terrains can be compared. Specifically, the erodibility efficiency is the percentage of the total available stream power expended to detach, remobilize, or transport a mass of sediment. Erodibility efficiencies were calculated for the (i) initial detachment, and for the (ii) remobilization and transport of sediment on the hillslopes and in the channels after wildfire in two different geological terrains.The initial detachment efficiencies for the main channel and tributary channel in the granitic terrain were 10 ± 9% and 5 ± 4% and were similar to those for the volcanic terrain, which were 5 ± 5% and 1 ± 1%. No initial detachment efficiency could be measured for the hillslopes in the granitic terrain because hillslope measurements were started after the first major rainstorm. The initial detachment efficiency in the volcanic terrain was 1.3 ± 0.41%. The average remobilization and transport efficiencies associated with flash floods in the channels also were similar in the granitic (0.18 ± 0.57%) and volcanic (0.11 ± 0.41%) terrains. On the hillslope the remobilization and transport efficiency was greater in the volcanic terrain (2.4%) than in the granitic terrain (0.65%). However, this may reflect the reduced sediment availability after the first major rainstorm (30-min maximum rainfall intensity  90 mm h^− 1) in the granitic terrain, while easily erodible fine colluvium remained on the hillslope after the first rainstorm (30-min maximum rainfall intensity = 7.2 mm h^− 1) in the volcanic terrain. The erosion response in channels and on hillslopes of the granitic and volcanic terrains was similar when compared using erodibility efficiencies.     

Climatic influences on Holocene variations in soil erosion rates on a small hill in the Mojave Desert

A detailed study of a small hill in NE Mojave Desert in eastern California was conducted to elucidate the effect of climate on the variations in soil erosion rates through Holocene. Field surveys and sampling were carried out to obtain information on topography, geomorphology, soil and vegetation conditions, seismic refraction, sediment deposition, and hillslope processes. Integration of this information allowed reconstruction of the hill topography at the end of the Pleistocene, deduction of the evolution of the hill from the end of the Pleistocene to the present, and estimation of total soil losses resulting from various hillslope processes. The estimates are consistent with the premise that early Holocene climate change resulted in vegetation change, soil destabilization, and topographic roughening. Current, very slow, hillslope transport rates (e.g., 5 mm ky⁻¹ by rodent burrowing, a presently important transport form) appear inconsistent with the inferred total soil loss rate (31 mm ky⁻¹). Packrat midden studies imply that the NE Mojave Desert experienced enhanced monsoonal precipitation in the early Holocene, presumably accentuating soil loss. Water erosion on one slope of the hill was simulated using Water Erosion Prediction Project (WEPP), a process-based erosion model, using 4 and 6 ky of precipitation input compatible with an appropriate monsoonal climate and the present climate, respectively. The WEPP-predicted soil losses for the chosen slope were compatible with inferred soil losses. Identification of two time periods within the Holocene with distinct erosion characteristics may provide new insight into the current state of Mojave Desert landform evolution.     

A model for post‐orogenic development of a mountain range and its foreland

Decaying mountain ranges often show a surprisingly dynamic pattern of landscape evolution. Although one might expect a simple, monotonic decline in relief over time, evidence from several inactive mountain ranges shows alternating sequences of deposition and erosion in the associated basins, suggesting variations in relief and exhumation rate in the ranges themselves. Examples include the Southern Rocky Mountains, the Pyrenees, the European Alps and the Atlas Mountains. In this paper, we explore the possible origins of post‐orogenic landscape dynamics using a simple mathematical model of a mountain range and an adjacent foreland basin. The analysis highlights the importance of mass balance. In particular, a switch from basin exhumation to renewed sedimentation requires either an increase in sediment influx from the range or a decrease in sediment outflux beyond the basin margin. Although it is widely understood that post‐orogenic changes in erosion and sediment flux can have multiple causes (including climate change, regional tectonic uplift or tilting, or exhumation of variable lithologies), an important implication of our analysis is that the impact of such changes must differ in sign or magnitude between the range and the basin to be recorded. This requirement places an important constraint on viable explanations for alternating sequences of deposition and erosion in a decaying mountain‐basin pair.     

Denudation rates and a topography-driven rainfall threshold in northern Chile: Multiple cosmogenic nuclide data and sediment yield budgets

The quantification of geomorphic process rates on the outcrop- and the orogen-scale is important to describe accurately the interaction between the relative effects of erosion, tectonics and climate on landscape evolution. We report single and paired cosmogenic nuclide (¹⁰Be, ²⁶Al and ²¹Ne) derived erosion rates and exposure ages on hillslope interfluves from the tectonically active western central Andes that show a distinct spatial variation. A positive correlation of erosion rates with elevation and present-day rainfall rates is observed. Erosion rates at lower altitudes–the hyperarid Coastal Cordillera and the Western Escarpment with the northern part of the Atacama Desert–are extremely low and of the order of 10–100 cm/My (nominal exposure ages 1–6 My). In contrast, erosion rates at higher altitudes–the semiarid Western Cordillera–range up to 4600 cm/My (nominal exposure ages 0.02–0.1 My). This latter average long-term bedrock erosion rate record, suggested to be coupled to an orographically controlled pattern of rainfall, is also reflected in the pattern of denudation rates derived from a short-term decadal record of limited sediment yield data. Specifically, denudation rates calculated from sediment flux data are of a similar order of magnitude as erosion rates deduced from long-lived cosmogenic nuclides from bedrock hillslope interfluves of the Western Cordillera. Nevertheless, the production and the supply of sediment from the western Andean slope are very limited.Analysis of multiple cosmogenic nuclides allows simultaneous determination of erosion rates and exposure ages but also reveals complex exposure histories of non-bedrock samples, such as boulders or amalgamated clast samples. Notably, this study shows that saturation of nuclides, usually assumed in studies where only a single nuclide is analyzed, is rather the exception than the rule, as revealed by erosion island plots. Constant erosion that started much later than the formation age of the rocks or episodic erosion by spalling can partially explain non-steady-state concentrations and more complicated exposure scenarios. Furthermore, the use of multiple nuclides with different half-lives allowed us to infer that no significant variations in long-term erosion rates have occurred and that at the Western Escarpment erosion rates have been low and constant for most of the late Neogene. Nevertheless, the time intervals necessary to reach steady-state concentrations for cosmogenic nuclides can be quite different from those needed for landscapes to reach steady state.     

Spatial relationship between topography and rock uplift patterns in asymmetric mountain ranges based on a stream erosion model

The spatial relationship between topography and rock uplift patterns in asymmetric mountain ranges was investigated using a stream erosion model in which the asymmetric rock uplift was given and erosion rates were proportional to the m-th power of the drainage area and the n-th power of the channel gradient. The model conditions were simple, and thus the effects of horizontal rock movement, diffusional processes, and erosion thresholds were neglected, and spatially uniform precipitation, lithology, and vegetation were assumed. In asymmetric mountain ranges, under realistic exponent conditions (m < n) and the above assumptions, the surface erosion rate is faster on the steeper side and slower on the gentler side. The topographic axis migrates away from the rock uplift axis toward the center of the mountain range owing to the contrast in erosion rates. This migration continues until the erosion is balanced with rock uplift. In a dynamic steady state, the topographic pattern is independent of the rock uplift rate as indicated by an analytical solution, and is prescribed by the rock uplift pattern and the exponents m and n. As the asymmetry of the rock uplift pattern increases, the topographic axis migrates a greater distance. The location of the topographic axis is related to the location of the rock uplift axis by a simple logarithmic function, for a wide range of m and n. The fit of the numerical results and the logarithmic function is particularly good when m = 0.5 and n = 1.0. If the rock uplift pattern in asymmetric mountain ranges is known, the value of n − 5m/4 can be constrained based on the logarithmic relation, assuming a dynamic steady state. On the other hand, if the value of n − 5m/4 is known in an asymmetric mountain range, the rock uplift pattern can be estimated directly from the topography. This relation was applied to the Suzuka Range in central Japan, and the value of n − 5m/4 was estimated for an assumed reverse fault motion.     

Effect of grass coverage on sediment yield in the hillslope-gully side erosion system

By scouring experiments, the changeable process and characteristics of sediment yield in the hillslope-gully side erosion system with different coverage degrees and spatial locations of grass were studied. Five grass coverage degrees of 0, 30%, 50%, 70%, 90%, three spatial locations of grass (upslope, mid-slope, low-slope) and two water inflow rates of 3.2 L/min, 5.2 L/min were applied to a 0.5 by 7 m soil bed in scouring experiments. Results showed that the sediment yield decreased with the increase of grass coverage degree at 3.2 L/min water inflow rate in scouring experiments and the sediment yield with different grass locations on the sloping surface was in the order of upper > middle > lower. At 5.2 L/min water inflow rate, the differences of sediment yield among various grass coverage degrees were increased, whereas the changeable tendency of sediment yield with different grass locations on the whole sloping surface was not very obvious. The proportion of sediment yield from the gully side increased in an exponential relationship with the increase of grass coverage degree. When the grass was located on the lower position of hillslope, the influence for accelerating gully erosion is the greatest.     

Patterns of rock fragment cover generated by tillage erosion

Intensively cultivated areas in the upper part of the Guadalentin catchment (southeast Spain) show a systematic spatial pattern of surface rock fragment cover (R_c). The objective of this paper is to quantify and to explain this spatial rock fragment cover pattern. Therefore, a map of an intensively cultivated area of 5 km² was digitised, and for each pixel total topographic curvature was calculated. Next, rock fragment cover was determined photographically at 35 sites with a range of total slope curvatures. A linear relation between total curvature and rock fragment cover was found, except for narrow concavities. It was hypothesised that this pattern can be explained by a significant net downslope movement of rock fragments and fine earth by tillage. The displacement distances of rock fragments by tillage with a duckfoot chisel were measured by monitoring the displacement of tracers (painted rock fragments and aluminium cubes) on 5 sites having different slopes. The rare of tillage erosion for one tillage pass with a duckfoot chisel, expressed by the diffusion constant (k), equals 282 kg/m for up and downslope tillage and only 139 kg/m for contour tillage. Nomograms indicate that mean denudation rates in almond groves due to tillage erosion (3 to 5 tillage passes per year) can easily amount to 1.5–2.6 mm/year for contour tillage and up to 3.6–5.9 mm/year for up- and downslope tillage for a field, 50 m long and having a slope of 20%. These figures are at least one order of magnitude larger than reported denudation rates caused by water erosion in similar environments. Hence tillage erosion contributes significantly to land degradation. The downslope soil flux induced by tillage not only causes considerable denudation on topographic convexities (hill tops and spurs) and upper field boundaries but also an important sediment accumulation in topographic concavities (hollows and valley bottoms) and at lower field boundaries. Kinetic sieving (i.e. the upward migration of rock fragments) by the tines of the duckfoot chisel also concentrates the largest rock fragments in the topsoil in such a way that a rock fragment mulch develops in narrow valleys and at the foot of the slopes. These results clearly indicate that tillage erosion is the main process responsible for the observed rock fragment cover pattern in the study area. Since the study area is representative for many parts of southern Spain where almond groves have expanded since 1970, the results have a wider application. They show to what extent intensive tillage of steep slopes has contributed to the increase in soil degradation, to changes in hillslope morphology (i.e. strong denudation of convexities, development of lynchets and rapid infilling of narrow valley bottoms) and to the development of rock fragment cover patterns which control the spatial variability of the hydrological and water erosion response within such landscapes.     

坡面草被覆盖对坡沟侵蚀产沙过程的影响

By scouring experiments, the changeable process and characteristics of sediment yield in the hillslope-gully side erosion system with different coverage degrees and spatial locations of grass were studied. Five grass coverage degrees of 0, 30%, 50%, 70%, 90%, three spatial locations of grass (upslope, mid-slope, low-slope) and two water inflow rates of 3.2 L/min, 5.2 L/min were applied to a 0.5 by 7 m soil bed in scouring experiments. Results showed that the sediment yield decreased with the increase of grass coverage degree at 3.2 L/min water inflow rate in scouring experiments and the sediment yield with different grass locations on the sloping surface was in the order of upper > middle > lower. At 5.2 L/min water inflow rate, the differences of sediment yield among various grass coverage degrees were increased, whereas the changeable tendency of sediment yield with different grass locations on the whole sloping surface was not very obvious. The proportion of sediment yield from the gully side increased in an exponential relationship with the increase of grass coverage degree. When the grass was located on the lower position of hillslope, the influence for accelerating gully erosion is the greatest. Foundation: National Basic Research Program of China, No.2007CB407201; National Key Technology R&D Program, No.2006BAB06B01-06; Science and Technique Development Foundation of YRIHR, No.200603 Author: Li Mian (1968–), Ph.D and Senior Engineer, specialized in soil erosion.     

Relief-rejuvenation and topographic length scales in a fluvial drainage basin, Napf area, Central Switzerland

This paper explores how, and to what extent, a phase of relief-rejuvenation modifies the mode of surface erosion in an approximately 63 km² drainage basin located at the northern border of the Swiss Alps (Luzern area). In the study area, the retreat of the Alpine glaciers at the end of the Last Glacial Maximum (LGM) caused base level to lower by approximately 80 m. The fluvial system adapted to the lowered base level by headward erosion. This is indicated by knickzones in the longitudinal stream profiles and by the continuous upstream narrowing of the width of the valley floor towards these knickzones. In the headwaters above these knickzones, processes are still to a significant extent controlled by the higher base level of the LGM. There, frequent exposure of bedrock in channels and especially on hillslopes implies that sediment flux is to a large extent limited by weathering rates. In the knickzones, however, exposure of bedrock in channels implies that sediment flux is supply-limited, and that erosion rates are controlled by stream power.The morphometric analysis reveals the existence of length scales in the topography that result from distinct geomorphic processes. Along the tributaries where the upstream sizes of the drainage basins exceed 100,000–200,000 m², the mode of sediment transport and erosion changes from predominantly hillslope processes (i.e., landsliding, creep of regolith, rock avalanches and to some extent debris flows) to processes in channels (fluvial processes and debris flows). This length scale reflects the minimum size of the contributing area for channelized processes to take over in the geomorphic development (i.e., threshold size of drainage basin). This threshold size depends on the ratio between production rates of sediment on hillslopes, and export rates of sediment by processes in channels. Consequently, in the headwaters, erosion rates and sediment flux, and hence landscape evolution rates, are to a large extent limited by weathering processes. In contrast, in the lower portion of the drainage basin that adjusts to the lowered base-level, rates of channelized erosion and relief formation are controlled mainly by stream power. Hence, this paper shows that base-level lowering, headward erosion and establishment of knickzones separate drainage basins in two segments with different controls on rates of surface erosion, sediment flux and relief formation.     

Soil erodibility and processes of water erosion on hillslope

The importance of the inherent resistance of soil to erosional processes, or soil erodibility, is generally recognized in hillslope and fluvial geomorphology, but the full implications of the dynamic soil properties that affect erodibility are seldom considered. In Canada, a wide spectrum of soils and erosional processes has stimulated much research related to soil erodibility. This paper aims to place this work in an international framework of research on water erosion processes, and to identify critical emerging research questions. It focuses particularly on experimental research on rill and interrill erosion using simulated rainfall and recently developed techniques that provide data at appropriate temporal and spatial scales, essential for event-based soil erosion prediction. Results show that many components of erosional response, such as partitioning between rill and interrill or surface and subsurface processes, threshold hydraulic conditions for rill incision, rill network configuration and hillslope sediment delivery, are strongly affected by spatially variable and temporally dynamic soil properties. This agrees with other recent studies, but contrasts markedly with long-held concepts of soil credibility as an essentially constant property for any soil type. Properties that determine erodibility, such as soil aggregation and shear strength, are strongly affected by climatic factors such as rainfall distribution and frost action, and show systematic seasonal variation. They can also change significantly over much shorter time scales with subtle variations in soil water conditions, organic composition, microbiological activity, age-hardening and the structural effect of applied stresses. Property changes between and during rainstorms can dramatically affect the incidence and intensity of rill and interrill erosion and, therefore, both short and long-term hillslope erosional response. Similar property changes, linked to climatic conditions, may also significantly influence the stability and resilience of plant species and vegetation systems. Full understanding of such changes is essential if current event-based soil erosion models such as WEPP and EUROSEM are to attain their full potential predictive precision. The complexity of the interacting processes involved may, however, ultimately make stochastic modelling more effective than physically based modelling in predicting hillslope response to erodibility dynamics.     

Fallout radionuclide tracers identify a switch in sediment sources and transport-limited sediment yield following wildfire in a eucalypt forest

Fire can alter sediment sources and transport rates in river basins, changing landforms and aquatic habitats and degrading downstream water quality. Variability in the response between environments, between fires, and with time since fire makes predicting the catchment-scale effect of individual fires difficult. This study applies the fallout radionuclides ¹³⁷Cs and ²¹⁰Pb_xs to trace the sources and transport of fine sediment through a river network following a wildfire of moderate to extreme severity in the 629-km² eucalypt-forested Nattai River water-supply catchment near Sydney, Australia. The tracer analysis showed that post-fire erosion caused a switch in fine (< 10 µm) sediment sources from 80% subsoil derived from gully and river bank erosion to 86% topsoil derived from hillslope surface erosion. The fine sediment phosphorus content increased 4–10 fold over pre-fire levels. Annual post-fire sediment yields estimated from suspended solids rating curves were 109–250 times higher than they would have been without fire. A large additional amount of sediment remained stored within the river network for at least four years, particularly in lower-gradient reaches. Analysis of a sediment core showed that surface erosion following a previous fire had supplied at least 29% of total catchment sediment yield over the past 36 years. It is concluded that wildfire can alter catchment sediment budgets in two ways. Firstly, a spatially-diffuse pulse of elevated erosion is associated with moderate or intense rainfall events in post-fire years. Secondly, pulses of elevated catchment sediment yield are driven by the timing and river sediment transport capacity of runoff events. Severe post-fire erosion and high interannual hydrologic variability can result in large sediment stores persisting within the river network for many years. Fallout radionuclide tracers are shown to be useful in quantifying fine sediment sources and transport dynamics following wildfire, and the contribution of wildfire to catchment sediment yield.     

Sediment budget analysis of slope–channel coupling and in-channel sediment storage in an upland catchment, southeastern Australia

Slope–channel coupling and in-channel sediment storage can be important factors that influence sediment delivery through catchments. Sediment budgets offer an appropriate means to assess the role of these factors by quantifying the various components in the catchment sediment transfer system. In this study a fine (< 63 µm) sediment budget was developed for a 1.64-km² gullied upland catchment in southeastern Australia. A process-based approach was adopted that involved detailed monitoring of hillslope and bank erosion, channel change, and suspended sediment output in conjunction with USLE-based hillslope erosion estimation and sediment source tracing using ¹³⁷Cs and ²¹⁰Pb_ex. The sediment budget developed from these datasets indicated channel banks accounted for an estimated 80% of total sediment inputs. Valley floor and in-channel sediment storage accounted for 53% of inputs, with the remaining 47% being discharged from the catchment outlet. Estimated hillslope sediment input to channels was low (5.7 t) for the study period compared to channel bank input (41.6 t). However an estimated 56% of eroded hillslope sediment reached channels, suggesting a greater level of coupling between the two subsystems than was apparent from comparison of sediment source inputs. Evidently the interpretation of variability in catchment sediment yield is largely dependent on the dynamics of sediment supply and storage in channels in response to patterns of rainfall and discharge. This was reflected in the sediment delivery ratios (SDR) for individual measurement intervals, which ranged from 1 to 153%. Bank sediment supply during low rainfall periods was reduced but ongoing from subaerial processes delivering sediment to channels, resulting in net accumulation on the channel bed with insufficient flow to transport this material to the catchment outlet. Following the higher flow period in spring of the first year of monitoring, the sediment supplied to channels during this interval was removed as well as an estimated 72% of the sediment accumulated on the channel bed since the start of the study period. Given the seasonal and drought-dependent variability in storage and delivery, the period of monitoring may have an important influence on the overall SDR. On the basis of these findings, this study highlights the potential importance of sediment dynamics in channels for determining contemporary sediment yields from small gullied upland catchments in southeastern Australia.     

Coupling relationships: Hillslope–fluvial linkages in the Hodder catchment, NW England

Variations in the coupling of sediment transfer between different parts of a fluvial catchment, e.g., hillslope to axial stream, can hamper understanding but are an integral part of the geomorphological record. Depositional environments respond to a combination of land use, climate, storms (floods), and autogenic conditioning. The distribution of sediment in the upland landscapes of NW England is out of equilibrium with contemporary climate and geomorphological processes; more a function of peri- and paraglacial mobilisation of glacigenic deposits. Soil and vegetation development after deglaciation have interrupted any progression toward sediment exhaustion with sediment release controlled largely by extrinsic perturbation, with late Holocene anthropogenic activity, climate and extreme hydrological events the likely candidates. This paper presents a new radiocarbon-dated Holocene geomorphological succession for the River Hodder (NW England), alongside evaluating new palaeoecological and geoarchaeological data to discern the impacts of human activity. These data show a late Holocene expansion in human occupation and use of the landscape since the Iron Age (700–0 cal. B.C.), with more substantial changes in the character and intensity of upland land use in the last 1300 years. The geomorphological responses in the uplands were the onset of considerable and widespread hillslope erosion (gullying) and associated alluvial fan development. Interpretation of the regional radiocarbon chronology limits gullying to four, more extensive and aggressive phases after 500 cal. B.C. The downstream alluvial system has responded with considerable valley floor deposition and lateral channel migration that augmented sediment supply by remobilising the existing floodplain terraces and led to the aggradation of a series of inset alluvial terraces. The timing of these changes between states of aggradation and incision in alluvial reaches reflects the increased connectivity between the hillslope and alluvial systems. Aspects of both the regional climate and land use histories are conducive to increasing discharge and sediment flux, but the region wide lowering of erosion thresholds appears a key driver conditioning these sediment-rich conditions and producing a landscape that was more susceptible to erosion under lower magnitude flows.