**Research article**
29 Nov 2019

**Research article** | 29 Nov 2019

# Estimating the disequilibrium in denudation rates due to divide migration at the scale of river basins

Timothée Sassolas-Serrayet Rodolphe Cattin Matthieu Ferry Vincent Godard and Martine Simoes

^{1},

^{1},

^{1},

^{2},

^{3}

**Timothée Sassolas-Serrayet et al.**Timothée Sassolas-Serrayet Rodolphe Cattin Matthieu Ferry Vincent Godard and Martine Simoes

^{1},

^{1},

^{1},

^{2},

^{3}

^{1}Géosciences Montpellier, Université de Montpellier – CNRS, Montpellier, France^{2}Aix Marseille Univ, CNRS, IRD, INRA, Coll France, CEREGE, Aix-en-Provence, France^{3}Université de Paris, Institut de physique du globe de Paris, CNRS, 75005 Paris, France

^{1}Géosciences Montpellier, Université de Montpellier – CNRS, Montpellier, France^{2}Aix Marseille Univ, CNRS, IRD, INRA, Coll France, CEREGE, Aix-en-Provence, France^{3}Université de Paris, Institut de physique du globe de Paris, CNRS, 75005 Paris, France

**Correspondence**: Timothée Sassolas-Serrayet (timothee.sassolas-serrayet@umontpellier.fr)

**Correspondence**: Timothée Sassolas-Serrayet (timothee.sassolas-serrayet@umontpellier.fr)

Received: 29 May 2019 – Discussion started: 26 Jun 2019 – Revised: 04 Oct 2019 – Accepted: 21 Oct 2019 – Published: 29 Nov 2019

Basin-averaged denudation rates may locally exhibit a wide dispersion, even in areas where the topographic steady state is supposedly achieved regionally. This dispersion is often attributed to the accuracy of the data or to some degree of natural variability of local erosion rates which can be related to stochastic processes such as landsliding. Another physical explanation of this dispersion is local and transient disequilibrium between tectonic forcing and erosion at the scale of catchments. Recent studies have shown that basin divide migration can potentially induce such perturbations, and they propose metrics to assess divide mobility based on cross-divide contrasts in headwater topographic features. Here, we use a set of landscape evolution models assuming spatially uniform uplift, rock strength and rainfall to assess the effect of divide mobility on basin-wide denudation rates. We propose using basin-averaged aggressivity metrics based on cross-divide contrasts (1) in channel *χ*, an integral function of position in the channel network; (2) in channel local gradient; and (3) in channel height, measured at a reference drainage area. From our simulations, we show that the metric based on differences in *χ* is the most reliable to diagnose local disequilibrium. The other metrics are more suitable for relatively active tectonic regions such as mountain belts, where contrasts in local gradient and elevation are more important. We find that the ratio of basin denudation associated with drainage migration to uplift can reach a factor of 2, regardless of the imposed uplift rate, erodibility, diffusivity coefficient or critical hillslope gradient. A comparison with field observations in the Great Smoky Mountains (southern Appalachians, USA) underlines the difficulty of using the metric based on *χ*, which depends on the – poorly constrained – elevation of the outlet of the investigated catchment. Regardless of the considered metrics, we show that observed dispersion is controlled by catchment size: a smaller basin may be more sensitive to divide migration and hence to disequilibrium. Our results thus highlight the relevance of divide stability analysis from digital elevation models as a fundamental preliminary step for basin-wide denudation rate studies based on cosmogenic radionuclide concentrations.

- Article
(13623 KB) -
Supplement
(6934 KB) - BibTeX
- EndNote

Topographic steady state, in which average topography is constant over time, is one of the key concepts of modern geomorphology (e.g. Gilbert, 1877; Hack, 1960; Montgomery, 2001). Though simple, this paradigm provides a useful framework to study landscape evolution related to tectonic and/or climatic forcing (e.g. Willett et al., 2001; Reinhart and Ellis, 2015), to spatial variations in rock strength (Perne et al., 2017) or to the geometry of active crustal structures (Lavé and Avouac, 2001; Stolar et al., 2007; Scherler et al., 2014; Le Roux-Mallouf et al., 2015). To define topographic steady state, the temporal and spatial scales of the processes involved are essential parameters. Compared to large-scale geodynamic processes operating over 1–100 Myr timescales, river incision and sediment transport are rapid processes driving landscapes to stable forms over this long timescale, whereas rapid climatic fluctuations during the Quaternary may prevent the occurrence of steady-state conditions in modern landscapes (Whipple, 2001).

The timescale of divide migration has received increasing attention in recent years. Although rivers exhibit a rapid adjustment to tectonic or
climatic changes to maintain their profiles, Whipple et al. (2017) show that
divides continue to migrate over time periods of 10^{6}–10^{7} years as
a response to the same changes. This suggests that long-term transience might
be pervasive in the planar structure of landscapes, even in the absence of
new variations in landscape characteristics or forcings (e.g. tectonic or
climate) (Hasbargen and Paola, 2000, 2003; Pelletier, 2004; Dahlquist et al., 2018). In addition to the influence of spatial variability of rock uplift rate, rock strength or rainfall (e.g. Reiners et al., 2003; Godard et al., 2006; Miller et al., 2013), this long timescale could also explain the persistence of spatial variations in denudation rates observed in tectonically inactive orogens which achieved regional-scale topographic steady state (Willett et al., 2014).

As an example, in the Great Smoky Mountains in the southern Appalachians,
uplift and erosion rates integrated over varying time periods from 10 kyr
to 100 Myr give a similar average magnitude of ca. 0.03 mm yr^{−1} (Matmon et al., 2003a, b; Portenga and Bierman, 2011). These results suggest a regional quasi-topographic steady state over the last ∼180 Myr, maintained by the isostatic response of the thickened crust since the
end of the Appalachian orogeny (Matmon et al., 2003a, b). Beyond this
average value, individual basin-wide denudation rates exhibit a strong
dispersion (up to a factor of 2, Fig. 1), which is not related to spatial
variation in rainfall or in erodibility of the substrate (Matmon et al.,
2003b). In a recent study, Willett et al. (2014) assess divide mobility from
the contrast in the channel head topographic metric *χ*, taken here as a
proxy for steady-state river profile elevation (Perron and Royden, 2012;
Royden and Perron, 2013), and propose an explanation in which a significant
part of the observed dispersion in denudation rates could be due to drainage
divide migration associated with contrasting erosion rates across divides.

