Journal cover Journal topic
Earth Surface Dynamics An interactive open-access journal of the European Geosciences Union
Journal topic
Earth Surf. Dynam., 7, 393-410, 2019
https://doi.org/10.5194/esurf-7-393-2019
Earth Surf. Dynam., 7, 393-410, 2019
https://doi.org/10.5194/esurf-7-393-2019

Research article 02 May 2019

Research article | 02 May 2019

# Cosmogenic 10Be in river sediment: where grain size matters and why

Cosmogenic 10Be in river sediment
Renee van Dongen1, Dirk Scherler1,2, Hella Wittmann1, and Friedhelm von Blanckenburg1,2 Renee van Dongen et al.
• 1GFZ German Research Centre for Geosciences, Earth Surface Geochemistry, Telegrafenberg, 14473 Potsdam, Germany
• 2Freie Universität Berlin, Institute for Geological Sciences, 12249 Berlin, Germany
Abstract

Concentrations of in-situ-produced cosmogenic 10Be in river sediment are widely used to estimate catchment-average denudation rates. Typically, the 10Be concentrations are measured in the sand fraction of river sediment. However, the grain size of bedload sediment in most bedrock rivers covers a much wider range. Where 10Be concentrations depend on grain size, denudation rate estimates based on the sand fraction alone are potentially biased. To date, knowledge about catchment attributes that may induce grain-size-dependent 10Be concentrations is incomplete or has only been investigated in modelling studies. Here we present an empirical study on the occurrence of grain-size-dependent 10Be concentrations and the potential controls of hillslope angle, precipitation, lithology, and abrasion. We first conducted a study focusing on the sole effect of precipitation in four granitic catchments located on a climate gradient in the Chilean Coastal Cordillera. We found that observed grain size dependencies of 10Be concentrations in the most-arid and most-humid catchments could be explained by the effect of precipitation on both the scouring depth of erosion processes and the depth of the mixed soil layer. Analysis of a global dataset of published 10Be concentrations in different grain sizes (n=73 catchments) – comprising catchments with contrasting hillslope angles, climate, lithology, and catchment size – revealed a similar pattern. Lower 10Be concentrations in coarse grains (defined as “negative grain size dependency”) emerge frequently in catchments which likely have thin soil and where deep-seated erosion processes (e.g. landslides) excavate grains over a larger depth interval. These catchments include steep (> 25) and humid catchments (> 2000 mm yr−1). Furthermore, we found that an additional cause of negative grain size dependencies may emerge in large catchments with weak lithologies and long sediment travel distances (> 2300–7000 m, depending on lithology) where abrasion may lead to a grain size distribution that is not representative for the entire catchment. The results of this study can be used to evaluate whether catchment-average denudation rates are likely to be biased in particular catchments.

1 Introduction

Catchment-average denudation rates are commonly estimated with in-situ-produced cosmogenic 10Be concentrations in river sediment (Bierman and Steig, 1996; Brown et al., 1995; Granger et al., 1996). 10Be is a rare isotope that is produced within quartz minerals by high-energy cosmic rays in the upper few metres of the Earth's surface (Gosse and Phillips, 2001). Its concentration records the time minerals were exposed to cosmic radiation, which is inversely proportional to denudation rates over timescales of 102–105 years (Bierman and Steig, 1996; Brown et al., 1995; Granger et al., 1996; Lal, 1991). Most studies use a sand fraction (0.1–2 mm) of river bedload sediment to estimate catchment-average denudation rates. However, bedload grain sizes found in bedrock rivers, where this method is frequently applied, are often much coarser (Fig. 1). The sand fraction provides a representative catchment-average denudation rate only if it is spatially and temporally representative for all erosion sources within the catchment (Bierman and Steig, 1996; von Blanckenburg, 2005; Brown et al., 1995; Gonzalez et al., 2017; Granger et al., 1996; Neilson et al., 2017; Willenbring et al., 2013). Evaluating this condition is challenging and requires a detailed understanding of the catchment and its erosion processes. Where 10Be concentrations differ amongst grain size fractions, using a non-representative grain size fraction could bias catchment-average denudation rates by a factor of 3 or more (Lukens et al., 2016). To date, there is no general consensus of what causes grain-size-dependent 10Be concentrations in a catchment. Some studies inferred that lower 10Be concentrations in coarse grains are caused by deep-seated erosion processes, such as landslides, which excavate material from greater depth where 10Be concentrations are lower (e.g. Aguilar et al., 2014; Belmont et al., 2007; Binnie et al., 2007; Brown et al., 1995; Puchol et al., 2014; Sosa Gonzalez et al., 2016a, b; Tofelde et al., 2018; West et al., 2014). In a recent study Tofelde et al. (2018) combined a detailed inventory of hillslope processes in a large semi-arid Andean catchment, with 10Be concentrations measured in the sand and gravel fraction of river sediment. They explained lower 10Be concentrations in the gravel compared to the sand by the scouring depth of erosion processes. However, other studies found that grain size reduction by abrasion during fluvial transport, or spatial variations in the provenance of different grain sizes can additionally account for grain-size-dependent 10Be concentrations (Carretier et al., 2009; Carretier and Regard, 2011; Lukens et al., 2016; Lupker et al., 2017; Matmon et al., 2003). This is particularly true for large and high-relief catchments that cover large elevation ranges because 10Be production rates depend on atmospheric depth. Carretier et al. (2015) conducted a comprehensive study on 10Be concentrations in different grain sizes on a precipitation gradient, sampling large Andean catchments. Despite significant contrasts in precipitation and, presumably, weathering and erosion processes, no systematic grain size dependency of 10Be concentrations controlled by precipitation emerged. A reason may be that as catchment size increases, more complexity in controlling factors is added that may blur potential trends in the data (Carretier et al., 2015; Portenga and Bierman, 2011). Hence, which types of catchments are sensitive to grain-size-dependent 10Be concentrations and biased catchment-average denudation rates remain elusive (e.g. Carretier et al., 2015), or they have only been addressed in modelling studies (e.g. Lukens et al., 2016).

Figure 1Grain size (GS) distributions of bedload sediment in rivers and grain sizes used for 10Be-derived catchment-average denudation rates (10Be GS, n=2735) (Codilean et al., 2018). Bedload grain size distributions were measured by pebble counts in bedrock rivers in Southern Italy and Sicily (Italy, n=3900) (Allen et al., 2015; Roda-Boluda et al., 2018) and the San Gabriel Mountains (SGM, n=5930) (DiBiase and Whipple, 2011; Scherler et al., 2016). Wolman pebble count fractions classified as < 2 mm are shown as 1 mm in the figure. Dashed line indicates 2 mm.

This paper presents the results of an empirical study in which we investigated the occurrence and cause of grain-size-dependent 10Be concentrations in river sediment. Our study consists of two parts. In the first part, we investigated the sole effect of precipitation in small granitic catchments in the Chilean Coastal Cordillera that differ mainly by mean annual precipitation. In the second part, we compiled and investigated a global dataset with previously published grain-size-dependent 10Be concentrations to include more catchment attributes in our analysis. In the following, we first provide a review of processes that control the grain size distribution and 10Be concentrations of river sediment to determine relevant catchment attributes for our analysis.

2 Why 10Be concentrations in river sediment can depend on grain size

## 2.1 Processes that control grain size

The grain size distribution of river sediment is a function of (1) weathering, (2) erosion processes at the hillslope, and (3) fluvial processes that change the grain size distribution. Chemical weathering at the hillslope converts bedrock into sediment of different grain sizes at a rate that is controlled by the properties of the parent rock, the climatic regime, and denudation rates at the surface (Sklar et al., 2017). The parent rock mineralogy sets rock dissolution rates and constrains the minimum size of individual minerals. Bedrock fractures provide water pathways and expose fresh bedrock to weathering (Lebedeva and Brantley, 2017; Oberlander, 1972; Ruxton and Berry, 1957). As chemical weathering rates are set primarily by the flux of water flowing through the regolith and, to some degree, also by temperature (Lasaga et al., 1994; Maher, 2010; White et al., 1999), there is a strong dependency of weathering and sediment production on the climatic regime (Dixon et al., 2016; Riebe et al., 2004). The presence of biota in humid climates can enhance the breakdown of rock fragments because microbes play an important role in chemical weathering processes, and growing roots can fracture bedrock (Drever, 1994; Ehrlich, 1998; Gabet and Mudd, 2010; Roering et al., 2010). The size reduction of bedrock fragments in the regolith depends, besides the chemical weathering rate, also on the time they spend in the regolith layer. This regolith residence time is controlled by the thickness of the regolith layer and the denudation rate, i.e. the rate of sediment removal from the surface (Anderson et al., 2007; Attal et al., 2015; Sklar et al., 2017).

Hillslope sediment is transported towards the river channel by a variety of erosion processes. Diffusive processes are considered to operate in slowly eroding soil-mantled landscapes and move relatively fine grains at or near the surface (Roering et al., 1999). In contrast, deep-seated erosion processes (e.g. landslides) are frequent in steep and rapidly eroding bedrock landscapes (e.g. Burbank et al., 1996; Hovius et al., 1997; Montgomery and Brandon, 2002). Especially when a critical threshold hillslope angle of ∼25–30 is exceeded, denudation rates are dominated by the frequency of landslides (Larsen and Montgomery, 2012; Montgomery and Brandon, 2002; Ouimet et al., 2009). Deep-seated erosion processes excavate sediment and bedrock fragments of any size from a greater depth interval (Casagli et al., 2003).

The type of erosion process is indirectly controlled by tectonic uplift rates, precipitation, and lithology. In a steady state landscape, denudation rates are set by tectonic uplift, which controls river incision and hillslope steepening (DiBiase et al., 2010; Scherler et al., 2014; Whipple and Tucker, 1999). Extreme precipitation events may over-saturate hillslopes and increase the susceptibility to hillslope failure (Chen et al., 2006; Gabet and Dunne, 2002). Some authors argue that deep-seated erosion processes are also important in arid landscapes (Aguilar et al., 2014). Finally, the bedrock strength and fracture abundance affects the susceptibility to hillslope failure (e.g. Clarke and Burbank, 2011; Perras and Diederichs, 2014) and constrains how much energy is needed for the detachment and transport of individual particles.

