No abstract

The stream power equation is commonly used to model river incision into bedrock. Although specific conditions allow an analytical approach, finite difference methods (FDMs) are most frequently used to solve this equation. FDMs inevitably suffer from numerical smearing which may affect their suitability for transient river incision modeling. We propose the use of a finite volume method (FVM) which is total variation diminishing (TVD) to simulate river incision in a more accurate way. The TVD_FVM is designed to simulate sharp discontinuities, making it very suitable to simulate river incision pulses. We show that the TVD_FVM is much better capable of preserving propagating knickpoints than FDMs, using Niagara Falls as an example. Comparison of numerical results obtained using the TVD_FVM with analytical solutions shows a very good agreement. Furthermore, the uncertainty associated with parameter calibration is dramatically reduced when the TVD_FVM is applied. The high accuracy of the TVD_FDM allows correct simulation of transient incision waves as a consequence of older uplift pulses. This implies that the TVD_FVM is much more suitable than FDMs to reconstruct regional uplift histories from current river profile morphology and to simulate river incision processes in general.

Landscape evolution models (LEMs) allow the study of earth surface responses to changing climatic and tectonic forcings. While much effort has been devoted to the development of LEMs that simulate a wide range of processes, the numerical accuracy of these models has received less attention. Most LEMs use first-order accurate numerical methods that suffer from substantial numerical diffusion. Numerical diffusion particularly affects the solution of the advection equation and thus the simulation of retreating landforms such as cliffs and river knickpoints. This has potential consequences for the integrated response of the simulated landscape. Here we test a higher-order flux-limiting finite volume method that is total variation diminishing (TVD-FVM) to solve the partial differential equations of river incision and tectonic displacement. We show that using the TVD-FVM to simulate river incision significantly influences the evolution of simulated landscapes and the spatial and temporal variability of catchment-wide erosion rates. Furthermore, a two-dimensional TVD-FVM accurately simulates the evolution of landscapes affected by lateral tectonic displacement, a process whose simulation was hitherto largely limited to LEMs with flexible spatial discretization. We implement the scheme in TTLEM (TopoToolbox Landscape Evolution Model), a spatially explicit, raster-based LEM for the study of fluvially eroding landscapes in TopoToolbox 2.Published by Copernicus Publications on behalf of the European Geosciences Union.

Abstract. Landscape evolution models can be used to assess the impact of rainfall variability on bedrock river incision over millennial timescales. However, isolating the role of rainfall variability remains difficult in natural environments, in part because environmental controls on river incision such as lithological heterogeneity are poorly constrained. In this study, we explore spatial differences in the rate of bedrock river incision in the Ecuadorian Andes using three different stream power models. A pronounced rainfall gradient due to orographic precipitation and high lithological heterogeneity enable us to explore the relative roles of these controls. First, we use an area-based stream power model to scrutinize the role of lithological heterogeneity in river incision rates. We show that lithological heterogeneity is key to predicting the spatial patterns of incision rates. Accounting for lithological heterogeneity reveals a nonlinear relationship between river steepness, a proxy for river incision, and denudation rates derived from cosmogenic radionuclide (CRNs). Second, we explore this nonlinearity using runoff-based and stochastic-threshold stream power models, combined with a hydrological dataset, to calculate spatial and temporal runoff variability. Statistical modeling suggests that the nonlinear relationship between river steepness and denudation rates can be attributed to a spatial runoff gradient and incision thresholds. Our findings have two main implications for the overall interpretation of CRN-derived denudation rates and the use of river incision models: (i) applying sophisticated stream power models to explain denudation rates at the landscape scale is only relevant when accounting for the confounding role of environmental factors such as lithology, and (ii) spatial patterns in runoff due to orographic precipitation in combination with incision thresholds explain part of the nonlinearity between river steepness and CRN-derived denudation rates. Our methodology can be used as a framework to study the coupling between river incision, lithological heterogeneity and climate at regional to continental scales.

No abstract

Abstract. Laboratory-scale experiments of erosion have demonstrated that landscapes have a natural (or intrinsic) response time to a change in precipitation rate. In the last few decades there has been growth in the development of numerical models that attempt to capture landscape evolution over long timescales. However, there is still an uncertainty regarding the validity of the basic assumptions of mass transport that are made in deriving these models. In this contribution we therefore return to a principal assumption of sediment transport within the mass balance for surface processes; we explore the sensitivity of the classic end-member landscape evolution models and the sediment fluxes they produce to a change in precipitation rates. One end-member model takes the mathematical form of a kinetic wave equation and is known as the stream power model, in which sediment is assumed to be transported immediately out of the model domain. The second end-member model is the transport model and it takes the form of a diffusion equation, assuming that the sediment flux is a function of the water flux and slope. We find that both of these end-member models have a response time that has a proportionality to the precipitation rate that follows a negative power law. However, for the stream power model the exponent on the water flux term must be less than one, and for the transport model the exponent must be greater than one, in order to match the observed concavity of natural systems. This difference in exponent means that the transport model generally responds more rapidly to an increase in precipitation rates, on the order of 10 5 years for post-perturbation sediment fluxes to return to within 50 % of their initial values, for theoretical landscapes with a scale of 100 × 100 km. Additionally from the same starting conditions, the amplitude of the sediment flux perturbation in the transport model is greater, with much larger sensitivity to catchment size. An important finding is that both models respond more quickly to a wetting event than a drying event, and we argue that this asymmetry in response time has significant implications for depositional stratigraphies. Finally, we evaluate the extent to which these constraints on response times and sediment fluxes from simple models help us understand the geological record of landscape response to rapid environmental changes in the past, such as the Paleocene-Eocene thermal maximum (PETM). In the Spanish Pyrenees, for instance, a relatively rapid (10 to 50 kyr) duration of the deposition of gravel is observed for a climatic shift that is thought to be towards increased precipitation rates. We suggest that the rapid response observed is more easily explained through a diffusive transport model because (1) the model has a faster response time, which is consistent with the documented stratigraphic data, (2) there is a high-amplitude spike in sediment flux, and (3) the assumption of instantaneous transport is difficult to justify for the transport of large grain sizes as an alluvi...

Abstract. Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics.

No abstract