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.     

Orogen‐scale drainage network evolution and response to erodibility changes: insights from numerical experiments

The continuous feedbacks among tectonics, surface processes, and climate are reflected in the distribution of catchments on active mountain ranges. Previous studies have shown a regularity of valley spacing across mountain ranges worldwide, but the origin of this geomorphological feature is currently not well known. In this work, we use a landscape evolution model to investigate the process of fluvial network organization and the evolution of regular ridge‐and‐valley patterns on simulated mountain ranges. In particular, we investigate the behavior of such patterns when subjected to a perturbation in landscape processes from a previous steady state, resulting from a sudden variation in the pattern of bedrock erodibility, from homogeneous to a gradient. We analyze the time evolution of the mean ratio λ' between the linear spacing of adjacent valleys and the half width of the mountain range. We show how a valley spacing ratio of ~0.5 is first achieved at steady state under uniform bedrock erodibility. After applying the gradient of bedrock erodibility across the landscape, we observe that λ' first increases and then decreases to a new steady‐state value that is smaller than the original value. A detailed analysis of the simulations, through observations of surface ‘snapshots’ at repeated time intervals, allows to gain some insight into the mechanisms governing this fluvial network reorganization process, driven by the migration of the main divide toward the side characterized by lower bedrock erodibility. On both sides of the range the new steady‐state valley spacing is obtained through mechanisms of catchment reorganization and competition between adjacent fluvial networks. In particular, catchment reorganization is characterized by the growth of smaller catchments between shrinking larger catchments on the side with lower erodibility, and the growth of larger catchments on the side with higher erodibility. Copyright © 2014 John Wiley & Sons, Ltd.     

Field evidence for the influence of weathering on rock erodibility and channel form in bedrock rivers

Erosion processes in bedrock‐floored rivers shape channel cross‐sectional geometry and the broader landscape. However, the influence of weathering on channel slope and geometry is not well understood. Weathering can produce variation in rock erodibility within channel cross‐sections. Recent numerical modeling results suggest that weathering may preferentially weaken rock on channel banks relative to the thalweg, strongly influencing channel form. Here, we present the first quantitative field study of differential weathering across channel cross‐sections. We hypothesize that average cross‐section erosion rate controls the magnitude of this contrast in weathering between the banks and the thalweg. Erosion rate, in turn, is moderated by the extent to which weathering processes increase bedrock erodibility. We test these hypotheses on tributaries to the Potomac River, Virginia, with inferred erosion rates from ~0.1 m/kyr to >0.8 m/kyr, with higher rates in knickpoints spawned by the migratory Great Falls knickzone. We selected nine channel cross‐sections on three tributaries spanning the full range of erosion rates, and at multiple flow heights we measured (1) rock compressive strength using a Schmidt hammer, (2) rock surface roughness using a contour gage combined with automated photograph analysis, and (3) crack density (crack length/area) at three cross‐sections on one channel. All cross‐sections showed significant (p < 0.01 for strength, p < 0.05 for roughness) increases in weathering by at least one metric with height above the thalweg. These results, assuming that the weathered state of rock is a proxy for erodibility, indicate that rock erodibility varies inversely with bedrock inundation frequency. Differences in weathering between the thalweg and the channel margins tend to decrease as inferred erosion rates increase, leading to variations in channel form related to the interplay of weathering and erosion rate. This observation is consistent with numerical modeling that predicts a strong influence of weathering‐related erodibility on channel morphology. Copyright © 2017 John Wiley & Sons, Ltd.     

Lateral bedrock erosion and valley formation in a heterogeneously layered landscape,Northeast Kansas

In this study, we present direct field measurements of modern lateral and vertical bedrock erosion during a 2-year study period, and optically stimulated luminescence (OSL) ages of fluvial material capping a flat bedrock surface at Kings Creek located in northeast Kansas, USA. These data provide insight into rates and mechanisms of bedrock erosion and valley-widening in a heterogeneously layered limestone-shale landscape. Lateral bedrock erosion outpaced vertical incision during our 2-year study period. Modern erosion rates, measured at erosion pins in limestone and shale bedrock reveal that shale erosion rate is a function of wetting and drying cycles, while limestone erosion rate is controlled by discharge and fracture spacing. Variability in fracture spacing amongst field sites controls the size of limestone block collapse into the stream, which either allowed continued lateral erosion following rapid detachment and transport of limestone blocks, or inhibited lateral erosion due to limestone blocks that protected the valley wall from further erosion. The OSL ages of fluvial material sourced from the strath terrace were older than any material previously dated at our study site and indicate that Kings Creek was actively aggrading and incising throughout the late Pleistocene. Coupling field measurements and observations with ages of fluvial terraces can be useful to investigate the timing and processes linked to how bedrock rivers erode laterally over time to form wide bedrock valleys.     