Once the sediment has reached the channel, processes like downstream abrasion, selective transport, and mixing of sediment sources control the grain size distribution at the sample location. If mixing within the channel is incomplete, single tributaries or local inputs of sediment (e.g. landslides) may dominate the grain size distribution (Binnie et al., 2006; Neilson et al., 2017; Niemi et al., 2005; Yanites et al., 2009). Downstream abrasion and selective transport result in a progressively smaller grain size distribution. Abrasion wears off the outer layers of clasts (Kodoma, 1994a, b; Sklar et al., 2006) and depends on the travel distance and velocity as well as the lithology of the clasts (Attal and Lavé, 2009). Selective transport preferentially deposits coarse grains when the transport capacity of water is low (Ferguson et al., 1996; Hoey and Ferguson, 1994) and thus further changes the grain size distribution.

## 2.2 Processes that control variations in 10Be concentrations

10Be concentrations in quartz grains depend on the 10Be production rate and the exposure time of the grain to cosmic rays (Gosse and Phillips, 2001). Processes that preferentially transport sediment from locations with contrasting 10Be production rates (i.e. different soil depths or elevations within the catchment) or exposure times (i.e. locations with different denudation rates) result in a larger variation in 10Be concentrations in the sediment. A certain sample location provides a spatially and temporally representative sediment sample when the sediment from different sources is sufficiently mixed (Binnie et al., 2006; Neilson et al., 2017; Niemi et al., 2005; Yanites et al., 2009). Because 10Be production rates decrease exponentially with depth (Gosse and Phillips, 2001), hillslope sediment that is excavated over a larger depth interval by landslides will obtain a larger variation in 10Be concentrations than sediment transported by diffusive processes near the surface (Aguilar et al., 2014; Belmont et al., 2007; Binnie et al., 2007; Brown et al., 1995; Puchol et al., 2014; Sosa Gonzalez et al., 2016a, b; Tofelde et al., 2018; West et al., 2014). In soil-mantled landscapes, bioturbation by burrowing animals and tree throw (Gabet et al., 2003) result in well-mixed surface layers with uniform 10Be concentrations (Brown et al., 1995; Granger et al., 1996; Schaller et al., 2018) (Fig. 2). These mixed soil layers are most likely to develop in humid and slowly eroding catchments, where biota is abundant. Eroded sediment from these layers is expected to have uniform 10Be concentrations. In rapidly eroding and arid landscapes, however, soils are typically very thin or absent, and the eroded sediment likely yield larger variations in 10Be concentrations (Fig. 2). Furthermore, fluvial processes can affect the grain size distribution at the sample location in a way that not all parts of a catchment are equally represented in different grain size fractions. For example, sediment provenance of grains from different elevations could play a role in catchments with heterogeneous rock types that produce different clast sizes or contain different quartz abundances (Bierman and Steig, 1996; Carretier et al., 2015). An unequal representation of elevations in different grain size fractions may also result from hydrodynamic sorting, downstream abrasion, and insufficient mixing of tributaries that drain different elevations (Carretier et al., 2009; Carretier and Regard, 2011; Lukens et al., 2016; Neilson et al., 2017). Combined with elevation-dependent 10Be production rates (and provided that denudation rates are constant), this could also result in grain-size-dependent 10Be concentrations (Carretier et al., 2015; Lukens et al., 2016; Matmon et al., 2003).

Figure 2The effect of hillslope steepness and precipitation on the depth of the mixed soil layer and 10Be concentrations as a function of depth. (a) Hillslope steepness and denudation rates control the thickness of the soil mantle by the removal of material from the top. A thick soil mantle likely develops in gently sloping and slowly eroding landscapes, whereas high denudation rates in steep landscapes prohibit the development of a thick soil mantle. (b) Precipitation provides water for chemical weathering. Humid landscapes likely develop a thick soil mantle, which may be absent in arid landscapes. Bioturbation in landscapes with thick soil mantles results in a well-mixed soil layer with a uniform 10Be concentration, which, in isotopic steady state, is equal to the surface concentration. In landscapes where a mixed soil layer is absent, 10Be concentrations decrease exponentially with depth.

## 2.3 Catchment attributes that potentially control grain-size-dependent 10Be concentrations

Based on the above review of processes that may influence grain-size-dependent 10Be concentrations, we selected the catchment attributes mean basin slope, mean travel distance, mean annual precipitation (MAP), and lithology as the main focus of our study. Here we consider mean basin slope as a topographic catchment attribute that controls denudation rates and the scouring depth of diffusive or deep-seated erosion processes. We selected MAP because of its effect on both weathering rates and the scouring depth of erosion processes and lithology because it affects chemical weathering rates, the grain size of individual minerals, and the susceptibility to hillslope failure. Finally, we selected mean travel distance of sediment as a metric for fluvial processes that are transport-dependent (e.g. abrasion and hydrodynamic sorting).

Table 1Sample location coordinates and catchment attributes of the research areas located in the Chilean Coastal Cordillera.

a Mean annual precipitation (MAP) estimates derived from the Global Precipitation Climatology Centre (GPCC) dataset (Meyer-Christoffer et al., 2015). b Total mean basin slope calculated with a 30 m digital elevation model (DEM). c Normalized channel steepness index.

3 Study area

For our case study, we selected four small catchments (< 10 km2) located in the Coastal Cordillera of central Chile (Table 1, Fig. 3). The Coastal Cordillera features a pronounced latitudinal climate and vegetation gradient, whereas the tectonic setting is rather uniform. The selected catchments are located in the National Park Pan de Azúcar (AZ) (∼26 S), the National Reserve Santa Gracia (SG) (∼30 S), the National Park La Campana (LC) (∼33 S) and the National Park Nahuelbuta (NA) (∼38 S). The catchments share a granodioritic lithology, though some minor variations in mineralogy exist between the sites (Oeser et al., 2018). The three northern catchments experience modern uplift rates of < 0.1 mm yr−1 (Melnick, 2016). The southernmost catchment is located in the Nahuelbuta Range, where uplift rates increased from 0.03–0.04 to > 0.2 mm yr−1 at 4±1.2 Ma (Glodny et al., 2008; Melnick et al., 2009). Because the sampled catchment is located upstream of a river channel knickpoint, it may not yet be influenced by the increased uplift rates (Crosby and Whipple, 2006). The climatic regime and mean annual precipitation (MAP) range from arid (MAP ∼13 mm yr−1) in Pan de Azúcar in the north to semi-arid in Santa Gracia (MAP ∼88 mm yr−1), Mediterranean in La Campana (MAP ∼358 mm yr−1), and temperate in Nahuelbuta (MAP ∼1213 mm yr−1) in the south (Meyer-Christoffer et al., 2015). This latitudinal increase in MAP results in an increase in vegetation density (Bernhard et al., 2018). The Normalized Difference Vegetation Index (NDVI) (Didan, 2015) varies from 0.1 in the northernmost catchment to 0.8 in the southernmost catchment (Fig. S1 in the Supplement). The increase in MAP is accompanied by an increase in chemical weathering rates measured in soil profiles located within or in proximity of the catchments (Oeser et al., 2018). The chemical depletion fraction (CDF), a measure to quantify chemical weathering (Riebe et al., 2003), increases from ∼0.1 in Pan de Azúcar (AZ) to ∼0.4–0.5 in Santa Gracia (SG) and ∼0.3–0.6 in La Campana (LC). Due to heterogeneities in bedrock samples collected in Nahuelbuta (NA), no reliable CDF could be assigned. 10Be depth profiles measured in two mid-slope soil profiles revealed an increasing thickness of the mixed soil layer, presumably due to bioturbation (Schaller et al., 2018). The depth of the mixed soil layer increases from ∼0–17.5 cm in Pan de Azúcar (AZ) to ∼25–45 cm in Santa Gracia (SG), ∼47.5–85 cm in La Campana, and ∼70 cm in Nahuelbuta (NA) (Schaller et al., 2018). By selecting small and low-relief catchments with similar lithology and a relatively uniform tectonic setting, we aim to explore the relationship between grain size and 10Be concentrations as controlled by precipitation.

Figure 3Research areas located on a precipitation gradient in the Chilean Coastal Cordillera. (a) The catchment locations (stars) on a mean annual precipitation (MAP) map from Climate Hazards Group InfraRed Precipitation (CHIRPS), underlain by a Shuttle Radar Topography Mission (SRTM) DEM-derived hillshade map. (b) Google Earth images showing the sample locations (stars), catchment outlines (white), channels with an upstream area of 0.2 km2 (blue lines), and the locations of two 10Be depth profiles measured by Schaller et al. (2018) in soil pits on a north-facing (N) and south-facing (S) slope (green dots).

4 Methods

## 4.1 Sampling and analytical methods

In each of the four Chilean catchments, we collected approximately 6 kg sand and pebbles from the active channel (Fig. 3) and conducted a Wolman pebble count (Wolman, 1954) with 1 m intervals to measure the grain size distribution at the sample locations (Fig. S2). We dried and sieved the samples in the laboratory to separate the grain size fractions 0.5–1, 1–2, 2–4, 4–8, 8–16, 16–32, and 32–64 mm. Before further processing, we crushed pebbles larger than 1 mm. To isolate pure quartz, we separated and purified the river sediment using standard physical and chemical separation methods (Kohl and Nishiizumi, 1992). We spiked between 10 and 20 g of pure quartz with 0.15 mg 9Be carrier, dissolved the quartz, and extracted beryllium following established protocols (e.g. von Blanckenburg et al., 2004). Accelerator mass spectrometry measurements were carried out at the University of Cologne, Germany. Reported 10Be∕9Be ratios have been normalized to the KN01-6-2 and KN01-5-3 standards, with nominal 10Be∕9Be ratios of $\mathrm{5.35}×{\mathrm{10}}^{-\mathrm{13}}$ and $\mathrm{6.32}×{\mathrm{10}}^{-\mathrm{12}}$, respectively. We calculated 10Be concentrations from 10Be∕9Be ratios and a blank correction was performed. We used MATLAB® and the CRONUS functions (Balco et al., 2008) with the time-independent (St) scaling scheme (Lal, 1991; Stone, 2000) and the sea level high latitude (SLHL) production rate of 4.01 atoms g−1 yr−1 (Borchers et al., 2016; Phillips et al., 2016) to calculate catchment-average denudation rate estimates from the 10Be concentrations