More recently, to characterize divide migrations, Forte and Whipple (2018)
introduced other metrics, referred to as “Gilbert metrics” (Gilbert, 1877),
based on the cross-divide contrast in channel local gradient and
height. This last study indeed focused on cross-divide contrasts in
headwater basin shape. Here, we propose extending these approaches by
modelling divide migration and by developing new metrics to assess divide
stability at the scale of the entire watershed, which are an expansion of
the aggressivity metric initially suggested by Willett et al. (2014). We use
these metrics to assess the effect of persistent divide mobility on
basin-averaged erosion rates at a timescale of 10^{4} years. We use numerical
landscape evolution models, taking into account both hillslope diffusion and
fluvial incision. For the sake of simplicity and to avoid the influence of
other factors such as topography, lithology, climate or vegetation, we
restrict our analysis to synthetic orogens with spatially uniform uplift,
rock strength and rainfall. After a brief presentation of our landscape
evolution model (LEM), we describe the methods developed to assess
basin-wide denudation rates and aggressivity metrics, such as average
cross-divide contrasts in channel *χ*, gradient and height. Next, we
investigate transient time and location of morphologic adjustments to divide
migrations. We explore the relevance and complementarity of tested relative
stability metrics between neighbouring basins. We then investigate the impact
of uplift rate, erodibility and hillslope processes on the dynamics of
divide migration and associated denudation rates. Finally, we apply our
approach to the basin-wide denudation rates dataset of Matmon (2003a, b) in
the case of the Great Smoky Mountains and propose new criteria to guide
future sampling strategies to assess basin-wide denudation rates from river
sands.

## 2.1 Landscape evolution model (LEM)

We use TTLEM (TopoToolbox Landscape Evolution Model) (Campforts et al., 2017), a landscape evolution model based on the MATLAB function library TopoToolbox 2 (Schwanghart and Scherler, 2014). This LEM uses a finite volume method (Campforts and Govers, 2015) to solve the following equation of mass conservation for rock or regolith subject to uplift and denudation:

where $\partial z/\partial t$ is the variation in elevation with time,
$(\partial z/\partial t{)}_{\mathrm{td}}$ is the change of elevation due to tectonic horizontal advection, *U* is the rock uplift rate, and *ρ*_{r}∕*ρ*_{s} is the density ratio between the bedrock and the regolith. We use a linear formulation of hillslope diffusion (Culling, 1963) limited by a critical slope *S*_{c}:

where *q*_{s} is the flux of soil or regolith material. When slope values exceed *S*_{c}, they are readjusted to the critical value by using a modified version of the excess topography algorithm (Blöthe et al., 2015). The diffusivity *D* gives the rate of soil or regolith material creep. Its magnitude ranges from 10^{−3} to 10^{−1} m^{2} yr^{−1} in natural settings and varies with soil thickness, lithology and
vegetation (Roering et al., 1999; Jungers et al., 2009; West et al., 2013;
Richardson et al., 2019). Hillslope diffusion is implemented in TTLEM using
an implicit scheme, which is unconditionally stable at large time steps
(Pelletier, 2008). A non-linear diffusion formulation (Perron, 2011) is also
implemented in TTLEM. However, we favoured the use of a linear diffusion with
a critical slope, which is more convenient for the time step used in our
simulations (5000 years) and the set of parameters considered (see Sect. 2.2). Due to the relatively coarse spatial resolution of our models
(90 m), any of these diffusion formulations generate negligible
topographic differences on the direct vicinity of crest lines (Roering et
al., 1999, Campforts et al., 2017) and do not affect our results (see Fig. S1 in the Supplement). Fluvial incision is calculated with a stream power law:

where *K* is the erodibility coefficient reflecting climate, hydraulic
roughness, sediment load and lithology. Its value ranges between
10^{−16} and 10^{0} m^{(1−2 m)} yr^{−1} (Kirby and Whipple, 2001; Harel et al., 2016). *A* is the upstream area. *x*_{Γ} is the along-stream distance from the outlet of the river. *m* and *n* are two parameters that are usually reported as a *m*∕*n* ratio ranging between 0.35 and 0.8. The river incision law is implemented in TTLEM using an explicit scheme based on a higher-order flux-limiting finite volume method (FVM) that is total variation diminishing (TVD-FVM) (see Campforts and Govers, 2015, and Campforts et al., 2017, for further details). Its main advantage is to eliminate numerical diffusion, which is present in most other schemes solving differential equations of river incision. This last point has a significant impact on the accuracy of basin-wide simulated denudation rates, making TTLEM a well-suited LEM for the purpose of this study.

## 2.2 Modelling approach and assumptions

### 2.2.1 Geometry and meshing

Since the computation is performed using a discretized land surface, smaller mesh sizes lead to detailed topography but lengthen the computation time and memory requirements. Hereinafter, we consider a reference square landscape model of 50 km side with a grid resolution of 90 m, which is a good compromise between computation time (3–5 h on a PC workstation) and the total amount of basins that can be studied (>1000). Our results are not affected when the grid resolution is 30 m nor when the model size is 100 km×100 km (see Fig. S2).

### 2.2.2 Boundary conditions

In order to isolate the effect of divide migrations on the variability of basin-wide denudation rates, we explore simple models with constant and spatially uniform uplift and precipitation rates, and we assume no horizontal advection, $(\partial z/\partial t{)}_{\mathrm{td}}=\mathrm{0}$. We use a Dirichlet boundary condition: simulation edges are not affected by uplift on a one pixel band to represent a stable base level for rivers. The model presents no initial topography, except for gaussian noise ranging between 0 and 50 m so as to initiate a random fluvial network.

### 2.2.3 Set of parameters

Firstly, we consider a reference model with parameters commonly used for
moderately active orogens: an uplift rate *U* of 0.1 mm yr^{−1}, a
diffusivity *D* of 10^{−2} m^{2} yr^{−1} (Roering et al., 1999), a
threshold slope *S*_{c} of 30^{∘} (Burbank et al., 1996; Montgomery and Brandon, 2002; Binnie et al., 2007), and a *m*∕*n* ratio of 0.5 with *m*=0.5 and *n*=1, an erodibility coefficient *K* of $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{5}}$ m^{(1−2 m)} yr^{−1}, a *ρ*_{r}∕*ρ*_{s} ratio of 1.3.

Secondly, all other parameters held constant, we investigate the specific
impact of uplift rate, erodibility and hillslope processes in other models
by varying *U*, *K*, *D* and *S*_{c} between 0.01 and 1 mm yr^{−1}, between $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{5}}$ m^{(1−2 m)} yr^{−1} and $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{6}}$ m^{(1−2 m)} yr^{−1}, between 10^{−3} and
10^{−1} m^{2} yr^{−1}, and between 20 and 40^{∘}, respectively.

In order to better constrain the variability of our results under similar conditions, we ran for each model five simulations using the same parameters but with different initial random topographies.

### 2.2.4 Timescale

The total duration of simulations is 10 Myr. The implicit scheme used to
simulate linear hillslope processes provides stable solutions regardless of
the time step. In contrast, the explicit scheme used to model fluvial
incision requires a time step that satisfies the Courant–Friedrich–Lewy
criterion. Hereinafter, we choose a time step Δ*t* of 5000 years for
hillslope diffusion. Our results are not affected when using a smaller Δ*t* (i.e. 1000 years) (see Fig. S1). Incision computation is nested in this time step and uses another time step that is automatically determined to assure model stability (Campforts et al., 2017).

