Modelling soil erosion with a downscaled landscape evolution model

The measurement and prediction of soil erosion is important for understanding both natural and disturbed landscape systems. In particular numerical models of soil erosion are important tools for managing landscapes as well as understanding how they have evolved over time. Over the last 40 years a variety of methods have been used to determine rates of soil loss from a landscape and these can be loosely categorized into empirical and physically based models. Alternatively, physically based landscape evolution models (LEMs) have been developed that provide information on soil erosion rates at much longer decadal or centennial scales, over large spatial scales and examine how they may respond to environmental and climatic changes. Both soil erosion LEMs are interested in similar outcomes (landscape development and sediment delivery) yet have quite different methodologies and parameterizations. This paper applies a LEM (the CAESAR model) for the first time at time and space scales where soil erosion models have largely been used. It tests the ability of the LEM to predict soil erosion on a 30 m experimental plot on a trial rehabilitated landform in the Northern Territory, Australia. It then continues to discuss the synergies and differences between soil erosion and LEMs. The results demonstrate that once calibrated for the site hydrology, predicted suspended sediment and bedload yields from CAESAR show a close correspondence in both volume and timing of field measured data. The model also predicts, at decadal scales, sediment loads close to that of field measured data. Findings indicate that the small‐scale drainage network that forms within these erosion plots is an important control on the timing and magnitude of sediment delivery. Therefore, it is important to use models that can alter the DEM to reflect changing topography and drainage network as well as having a greater emphasis on channel processes. Copyright © 2012 John Wiley & Sons, Ltd. and Commonwealth of Australia     

The impact of depression removal on catchment geomorphology,soil erosion and landscape evolution

The current generation of landscape evolution models use a digital elevation model for landscape representation. These programs also contain a hydrological model that defines overland flow with the drainage network routed to an outlet. One of the issues with landscape evolution modelling is the hydrological correctness of the digital elevation model used for the simulations. Despite the wide use and increased quality of digital elevation models, data pits and depressions in the elevation data are a common feature and their removal will remain a necessary step for many data sets. This study examines whether a digital elevation model can be hydrologically correct (i.e. all depressions removed so that all water can run downslope) before use in a landscape evolution model and what effect depression removal has on long‐term geomorphology and hydrology. The impact on sediment transport rates is also examined. The study was conducted using a field catchment and a proposed landform for a post‐mining landscape. The results show that there is little difference in catchment geomorphology and hydrology for the non‐depression removed and depression removed data sets. The non‐depression removed and depression removed digital elevation models were also evaluated as input to a landscape evolution model for a 50 000 year simulation period. The results show that after 1000 years there is little difference between the data sets, although sediment transport rates did vary considerably early on in the simulation. Copyright © 2007 John Wiley & Sons, Ltd.     

Linking landslides,hillslope erosion,and landscape evolution

This Virtual Issue highlights 10 recent innovative, unconventional, or otherwise significant contributions to Earth Surface Processes and Landforms that help advance the state‐of‐the‐art in research on linkages between landslides, hillslope erosion, and landscape evolution. The selected studies address this feedback within a temporal spectrum that ranges from the event to the millennial scale, thus underscoring the importance of detailed field observations, high‐resolution digital topographic data, and geochronological methods for increasing our capability of quantifying landslide processes and hillslope erosion. Copyright © 2009 John Wiley & Sons, Ltd.     

Quantifying erosional equilibrium across a slowly eroding,soil mantled landscape