## 4.2 Global compilation

We compiled data from previously published studies that measured 10Be concentrations in different grain size fractions sampled at the same location. Because we are interested in small- to medium-sized bedrock catchments, and to reduce the effect of long-term floodplain sediment storage in large basins, we discarded basins with an area of > 5000 km2. We also removed studies that only measured sand fractions (< 2 mm) as weathering and erosion processes affecting these sand-sized grain size fractions may be similar. Hence, we only selected studies measuring at least one sand fraction (mean grain size < 2 mm) and at least one coarser grain size fraction (mean grain size > 2 mm) (Aguilar et al., 2014; Belmont et al., 2007; Brown et al., 1995; Carretier et al., 2015; Clapp et al., 2002; Derrieux et al., 2014; Heimsath et al., 2009; Matmon et al., 2003; Palumbo et al., 2010; Puchol et al., 2014; Reinhardt et al., 2007; Stock et al., 2009; Sullivan, 2007; Tofelde et al., 2018; Granger et al. 1996; Safran et al., 2005). From each selected catchment we compiled the reported grain size classes, the corresponding 10Be concentrations (± analytical errors), and the sample location coordinates (Table S3). For studies that reported a grain size fraction as “larger than”, we assumed that the upper grain size limit corresponds to twice the lower limit (e.g. reported: > 2 mm, data compilation: 2–4 mm). We acknowledge that this range might be incorrect, but a fixed grain size range was required for proper data analysis. We transformed the measured grain sizes to φ-based grain size classes, which is the negative logarithm to the base 2 of the grain size diameter (Krumbein, 1934, 1938). The range of grain sizes we investigated (0.063 to 200 mm) corresponds to φ values of −4 to 7.64. To compare data from different study areas with different 10Be production rates, we normalized the 10Be concentrations (± analytical uncertainties) by the arithmetic mean concentration of all samples from the same catchment.

To assess the influence of the identified catchment attributes (mean basin slope, mean travel distance, mean annual precipitation, and lithology) (Sect. 2.3) on grain-size trends in the global compilation, we used a 90 m resolution SRTM DEM (Jarvis et al., 2008). We obtained upstream areas based on the published sample coordinates and using the flow routing tools of the TopoToolbox v2 (Schwanghart and Scherler, 2014). We calculated the following topographic parameters: catchment area, mean basin slope, total relief (maximum elevation minus minimum elevation), and the mean travel distance of sediment to the sample location, which is calculated as the arithmetic mean travel distance of all pixels in the catchment to the sample location). The agreement between the published and recalculated topographic parameters is good, and minor deviations likely result from differences in DEM resolution (Fig. S3). We obtained an estimate of mean annual precipitation (MAP) in each catchment using the 0.25 resolution gridded precipitation dataset from the Global Precipitation Climatology Centre (Meyer-Christoffer et al., 2015). To classify catchment lithology we used the Global Lithological Map (GLiM; Hartmann and Moosdorf, 2012) together with the lithology reported in the original publications. We defined four different lithological classes: sedimentary, magmatic, metamorphic, and mixed (> 3 different rock types in a catchment).

Next, we used Sternberg's law to estimate the extent of abrasion of bedload sediment during fluvial transport to define a travel distance threshold after which abrasion becomes significant:

$\begin{array}{}\text{(1)}& D\left(L\right)={D}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}{e}^{-\mathit{\alpha }L}.\end{array}$

Using Eq. (1), we calculated the grain size D at the sample location, which is derived from an initial grain size D0 at the source, that travelled distance L and decreased in size at a rate given by the reduction coefficient α (Kodama, 1994; Lewin and Brewer, 2002; Sklar et al., 2006; Sklar and Dietrich, 2008). The reduction coefficient depends on both grain velocity and lithology. Rocks with low tensile strength reduce faster in size during transport than rocks with high tensile strength (Attal and Lavé, 2009). We chose the reduction coefficients based on literature values for field settings (sedimentary rocks: α=0.0003 m−1, magmatic rocks: α=0.0002 m−1, metamorphic rocks: α=0.0001 m−1), which are typically higher than experimental studies due to different particle collision dynamics and the lack of weathering in experimental studies (Sklar et al., 2006). We considered the effect of abrasion to be negligible when a grain size at the sample location (D) falls in the same phi-grain-size class as at its erosion source (D0). As an example, for abrasion to be significant, a grain size of 2 mm at the erosion source must be smaller than 1 mm at the sample location to fall in a lower phi-grain-size class. This results in abrasion thresholds for sedimentary, magmatic, and metamorphic rocks of 2300, 3500, and 7000 m, respectively. For catchments underlain by mixed lithologies, the abrasion threshold lies between 2300 and 7000 m (Attal and Lavé, 2009).

We quantified the relationship between grain size and 10Be concentrations by calculating a “grain size dependency” for each sample set (Fig. 4). This is the slope of a linear fit through the 10Be concentrations of different grain size classes. To account for uncertainties in 10Be concentrations and for grain size ranges, we used a Monte Carlo approach (n=10 000) to randomly select a point between the mean (± analytical error) 10Be concentrations and the analysed grain size range. We thus obtained a mean (± standard deviation) grain size dependency for each catchment. A positive grain size dependency indicates higher 10Be concentrations in coarser grains, and vice versa.

Figure 4Schematic showing the concept of measuring grain size dependency. A random point was selected within each grain size range and from the corresponding 10Be concentrations (± analytical error) (boxes). The slope of a line fitted to the randomly selected points represents the grain size dependency. We used a Monte Carlo simulation of 1000 runs to account for the width of the grain size range and the analytical errors in 10Be concentrations. This yields a mean grain size dependency with error bars.

Next, we used the Kolmogorov–Smirnov test (KS test, 5 % significance interval) to test whether a particular mean basin slope, MAP, or sediment travel distance class showed a significantly different distribution of grain size dependencies (Kolmogorov, 1933; Smirnov, 1939). Finally, we calculated linear regression statistics between the grain size dependency values and the catchment attributes mean basin slope, MAP, and mean travel distance and applied a multivariate linear regression model including the effect of all 3 catchment attributes. We did this for the entire dataset and for each individual lithology. As part of the multivariate statistics, we calculated the relative importance (RI) of all catchment attributes, using the approach proposed by Lindeman, Merenda and Gold (LGM) (Lindeman et al., 1980) of the “Relaimpo” R studio package (Grömping, 2006). This provides the percentage contribution of each catchment attribute to the multivariate regression model R2.

Table 2Cosmogenic nuclide samples from the Chilean Coastal Cordillera. IGSN number, analysed quartz mass, 9Be carrier mass, 10Be∕9Be ratio (±1σ), 10Be concentrations (±2σ analytical error), spallation (Psp) and muogenic (Pmu) production rates, and calculated denudation rates (±2σ).

5 Results

## 5.1 Chilean Coastal Cordillera

The measured 10Be concentrations in the most arid catchment (AZ) range from 2.8 to 4.6×105 atoms (g quartz)−1, resulting in catchment-average denudation rates of 5.8±0.7 to 10.1±1.1 mm kyr−1 (Table 2). In the semi-arid catchment (SG), the 10Be concentrations range from 3.6 to 5.2×105 atoms (g quartz)−1, which corresponds to catchment-average denudation rates of 7.5±0.8 to 11.0±1.4 mm kyr−1 (Table 2). The 10Be concentrations in the Mediterranean catchment (LC) are a factor 10 lower compared to the other catchments and range from 0.2 to 0.6×105 atoms (g quartz)−1, which results in catchment-average denudation rates of 103.7±12.4 to 384.1±54.5 mm kyr−1 (Table 2). The temperate catchment (NA) yielded 10Be concentrations ranging from 1.8 to 2.9×105 atoms (g quartz)−1, resulting in catchment-average denudation rates of 24.0±2.6 to 40.2±4.5 mm kyr−1 (Table 2). Only the arid (AZ) and Mediterranean (LC) catchments show a consistent, but noisy, trend between 10Be concentrations and grain sizes. In the arid catchment (AZ), 10Be concentrations are decreasing with increasing grain size. The 2σ variability of 10Be concentrations measured in all grain size fractions deviates ±18 % from the mean (Fig. 5). In the Mediterranean catchment (LC), the 10Be concentrations of all grain size fractions vary up to ±40 % from the mean and display a noisy but positive grain size dependency, i.e. increasing 10Be concentrations with increasing grain size (Fig. 5). In both the semi-arid (SG) and temperate catchments (NA), the 2σ variability in 10Be concentrations is low (±12 % and ±14 %, respectively) and rather unsystematic (Fig. 5). The smallest grain size fractions (0.5–4 mm) in the semi-arid catchment (SG) show a decreasing trend, but this trend increases again for coarser grain size fractions (4–32 mm). In the temperate catchment (NA), 10Be concentrations are uniform in the five smallest grain size fractions (0.5–16 mm), but this trend breaks down at the two largest grain size fractions (16–64 mm), which have lower 10Be concentrations.

Figure 5Normalized 10Be concentrations (±2σ analytical error) measured in seven different grain size classes. The 10Be concentrations are normalized to the arithmetic mean of all grain size fractions within a catchment. The red line indicates the normalized catchment-average 10Be concentration.

## 5.2 Global compilation

The global compilation includes 73 catchments covering a wide range of different hillslope angles, sediment travel distances, MAP, and lithologies (Fig. 6). Figure 7 shows the data of all catchments, classified in four slope classes and colour coded by lithology. Each box represents the normalized 10Be concentrations (± analytical uncertainties) and the grain size range of a single sample. Uncertainties in 10Be concentrations tend to be larger for samples from steeper catchments (> 10), which may be related to higher denudation rates and therefore lower 10Be concentrations. Generally, uncertainties are larger for low 10Be∕9Be ratios. In catchments with mean basin hillslope angles < 10, 10Be concentrations are relatively similar across all grain size classes. In steeper hillslope classes, coarse grains reveal lower 10Be concentrations compared to fine grains, with the largest deviations in catchments with hillslope angles > 25 (Fig. 7). We discern no pattern related to lithology from this figure but we emphasize that magmatic lithologies are more abundant in shallow sloping catchments, whereas metamorphic lithologies are more abundant in steep catchments.

