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     

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.     

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.     

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.     

Postglacial to Holocene landscape evolution and process rates in steep alpine catchments

Climate change and high magnitude mass wasting events pose adverse societal effects and hazards, especially in alpine regions. Quantification of such geomorphic processes and their rates is therefore critical but is often hampered by the lack of appropriate techniques and the various spatiotemporal scales involved in these studies. Here we exploit both in situ cosmogenic beryllium-10 (¹⁰Be) and carbon-14 (¹⁴C) nuclide concentrations for deducing exposure ages and tracing of sediment through small alpine debris flow catchments in central Switzerland. The sediment cascade and modern processes we track from the source areas, through debris flow torrents to their final export out into sink regions with cosmogenic nuclides over an unprecedented five-year time series with seasonal resolution. Data from a seismic survey and a 90 m core revealed a glacially overdeepened basin, filled with glacial and paraglacial sediments. Surface exposure dating of fan boulders and radiocarbon ages constrain the valley fill from the last deglaciation until the Holocene and show that most of the fan existed in early Holocene times already. Current fan processes are controlled by episodic debris flow activity, snow (firn) and rock avalanches. Field investigations, digital elevation models (DEMs) of difference and geomorphic analysis agree with sediment fingerprinting with cosmogenic nuclides, highlighting that the bulk of material exported today at the outlet of the subcatchments derives from the lower fans. Cosmogenic nuclide concentrations steadily decrease from headwater sources to distal fan channels due to the incorporation of material with lower nuclide concentrations. Further downstream the admixture of sediment from catchments with less frequent debris flow activity can dilute the cosmogenic nuclide signals from debris flow dominated catchments but may also reach thresholds where buffering is limited. Consequently, careful assessment of boundary conditions and driving forces is required when apparent denudation rates derived from cosmogenic nuclide analysis are upscaled to larger regions. © 2018 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.     

A 13‐year record of erosion on badland sites in the Karoo,South Africa

Land degradation in South Africa has been of concern for more than 100 years with both climate change and inappropriate land management (overgrazing) being proposed as primary drivers. However, there are few quantitative studies of degradation and, in particular, few of erosion by water. Badlands, taken here to be the landform which results from extreme erosion, have been notably neglected. We report on 13 consecutive years of erosion pin measurements of badland erosion on 10 study sites in the Sneeuberg uplands of the eastern Karoo in South Africa. The study sites are on Holocene colluvium which mantles footslopes. They have been subject to overgrazing for at least 100 years, c. 1850–1950. Currently they are lightly grazed by sheep. The area receives about 500 mm rainfall per year. The sites are remote, with only informal, farmer‐operated, daily raingauges nearby. The nearest sub‐daily raingauge is c. 55 km distant. Also we report on an analysis of the erosion pin data which focuses on establishing the origins and context of the badlands, including the relationship between study sites and adjacent valley‐bottom gully systems; compare erosion rates on our study sites with rates determined by erosion pins on other badland sites; and discuss the implications of these erosion rates for landscape development and off‐site impacts. Net erosion rates on the study sites are relatively high compared with global badland rates and range from 3.1 to 8.5 mm yr^‐1 which may be extrapolated to 53 to 145 t ha yr^‐1 (using a measured bulk density of 1.7 g cm^‐3). However, comparisons with badland sites elsewhere are difficult because of different measuring methodologies, lithologies, climate and dominant processes. Erosion rates on the study sites are strongly influenced by rainfall amounts and, in particular, by daily rainfall events which exceed ~10 mm: this is the threshold intensity at which runoff has been observed to commence on badlands. Of significance, but of lesser influence, is weathering, mainly by wetting and drying: this prepares bare surfaces for erosion. However, questions remain regarding the role of site characteristics, and of processes at each site, in determining between‐site differences in erosion rate. Crude extrapolation of current rates of erosion, in conjunction with depths of incision into the badlands, suggests that badland development started around 200 years ago, probably as a response to the introduction of European‐style stock farming which resulted in overgrazing. We assume, but cannot quantify, the additional influence of periods of drought and burning in the erosional history of the area. Intermittent connection of these badlands to valley‐bottom gullies and therefore to small farm dams and ultimately to large water storage reservoirs increases their impact on local water resources. Copyright © 2015 John Wiley & Sons, Ltd.     

