Numerical modelling landscape and sediment flux response to precipitation rate change

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 a growth in the development of numerical models that attempt to capture landscape evolution over long time-scales. Recently, a sub-set of these numerical models have been used to invert river profiles for past tectonic conditions even during variable climatic conditions. However, there is still an uncertainty over validity of the basic assumption of mass transport that are made in deriving these models. In this contribution we therefore return to a principle assumption of sediment transport within the mass balance for surface processes, and explore the sensitivity of the classic end-member landscape evolution models to 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, where sediment is assumed to be transported immediately out of the model domain. The second end-member model takes the form of a diffusion equation, and assumes 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. For the stream power model the exponent on the water flux term must be less than one, and for the sediment 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 sediment transport model responds more rapidly to an increase in precipitation rates, on the order of 105 years for a landscape with a scale of 105 m. In nature, landscape response times to a rapid environmental change have been estimated for events such as the Paleocene-Eocene thermal maximum (PETM). In the Spanish Pyrenees, a relatively rapid, 20 to 100 kyr, duration of deposition of gravel during the PETM is observed for a climatic shift that is thought to be towards increased precipitation rates. We suggest the rapid response observed is more easily explained through a diffusive sediment transport model, as (1) this model has a faster response time, consistent with the documented stratigraphic data, and (2) the assumption of instantaneous transport is difficult to justify for the transport of large grain sizes as an alluvial bed-load.