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     

Testing of the SIBERIA landscape evolution model using the Tin Camp Creek,Northern Territory,Australia, field catchment

The SIBERIA landscape evolution model was used to simulate the geomorphic development of the Tin Camp Creek natural catchment over geological time. Measured hydrology, erosion and geomorphic data were used to calibrate the SIBERIA model, which was then used to make independent predictions of the landform geomorphology of the study site. The catchment, located in the Northern Territory, Australia is relatively untouched by Europeans so the hydrological and erosion processes that shaped the area can be assumed to be the same today as they have been in the past, subject to the caveats regarding long‐term climate fluctuation. A qualitative, or visual comparison between the natural and simulated catchments indicates that SIBERIA can match hillslope length and hillslope profile of the natural catchments. A comparison of geomorphic and hydrological statistics such as the hypsometric curve, width function, cumulative area distribution and area–slope relationship indicates that SIBERIA can model the geomorphology of the selected Tin Camp Creek catchments. Copyright 2002 © Environmental Research Institute of the Supervising Scientist, Commonwealth of Australia.     

A catchment scale evaluation of the SIBERIA and CAESAR landscape evolution models

Landscape evolution models provide a way to determine erosion rates and landscape stability over times scales from tens to thousands of years. The SIBERIA and CAESAR landscape evolution models both have the capability to simulate catchment–wide erosion and deposition over these time scales. They are both cellular, operate over a digital elevation model of the landscape, and represent fluvial and slope processes. However, they were initially developed to solve research questions at different time and space scales and subsequently the perspective, detail and process representation vary considerably between the models. Notably, CAESAR simulates individual events with a greater emphasis on fluvial processes whereas SIBERIA averages erosion rates across annual time scales. This paper describes how both models are applied to Tin Camp Creek, Northern Territory, Australia, where soil erosion rates have been closely monitored over the last 10 years. Results simulating 10 000 years of erosion are similar, yet also pick up subtle differences that indicate the relative strengths and weaknesses of the two models. The results from both the SIBERIA and CAESAR models compare well with independent field data determined for the site over different time scales. Representative hillslope cross‐sections are very similar between the models. Geomorphologically there was little difference between the modelled catchments after 1000 years but significant differences were revealed at longer simulation times. Importantly, both models show that they are sensitive to input parameters and that hydrology and erosion parameter derivation has long‐term implications for sediment transport prediction. Therefore selection of input parameters is critical. This study also provides a good example of how different models may be better suited to different applications or research questions. Copyright © 2010 John Wiley & Sons, Ltd and Commonwealth of Australia     

The use of a landscape simulator in the validation of the SIBERIA landscape evolution model: transient landforms

SIBERIA is a physically based model for the geomorphic evolution of landforms. It is essential that the SIBERIA model be tested or validated against controlled landform development. Previous studies have demonstrated that SIBERIA is able to simulate declining equilibrium landforms and in this paper we examine SIBERIA's ability to simulate landforms as they evolve to their declining equilibrium form. These landscapes are termed transient landforms. Landscapes generated by SIBERIA were compared to those produced by a physical model (experimental model landforms) at stages of evolution. Comparison of the experimental landscapes with the simulated landscapes using total mass, hypsometric curve, width function, cumulative area distribution and area–slope demonstrate that SIBERIA can simulate the experimental model landscape during development (i.e. transient landscapes). Copyright © 2002 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.     

Modelling landscape evolution

Geomorphology is currently in a period of resurgence as we seek to explain the diversity, origins and dynamics of terrain on the Earth and other planets in an era of increased environmental awareness. Yet there is a great deal we still do not know about the physics and chemistry of the processes that weaken rock and transport mass across a planet's surface. Discovering and refining the relevant geomorphic transport functions requires a combination of careful field measurements, lab experiments, and use of longer‐term natural experiments to test current theory and develop new understandings. Landscape evolution models have an important role to play in sharpening our thinking, guiding us toward the right observables, and mapping out the logical consequences of transport laws, both alone and in combination with other salient processes. Improved quantitative characterization of terrain and process, and an ever‐improving theory that describes the continual modification of topography by the many and varied processes that shape it, together with improved observation and qualitative and quantitative modelling of geology, vegetation and erosion processes, will provide insights into the mechanisms that control catchment form and function. This paper reviews landscape theory – in the form of numerical models of drainage basin evolution and the current knowledge gaps and future computing challenges that exist. Copyright © 2010 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.     