Figure 6Sample locations and catchment attributes of all catchments in the global compilation (n=73). (a) Coarsest phi-grain-size classes measured in each study (the smallest grain size was always a sand fraction < 2 mm). (b) Mean travel distance of sediment, calculated as the arithmetic mean of each grid cell's travel distance towards the catchment outlet. (c) Mean basin slope of each catchment, calculated as the arithmetic mean of the hillslope angles at each grid cell. (d) Mean annual precipitation (MAP) in each catchment, derived from the Global Precipitation Climatology Centre (GPCC) dataset (Meyer-Christoffer et al., 2015).

Figure 8 (and Fig. S4 for plots separated by lithology) shows the grain size dependencies of individual catchments, resulting from the slope of a linear fit to the 10Be concentrations of all grain size classes (see Methods section and Fig. 4). Overall, we observe more sample sets that display significantly (i.e. error bar does not overlap with 0) negative (56.2 %) trends in grain-size-dependent 10Be concentrations than positive (32.8 %; Fig. 8). A percentage of 11.0 % of the sample sets have grain size dependencies that are not significantly different from zero, and thus reveal no grain size dependency. Furthermore, negative grain size dependencies are typically stronger (i.e. higher absolute differences between grain sizes) than positive grain size dependencies.

The calculated grain size dependencies reveal a significant break point at a mean hillslope angle of ∼15 (KS test, Fig. 8a). Catchments with mean hillslope angles < 15 reveal a distribution with predominantly weak grain size dependencies. Steep catchments with hillslope angles > 15 show a wider distribution with predominantly negative grain size dependencies (62.3 % significantly negative). A percentage of 70.0 % of the catchments that exceed the threshold hillslope (> 25) have significantly negative grain size dependencies. Our analysis of sediment travel distance shows that the amount and magnitude of negative grain size dependencies slightly increase at longer sediment travel distances (Fig. 8b). However, catchments that exceeded the abrasion threshold (sedimentary: 2300 m, magmatic: 3500 m, metamorphic: 7000 m, mixed: 2300–7000 m) show no significantly different grain size dependency distribution based on the KS test. Finally, the data suggest a slightly increasing amount and magnitude of negative grain size dependencies with increasing MAP. Humid catchments (MAP > 2000 mm yr−1) reveal a distribution of predominantly (90 %) significantly negative grain size dependencies, which is significantly different (KS test) from catchments with MAP < 2000 mm yr−1 (Fig. 8c). However, only a low number of catchments with MAP > 2000 mm yr−1 compose the distribution. Catchments underlain by sedimentary and metamorphic rocks show the most significant negative grain size dependencies (66.7 % and 65.4 %, respectively), followed by catchments underlain by mixed lithologies (50.0 %). The number of significantly negative grain size dependencies is lowest for catchments underlain by magmatic lithologies (37.5 %). None of the lithologies revealed a significantly different grain size dependency distribution based on the KS test.

Figure 7All global compilation samples divided into mean basin slope classes and colour coded by lithology. The boxes indicate normalized mean 10Be concentrations (± analytical errors) and the phi-grain-size range. The 10Be concentrations are normalized by the arithmetic mean 10Be concentration of all samples from the same catchment.

Linear regressions of grain size dependencies as a function of mean basin slope revealed significantly negative trends for all lithologies combined (p=0.002) and for metamorphic catchments (p=0.017) but not for the other lithologies alone (Table S4, Fig. S4). MAP showed a significantly negative relationship with grain size dependencies for all lithologies combined (p=0.007), and for catchments underlain by magmatic (p=0.006) lithology. However, the trend for magmatic catchments mainly results from one negative data point at higher MAP. No significant linear trends emerged between mean travel distance and grain size dependencies for any of the lithologies.

When considering the combined influence of mean basin slope, MAP, and mean travel distance with a multivariate linear model, we found that the variance of all lithologies combined is significantly described (p=0.004) by two out of all three factors, but that the explained variance is low (R2=0.190, Table S5). Most of the variance is related to mean basin slope (relative importance, RI =9.1 %), followed by MAP (RI =7.6 %), whereas mean travel distance revealed no significant contribution (Table S5, Fig. 9). Furthermore, multivariate models yielded significant results when considering only magmatic (p=0.031, R2=0.552) and metamorphic catchments (p=0.077, R2=0.276). In magmatic catchments, a large proportion of the variability is ascribed to MAP (RI =51.4 %), which results from the above mentioned negative data point at higher MAP. Finally, most of the variability in magmatic catchments was significantly described by mean basin slope (RI =22.0 %). Multivariate statistics yielded insignificant results for mixed and sedimentary lithologies, possibly due to too few catchments to disclose unambiguous trends.

Figure 8Grain size dependencies of all catchments in the global compilation (n=73) as a function of (a) mean basin slope, (b) mean travel distance, and (c) mean annual precipitation (MAP). Coloured symbols depict lithological classes. Grain size dependencies are derived from the slope of a linear fit to the normalized 10Be concentrations and grain sizes from a sample set, as described in Fig. 4. Dashed lines indicate the threshold hillslope (a) and abrasion thresholds (b) mentioned in the text. Global compilation statistics are provided in Table S4 and Fig. S4.

6 Discussion

## 6.1 Grain-size-dependent 10Be concentrations in the Chilean Coastal Cordillera

The sampled catchments on the climatic gradient in the Chilean Coastal Cordillera only show a systematic trend of 10Be concentrations with grain size in the arid (AZ) and Mediterranean catchments (LC). In both catchments, the 10Be concentrations of river sediment correspond to concentrations measured in the subsurface of the soil profiles (Fig. 10; Schaller et al., 2018). Because the difference between 10Be production rates of the catchment on average and at the soil profiles is small (< 10 %), we can compare measured 10Be concentrations directly. In the arid catchment (AZ), both the negative grain size dependencies and the fact that 10Be concentrations correspond to concentrations at ∼1 m depth in the soil profiles suggest that erosion processes (e.g. rock falls, landslides, gully head retreat) – which excavate sediment from intermediate to greater depth during rare precipitation events or earthquakes (e.g. Mather et al., 2014; Pinto et al., 2008) – may occur in this catchment. All of the measured 10Be concentrations in river sediment in the Mediterranean catchment (LC) are considerably lower than the concentrations measured at the surface of soil profiles in a close proximity to the catchment (Fig. 10; Schaller et al., 2018). This suggests that the catchment experiences faster erosion processes compared to the location of the soil pit, which is confirmed by debris flow scars observed at high elevation in the catchment (Fig. S5). Deep-seated erosion processes and insufficient mixing in a small-sized catchment may make a sample non-representative for the entire catchment (Niemi et al., 2005; Yanites et al., 2009). However, the noisy, but overall positive, grain size dependency in the Mediterranean catchment (LC) contradicts this hypothesis (Fig. S5), as debris flows would presumably excavate coarse grains from greater depth. Higher 10Be production rates at the elevation where debris flows originate, and the condition that coarse grains only originate from that area, cannot account for the positive grain size dependency alone (Fig. S5). Without being able to clarify this issue, the lower 10Be concentrations of river sediment, combined with the observed greater scatter in the positive grain size dependency, may hint at selective transport and longer residence times of coarse grains at higher elevations.

Figure 9Relative importance of mean basin slope, MAP, and mean travel distance to the multivariate linear regression R2 value. Results are given for all lithologies combined (grey) and differentiated by lithology (colours). Multivariate linear regression results are provided in Table S5.

Figure 10Comparison of 10Be concentrations measured in (a) river sediment, and (b) north- and south-facing soil profiles (Schaller et al., 2018), from the same catchments in the Chilean Coastal Cordillera. MDN and MDS are the soil mixing depths of the north- and south-facing hillslopes, respectively. Note the reduced x axis range of the Mediterranean catchment (LC). Shaded areas show range of 10Be concentrations in river sediment for comparison with soil profiles.

The 10Be concentrations in river sediment from the semi-arid (SG) and temperate (NA) catchments show little variations and are similar to concentrations measured near the surface in soil pits (Fig. 10; Schaller et al., 2018). Within the temperate catchment (NA), the uniform 10Be concentrations in grains < 16 mm suggests that these originate from the ∼70 cm thick mixed soil layer, whereas the lower 10Be concentrations in grains > 16 mm suggest these may be derived from below the mixed layer (Fig. 10; Schaller et al., 2018). In the semi-arid (SG) catchment, the measured samples from the channel show similar 10Be concentrations compared to those measured in the mixed soil layer of the north-facing hillslope and higher 10Be concentrations compared to the mixed layer of the south-facing hillslope (Fig. 10; Schaller et al., 2018). This suggests that grains are unlikely to be derived from greater depth, where 10Be concentrations are lower.

We propose that the existing or missing trends in the arid (AZ), semi-arid (SG), and temperate (NA) catchments are mainly related to differences in precipitation and the excavation depth of the erosion processes. These catchments show minor variations in mean basin slope, and hence we do not expect big differences in erosion processes due to changes in slope alone. Furthermore, the limited relief of these catchments excludes differences in 10Be production rates and local sediment sources to influence observed differences in 10Be concentrations. However, steeper hillslope angles and higher total relief may have overruled the effect of precipitation in the La Campana catchment. We do not expect a control related to the different catchment sizes in any of the catchments because granitic rocks have a low abrasion breakdown rate (Attal and Lavé, 2009) and the mean travel distances were small (< 1 km).

In summary, we think our new samples from the Chilean Coastal Cordillera suggest an influence of MAP on grain-size-dependent 10Be concentrations only in the most-arid and most-humid catchments by its effect on the thickness of the mixed soil layer and the scouring depth of erosion processes that transport larger grains from below the mixed soil layer.

## 6.2 Grain-size-dependent 10Be concentrations in the global compilation

### 6.2.1 Mean basin slope