## 2.3 Basin-wide denudation rates and aggressivity metrics

### 2.3.1 Basin-wide denudation rates

We derive basins from the synthetic DEMs (digital elevation models) using an
accumulation map computed with a single flow direction algorithm implemented
in TopoToolbox (Schwanghart and Scherler, 2014). Next, we calculate for each
basin the variation in average elevation over a time interval of 10 kyr. The
drainage network migrates during the simulation, so we only survey the
basins that keep the same outlet location during this time interval.
Furthermore, due to divide mobility, the geometry of watersheds can also
change. Hence, we measure the average difference in elevation inside the
basin perimeter after 10 kyr. Here we only assess the surface uplift *U*_{s} (England and Molnar, 1990). To approximate the denudation rates *E* for each basin, we sum the surface uplift *U*_{s} with the rock uplift rate *U* and divide the result by the time interval. By considering the relatively small period over which we integrate denudation (10 kyr), we then assume that these approximations have a negligible impact on the results. If the basin is in a topographic steady state, *U*_{s} is equal to zero and *E* is equal to the background uplift rate. Thus, a positive (negative) value of *U*_{s} traduce a deficit (an excess) of denudation. Calculated that way, *E* is sensitive to divide migration but also to transient features like knickpoints that migrate along the river network. In our simulations, knickpoints may develop due to (1) the dissection of the initial flat surface or (2) discrete drainage captures (see Sect. 3.1). We use the “knickpointfinder” algorithm implemented in TopoToolbox (Schwanghart and Scherler, 2014) to identify the affected basins.

### 2.3.2 From cross-divide metrics to basin-averaged aggressivity metrics

Most recent studies have focused on the relationship between drainage divide
mobility and headwater cross-divide contrast in either *χ*, gradient,
height or local relief values (e.g. Whipple et al., 2017; Forte and
Whipple, 2018). Here, in line with Willett et al. (2014; see the Supplement
therein) we focus on the specific influence of divide migration on
denudation rates at the scale of the entire stream basin. Our approach aims
to integrate cross-divide contrasts in drainage network properties along the
entire basin perimeter. We then obtain basin-averaged aggressivity metrics
that determine if a watershed is either growing or shrinking (Willett et al.,
2014).

First, we assess *χ*, local topographic gradient *G* and height *H* of
the drainage network at a reference drainage area *A*_{ref} (Fig. 2).
Ideally, *A*_{ref} must be equal to the area at which channelization occurs (Forte and Whipple, 2018). However, it is challenging to locate the accurate position of channel heads (Clubb et al., 2016). Hence, we use a constant value of *A*_{ref} set to 1 km^{2}. The parameter *χ* is an integral function of position along the channel network (Perron and
Royden, 2012) described by the equation:

where *A*(*x*) is the upstream drainage area at location *x*, *A*_{0} is
an arbitrary scaling area set to 1 km^{2}. The *m* over *n* ratio refers here to the reference concavity of an equilibrated river profile. Its value is set to 0.5 in accordance with the model parameters. For each independent drainage network, we integrate *χ* from the outlet *x*_{b}, located at the model boundary (<1 m high), to the channel heads. Local gradient is determined for each DEM pixel from its eight connected neighbours. Height is simply extracted from the DEM.

Then, we calculate the difference in metrics (Δ*χ*, Δ*G* and Δ*H*) across the segments of divide shared by two reference basins. Finally, the aggressivity metric is obtained by averaging these
cross-divide differences along the perimeter of each sampled basin (Fig. 2). This way, the sign of the aggressivity metric in a basin corresponds to
the difference in the averaged value of considered metric difference
(Δ*χ*, Δ*G* and Δ*H*) in this basin with respect to
its neighbours. This method has the advantage of weighting individual
divide segments by the number of pixels they contain and then providing a
robust assessment of the basin aggressivity. Aggressivity metrics based on *χ*, *G* and *H* are hereafter referred to as Δ*χ*_{av},
Δ*G*_{av} and Δ*H*_{av}, respectively. However, due to topology issues, some parts of the perimeter of the sampled basins may be not shared by two reference basins (Fig. 2). We quantify this incompleteness by assessing the ratio of documented pixels over the total amount of pixel along the basin perimeter. We refer to this ratio as the “confidence index” CI, assuming that a higher CI is associated with a more robust basin aggressivity assessment.

## 3.1 Evolution of reference model

A detailed analysis of the DEM suggests that during the initial phase, the flat initial surface (Fig. 3a) is progressively uplifted to form a plateau. At the same time the edges of this plateau are gradually regressively eroded by drainage networks that spread from the base level toward the centre of the model (Fig. 3b and c). This transient landscape is completely dissected after 2 Myr. From this time and until the end of the simulation, landscape changes are mainly due to competition between watersheds, resulting in continuous divide migrations with decreasing intensity as the model is moving toward a total topographic equilibrium (Fig. 3d–f; video no. 1 in the Supplement).

To define the time period of regional steady state, we measure the average elevation, the maximum elevation and the average denudation rate over the entire model for each time step (Fig. 4a). We identify two distinct stages during the evolution of our reference simulation. During the first million years, due to long wavelength topographic building, the calculated landscapes are far from steady state. This leads to a major increase in the mean elevation from ca. 25 to ca. 75 m. In a second stage, this trend reverses and the mean elevation decreases asymptotically toward ca. 60 m until the end of the simulation.

The evolution of the maximum elevation follows the same pattern but can be affected by temporal changes in the location and altitude of highest peaks. The maximum elevation increases between ca. 50 and ca. 250 m over the first 3 Myr (Fig. 4a) then decreases progressively to remain at ca. 200 m during the rest of the simulation.

We compute the average denudation rate from the rock uplift rate and from average elevation change over the entire model between two time steps:

where (Δ*z*∕Δ*t*)_{av} is the average surface uplift over the entire model on a time-step Δ*t*, *U* is the imposed uniform uplift rate (0.1 mm yr^{−1}) and *E*_{av} is the average “real” denudation rate. During the first 0.25 Myr, the mean denudation rate falls abruptly from ca. 0.6 mm yr^{−1} to nearly 0 mm yr^{−1} as a consequence of diffusion over the initial flat topography. After that time and until the first 1 Myr, the mean denudation rate increases but remains lower than the uplift rate, leading to the increase in average elevation over this time period. In the following 1 Myr, *E*_{av} exceeds the uplift rate to reach up to 0.104 mm yr^{−1} before it gently decreases to 0.1 mm yr^{−1} until the end of the simulation. This shows that topography tends to – but never reaches – a strict steady state over the simulation time. Abrupt changes in *E*_{av} after ca. 2.5, 3.5, 4, 5 and 9.5 Myr (red circles in Fig. 4b) are related to major local captures in the drainage network, which can be observed during the model evolution (red circles in Fig. 3e and f and video no. 1 in the Supplement).

