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.     

Modelling centennial sediment waves in an eroding landscape – catchment complexity

Sediment flux dynamics in fluvial systems have often been related to changes in external drivers of topography, climate or land cover. It is well known that these dynamics are non‐linear. Recently, model simulations of fluvial activity and landscape evolution have suggested that self‐organization in landscapes can also cause internal complexity in the sedimentary record. In this contribution one particular case of self‐organization is explored in the Sabinal field study area, Spain, where several dynamic zones of sedimentation and incision are observed along the current river bed. Whether these zones can be caused by internal complexity was tested with landscape evolution model (LEM) LAPSUS (Landscape Process Modelling at Multi‐dimensions and Scales). During various 500 year simulations, zones of sedimentation appear to move upstream and downstream in eroding river channels (‘waves’). These waves are visualized and characterized for a range of model settings under constant external forcing, and the self‐organizing process behind their occurrence is analysed. Results indicate that this process is not necessarily related to simplifications in the model and is more generic than the process of bed‐armouring that has recently been recognized as a cause for complexity in LEM simulations. We conclude that autogenic sediment waves are the result of the spatial propagation in time of feedbacks in local transport limited (deposition) and detachment limited (erosion) conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

The shape of patterns to come: from initial formation to long‐term evolution

Many studies focus on the emergence and development of rhythmic landscape patterns. In this contribution we explore the different behaviors found as patterns evolve; the trajectories that patterns exhibit as they transit from infinitesimal‐amplitude perturbation to a statistically steady state (or in some cases to continued statistical evolution). The variety of behaviors observed, either through field and laboratory experiments or numerical modeling, can be reduced to four classes: (a) simple stabilization where predictions based on the initial growth of small perturbations corresponds with the characteristics of patterns observed in nature; (b) significant pattern coarsening en route to saturated wavelength, where non‐linear interactions between finite‐amplitude pattern elements change the geometric properties of a pattern as it approaches steady‐state; (c) perpetual coarsening where the wavelength associated with the emerging pattern continues to grow over time and is only limited by physical boundaries or external constrains; (d) slow evolution toward a different attractor, a novel behavior observed in numerical modeling that involves profound temporal changes in pattern characteristics. Within these classes we also observe generalizable non‐linear behaviors: dependence on initial conditions, the emergence of pattern‐scale variables such as pattern defects, and the presence of multiple stable states. Predicting the shape of patterns to come remains a challenge – one that we suggest requires a range of modeling approaches to address both initial instabilities and the emergent properties of evolving patterns, which involve disparate forms of non‐linear interactions. Consideration of generic system behaviors at the pattern scale could enhance future pattern formation studies, facilitating appropriate pairings of analysis approaches and pattern‐evolution modes. Copyright © 2013 John Wiley & Sons, Ltd.     

Natural experiments in landscape evolution

A natural experiment in landscape evolution is a case study of landform development in which only one element varies significantly, and for which the driving forces, initial conditions, and/or boundary conditions are well constrained. Natural experiments provide a means of testing landscape evolution theory on the large space and time scales to which that theory applies. Natural experiments can involve either steady or transient conditions. Cases with steady conditions allow one to test predictions about the relationships among topography, erosion rates, and various attributes related to climate and material properties. Transient cases are valuable for distinguishing between models whose predictions might be similar, and therefore indistinguishable, under steady conditions. Essential ingredients of a natural experiment include minimal variation in all but one factor, good constraints on timing and/or rates, well‐characterized processes, and high quality topographic data. Other useful ingredients include information about intermediate topographic states (such as a former valley profile revealed by strath terraces), and knowledge of the time history of erosion rates. In order to deepen our understanding of the physics and chemistry of long‐term landscape evolution, there is a pressing need to identify natural experiments and develop the necessary databases to take advantage of them. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

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     

Seismic response of self‐centring hysteretic SDOF systems