The textbook concept of an equilibrium landscape, which posits that soil production and erosion are balanced and equal channel incision, is rarely quantified for natural systems. In contrast to mountainous, rapidly eroding terrain, low relief and slow-eroding landscapes are poorly studied despite being widespread and densely inhabited. We use three field sites along a climosequence in South Africa to quantify very slow (2-5 m/My) soil production rates that do not vary across hillslopes or with climate. We show these rates to be indistinguishable from spatially invariant catchment-average erosion rates while soil depth and chemical weathering increase strongly with rainfall across our sites. Our analyses imply landscape-scale equilibrium although the dominant means of denudation varies from physical weathering in dry climates to chemical weathering in wet climates. In the two wetter sites, chemical weathering is so significant that clay translocates both vertically in soil columns and horizontally down hillslope catenas, resulting in particle size variation and the accumulation of buried stone lines at the clay-rich depth. We infer hundred-thousand-year residence times of these stone lines and suggest that bioturbation by termites plays a key role in exhuming sediment into the mobile soil layer from significant depths below the clay layer. Our results suggest how tradeoffs in physical and chemical weathering, potentially modulated by biological processes, shape slowly eroding, equilibrium landscapes. © 2019 John Wiley & Sons, Ltd.     

Stochastic processes of soil production and transport: erosion rates,topographic variation and cosmogenic nuclides in the Oregon Coast Range

Landscapes in areas of active uplift and erosion can only remain soil‐mantled if the local production of soil equals or exceeds the local erosion rate. The soil production rate varies with soil depth, hence local variation in soil depth may provide clues about spatial variation in erosion rates. If uplift and the consequent erosion rates are sufficiently uniform in space and time, then there will be tendency toward equilibrium landforms shaped by the erosional processes. Soil mantle thickness would adjust such that soil production matched the erosion. Previous work in the Oregon Coast Range suggested that there may be a tendency locally toward equilibrium between hillslope erosion and sediment yield. Here results from a new methodology based on cosmogenic radionuclide accumulation in bedrock minerals at the base of the soil column are reported. We quantify how soil production varies with soil thickness in the southern Oregon Coast Range and explore further the issue of landscape equilibrium. Apparent soil production is determined to be an inverse exponential function of soil depth, with a maximum inferred production rate of 268 m Ma^?1 occurring under zero soil depth. This rate depends, however, on the degree of weathering of the underlying bedrock. The stochastic and large‐scale nature of soil production by biogenic processes leads to large temporal and spatial variations in soil depth; the spatial variation of soil depth neither supports nor rejects equilibrium morphology. Our observed catchment‐averaged erosion rate of 117 m Ma^?1 is, however, similar to that estimated for the region by others, and to soil production rates under thin and intermediate soils typical for the steep ridges. We suggest that portions of the Oregon Coast Range may be eroding at roughly the same rate, but that local competition between drainage networks and episodic erosional events leads to landforms that are out of equilibrium locally and have a spatially varying soil mantle. Copyright © 2001 John Wiley & Sons, Ltd.     

Hillslope erosion measurement—a simple approach to a complex process

The measurement of hillslope erosion can be a difficult, costly and time‐consuming activity. Many techniques are available, ranging from using environmental tracers, to LiDAR. Erosion measurements using erosion pins are assessed and compared with regional scale erosion data, hillslope data obtained using ¹³⁷Cs and erosion modelling results. The pins produced erosion rates which are within the range determined using ¹³⁷Cs and model data but above that of regional denudation rates. Our findings demonstrate that inexpensive erosion pins can provide reliable data on hillslope erosion. © 2015 Commonwealth of Australia. Hydrological Processes © 2015 John Wiley & Sons Ltd.     

Hillslope erosion in a grassland environment: Calibration and evaluation of the SIBERIA landscape evolution model

Field measurement and modelling of soil erosion provides insights into landscape systems as well as the potential for enhanced landscape management. There are a number of field and numerical methods by which soil erosion and deposition can be quantified. Here we examine the capability of the SIBERIA landscape evolution model to quantify short-term erosion and deposition on a well-managed cattle grazing landscape on the east coast of Australia. The model is calibrated by two methods (1) a geomorphological approach using a site digital elevation model (DEM) and soil data and (2) a laboratory-scale flume. The two calibration processes resulted in similar model input parameters and estimated erosion rates of 3.1 t ha⁻¹ year⁻¹ and 4.4 t ha⁻¹ year⁻¹, respectively. These were found to closely match erosion rates estimated using the environmental tracer ¹³⁷Cs (2.7–4.8 t ha⁻¹ year⁻¹). However, erosion and deposition estimated at individual points along the hillslope was not well correlated with ¹³⁷Cs at the same position due to the temporal averaging of the model and microtopography. Sensitivity analysis showed the model was more sensitive to parameterisation than sub-DEM-scale topography. This places confidence in the model's ability to estimate erosion and deposition across an entire hillslope and catchment on decadal time scales. We also highlight the robustness and flexibility of the calibration methods.     