The effect of mean basin slope on grain-size-dependent 10Be concentrations is apparent as weak grain size dependencies in gently sloping catchments, and predominantly negative grain size trends in steep catchments (Figs. 7, 8a). Mean basin slope may control grain-size-dependent 10Be concentrations through its effect on the thickness of soils and the scouring depth of erosion processes. In gently sloping catchments, denudation rates are typically low (e.g. Portenga and Bierman, 2011) and well-mixed soil layers with uniform 10Be concentrations can develop. Diffusive erosion processes transport sediment from near the surface, which results in uniform 10Be concentrations. In contrast, in steep landscapes, denudation rates are usually high and soils are thin or absent if denudation rates exceed the soil production limit (∼170 mm kyr−1; Dixon and von Blanckenburg, 2012). Such catchments are typically dominated by deep-seated hillslope processes (Hovius et al., 1997). Negative grain size dependencies thus occur because coarse grains are excavated from greater depth, where 10Be concentrations are lower. The highest percentage of negative grain size dependencies are found in catchments steeper than 25. In these catchments, many hillslopes have likely reached the threshold hillslope angle of ∼25–30 at which denudation rates are dominated by the frequency of landslides (Larsen and Montgomery, 2012; Montgomery and Brandon, 2002; Ouimet et al., 2009). Linear regression models revealed a stronger control of mean basin slope on grain-size-dependent 10Be concentrations than MAP and mean travel distance; however, the R2 values of the regression models were low (Table S5, Fig. 9). This conforms with previous studies that also found negative grain size dependencies which emerged from a transition of transport-limited to detachment-limited erosion processes, and therefore deep-seated erosion processes (Binnie et al., 2007; Brown et al., 1995; Lukens et al., 2016; Reinhardt et al., 2007; Sosa Gonzalez et al., 2016a, b; Tofelde et al., 2018). It is notable that the most-negative grain size dependencies occur in catchments underlain by sedimentary rocks (Table S4 and Fig. S4). This may be due to lower rock mass strength of sedimentary rocks, which partly stems from the presence of bedding planes, making them more susceptible to hillslope failure (e.g. Clarke and Burbank, 2011; Perras and Diederichs, 2014).

### 6.2.2 Sediment travel distance

Our results revealed a weak negative control of sediment travel distance on grain size dependencies; however, no significant relationships were found. The negative control is strongest for sedimentary catchments in which negative grain size dependencies appear to be more frequent in catchments with long sediment travel distances, and the most negative grain size dependencies appear when travel distances exceeded the abrasion threshold (Fig. 8b). Possibly the lower rock strength of sedimentary rocks promotes the breakdown into smaller particles and increases the grain's sensitivity to abrasion (Attal and Lavé, 2009; Sklar and Dietrich, 2001). Due to abrasion, distant erosion sources may be over-represented in finer grain size fractions and under-represented in coarser ones (Lukens et al., 2016). As travel distance scales with elevation (Fig. S6), and therefore 10Be production rates, sediment from high elevations may have inherently higher nuclide concentrations (Lal, 1991). In contrast, coarse grains, which experienced less abrasion, may originate from lower elevations, with lower 10Be production rates. This elevation dependence of certain grain size fractions may induce a negative grain size-dependency. Secondly, if abrasion were to reduce river sediment of decimetre or metre scales to sand size, the centre of such clasts would have lower concentrations (Carretier and Regard, 2011; Lupker et al., 2017). However, the associated travel distance has to be considerably longer, and the initial clast must be large. For example, abrasion of an initial 25 cm sized granitic cobble over a distance of ∼8 km would result in a size reduction of 10 cm and expose a centre with a 10Be concentration that is only 8.5 % lower compared to the outer layers (Balco et al., 2008; Sklar et al., 2006). The by-product of abrasion, which typically is of silt or clay size (Sklar et al., 2006), unlikely affects the measured 10Be concentrations, as it is finer than the grain size classes typically analysed (Lukens et al., 2016). We did not observe a control of sediment travel distance in catchments with mixed lithologies. The provenance of distinct grain sizes from different lithologies has not resulted in a dominantly positive or negative grain size dependency, possibly because the spatial arrangement of different lithologies in a landscape is not necessarily elevation-dependent, or because these lithologies yield minor differences in grain sizes.

### 6.2.3 Mean annual precipitation

The global compilation suggested an additional control of MAP on grain-size-dependent 10Be concentrations. The amount and magnitude of negative grain size trends seem to increase with increasing MAP. The highest percentage of negative grain size dependencies is found in humid catchments (> 2000 mm yr−1). However, this trend is related to a low number of catchments. Negative grain size dependencies at higher MAP values could be related to higher denudation rates and increasing depth of erosion processes (e.g. precipitation-induced landslides; Chang et al., 2007; Chen et al., 2006; Lin et al., 2008). This differs from our interpretation of the results from the Chilean Coastal Cordillera in which we emphasize the control of MAP on the thickness of the mixed soil layer. The discrepancy with the global compilation may result from the additional effect of hillslope angle, which also influences the thickness of the soil mantle and the depth of erosion processes (Heimsath et al., 2009).

## 6.3 Implications

Our results and the above discussion suggest that grain size trends in 10Be concentrations are best explained by the effects of hillslope angle and MAP on the presence and thickness of mixed soil layers and the scouring depth of erosion processes. In large catchments, an additional effect may emerge by abrasion during fluvial transport, which could induce a non-representative grain size distribution. At present, however, it is difficult to quantify the relative roles of hillslope angle, precipitation, travel distance, and lithology, because these parameters tend to be partly correlated. For example, high and steep topography is often associated with high amounts of orographic precipitation, and long travel distances are associated with high total relief (Fig. S6).

In any case, the presumed role that soils and different hillslope erosion processes play for grain-size-dependent 10Be concentrations is likely not linearly related to variables like mean hillslope angle or mean annual precipitation. Instead, our results are consistent with the presence of thresholds. Landslides likely become important when hillslope angles exceed a critical threshold (Burbank et al., 1996) and once precipitation is high enough to sustain vegetation and soils, diffusive processes may dominate gently sloping and soil-mantled landscapes. Such a threshold control on the occurrence of grain-size-dependent 10Be concentrations may be the reason why our linear regression statistics yielded mostly insignificant results or low R2 values (Fig. 9 and Table S5). More data may allow for better constraining of the controls and relative importance of these factors in the future. It additionally highlights the importance of systematic studies on single factors, like our study on the sole effect of MAP in the Chilean Coastal Cordillera.

We evaluated the likelihood of grain-size-dependent 10Be concentrations and a potential bias in previously published 10Be-derived catchment-average denudation rates by comparing our findings with a recently published global compilation (Codilean et al., 2018). Out of 2537 different catchments with an area < 5000 km2, 55.7 % have hillslope angles > 15, where our data first show significant grain size effects, and 23.3 % have hillslope angles > 25. When considering sediment travel distances, using the relationship between catchment area and sediment travel distance that emerged from our global compilation (R2=0.99; Fig. S6), about 61.9 %, 49.8 %, and 29.2 % of the catchments exceeded the sediment travel distances of 2300, 3500, and 7000 m, respectively. Finally, 11.5 % of the catchments have MAP > 2000 mm yr−1, based on GPCC-derived MAP (Meyer-Christoffer et al., 2015) at the sample location. Therefore, previously published catchment-average denudation rates may more frequently be biased as a result of steep hillslopes and long sediment travel distance and less frequently by the influence of MAP. When considering a combined effect of all controlling factors in each catchment (slope > 25, sediment travel distance > 7000 m, and MAP > 2000 mm yr−1), 49.1 % of the catchments are predicted to be devoid of grain size dependencies of 10Be concentrations and biased catchment-average denudation rates, whereas 50.9 % might contain a bias because one or more of the controlling factors has exceeded the threshold values that emerged from our study.

7 Conclusions

In this paper, we used a field study in Chile and a global compilation of previously published data to assess in what types of catchments grain-size-dependent 10Be concentrations may lead to biased estimates in catchment-average denudation rates. Our results suggest that mean basin slope and MAP control grain-size-dependent 10Be concentration through their effect on the presence and thickness of a mixed soil layer and the depth of erosion processes. Hillslope steepness appears to exert the most important influence on grain-size-dependent 10Be concentrations. Our global compilation results show that the influence of MAP is limited to humid catchments (> 2000 mm yr−1), whereas our case study in Chile suggests an additional control in arid catchments (< 100 mm yr−1). Furthermore, grain-size-dependent 10Be concentrations may occur in large catchments with long sediment travel distances (> 2300 to > 7000 m, depending on lithology), where abrasion may induce non-representative grain size distributions, but this control is less apparent in the current data. We suggest that due to the presence of thresholds, the catchment steepness, MAP, and sediment travel distance are non-linearly related to grain-size-dependent 10Be concentrations, which complicates efforts to disentangle and quantify their relative roles. The results of our study can be used to evaluate whether catchment-average denudation rates may be biased in particular catchments.

Data availability
Data availability.

All tables (S1–S5) and figures (S1–S6) in the Supplement are available in the data supplement of van Dongen et al. (2019): https://doi.org/10.5880/GFZ.3.3.2019.002. These data are freely available under the Creative Commons Attribution 4.0 International (CC BY 4.0) open-access license at GFZ data services. When using these data please cite this paper.

Sample availability
Sample availability.

The metadata of all samples in Table 2 can be accessed via http://igsn.org/.

Author contributions
Author contributions.

RvD carried out fieldwork, laboratory work, and data evaluation. DS conceived the study and was the main advisor during fieldwork, data evaluation, and article writing. HW was responsible for the cosmogenic-10Be laboratory training and laboratory supervision. FvB and HW were available for extensive discussion during data evaluation. RvD prepared the manuscript with contributions and edits from all co-authors.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We acknowledge support from the German Science Foundation (DFG) priority research program SPP-1803 “EarthShape: Earth Surface Shaping by Biota” (grant SCHE 1676/4-1 to Dirk Scherler). We are grateful to the Chilean National Park Service (CONAF) for providing access to the sample locations and on-site support of our research. We also thank Luca Mao, Ricardo Carrillo, and Melle Koelewijn for their support during fieldwork; Steven A. Binnie and Stefan Heinze from Cologne University for conducting AMS measurements; and Michael J. Henehan for his help with statistical analysis. We thank Amanda Schmidt, Stefanie Tofelde, and the two anonymous reviewers for their constructive reviews.

The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.

Review statement
Review statement.

This paper was edited by Michele Koppes and reviewed by Amanda Schmidt and two anonymous referees.

References

