Theoretical cosmogenic nuclide concentration in river bed load clasts: Does it depend on clast size?

Terrestrial Cosmogenic Nuclides (TCNs) have been widely used to date the exposure of alluvial surfaces and to estimate catchment-scale erosion rates. However, TCN concentration differences in samples of different grain sizes remain to be fully understood. In order to explore the possibility that river processes generate such differences, we develop a numerical model to calculate along-stream clast-scale TCN concentrations. Using the hillslope model, there is a progressive detachment of successive clasts of specific sizes followed by their instantaneous fall into the river. In the river, transport velocity and TCN concentration evolution in a clast depend 1) on the probability of being trapped in the sediment mixing layer of the river or within an adjacent terrace; 2) on its size which decreases downstream by attrition. The size-dependent transport law corresponds to the partial transport state in a river. We model the distribution of TCN concentrations in different clast size fractions in the 0–5 cm radius range for catchments in steady-state erosion, and for catchments experiencing sedimentation.We propose that clast attrition tends to increase the variance of TCN concentrations of the small clast size fractions because these fractions incorporate initially big clasts that travelled a long distance in addition to small clasts contributed near the outlet. We obtained numerous clast size–TCN concentration correlations, positive or negative, the significance of which depends on the initial clast size distribution, hillslope erosion rate, river length and lithology. For an equilibrium catchment, even large, we found that the addition of TCN concentration acquired during river transport is negligible compared to TCN concentration acquired on a hillslope, although a clast size–TCN concentration relationship can result from or be modified by clast attrition. On the contrary, aggrading catchments may show a significant clast size-dependent TCN concentration increase during river transport. This may introduce a small bias in the TCN-derived catchment erosion rate, but it could be used positively to quantify the mean transport velocity of clasts of different sizes over thousands of years. In addition, the lack of correlation between TCN concentration and clast size does not imply that the mean transport velocity is the same for all clast size fractions. Overall, our study provides an alternative explanation for observed clast size-dependent TCN concentrations and brings to the fore the need for measuring TCN concentration in larger clast size fractions than is usually done. To see if the byproducts of abrasion dilute or increase the TCN concentration of sand, all products should be included in a future study.     

Improved discrimination of subglacial and periglacial erosion using 10Be concentration measurements in subglacial and supraglacial sediment load of the Bossons glacier (Mont Blanc massif,France)

Deciphering the complex interplays between climate, uplift and erosion is not straightforward and estimating present‐day erosion rates can provide useful insights. Glaciers are thought to be powerful erosional agents, but most published ‘glacial’ erosion rates combine periglacial, subglacial and proglacial erosion processes. Within a glaciated catchment, sediments found in subglacial streams originate either from glacial erosion of substratum or from the rock walls above the glacier that contribute to the supraglacial load. Terrestrial cosmogenic nuclides (TCN) are produced by interactions between cosmic ray particles and element targets at the surface of the Earth, but their concentration becomes negligible under 15 m of ice. Measuring TCN concentrations in quartz sand sampled in subglacial streams and in supraglacial channels is statistically compliant with stochastic processes (e.g. rockfalls) and may be used to discriminate subglacial and periglacial erosion. Results for two subglacial streams of the Bossons glacier (Mont Blanc massif, France) show that the proportion of sediments originating from glacially eroded bedrock is not constant: it varies from 50% to 90% (n = 6). The difference between the two streams is probably linked to the presence or absence of supraglacial channels and sinkholes, which are common features of alpine glaciers. Therefore, most of the published mean catchment glacial erosion rates should not be directly interpreted as subglacial erosion rates. In the case of catchments with efficient periglacial erosion and particularly rockfalls, the proportion of sediments in the subglacial stream originating from the supraglacial load could be considerable and the subglacial erosion rate overestimated. Here, we estimate warm‐based subglacial and periglacial erosion rates to be of the same order of magnitude: 0.39 ± 0.33 and 0.29 ± 0.17 mm a^?1, respectively. Copyright © 2015 John Wiley & Sons, Ltd.     

Interpreting erosion rates from cosmogenic radionuclide concentrations measured in rapidly eroding terrain