Exploring effects of rainfall intensity and duration on soil erosion at the catchment scale using openLISEM: Prado catchment,SE Spain

In semi‐arid areas, high‐intensity rainfall events are often held responsible for the main part of soil erosion. Long‐term landscape evolution models usually use average annual rainfall as input, making the evaluation of single events impossible. Event‐based soil erosion models are better suited for this purpose but cannot be used to simulate longer timescales and are usually applied to plots or small catchments. In this study, the openLISEM event‐based erosion model was applied to the medium‐sized (～50 km²) Prado catchment in SE Spain. Our aim was to (i) test the model's performance for medium‐sized catchments, (ii) test the ability to simulate four selected typical Mediterranean rainfall events of different magnitude and (iii) explore the relative contribution of these different storms to soil erosion using scenarios of future climate variability. Results show that because of large differences in the hydrologic response between storms of different magnitudes, each event needed to be calibrated separately. The relation between rainfall event characteristics and the calibration factors might help in determining optimal calibration values if event characteristics are known. Calibration of the model features some drawbacks for large catchments due to spatial variability in K_sat values. Scenario calculations show that although ～50% of soil erosion occurs as a result of high frequency, low‐intensity rainfall events, large‐magnitude, low‐frequency events potentially contribute significantly to total soil erosion. The results illustrate the need to incorporate temporal variability in rainfall magnitude–frequency distributions in landscape evolution models. Copyright © 2011 John Wiley & Sons, Ltd.     

An evaluation of an enhanced soil erosion and landscape evolution model: a case study assessment of the former Nabarlek uranium mine,Northern Territory,Australia

The assessment of post‐mining landscapes as case studies is an important part of the evaluation of current rehabilitation practices. A necessary part of this assessment is to predict the surface stability of the landform using erosion and landform evolution modelling techniques. In the short term, erosion on a rehabilitated mine site can lead to increased sediment loads and transport of other mine related contaminants in downstream waterways. It is well recognized that in many mine areas the erodibility of surface materials can, and does, vary. This is a particularly significant issue on mine sites, where the surface conditions may range from areas of undisturbed natural surface materials, waste rock dumps constructed with materials exhumed from the sub‐surface, and other areas that have a mix of waste rock and soil to enhance the growth of vegetation. A further significant issue is that when the subsurface materials are exposed to surface conditions they can weather rapidly, changing their erodibility. This paper uses a new version of the SIBERIA landscape evolution and soil erosion model to evaluate the former Nabarlek uranium mine site in the Northern Territory, Australia. This new version of SIBERIA uses spatially variable erosion and hydrology parameters across the study domain to represent different erodibilities of surface materials, thus allowing better representation of catchment heterogeneity. The results demonstrate that the model predicts erosion rates similar to that of other modelled results and independent field data, providing confidence in the model and its parameterization. The tailings, deposited in the mined out pit and capped with waste rock, appear to be safely encapsulated for the modelled period. Copyright © 2008 John Wiley & Sons, Ltd and Commonwealth of Australia (Department of the Environment and Water Resources Supervising Scientist).     

Disaggregating soil erosion processes within an evolving experimental landscape

Soil‐mantled landscapes subjected to rainfall, runoff events, and downstream base level adjustments will erode and evolve in time and space. Yet the precise mechanisms for soil erosion also will vary, and such variations may not be adequately captured by soil erosion prediction technology. This study sought to monitor erosion processes within an experimental landscape filled with packed homogenous soil, which was exogenically forced by rainfall and base level adjustments, and to define the temporal and spatial variation of the erosion regimes. Close‐range photogrammetry and terrain analysis were employed as the primary methods to discriminate these erosion regimes. Results show that (1) four distinct erosion regimes can be identified (raindrop impact, sheet flow, rill, and gully), and these regimes conformed to an expected trajectory of landscape evolution; (2) as the landscape evolved, the erosion regimes varied in areal coverage and in relative contribution to total sediment efflux measured at the outlet of the catchment; and (3) the sheet flow and rill erosion regimes dominated the contributions to total soil loss. Disaggregating the soil erosion processes greatly facilitated identifying and mapping each regime in time and space. Such information has important implications for improving soil erosion prediction technology, for assessing landscape degradation by soil erosion, for mapping regions vulnerable to future erosion, and for mitigating soil losses and managing soil resources. Copyright © 2017 John Wiley & Sons, Ltd.     