Distance,time and scale in soil erosion processes

This commentary brings together, as a virtual Special Issue, a number of recent papers in Earth Surface Processes and Landforms that are all related to issues of scale in soil erosion. Empirical concepts that were developed in the 1940s are now in need of re‐thinking, and papers are increasingly exploring, through modelling and measurement, appropriate ways to recognize the mechanisms that connect processes across time and space scales. Issues include a more nuanced approach to selective transportation, responses to variability in surface and sub‐surface conditions and the need to analyse measurements in ways that can be transferred between sites and storms. Copyright © 2010 John Wiley & Sons, Ltd.     

Integrating in situ measurements of an index of connectivity to assess soil erosion processes in vineyards

ABSTRACT

In this research, we survey soil erosion processes using an index of connectivity and a non-invasive and long-term assessment in situ technique: the improved stock unearthing method (ISUM), for which the vineyard of Castilla La Mancha under tillage management was selected. Our results show, that in 10 years, the total average soil surface level decreased by – 1.6 cm and the total soil mobilization was up to 17.7 Mg ha^?1 year^?1. The surrounding lands of the survey plot showed meagre connectivity; however, smaller linear features with higher values were found as possible locations for potential rill generation. The survey plot is traversed by one of these linear features. In this inter-row survey, we found relatively low connectivity values. As a conclusion, we confirm that both methods can be useful to assess soil erosion processes in vineyards and detect areas that could increase the desertification as a consequence of non-sustainable soil erosion rates.     

Tracing the temporal evolution of soil redistribution rates in an agricultural landscape using 239+240Pu and 10Be

Two principal groups of processes shape mass fluxes from and into a soil: vertical profile development and lateral soil redistribution. Periods having predominantly progressive soil forming processes (soil profile development) alternate with periods having predominantly regressive processes (erosion). As a result, short-term soil redistribution – years to decades – can differ substantially from long-term soil redistribution; i.e. centuries to millennia. However, the quantification of these processes is difficult and consequently their rates are poorly understood. To assess the competing roles of erosion and deposition we determined short- and long-term soil redistribution rates in a formerly glaciated area of the Uckermark, northeast Germany. We compared short-term erosion or accumulation rates using plutonium-239 and -240 (^239+240Pu) and long-term rates using both in situ and meteoric cosmogenic beryllium-10 (¹⁰Be). Three characteristic process domains have been analysed in detail: a flat landscape position having no erosion/deposition, an erosion-dominated mid-slope, and a deposition-dominated lower-slope site. We show that the short-term mass erosion and accumulation rates are about one order of magnitude higher than long-term redistribution rates. Both, in situ and meteoric ¹⁰Be provide comparable results. Depth functions, and therefore not only an average value of the topsoil, give the most meaningful rates. The long-term soil redistribution rates were in the range of −2.1 t ha^-1 yr^-1 (erosion) and +0.26 t ha^-1 yr^-1 (accumulation) whereas the short-term erosion rates indicated strong erosion of up to 25 t ha^-1 yr^-1 and accumulation of 7.6 t ha^-1 yr^-1. Our multi-isotope method identifies periods of erosion and deposition, confirming the ‘time-split approach’ of distinct different phases (progressive/regressive) in soil evolution. With such an approach, temporally-changing processes can be disentangled, which allows the identification of both the dimensions of and the increase in soil erosion due to human influence. © 2019 John Wiley & Sons, Ltd.     

Quantifying hillslope erosion rates and processes for a coastal California landscape over varying timescales