A combination of numerical analysis and ¹⁰Be concentrations measured in sediment samples from the high‐relief Torrente catchment, southern Spain, allows us to investigate the sampling requirements for determining erosion rates using cosmogenic nuclides in high‐relief, landslide‐dominated terrain. We use simple modelling to quantify the effect of particle spalling and/or landsliding on erosion rates determined using a cosmogenic in‐situ produced isotope. Analytical results show that the cosmogenic nuclide concentration of a surface experiencing regular detachment of a grain or block may be considered to be in steady state, and ‘in‐situ’ erosion rates estimated, when an appropriate number of spatially independent samples are amalgamated. We present equations that enable calculation of the number of bedrock samples that must be amalgamated for the estimation of mean erosion rates on an outcrop experiencing regular detachment of a grain or chip of thickness L every T years. Our findings confirm that mean catchment erosion rates may be reliably estimated from ¹⁰Be concentrations in fluvial sediment in high‐relief rapidly eroding terrain. These catchment‐wide integrated erosion rates can be calculated where erosion is primarily accomplished through shallow (<3 m) spalling processes; where deep‐seated (>3 m) landslides are the dominant mode of erosion only minimum erosion rates can be determined. Lastly, we present erosion rate measurements from the Torrente catchment that reveal variation of two orders of magnitude (0·03–1·6 m ka^?1) quantifying the high degree of spatial variation in erosion rates expected within rapidly uplifting catchments. Copyright © 2006 John Wiley & Sons, Ltd.     

Very slow erosion rates and landscape preservation across the southwestern slope of the Ladakh Range,India

Erosion rates are key to quantifying the timescales over which different topographic and geomorphic domains develop in mountain landscapes. Geomorphic and terrestrial cosmogenic nuclide (TCN) methods were used to determine erosion rates of the arid, tectonically quiescent Ladakh Range, northern India. Five different geomorphic domains are identified and erosion rates are determined for three of the domains using TCN ¹⁰Be concentrations. Along the range divide between 5600 and 5700 m above sea level (asl), bedrock tors in the periglacial domain are eroding at 5.0 ± 0.5 to 13.1 ± 1.2 meters per million years (m/m.y.)., principally by frost shattering. At lower elevation in the unglaciated domain, erosion rates for tributary catchments vary between 0.8 ± 0.1 and 2.0 ± 0.3 m/m.y. Bedrock along interfluvial ridge crests between 3900 and 5100 m asl that separate these tributary catchments yield erosion rates <0.7 ± 0.1 m/m.y. and the dominant form of bedrock erosion is chemical weathering and grusification. Erosion rates are fastest where glaciers conditioned hillslopes above 5100 m asl by over‐steepening slopes and glacial debris is being evacuated by the fluvial network. For range divide tors, the long‐term duration of the erosion rate is considered to be 40–120 ky. By evaluating measured ¹⁰Be concentrations in tors along a model ¹⁰Be production curve, an average of ~24 cm is lost instantaneously every ~40 ky. Small (<4 km²) unglaciated tributary catchments and their interfluve bedrock have received very little precipitation since ~300 ka and the long‐term duration of their erosion rates is 300–750 ky and >850 ky, respectively. These results highlight the persistence of very slow erosion in different geomorphic domains across the southwestern slope of the Ladakh Range, which on the scale of the orogen records spatial changes in the locus of deformation and the development of an orogenic rain shadow north of the Greater Himalaya. Copyright © 2014 John Wiley & Sons, Ltd.     

Constant denudation rates in a high alpine catchment for the last 6 kyrs

Terrestrial cosmogenic nuclides (TCN) have widely been used as proxies in determining denudation rates in catchments. Most studies were limited to samples from modern active streams, thus little is known about the magnitude and causes of TCN variability on millennial time scales. In this work we present a 6 kyrs long, high resolution record of ¹⁰Be concentrations (n = 18), which were measured in sediment cores from an alluvial fan delta at the outlet of the Fedoz Valley in the Swiss Alps. This record is paired with a 3‐year time series (n = 4) of ¹⁰Be measured in sediment from the active stream currently feeding this fan delta. The temporal trend in the ¹⁰Be concentrations after correction for postdepositional production of ¹⁰Be was found to be overall constant and in good agreement with the modern river ¹⁰Be concentration. The calculated mean catchment‐wide denudation rate amounts to 0.73 ± 0.18 mm yr^?1. This fairly constant level of ¹⁰Be concentrations can be caused by a constant denudation rate over time within the catchment or alternatively by a buffered signal. In this contribution we suggest that the large alluvial floodplain in the Fedoz Valley may act as an efficient buffer on Holocene time scales in which sediments with different ¹⁰Be signatures are mixed. Therefore, presumable variations in the ¹⁰Be signals derived from changes in denudation under a fluctuating Holocene climate are only poorly transferred to the catchment outlet and not recorded in the ¹⁰Be record. However, despite the absence of high frequency signals, we propose that the buffered and averaged ¹⁰Be signal could be meaningfully and faithfully interpreted in terms of long‐term catchment‐averaged denudation rate. Our study suggests that alluvial buffers play an important role in regulating the ¹⁰Be signal exported by some alpine settings that needs to be taken into account and further investigated. Copyright © 2016 John Wiley & Sons, Ltd.     