Amphitheatre-headed canyons of Southern Utah: Stratigraphic control of canyon morphology

Amphitheatre-headed canyons are common on Earth and Mars and researchers have long sought to draw inferences about canyon-forming processes from the morphology of canyon heads and associated knickpoints, often suggesting that amphitheatre heads indicate erosion by groundwater seepage erosion. However, the conditions and processes that lead to amphitheatre-headed canyon formation have been debated for many years. We consider two hypotheses that attribute the amphitheatre-headed canyon formation to fluvial erosion of strong-over-weak stratigraphy or, alternatively, groundwater spring discharge and seepage erosion. A spatial analysis of canyon-form distribution with respect to local stratigraphy along the Escalante River and on Tarantula Mesa, Utah indicates that canyon form is most closely related to variations in local sedimentary rock strata, rather than inferred groundwater spring intensity. Lateral facies variations that affect the continuity of strong layers can induce or disrupt the formation of amphitheatres. Furthermore, we find that amphitheatre retreat rate is dictated by the interaction of fluvial processes downstream of the amphitheatre headwalls and stratigraphy, rather than waterfall and groundwater processes that likely importantly influence headwall form. We conclude that fluvial erosion of strong-over-weak stratigraphic layering alone is sufficient to form amphitheatres at knickpoints and canyon heads. Thus, we re-affirm that formation process should not be inferred from canyon-head morphology, particularly where a strong-over-weak layering is known or plausible. © 2020 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.     

Developing,choosing and using landscape evolution models to inform field‐based landscape reconstruction studies

Landscape evolution models (LEMs) are an increasingly popular resource for geomorphologists as they can operate as virtual laboratories where the implications of hypotheses about processes over human to geological timescales can be visualized at spatial scales from catchments to mountain ranges. Hypothetical studies for idealized landscapes have dominated, although model testing in real landscapes has also been undertaken. So far however, numerical landscape evolution models have rarely been used to aid field‐based reconstructions of the geomorphic evolution of actual landscapes. To help make this use more common, we review numerical landscape evolution models from the point of view of model use in field reconstruction studies. We first give a broad overview of the main assumptions and choices made in many LEMs to help prospective users select models appropriate to their field situation. We then summarize for various timescales which data are typically available and which models are appropriate. Finally, we provide guidance on how to set up a model study as a function of available data and the type of research question. Copyright © 2017 John Wiley & Sons, Ltd.     

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 legacy of initial conditions in landscape evolution

Landscapes subject to constant forcing tend to evolve toward equilibrium states in which individual landforms have similar characteristics. Yet, even in landscapes at or near equilibrium, there can be significant variability among individual landforms. Furthermore, sites subject to similar processes and conditions can have different mean landform characteristics. This variability is often ascribed to on‐going transient evolution, or to heterogeneity in processes, material properties, forcing, or boundary conditions. Three surprising outcomes of landform evolution models suggest, however, that such variability could arise in equilibrium landscapes without any heterogeneity in the physical processes shaping the topography. First, homogeneous systems subjected to constant forcing can generate a heterogeneous distribution of equilibrium landforms. Second, even simple non‐linear systems can have multiple stable equilibrium states. Third, evolving landscapes can exhibit path dependence and hysteresis. We show how these three mechanisms can produce variability in landforms that arises from the characteristics of the initial topographic surface rather than from heterogeneity in geomorphic processes. Numerical experiments on the formation of low‐order fluvial valleys and transportational cyclic steps in erodible streambeds illustrate why it is important to consider the influence of initial conditions when comparing models with natural topography, estimating the uncertainty of model predictions, and studying how landscapes respond to disturbances. Copyright © 2011 John Wiley & Sons, Ltd.     

Long‐term landscape evolution: linking tectonics and surface processes