The Earth's surface erodes by processes that occur over different spatial and temporal scales. Both continuous, low‐magnitude processes as well as infrequent, high‐magnitude events drive erosion of hilly soil‐mantled landscapes. To determine the potential variability of erosion rates we applied three independent, field‐based methods to a well‐studied catchment in the Marin Headlands of northern California. We present short‐term, basin‐wide erosion rates determined by measuring pond sediment volume (40 years) and measured activities of the fallout nuclides ¹³⁷Cs and ²¹⁰Pb (40–50 years) for comparison with long‐term (>10 ka) rates previously determined from in situ‐produced cosmogenic ¹⁰Be and ²⁶Al analyses. In addition to determining basin‐averaged rates, ¹³⁷Cs and ²¹⁰Pb enable us to calculate point‐specific erosion rates and use these rates to infer dominant erosion processes across the landscape. When examined in the context of established geomorphic transport laws, the correlations between point rates of soil loss from ¹³⁷Cs and ²¹⁰Pb inventories and landscape morphometry (i.e. topographic curvature and upslope drainage area) demonstrate that slope‐driven processes dominate on convex areas while overland flow processes dominate in concave hollows and channels. We show a good agreement in erosion rates determined by three independent methods: equivalent denudation rates of 143 ± 41 m Ma^?1 from pond sediment volume, 136 ± 36 m Ma^?1 from the combination of ¹³⁷Cs and ²¹⁰Pb, and 102 ± 25 m Ma^?1 from ¹⁰Be and ²⁶Al. Such agreement suggests that erosion of this landscape is not dominated by extreme events; rather, the rates and processes observed today are indicative of those operating for at least the past 10 000 years. Copyright © 2006 John Wiley & Sons, Ltd.     

Rill development and soil erosion: a laboratory study of slope and rainfall intensity

A total of 15 rainfall simulation experiments were conducted in a 1 m by 2 m box varying slope (10, 20, 30%) and rainfall intensity (60, 90, 120 mm h^?1). The experiments were performed to study how rill networks initiate and evolve over time under controlled conditions with regard to the treatment variables considered, and to allow for input in a computer simulation model. Runoff and sediment yield samples were collected. Digital elevation models were calculated by means of photogrammetry for several time steps of most experiments. The soil used in the experiments was a basal till derived Cambisol typical for the Swiss Plateau. While significant differences were found for sediment yield, runoff did not vary significantly with treatment combinations. Increasing rainfall intensity had a larger effect on sediment yield than increasing slope. Rill density and energy expenditure decreased with time, suggesting that energy expenditure was a useful parameter to describe the emergence of rill network at the laboratory scale. Copyright © 2010 John Wiley & Sons, Ltd.     

The increase of rainfall erosivity and initial soil erosion processes due to rainfall acidification

The drastic growth of population in highly industrialized urban areas, as well as fossil fuel use, is increasing levels of airborne pollutants and enhancing acid rain. In rapidly developing countries such as Iran, the occurrence of acid rain has also increased. Acid rain is a driving factor of erosion due to the destructive effects on biota and aggregate stability; however, little is known about its impact on specific rates of erosion at the pedon scale. Thus, the present study aimed to investigate the effect of acid rain at pH levels of 5.25, 4.25, and 3.75 for rainfall intensities of 40, 60, and 80 mm h^?1 on initial soil erosion processes under dry and saturated soil conditions using rainfall simulations. The results were compared using a two‐way ANOVA and Duncan tests and showed that initial soil erosion rates with acidic rain and non‐acidic rain under dry soil conditions were significantly different. The highest levels of soil particle loss due to splash effects in all rainfall intensities were observed with the most acidic rain (pH = 3.75), reaching maximum values of 16 g m^?2 min^?1. The lowest levels of particle losses were observed in the control plot where non‐acidic rain was used, with values ranging from 3.8 to 8.1 g m^?2 min^?1. Similarly, under saturated soil conditions, the lowest level of soil particle loss was observed in the control plot, and the highest peaks of soil loss were observed for the most acidic rains (pH = 3.75 and pH = 4.25), reaching maximum average values of 40 g m^?2 min^?1. However, for saturated soils with acidic water but with non‐acidic rain, the highest soil particle loss was observed for the control plot for all the rainfall intensities. In conclusion, acidic rain has a negative impact on soils, which can be more intense with a concomitant increase in rainfall intensity. Rapid solutions, therefore, need to be found to reduce the emission of pollutants into the air, otherwise, rainfall erosivity may drastically increase.     