Shoreline erosion: a cautionary note in using small farm dams to determine catchment erosion rates

Data from 10 small farm dams in SE Australia show that shoreline erosion due to farm livestock access to the dams can account for a significant proportion (up to 85%) of sediment contained in the dam. The volumes of sediment resulting from such shoreline erosion may be of the same order as the volumes produced by gully erosion in the dams' catchments, prompting caution in using farm dams to which livestock have access to determine small catchment erosion rates. Other issues, related to the trap efficiency, also mean that erosion estimates based on farm dam sedimentation should be treated with caution. © 1998 John Wiley & Sons, Ltd.     

Dynamics of soil erosion rates and controlling factors in the Northern Ethiopian Highlands – towards a sediment budget

This paper analyses the factors that control rates and extent of soil erosion processes in the 199 ha May Zegzeg catchment near Hagere Selam in the Tigray Highlands (Northern Ethiopia). This catchment, characterized by high elevations (2100–2650 m a.s.l.) and a subhorizontal structural relief, is typical for the Northern Ethiopian Highlands. Soil loss rates due to various erosion processes, as well as sediment yield rates and rates of sediment deposition within the catchment (essentially induced by recent soil conservation activities), were measured using a range of geomorphological methods. The area‐weighted average rate of soil erosion by water in the catchment, measured over four years (1998–2001), is 14·8 t ha^?1 y^?1, which accounts for 98% of the change in potential energy of the landscape. Considering these soil loss rates by water, 28% is due to gully erosion. Other geomorphic processes, such as tillage erosion and rock fragment displacement by gravity and livestock trampling, are also important, either within certain land units, or for their impact on agricultural productivity. Estimated mean sediment deposition rate within the catchment equals 9·2 t ha^?1 y^?1. Calculated sediment yield (5·6 t ha^?1 y^?1) is similar to sediment yield measured in nearby catchments. Seventy‐four percent of total soil loss by sheet and rill erosion is trapped in exclosures and behind stone bunds. The anthropogenic factor is dominant in controlling present‐day erosion processes in the Northern Ethiopian Highlands. Human activities have led to an overall increase in erosion process intensities, but, through targeted interventions, rural society is now well on the way to control and reverse the degradation processes, as can be demonstrated through the sediment budget. Copyright © 2007 John Wiley & Sons, Ltd.     

Complex multiple cosmogenic nuclide concentration and histories in the arid Rio Lluta catchment,northern Chile

Terrestrial cosmogenic nuclide (TCN) concentrations measured in river sediments can be used to estimate catchment‐wide denudation rates. By investigating multiple TCN the steadiness of sediment generation, transport and depositional processes can be tested. Measurements of ¹⁰Be, ²¹Ne and ²⁶Al from the hyper‐ to semi‐arid Rio Lluta catchment, northern Chile, yield average single denudation rates ranging from 12 to 75 m Myr^–1 throughout the catchment. Paired nuclide analysis reveals complex exposure histories for most of the samples and thus the single nuclide estimates do not exclusively represent catchment‐wide denudation rates. The lower range of single nuclide denudation rates (12–17 m Myr^–1), established with the noble gas ²¹Ne, is in accordance with palaeodenudation rates derived from ²¹Ne/¹⁰Be and ²⁶Al/¹⁰Be ratio analysis. Since this denudation rate range is measured throughout the system, it is suggested that a headwater signal is transported downstream but modulated by a complex admixture of sediment that has been stored and buried at proximal hillslope or terrace deposits, which are released during high discharge events. That is best evidenced by the stable nuclide ²¹Ne, which preserves the nuclide concentration even during storage intervals. The catchment‐wide single ²¹Ne denudation rates and the palaeodenuation rates contrast with previous TCN‐derived erosion rates from bedrock exposures at hillslope interfluves by being at least one order of magnitude higher, especially in the lower river course. These results support earlier studies that identified a coupling of erosional processes in the Western Cordillera contrasting with decoupled processes in the Western Escarpment and in the Coastal Cordillera. Copyright © 2008 John Wiley & Sons, Ltd.     