Landscapes on the edge: examining the role of climatic interactions in shaping coastal watersheds using a coastal–terrestrial landscape evolution model

Incised coastal gullies (ICGs) are dynamic features found at the terrestrial‐coastal interface. Their geomorphic evolution is driven by the interactions between processes of fluvial knickpoint migration and coastal cliff erosion. Under scenarios of future climate change the frequency and magnitude of the climatological drivers of both terrestrial (fluvial and hillslope) and coastal (cliff erosion) processes are likely to change, with an adjunct impact on these types of coastal features. Here we explore the response of an incised coastal gully to changes in both terrestrial and coastal climate in order to elucidate the key process interactions which drive ICG evolution. We modify an extant landscape evolution model, CHILD, to incorporate processes of soft‐cliff erosion. This modified version, termed the Coastal‐Terrestrial‐CHILD (CT‐CHILD) model, is then employed to explore the interactions between changing terrestrial and coastal driving forces on the future evolution of an ICG found on the south‐west Isle of Wight, UK. It was found that the magnitude and frequency of storm events will play a key role in determining the future trajectory of ICGs, highlighting a need to understand the role of event sequencing in future projections of landscape evolution. Furthermore, synergistic (positive) and antagonistic (negative) interactions were identified between coastal and terrestrial parameters, such as wave height intensity and precipitation duration, which act to modulate the impact of changes in any one parameter. Of note was the role played by wave height intensity in driving coastal erosion, which was found to play a more important role than sea‐level rise in determining rates of coastal erosion. This highlights the need for a greater focus on wave height in studies of soft‐cliff erosion. Copyright © 2014 John Wiley & Sons, Ltd.     

Impact of hillslope vegetation restoration on gully erosion induced sediment yield

As a result of serious soil erosion on the Loess Pla-teau of China, about 1.6 billion tons of silt discharge into the downstream and 0.4 billion tons deposit on the riverbed every year, causing serious threat to the life and property of the local people on both banks of the lower Yellow River[1]. Since the 1950s, the Chinese government has initiated the work on soil and water conservation and environmental management on the Loess Plateau and formulated the guiding principle of hillslope and …     

Process–form relationships in Southern Italian badlands: erosion rates and implications for landform evolution

Characteristic badlands are incised into Plio‐Pleistocene clays in Basilicata, southern Italy, creating steep, scarp slopes with knife‐edge ridges (calanchi) and small dome‐shaped forms (biancane). Erosion pin data for the period 1997–2003 give mean annual erosion rates for dome‐shaped biancane in the range 9–19 mm a^?1, while rates for the calanchi scarps are lower, at 7–10 mm a^?1. The erosion pin data also show a non‐linear relationship with slope angle. Maximum erosion rates coincide with a slope angle of 35°, within an envelope defined by combining the theoretical effects of both rainsplash and surface weathering. Monitoring of surface changes and erosion rates for two 0·5 m² cleared swathes on biancane forms reveals a complex relationship between weathering and erosion. Characteristic forms can develop from large blocks of intact clay bedrock over a time period of less than 30 a. The implications of the measured erosion rates for the landform association of mountain front/pediment/domed inselberg are explored. Copyright © 2005 John Wiley & Sons, Ltd.     

Postglacial topographic evolution of glaciated valleys: a stochastic landscape evolution model