Interrill soil erosion processes and their interaction on low slopes

Soil erosion by water is mostly the result of rainfall‐driven and runoff‐driven processes taking place simultaneously during a storm event. However, the effect of interaction between these two erosion processes has received limited attention. Most laboratory experiments indicate that the rate of erosion in a rain‐impacted flow is greater than for un‐impacted flows of similar depth and velocity; however, negative interaction between the two processes has also been reported. There is no provision for any such interaction in any of the current erosion models. This paper reports on the results of a number of exact experiments on three soil types carried out in the flume of Griffith University's large rainfall simulator to study interaction between rain and runoff processes. The results show that interaction is generally positive under approximately steady state condition and there is very limited sign of negative interaction reported by others. Results provide strong evidence that raindrops continuously peel fine sediment from larger stable aggregates. This mechanism could be the reason for positive interaction during simultaneous rainfall and flow driven erosion in well aggregated soils as a result of increased fine particles in the eroded sediment. Strong positive interaction between rain and runoff erosion also occurs for medium to large aggregates. This strongly suggests that mechanisms that are not well understood are operational. It is quite possible that particle movement can be stimulated by rolling or creeping in a size‐selective manner. Indeed, such additional mechanisms may well be largely responsible for the positive interaction observed between rain and surface flow. Copyright © 2006 John Wiley & Sons, Ltd.     

A conceptual model for determining soil erosion by water

Current estimates of rates of soil erosion by water derived from plots are incompatible with estimates of long‐term lowering of large drainage basins. Traditional arguments to reconcile these two disparate rates are ?awed. The ?ux of sediment leaving a speci?ed area cannot be converted to a yield simply by dividing by the area, because there is no simple relationship between ?ux and area. Here, we develop an approach to the determination of erosion rates that is based upon the entrainment rates and travel distances of individual particles. The limited available empirical data is consistent with the predictions of this approach. Parameterization of the equations to take account of such factors as gradient and sediment supply is required to proceed from the conceptual framework to quantitative measurements of erosion. However, our conceptual model solves the apparent paradox of the sediment delivery ratio, resolves recent discussion about the validity of erosion rates made using USLE erosion plots, and potentially can reconcile erosion rates with known lifespans of continents. Our results imply that previous estimates of soil erosion are fallacious. Copyright © 2004 John Wiley & Sons, Ltd.     

Impact of vegetation on erosion: Insights from the calibration and test of a landscape evolution model in alpine badland catchments

While it is well recognized that vegetation can affect erosion, sediment yield and, over longer timescales, landform evolution, the nature of this interaction and how it should be modeled is not obvious and may depend on the study site. In order to develop quantitative insight into the magnitude and nature of the influence of vegetation on catchment erosion, we build a landscape evolution model to simulate erosion in badlands, then calibrate and evaluate it against sediment yield data for two catchments with contrasting vegetation cover. The model couples hillslope gravitational transport and stream alluvium transport. Results indicate that hillslope transport processes depend strongly on the vegetation cover, whereas stream transport processes do not seem to be affected by the presence of vegetation. The model performance in prediction is found to be higher for the denuded catchment than for the reforested one. Moreover, we find that vegetation acts on erosion mostly by reducing soil erodibility rather than by reducing surface runoff. Finally, the methodology we propose can be a useful tool to evaluate the efficiency of previous revegetation operations and to provide guidance for future restoration work. © 2019 John Wiley & Sons, Ltd.     

Sediment transport capacity and erosion processes: model concepts and reality