Based on these results, we will consider that a regional topographic
steady state is reached between 1.5 and 2 Ma, when the plateau relict
topography is totally eroded and *E*_{av} begins to decrease (Figs. 3 and 4). This time is consistent with the time required to reach topographic
steady state proposed from models with constant uplift rate and no
horizontal advection (Willett et al., 2001).

## 3.2 Basin-wide denudation rates variability

We calculate basin-wide denudation rates *E* upstream of each stable drainage network confluence after 2.5, 5 and 10 Myr of simulation (Fig. 5a–c, respectively). Regardless of the duration, we observe a significant variability in the calculated denudation rates depending on basin size. As exposed by Forte and Whipple (2018), the erosion rate contrasts across divides are spatially limited to areas very near the divides. Thus, the variability is maximum for small basins (ca. 1 km^{2}) and decreases with increasing basin area. In our approach, small basins are nested in larger ones. Hence, these results can be related to the averaging of denudation rates along the drainage network, in agreement with the measurements of Matmon et al. (2003b). This variability also decreases with time (Fig. 5a–c). For basins with an excess of denudation relative to the uplift rate *U*, the *E*∕*U* ratio can reach up to 2.5 after 2.5 Myr but only 2 after 5 Myr and 1.7 after 10 Myr. Basins with a denudation excess that stand out of the general trend at 10 Ma (Fig. 5c) are associated with a capture event visible in Fig. 4b. For basins with a deficit of denudation, the evolution of the ratio is less obvious. It can be lower than 0.5 after 2.5 Myr, but it increases slightly to 0.6 until 10 Myr. These results reflect a significant spatial variability of the difference between basin-wide denudation rates and uplift rate. To assess more accurately the temporal evolution of this variability, we calculate *E* every 0.5 Myr for three distinct categories of basin sizes: 1–2, 10–20 and 100–200 km^{2}. We then estimate the mean absolute deviation (MAD) from the uplift rate by considering separately basins with a denudation in excess or in deficit of uplift rate (Fig. 5d). Until 1.5 Ma, basins are located on the plateau where denudation rate is null. This leads to a low MAD for basins with a denudation deficit and to the absence of basins with a denudation excess. After 1.5 Ma, basins in deficit exhibit an increase in MAD from nearly −0.15 to −0.04 mm yr^{−1}, regardless of the area class considered. For basins in excess, the MAD value decreases through time, depending on drainage area: from ca. 0.25 to ca. 0.07 mm yr^{−1} for basins with an area of 1–2 km^{2}, from
ca. 0.2 to ca. 0.07 mm yr^{−1} for basins of 10–20 km^{2} and from ca. 0.7 to ca. 0.04 mm yr^{−1} for the largest basins. We see a coherent evolution of this difference over the simulation time, consistent with the model progression toward topographic equilibrium.

The spatial variability of the denudation rates is neither homogeneous nor
randomly distributed (Fig. 6a). The location of drainage basins with
denudation rates far from the equilibrium value of 0.1 mm yr^{−1}
coincides with migrating drainage divides (Fig. 3d) and with cross-divide
contrasts in channel *χ*, gradient and height (Fig. 6b–d). Following
Willett et al. (2014) and Forte and Whipple (2018), the divide migrations
predicted by these contrasts are consistent with the direction of divide
mobility obtained from our model. One may note that the higher the contrast
in these parameters across the divide, the higher the deviation of the
denudation rate from the uplift rate, and therefore from topographic
equilibrium. None of the sampled basins in this dataset contain a
knickpoint. Thus, these results based on simulations assuming uniform and
constant properties as well as constant boundary conditions confirm that the
dispersion observed in denudation rates is primarily controlled by divide
migration. Basins that expand (shrink) show higher (lower) denudation rates
compared to uplift rate, and they are hereafter referred to as aggressors
(victims), following the terminology adopted by Willett et al. (2014).

## 3.3 Deviation of denudation rates from the uplift rate and basin aggressivity

Willett et al. (2014) showed that the basin-averaged cross-divide contrast
in *χ*, could be used to deduce an aggressivity metric for basins. We
extend this basin-scale approach to the Gilbert metrics recently proposed
by Forte and Whipple (2018) including cross-divide contrast in headwater
gradient and elevation.

We here assess the relationship between the *E*∕*U* ratio and these
aggressivity metrics. First, to exclude variability related to both basin
area and time, we focus on a single class of basins with a size of 2–4 km^{2} gathered from five computed reference models after a simulation duration of 2.5 Myr. Denudation rates may be affected by
knickpoints, which are a source of transient perturbation at the scale of
the catchment. Therefore, in order to focus only on perturbations associated
with drainage divide dynamics, basins that contain knickpoints are ignored.
In agreement with cross-divide metrics tested by Forte and Whipple (2018),
graphs in Fig. 7 must be divided into four quadrants. Aggressor (victim)
basins have negative (positive) Δ*χ*_{av} and Δ*H*_{av} values and conversely a positive (negative) Δ*G*_{av} value (Fig. 2). Theoretically, aggressor (victim) basins have higher (lower) denudation rates than the underlying uplift rate. This result is verified for ca. 81 %, 52 % and 81 % of basins
for Δ*χ*_{av}, Δ*G*_{av} and Δ*H*_{av}, respectively. For this limited dataset, the
evolution between *E*∕*U* and both Δ*χ*_{av} and Δ*H*_{av} may be defined by a linear relationship (Fig. 7b). Compared to other metrics, Δ*G*_{av} is less sensitive to drainage migration and shows a more scattered distribution.

In natural settings, the stage of evolution of landscapes cannot be easily
defined and the total amount of basins with a specific size may be limited.
The large dataset from our modelling can provide further insights by
gathering the results obtained every 0.5 Myr for seven classes of basin
areas expanding geometrically with a multiplying factor of 2 from 1–2 to
64–128 km^{2} (Fig. 7b). Basins that contain knickpoints are
discarded from the analysis. When all classes of drainage areas are combined
together, we still obtain a clear relationship between aggressivity metrics
and *E*∕*U*, with 77 %, 56 % and 78 % of basins lying in aggressor
or victim quadrants for Δ*χ*_{av}, Δ*G*_{av} and Δ*H*_{av}, respectively (Fig. 7b). Our results highlight the major control of basin size on the dispersion *E*∕*U*. Part of the variability intrinsic for each class of basin area may in turn be explained by heterogeneities in aggressivity between different parts of a basin. Figure 7c shows that this dispersion is related to the standard deviation of
aggressivity metrics, Δ*χ*_{std}, Δ*G*_{std}
and Δ*H*_{std}. In other words, basins where different divide segments migrate at different rates or in different directions are more scattered. The lower the confidence index, the more scattered the results (Fig. 7d). Thus, some dispersion may come from approximations due to undocumented divide segments performed when averaging metric differences between reference basins (Fig. 2). One may note two different trends for victim and aggressor basins. Aggressors show a more scattered distribution for Δ*χ*_{av} and Δ*G*_{av} metrics. When compared to victims, these basins have hillslopes closer to the critical value *S*_{c} (Fig. S3). Hence, the dispersion may be explained by the non-linear relationship existing between denudation rates and basin slope (Montgomery and Brandon, 2002; Binnie et al., 2007).

