Constraining the stream power law: a novel approach combining a landscape evolution model and an inversion method

In the past few decades, many studies have been dedicated to the understanding of the interactions between tectonics and erosion, in many instances through the use of numerical models of landscape evolution. Among the numerous parameterizations that have been developed to predict river channel evolution, the stream power law, which links erosion rate to drainage area and slope, remains the most widely used. Despite its simple formulation, its power lies in its capacity to reproduce many of the characteristic features of natural systems (the concavity of river profile, the propagation of knickpoints, etc.). However, the three main coefficients that are needed to relate erosion rate to slope and drainage area in the stream power law remain poorly constrained. In this study, we present a novel approach to constrain the stream power law coefficients under the detachment-limited mode by combining a highly efficient landscape evolution model, FastScape, which solves the stream power law under arbitrary geometries and boundary conditions and an inversion algorithm, the neighborhood algorithm. A misfit function is built by comparing topographic data of a reference landscape supposedly at steady state and the same landscape subject to both uplift and erosion over one time step. By applying the method to a synthetic landscape, we show that different landscape characteristics can be retrieved, such as the concavity of river profiles and the steepness index. When applied on a real catchment (in the Whataroa region of the South Island in New Zealand), this approach provides well-resolved constraints on the concavity of river profiles and the distribution of uplift as a function of distance to the Alpine Fault, the main active structure in the area.