The retreat of valley glaciers has a dramatic effect on the stability of glaciated valleys and exerts a prolonged influence on the subsequent fluvial sediment transport regime. We have studied the evolution of an idealized glaciated valley during the period following retreat of ice using a numerical model. The model incorporates a stochastic process to represent deep‐seated landsliding, non‐linear diffusion to represent shallow landsliding and an approximation of the Bagnold relation to represent fluvial sediment transport. It was calibrated using field data from several recent surveys within British Columbia, Canada. We present ensemble model results and compare them with results from a deterministic linear‐diffusion model to show that explicit representation of large landslides is necessary to reproduce the morphology and channel network structure of a typical postglacial valley. Our model predicts a rapid rate of fluvial sediment transport following deglaciation with a subsequent gradual decline, similar to that inferred for Holocene time. We also describe how changes in the model parameters affect the estimated magnitude and duration of the paraglacial sediment pulse. Copyright © 2005 John Wiley & Sons, Ltd.     

Integrating the LISFLOOD‐FP 2D hydrodynamic model with the CAESAR model: implications for modelling landscape evolution

Landscape evolution models (LEMs) simulate the geomorphic development of river basins over long time periods and large space scales (100s–1000s of years, 100s of km²). Due to these scales they have been developed with simple steady flow models that enable long time steps (e.g. years) to be modelled, but not shorter term hydrodynamic effects (e.g. the passage of a flood wave). Nonsteady flow models that incorporate these hydrodynamic effects typically require far shorter time steps (seconds or less) and use more expensive numerical solutions hindering their inclusion in LEMs. The recently developed LISFLOOD‐FP simplified 2D flow model addresses this issue by solving a reduced form of the shallow water equations using a very simple numerical scheme, thus generating a significant increase in computational efficiency over previous hydrodynamic methods. This leads to potential convergence of computational cost between LEMs and hydrodynamic models, and presents an opportunity to combine such schemes. This paper outlines how two such models (the LEM CAESAR and the hydrodynamic model LISFLOOD‐FP) were merged to create the new CAESAR‐Lisflood model, and through a series of preliminary tests shows that using a hydrodynamic model to route flow in an LEM affords many advantages. The new model is fast, computationally efficient and has a stronger physical basis than a previous version of the CAESAR model. For the first time it allows hydrodynamic effects (tidal flows, lake filling, alluvial fans blocking valley floor) to be represented in an LEM, as well as producing noticeably different results to steady flow models. This suggests that the simplification of using steady flow in existing LEMs may bias their findings significantly. Copyright © 2013 John Wiley & Sons, Ltd.     

Ecohydrological evolution model on riparian vegetation in hyperarid regions and its validation in the lower reach of Tarim River

The evolution of the ecohydrological system driven by external climatic forcing and internal feedbacks between vegetation and hydrology, which is more remarkable in arid and semiarid regions, has attracted substantial research attention in recent years. To examine critically the state‐of‐the‐art assumptions and dynamic equations used in the evolution study of an ecohydrological system, the rule of proceeding from simplicity to complexity should be followed. The riparian vegetation ecohydrological system in hyperarid regions (e.g. the lower Tarim River) can serve as a starting point given its simplicity, which has been seldom examined before in terms of system evolution. Further, the water transfer practice from 2000 to 2006 in the lower Tarim River serves as a valuable prototype experiment for model validation. This is because the remarkable changes in groundwater and vegetation in the area have taken place within a shorter period and thus can be easily observed. In the present study, the ecohydrological evolution model on riparian vegetation (ERV model) in hyperarid regions was proposed by coupling groundwater movement and vegetation dynamics. In the ERV model, the groundwater table serves as a critical feedback variable that determines the vegetation dynamics (colonization and mortality) and is determined by vegetation transpiration other than groundwater movement. The monitored groundwater table by wells and satellite‐observed vegetation coverage from the Moderate Resolution Imaging Spectroradiometer are used for model validation. The simulation results show the good performance of the ERV model with uncalibrated parameters. It was also calibrated manually using a multiobjective method, and the fine‐tuned parameters are close to the uncalibrated ones, indicating the robustness of the model. The analysis shows further that the increased evapotranspiration is substantially due to the water transfer and thus the vegetation growth, which implies the importance of ecohydrological coupling for long‐term hydrological modelling. Copyright © 2011 John Wiley & Sons, Ltd.     