Impact of climate on landscape form,sediment transfer and the sedimentary record

The relationship between climate, landscape connectivity and sediment export from mountain ranges is key to understanding the propagation of erosion signals downstream into sedimentary basins. We explore the role of connectivity in modulating the composition of sediment exported from the Frontal Cordillera of the south-central Argentine Andes by comparing three adjacent and apparently similar semi-glaciated catchment-fan systems within the context of an along-strike precipitation gradient. We first identify that the bedrock exposed in the upper, previously glaciated reaches of the cordillera is under-represented in the lithological composition of gravels on each of three alluvial fans. There is little evidence for abrasion or preferential weathering of sediment sourced from the upper cordillera, suggesting that the observed bias can only be explained by sediment storage in these glacially widened and flattened valleys of the upper cordillera (as revealed by channel steepness mapping). A detailed analysis of the morphology of sedimentary deposits within the catchments reveals catchment-wide trends in either main valley incision or aggradation, linked to differences in hillslope–channel connectivity and precipitation. We observe that drier catchments have poor hillslope–channel connectivity and that gravels exported from dry catchments have a lithological composition depleted in clasts sourced from the upper cordillera. Conversely, the catchment with the highest maximum precipitation rate exhibits a high degree of connectivity between its sediment sources and the main river network, leading to the export of a greater proportion of upper cordillera gravel as well as a greater volume of sand. Finally, given a clear spatial correlation between the resistance of bedrock to erosion, mountain range elevation and its covariant, precipitation, we highlight how connectivity in these semi-glaciated landscapes can be preconditioned by the spatial distribution of bedrock lithology. These findings give insight into the extent to which sedimentary archives record source erosion patterns through time.     

Exploring suspended sediment delivery dynamics of two Mediterranean nested catchments

This paper investigates suspended sediment transport and dynamics of two nested agricultural lowland Mediterranean catchments with a difference of two orders of magnitude in the surface area (i.e., 1 and 264 km²). The effects of the drainage catchment area over the specific suspended sediment yield are assessed by using the nested approach over various timeframes. A detailed analysis of the rainfall–runoff–sediment transport relationships during the 2‐year study period shows that the hydrological and sedimentological responses were extremely variable for both catchments. Very low or no correlations were observed between the rainfall intensity and the selected hydrological variables and sediment loads. However, remarkable or high correlations were obtained between the rainfall intensity and the maximum and average suspended sediment concentrations, indicating that rainfall per unit time has little control on the hydrological response, but that, simultaneously, its high‐erosive power triggers sediment production, increasing the sedimentary response of the catchments. This study also illustrates how sediment is mainly transported during floods, producing predominantly clockwise hysteretic loops. Moreover, the small headwater catchment exerts a reduced (or even negligible) effect over the hydro‐sedimentary response of the larger downstream catchment, caused by the reduced sediment availability in a landscape with an inherent disconnection of the sediment pathways.     

The relationship between dissolved organic carbon in stream water and soil organic carbon pools at different spatial scales

The relationship between stream water DOC concentrations and soil organic C pools was investigated at a range of spatial scales in subcatchments of the River Dee system in north‐east Scotland. Catchment percentage peat cover and soil C pools, calculated using local, national and international soils databases, were related to mean DOC concentrations in streams draining small‐ (<5 km²), medium‐ (12–38 km²) and large‐scale (56–150 km²) catchments. The results show that, whilst soil C pool is a good predictor of stream water DOC concentration at all three scales, the strongest relationships were found in the small‐scale catchments. In addition, in both the small‐ and large‐scale catchments, percentage peat cover was as a good predictor of stream water DOC concentration as catchment soil C pool. The data also showed that, for a given soil C pool, streams draining lowland (<700 m) catchments had higher DOC concentrations than those draining upland (>700 m) catchments, suggesting that disturbance and land use may have a small effect on DOC concentration. Our results therefore suggest that the relationship between stream water DOC concentration and catchment soil C pools exists at a range of spatial scales and this relationship appears to be sufficiently robust to be used to predict the effects of changes in catchment soil C storage on stream water DOC concentration. Copyright © 1999 John Wiley & Sons, Ltd.     