The seismic response of single‐degree‐of‐freedom (SDOF) systems incorporating flag‐shaped hysteretic structural behaviour, with self‐centring capability, is investigated numerically. For a SDOF system with a given initial period and strength level, the flag‐shaped hysteretic behaviour is fully defined by a post‐yielding stiffness parameter and an energy‐dissipation parameter. A comprehensive parametric study was conducted to determine the influence of these parameters on SDOF structural response, in terms of displacement ductility, absolute acceleration and absorbed energy. This parametric study was conducted using an ensemble of 20 historical earthquake records corresponding to ordinary ground motions having a probability of exceedence of 10% in 50 years, in California. The responses of the flag‐shaped hysteretic SDOF systems are compared against the responses of similar bilinear elasto‐plastic hysteretic SDOF systems. In this study the elasto‐plastic hysteretic SDOF systems are assigned parameters representative of steel moment resisting frames (MRFs) with post‐Northridge welded beam‐to‐column connections. In turn, the flag‐shaped hysteretic SDOF systems are representative of steel MRFs with newly proposed post‐tensioned energy‐dissipating connections. Building structures with initial periods ranging from 0.1 to 2.0s and having various strength levels are considered. It is shown that a flag‐shaped hysteretic SDOF system of equal or lesser strength can always be found to match or better the response of an elasto‐plastic hysteretic SDOF system in terms of displacement ductility and without incurring any residual drift from the seismic event. Copyright © 2002 John Wiley & Sons, Ltd.     

Possible evidence for underlying non‐linear dynamics in steep‐faced glaciodeltaic progradational ­successions

Steep‐faced glaciodeltaic progradational successions are often studied in order to reconstruct the behaviour of the glacial feeder system, or changes in the sediment sink. This paper analyses the magnitude and frequency of depositional events associated with steep‐faced glacier‐fed progradational successions recorded in Scandinavia and Ireland. The successions exhibit depositional patterns that may be interpreted as a function of underlying non‐linear dynamics. A number of the sequences display fractal scaling in the frequency and thickness of foreset units. Other successions demonstrate chaotic patterns and strong relationships between delta‐front angle and bed thicknesses, suggesting that the progradation of such sequences is self‐organized, and to an extent occurs independently of forcing by the feeder system that provides sediment to the delta front. These patterns of sedimentation appear to be a function of the steepness of the delta front and/or the textural characteristics of the sediment. This paper provides further evidence for the simultaneous presence of order and chaos in Earth surface processes and calls into question the extent to which palaeoenvironmental reconstructions may be made from steep‐faced progradational successions. Copyright © 2000 John Wiley & Sons, Ltd.     

Reciprocal theorems in structural dynamics including initial conditions

For completeness purposes, as well as for practical reasons, this work investigates the well‐known Betti–Rayleigh reciprocal theorem for structural dynamics in a way which includes the effect of initial conditions. It then presents a natural and consistent way for introducing the Duhamel integral for the transient response of a dynamical system through the aforementioned reciprocal theorem. Copyright © 2005 John Wiley & Sons, Ltd.     

Complexity and non‐linearity in earth surface processes – concepts,methods and applications

Complexity has long been recognized and is increasingly becoming mainstream in geomorphology. However, the relative novelty of various concepts and techniques associated to it means that ambiguity continues to surround complexity. In this commentary, we present and discuss a variety of recent contributions that have the potential to help clarify issues and advance the use of complexity in geomorphology. Copyright © 2015 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.     

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.     

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 interaction between hydrology and geomorphology in a landscape simulator experiment

An experimental landscape simulator has been developed which uses a rainfall simulator to create overland flow and erosion. The simulator uses rainfall sprinklers that eliminate rainsplash and an artificial soil which has little cohesion. Experimental landscapes developed in the simulator evolved according to Howard's headward growth model. Elements of Glock's model could be identified during evolution (i.e. initiation and maximum extension), but other stages of this model were not observed (i.e. extension and integration). The Horton concept of cross‐grading and micropiracy and stream piracy was not observed despite the dominance of overland flow, nor the groundwater headward growth mechanism proposed by Dunne, the latter due to experimental design, which eliminated any perched groundwater table. The experimental apparatus produced model landscapes that are scaled‐down analogues of real world processes. Copyright © 2001 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.     

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.     

Linear and non‐linear scaling of the Yangtze River flow

Power law correlation properties of sign and magnitude series have been studied based on the series of observation records of flow of the River Yangtze. The results obtained give improved insight into and understanding of the linear and non‐linear processes of the water cycle. With the newly developed Delayed Vector Variance method and the surrogate test, the documented linkage between the sign series and the linear process, and that between the magnitude series and non‐linear process can be verified. The spectra estimated by detrended fluctuation analysis method show different properties of intra‐annual and inter‐annual correlations in both sign and magnitude series. The linear process behaves as an 1/f noise at a time scale less than about 60 days, but shows features of anti‐persistence in terms of long‐term fluctuation. The magnitudes are clustered in three ways mainly caused by non‐linear processes, i.e. periodic clustering, strong short‐term clustering of 1/f noise at time scales less than 20 days, and long‐term clustering with weak persistence. Copyright © 2007 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.