Aguilar, G., Carretier, S., Regard, V., Vassallo, R., Riquelme, R., and Martinod, J.: Grain size-dependent 10Be concentrations in alluvial stream sediment of the Huasco Valley, a semi-arid Andes region, Quat. Geochronol., 19, 163–172, https://doi.org/10.1016/j.quageo.2013.01.011, 2014.

Allen, P. A., Armitage, J. J., Whittaker, A. C., Michael, N. A., Roda-boluda, D., and Arcy, M. D.: Fragmentation Model of the Grain Size Mix of Sediment Supplied to Basins, J. Geol., 123, 405–427, https://doi.org/10.1086/683113, 2015.

Anderson, S. P., von Blanckenburg, F., and White, A.: Physical and Chemical Controls on the Chritical Zone, Elements, 3, 315–319, https://doi.org/10.2113/gselements.3.5.315, 2007.

Attal, M. and Lavé, J.: Pebble abrasion during fluvial transport: Experimental results and implications for the evolution of the sediment load along rivers, J. Geophys. Res.-Earth, 114, 1–22, https://doi.org/10.1029/2009JF001328, 2009.

Attal, M., Mudd, S. M., Hurst, M. D., Weinman, B., Yoo, K., and Naylor, M.: Impact of change in erosion rate and landscape steepness on hillslope and fluvial sediments grain size in the Feather River basin (Sierra Nevada, California), Earth Surf. Dynam., 3, 201–222, https://doi.org/10.5194/esurf-3-201-2015, 2015.

Balco, G., Stone, J. O., Lifton, N. A. and Dunai, T. J.: A complete and easily accessible means of calculating surface exposure ages or erosion rates from 10Be and 26Al measurements, Quat. Geochronol., 3, 174–195, https://doi.org/10.1016/j.quageo.2007.12.001, 2008.

Belmont, P., Pazzaglia, F. J., and Gosse, J. C.: Cosmogenic 10Be as a tracer for hillslope and channel sediment dynamics in the Clearwater River, western Washington State, Earth Planet. Sc. Lett., 264, 123–135, https://doi.org/10.1016/j.epsl.2007.09.013, 2007.

Bernhard, N., Moskwa, L. M., Schmidt, K., Oeser, R. A., Aburto, F., Bader, M. Y., Baumann, K., von Blanckenburg, F., Boy, J., van den Brink, L., Brucker, E., Büdel, B., Canessa, R., Dippold, M. A., Ehlers, T. A., Fuentes, J. P., Godoy, R., Jung, P., Karsten, U., Köster, M., Kuzyakov, Y., Leinweber, P., Neidhardt, H., Matus, F., Mueller, C. W., Oelmann, Y., Oses, R., Osses, P., Paulino, L., Samolov, E., Schaller, M., Schmid, M., Spielvogel, S., Spohn, M., Stock, S., Stroncik, N., Tielbörger, K., Übernickel, K., Scholten, T., Seguel, O., Wagner, D., and Kühn, P.: Pedogenic and microbial interrelations to regional climate and local topography: new insights from a climate gradient (arid to humid) along the Coastal Cordillera of Chile, Catena, 170, 335–355, https://doi.org/10.1016/j.catena.2018.06.018, 2018.

Bierman, P. R. and Steig, E. J.: Estimating rates of denudation using cosmogenic isotope abundances in sediment, Earth Surf. Proc. Land., 21, 125–139, https://doi.org/10.1002/(SICI)1096-9837(199602)21:2<125::AID-ESP511>3.0.CO;2-8, 1996.

Binnie, S. A., Phillips, W. M., Summerfield, M. A., and Keith Fifield, L.: Sediment mixing and basin-wide cosmogenic nuclide analysis in rapidly eroding mountainous environments, Quat. Geochronol., 1, 4–14, https://doi.org/10.1016/j.quageo.2006.06.013, 2006.

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, https://doi.org/10.1130/G23641A.1, 2007.

Borchers, B., Marrero, S., Balco, G., Caffee, M., Goehring, B., Lifton, N., Nishiizumi, K., Phillips, F., Schaefer, J., and Stone, J.: Geological calibration of spallation production rates in the CRONUS- Earth project, Quat. Geochronol., 31, 188–198, https://doi.org/10.1016/j.quageo.2015.01.009, 2016.

Brown, E. T., Stallard, R. F., Larsen, M. C., Raisbeck, G. M., and Yiou, F.: Denudation rates determined from the accumulation of in situ-produced 10Be in the luquillo experimental forest, Puerto Rico, Earth Planet. Sc. Lett., 129, 193–202, https://doi.org/10.1016/0012-821X(94)00249-X, 1995.

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.

Carretier, S. and Regard, V.: Is it possible to quantify pebble abrasion and velocity in rivers using terrestrial cosmogenic nuclides?, J. Geophys. Res.-Earth, 116, 1–17, https://doi.org/10.1029/2011JF001968, 2011.

Carretier, S., Regard, V., and Soual, C.: Theoretical cosmogenic nuclide concentration in river bed load clasts: Does it depend on clast size?, Quat. Geochronol., 4, 108–123, https://doi.org/10.1016/j.quageo.2008.11.004, 2009.

Carretier, S., Regard, V., Vassallo, R., Aguilar, G., Martinod, J., Riquelme, R., Christophoul, F., Charrier, R., Gayer, E., Farías, M., Audin, L., and Lagane, C.: Differences in 10Be concentrations between river sand, gravel and pebbles along the western side of the central Andes, Quat. Geochronol., 27, 33–51, https://doi.org/10.1016/j.quageo.2014.12.002, 2015.

Casagli, N., Ermini, L., and Rosati, G.: Determining grain size distribution of the material composing landslide dams in the Northern Apennines: Sampling and processing methods, Eng. Geol., 69, 83–97, https://doi.org/10.1016/S0013-7952(02)00249-1, 2003.

Chang, K. T., Chiang, S. H., and Hsu, M. L.: Modeling typhoon- and earthquake-induced landslides in a mountainous watershed using logistic regression, Geomorphology, 89, 335–347, https://doi.org/10.1016/j.geomorph.2006.12.011, 2007.

Chen, H., Dadson, S., and Chi, Y. G.: Recent rainfall-induced landslides and debris flow in northern Taiwan, Geomorphology, 77, 112–125, https://doi.org/10.1016/j.geomorph.2006.01.002, 2006.

Clapp, E. M., Bierman, P. R., and Caffee, M.: Using 10Be and 26Al to determine sediment generation rates and identify sediment source areas in an arid region drainage basin, Geomorphology, 45, 89–104, https://doi.org/10.1016/S0169-555X(01)00191-X, 2002.

Clarke, B. A. and Burbank, D. W.: Quantifying bedrock – fracture patterns within the shallow subsurface?: Implications for rock mass strength, bedrock landslides, and erodibility, J. Geophys. Res., 116, F04009, https://doi.org/10.1029/2011JF001987, 2011.

Codilean, A. T., Munack, H., Cohen, T. J., Saktura, W. M., Gray, A., and Mudd, S. M.: OCTOPUS: an open cosmogenic isotope and luminescence database, Earth Syst. Sci. Data, 10, 2123–2139, https://doi.org/10.5194/essd-10-2123-2018, 2018.

Crosby, B. T. and Whipple, K. X.: Knickpoint initiation and distribution within fluvial networks: 236 waterfalls in the Waipaoa River, North Island, New Zealand, Geomorphology, 82, 16–38, https://doi.org/10.1016/j.geomorph.2005.08.023, 2006.

Derrieux, F., Siame, L. L., Bourlès, D. L., Chen, R. F., Braucher, R., Léanni, L., Lee, J. C., Chu, H. T., and Byrne, T. B.: How fast is the denudation of the Taiwan mountain belt? Perspectives from in situ cosmogenic 10Be, J. Asian Earth Sci., 88, 230–245, https://doi.org/10.1016/j.jseaes.2014.03.012, 2014.

DiBiase, R. A. and Whipple, K. X.: The influence of erosion thresholds and runoff variability on the relationships among topography, climate, and erosion rate, J. Geophys. Res., 116, F04036, https://doi.org/10.1029/2011JF002095, 2011.

DiBiase, R. A., Whipple, K. X., Heimsath, A. M., and Ouimet, W. B.: Landscape form and millenial erosion rates in the San Gabriel Mountains, Ca, Earth Planet. Sc. Lett., 289, 134–144, https://doi.org/10.1016/j.epsl.2009.10.036, 2010.

Didan, K.: MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid V006 [Data set], https://doi.org/10.5067/MODIS/MOD13Q1.006, 2015.

Dixon, J. L. and von Blanckenburg, F.: Soils as pacemakers and limiters of global silicate weathering, C. R. Geosci., 344, 597–609, https://doi.org/10.1016/j.crte.2012.10.012, 2012.

Dixon, J. L., Chadwick, O. A., and Vitousek, P. M.: Climate-driven threshold for chemical weathering in postglacial soils of New Zealand, J. Geophys. Res.-Earth, 121, 1619–1634, https://doi.org/10.1002/2016JF003864, 2016.

Drever, J. I.: The effect of land plants on weathering rates of silicate minerals, Geochim. Cosmochim. Ac., 58, 2325–2332, https://doi.org/10.1016/0016-7037(94)90013-2, 1994.

Ehrlich, H. L.: Geomicrobiology: Its significance for geology, Earth Sci. Rev., 45, 45–60, https://doi.org/10.1016/S0012-8252(98)00034-8, 1998.

Ferguson, R., Hoey, T., Wathen, S., and Werritty, A.: Field evidence for rapid downstream fining of river gravels through selective transport, Geology, 24, 179–182, https://doi.org/10.1130/0091-7613(1996)024<0179:FEFRDF>2.3.CO;2, 1996.

Gabet, E. J. and Dunne, T.: Landslides on coastal sage-scrub and grassland hillslopes in a severe El Niño winter: The effects of vegetation conversion on sediment delivery, B. Geol. Soc. Am., 114, 983–990, https://doi.org/10.1130/0016-7606(2002)114<0983:LOCSSA>2.0.CO;2, 2002.

Gabet, E. J. and Mudd, S. M.: Bedrock erosion by root fracture and tree throw: A coupled biogeomorphic model to explore the humped soil production function and the persistence of hillslope soils, J. Geophys. Res.-Earth, 115, 1–14, https://doi.org/10.1029/2009JF001526, 2010.