Using 137Cs measurements to establish catchment sediment budgets and explore scale effects

Growing awareness of the wider environmental significance of fine sediment transport by rivers and associated sediment problems linked to sediment–water quality interactions, nutrient and contaminant transfer, and the degradation of aquatic habitats has resulted in the need for an improved understanding of the mobilization and transfer of sediment in catchments to support the development of effective sediment management strategies. The sediment budget provides a key integrating concept for assembling information on the internal functioning of a catchment in terms of its sediment dynamics by providing information on the mobilization, transfer, storage and output of sediment. One key feature of a catchment sediment budget is the relationship between the sediment yield at the catchment outlet and rates of sediment mobilization and transfer within the catchment, which is commonly represented by the sediment delivery ratio. To date, most attempts to derive estimates of this ratio have been based on a comparison of the measured sediment yield from a catchment with an estimate of the erosion occurring within the catchment, derived from an erosion prediction procedure, such as the Universal Soil Loss Equation (USLE) or its revised version, RUSLE. There is a need to obtain more direct and spatially distributed evidence of the erosion rates occurring within a catchment and to characterize the links between sediment mobilization, transfer, storage and output more explicitly. In this context, fallout radionuclides have proved particularly useful as sediment tracers. This paper reports the results of a study aimed at exploring the use of caesium‐137 (¹³⁷Cs) measurements to establish sediment budgets for three catchments of different sizes and contrasting land use located in Calabria, southern Italy. Long‐term measurements of sediment output were available for the catchments, and, by using the estimates of gross and net rates of soil loss within the catchments provided by ¹³⁷Cs measurements, it was possible to establish the key components of the sediment budget for each catchment. By documenting the sediment budgets of three catchments of different sizes, the study provides a basis for exploring the effects of scale on catchment sediment budgets and, in particular, the increasing importance of catchment storage as the size of the catchment increases. The results of this study demonstrate a reduction in the sediment delivery ratio from 98 to 2% as catchment area increases from 1·47 ha to 31·2 km². Copyright © 2010 John Wiley & Sons, Ltd.     

Improving bank erosion modelling at catchment scale by incorporating temporal and spatial variability

Bank erosion can contribute a significant portion of the sediment budget within temperate catchments, yet few catchment scale models include an explicit representation of bank erosion processes. Furthermore, representation is often simplistic resulting in an inability to capture realistic spatial and temporal variability in simulated bank erosion. In this study, the sediment component of the catchment scale model SHETRAN is developed to incorporate key factors influencing the spatio‐temporal rate of bank erosion, due to the effects of channel sinuosity and channel bank vegetation. The model is applied to the Eden catchment, north‐west England, and validated using data derived from a GIS methodology. The developed model simulates magnitudes of total catchment annual bank erosion (617–4063 t y^‐1) within the range of observed values (211–4426 t yr^‐1). In addition, the model provides both greater inter‐annual and spatial variability of bank eroded sediment generation when compared with the basic model, and indicates a potential 61% increase of bank eroded sediment as a result of temporal flood clustering. The approach developed within this study can be used within a number of distributed hydrologic models and has general applicability to temperate catchments, yet further development of model representation of bank erosion processes is required. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

The ‘humped’ soil production function: eroding Arnhem Land,Australia