Modified MMF (Morgan–Morgan–Finney) model for evaluating effects of crops and vegetation cover on soil erosion

Modifications are made to the revised Morgan–Morgan–Finney erosion prediction model to enable the effects of vegetation cover to be expressed through measurable plant parameters. Given the potential role of vegetation in controlling water pollution by trapping clay particles in the landscape, changes are also made to the way the model deals with sediment deposition and to allow the model to incorporate particle‐size selectivity in the processes of erosion, transport and deposition. Vegetation effects are described in relation to percentage canopy cover, percentage ground cover, plant height, effective hydrological depth, density of plant stems and stem diameter. Deposition is modelled through a particle fall number, which takes account of particle settling velocity, flow velocity, flow depth and slope length. The detachment, transport and deposition of soil particles are simulated separately for clay, silt and sand. Average linear sensitivity analysis shows that the revised model behaves rationally. For bare soil conditions soil loss predictions are most sensitive to changes in rainfall and soil parameters, but with a vegetation cover plant parameters become more important than soil parameters. Tests with the model using field measurements under a range of slope, soil and crop covers from Bedfordshire and Cambridgeshire, UK, give good predictions of mean annual soil loss. Regression analysis of predicted against observed values yields an intercept value close to zero and a line slope close to 1·0, with a coefficient of efficiency of 0·81 over a range of values from zero to 38·6 t ha^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Climate change and soil erosion in a small alpine lake basin on the Loess Plateau,China

Multi‐proxy indices retrieved from sediments in Lake Chaonaqiu, an alpine lake on the western Loess Plateau (LP) of China, were used to reconstruct a precipitation history over the last ~300 years. The results correlate well with records from tree rings and historical documents in neighboring regions. We show that the lake oscillated between two states, i.e. wetter climatic conditions, which favored denser vegetation cover, and promoted weaker catchment soil erosion; and drier climatic conditions, which lead to less vegetation coverage, correlate with stronger surface soil erosion. Several intensive soil erosion events were identified in the sediment cores, and most of these occurred during decadal/multi‐decadal dry periods, and correlate well with flood events documented in historical literature. The results of this study show that soil erosion by flood events is particularly intense during dry periods, and further highlights the role of vegetation cover in the conservation of water and soil in small lake basins on the Chinese LP. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluation of the LISEM soil erosion model in two catchments in the East African Highlands

Under increasing population pressure, soil erosion has become a threat in the East African Highlands, and erosion modelling can be useful to quantify this threat. To test its applicability for this region, the LISEM soil erosion model was applied to two small catchments, one in the Usumbara Mountains, Tanzania, and the other on the slopes of Mount Kenya. Input data for the model were collected in both catchments, as were data on runoff and erosion that were used for calibration and validation of the model. LISEM was first calibrated on catchment outlet data, and afterwards simulated spatial patterns of erosion were compared to available erosion data. The results showed that LISEM can, after calibration, give good discharge predictions for some events, but not for all. However, LISEM generally overpredicted soil loss from the catchments. Comparison with observed erosion patterns did not show overprediction, but according to the model, erosion was more widespread than was observed. There are several reasons for these discrepancies. First, it is difficult to obtain enough accurate data to run the model, such as accurate maps, rainfall data and soil and plant characteristics. Second, it is also difficult to obtain accurate data to evaluate the performance of the model, either for the catchment outlet or spatially, therefore observed erosion rates are also uncertain. Third, the model could not deal correctly with complex events, i.e. those having double rainfall peaks, and might also have difficulties with catchment characteristics such as soil type and the complexity of land use. Finally, LISEM could not deal with events in which throughflow or baseflow played a role, which was to be expected since those processes are not simulated by LISEM. Nevertheless, LISEM could be calibrated to give good discharge predictions for some events, and also gave reasonable results when compared to data obtained from erosion plots. Furthermore, only complex, distributed, storm‐based models such as LISEM can give spatial predictions for single storms. Therefore, it is concluded that if the aim is spatial prediction on an event basis, there is no alternative to complex erosion models such as LISEM, but if the aim is to predict average annual erosion, the data‐demanding, physically based LISEM erosion model may not be the most appropriate model. Copyright © 2005 John Wiley & Sons, Ltd.     