## 4.1 Sensitivity tests

The reference model involves various parameters related to uplift, fluvial
incision and hillslope denudation. A systematic analysis of trade-offs
between all parameters is out of the scope of this article. In this
section, we assess the sensitivity of the results to both tectonic and
erosion processes by studying the specific impact of uplift *U*,
erodibility *K*, diffusivity *D* and critical hillslope gradient *S*_{c} taken separately. Varying these parameters may change the simulation time required to erode the plateau associated with the initial boundary conditions. In this section, to reduce sensitivity dependence on these initial conditions, we only consider results obtained between 5 and 10 Ma.

### 4.1.1 Sensitivity to uplift rate

We test rock uplift rates of 0.01, 0.1 (hereafter called reference model) and 1 mm yr^{−1} to cover the range of a large variety of geodynamic settings (Champagnac et al., 2012). It is well known that a river responds to a fall in base level (due to changes in rock uplift rate or other forcing) by cutting downward into its bed, deepening and widening its active channel. In our simulations, changes in uplift rate lead to variations in the density of the drainage network. Compared to the reference model, an uplift rate of 1 mm yr^{−1} (0.01 mm yr^{−1}) results in a decrease (increase) of drainage density. These results are consistent with previous studies that show an inverse relationship between drainage density and erosion rates in equilibrium topography when using a threshold slope for diffusion processes (Tucker and Bras, 1998; Clubb et al., 2016). An increase in uplift rate favours river entrenchment leading to an increase in the range of Δ*G*_{av} and Δ*H*_{av} (Fig. 8). Hence, these two Gilbert metrics appear to be well suited to diagnose local disequilibrium for higher uplift rates (i.e. ≥1 mm yr^{−1}). Conversely an increase in uplift rate induces a lower range of values for Δ *χ*_{av}. This last observation is explained by the decrease in drainage density and associated stream length.

Maximum variability of *E*∕*U* reaches a factor of 2, regardless of the assumed uplift rate between 1 and 0.01 mm yr^{−1}. The observed small differences suggest that limited uplift rates promote diffusive processes (see Sect. 4.1.3).

### 4.1.2 Influence of erodibility

Fluvial erosion is proportional to the erodibility coefficient *K* that may
reflect, among others, rock strength and climate. We let this parameter vary
between $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{6}}$ and $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{5}}$ m^{(1−2 m)} yr^{−1}. As expected from (Eqs. 1 and 3), we find that erodibility and uplift rates have opposite effects. Lower (higher) values of erodibility lead to higher (lower) average topography. Thus, an increase (decrease) in erodibility decreases (increases) the range of all aggressivity metrics (Fig. 9). Lower values of erodibility also increase the range of the *E*∕*U* ratio. Models with higher (lower) erodibility reach a quasi-topographic steady state earlier (at a later stage). Hence, differences in the variability of *E*∕*U* may be related to different stages of evolution for each model over the period we consider (5 to 10 Ma) (Fig. 5d).

### 4.1.3 Influence of hillslope processes

Hillslope denudation is proportional to the diffusivity coefficient *D* and
depends on the critical slope *S*_{c} (Eq. 2). To test the effect of hillslope processes, we let *D* vary between 10^{−3} and 10^{−1} m^{2} yr^{−1}. Compared to the reference model, we find no differences in the case of a lower diffusivity (i.e. 0.001 m^{2} yr^{−1}) (Fig. 10c). In contrast, for models with higher diffusivity coefficient (i.e. 0.1 m^{2} yr^{−1}), this parameter has a significant effect on both the range of *E*∕*U* and the aggressivity metric Δ*G*_{av} (Fig. 10a). This result is consistent with the observations described in Sect. 4.1.1. It derives from a stronger impact of diffusive processes, which decrease local slopes in the vicinity of divides. In our modelling, the local slopes remain lower than the fixed critical value. Then, assuming a critical slope between 20 and 40^{∘}, we find that *S*_{c} does not affect significantly the relationship between the *E*∕*U* ratio and the studied metrics (Fig. 11).

Altogether, these sensitivity tests demonstrate the robustness of our findings. Regardless of the tested parameter values, we observe a relationship between aggressivity metrics and deviation of denudation rates from uplift rates. Thus, aggressivity metrics are, to the 1st-order, reliable metrics to assess the effect of divide mobility on basin-wide denudation rates inferred from simulations. In the following section, we apply this approach to field observations and discuss the consequences for sampling and interpretation.

## 4.2 Implications for the interpretation of basin-wide denudation rates

Over the last decades, measurements of cosmogenic radionuclide (CRN) concentrations in alluvial sediments (see Granger et al., 2013, and references therein), of suspended sediments (Gabet et al., 2008) and of detrital thermochronology (Huntington and Hodges, 2006) have become common practices to assess basin-wide denudation rates. However, their interpretation remains debated, even in settings where topographic steady state is supposedly achieved regionally.

### 4.2.1 Application to the Great Smoky Mountains

As previously mentioned (Matmon et al., 2003a, b), while the Great Smoky
Mountains in the southern Appalachians are expected to be in a
quasi-topographic steady state, basin-wide denudation rates show a strong
dispersion up to a factor of 2 in comparison to the estimated uplift rate
(ca. 0.03 mm yr^{−1}; see Fig. 1). We use the data associated with 40 basins originally sampled by Matmon et al. (2003a, b) and for which
denudation rates were recalculated by Portenga and Bierman (2011).
Following our method, we calculate the three basin-averaged aggressivity
metrics Δ*χ*_{av}, Δ*G*_{av} and Δ*H*_{av} associated with these 40 catchments (Fig. 12; see also Fig. S3). The calculation of *χ* requires us to define the elevation of the catchment outlets *H*_{b} and the *m*∕*n* ratio (Eq. 4). As underlined by Forte and Whipple (2018), the choice of the “correct” outlet elevation is non-trivial in natural settings. We first consider a local base level given by the Tennessee River. To test the relevance of this choice, we also use a base level located at a fixed arbitrary elevation *H*_{b}=400 m. We assume the same *m*∕*n* ratio value of 0.45 as used by Willett et al. (2014) for the Great Smoky Mountains. For all calculated metrics, the majority (ca. 58 % for Δ*G*_{av} and ca. 66 % for Δ*H*_{av}) of the basins is located in the expected quadrants (see Fig. 7). However, more attention must be given to the results based on Δ*χ*_{av}. For this metric, ca. 58 % of the analysed basins lie in the expected quadrant when we consider the Tennessee River as the local base level versus ca. 68 % for *H*_{b}=400 m (Fig. 12b). Although the overall results are similar, we show that the choice of a different base level *H*_{b} leads to significant variations in Δ*χ*_{av} for individual basins. This highlights the main weakness of the Δ*χ*_{av} metric, which is highly sensitive to the choice of the proper base level *H*_{b}. Nevertheless, our results confirm the findings by Willett et al. (2014), suggesting that a significant part of the data variance observed in the Matmon et al. (2003a, b) can be explained by divide migration (Fig. 12), raising this possible explanation for the variability of most natural datasets. One may note that the southern Appalachians exhibit migrating knickpoints that can locally affect denudation rates (Gallen et al., 2011, 2013). This last point can also explain part of the observed variability in this dataset but this specific impact is beyond the scope of the present study.