We report erosion rates and processes, determined from in situ‐produced beryllium‐10 (¹⁰Be) and aluminum‐26 (²⁶Al), across a soil‐mantled landscape of Arnhem Land, northern Australia. Soil production rates peak under a soil thickness of about 35 cm and we observe no soil thicknesses between exposed bedrock and this thickness. These results thus quantify a well‐defined ‘humped’ soil‐production function, in contrast to functions reported for other landscapes. We compare this function to a previously reported exponential decline of soil production rates with increasing soil thickness across the passive margin exposed in the Bega Valley, south‐eastern Australia, and found remarkable similarities in rates. The critical difference in this work was that the Arnhem Land landscapes were either bedrock or mantled with soils greater than about 35 cm deep, with peak soil production rates of about 20 m/Ma under 35–40 cm of soil, thus supporting previous theory and modeling results for a humped soil production function. We also show how coupling point‐specific with catchment‐averaged erosion rate measurements lead to a better understanding of landscape denudation. Specifically, we report a nested sampling scheme where we quantify average erosion rates from the first‐order, upland catchments to the main, sixth‐order channel of Tin Camp Creek. The low (～5 m/Ma) rates from the main channel sediments reflect contributions from the slowly eroding stony highlands, while the channels draining our study area reflect local soil production rates (～10 m/Ma off the rocky ridge; ～20 m/Ma from the soil mantled regions). Quantifying such rates and processes help determine spatial variations of soil thickness as well as helping to predict the sustainability of the Earth's soil resource under different erosional regimes. Copyright © 2009 John Wiley & Sons, Ltd.     

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

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

Short‐term versus medium‐term monitoring for detecting gully‐erosion variability in a Mediterranean environment

This study investigates how medium‐term gully‐development data differ from short‐term data, and which factors influence their spatial and temporal variability at nine selected actively retreating bank gullies situated in four Spanish basin landscapes. Small‐format aerial photographs using unmanned, remote‐controlled platforms were taken at the gully sites in short‐term intervals of one to two years over medium‐term periods of seven to 13 years and gully change during each period was determined using stereophotogrammetry and a geographic information system. Results show a high variability of annual gully retreat rates both between gullies and between observation periods. The mean linear headcut retreat rates range between 0·02 and 0·26 m a^–1. Gully area loss was between 0·8 and 22 m² a^–1 and gully volume loss between 0·5 to 100 m³ a^–1, of which sidewall erosion may play a considerable part. A non‐linear relationship between catchment area and medium‐term gully headcut volume change was found for these gullies. The short‐term changes observed at the individual gullies show very high variability: on average, the maximum headcut volume change observed in 7–13 years was 14·3 times larger than the minimum change. Dependency on precipitation varies but is clearly higher for headcuts than sidewalls, especially in smaller and less disturbed catchments. The varying influences of land use and human activities with their positive or negative effects on runoff production and connectivity play a dominant role in these study areas, both for short‐term variability and medium‐term difference in gully development. The study proves the value of capturing spatially continuous, high‐resolution three‐dimensional data using small‐format aerial photography for detailed gully monitoring. Results confirm that short‐term data are not representative of longer‐term gully development and demonstrate the necessity for medium‐ to long‐term monitoring. However, short‐term data are still required to understand the processes – particularly human activity at varying time scales – causing fluctuations in gully erosion rates. Copyright © 2011 John Wiley & Sons, Ltd.     

Which catchment characteristics control the temporal dependence structure of daily river flows?

Hydrological classification systems seek to provide information about the dominant processes in the catchment to enable information to be transferred between catchments. Currently, there is no widely agreed‐upon system for classifying river catchments. This paper develops a novel approach to classifying catchments based on the temporal dependence structure of daily mean river flow time series, applied to 116 near‐natural ‘benchmark’ catchments in the UK. The classification system is validated using 49 independent catchments. Temporal dependence in river flow data is driven by the flow pathways, connectivity and storage within the catchment and can thus be used to assess the influence catchment characteristics have on moderating the precipitation‐to‐flow relationship. Semi‐variograms were computed for the 116 benchmark catchments to provide a robust and efficient way of characterising temporal dependence. Cluster analysis was performed on the semi‐variograms, resulting in four distinct clusters. The influence of a wide range of catchment characteristics on the semi‐variogram shape was investigated, including: elevation, land cover, physiographic characteristics, soil type and geology. Geology, depth to gleyed layer in soils, slope of the catchment and the percentage of arable land were significantly different between the clusters. These characteristics drive the temporal dependence structure by influencing the rate at which water moves through the catchment and/or the storage in the catchment. Quadratic discriminant analysis was used to show that a model with five catchment characteristics is able to predict the temporal dependence structure for un‐gauged catchments. This method could form the basis for future regionalisation strategies, as a way of transferring information on the precipitation‐to‐flow relationship between gauged and un‐gauged catchments. © 2014 The Authors. Hydrological Processes by published by John Wiley & Sons, Ltd.     

Modelling mean annual sediment yield using a distributed approach