Research in landscape evolution over millions to tens of millions of years slowed considerably in the mid‐20th century, when Davisian and other approaches to geomorphology were replaced by functional, morphometric and ultimately process‐based approaches. Hack's scheme of dynamic equilibrium in landscape evolution was perhaps the major theoretical contribution to long‐term landscape evolution between the 1950s and about 1990, but it essentially ‘looked back’ to Davis for its springboard to a viewpoint contrary to that of Davis, as did less widely known schemes, such as Crickmay's hypothesis of unequal activity. Since about 1990, the field of long‐term landscape evolution has blossomed again, stimulated by the plate tectonics revolution and its re‐forging of the link between tectonics and topography, and by the development of numerical models that explore the links between tectonic processes and surface processes. This numerical modelling of landscape evolution has been built around formulation of bedrock river processes and slope processes, and has mostly focused on high‐elevation passive continental margins and convergent zones; these models now routinely include flexural and denudational isostasy. Major breakthroughs in analytical and geochronological techniques have been of profound relevance to all of the above. Low‐temperature thermochronology, and in particular apatite fission track analysis and (U–Th)/He analysis in apatite, have enabled rates of rock uplift and denudational exhumation from relatively shallow crustal depths (up to about 4 km) to be determined directly from, in effect, rock hand specimens. In a few situations, (U–Th)/He analysis has been used to determine the antiquity of major, long‐wavelength topography. Cosmogenic isotope analysis has enabled the determination of the ‘ages’ of bedrock and sedimentary surfaces, and/or the rates of denudation of these surfaces. These latter advances represent in some ways a ‘holy grail’ in geomorphology in that they enable determination of ‘dates and rates’ of geomorphological processes directly from rock surfaces. The increasing availability of analytical techniques such as cosmogenic isotope analysis should mean that much larger data sets become possible and lead to more sophisticated analyses, such as probability density functions (PDFs) of cosmogenic ages and even of cosmogenic isotope concentrations (CICs). PDFs of isotope concentrations must be a function of catchment area geomorphology (including tectonics) and it is at least theoretically possible to infer aspects of source area geomorphology and geomorphological processes from PDFs of CICs in sediments (‘detrital CICs’). Thus it may be possible to use PDFs of detrital CICs in basin sediments as a tool to infer aspects of the sediments' source area geomorphology and tectonics, complementing the standard sedimentological textural and compositional approaches to such issues. One of the most stimulating of recent conceptual advances has followed the considerations of the relationships between tectonics, climate and surface processes and especially the recognition of the importance of denudational isostasy in driving rock uplift (i.e. in driving tectonics and crustal processes). Attention has been focused very directly on surface processes and on the ways in which they may ‘drive’ rock uplift and thus even influence sub‐surface crustal conditions, such as pressure and temperature. Consequently, the broader geoscience communities are looking to geomorphologists to provide more detailed information on rates and processes of bedrock channel incision, as well as on catchment responses to such bedrock channel processes. More sophisticated numerical models of processes in bedrock channels and on their flanking hillslopes are required. In current numerical models of long‐term evolution of hillslopes and interfluves, for example, the simple dependency on slope of both the fluvial and hillslope components of these models means that a Davisian‐type of landscape evolution characterized by slope lowering is inevitably ‘confirmed’ by the models. In numerical modelling, the next advances will require better parameterized algorithms for hillslope processes, and more sophisticated formulations of bedrock channel incision processes, incorporating, for example, the effects of sediment shielding of the bed. Such increasing sophistication must be matched by careful assessment and testing of model outputs using pre‐established criteria and tests. Confirmation by these more sophisticated Davisian‐type numerical models of slope lowering under conditions of tectonic stability (no active rock uplift), and of constant slope angle and steady‐state landscape under conditions of ongoing rock uplift, will indicate that the Davis and Hack models are not mutually exclusive. A Hack‐type model (or a variant of it, incorporating slope adjustment to rock strength rather than to regolith strength) will apply to active settings where there is sufficient stream power and/or sediment flux for channels to incise at the rate of rock uplift. Post‐orogenic settings of decreased (or zero) active rock uplift would be characterized by a Davisian scheme of declining slope angles and non‐steady‐state (or transient) landscapes. Such post‐orogenic landscapes deserve much more attention than they have received of late, not least because the intriguing questions they pose about the preservation of ancient landscapes were hinted at in passing in the 1960s and have recently re‐surfaced. As we begin to ask again some of the grand questions that lay at the heart of geomorphology in its earliest days, large‐scale geomorphology is on the threshold of another ‘golden’ era to match that of the first half of the 20th century, when cyclical approaches underpinned virtually all geomorphological work. Copyright © 2007 John Wiley & Sons, Ltd.     