Based on both our simulations and this field dataset, we propose favouring
the use of Δ*χ*_{av} and Δ*H*_{av}. Among the tested metrics, Δ*G*_{av} appears the least sensitive to disequilibrium, excepted in active mountain belts with rock uplift *U*≥1 mm yr^{−1}.

### 4.2.2 Assessment of topographic disequilibrium

Topographic steady state is a very convenient assumption and concept to deduce the uplift pattern in mountains ranges from denudation rates, and thus to obtain significant information on the geometry of active structures and on orogen dynamics (Lavé and Avouac, 2001; Godard et al., 2014; Scherler et al., 2014; Le Roux-Mallouf et al., 2015). However, this assumption is seldom verified at the scale of sampled watersheds.

On the basis of our modelling, we show that the competition between low-order basins has a significant impact on basin-wide denudation rates. The proposed approach provides a new tool to assess the potential deviation from topographic steady state based on aggressivity metrics and drainage area, which can both be inferred from a simple DEM: the closer to zero the aggressivity metrics and the lower the standard deviation of cross-divide metrics, the more representative of uplift rate the measured denudation rates.

### 4.2.3 Improvement of sampling strategy

Basin-wide denudation rates obtained from CRN concentration measurements, suspended sediments or detrital thermochronology depend on many parameters including lithology, ice cover, rainfall, landslide activity or tectonic uplift (Vance et al., 2003; Bierman and Nichols, 2004; Wittmann et al., 2007; Yanites et al., 2009; Norton et al., 2010; Godard et al., 2012; Whipp and Elhers, 2019). Hence, to unravel the influence of tectonics from other processes, a specific sampling strategy is usually recommended: (1) to sample catchments with homogeneous lithologies to limit the effect of spatial variations in the abundance of target minerals in bedrock formations; (2) to select catchments with no ice cover (past or present) because the input of glacier-derived sediments can significantly complicate the interpretation of CRN concentrations; (3) to choose areas with spatially uniform rainfall distribution; and (4) to consider watersheds where the relative contribution of landslides to long-term landscape evolution is low. Unfortunately, these different criteria imply selection of watersheds with variable sizes. The first three criteria favour the sampling of small catchments, whereas the last one requires basins large enough to be less affected by landslides.

Our approach suggests the need to pre-assess targeted basins for their potential divide mobility before sampling for CRN concentration measurements. If the objective is to quantify the background uplift rate, one should sample basins that satisfy the conditions we previously described in the current section and also display an aggressivity close to zero and with the smallest associated standard deviation. Conversely, to quantify the specific denudation rate associated with the migration of drainage divides, small aggressor or victim basins should be favoured.

Based on our simulations, a relationship between the maximum of erosion
variability (0.5 and 99.5 percentiles, respectively) due to divide mobility $\left[\right(E-U)/U{]}_{\mathrm{max}}$ and the catchment size *A* can be derived (Fig. 13). Our results suggest a logarithm dependence between these two parameters, regardless of the assumed *U*, *K*, *D* and *S*_{c}:

with *c*_{1} and *c*_{2} being two parameters that depend on the balance between
erosion processes, uplift rate and state of evolution of the landscape.

Calculations from a landscape evolution model assuming spatially uniform uplift, rock strength and rainfall confirm that the concept of topographic steady state is relevant at the scale of entire mountain belts, but this represents an oversimplification at the scale of individual watersheds. Our simulations underline the role of divide mobility on deviations from equilibrium, which can lead to significant differences between tectonic uplift rate and basin-wide denudation rates even if an overall topographic steady state is achieved at large scale.

To better assess these deviations, we propose new basin-averaged
aggressivity metrics – Δ*χ*_{av}, Δ*G*_{av} and Δ*H*_{av} – based on the approach by Willett et al. (2014) and Forte and Whipple (2018). They include mean cross-divide contrasts in channel *χ*, local gradient and height. From our calculations, Δ*χ*_{av} is the most reliable aggressivity metric to assess local disequilibrium, but it is highly dependent on the chosen base level, which remains hard to constrain. Gilbert metrics Δ*G*_{av} and Δ*H*_{av} are more suitable for relatively high uplift rate (i.e. ≥1 mm yr^{−1}). Altogether, our metrics reveal that deviation of denudation rates from uplift rate related to divide migrations depends on both basin aggressivity and basin area. This last parameter has a key control on the dispersion in *E*∕*U*, which can reach a factor of 2, regardless of the imposed uplift rate (here 0.01–1 mm yr^{−1}), erodibility (here $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{6}}$–$\mathrm{5}\times {\mathrm{10}}^{-\mathrm{5}}$ m^{(1−2 m)} yr^{−1}), diffusivity (here 10^{−3}–10^{−1} m^{2} yr^{−1}) or hillslope gradient (here 20–40^{∘}). By comparing our results to CRN measurements from the Great Smoky Mountains (Matmon et al., 2003a, b), we show that this approach can be used to improve field sampling strategies and provides a new tool to derive a minimal uncertainty in basin-wide denudation rates due to topographic disequilibrium.

For the sake of simplicity our models involve spatially homogenous and time
invariant parameters. Additional simulations are now needed to test this
approach in more complex settings, including spatial and temporal
variability in climate and tectonic forcing or parameters like stream power
equation exponents *n* and *m*.

The data that support the findings of this study are available from the corresponding author on request.

The supplement related to this article is available online at: https://doi.org/10.5194/esurf-7-1041-2019-supplement.

RC and MF initiated this study. TSS performed the simulations and topographic analyses. All authors contributed to the writing of the paper.

The authors declare that they have no conflict of interest.

We are greatly indebted to referees Fiona Clubb and Adam Forte for providing constructive reviews that significantly improved the quality of the article. We thank Wolfgang Schwanghart and Benjamin Campforts for providing the TopoToolbox and TTLEM codes to analyse and to simulate landscape evolution. Through the work of Martine Simoes, this study contributes to the IdEx Université de Paris ANR-18-IDEX-0001 and is IPGP contribution no. 4082.

Timothée Sassolas-Serrayet's PhD is supported by a fellowship from the French Ministry for Higher Education. This research has been supported by the Agence Nationale de la Recherche (project ANR-18-CE01-0017 (Topo-Extreme)).