Effects of spatially structured vegetation patterns on hillslope erosion in a semiarid Mediterranean environment: a simulation study

A general trend of decreasing soil loss rates with increasing vegetation cover fraction is widely accepted. Field observations and experimental work, however, show that the form of the cover‐erosion function can vary considerably, in particular for low cover conditions that prevail on arid and semiarid hillslopes. In this paper the structured spatial distribution of the vegetation cover and associated soil attributes is proposed as one of the possible causes of variation in cover–erosion relationships, in particular in dryland environments where patchy vegetation covers are common. A simulation approach was used to test the hypothesis that hillslope discharge and soil loss could be affected by variation in the spatial correlation structure of coupled vegetation cover and soil patterns alone. The Limburg Soil Erosion Model (LISEM) was parameterized and verified for a small catchment with discontinuous vegetation cover at Rambla Honda, SE Spain. Using the same parameter sets LISEM was subsequently used to simulate water and sediment fluxes on 1 ha hypothetical hillslopes with simulated spatial distributions of vegetation and soil parameters. Storms of constant rainfall intensity in the range of 30–70 mm h^?1 and 10–30 min duration were applied. To quantify the effect of the spatial correlation structure of the vegetation and soil patterns, predicted discharge and soil loss rates from hillslopes with spatially structured distributions of vegetation and soil parameters were compared with those from hillslopes with spatially uniform distributions. The results showed that the spatial organization of bare and vegetated surfaces alone can have a substantial impact on predicted storm discharge and erosion. In general, water and sediment yields from hillslopes with spatially structured distributions of vegetation and soil parameters were greater than from identical hillslopes with spatially uniform distributions. Within a storm the effect of spatially structured vegetation and soil patterns was observed to be highly dynamic, and to depend on rainfall intensity and slope gradient. Copyright © 2005 John Wiley & Sons, Ltd.     

The effect of landform variation on vegetation patterning and related sediment dynamics

Semi‐arid ecosystems are often spatially self‐organized in typical patterns of vegetation bands with high plant cover interspersed with bare soil areas, also known as ‘tiger bush’. In modelling studies, most often, straight planar slopes were used to analyse vegetation patterning. The effect of slope steepness has been investigated widely, and some studies investigated the effects of microtopography and hillslope orientation. However, at the larger catchment scale, the overall form of the landscape may affect vegetation patterning and these more complex landscapes are much more prevalent than straight slopes. Hence, our objective was to determine the effect of landform variation on vegetation patterning and sediment dynamics. We linked two well‐established models that simulate (a) plant growth, death and dispersal of vegetation, and (b) erosion and sedimentation dynamics. The model was tested on a straight planar hillslope and then applied to (i) a set of simple synthetic topographies with varying curvature and (ii) three more complex, real‐world landscapes of distinct morphology. Results show banded vegetation patterning on all synthetic topographies, always perpendicular to the slope gradient. Interestingly, we also found that movement of bands – a debated phenomenon – seems to be dependent on curvature. Vegetation banding was simulated on the slopes of the alluvial fan and along the valley slopes of the dissected and rolling landscapes. In all landscapes, local valleys developed a full vegetation cover induced by water concentration, which is consistent with observations worldwide. Finally, banded vegetation patterns were found to reduce erosion significantly as compared to other vegetation configurations. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.