Long-term landscape evolution in Australia

The main features of the Australian physical landscape are of the order of 10⁷-10⁸ years old. This contradicts the widely held view that little of the Earth's topography predates the Quaternary and that erosion cycles are carried to planation within tens of millions of years. Much of the Australian landscape must have developed over similar timescales to that of the tectonic evolution of the continent itself. The study of the geomorphology of such ancient terrains may therefore be seriously deficient unless it is considered within the context of continental-scale tectonic development. Application of this approach shows that there are strong links between the geomorphology of Australia and plate movements, ocean spreading, plate convergence, tectonostratigraphic terranes, orogenesis and epeirogenesis. The most important factor contributing to the survival of ancient landscapes in Australia is the low rate of denudation which the continent has experienced during the Mesozoic and Cenozoic. This is largely a consequence of orogenic stability, although the absence of significant Quaternary glaciation may also be of importance. However, in order for landforms to have survived over such timespans, denudation must not only have been low, but must also have been highly localized over space and time. This has been the case both on a regional scale, with long-term denudation rates of 0-2 m Ma^?1 in central Australia contrasting with higher rates along the continental margins, and on a local scale, with denudation confined to valleys, leaving divides and interfluves almost unscathed.     

Coupled numerical–analytical approach to landscape evolution modeling

The Earth's topography is shaped by surface processes that operate on various scales. In particular, river processes control landscape dynamics over large length scales, whereas hillslope processes control the dynamics over smaller length scales. This scale separation challenges numerical treatments of landscape evolution that use space discretization. Large grid spacing cannot account for the dynamics of water divides that control drainage area competition, and erosion rate and slope distribution. Small grid spacing that properly accounts for divide dynamics is computationally inefficient when studying large domains. Here we propose a new approach for landscape evolution modeling that couples irregular grid‐based numerical solutions for the large‐scale fluvial dynamics and continuum‐based analytical solutions for the small‐scale fluvial and hillslope dynamics. The new approach is implemented in the landscape evolution model DAC (divide and capture). The geometrical and topological characteristics of DAC's landscapes show compatibility with those of natural landscapes. A comparative study shows that, even with large grid spacing, DAC predictions fit well an analytical solution for divide migration in the presence of horizontal advection of topography. In addition, DAC is used to study some outstanding problems in landscape evolution. (i) The time to steady‐state is investigated and simulations show that steady‐state requires much more time to achieve than predicted by fixed area calculations, due to divides migration and persistent reorganization of low‐order streams. (ii) Large‐scale stream captures in a strike‐slip environment are studied and show a distinct pattern of erosion rates that can be used to identify recent capture events. (iii) Three tectono‐climatic mechanisms that can lead to asymmetric mountains are studied. Each of the mechanisms produces a distinct morphology and erosion rate distribution. Application to the Southern Alps of New Zealand suggests that tectonic advection, precipitation gradients and non‐uniform tectonic uplift act together to shape the first‐order topography of this mountain range. Copyright © 2014 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.     

Late quaternary speleogenesis and landscape evolution in the northern Apennine evaporite areas

Gypsum beds host the majority of the caves in the north‐eastern flank of the Apennines, in the Emilia Romagna region (Italy). More than six hundred of these caves have been surveyed, including the longest known epigenic gypsum cave systems in the world (Spipola‐Acquafredda, ~11 km). Although this area has been intensively studied from a geological point of view, the age of the caves has never been investigated in detail. The rapid dissolution of gypsum and uplift history of the area have led to the long‐held view that speleogenesis commenced only during the last 130 000 years. Epigenic caves only form when the surface drainage system efficiently conveys water into the underground. In the study area, this was achieved after the dismantling of most of the impervious sediments covering the gypsum and the development of protovalleys and sinkholes. The time necessary for these processes can by constrained by understanding when caves were first formed. The minimum age of karst voids can be indirectly estimated by dating the infilling sediments. U–Th dating of carbonate speleothems growing in gypsum caves has been applied to 20 samples from 14 different caves from the Spipola‐Acquafredda, Monte Tondo‐Re Tiberio, Stella‐Rio Basino, Monte Mauro, and Castelnuovo systems. The results show that: (i) caves have been forming since at least ~600 kyr ago; (ii) the peak of speleogenesis was reached during relatively cold climate stages, when rivers formed terraces at the surface and aggradation caused paragenesis in the stable cave levels; (iii) ~200 000 years were necessary for the dismantling of most of the sediments covering the karstifiable gypsum and the development of a surface mature drainage network. Besides providing a significant contribution to the understanding of evaporite karst evolution in the Apennines, this study refines our knowledge on the timescale of geomorphological processes in a region affected by rapid uplifting. Copyright © 2016 John Wiley & Sons, Ltd.     