Predicting the impact of linear landscape elements on surface runoff,soil erosion,and sedimentation in the Wahnbach catchment,Germany

Model predictions concerning the endangerment of on‐site and off‐site damages due to runoff, soil erosion and sedimentation under alternative design and operation policies are of particular importance in recent catchment planning and management. By using the raster‐based model approach, linear landscape elements, such as streets and roads, and their impacts on flow paths are often neglected. Therefore, the aim of this study was to analyse the effects of linear landscape elements on patterns of soil erosion, sediment transport and sedimentation. To accomplish this, roads are considered while determining flow paths. Simulations in the well‐investigated catchment of the Wahnbach River (54 km²) in a low mountain range in Germany were carried out using a combination of different models for hydrology and soil erosion. Although the study focuses on the catchment scale of the Wahnbach River, detailed investigations concerning the sub‐catchment scale (21 ha) were also conducted. The simulation results show that these spatial structures mainly affect the pattern of soil erosion and sedimentation. On the sub‐catchment scale, improved identification of active zones for sediment dynamic becomes possible. On the catchment scale, the predicted runoff is about 20% higher, and sediment outputs were four times larger than predicted when roads were considered. Soil erosion increases by 37% whereas sedimentation is reduced by 29%. The model improvement could not be evaluated on the catchment scale because of the high variability and heterogeneity of land use and soils, but road impacts could be explained by simulations on the sub‐catchment scale. It can be concluded that runoff concentration due to rerouted flow paths leads to lower non‐concentrated and higher concentric‐linear surface runoff. Thus, infiltration losses decline and surface runoff and soil erosion increase because sedimentation is reduced. Further, runoff concentration can cause soil erosion hot spots. In the model concept used in this study, buffering of runoff and sediments on the upslope side of roads and in local depressions adjacent to roads cannot be simulated. Flow paths will only be rerouted because of road impacts, but the temporal ponding of water is not simulated. Therefore, the drastic increase of predicted sediment output due to road impact does not seem to be reliable. However, results indicate that the consideration of roads when determining flow paths enabled more detailed simulations of surface runoff, soil erosion and sedimentation. Thus, progress in model‐based decision‐making support for river catchment planning and management can be achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

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

A qualitative and quantitative evaluation of experimental model catchment evolution

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.     

Linking uplift,erosion, and sedimentation using landscape evolution models: Madagascar since the Late Cretaceous

We present a study to estimate the large-scale landscape history of a continental margin, by establishing a source-to-sink volume balance between the eroding onshore areas and the offshore basins. Assuming erosion as the primary process for sediment production, we strive to constrain a numerical model of landscape evolution that balances the volumes of eroded materials from the continent and that deposited in the corresponding basins, with a ratio imposed for loss of erosion products. We use this approach to investigate the landscape history of Madagascar since the Late Cretaceous. The uplift history prescribed in the model is inferred from elevations of planation surfaces formed at various ages. By fitting the volumes of terrigenous sediments in the Morondava Basin along the west coast and the current elevation of the island, the landscape evolution model is optimized by constraining the erosion law parameters and ratios of sediment loss. The results include a best-fit landscape evolution model, which features two major periods of uplift and erosion during the Late Cretaceous and the middle to late Cenozoic. The model supports suggestions from previous studies that most of the high topography of the island was constructed since the middle to late Miocene, and on the central plateau the erosion has not reached an equilibrium with the high uplift rates in the late Cenozoic. Our models also indicate that over the geological time scale a significant portion of materials eroded from Madagascar was not archived in the offshore basin, possibly consumed by chemical weathering, the intensity of which might have varied with climate.     

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.     

Modeling soil erosion using a spatially distributed model in a karst catchment of northwest Guangxi,China