Gabet, E. J., Reichman, O. J., and Seabloom, E. W.: The Effects of Bioturbation on Soil Processes and Sediment Transport, Annu. Rev. Earth Pl. Sc., 31, 249–273, https://doi.org/10.1146/annurev.earth.31.100901.141314, 2003.

Glodny, J., Graaefe, K., and Rosenau, M.: Mesozoic to Quaternary continental margin dynamics in South-Central Chile (36–42  S): the apatite and zircon fission track perspective, Int. J. Earth Sci., 97, 1271–1291, https://doi.org/10.1007/s00531-007-0203-1, 2008.

Gonzalez, V. S., Schmidt, A. H., Bierman, P. R., and Rood, D. H.: Spatial and temporal replicability of meteoric and in situ 10Be concentrations in fluvial sediment, Earth Surf. Proc. Land., 42, 2570–2584, https://doi.org/10.1002/esp.4205, 2017.

Gosse, J. C. and Phillips, F. M.: Terrestrial in situ cosmogenic nuclides:theory and application, Quaternary Sci. Rev., 20, 1475–1560, https://doi.org/10.1016/S0277-3791(00)00171-2, 2001.

Granger, D., Kirchner, J., and Finkel, R.: Spatially averaged long-term erosion rates measured from in situ-produced cosmogenic nuclides in alluvial sediment, J. Geol., 104, 249–257, https://doi.org/10.1086/629823, 1996.

Grömping, U.: R package relaimpo: relative importance for linear regression, J. Stat. Softw., 17, 139–147, https://doi.org/10.18637/jss.v017.i01, 2006.

Hartmann, J. and Moosdorf, N.: The new global lithological map database GLiM: A representation of rock properties at the Earth surface, Geochem. Geophy. Geosy., 13, 1–37, https://doi.org/10.1029/2012GC004370, 2012.

Heimsath, A. M., Fink, D., and Hancock, G. R.: The “humped” soil production function?: eroding Arnhem Land, Australia, Earth Surf. Proc. Land., 1684, 1674–1684, https://doi.org/10.1002/esp.1859, 2009.

Hoey, T. B. and Ferguson, R.: Numerical simulation of downstream fining by selective transport in gravel bed rivers: Model development and illustration, Water Resour., 30, 2251–2260, https://doi.org/10.1029/94WR00556, 1994.

Hovius, N., Stark, C. A., and Allen, P. A.: Sediment flux from a mountain belt derived by landsliding mapping, Geology, 25, 231–234, https://doi.org/10.1130/0091-7613(1997)025<0231:SFFAMB>2.3.CO;2, 1997.

Jarvis, A., Reuter, H. I., Nelson, A., and Guevara, E.: Hole-filled SRTM for the globe Version 4, available from the CGIAR-CSI SRTM 90m Database, available at: http://srtm.csi.cgiar.org (last access: 15 October 2018), 2008.

Kodama, Y: Downstream changes in the lithology and grain size of fluvial gravels, the Watarase River, Japan; evidence of the role of abrasion in downstream fining. J. Sediment. Res., 64, 68–75, https://doi.org/10.1306/D4267D0C-2B26-11D7-8648000102C1865D, 1994a.

Kodama, Y: Experimental Study of Abrasion and Its Role in Producing Downstream Fining in Gravel-Bed Rivers, J. Sediment. Res., 64, 76–85, https://doi.org/10.2110/jsr.64.76, 1994b.

Kohl, C. P. and Nishiizumi, K.: Chemical isolation of quartz for measurement of in-situ -produced cosmogenic nuclides, Geochim. Cosmochim. Ac., 56, 3583–3587, https://doi.org/10.1016/0016-7037(92)90401-4, 1992.

Kolmogorov, A.: Sulla determinazione empirica di una legge di distribuzione, G. dell'Instituto Ital. degli Attuari, 4, 83–91, 1933.

Krumbein, W. C.: Size frequency distributions of sediments, J. Sediment. Petrol., 4, 65–77, 1934.

Krumbein, W. C.: Size frequency distributions of sediments and the normal phi curve, J. Sediment. Petrol., 8, 84–90, 1938.

Lal, D.: Cosmic ray labeling of erosion surfaces?: in situ nuclide production rates and erosion models, Earth Planet. Sc. Lett., 104, 424–439, 1991.

Larsen, I. J. and Montgomery, D. R.: Landslide erosion coupled to tectonics and river incision, Nat. Geosci., 5, 468–473, https://doi.org/10.1038/ngeo1479, 2012.

Lasaga, A. C., Soler, J. M., Ganor, J., Burch, T. E., and Nagy, K. L.: Chemical weathering rate laws and global geochemical cycles, Geochim. Cosmochim. Ac., 58, 2361–2386, https://doi.org/10.1016/0016-7037(94)90016-7, 1994.

Lebedeva, M. I. and Brantley, S. L.: Weathering and erosion of fractured bedrock systems, Earth Surf. Proc. Land., 42, 2090–2108, https://doi.org/10.1002/esp.4177, 2017.

Lewin, J. and Brewer, P. A.: Laboratory Simulation of Clast Abrasion, Earth Surf. Proc. Land., 27, 145–164, https://doi.org/10.1002/esp.306, 2002.

Lin, G. W., Chen, H., Chen, Y. H., and Horng, M. J.: Influence of typhoons and earthquakes on rainfall-induced landslides and suspended sediments discharge, Eng. Geol., 97, 32–41, https://doi.org/10.1016/j.enggeo.2007.12.001, 2008.

Lindeman, R. H., Merenda, P. F., and Gold, R. Z.: Introduction to Bivariate and Multivariate Analysis, Scott, Foresman, Glenview, IL, 1, 1980.

Lukens, C. E., Riebe, C. S., Sklar, L. S., and Shuster, D. L.: Grain size bias in cosmogenic nuclide studies of stream sediment in steep terrain, J. Geophys. Res.-Earth, 121, 978–999, https://doi.org/10.1002/2016JF003859, 2016.

Lupker, M., Lavé, J., France-Lanord, C., Christl, M., Bourlès, D., Carcaillet, J., Maden, C., Wieler, R., Rahman, M., Bezbaruah, D., and Xiaohan, L.: 10Be systematics in the Tsangpo-Brahmaputra catchment: the cosmogenic nuclide legacy of the eastern Himalayan syntaxis, Earth Surf. Dynam., 5, 429–449, https://doi.org/10.5194/esurf-5-429-2017, 2017.

Maher, K.: The dependence of chemical weathering rates on fluid residence time, Earth Planet. Sc. Lett., 294, 101–110, https://doi.org/10.1016/j.epsl.2010.03.010, 2010.

Mather, A. E., Hartley, A. J., and Griffiths, J. S.: The giant coastal landslides of Northern Chile: Tectonic and climate interactions on a classic convergent plate margin, Earth Planet. Sc. Lett., 388, 249–256, https://doi.org/10.1016/j.epsl.2013.10.019, 2014.

Matmon, A., Bierman, P. R., Larsen, J., Southworth, S., Pavich, M., Finkel, R., and Caffee, M.: Erosion of an ancient mountain range, the Great Smoky Mountains, North Carolina and Tennessee, Am. J. Sci., 303, 817–855, https://doi.org/10.2475/ajs.303.9.817, 2003.

Melnick, D.: Rise of the central Andean coast by earthquakes straddling the Moho, Nat. Geosci., 9, 1–8, https://doi.org/10.1038/ngeo2683, 2016.

Melnick, D., Bookhagen, B., Strecker, M. R., and Echtler, H. P.: Segmentation of megathrust rupture zones from fore-arc deformation patterns over hundreds to millions of years, Arauco peninsula, Chile, J. Geophys. Res.-Sol. Ea., 114, 1–23, https://doi.org/10.1029/2008JB005788, 2009.

Meyer-Christoffer, A., Becker, A., Finger, P., Rudolf, B., Schneider, U., and Ziese, M.: GPCC Climatology Version 2015 at 0.25: Monthly Land-Surface Precipitation Climatology for Every Month and the Total Year from Rain-Gauges built on GTS-based and Historic Data, https://doi.org/10.5676/DWD_GPCC/CLIM_M_V2015_025, 2015.

Montgomery, D. R. and Brandon, M. T.: Topographic controls on erosion rates in tectonically active mountain ranges, Earth Planet. Sc. Lett., 201, 481–489, https://doi.org/10.1016/S0012-821X(02)00725-2, 2002.

Neilson, T. B., Schmidt, A. H., Bierman, P. R., Rood, D. H., and Sosa Gonzalez, V.: Efficacy of in situ and meteoric 10Be mixing in fluvial sediment collected from small catchments in China, Chem. Geol., 471, 119–130, https://doi.org/10.1016/j.chemgeo.2017.09.024, 2017.

Niemi, N. A., Oskin, M., Burbank, D. W., Heimsath, A. M., and Gabet, E. J.: Effects of bedrock landslides on cosmogenically determined erosion rates, Earth Planet. Sc. Lett., 237, 480–498, https://doi.org/10.1016/j.epsl.2005.07.009, 2005.

Oberlander, T. M.: Morphogenesis of Granitic Boulder Slopes in the Mojave Desert, California, J. Geol., 80, 1–20, 1972.

Oeser, R. A., Stroncik, N., Moskwa, L., Bernhard, N., Schaller, M., Canessa, R., Brink, L. Van Den, Köster, M., Brucker, E., Stock, S., Pablo, J., Godoy, R., Javier, F., Oses, R., Osses, P., Paulino, L., Seguel, O., Bader, M. Y., Boy, J., Dippold, M. A., Ehlers, T. A., Kühn, P., Kuzyakov, Y., Leinweber, P., Scholten, T., Spielvogel, S., Spohn, M., Übernickel, K., Tielbörger, K., Wagner, D., and Blanckenburg, F. Von: Chemistry and microbiology of the Critical Zone along a steep climate and vegetation gradient in the Chilean Coastal Cordillera, Catena, 170, 183–203, https://doi.org/10.1016/j.catena.2018.06.002, 2018.