In this paper a spatially distributed model for the calculation of sediment delivery to river channels is presented (SEDEM: SEdiment DElivery Model). The model consists of two components: (1) the calculation of a spatial pattern of mean annual soil erosion rates in the catchment using a RUSLE (Revised Soil Erosion Equation) approach; and (2) the routing of the eroded sediment to the river channel network taking into account the transport capacity of each spatial unit. If the amount of routed sediment exceeds the local transport capacity, sediment deposition occurs. An existing dataset on sediment yield for 24 catchments in central Belgium was used to calibrate the transport capacity parameters of the model. A validation of the model results shows that the sediment yield for small and medium sized catchments (10–5000 ha) can be predicted with an average accuracy of 41 per cent. The predicted sediment yield values with SEDEM are significantly more accurate than the predictions using a lumped regression model. Moreover a spatially distributed approach allows simulation of the effect of different land use scenarios and soil conservation techniques. Copyright © 2001 John Wiley & Sons, Ltd.     

Generation of realistic synthetic catchments to explore fine continental surface processes

Understanding, analysing, and predicting the erosion mechanisms and sedimentary flows produced by catchments plays a key role in environmental conservation and restoration management and policies. Numerical case-testing studies are generally undertaken to analyse the sensitivity of flood and soil erosion processes to the physical characteristics of catchments. Most analyses are conducted on simple virtual catchments with physical characteristics that, unlike real catchments, are perfectly controlled. Virtual catchments generally correspond to V-shaped valley catchments. However, although these catchments are suitable for methodical analysis of the results, they do not provide a realistic representation of the spatial structures of the landscape and field conditions. They can, therefore, lead to potential modelling errors and can make it difficult to extend or generalize their results. Our proposed method bridges the gap between real and traditional virtual catchments by creating realistic virtual catchments with perfectly controllable physical characteristics. Our approach represents a real alternative to traditional test case procedures and provides a new framework for geomorphological and hydrological communities. It combines a field procedural generation approach, geographic information system processing procedures, and the CAESAR-Lisflood landscape evolution model. We illustrate how each of these components acts in the process of generating virtual catchments. Five physical parameters were adjusted and tested for each virtual catchment: drainage density, hypsometric integral, mean slope of the main channel, granulometry, and land use. One of our virtual catchments is compared with a real catchment and a virtual catchment produced by a standard method. This comparison indicates that our approach can produce more realistic virtual catchments than those produced by more traditional methods, while a high degree of controllability is maintained. This new method of generating virtual catchments therefore offers significant research potential to identify the impacts of the physical characteristics of catchments on hydro-sedimentary dynamics and responses.     

Estimation of sub-annual inter-catchment groundwater flow using short-term water balance method

Predicting inter-catchment groundwater flow (IGF) is essential because IGF greatly affects stream water discharge and water chemistry. However, methods for estimating sub-annual IGF and clarifying its mechanisms using minimal data are limited. Thus, we quantified the sub-annual IGF and elucidated its driving factors using the short-term water balance method (STWB) for three forest headwater catchments in Japan (named here catchment A, B and As). Our previous study using the chloride mass balance indicated that annual IGF of catchment A (49.0 ha) can be negligible. Therefore, we calculated the daily evapotranspiration (ET) rate using the Priestley–Taylor expression and the 5-year water balance in catchment A (2010–2014). The sub-annual IGF of the three catchments was then calculated by subtracting the ET rate from the difference between rainfall and stream discharge during the sub-annual water balance periods selected using the STWB. The IGF rates of catchment B (7.0 ha), which is adjacent to catchment A, were positive in most cases, indicating that more groundwater flowed out of the catchment than into it, and exhibited positive linear relationships with rainfall and stream discharge. This suggested that as the catchments became wetter, more groundwater flowed out of catchment B. Conversely, the IGF rates of catchment As (5.3 ha), included in catchment A, were negative in most cases, indicating that more groundwater flowed into the catchment than out from it, and exhibited negative linear relationships with rainfall and stream discharge. Given the topography of the catchments studied, infiltration into the bedrock was the probable reason for the IGF outflow from catchment B. We hypothesized that in catchment As, the discrepancy between the actual hydrological boundary and the surface topographic boundary could have caused an IGF inflow. This study provides a useful tool for determining an IGF model structure to be incorporated into rainfall-runoff models.