Sediment transport capacity, T_c, defined as the maximum amount of sediment that a flow can carry, is the basic concept in determining detachment and deposition processes in current process-based erosion models. Although defined conceptually and used extensively in modelling erosion, T_c was rarely measured. Recently, a series of laboratory studies designed to quantify effects of surface hydrologic conditions on erosion processes produced data sets feasible to evaluate the concept of T_c. A dual-box system, consisting of 1·8 m long sediment feeder box and a 5 m long test box, was used. Depending on the relative magnitudes of sediment delivery from feeder and test boxes, five scenarios are proposed ranging from deposition-dominated to transport-dominated sediment regimes. Results showed that at 5 per cent slope under seepage or 10 per cent slope under drainage conditions, the runoff from the feeder box caused in the additional sediment transport in the test box, indicating a transport-dominated sediment regime. At 5 per cent slope under drainage conditions, deposition occurred at low rainfall intensities. Increases in slope steepness, rainfall intensity and soil erodibility shifted the dominant erosion process from deposition to transport. Erosion process concepts from the Meyer–Wishmeier, Foster–Meyer and Rose models were compared with the experimental data, and the Rose model was found to best describe processes occurring during rain. A process-based erosion model needs to have components that can represent surface conditions and physical processes and their dynamic interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

Model‐based evaluation of impact of soil redistribution on soil organic carbon stocks in a temperate hedgerow landscape

Agroforestry systems are promoted for providing a number of ecosystem services and environmental benefits, including soil protection and carbon sequestration. This study proposes a modelling approach to quantify the impact of soil redistribution on soil organic carbon (SOC) storage in a temperate hedgerow landscape. Evolution of SOC stocks at the landscape scale was examined by simulating vertical and horizontal SOC transfers in the 0–105 cm soil layer due to soil redistribution by tillage and water processes. A spatially explicit SOC dynamics model (adapted from RothC‐26.3) was used, coupled with a soil‐redistribution model (LandSoil). SOC dynamics were simulated over 90 years in an agricultural hedgerow landscape dedicated to dairy farming, with a mix of cropping and grasslands. Climate and land use were simulated considering business‐as‐usual scenarios derived from existing information on the study area. A net decrease in SOC stocks was predicted at the end of the simulation period. Soil redistribution induced a net SOC loss equivalent to 2 kg C ha^?1 yr^?1 because of soil exportation out of the study site and an increase in SOC mineralization. Hedgerows and woods were the only land use in which soil redistribution induced net SOC storage. Soil tillage was the main process that induced soil redistribution within cultivated fields. Soil exportation out of the study area was due to erosion by water, but remained low because of the protective role of the hedgerow network. These soil transfers redistributed SOC stocks in the landscape, mostly within cultivated fields. Copyright © 2016 John Wiley & Sons, Ltd.     

The extent of soil erosion in upland England and Wales

Following recommendations by the 19th Royal Commission on Environmental Pollution, the area, causes and rates of upland soil erosion in England and Wales were investigated between 1997 and 1999. This paper describes the methods and results of the field survey of 1999 in which the extent of eroded ground was determined. 2. The area of degraded soil and the volume of eroded material were both determined from the dimensions of individual erosion features at 399 field sites located on an orthogonal grid across the uplands. Using measurements of individual erosion features, degraded soil extent in upland England and Wales was estimated at almost 25 000 ha, 2·46 per cent of the total upland area surveyed. Half this eroded area was revegetated and no longer subject to continued accelerated soil loss in 1999. The total volume of eroded material was estimated at 0·284 km³. Although deposition of eroded material occurred within 20 per cent of eroded field sites, the total volume of redeposited material was less than 1 per cent of the total volume of eroded soil. 3. Erosion was more extensive on peat soils than on dry, wet mineral or wet peaty mineral soils. In addition, the higher incidence of erosion at high altitudes and on low slopes reinforced the relationship between erosion and areas of peat formation. Copyright © 2002 John Wiley & Sons, Ltd.