The influence of bedrock orientation on the landscape evolution,surface morphology and denudation (10Be) at the Niesen,Switzerland

Landscape evolution and surface morphology in mountainous settings are a function of the relative importance between sediment transport processes acting on hillslopes and in channels, modulated by climate variables. The Niesen nappe in the Swiss Penninic Prealps presents a unique setting in which opposite facing flanks host basins underlain by identical lithologies, but contrasting litho‐tectonic architectures where lithologies either dip parallel to the topographic slope or in the opposite direction (i.e. dip slope and non‐dip slope). The north‐western facing Diemtigen flank represents such a dip slope situation and is characterized by a gentle topography, low hillslope gradients, poorly dissected channels, and it hosts large landslides. In contrast, the south‐eastern facing Frutigen side can be described as non‐dip slope flank with deeply incised bedrock channels, high mean hillslope gradients and high relief topography. Results from morphometric analysis reveal that noticeable differences in morphometric parameters can be related to the contrasts in the relative importance of the internal hillslope‐channel system between both valley flanks. While the contrasting dip‐orientations of the underlying flysch bedrock has promoted hillslope and channelized processes to contrasting extents and particularly the occurrence of large landslides on the dip slope flank, the flank averaged beryllium‐10 (¹⁰Be)‐derived denudation rates are very similar and range between 0.20 and 0.26 mm yr^?1. In addition, our denudation rates offer no direct relationship to basin's slope, area, steepness or concavity index, but reveal a positive correlation to mean basin elevation that we interpret as having been controlled by climatically driven factors such as frost‐induced processes and orographic precipitation. Our findings illustrate that while the landscape properties in this part of the northern Alpine border can mainly be related to the tectonic architecture of the underlying bedrock, the denudation rates have a strong orographic control through elevation dependent mean annual temperature and precipitation. Copyright © 2013 John Wiley & Sons, Ltd.     

Topographic signatures of progressive glacial landscape transformation

The Pleistocene glaciations left a distinct topographic footprint in mountain ranges worldwide. The geometric signature of glacial topography has been quantified in various ways, but the temporal development of landscape metrics has not been traced in a landscape evolution model so far. However, such information is needed to interpret the degree of glacial imprint in terms of the integrated signal of temporal and spatial variations in erosion as a function of glacial occupation time. We apply a surface process model for cold-climate conditions to an initially fluvial mountain range. By exploring evolving topographic patterns in model time series, we determine locations where topographic changes reach a maximum and where the initial landscape persists. The signal of glacial erosion, expressed by the overdeepening of valleys and the steepening of valley flanks, develops first at the glacier front and migrates upstream with ongoing glacial erosion. This leads to an increase of mean channel slope and its variance. Above steep flanks and head-walls, however, the observed mean channel slope remains similar to the mean channel slope of the initial fluvial topography. This leads to a characteristic turning point in the channel slope–elevation distribution above the equilibrium line altitude, where a transition from increasing to decreasing channel slope with elevation occurs. We identify this turning point and a high channel slope variance as diagnostic features to quantify glacial imprint. Such features are abundant in glacially imprinted mid-latitude mountain ranges such as the Eastern Alps. By analysing differently glaciated parts of the mountain range, we observe a decreasing clarity of this diagnostic morphometric property with decreasing glacial occupation. However, catchments of the unglaciated eastern fringe of the Alps also feature turning points in their channel slope–elevation distributions, but in contrast to the glaciated domain, the variance of channel slope is small at all elevation levels.     

The impact of different rainfall products on landscape modelling simulations

