We present a statistical technique for analyzing longitudinal channel profiles. Our technique is based on the integral approach to channel analysis: Drainage area is integrated over flow distance to produce a transformed coordinate, χ, which has dimensions of length. Assuming that profile geometry is conditioned by the stream power law, defined as E = KA^mSⁿ where E is erosion rate, K is erodibility, A is drainage area, S is channel gradient, and m and n are constants, the slope of a transformed profile in χ‐elevation space should reflect the ratio of erosion rate to channel erodibility raised to a power 1/n; this quantity is often referred to as the channel steepness and represents channel slope normalized for drainage area. Our technique tests all possible contiguous segments in the channel network to identify the most likely locations where channel steepness changes and also identifies the most likely m/n ratio. The technique identifies locations where either erodibility or erosion rates are most likely to be changing. Tests on a simulated landscape demonstrate that the technique can accurately retrieve both the m/n ratio and the correct number and location of segments eroding at different rates where model assumptions apply. Tests on natural landscapes illustrate how the method can distinguish between spurious channel convexities due to incorrect selection of the m/n ratio from those which are candidates for changing erodibility or erosion rates. We also show how, given erosion or uplift rate constraints, the method can be used to constrain the slope exponent, n.

The concavity index, , describes how quickly river channel gradient declines downstream. It is used in calculations of normalized channel steepness index, , a metric for comparing the relative steepness of channels with different drainage area. It is also used in calculating a transformed longitudinal coordinate, , which has been employed to search for migrating drainage divides. A value of 0.45 is typically assumed in studies. Here we quantify the variability in across multiple landscapes distributed across the globe. We describe the degree to which both the spatial distribution and magnitude of and can be distorted if is assumed rather than constrained. Differences between constrained and assumed of 0.1 or less are unlikely to affect the spatial distribution and relative magnitude of values, but larger differences can change the spatial distribution of and in extreme cases invert differences in relative steepness: relatively steep reaches can appear relatively gentle as quantified by . These inversions are function of the range of drainage area in the considered watersheds. We also demonstrate that the coordinate, and therefore the detection of migrating drainage divides, is sensitive to varying values of . The median of most likely across a wide range of mountainous and upland environments is 0.425. This wide range of variability suggests workers should not assume any value for , but should instead calculate a representative for the landscape of interest, and exclude basins for which this value is a poor fit.

Drainage density is a fundamental landscape metric describing the extent of the fluvial network. We compare the relationship between drainage density (D_d) and erosion rate (E) using the Channel-Hillslope Integrated Landscape Development (CHILD) numerical model. We find that varying the channel slope exponent (n) in detachment-limited fluvial incision models controls the relationship between D_d and E, with n > 1 resulting in increasing D_d with E if all other parameters are held constant. This result is consistent when modeling both linear and nonlinear hillslope sediment flux. We also test the relationship between D_d and E in five soil-mantled landscapes throughout the USA: Feather River, CA; San Gabriel Mountains, CA; Boulder Creek, CO; Guadalupe Mountains, NM; and Bitterroot National Forest, ID. For two of these field sites we compare D_d to cosmogenic radionuclide (CRN)-derived erosion rates, and for each site we use mean hilltop curvature as a proxy for erosion rate where CRN-derived erosion rates are not available. We find that there is a significant positive relationship between D_d, E, and hilltop curvature across every site, with the exception of the San Gabriel Mountains, CA. This relationship is consistent with an n exponent greater than 1, suggesting that at higher erosion rates, the transition between advective and diffusive processes occurs at smaller contributing areas in soil-mantled landscapes.

We report a new program for calculating catchment-averaged denudation rates from cosmogenic nuclide concentrations. The method (Catchment-Averaged denudatIon Rates from cosmogenic Nuclides: CAIRN) bundles previously reported production scaling and topographic shielding algorithms. In addition, it calculates production and shielding on a pixel-by-pixel basis. We explore the effect of sampling frequency across both azimuth (Δθ) and altitude (Δϕ) angles for topographic shielding and show that in high relief terrain a relatively high sampling frequency is required, with a good balance achieved between accuracy and computational expense at Δθ = 8° and Δϕ = 5°. CAIRN includes both internal and external uncertainty analysis, and is packaged in freely available software in order to facilitate easily reproducible denudation rate estimates. CAIRN calculates denudation rates but also automates catchment averaging of shielding and production, and thus can be used to provide reproducible input parameters for the CRONUS family of online calculators.

Hillslope length is a fundamental attribute of landscapes, intrinsically linked to drainage density, landslide hazard, biogeochemical cycling and hillslope sediment transport. Existing methods to estimate catchment average hillslope lengths include inversion of drainage density or identification of a break in slope–area scaling, where the hillslope domain transitions into the fluvial domain. Here we implement a technique which models flow from point sources on hilltops across pixels in a digital elevation model (DEM), based on flow directions calculated using pixel aspect, until reaching the channel network, defined using recently developed channel extraction algorithms. Through comparisons between these measurement techniques, we show that estimating hillslope length from plots of topographic slope versus drainage area, or by inverting measures of drainage density, systematically underestimates hillslope length. In addition, hillslope lengths estimated by slope–area scaling breaks show large variations between catchments of similar morphology and area. We then use hillslope length–relief structure of landscapes to explore nature of sediment flux operating on a landscape. Distinct topographic forms are predicted for end‐member sediment flux laws which constrain sediment transport on hillslopes as being linearly or nonlinearly dependent on hillslope gradient. Because our method extracts hillslope profiles originating from every ridgetop pixel in a DEM, we show that the resulting population of hillslope length–relief measurements can be used to differentiate between linear and nonlinear sediment transport laws in soil mantled landscapes. We find that across a broad range of sites across the continental United States, topography is consistent with a sediment flux law in which transport is nonlinearly proportional to topographic gradient. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.