Reliable assessment of the spatial distribution of soil erosion is important for making land management decisions, but it has not been thoroughly evaluated in karst geo‐environments. The objective of this study was to modify a physically based, spatially distributed erosion model, the revised Morgan, Morgan and Finney (RMMF) model, to estimate the superficial (as opposed to subsurface creep) soil erosion rates and their spatial patterns in a 1022 ha karst catchment in northwest Guangxi, China. Model parameters were calculated using local data in a raster geographic information system (GIS) framework. The cumulative runoff on each grid cell, as an input to the RMMF model for erosion computations, was computed using a combined flow algorithm that allowed for flow into multiple cells with a transfer grid considering infiltration and runoff seepage to the subsurface. The predicted spatial distributions of soil erosion rates were analyzed relative to land uses and slope zones. Results showed that the simulated effective runoff and annual soil erosion rates of hillslopes agreed well with the field observations and previous quantified redistribution rates with caesium‐137 (¹³⁷Cs). The estimated average effective runoff and annual erosion rate on hillslopes of the study catchment were 18 mm and 0.27 Mg ha^?1 yr^?1 during 2006–2007. Human disturbances played an important role in accelerating soil erosion rates with the average values ranged from 0.1 to 3.02 Mg ha^?1 yr^?1 for different land uses. The study indicated that the modified model was effective to predict superficial soil erosion rates in karst regions and the spatial distribution results could provide useful information for developing local soil and water conservation plans. Copyright © 2014 John Wiley & Sons, Ltd.     

Initial hydro‐geomorphic development and rill network evolution in an artificial catchment

The formation of erosion rills and gullies is a critical step in land surface development, but possibilities to study initial unaffected surface development under natural conditions and with well‐defined initial and boundary conditions are rare. The objective of this study was to characterize rill network development from ’point zero’ in the artificially‐created catchment ‘Hühnerwasser’. To ensure unaffected development, the study was largely restricted to the analysis of remotely‐sensed data. We analyzed a series of photogrammetry‐based digital elevation models (DEMs) for 10 points in time, over a period of five years and beginning with the initial state. The evolving erosion rill network was quantitatively described based on mapping from aerial photographs. DEMs and rill network maps were combined to specifically analyze the development of morphometry for different parts of the network and to characterize energy dissipation and connectivity. The restriction to remote‐sensing data did not allow for analyzing specific processes governing rill network development, nevertheless, two major development phases could be characterized. We observed a phase of growth of the rill network along with variations in drainage patterns during the first two years of development and a subsequent phase of reduction of its area along with comparably stable patterns. Region‐specific analysis of morphometry indicates that, besides effects of changing sediment characteristics and vegetation cover development, locally evolving hydro‐geomorphic feedback cycles influenced this development. Results show an increasing similarity of overall statistical characteristics (e.g. drainage density) for two parts of the catchment, but a persistent influence of initial conditions on specific rill geometry. The observed development towards higher orderliness and increased connectivity is consistent with experiments and concepts on drainage network evolution across scales; however, we did not observe major influences of rill piracy and cross grading or a reduction of energy dissipation with network development. Copyright © 2013 John Wiley & Sons, Ltd.     

The effect of storm movement on infiltration,runoff and soil erosion in a semi-arid catchment

Rainfall characteristics are key factors influencing infiltration and runoff generation in catchment hydrology, particularly for arid and semiarid catchments. Although the effect of storm movement on rainfall-runoff processes has been evaluated and emphasized since the 1960s, the effect on the infiltration process has barely been considered. In this study, a physically based distributed hydrological model (InHM) was applied to a typical semi-arid catchment (Shejiagou, 4.26 km²) located in the Loess Plateau, China, to investigate the effect of storm movement on infiltration, runoff and soil erosion at the catchment scale. Simulations of 84 scenarios of storm movement were conducted, including storms moving across the catchment in both the upstream and downstream directions along the main channel, while in each direction considering four storm moving speeds, three rainfall depths and two storm ranges. The simulation results showed that, on both the hillslopes facing downstream (facing south) and in the main channel, the duration of the overland flow process under the upstream-moving storms was longer than that under the downstream-moving storms. Thus, the duration and volume of infiltration under upstream-moving storms were larger in these areas. For the Shejiagou catchment, as there are more hillslopes facing downstream, more infiltration occurred under the upstream-moving storms than the downstream-moving storms. Therefore, downstream-moving storms generated up to 69% larger total runoff and up to 351% more soil loss in the catchment than upstream-moving storms. The difference in infiltration between the storms moving upstream and downstream decreased as the storm moving speed increased. The relative difference in total runoff and sediment yield between the storms moving upstream and downstream decreased with increasing rainfall depth and storm speed. The results of this study revealed that the infiltration differences under moving storms largely influenced the total runoff and sediment yield at the catchment scale, which is of importance in runoff prediction and flood management. The infiltration differences may be a potential factor leading to different groundwater, vegetation cover and ecology conditions for the different sides of the hillslopes.