This paper was edited by Simon Mudd and reviewed by Adam Forte and Fiona Clubb.

Bierman, P. and Nichols, K. K.: Rock to sediment slope to sea with ^{10}Be rates of landscape change, Annu. Rev. Earth Planet. Sci., 32,
215–255, 2004.

Binnie, S. A., Phillips, W. M., Summerfield, M. A., and Fifield, L. K.: Tectonic uplift, threshold hillslopes, and denudation rates in a developing mountain range, Geology, 35, 743–746, 2007.

Blöthe, J. H., Korup, O., and Schwanghart, W.: Large landslides lie low: Excess topography in the Himalaya-Karakoram ranges, Geology, 43, 523–526, https://doi.org/10.1130/G36527.1, 2015.

Burbank, D. W., Leland, J., Fielding, E., Anderson, R. S., Brozovic, N., Reid, M. R., and Duncan, C.: Bedrock incision, rock uplift and threshold hillslopes in the northwestern Himalayas, Nature, 379, 505–510, https://doi.org/10.1038/379505a0, 1996.

Campforts, B. and Govers, G.: Keeping the edge: A numerical method that avoids knickpoint smearing when solving the stream power law, J. Geophys. Res.-Ea. Surf., 120, 1189–1205, 2015.

Campforts, B., Schwanghart, W., and Govers, G.: Accurate simulation of transient landscape evolution by eliminating numerical diffusion: the TTLEM 1.0 model, Earth Surf. Dynam., 5, 47–66, https://doi.org/10.5194/esurf-5-47-2017, 2017.

Champagnac, J. D., Molnar, P., Sue, C., and Herman, F.: Tectonics, climate and mountain topography, J. Geophys. Res., 117, B02403, https://doi.org/10.1029/2011JB008348, 2012.

Clubb, F. J., Mudd, S. M., Attal, M., Milodowski, D. T., and Grieve, S. W. D.: The relationship between drainage density, erosion rate, and hilltop curvature: Implications for sediment transport processes, J. Geophys. Res.-Ea. Surf., 121, 1724–1745, https://doi.org/10.1002/2015JF003747, 2016.

Culling, W. E. H.: Soil Creep and the Development of HillsideSlopes, J. Geol., 71, 127–161, https://doi.org/10.1086/626891, 1963.

Dahlquist, M. P., West, A. J., and Li, G.: Landslide-driven drainage divide migration, Geology, 46, 403–406, https://doi.org/10.1130/g39916.1, 2018.

England, P. and Molnar, P.: Surface uplift, uplift of rocks, and exhumation of rocks, Geology, 18, 1173–1177, 1990.

Forte, A. M. and Whipple, K. X.: Criteria and tools for determining drainage divide stability, Earth Planet. Sc. Lett., 493, 102–117, 2018.

Gabet, E. J., Burbank, D. W., Pratt-Sitaula, B., Putkonen, J., and Bookhagen, B.: Modern erosion rates in the High Himalayas of Nepal, Earth Planet. Sc. Lett., 267, 482–494, https://doi.org/10.1016/j.epsl.2007.11.059, 2008.

Gallen, S. F., Wegmann, K. W., Frankel, K. L., Hughes, S., Lewis, R. Q., Lyons, N., Paris, P., Ross, K., Bauer, J. B., and Witt, A. C.: Hillslope response to knickpoint migration in the Southern Appalachians: implications for the evolution of post-orogenic landscapes, Earth Surf. Proc. Land., 36, 1254–1267, 2011.

Gallen, S. F., Wegmann, K. W., and Bohnenstiehl, D. R.: Miocene rejuvenation of topographic relief in the southern Appalachians, GSA Today, 23, 4–10, https://doi.org/10.1130/GSATG163A.1, 2013.

Gilbert G. K.: Report on the Geology of the Henry Mountains, USGS Report, US Government Printing Office, Washington, D.C., 1877.

Godard, V., Lavé, J., and Cattin, R.: Numerical modelling of erosion processes in the Himalayas of Nepal: Effects of spatial variations of rock strength and precipitation, Geol. Soc. Lond. Spec. Publ., 253, 341–358, 2006.

Godard, V., Burbank, D. W., Bourlès, D. L., Bookhagen, B., Braucher, R.
and Fisher, G. B.: Impact of glacial erosion on ^{10}Be concentrations in fluvial sediments of the Marsyandi catchment, central Nepal, J. Geophys. Res.-Ea. Surf., 117, F03013, https://doi.org/10.1029/2011JF002230, 2012.

Godard, V., Bourlès, D. L., Spinabella, F., Burbank, D. W., Bookhagen, B., Fisher, G. B., Moulin, A., and Léanni, L.: Dominance of tectonics over climate in Himalayan denudation, Geology, 42, 243–246, 2014.

Granger, D. E., Lifton, N. A., and Willenbring, J. K.: A cosmic trip: 25 years of cosmogenic nuclides in geology, GSA Bull., 125, 1379–1402, 2013.

Hack, J. T.: Interpretation of erosional topography in humid temperate regions, Am. J. Sci., 258-A, 80–97, 1960.

Harel, M.-A., Mudd, S. M., and Attal, M.: Global analysis of the stream power
law parameters based on worldwide ^{10}Be denudation rates, Geomorphology, 268, 184–196, https://doi.org/10.1016/j.geomorph.2016.05.035, 2016.

Hasbargen, L. E. and Paola, C.: Landscape instability in an experimental drainage basin, Geology, 28, 1067–1070, 2000.

Hasbargen, L. E. and Paola, C.: How predictable is local erosion rates in eroding landscapes?, Predicti. Geomorphol., 135, 231–240, https://doi.org/10.1029/135GM16, 2003.

Huntington, K. W., and Hodges, K. V.: A comparative study of detrital mineral and bedrock age-elevation methods for estimating erosion rates, J. Geophys. Res.-Ea. Surf., 111, 1–11, https://doi.org/10.1029/2005JF000454, 2006.

Jungers, M. C., Bierman, P. R., Matmon, A., Nichols, K., Larsen, J., and
Finkel, R.: Tracing hillslope sediment production and transport with in situ
and meteoric ^{10}Be, J. Geophys. Res.-Ea. Surf., 114, F04020, https://doi.org/10.1029/2008JF001086, 2009.

Kirby, E. and Whipple, K.: Quantifying differential rock-uplift rates via stream profile analysis, Geology, 29, 415–418, 2001.

Lavé, J. and Avouac, J.-P.: Fluvial incision and tectonic uplift across the Himalayas of central Nepal, J. Geophys. Res., 106, 26561–26591, 2001.

Le Roux-Mallouf, R., Godard, V., Cattin, R., Ferry, M., Gyeltshen, J., Ritz, J. F., Drupka, D., Guillou, V., Arnold, M., Aumaître, G., and Bourles, D. L.: Evidence for a wide and gently dipping Main Himalayan Thrust in western Bhutan, Geophys. Res. Lett., 42, 3257–3265, 2015.