Ouimet, W. B., Whipple, K. X., and Granger, D. E.: Beyond threshold hillslopes: Channel adjustment to base-level fall in tectonically active mountain ranges, Geology, 37, 579–582, https://doi.org/10.1130/G30013A.1, 2009.

Palumbo, L., Hetzel, R., Tao, M., and Li, X.: Topographic and lithologic control on catchment-wide denudation rates derived from cosmogenic 10Be in two mountain ranges at the margin of NE Tibet, Geomorphology, 117, 130–142, https://doi.org/10.1016/j.geomorph.2009.11.019, 2010.

Perras, M. A. and Diederichs, M. S.: A Review of the Tensile Strength of Rock: Concepts and Testing, Geotech. Geol. Eng., 32, 525–546, https://doi.org/10.1007/s10706-014-9732-0, 2014.

Phillips, F. M., Argento, D. C., Balco, G., Caffee, M. W., Clem, J., Dunai, T. J., Finkel, R., Goehring, B., Gosse, J. C., Hudson, A. M., Jull, A. J. T., Kelly, M. A., Kurz, M., Lal, D., Lifton, N., Marrero, S. M., Nishiizumi, K., Reedy, R. C., Schaefer, J., Stone, J. O. H., Swanson, T., and Zreda, M. G.: The CRONUS-Earth Project: A synthesis, Quat. Geochronol., 31, 119–154, https://doi.org/10.1016/j.quageo.2015.09.006, 2016.

Pinto, L., Hérail, G., Sepúlveda, S. A., and Krop, P.: A Neogene giant landslide in Tarapacá, northern Chile: A signal of instability of the westernmost Altiplano and palaeoseismicity effects, Geomorphology, 102, 532–541, https://doi.org/10.1016/j.geomorph.2008.05.044, 2008.

Portenga, E. W. and Bierman, P. R.: Understanding earth's eroding surface with 10Be, GSA Today, 21, 4–10, https://doi.org/10.1130/G111A.1, 2011.

Puchol, N., Lavé, J., Lupker, M., Blard, P. H., Gallo, F., and France-Lanord, C.: Grain-size dependent concentration of cosmogenic 10Be and erosion dynamics in a landslide-dominated Himalayan watershed, Geomorphology, 224, 55–68, https://doi.org/10.1016/j.geomorph.2014.06.019, 2014.

Reinhardt, L. J., Hoey, T. B., Barrows, T. T., Dempster, T. J., Bishop, P., and Fifield, L. K.: Interpreting erosion rates from cosmogenic radionuclide concentrations measured in rapidly eroding terrain, Earth Surf. Proc. Land., 32, 390–406, https://doi.org/10.1002/esp.1415, 2007.

Riebe, C. S., Kirchner, J. W., and Finkel, R. C.: Long-term rates of chemical weathering and physical erosion from cosmogenic nuclides and geochemical mass balance, Geochim. Cosmochim. Ac., 67, 4411–4427, https://doi.org/10.1016/S0016-7037(03)00382-X, 2003.

Riebe, C. S., Kirchner, J. W., and Finkel, R. C.: Erosional and climatic effects on long-term chemical weathering rates in granitic landscapes spanning diverse climate regimes, Earth Planet. Sc. Lett., 224, 547–562, https://doi.org/10.1016/j.epsl.2004.05.019, 2004.

Roda-Boluda, D. C., D'Arcy, M., McDonald, J., and Whittaker, A. C.: Lithological controls on hillslope sediment supply: insights from landslide activity and grain size distributions, Earth Surf. Proc. Land., 43, 956–977, https://doi.org/10.1002/esp.4281, 2018.

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, https://doi.org/10.1029/1998WR900090, 1999.

Roering, J. J., Marshall, J., Booth, A. M., Mort, M., and Jin, Q.: Evidence for biotic controls on topography and soil production, Earth Planet. Sc. Lett., 298, 183–190, https://doi.org/10.1016/j.epsl.2010.07.040, 2010.

Ruxton, B. P. and Berry, L.: Weathering of granite and associated erosional features in Hong Kong, Bull. Geol. Soc. Am., 68, 1263–1292, https://doi.org/10.1130/0016-7606(1957)68[1263:WOGAAE]2.0.CO;2, 1957.

Safran, E. B., Bierman, P. R., Aalto, R., Dunne, T., Whipple, K. X., and Caffee, M.: Erosion rates driven by channel network incision in the Bolivian Andes. Earth Surface Processes and Landforms, 30, 1007–1024, https://doi.org/10.1002/esp.1259, 2005.

Schaller, M., Ehlers, T. A., Lang, K. A. H., Schmid, M., and Fuentes-Espoz, J. P.: Addressing the contribution of climate and vegetation cover on hillslope denudation, Chilean Coastal Cordillera (26–38S), Earth Planet. Sc. Lett., 489, 111–122, https://doi.org/10.1016/j.epsl.2018.02.026, 2018.

Scherler, D., Lamb, M. P., Rhodes, E. J., and Avouac, J. P.: Climate-change versus landslide origin of fill terraces in a rapidly eroding bedrock landscape: San Gabriel River, California, Bull. Geol. Soc. Am., 128, 1228–1248, https://doi.org/10.1130/B31356.1, 2016.

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.

Sklar, L. S. and Dietrich, W. E.: Sediment and rock strength controls on river incision into bedrock, Geology, 29(12), 1087–1090, https://doi.org/10.1130/0091-7613(2001)029<1087:SARSCO>2.0.CO;2, 2001.

Sklar, L. S. and Dietrich, W. E.: Implications of the Saltation-Abrasion Bedrock Incision Model for Steady-State River Longitudinal Profile Relief and Concavity, Earth Surf. Proc. Land., 33, 1129–1151, https://doi.org/10.1002/esp.1689, 2008.

Sklar, L. S., Dietrich, W. E., Foufoula-Georgiou, E., Lashermes, B., and Bellugi, D.: Do gravel bed river size distributions record channel network structure?, Water Resour. Res., 42, 1–22, https://doi.org/10.1029/2006WR005035, 2006.

Sklar, L. S., Riebe, C. S., Marshall, J. A., Genetti, J., Leclere, S., Lukens, C. L., and Merces, V.: The problem of predicting the size distribution of sediment supplied by hillslopes to rivers, Geomorphology, 277, 31–49, https://doi.org/10.1016/j.geomorph.2016.05.005, 2017.

Smirnov, N. V.: On the estimation of the discrepancy between empirical curves of distribution for two independent samples, Bull. Math. Univ. Moscou, 2, 3–16, 1939.

Sosa Gonzalez, V., Bierman, P. R., Fernandes, N. F., and Rood, D. H.: Long-term background denudation rates of southern and southeastern Brazilian watersheds estimated with cosmogenic 10 Be, Geomorphology, 268, 54–63, 2016a.

Sosa Gonzalez, V., Bierman, P. R., Nichols, K. K., and Rood, D. H.: Long-term erosion rates of Panamanian drainage basins determined using in situ 10 Be, Geomorphol. Tech., 275, 1–15, 2016b.

Stock, G. M., Frankel, K. L., Ehlers, T. a., Schaller, M., Briggs, S. M., and Finkel, R. C.: Spatial and temporal variations in denudation of the Wasatch Mountains, Utah, USA, Lithosphere, 1, 34–40, https://doi.org/10.1130/L15.1, 2009.

Stone, J. O.: Air pressure and cosmogenic isotope production, J. Geophys. Res., 105, 753–759, https://doi.org/10.1029/2000JB900181, 2000.

Sullivan, C. L.: 10Be Erosion Rates and Landscape Evolution of the Blue Ridge Escarpment, Southern Appalachian Mountains, MS thesis, The University of Vermont, 2007.

Tofelde, S., Duesing, W., Schildgen, T. F., Wickert, A. D., Wittmann, H., Alonso, R. N., and Strecker, M.: Effects of deep-seated versus shallow hillslope processes on cosmogenic10Be concentrations in fluvial sand and gravel, Earth Surf. Proc. Land., 39, 3086–3098, https://doi.org/10.1002/esp.4471, 2018.

von Blanckenburg, F.: The control mechanisms of erosion and weathering at basin scale from cosmogenic nuclides in river sediment, Earth Planet. Sc. Lett., 237, 462–479, https://doi.org/10.1016/j.epsl.2005.06.030, 2005.

von Blanckenburg, F., Hewawasam, T., and Kubik, P. W.: Cosmogenic nuclide evidence for low weathering and denudation in the wet, tropical highlands of Sri Lanka, J. Geophys. Res., 109, F03008, https://doi.org/10.1029/2003JF000049, 2004.

van Dongen, R., Scherler, D., Wittmann, H., and von Blanckenburg, F.: Data supplement to: Cosmogenic 10Be in river sediment: where grain size matters and why, GFZ Ger. Res. Cent. Geosci., https://doi.org/10.5880/GFZ.3.3.2019.002, 2019.

West, A. J., Hetzel, R., Li, G., Jin, Z., Zhang, F., Hilton, R. G., and Densmore, A. L.: Dilution of 10 Be in detrital quartz by earthquake-induced landslides?: Implications for determining denudation rates and potential to provide insights into landslide sediment dynamics, Earth Planet. Sc. Lett., 396, 143–153, https://doi.org/10.1016/j.epsl.2014.03.058, 2014.

Whipple, K. X. and Tucker, G. E.: Dynamics of the stream-power river incision model: Implications for height limits of mountain ranges, landscape response timescales and research needs, J. Geophys. Res., 104, 661–674, 1999.

White, A. F., Blum, A. E., Bullen, T. D., Vivit, D. V., Schulz, M., and Fitzpatrick, J.: The effect of temperature on experimental and natural chemical weathering rates of granitoid rocks, Geochim. Cosmochim. Ac., 63, 3277–3291, 1999.

Willenbring, J. K., Gasparini, N. M., Crosby, B. T., and Brocard, G.: What does a mean mean? The temporal evolution of detrital cosmogenic denudation rates in a transient landscape, Geology, 41, 1215–1218, https://doi.org/10.1130/G34746.1, 2013.

Wolman, M. G.: A method of sampling coarse river-bed material, Trans. Am. Geophys. Union, 35, 951–956, https://doi.org/10.1029/TR035i006p00951, 1954.

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.-Earth, 114, 114, F01007, https://doi.org/10.1029/2008JF001088, 2009.