Rainfall products can contain significantly different spatiotemporal estimates, depending on their underlying data and final constructed resolution. Commonly used products, such as rain gauges, rain gauge networks, and weather radar, differ in their information content regarding intensities, spatial variability, and natural climatic variability, therefore producing different estimates. Landscape evolution models (LEMs) simulate the geomorphic changes in landscapes, and current models can simulate timeframes from event level to millions of years and some use rainfall inputs to drive them. However, the impact of different rainfall products on LEM outputs has never been considered. This study uses the STREAP rainfall generator, calibrated using commonly used rainfall observation products, to produce longer rainfall records than the observations to drive the CAESAR-Lisflood LEM to examine how differences in rainfall products affect simulated landscapes. The results show that the simulation of changes to basin geomorphology is sensitive to the differences between rainfall products, with these differences expressed linearly in discharges but non-linearly in sediment yields. Furthermore, when applied over a 1500-year period, large differences in the simulated long profiles were observed, with the simulations producing greater sediment yields showing erosion extending further downstream. This suggests that the choice of rainfall product to drive LEMs has a large impact on the final simulated landscapes. The combination of rainfall generator model and LEMs represents a potentially powerful method for assessing the impacts of rainfall product differences on landscapes and their short- and long-term evolution. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Complexities of landscape evolution during incision through layered stratigraphy with contrasts in rock strength

Variation in the erodibility of rock units has long been recognized as an important determinant of landscape evolution but has been little studied in landscape evolution models. We use a modified version of the Channel‐Hillslope Integrated Landscape Development (CHILD) model, which explicitly allows for variations in rock strength, to reveal and explore the remarkably rich, complex behavior induced by rock erodibility variations in even very simple geologic settings with invariant climate and tectonics. We study the importance of relative contrasts in erodibility between just two units, the order of these units (whether hard rocks overlie soft or soft rocks overlie hard) and the orientation of the contact between the two units. We emphasize the spatial and temporal evolution of erosion rates, which have important implications for basin analysis, detrital mineral records, and the interpretation of cosmogenic isotope concentrations in detrital samples. Results of the landscape evolution modeling indicate that the stratigraphic order of units in terms of erodibility, the gross orientation of the contact (i.e. dipping away or toward the outlet of the landscape) and the contact dip angle all have measurable effects on landscape evolution, including significant spatial and temporal variations in erosion rates. Steady‐state denudation conditions are unlikely to develop in landscapes with significant contrasts in rock strength in horizontal to moderately tilted rock layers, at least at the scale of the entire landscape. Additionally, our results demonstrate that there is no general relation between rock erodibility and erosion rates in natural settings. Although rock erodibility directly controls the erosion rate constant in our models, it is not uncommon for higher erosion rates to occur in the harder, less erodible rock. Indeed erosion rates may be either greater or less than the rock uplift rate (invariant in time and space in our models) in both hard and soft rocks, depending on the local geology, topography, and the pattern of landscape evolution. Copyright © 2016 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.     

An empirical model of event scale cohesive bank profile evolution

The highly stochastic nature of riverbank erosion has driven the need for spatially explicit empirical models. Detailed bank profile surveys along a meander bend of the Brandywine Creek in Pennsylvania, USA, before and after 28 high flow events over a 2·5 year period are used to develop an empirical model of cohesive bank profile erosion. Two hundred and thirty‐six bank erosion observations are classified as hydraulic erosion or subaerial erosion. Threshold conditions required to initiate bank erosion cannot be defined based on field measurements. Using the near‐bank velocity and the number of freeze–thaw cycles as predictors, regression equations are derived for hydraulic erosion that specify the length, thickness, and location on the bank face of eroded blocks. An empirical discriminant function defines the critical geometry of overhang failures, and the volumes removed by overhang failures are computed using another regression equation. All the regression equations are significant, but have low correlation coefficients, suggesting that cohesive bank erosion has a strong stochastic component. Individual events typically remove small masses of soil (average volume 0·084 m³/m) a few centimeters thick (median = 0·057 m) and a few decimeters in length (median = 0·50 m) from the lower third of the bank. Hydraulic erosion is responsible for 87% of all erosion. When applied to three survey sites not used in its development, the profile model predicts the total volume of erosion with errors of 23%, 5% and 1%. Twenty‐four percent of computed erosion volumes for single events are within 50% of observed volumes at these three sites. Extending the approach to decadal timescales and to entire bends will require three‐dimensional observations of bank failure, and spatially and temporally explicit methods to account for the influence of individual large trees on bank failures and near‐bank hydraulic processes. Copyright © 2009 John Wiley & Sons, Ltd.