Matmon, A., Bierman, P. R., Larsen, J., Southworth, S. Pavich, M., Kinkel, R., and Caffee, M.: Erosion of an ancient mountain range, the Great Smoky Mountains, North Carolina and Tennessee, Am. J. Sci., 303, 817–855, 2003a.

Matmon, A., Bierman, P. R., Larsen, J., Southworth, S., Pavich, M., and Caffee, M.: Temporally and spatially uniform rates of erosion in the southern Appalachian Great Smoky Mountains, Geology, 31, 155–158, 2003b.

Miller, S. R., Sak, P. B., Kirby, E., and Bierman, P. R.: Neogene rejuvenation of central Appalachian topography: Evidence for differencial rock uplift from stream profiles and erosion rates, Earth Planet. Sc. Lett., 369–370, 1–12, https://doi.org/10.1016/j.epsl.2013.04.007, 2013.

Montgomery, D. R.: Slope distributions, threshold hillslopes, and steady-state topography, Am. J. Sci., 301, 432–454, 2001.

Montgomery, D. R. and Brandon, M. T.: Topographic controls on erosion rates in tectonically active mountain ranges, Earth Planet. Sc. Lett., 201, 481–489, 2002.

Norton, K. P., von Blanckenburg, F., and Kubik, P. W.: Cosmogenic nuclide-derived rates of diffusive and episodic erosion in the glacially sculpted upper Rhone Valley, Swiss Alps, Earth Surf. Proc. Land., 35, 651–662, 2010.

Pelletier, J. D.: Persistent drainage migration in a numerical landscape evolution model, Geophys. Res. Lett., 31, L20501, https://doi.org/10.1029/2004GL020802, 2004.

Pelletier, J. D.: Quantitative Modeling of Earth Surface Processes, Cambridge University Press, Cambridge, UK, 2008.

Perne, M., Covington, M. D., Thaler, E. A., and Myre, J. M.: Steady state, erosional continuity, and the topography of landscapes developed in layered rocks, Earth Surf. Dynam., 5, 85–100, https://doi.org/10.5194/esurf-5-85-2017, 2017.

Perron, J. T.: Numerical methods for nonlinear hillslope transport laws, J. Geophys. Res.-Ea. Surf. 116, F02021, https://doi.org/10.1029/2010JF001801, 2011.

Perron, J. T. and Royden, L.: An integral approach to bedrock river profile analysis, Earth Surf. Proc. Land., 38, 570–576, 2012.

Portenga, E. W. and Bierman, P. R.: Understanding Earth's eroding surface with ^{10}Be, GSA Today, 21, 4–10, 2011.

Reiners, P. W., Ehlers, T. A., Mitchell, S. G., and Montgomery, D. R.: Coupled spatial variations in precipitation and long-term erosion rates across the Washington Cascades, Nature, 426, 645–647, 2003.

Reinhardt, L. and Ellis, M. A.: The emergence of topographic steady state in a perpetually dynamic self-organized critical landscape, Water Resour. Res., 51, 4986–5003, 2015.

Richardson, P. W., Perron, J. T., and Schurr, N. D.: Influences of climate and life on hillslope sediment transport, Geology, 47, 1–4, https://doi.org/10.1130/G45305.1, 2019.

Roering, J. J., Kirchner, J. W., and Dietrich, W. E.: Evidence for nonlinear, diffusive sediment transport on hillslopes and implications for landscape morphology, Water Resour. Res., 35, 853–870, 1999.

Royden, L. and Perron, J. T.: Solutions of the stream power equation and application to the evolution of river longitudinal profiles, J. Geophys. Res.-Ea. Surf., 118, 497–518, 2013.

Scherler, D., Bookhagen, B., and Strecker, M. R.: Tectonic control on
^{10}Be-derived erosion rates in the Garhwal Himalaya, India, J.
Geophys. Res.-Ea. Surf., 119, 83–105, 2014.

Schwanghart, W. and Scherler, D.: Short Communication: TopoToolbox 2 – MATLAB-based software for topographic analysis and modeling in Earth surface sciences, Earth Surf. Dynam., 2, 1–7, https://doi.org/10.5194/esurf-2-1-2014, 2014.

Stolar, D. B., Willett, S. D., and Montgomery, D. R.: Characterization of topographic steady state in Taiwan, Earth Planet. Sc. Lett., 261, 421–431, 2007.

Tucker, G. E. and Bras, R. L.: Hillslope processes, drainage density, and landscape morphology, Water Resour. Res., 34, 2751–2764, 1998.

Vance, D., Bickle, M., Ivy-Ochs, S., and Kubik, P. W.: Erosion and exhumation in the Himalaya from cosmogenic isotope inventories of river sediment, Earth Planet. Sc. Lett., 206, 273–288, 2003.

West, N., Kirby, E., Bierman, P., Slingerland, R., Ma, L., Rood, D., and
Brantley, S.: Regolith production and transport at the Susquehanna Shale Hills Critical Zone Observatory, part 2: insights from meteoric ^{10}Be, J. Geophys. Res.-Ea. Surf., 118, 1877–1896, 2013.

Whipp, D. M. and Elhers T. A.: Quantifying landslide frequency and sediment residence time in Nepal Himalaya, Sci. Adv., 5, eaav3482, https://doi.org/10.1126/sciadv.aav3482, 2019.

Whipple, K. X.: Fluvial landscape response time: How plausible is steady-state denudation?, Am. J. Sci., 301, 313–325, 2001.

Whipple, K. X., Forte, A. M., DiBiase, R. A., Gasparini, N. M., and Ouimet, W. B.: Timescales of landscape response to divide migration and drainage capture: Implications for the role of divide mobility in landscape evolution, J. Geophys. Res.-Ea. Surf., 122, 248–273, 2017.

Willett, S. D., Slingerland, R., and Hovius, N.: Uplift, shortening, and steady state topography in active mountain belts, Am. J. Sci., 301, 455–485, 2001.

Willett, S. D., McCoy, S. W., Perron, J. T., Goren, L., and Chen, C. Y.: Dynamic reorganization of river basins, Science, 343, 1248765, https://doi.org/10.1126/science.1248765, 2014.

Wittmann, H., von Blanckenburg, F., Kruesmann, T., Norton, K. P., and Kubik, P. W.: Relation between rock uplift and denudation from cosmogenic nuclides in river sediment in the Central Alps of Switzerland, J. Geophys. Res.-Ea. Surf., 112, F04010, https://doi.org/10.1029/2006JF000729, 2007.

Yanites, B. J., Tucker, G. E., and Anderson, R. S.: Numerical and analytical models of cosmogenic radionuclide dynamics in landslide-dominated drainage basins, J. Geophys. Res.-Ea. Surf., 114, F01007, https://doi.org/10.1029/2008JF001088, 2009.