Articles | Volume 10, issue 5
Earth Surf. Dynam., 10, 909–928, 2022
Earth Surf. Dynam., 10, 909–928, 2022
Research article
19 Sep 2022
Research article | 19 Sep 2022

Quantification of post-glacier bedrock surface erosion in the European Alps using 10Be and optically stimulated luminescence exposure dating

Quantification of post-glacier bedrock surface erosion in the European Alps using 10Be and optically stimulated luminescence exposure dating
Joanne Elkadi1, Benjamin Lehmann2, Georgina E. King1, Olivia Steinemann3, Susan Ivy-Ochs3, Marcus Christl3, and Frédéric Herman1 Joanne Elkadi et al.
  • 1Institute of Earth Surface Dynamics, University of Lausanne, 1015 Lausanne, Switzerland
  • 2INSTAAR and Department of Geological Sciences, University of Colorado Boulder, Boulder, CO 80309, USA
  • 3Laboratory of Ion Beam Physics, ETH Zürich, Otto-Stern-Weg 5, 8093 Zurich, Switzerland

Correspondence: Joanne Elkadi (


The retreat of glaciers since the Last Glacial Maximum in the European Alps has left an imprint on topography through various erosional processes. However, few methods are currently capable of resolving these mechanisms on Late Glacial to Holocene timescales. Quantifying the relative contributions of mountain erosion, during these different climate cycles, is useful for understanding long-term landscape evolution and the links between global climate and erosion. Here, we combine three optically stimulated luminescence (OSL) exposure dating signals with 10Be surface exposure dating to constrain the post-glacier erosion rates of bedrock samples adjacent to the Gorner Glacier in the European Alps. The results reveal erosion rates of the order of 10−2 to 10−1mm a−1, in general agreement with other studies in the region, as well as a strong negative correlation between erosion rate and elevation, suggesting that frost crack weathering is perhaps not the dominant form of post-glacier weathering. Finally, a global compilation of both subglacial and periglacial erosion rates shows that periglacial erosion rates could be greater than previously thought. Yet subglacial erosion remains higher, implying that it continues to have the stronger influence on shaping landscapes. Therefore, with a changing climate, periglacial erosion rates are likely to remain transient. These insights could lead to important implications for landscape evolution models.

1 Introduction

The interplay between erosion and climate has sparked debates and inspired research aimed at better understanding the efficacy of various erosion mechanisms in long-term landscape evolution as well as the role that climate and its variability play in setting these erosion rates (e.g. Zhang et al., 2001; Molnar, 2004; Willenbring and von Blanckenburg, 2010; Lupker et al., 2013; Cogez et al., 2015; Herman et al., 2013; Herman and Champagnac, 2016; Willenbring and Jerolmack; 2016). Globally, continental topography has been shaped partly through erosional processes associated with rivers, glaciers, soils, rockfall and weathering. For high mountain environments specifically, the strong imprint of glacial and non-glacial erosion is observed at mid to high latitudes, but their specific mechanisms and respective impacts on the topography remain convoluted (e.g. André, 2002a; Ballantyne, 2002; Koppes and Montgomery, 2009). Here, non-glacial erosion refers broadly to any erosion occurring in a glacial environment that is not related to subglacial erosion. It is necessary to better quantify this to develop our knowledge of the influence of mountain erosion on the global feedback loop that exists between climate and erosion during glacial and interglacial times.

Currently in alpine environments, glacial erosion and its associated processes are thought to play a dominant role, and thus extensive research has addressed its quantification as well as the timing of deglaciations (e.g. Hallet et al., 1996; Montgomery, 2002; Ivy-Ochs and Briner, 2014; Herman et al., 2015; 2018; Wirsig, 2016; Wirsig et al., 2016, 2017; Ruszkiczay-Rüdiger et al., 2021; Steinemann et al., 2021). In contrast, studies exploring periglacial erosion or erosion during interglacial times have mainly investigated local fluvial incision (e.g. Korup and Schlunegger, 2007; Valla et al., 2010; Rolland et al., 2017) or catchment-wide erosion rates. A small number of studies have successfully quantified periglacial and interglacial erosion rates from bedrock surfaces through novel techniques (e.g. Kirkbride and Bell, 2010; Sohbati et al., 2018; Smedley et al., 2021), which is discussed in further detail in Sect. 1.1. Nevertheless, disentangling the relative contributions of the various erosional processes remains challenging (e.g. Hallet et al., 1996; Delmas et al., 2009; O'Farrell et al., 2009; Guillon et al., 2015; Cook et al., 2020).

We extend this dataset by applying the recently developed approach from Lehmann et al. (2019), which combines beryllium-10 (10Be) terrestrial cosmogenic nuclide (TCN) dating with optically stimulated luminescence (OSL) surface exposure dating to investigate bedrock post-glacier erosion rates (i.e. erosion since glacier retreat). This was done for six samples down a vertical transect adjacent to the Gorner Glacier near Zermatt, Switzerland. Then, we examine any potential trends between elevation or slope with our erosion rate results and find a strong negative correlation between erosion rate and elevation but no correlation between erosion rate and surface slope. Finally, the post-glacier erosion rates from this study are combined with global studies of both subglacial and periglacial erosion rates to reveal that periglacial erosion rates could be more comparable with subglacial erosion rates than anticipated.

1.1 Measuring erosion rates in deglaciated environments

At present, there exists a wide range of analytical techniques capable of quantifying bedrock erosion rates across different time intervals (please refer to Moses et al., 2014, and Turowski and Cook, 2017, for in-depth reviews). For timescales of the order of seconds to decades, these methods can include remote sensing (e.g. photogrammetry on both small and large spatial scales; Inkpen et al., 2000; Dornbusch et al., 2008) or measurements relative to anthropogenic reference points (e.g. lettering on gravestones; Inkpen and Jackson, 2000). On the other hand, studies targeting longer timescales (> 103 years) have measured relative to natural reference points (e.g. resistant quartz veins; Dahl, 1967; André, 2002b; Nicholson, 2008), exploited the half-lives of different cosmogenic nuclides (e.g. Nishiizumi et al., 1986; Bierman and Caffee, 2002; Balco et al., 2008) or used thermochronometry (e.g. Reiners and Brandon, 2006; Herman and King, 2018). Unfortunately, there is a lack of available methodological approaches to quantify bedrock erosion rates across the intermediate time interval, which has recently driven focussed research to devise new methods capable of doing so (Sohbati et al., 2018; Brown and Moon, 2019; Lehmann et al., 2019).

A small number of studies worldwide have already attempted to calculate periglacial rock surface erosion rates and have yielded a wide range of results. These include an investigation using TCN in the western US mountain ranges that estimated the maximum surface erosion rates of alpine bedrock summits at 7.6 × 10−3mm a−1 (Small et al., 1997), in contrast to another TCN study in the Nepal high Himalayas which instead found erosion rates of 8 × 10−2–2 × 10−1mm a−1 (Heimsath and McGlynn, 2007). In Europe, using reference quartz veins in Norway found erosion rates of 5.5 × 10−4mm a−1 (Nicholson, 2008), while incorporating the edge roundness of boulders in Scotland produced erosion rates of 3.3 × 10−3mm a−1 (Kirkbride and Bell, 2010). Recent studies combining TCN and OSL surface exposure dating in the Eastern Pamirs, China, and the Mont Blanc Massif, France, revealed bedrock surface erosion rates of < 3.8 × 10−5 and 1.72 × 10−3mm a−1 (Sohbati et al., 2018) and 3.53 × 10−3–4.3 mm a−1 (Lehmann et al., 2019, 2020), respectively. Smedley et al. (2021) also combined TCN and OSL surface exposure dating in north-western Scotland to derive interglacial erosion rates over the last 4 ka that were consistent with local independent erosion rate estimates (Kirkbride and Bell, 2010) and further inferred that some of their results could be explained by climatic fluctuations that are known to have occurred over that time period. Finally, a global compilation calculated by Portenga and Bierman (2011) gave an erosion rate of 1.2 × 10−2mm a−1 by averaging the results from studies that applied 10Be to bedrock surfaces. Here, we use a newly developed approach (Lehmann et al., 2019) that combines two surface exposure dating methods – 10Be and OSL – to investigate bedrock post-glacier erosion rates and onset times. In this case, the definition of erosion will be the removal of bedrock surface material. Bedrock surfaces offer great potential for the quantification of post-glacial erosion as (1) they are almost instantaneously exposed to the atmosphere once the ice retreats and thus immediately begin to record any changes to the surface associated with post-glacial erosion and (2) they are relatively durable in nature, rendering them capable of recording long-term erosional histories.

TCNs are formed at or near the Earth's surface within specific target minerals as a result of the Earth's constant bombardment by high-energy cosmic rays (Dunai, 2010; Gosse and Phillips, 2001). Consequently, following exposure, the concentration of nuclides measured in bedrock can be converted into an apparent exposure age. In this study, we focus on measuring 10Be which is found in quartz. In contrast, OSL is a trapped charge dating technique where a mineral, such as quartz or feldspar, emits light upon light stimulation due to electrons trapped in defects in the mineral's crystal lattice (Huntley et al., 1989; Aitken, 1998). The intensity of the light emitted is an indication of the concentration of trapped electrons. In recent years, the application of OSL to rock surface exposure dating has proven successful in a variety of settings (e.g. Sohbati et al., 2015; Lehmann et al., 2018; Liu et al., 2019) and is based on the principle that, for an exposed surface, the Sun's energy is sufficient to naturally reduce the surface luminescence signal to zero (e.g. Sohbati et al., 2011, 2012). This phenomenon is termed “bleaching”. Due to the attenuation of light, this bleaching effect decreases exponentially with depth into the rock until it becomes negligible (Habermann et al., 2000; Polikreti et al., 2002, 2003; Laskaris and Liritzis, 2011). Nonetheless, studies have shown that this depth of bleaching increases with exposure time (Habermann et al., 2000; Polikreti et al., 2002, 2003; Laskaris and Liritzis, 2011; Sohbati et al., 2011, 2012; Lehmann et al., 2018; Gliganic et al., 2019; Sellwood et al., 2019) and, after calibration to account for rock-specific bleaching rates, this bleaching depth can be translated directly into an apparent exposure age for surfaces which have not been affected by erosion (e.g. Lehmann et al., 2018; Sohbati et al., 2018). Recent luminescence instrument developments (Lapp et al., 2015) have facilitated the measurement of rock slices without requiring further mineral separation. Multiple luminescence signals can now be measured from the same slice to obtain the maximum amount of information (e.g. Jenkins et al., 2018; Luo et al., 2018; Meyer et al., 2018; Elkadi et al., 2021; Smedley et al., 2021). Although both TCN and OSL surface exposure dating are influenced by exposure time, they are also affected by surface erosion, and if this is not accounted for, then it can lead to an underestimation of exposure ages (e.g. Gosse and Phillips, 2001; Lehmann et al., 2019, 2020). However, the two methods have different sensitivities to erosion (Sohbati et al., 2018; Lehmann et al., 2019) since TCNs are formed several metres below a rock surface (Lal, 1991), whereas bleaching fronts in OSL depth profiles are created only in the top millimetres to centimetres (Vafiadou et al., 2007; Sohbati et al., 2011, 2012; Freiesleben et al., 2015). By using these two techniques in conjunction, this difference in sensitivity can be exploited to calculate the surface erosion of bedrock.

1.2 Study area

The post-glacier erosion rates were calculated from the flanks of the Gorner Glacier, located near the village of Zermatt, Switzerland. This area was chosen due to its well-constrained glacial history, consisting of a rich collection of geological maps and aerial photos as well as human observations, but also as a result of its proximity to the only other study that applied this method in the western Alps (Lehmann et al., 2019, 2020), allowing for direct comparisons.

Figure 1Study area and sampling sites on (a) regional, (b) glacier basin and (c) local scales. Panel (b) further depicts the glacier retreat history obtained from the Global Land Ice Measurements from Space initiative (GLIMS) (GLIMS Consortium, 2005a–f; Raup et al., 2007). The photo in the background is from the Swiss Federal Office of Topography and was taken in 2015. Panel (c) illustrates the difference in surface preservation with elevation and thus exposure age.

Table 1Summary of samples and measurements. All errors correspond to 1σ and encompass propagated uncertainties from the AMS measurements, blank correction and the local production rate.

a Estimated using old geological maps and aerial photos provided by the Swiss Federal Office of Topography. b The 10Be/9Be blank ratio used for the correction is 3.2 ± 1.7 × 10−15.

Download Print Version | Download XLSX

Bedrock material was collected from six sampling sites down a vertical transect, with sample lithologies consisting of hornfels, schist and gneiss (Table 1). Geomorphological reconstructions (Bini et al., 2009) suggest that, aside from the highest-elevation sample, the transect was covered in ice during the Last Glacial Maximum (LGM) and has been de-glaciated since. The three lower-elevation samples (GG04, GG05 and GG06) had additional exposure age information from old maps and photos acquired from the Swiss Federal Office of Topography. The three uppermost samples exhibited significant weathering, whereas the three lower-elevation samples had extremely well-preserved glacial morphologies and striations (Figs. 1 and S1 in the Supplement). Samples were collected from bedrock using a combination of a hammer, chisel and Husqvarna K760 power cutter with a diamond blade. Between two and four blocks with dimensions of  15 cm× 15 cm× 10 cm were extracted at each site, allowing for a sufficient amount of material for both OSL and 10Be surface exposure dating.

Table 2Protocol used for measuring the luminescence signals in the rock slices.

Value represents the wavelength that the emission is centred on.

Download Print Version | Download XLSX

2 Methodology

2.110Be sample preparation, measurement and age calculation

Sample preparation for 10Be dating was undertaken at ETH Zürich, Switzerland, and began with isolating quartz from the bulk rock. To do this, the uppermost 4 cm of each sample was crushed and sieved to obtain grain sizes between 250 and 1000 µm before being passed through a Frantz magnetic separator and subsequently treated with HCl and a low-concentration HF solution. Once pure quartz was isolated and dried, extraction of 10Be followed the well-established procedure outlined in Kohl and Nishiizumi (1992) and Ivy-Ochs (1996). First, the quartz grains were spiked with 0.25 mg of a 9Be carrier and dissolved with suprapure HF (40 %). Samples were then purified using two ion exchange resins to remove unwanted anions and cations followed by the final selective pH precipitation of Be(OH)2. 10Be/9Be ratios of the samples were measured using the 500 kV TANDY system at the accelerator mass spectrometry facility of ETH Zürich, Switzerland, using the in-house standard S2007N (Christl et al., 2013) calibrated against the 07KNSTD standard (Nishiizumi et al., 2007).

Surface exposure ages were determined using a modified version of the CREp online calculator (Martin et al., 2017; Lehmann et al., 2019) allowing for the application of a step-wise erosion rate correction. The calculations were done with blank-corrected 10Be/9Be ratios (full chemistry long-term laboratory procedural blank of 10Be/9Be (3.2 ± 1.7) × 10−15), the rescaled sea-level high-latitude (SLHL) Chironico landslide local production rate (4.16 ± 0.10 atomsgqtz-1a-1; Claude et al., 2014) corrected for each sample's latitude, longitude and elevation, the Lifton–Sato–Dunai (LSD) scaling scheme (Lifton et al., 2014), the ERA40 atmospheric model (Uppala et al., 2005) and the Lifton VDM 2016 geomagnetic database (Pavon-Carrasco et al., 2014, for ages between 0 and 14 ka; GLOPIS-75, Laj et al., 2004, for ages between 14 and 74 ka). Further information regarding the input data can be found in Table 1.

2.2 OSL

2.2.1 Sample preparation and measurement

Luminescence sample preparation followed the method described in Lehmann et al. (2018) and was carried out at the University of Lausanne, Switzerland. Following collection, samples were immediately placed in black, light-obstructing bags, and all ensuing laboratory work was done under subdued, red-light conditions. Laboratory work began with coring the samples using a water-cooled Husqvarna DM220 drill to extract multiple cores per sample, each with a diameter of 10 mm. Attention was given to drill as far from the edges as possible to avoid any potential signal resetting that may have occurred during fieldwork, and from areas with minimal lichen cover and red, iron-oxide staining that would have otherwise impeded light penetration and impacted the luminescence signal. The cores were then cut into thin slices ( 0.7 mm thickness) using a Buehler IsoMet low-speed saw mounted with a 0.3 mm-thick diamond blade and in the presence of a lubricant. The exact thickness of each slice was measured using a TESA Digitcal Caliper to allow precise reconstructions of the luminescence signal with depth. Example disc slices can be found in Fig. S2.

Table 3Symbols used in the luminescence surface exposure dating model.

Download Print Version | Download XLSX

Luminescence measurements followed the protocol outlined in Table 2, with low heating rates, extended isothermal holding times of samples prior to stimulation (1 Cs-1 and 100 s, respectively; Jenkins et al., 2018) and fragments of the slices placed in metal cups during measurement (Elkadi et al., 2021). Three signals were measured per sample – the IRSL50 and post-IR IRSL225 signals from feldspar and the post-IR IRSL225 OSL125 (hereinafter referred to as OSL125) signal predominantly from quartz. This was done to extract the maximum amount of information possible from each sample as the signals vary in bleaching rates. Since quartz and feldspar minerals are best stimulated using different wavelengths, different filters were required to isolate the luminescence emissions from the simulation wavelengths (Table 2). All measurements were performed using Risø TL-DA 20 TL-OSL readers (Bøtter-Jensen et al., 2010) equipped with a DASH head (Lapp et al., 2015) and a 90Sr/90Y beta source. The environmental dose rate (D˙) was calculated from U, Th, K and Rb concentrations of the bulk rock sample determined using ICPMS at Actlabs, Canada, and the DRAC online calculator (Durcan et al., 2015). All the luminescence signals were subsequently screened using three acceptance criteria: (1) maximum error of the test dose signal (Tn< 15 %, (2) Tn greater than 3σ above the background signal and (3) monotonic signal decay indicative of good heating (Elkadi et al., 2021). Any slices which did not meet these criteria were excluded from further analysis as their results were considered not reproducible.

2.2.2 Constraining the surface exposure dating model

The evolution of a luminescence signal L(x,t,r) (dimensionless) into a rock surface for a given depth x (mm), time t (year) and recombination centre distance r (dimensionless) can be modelled using the differential equation below (Lehmann et al., 2019):

(1) d L ( x , t , r ) d t = D ˙ D 0 [ 1 - L ( x , t , r ) ] - L ( x , t , r ) σ φ 0 e - μ x - L ( x , t , r ) s e - ρ - 1 3 r + ε ˙ ( t ) d L ( x , t , r ) d x .

This equation describes the four competing processes occurring following a surface's exposure to daylight: (1) electron trapping as a result of ambient radiation, (2) optical electron detrapping due to daylight exposure (bleaching), (3) athermal electron detrapping of the IRSL signal, most likely from quantum mechanical tunnelling in feldspars (Huntley, 2006; Kars et al., 2008), and (4) surface erosion. Definitions of the symbols can be found in Table 3, and we refer the reader to Lehmann et al. (2019) for further descriptions. A D0 value of 500 Gy (Lehmann et al., 2019) was selected for all samples following sensitivity tests that revealed the negligible effect of D0 on the final modelling results, even when varied by orders of magnitude.

The term describing optical electron detrapping contains two unknown parameters – σφ0 and μ – which have been shown to vary greatly across different lithologies, minerals and locations (e.g. Sohbati et al., 2012; Lehmann et al., 2018; Ou et al., 2018). Constraining the values of σφ0 and μ is one of the biggest challenges in luminescence surface exposure dating. One method of doing so is by calibration from the luminescence profiles of independently known exposure age samples, provided that the calibration and unknown age samples are from the same region and preferably of similar mineralogical composition and orientation (Meyer et al., 2018; Gliganic et al., 2019; Fuhrmann et al., 2022). Previous calibrations have involved the use of historical records (Lehmann et al., 2018), road cut outcrops (Sohbati et al., 2012; Smedley et al., 2021) or the creation of a freshly exposed surface that can be resampled at a later date (Gliganic et al., 2019). In this study, we created sample-specific calibration samples by exposing fresh surfaces for  1 year at each sample site and subsequently exploited the luminescence signal formed within this year of exposure to calculate the unknown σφ0 and μ values for all three lithologies and luminescence signals. At one site, we were able to collect calibration samples in two different orientations to improve our understanding regarding the influence of orientation on a luminescence profile with depth (Supplement).

The data were inverted using a Monte Carlo approach to constrain σφ0, μ and t. Since each unknown age sample has a site-specific calibration sample, the calibration and unknown age samples were solved for simultaneously using the same σφ0 and μ values. To do this, for each sample, at first a luminescence profile with depth for the unknown age surface was generated using random values of σφ0, μ and t and compared with the observed measured profile using a misfit function as follows:

(2) misfit unknown = i = 1 n 1 a L n T n meas ( i ) - L n T n pred ( i ) ,

where n is the number of rock slices in a sample, a is the standard deviation of the plateau in the luminescence–depth profile determined qualitatively and LnTn is the luminescence signal for each rock slice where (LnTn)meas(i) represents the luminescence signal measured in the sample and (LnTn)pred(i) is the luminescence signal predicted using the random parameter values and Eq. (1).

This misfit calculation was then also done for the known-age calibration surface with the same values of σφ0 and μ, although the exact number of days of exposure ( 1 year) was used instead of the randomly generated t value applied for the predicted unknown age profile. Next, the sum of the misfits (misfitcombined) generated from the two profiles was used to estimate a likelihood value using

(3) L = exp - 1 2 misfit combined .

Finally, a rejection algorithm of likelihood <R was applied, where R was a randomly generated value between 0 and 1. A probability density function of σφ0, μ and t was constructed from the values that were retained. To ensure that the parameter space was sufficiently explored, we ran 1.25 × 108 trials during the Monte Carlo search for each individual sample, with σφ0 values between 10−7 and 10−5s−1, μ between 1 and 3 mm−1 and t between 1 and 200 years.

2.3 Estimating erosion rates

Since the 10Be concentrations and OSL depth profiles in a rock surface are both influenced by exposure and surface erosion, an erosion history can be inferred by jointly inverting the 10Be and OSL data, as described in Lehmann et al. (2019). While the effects of complex, stochastic erosion histories have been investigated (Brown and Moon, 2019), here we assume a simple, step-wise erosion history where, at a specific time in the past, the surface goes from experiencing no erosion to an instantaneous onset of fixed rate of erosion. The inversion method tests random pairs of erosion (ε˙) and erosion onset times (ts) in log space to find the pairs representative of erosion histories which are most likely to successfully reproduce the 10Be and OSL data measured from the bedrock surfaces. In this study, we tested 104 pairs of ε˙ and ts with a range of possible ε˙ values from 10−6 to 10−2mm a−1 and ts values from 10−1 to 104a. The σφ0 and μ values used in the inversion were the median values found in Tables S2–S4.

3 Results

3.110Be ages

The 10Be apparent exposure ages, assuming zero correction for erosion, showed no trend with elevation (Table 1). The highest-elevation sample (GG01) is younger than suggested from geomorphic ice thickness reconstructions (Bini et al., 2009), which could reflect periglacial, rather than glacial, exposure (Gallach et al., 2018, 2020). However, post-glacial erosion of the surface cannot be eliminated entirely as there exist large uncertainties in LGM ice thickness reconstructions due to discrepancies in results derived from geomorphological observations and models/simulations. These differences are up to 800 m in some areas of the European Alps (Becker et al., 2016, 2017). On the other hand, samples GG02 and GG03 yield ages commensurate with the decay of the Egesen stadial glaciers, which has been dated to 13.0–11.5 ka (e.g. Ivy-Ochs et al., 2009; Protin et al., 2019), and these results are in agreement with a set of exposure ages calculated from polished bedrock samples at the nearby Triftje Glacier (Kronig et al., 2017). This provides valuable information on Younger Dryas ice thicknesses in the European Alps and interestingly is similar to findings in Lehmann et al. (2020) nearby. The 10Be ages for the three lower-elevation samples (GG04, GG05 and GG06) of 0.35–3.52 ka contrast with information obtained from the Swiss Federal Office of Topography's old geological maps and aerial photos, which show these surfaces were only exposed 22–40 years ago. This implies that the samples suffer from inheritance, which is noteworthy because it would mean that the Gorner Glacier advanced at one point during the Holocene but did not erode the  3 m necessary to reset the 10Be signal, as one may expect. The occurrence of inheritance in the three lower-elevation samples reveals the complicated exposure history these surfaces have experienced, reinforced by studies across the European Alps which imply repeated oscillations in glacier extent during the Holocene following the Egesen stadial (e.g. Hormes et al., 2001; Ivy-Ochs et al., 2009; Goehring et al., 2011; Kronig et al., 2017; Protin et al., 2019). While glaciers in the European Alps were likely smaller than the present day during the middle Holocene (e.g. Solomina et al., 2015), evidence for subsequent re-advances in glacier extent at the Gorner Glacier during the Löbben period and Little Ice Age has been determined (Holzhauser, 1995, 2010) as well as at various sites across the European Alps (e.g. Holzhauser, 1995, 2010; Ivy-Ochs et al., 2009; Schimmelpfennig et al., 2014; Kronig et al., 2017; Protin et al., 2019).

Figure 2Probability distribution inversion results for the unknown IRSL50 parameters σφ0, μ and t (a–c) and luminescence–depth profiles (d) from sample GG06. In the luminescence–depth profile, the blue dots represent IRSL50 luminescence measurements for the unknown exposure age sample, and the green dots are for the known age calibration sample that was exposed for  1 year. The dashed lines are the corresponding model fits, using the median σφ0 and μ values and their respective exposure times. Measurement errors are derived from the square root of the luminescence counts.


3.2 OSL unknown parameters

For all the samples, at least three cores were measured for both the unknown age and calibration samples, and visual assessment of the luminescence profiles with depth confirmed that the surfaces had recorded only one exposure event. Furthermore, as expected, all unknown age sample bleaching consistently penetrated to a deeper depth when compared with their corresponding calibration sample. The results from investigating the effect of calibration sample orientation revealed that, across the three luminescence signals, the σφ0 and μ values from the two different orientations overlap within 2σ (Supplement). This suggests that the effects of sampling a calibration sample in a different orientation to the unknown age sample are minimal.

Table 4Summary of inferred erosion histories across the three luminescence signals.

Download Print Version | Download XLSX

For all three luminescence signals, the inversion produced probability density functions of the unknown luminescence parameters, with the best-suited σφ0, μ and t values for each sample summarised in Tables S2–S4. Figure 2 shows an example result of the IRSL50 unknown parameter inversions taken from sample GG06. Overall, the μ values ranged from 0.59 to 2.45, 1.24 to 2.55 and 1.03 to 2.66 mm−1 for the IRSL50, OSL125 and post-IR IRSL225 signals, respectively. For the σφ0 values, the values ranged from 9.17 × 10−7 to 1.82 × 10−6s−1 for IRSL50, 1.33 × 10−7 to 1.50 × 10−6s−1 for OSL125 and 1.07 × 10−7 to 7.34 × 10−6s−1 for post-IR IRSL225. As shown in Fig. S3, the σφ0 solutions are all of comparable magnitude and overlap within 1σ. This is promising as the σφ0 parameter is region- and mineral-dependent (Sohbati et al., 2011), and it is expected that samples from the same location will share similar values. In contrast, the μ parameter results vary more than anticipated. Although, for each sample, aside from sample GG02, the μ values from the three signals all overlap within 1σ, samples of the same lithology down the transect do not have overlapping results, indicating no trend between μ and lithology for this sample set. We speculate that the observed spread reflects mineralogical variations (e.g. Meyer et al., 2018). Finally, the apparent OSL exposure ages (t) calculated were orders of magnitude lower than suggested by their setting and 10Be results (Tables S2–S4).

Figure 3Luminescence–depth profiles and probability distribution inversion results for sample GG01 across the three luminescence signals used in this study – IRSL50 (a, b), OSL125 (c, d) and post-IR IRSL225 (e, f) – which were inverted independently. For the luminescence profiles (a, c, e), the blue dots represent the luminescence measurements at that particular depth. Measurement errors are derived from the square root of the luminescence counts. The dashed black line represents the reference profile expected when using the 10Be exposure age, uncorrected for erosion, and the red lines are the inverted solutions using the erosion model and the values of ε˙ and ts deemed most likely to fit the data (likelihood > 0.95). The white zones (panels b, d and f) represent the pairs of ε˙ and ts which cannot predict the 10Be concentration of the sample.


Figure 4IRSL50 luminescence–depth profiles and probability distribution inversion results for samples GG02 (a, b) and GG04 (c, d). For the luminescence profiles (a, c), the blue dots represent IRSL50 luminescence measurements at that particular depth. Measurement errors are derived from the square root of the luminescence counts. The dashed black line represents the reference profile expected when using the 10Be exposure age, uncorrected for erosion, and the red lines are the inverted solutions using the erosion model and the values of ε˙ and ts deemed most likely to fit the data (likelihood > 0.95). The luminescence–depth profile for GG04 (c) also includes a reference profile (dashed green line) when using an exposure age informed from archives of old geological maps and aerial photos. This profile is overlain by the other profiles generated using the erosion model, confirming the information gleaned from the corresponding probability distribution plot that this surface has experienced very low erosion rates. The probability distributions highlight the difference between a sample in steady (b) or transient (d) state with erosion. In panel (b), the white zone represents the pairs of ε˙ and ts which cannot predict the 10Be concentration of sample GG02.


3.3 Erosion histories

The inversion outcomes for ε˙ and ts for all three luminescence signals, following the method described in Sect. 2.4, are reported in Table 4. The inverted steady-state erosion rates calculated across the three signals are generally consistent – an example from sample GG01 is shown in Fig. 3. Across all the samples, for the IRSL50 signal, the rates varied from 3.43 × 10−2 to 0.24 mm a−1, OSL125 from 3.13 × 10−2 to 0.18 mm a−1 and post-IR IRSL225 from 1.12 × 10−2 to 0.16 mm a−1 (Table 4). Since the three lower-elevation samples suffered from inheritance rendering the 10Be ages unusable, exposure age information from the historical maps and photos was employed (Table 1) for the inversion of these surfaces' post-glacier erosion rates using a slightly altered version of the inversion code. Of the six samples, the majority had luminescence profiles in steady state with erosion, thus allowing for the extraction of ε˙ and ts values. Conversely, some samples (e.g. GG04) or signals (e.g. OSL125 for GG02) reflected a transient state whereby a wide range of ε˙ and ts combinations was able to explain the data, and so these samples were excluded from further analysis.

Luminescence–depth profiles and probability density plots of ε˙ and ts were generated for each sample, and the IRSL50 results from samples GG02 and GG04 are shown as examples in Fig. 4. Each luminescence–depth plot includes the experimental data measured from the samples as well as a reference profile (dashed black line) plotted using solely the 10Be age without correcting for erosion. As seen in Fig. 4, there is a clear mismatch between the depth of the experimental data when compared with that of the 10Be reference plot – in the case of sample GG02 it translates to a difference in depth of 12 mm. As they both record the same exposure event, this discrepancy in depth confirms the presence of surface erosion and explains the unexpectedly low apparent OSL exposure ages mentioned in Sect. 3.2. Erosion removes material from the surface and therefore alters the depth of the luminescence profiles in samples so that they are shallower than what would be observed with a non-eroding profile. Previous studies (Lehmann et al., 2018, 2019; Sohbati et al., 2018; Smedley, 2021) have also reported underestimations in OSL apparent ages as a result of high erosion rates (> 1 mm ka−1). When plotting a profile using Eq. (1) and the most likely solutions of ε˙ and ts (likelihood > 0.95), it is immediately clear that these are a better fit to the experimental data (red lines).

Combining the inferred ε˙ at steady state with the corresponding minimum ts, we can calculate the minimum amount of material that has been removed from the surfaces as a result of erosion. For all three signals, the results indicate that the higher-elevation samples have had more material removed than the lower-elevation samples. For example, the IRSL50 data suggest that the highest-elevation samples, GG01 and GG02, have had at least 7 and 20 mm removed, respectively, as opposed to the lower-elevation samples, GG05 and GG06, which have lost 4 and 1 mm. This is also supported by the natural texture of the sites.

4 Discussion

4.1 Steady-state/transient erosion histories

As mentioned in Sect. 3.3, while the majority of the sample and signal outputs were in steady state, some transient erosion states were also observed. In the case of sample GG04, where all three signals consistently indicate a system in a transient state of erosion, this likely reflects a response to a localised stochastic erosion process (e.g. surface spallation) that has removed sufficient surface material to place all the luminescence–depth profiles in disequilibrium. In some cases, however, within the same sample there existed discrepancies between the three luminescence signals – some signals were in steady state and others in a transient state, a phenomenon that has also been observed by Smedley et al. (2021). Here, this was seen for the IRSL50 and post-IR IRSL225 signals in sample GG03 and the OSL125 signal in sample GG02. One potential explanation lies in the luminescence bleaching depths (Fig. S4). As expected, for all the samples, the IRSL50 signal is always bleached deepest, which has also been reported in previous studies (e.g. Smedley et al., 2021; Fuhrmann et al., 2022). On the other hand, there is no clear pattern with regards to the OSL125 and post-IR IRSL225 signals. In some samples, the OSL125 signal is bleached more deeply (e.g. GG05) than the post-IR IRSL225, yet in other samples the opposite trend occurs (e.g. GG02). Since, for these samples, both signals are generally more difficult to bleach than the IRSL50 signal, this increases their sensitivity to erosion and thus renders them prone to transient states of erosion, which could explain the presence of both transient and steady-state erosion systems across signals within the same sample. Unfortunately, the reason for the transient IRSL50 signal in sample GG03 remains unclear.

Figure 5Distribution of inverted bedrock surface erosion rates ε˙ with elevation (a) and surface slope (b). The samples presented in this study are shown in blue (dark blue is the average values, light blue is the individual values from all three luminescence signals) alongside results from a nearby study at the Mont Blanc Massif in green (Lehmann et al., 2020).


4.2 Dominant influences on post-glacier erosion rates

Several factors, often working in combination with each other, modulate bedrock surface erosion rates. These include, but are not limited to, temperature, elevation and surface slope. A global compilation of 10Be erosion rate measurements from bedrock surfaces, integrated over 103–106 years across various tectonic settings, climate zones and lithologies, failed to find a parameter that strongly dictates outcrop erosion rates (Portenga and Bierman, 2011). This contrasted with the results from drainage basin erosion rates, where mean basin slope was revealed to be the most dominant factor (Portenga and Bierman, 2011). Lithology is known to sometimes play a dominant role in modulating the bedrock surface erosion (e.g. Twidale, 1982; Ford and Williams, 1989; Moses et al., 2014). This is because the degradation of rocks (weathering) by either chemical or physical means is what subsequently provides material for transport (erosion), and the rate of this breakdown can be primarily controlled by rock lithology. For example, studies in northern Europe (André, 2002b; Nicholson, 2008) calculated post-glacier erosion rates using reference surfaces and suggested that, in some cases, lithology and/or biotic influences have a greater influence on the breakdown of crystalline rocks than environmental or climatic factors. However, in this study, there does not appear to be a relationship between the erosion rates calculated and the lithology of the samples (Fig. S5). This could be explained by the metamorphic nature of the rocks (hornfels, schist and gneiss) that renders the surfaces more resistant than other lithologies. Since lithology is not a dominant influence on weathering in our study area, the erosion rate must be controlled by other environmental or climatic factors.

To further investigate the potential influences on erosion rate in our study area, the inferred post-glacier erosion rates were plotted against elevation and surface slope (Fig. 5). When looking at the signals individually, the OSL125 and post-IR IRSL225 results reveal an anti-correlation between post-glacier erosion rates and elevation, whereas no trend is observed in the IRSL50 data (Table 4). Although the luminescence signals target different minerals and traps, they are all still from the same sample and thus have experienced the same history. Based on this, an average of the three signals was calculated for each site to generate one post-glacier erosion rate value. The trend in erosion rate and elevation observed for the OSL125 and post-IR IRSL225 data is maintained when analysing these erosion averages down the transect (Table 4, Fig. 5). Overall, the results exhibit a strong negative correlation between average erosion rate and elevation (r2= 0.95) but no correlation between erosion rate and surface slope (r2= 0.03). This trend between erosion rate and elevation is in agreement with a study undertaken at the Mont Blanc Massif nearby, which also found a negative correlation between erosion rate and elevation (r2= 0.53) that was stronger than the positive correlation between erosion rate and slope (r2= 0.22) (Lehmann et al., 2020).

The anti-correlation between erosion rate and elevation at these two sites is surprising since surfaces at these elevations in mountain environments are typically exposed to frost crack weathering. This occurs when rocks are subjected to temperatures between 3 and 8 C, termed the frost crack window (Anderson, 1998), and we would therefore expect increasing erosion rates with elevation. A similar observation was made by Small et al. (1997), who found that the bedrock erosion rates from their study, located in an alpine setting, were surprisingly similar to values reported from other environments (excluding arid settings) even though frost crack weathering should be more present. While the presence of snow can help protect the bedrock by maintaining the bedrock surface temperature at around 0 C, the lack of correlation between slope and erosion rate for this study site, and the weak correlation at the Mont Blanc Massif, imply that frost crack weathering is perhaps not a dominant form of post-glacier erosion in these areas. This is further supported when visually observing the sampling sites (Fig. S1), as, aside from sample GG04, there is no clear evidence of rockfall scars or surface spallation at the other sites, suggesting that bedrock erosion is most likely occurring through continuous grain-by-grain erosion.

Potential explanations for the apparent trend in erosion rate with elevation could include (1) wind erosion as a result of katabatic winds coming from the glacier surface (Oerlemans and Grisogono, 2002), resulting in increased exposure of surfaces near the glacier to wind erosion compared with surfaces higher up the valley sides, (2) the accumulation of water at lower elevations downslope due to gravity, facilitating erosion mechanisms that require the presence of water, (3) increased precipitation at the lower-elevation sites, following the Clausius–Clapeyron relationship which estimates a  25 % increase in the water-holding capacity of the atmosphere compared with the highest-elevation site (Supplement), resulting in greater chemical weathering and subsequent erosion, or (4) observed patterns of glacial erosion in valley profiles due to quarrying and/or abrasion both scaling with ice-sliding velocity (Harbor, 1992; Fabel et al., 2004; Goehring et al., 2011; Wirsig, 2016; Herman et al., 2021) and the subsequent damaging effect of the variation in ice load on the underlying bedrock (e.g. Leith et al., 2014). To further explain (4), higher glacial erosion rates were likely present at lower elevations (due to the fastest ice-sliding velocities) and therefore could have inflicted more damage on the underlying bedrock than surfaces at higher elevations experiencing lower ice-sliding velocities. This would have weakened the lower-elevation surfaces to a greater extent, rendering them more susceptible to post-glacier erosion mechanisms following ice retreat.

Although the two study areas both observe a negative correlation, it is clear from Fig. 5 that the decrease in erosion rate with elevation is more pronounced at the Mont Blanc Massif than in this study area. This difference is likely due to local variations influencing the dominant post-glacier erosional mechanisms. While a definitive explanation for this is still unclear, several possibilities exist. One option is differences in lithology – all the Mont Blanc Massif samples are from the same lithology (granite) and of igneous origin, whereas our sites cover three different lithologies that are all metamorphic. Alternatively, this observation could be due to elevation, since the samples in this study were collected at a higher elevation than the Mont Blanc Massif samples. Finally, the contrast in slope might reflect a potential relationship between erosion rate and exposure time as, due to the stochastic nature of weathering, surfaces exposed for shorter periods of time have the potential to derive higher erosion rates than actual long-term averages (e.g. Koppes and Montgomery, 2009; Ganti et al., 2016; Lehmann et al., 2020; Smedley et al., 2021). Here, the Mont Blanc Massif samples (with the more pronounced relationship) have a greater difference in exposure times (between  20 ka and a few years) than the Gorner Glacier samples (Table 1). Sample GG01 from the Gorner samples does not conform to this potential relationship, since it has a younger exposure age, but this might be due to a periglacial erosional influence on its exposure age. However, even though the slopes of the trends differ, it is encouraging that the erosion rates from the two studies are comparable and that they both present a negative correlation between elevation and erosion rate.

4.3 Comparison with other bedrock erosion studies

Attempts to quantify bedrock surface erosion rates have been made worldwide using a variety of techniques integrated across different timescales. TCN methods are generally representative of long-term averages (103–107 years) (Small et al., 1997; Heimsath and McGlynn, 2007; Portenga and Bierman, 2011), whereas other techniques exist that work on shorter timescales (centennial to millennial), such as comparisons with reference surfaces (André, 2002b; Nicholson, 2008), using the effective radii of curvatures of glacial and landslide boulders as a proxy for erosion (Kirkbride and Bell, 2010) or OSL applications (Sohbati et al., 2018; Lehmann et al., 2019, 2020; Smedley et al., 2021).

Figure 6Summary of the results of studies worldwide which have attempted to calculate subglacial erosion rates (blue) and periglacial erosion rates (red). Results involving bedrock of sedimentary lithology were excluded as they are not comparable with our study area. The arrows signify “greater than” or “less than”. The grey dashed lines are the mean subglacial and periglacial erosion rates, accompanied by the respective mean and standard deviation values. a Centennial timescale using 10Be of surficial bedrock surfaces and a* orbital timescale using 10Be depth profiles in a bedrock core. b Measured at marginal locations using 10Be and 14C and at riegel formations using b* 10Be and b** 14C. c Measured on abraded bedrock using 10Be and then c* calculated an approximate basin-wide erosion rate (incorporating both abrasion and quarrying). d Reported results from two separate glaciers. * Guillon et al. (2015) applied a TCN technique but in sediment flux measurements, as opposed to bedrock.


In the western US mountain ranges, in situ TCN 10Be and 26Al were used to estimate maximum mean bedrock erosion rates of 7.6 × 10−3mm a−1 from high alpine summit surfaces that showed no evidence of past glaciations (Small et al., 1997), while studies in north-western Scotland, using boulder radii measurements of glacial deposits, and southern Norway, using reference quartz veins from ice-scoured bedrock surfaces, calculated erosion rates of 3.3 × 10−3mm a−1 (Kirkbride and Bell, 2010) and 5.5 × 10−4mm a−1 (Nicholson, 2008), respectively. Furthermore, a recent study investigating landslide and glacier erratic boulder erosion rates in the Eastern Pamirs of China using OSL depth profiles found minimum steady-state erosion rates of < 3.8 × 10−5 and 1.72 × 10−3mm a−1 (Sohbati et al., 2018). These rates are up to 4 orders of magnitude less than those presented in this study. However, our values are in general agreement with erosion rates reported from studies with climates broadly similar to our study area. This includes results from the Nepal high Himalayas, a study which measured bedrock TCN 10Be and 26Al in valley ridge crests and sidewalls and reported erosion rates of 8 × 10−2–0.2 mm a−1 (Heimsath and McGlynn, 2007). In Europe, André (2002b) used quartz veins, quartzite layers and microcline phenocrysts as reference surfaces from roches moutonnées and glaciofluvially scoured outcrops to calculate a bedrock surface erosion rate of 0.7–5 mm a−1, and a study in the Mont Blanc Massif, applying the same technique as this study to previously glaciated bedrock surfaces, found post-glacial erosion rates of 3.53 × 10−3–4.3 mm a−1 (Lehmann et al., 2019, 2020). Aside from these periglacial studies, Smedley et al. (2021) combined 10Be measurements in sandstone boulders from a rock avalanche with three luminescence signals (IRSL50, post-IR IRSL150 and post-IR IRSL225) to determine interglacial erosion rates in north-western Scotland over the last  4.5 kyr. While the erosion rates derived were in a transient state, they were still comparable with the erosion rates calculated in this study, although it must be noted that only the lower range of values (6–14 mm a−1) from Smedley et al. (2021) can be realistically maintained over the timescale investigated. In addition, our results are in agreement with a compilation of 10Be bedrock measurements which calculated an average global outcrop erosion rate of 1.2 × 10−2mm a−1 (Portenga and Bierman, 2011).

Surprisingly, when comparing the results of this study, and previously calculated periglacial erosion rates, with estimations of subglacial erosion rates, the results demonstrate comparable orders of magnitude. The notion that subglacial and periglacial erosion rates are more similar than previously thought has been suggested previously (O'Farrell et al., 2009; Guillon et al., 2015). A summary of worldwide subglacial and periglacial erosion rates is displayed in Fig. 6, although readers should remain cautious that the studies in this compilation are integrated over different timescales depending on the method itself, which could introduce bias into settings where erosion is stochastic.

Erosion at a glacier bed primarily occurs through abrasion and plucking. Although the latter is thought to be more dominant, it remains difficult to uncouple the two processes with certainty when estimating erosion rates beneath a glacier. Calculations so far have returned rates that differ by several orders of magnitude (Hallet et al., 1996; Koppes and Montgomery, 2009; Koppes et al., 2015; Herman et al., 2015, 2021), mostly resulting from the sliding velocities of glaciers. It must be noted that the various methods are not only influenced by integration over different timescales, but, for previously glaciated surfaces, also on the duration of ice cover in the corresponding study areas. In the European Alps, direct measurements from quartz veins at a glacier snout in the Swiss Alps presented abrasion rates of 0.9–3.75 mm a−1 (Embleton and King, 1975), while measurements on marble plates cemented to the glacier bed of the Glacier d'Argentière, France, gave a rate of 36 mm a−1 (Boulton, 1979).

Bedrock subglacial rates have also been determined for formerly glaciated surfaces by exploiting the difference in half-lives between TCN 10Be and 14C in bedrock. In these studies, sampling deliberately targeted surfaces which displayed no apparent signs of quarrying to attempt to isolate abrasion rates, and the results produced values between 0.02 and > 5 mm a−1 (Goehring et al., 2011; Wirsig, 2016; Wirsig et al., 2017; Steinemann et al., 2021). The large range is due to differences in sample locations – for example, in Steinemann et al. (2021) the lower erosion rates were taken from marginal positions and the higher rates from the glacial trough. From a more global perspective, the application of 36Cl in Washington, USA, found subglacial erosion rates of 0.09–0.35 mm a−1 (Briner and Swanson, 1998), and bedrock TCN measurements in Greenland using 10Be yielded rates of 0.39–1.1 mm a−1 (Young et al., 2016). In the study done by Young et al. (2016), the authors suggest that their results predominantly reflect subglacial abrasion, due to their sampling strategy, and proceed to estimate a likely total basin-wide erosion rate of 1–1.8 mm a−1 by assuming that 30 %–60 % of a glacier's bedrock erosion budget is attributed to quarrying. An advantage of applying TCN-based measurements to calculate subglacial erosion rates is that this allows for the calculation of erosion rates at multiple points, which could subsequently reveal any potential spatial variability in subglacial erosion rates.

Alternatively, in presently glaciated areas, contemporary sediment volume measurements at, or beyond, a glacier terminus can be coupled with ice velocities to provide insight into glacial erosion rates on a greater catchment scale. In the European Alps, studies applying this have produced values between 0.1 and 1 mm a−1 (Hallet et al., 1996), but findings around the globe have occasionally reported higher values (e.g. Koppes and Montgomery, 2009; Koppes et al., 2015; Cook et al., 2020). For example, in New Zealand, glacier sliding velocities were mapped using remote sensing and combined with sediment flux measurements over 5 months at a glacier front to produce a glacial erosion value of 10 mm a−1 (Herman et al., 2015), while suspended sediment load measurements collected at the Leverett Glacier in Greenland over two melt seasons in 2009–2010 produced a subglacial erosion rate of 4.6 ± 2.6 mm a−1 (Cowton et al., 2012). A full compilation of glacier erosion rates, calculations and methods can be found in Herman et al. (2021). Results from sediment yield studies and the ensuing interpretations of subglacial erosion rates should be treated with caution as there are elements of the method which introduce potential bias. Nevertheless, the challenge of accessibility to the ice–bed interface beneath a glacier renders it difficult to estimate subglacial erosion rates by other means, and sediment flux measurements are often the only data available.

5 Conclusion

This study demonstrates the value of combining 10Be and OSL surface dating techniques for quantifying post-glacier bedrock erosion rate histories across timescales of the order of 101 to 104 years. We extended the method introduced by Lehmann et al. (2019, 2020) by measuring three OSL signals (IRSL50, post-IR IRSL225 and OSL125) for the samples in this study. The results show that using multiple OSL signals can not only yield additional constraints for the method, but also provide information in the absence of other data – for example, the IRSL signals in sample GG03 were not in steady state with erosion and therefore could not be used to calculate an erosion rate, but the OSL125 signal could be used.

Averaging the erosion rate results for the three signals at each sample site resulted in post-glacier erosion rates that vary from 9.72 × 10−2 to 1.51 × 10−1mm a−1. The magnitudes of the erosion rates found here at the Gorner Glacier are in agreement with a nearby study at the Mont Blanc Massif (Lehmann et al., 2019, 2020). Plotting the post-glacier erosion rates against elevation and surface slope for the samples in this study indicates a strong anti-correlation of erosion rate with elevation and no correlation between erosion rate and slope. This is in broad agreement with the results from the Mont Blanc Massif; however, the trends there are more pronounced. We suspect this is a result of local differences, such as lithology and/or elevation, influencing the dominant post-glacier erosion mechanisms present or the reflection of a potential relationship between erosion rate and exposure time. Finally, a global compilation of both subglacial and periglacial erosion rates reveals rates that are more comparable than previously assumed, although subglacial erosion rates remain higher. Nevertheless, this could lead to important implications for landscape evolution models and assessing the impact of Quaternary climate on mountain erosion.

Code availability

The code used in this paper is available here: (Lehmann, 2020).

Data availability

The dataset can be accessed using the following DOI: (Elkadi et al., 2022).


The supplement related to this article is available online at:

Author contributions

FH and JE developed the project. JE organised and undertook the fieldwork with assistance from FH and BL. JE did the OSL and 10Be preparation and analysis and performed the modelling. OS, SI and MC aided in the 10Be dating preparation and analysis. JE, GK and FH interpreted the data and contributed to writing the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors would like to thank Sebastián Vivero Andrade, Dilan Rech, Chloé Bouscary, Aurélien Ballu, Luca Malatesta and Günther Prasicek for fieldwork support as well as the three reviewers for the invaluable feedback.

Financial support

This research has been supported by the Bundesbehörden der Schweizerischen Eidgenossenschaft (grant no. 2017.1136).

Review statement

This paper was edited by Arjen Stroeven and reviewed by Derek Fabel, Rachel Smedley, and Pierre Valla.


Aitken, M. J.: An Introduction to Optical Dating: The Dating of Quaternary Sediments by the Use of Photon-Stimulated Luminescence, Oxford University Press, ISBN 978-0-198-54092-2, 1998. 

Anderson, R. S.: Near-Surface Thermal Profiles in Alpine Bedrock: Implications for the Frost Weathering of Rock, Arct. Antarct. Alp. Res., 30, 362–372, 1998. 

André, M.-F.: Do periglacial landscapes evolve under periglacial conditions?, Geomorphology, 52, 149–164, 2002a. 

André, M.-F.: Rates of postglacial rock weathering on glacially scoured outcrops (Abisko-Riksgränsen area, 68 N), Geogr. Ann. A, 84, 139–150, 2002b. 

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, 2008. 

Ballantyne, C. K.: Paraglacial geomorphology, Quaternary Sci. Rev., 21, 1935–2017, 2002. 

Balter-Kennedy, A., Young, N. E., Briner, J. P., Graham, B. L., and Schaefer, J. M.: Centennial- and Orbital-Scale Erosion Beneath the Greenland Ice Sheet Near Jakobshavn Isbræ, J. Geophys. Res.-Earth, 126, e2021JF006429,, 2021. 

Bierman, P. R. and Caffee, M.: Cosmogenic exposure and erosion history of Australian bedrock landforms, Geol. Soc. Am. Bull., 114, 787–803, 2002. 

Bini, A., Buoncristiani, J.-F., Couterrand, S., Ellwanger, D., Felber, M., Florineth, D., Graf, H. R., Keller, O., Kelly, M., Schlüchter, C., and Schoeneich, P.: Die Schweiz während des letzteiszeitlichen Maximums (LGM), Bundesamt für Landestopografie swisstopo, 2009. 

Becker, P., Seguinot, J., Jouvet, G., and Funk, M.: Last Glacial Maximum precipitation pattern in the Alps inferred from glacier modelling, Geogr. Helv., 71, 173–187, 2016. 

Becker, P., Funk, M., Schlüchter, C., and Hutter, K.: A study of the Würm glaciation focused on the Valais region (Alps), Geogr. Helv., 72, 421–442, 2017. 

Bøtter-Jensen, L., Thomsen, K. J., and Jain, M.: Review of optically stimulated luminescence (OSL) instrumental developments for retrospective dosimetry, Radiat. Meas., 45, 253–257, 2010. 

Boulton, G. S.: Processes of glacier erosion on different substrata, J. Glaciol., 23, 15–38, 1979. 

Briner, J. P. and Swanson, T. W.: Using inherited cosmogenic 36Cl to constrain glacial erosion rates of the Cordilleran ice sheet, Geology, 26, 3–6, 1998. 

Brown, N. D. and Moon, S.: Revisiting erosion rate estimates from luminescence profiles in exposed bedrock surfaces using stochastic erosion simulations, Earth Planet. Sc. Lett., 528, 115842, 2019. 

Christl, M., Vockenhuber, C., Kubik, P. W., Wacker, L., Lachner, J., Alfimov, V., and Synal, H. A.: The ETH Zurich AMS facilities: performance parameters and reference materials, Nucl. Instrum. Meth. B, 294, 29–38, 2013. 

Claude, A., Ivy-Ochs, S., Kober, F., Antognini, M., Salcher, B., and Kubik, P. W.: The Chironico landslide (Valle Leventina, southern Swiss Alps): age and evolution, Swiss J. Geosci., 107, 273–291, 2014. 

Cogez, A., Meynadier, L., Allègre, C., Limmois, D., Herman, F., and Gaillardet, J.: Constraints on the role of tectonic and climate on erosion revealed by two time series analysis of marine cores around New Zealand, Earth Planet. Sc. Lett., 410, 174–185, 2015. 

Cook, S. J., Swift, D. A., Kirkbride, M. P., Knight, P. G., and Waller, R. I.: The empirical basis for modelling glacial erosion rates, Nat. Commun., 11, 759,, 2020. 

Cowton, T., Nienow, P., Batholomew, I., Sole, A., and Mair, D.: Rapid erosion beneath the Greenland ice sheet, Geology, 40, 343–346, 2012. 

Dahl, R.: Post-glacial micro-weathering of bedrock surfaces in the Narvik district of Norway, Geogr. Ann. A, 49, 155–166, 1967. 

Delmas, M., Calvet, M., and Gunnell, Y.: Variability of Quaternary glacial erosion rates – A global perspective with special reference to the Eastern Pyrenees, Quaternary Sci. Rev., 28, 484–498, 2009. 

Dornbusch, U., Moses, C., Robinson, D. A., and Williams, R.: Soft copy photogrammetry to measure shore platform erosion on decadal time scales, J. Coast. Conserv., 11, 193–200, 2008. 

Dunai, T. J.: Cosmogenic Nuclides: Principles, Concepts and Applications in the Earth Surface Sciences, Cambridge University Press, Cambridge, ISBN 978-0-521-87380-2, 2010. 

Durcan, J. A., King, G. E., and Duller, G. A. T.: DRAC: Dose Rate and Age Calculator for trapped charge dating, Quat. Geochronol., 28, 54–61, 2015. 

Elkadi, J., King, G. E., Lehmann, B., and Herman, F.: Reducing variability in OSL rock surface dating profiles, Quat. Geochronol., 64, 101169,, 2021. 

Elkadi, J., Lehmann, B., King, G., Steinemann, O., Ivy-Ochs, S., Christl, M., and Herman, F.: Dataset: Quantification of post-glacier bedrock surface erosion in the European Alps using 10Be and optically stimulated luminescence exposure dating, Zenodo [data set],, 2022. 

Embleton, C. and King, C. A. M.: Glacial and Periglacial Geomorphology, Edward Arnold, London, ISBN 978-0-470-23895-0, 1975. 

Fabel, D., Harbor, J., Dahms, D., James, A., Elmore, D., Horn, L., Daley, K., and Steele, C.: Spatial patterns of glacial erosion at a valley scale derived from terrestrial cosmogenic 10Be and 26Al concentrations in rock, Ann. Assoc. Am. Geogr., 94, 241–255, 2004. 

Ford, D. and Williams, P.: Karst Geomorphology and Hydrology, Unwin Hyman, London, ISBN 978-94-011-7780-1, 601 pp., 1989. 

Freiesleben, T., Sohbati, R., Murray, A. S., Jain, M., al Khasawneh, S., Hvidt, S., and Jakobsen, B.: Mathematical model quantifies multiple daylight exposure and burial events for rock surfaces using luminescence dating, Radiat. Meas., 81, 16–22, 2015. 

Fuhrmann, S., Meyer, M. C., Gliganic, L. A., and Obleitner, F.: Testing the effects of aspect and total insolation on luminescence depth profiles for rock surface exposure dating, Radiat. Meas., 153, 106732, 2022. 

Gallach, X., Ravanel, L., Egli, M., Brandova, D., Schaepman, M., Christl, M., Gruber, S., Deline, P., Carcaillet, J., and Pallandre, F.: Timing of rockfalls in the Mont Blanc massif (Western Alps): evidence from surface exposure dating with cosmogenic 10Be, Landslides, 15, 1991–2000, 2018. 

Gallach, X., Carcaillet, J., Ravanel, L., Deline, P., Ogier, C., Rossi, M., Malet, E., and Garcia-Sellés, D.: Climatic and structural controls on Late-glacial and Holocene rockfall occurrence in high-elevated rock walls of the Mont Blanc massif (western Alps), Earth Surf. Processes, 45, 3071–3091,, 2020. 

Ganti, V., von Hagke, C., Scherler, D., Lamb, M. P., Fischer, W. W., and Avouac, J. P.: Timescale bias in erosion rates of glaciated landscapes, Science Advances, 2, e1600204,, 2016. 

Gliganic, L. A., Meyer, M. C., Sohbati, R., Jain, M., and Barrett, S.: OSL surface exposure dating of a lithic quarry in Tibet: laboratory validation and application, Quat. Geochronol., 49, 199–204, 2019. 

GLIMS Consortium: GLIMS Glacier Database, Version 1, [Analysis_IDs 85831–86727, Paul, F. (submitter); Paul, F. (analyst(s))], NASA National Snow and Ice Data Center, Distributed Active Archive Center, Boulder, Colorado, USA, 2005a. 

GLIMS Consortium: GLIMS Glacier Database, Version 1, [Analysis_IDs 165649–167068, Fischer, M. (submitter); Fischer, M. (analyst(s))]. NASA National Snow and Ice Data Center, Distributed Active Archive Center, Boulder, Colorado, USA, 2005b. 

GLIMS Consortium: GLIMS Glacier Database, Version 1, [Analysis_IDs 328449–335984, Paul, F. (submitter); Paul, F. (analyst(s))], NASA National Snow and Ice Data Center, Distributed Active Archive Center, Boulder, Colorado, USA, 2005c. 

GLIMS Consortium: GLIMS Glacier Database, Version 1, [Analysis_IDs 336896–338862, Paul, F. (submitter); Maisch, M. (analyst(s))], NASA National Snow and Ice Data Center, Distributed Active Archive Center, Boulder, Colorado, USA, 2005d. 

GLIMS Consortium: GLIMS Glacier Database, Version 1. [Analysis_IDs 338863–340868, Fischer, M. (submitter); Fischer, M. (analyst(s))], NASA National Snow and Ice Data Center, Distributed Active Archive Center, Boulder, Colorado, USA, 2005e. 

GLIMS Consortium: GLIMS Glacier Database, Version 1, [Analysis_IDs 341659–343718, Maisch, M. (submitter); Benz, Ch.; Wipf, A. (analyst(s))], NASA National Snow and Ice Data Center, Distributed Active Archive Center, Boulder, Colorado, USA, 2005f. 

Goehring, B. M., Schaefer, J. M., Schluechter, C., Lifton, N. A., Finkel, R. C., Timothy Jull, A. J., Akçar, N., and Alley, R. B.: The Rhone Glacier was smaller than today for most of the Holocene, Geology, 39, 679–692, 2011. 

Gosse, J. C. and Phillips, F. M.: Terrestrial in situ cosmogenic nuclides: theory and application, Quaternary Sci. Rev., 20, 1475–1560, 2001. 

Guillon, H., Mugnier, J.-L., Buoncristiani, J.-F., Carcaillet, J., Godon, C., Prud'homme, C., van der Beek, P., and Vassallo, R.: Improved discrimination of subglacial and periglacial erosion using 10Be concentration measurements in subglacial and supraglacial sediment load of the Bossons glacier (Mont Blanc massif, France), Earth Surf. Processes, 40, 1202–1215,, 2015. 

Habermann, J., Schilles, T., Kalchgruber, R., and Wagner, G. A.: Steps towards surface dating using luminescence, Radiat. Meas., 32, 847–851, 2000. 

Hallet, B., Hunter, L., and Bogen, J.: Rates of erosion and sediment evacuation by glaciers: A review of field data and their implications, Global Planet. Change, 12, 213–235, 1996. 

Harbor, J. M.: Numerical modeling of the development of U-shaped valleys by glacial erosion, Geol. Soc. Am. Bull., 104, 1364–1375, 1992. 

Heimsath, A. M. and McGlynn, R.: Quantifying periglacial erosion in the Nepal high Himalaya, Geomorphology, 97, 5–23, 2007. 

Herman, F. and Champagnac, J.-D.: Plio-Pleistocene increase of erosion rates in mountain belts in response to climate change, Terra Nova, 28, 2–10, 2016. 

Herman, F. and King, G. E.: Luminescence Thermochronometry: Investigating the Link between Mountain Erosion, Tectonics and Climate, Elements, 14, 33–38, 2018. 

Herman, F., Seward, D., Valla, P. G., Carter, A., Kohn, B., Willet, S. D., and Ehlers, T. A.: Worldwide acceleration of mountain erosion under a cooling climate, Nature, 504, 423–426, 2013. 

Herman, F., Beyssac, O., Brughelli, M., Lane, S. T., Leprince, S., Adatte, T., Lin, J. Y. Y., Avouac, J.-P., and Cox, S. C.: Erosion by an Alpine glacier, Science, 350, 193–195, 2015. 

Herman, F., Braun, J., Deal, E., and Prasicek, G.: The response time of glacial erosion, J. Geophys. Res.-Earth, 123, 801–817, 2018. 

Herman, F., De Doncker, F., Delaney, I., Prasicek, G., and Koppes, M.: The impact of glaciers on mountain erosion, Nat. Rev. Earth Environ., 2, 422–435, 2021. 

Holzhauser, H.: Gletscherschwankungen innerhalb der letzten 3200 Jahre am Beispiel des grossen Aletsch- und des Gornergletschers, Neue Ergebniss, Gletscher im ständigen Wandel, Schweizerische Gletscherkommission, ISBN 978-3-7281-2208-7, 101–123, 1995. 

Holzhauser, H.: Zur Geschichte des Gornergletschers ein Puzzle aus historischen Dokumenten und fossilen Hölzern aus dem Gletschervorfeld, Geographisches Institut der Universität Bern, Bern, 2010. 

Hormes, A., Müller, B. U., and Schlüchter, C.: The Alps with little ice: evidence for eight Holocene phases of reduced glacier extent in the Swiss Alps, Holocene, 11, 255–265, 2001. 

Humphrey, N. F. and Raymond, C.: Hydrology, erosion and sediment production in a surging glacier: Variegated Glacier, Alaska, 1982–1983, J. Glaciol., 40, 539–552, 1994. 

Huntley, D. J.: An explanation of the power-law decay of luminescence, J. Phys.-Condens. Mat., 18, 1359–1365, 2006. 

Huntley, D. J., McMullan, W. G., Godfrey-Smith, D. I., and Thewalt, M. LW.: Time-dependent recombination spectra arising from optical ejection of trapped charges in feldspars, J. Lumin., 44, 41–46, 1989. 

Inkpen, R. J. and Jackson, J.: Contrasting weathering rates in coastal, urban and rural areas in southern Britain: preliminary investigations using gravestones, Earth Surf. Processes, 25, 229–238, 2000. 

Inkpen, R. J., Collier, P., and Fontana, D. J. L.: Close-range photogrammetric analysis of rock surfaces, Z. Geomorphol., 120, 67–81, 2000. 

Ivy-Ochs, S.: The Dating of Rock Surfaces Using in Situ Produced 10Be, 26Al and 36Cl, with examples from Antarctica and the Swiss Alps, PhD thesis, ETH, Zurich, 1996. 

Ivy-Ochs, S. and Briner, J. P.: Dating Disappearing Ice with Cosmogenic Nuclides, Elements, 10, 351–356, 2014. 

Ivy-Ochs, S., Kerschner, H., Maisch, M., Christl, M., Kubik, P. W., and Schlüchter, C.: Latest Pleistocene and Holocene glacier variations in the European Alps, Quaternary Sci. Rev., 28, 21–22, 2009. 

Jenkins, G. T. H., Duller, G. A. T., Roberts, H. M., Chiverrell, and R. C., Glasser, N. F.: A new approach for luminescence dating glaciofluvial deposits- high precision optical dating of cobbles, Quaternary Sci. Rev., 192, 263–273, 2018. 

Kars, R. H., Wallinga, J., and Cohen, K. M.: A new approach towards anomalous fading correction for feldspar IRSL dating- tests on samples in field saturation, Radiat. Meas., 43, 786–790, 2008. 

Kirkbride, M. P. and Bell, C. M.: Edge roundness of boulders of Torridonian Sandstone (northwest Scotland): Applications for relative dating and implications for warm and cold climate weathering rates, Boreas, 39, 187–198, 2010. 

Kohl, C. P. and Nishiizumi, K.: Chemical isolation of quartz for measurement of in-situ produced cosmogenic nuclides, Geochim. Cosmochim. Ac., 56, 3583–3587, 1992. 

Koppes, M. and Montgomery, D.: The relative efficacy of fluvial and glacial erosion over modern to orogenic timescales, Nat. Geosci., 2, 644–647, 2009. 

Koppes, M., Hallet, B., Rignot, E., Mouginot, J., Smith Wellner, J., and Boldt, K.: Observed latitudinal variations in erosion as a function of glacier dynamics, Nature, 526, 100–103, 2015. 

Korup, O. and Schlunegger, F.: Bedrock landsliding, river incision, and transcience of geomorphic hillslope-channel coupling, J. Geophys. Res., 12, F03027,, 2007. 

Kronig, O., Ivy-Ochs, S., Hajdas, I., Christl, M., Wirsig, C., and Schlüchter, C.: Holocene evolution of the Triftje- and the Oberseegletscher (Swiss Alps) constrained with 10Be exposure and radiocarbon dating, Swiss J. Geosci., 111, 117–131, 2017. 

Laj, C., Kissel, C., and Beer, J.: High resolution global Paleointensity Stack since 75 kyr (GLOPIS-75) calibrated to Absolute values, Timescales Paleomagnetic Field, Geoph. Monog. Series, 145, 255–265, 2004. 

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

Lapp, T., Kook, M. H., Murray, A. S., Thomsen, K. J., Buylaert, J.-P., and Jain, M.: A new luminescence detection and stimulation head for the Risø TL/OSL reader, Radiat. Meas., 81, 178–184, 2015. 

Laskaris, N. and Liritzis, I.: A new mathematical approximation of sunlight attenuation in rocks for surface luminescence dating, J. Lumin., 131, 1874–1884, 2011. 

Lehmann, B.: OSL surface exposure calibration – new codes, GitHub [code], (last access: 2 September 2022), 2020. 

Lehmann, B., Valla, P. G., King, G. E., and Herman, F.: Investigation of OSL surface exposure dating to reconstruct post-LIA glacier fluctuations in the French Alps (Mer de Glace, Mont Blanc massif), Quat. Geochronol., 44, 64–74, 2018. 

Lehmann, B., Herman, F., Valla, P. G., King, G. E., and Biswas, R. H.: Evaluating post-glacial bedrock erosion and surface exposure duration by coupling in situ optically stimulated luminescence and 10Be dating, Earth Surf. Dynam., 7, 633–662,, 2019. 

Lehmann, B., Herman, F., Valla, P. G., King, G. E., Biswas, R. H., Ivy-Ochs, S., Steinemann, O., and Christl, M.: Postglacial erosion of bedrock surfaces and deglaciation timing: new insights from the Mont Blanc massif (Western Alps), Geology, 48, 139–144, 2020. 

Leith, K., Moore, J. R., Amann, F., and Loew, S.: Sub-glacial extensional fracture development and implications for Alpine valley evolution, J. Geophys. Res.-Earth, 119, 62–81, 2014. 

Lifton, N. A., Sato, T., and Dunai, T. J.: Scaling in situ cosmogenic nuclide production rates using analytical approximations to atmospheric cosmic-ray fluxes, Earth Planet. Sc. Lett., 386, 149–160, 2014. 

Liu, J., Cui, F., Murray, A. S., Sohbati, R., Jain, M., Gao, H., Li, W., Li, C., Li, P., Zhou, T., and Chen, J.: Resetting of the luminescence signal in modern riverbed cobbles along the course of the Shiyang River, China, Quat. Geochronol., 49, 184–190, 2019. . 

Lupker, M., France-Lanord, C., Galy, V., Lavé, J., and Kudrass, H.: Increasing chemical weathering in the Himalayan system since the Last Glacial Maximum, Earth Planet. Sc. Lett., 365, 243–252, 2013. 

Luo, M., Chen, J., Liu, J., Qin, J., Owen, L.A., Han, F., Yang, H., Wang, H., Zhang, B., Yin, J., and Li, Y.: A test of rock surface luminescence dating using Glaciofluvial boulders from the Chinese Padmir, Radiat. Meas., 120, 290–297, 2018. 

Martin, L. C. P., Blard, P.-H., Balco, G., Lave, J., Delunel, R., Lifton, N., and Laurent, V.: The CRep program and the ICE-D production rate calibration database: A fully parameterizable and updated online tool to compute cosmic ray exposure ages, Quat. Geochronol., 38, 25–49, 2017. 

Meyer, M. C., Gliganic, L. A., Jain, M., Sohbati, R., and Schmidmair, D.: Lithological controls on light penetration into rock surfaces- implications for OSL and IRSL surface exposure dating, Radiat. Meas., 120, 298–304, 2018. 

Molnar, P.: Late Cenozoic increase in accumulation rates of terrestrial sediment: How Might Climate Change Have Affected Erosion Rates?, Annu. Rev. Earth Pl. Sc., 32, 67–89, 2004. 

Montgomery, D. R.: Valley formation by fluvial and glacial erosion, Geology, 30, 1047–1050, 2002. 

Moses, C., Robinson, D., and Barlow, J.: Methods for measuring rock surface weathering and erosion: a critical review, Earth Sci. Rev., 135, 141–161, 2014. 

Nicholson, D. T.: Rock control on microweathering of bedrock surfaces in a periglacial environment, Geomorphology, 101, 655–665, 2008. 

Nishiizumi, K., Lal, D., Klein, J., Middleton, R., and Arnold, J.: Production of 10Be and 26Al by cosmic rays in terrestrial quartz in-situ and implications for erosion rates, Nature, 319, 134–136, 1986. 

Nishiizumi, K., Imamura, M., Caffee, M. W., Southon, J. R., Finkel, R. C., and McAninch, J.: Absolute calibration of 10Be AMS standards, Nucl. Instrum. Meth. B, 258, 403–413, 2007. 

Oerlemans, J. and Grisogono, B.: Glacier winds and parametrisation of the related surface heat fluxes, Tellus A, 54, 440–452, 2002. 

O'Farrell, C. R., Heimsath, A. M., Lawson, D. E., Jorgensen, L. M., Evenson, E. B., Larson, G., and Denner, J.: Quantifying periglacial erosion: insights on a glacial sediment budget, Matanuska Glacier, Alaska, Earth Surf. Processes, 34, 2008–2022, 2009. 

Ou, X. J., Roberts, H. M., Duller, G. A. T., Gunn, M. D., and Perkins, W. T.: Attenuation of light in different rock types and implications for rock surface luminescence dating, Radiat. Meas., 120, 305–311, 2018. 

Pavon-Carrasco, F. J., Osete, M. L., Torta, J. M., and De Santis, A.: A geomagnetic field model for the Holocene based on archaeomagnetic and lava flow data, Earth Planet. Sc. Lett., 388, 98–109, 2014. 

Polikreti, K., Michael, C., and Maniatis, Y.: Authenticating marble sculptures with thermoluminescence, Ancient TL, 20, 11–18, 2002. 

Polikreti, K., Michael, C. T., and Maniatis, Y.: Thermoluminescence characteristics of marble and dating of freshly excavated marble objects, Radiat. Meas., 37, 87–94, 2003. 

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

Protin, M., Schimmelpfennig, I., Mugnier, J.-L., Ravanel, L., Le Roy, M., Deline, P., Favier, V., and Buoncristiani, J.-F.: Climatic reconstructions for the Younger Dryas/Early Holocene transition and the Litte Ice Age based on paleo-extents of Argentière glacier (French Alps), Quaternary Sci. Rev., 221, 105863,, 2019. 

Raup, B. H., Racoviteanu, A., Khalsa, S. J. S., Helm, C., Armstrong, R., and Arnaud, Y.: The GLIMS Geospatial Glacier Database: a new tool for studying glacier change, Global Planet. Change, 56, 101–110, 2007. 

Reiners, P. W. and Brandon, M. T.: Using Thermochronology to Understand Orogenic Erosion, Annu. Rev. Earth Pl. Sc., 34, 419–466, 2006. 

Riihimaki, C. A., MacGregor, K. R., Anderson, R. S., Anderson, S. P., and Loso, M. G.: Sediment evacuation and glacial erosion rates at a small alpine glacier, J. Geophys. Res.-Earth, 110, F03003,, 2005. 

Rolland, Y., Petit, C., Saillard, M., Braucher, R., Bourlés, D., Darnault, R., Cassol, D., and ASTER Team: Innner gorges incision history: A proxy for deglaciation? Insights from Cosmic Ray Exposure dating (10Be and 36Cl) or river-polished surfaces (Tinée River, SW Alps, France), Earth Planet. Sc. Lett., 457, 271–278, 2017. 

Ruszkiczay-Rüdiger, Z., Kern, Z., Urdea, P., Madarász, B., Braucher, R., and ASTER Team: Limited glacial erosion during the last glaciation in mid-latitude cirques (Retezat Mts, Southern Carpathians, Romania), Geomorphology, 384, 107719,, 2021. 

Schimmelpfennig, I., Schaefer, J. M., Akcar, N., Koffman, T., Ivy-Ochs, S., Schwartz, R., Finkel, R. C., Zimmerman, S., and Schlüchter, C.: A chronology of Holocene and Little Ice Age glacier culminations of the Steingletscher, Central Alps, Switzerland based on high-sensitivity beryllium-10 moraine dating, Earth Planet. Sc. Lett., 393, 220–230, 2014. 

Sellwood, E. L., Guralnik, B., Kook, M., Prasad, A. K., Sohbati, R., Hippe, K., Wallinga, J., and Jain, M.: Optical bleaching front in bedrock revealed by spatially-resolved infrared photoluminescence, Sci. Rep., 9, 2611,, 2019. 

Small, E. E., Anderson, R. S., Repka, J. L., and Finkel, R.: Erosion rates of alpine bedrock summit surfaces deduced from in situ 10Be and 26Al, Earth Planet. Sc. Lett., 150, 413–425, 1997. 

Smedley, R. K., Small, D., Jones, R. S., Brough, S., Bradley, J., and Jenkins, G. T. H.: Erosion rates in a wet, temperate climate derived from rock luminescence techniques, Geochronology, 3, 525–543, 2021. 

Sohbati, R., Murray, A. S., Jain, M., Buylaert, J.-P., and Thomsen, K.: Investigating the resetting of OSL signals in rock surfaces, Geochronometria, 38, 249–258, 2011. 

Sohbati, R., Murray, A. S., Chapot, M. S., Jain, M., and Pederson, J.: Optically stimulated luminescence (OSL) as a chronometer for surface exposure dating, J. Geophys. Res.-Sol. Ea., 117, B09202,, 2012. 

Sohbati, R., Murray, A. S., Jain, M., and Avner, U.: Age of a prehistoric “Rodedian” cult site constrained by sediment and rock surface dating techniques, Quat. Geochronol., 30, 90–99, 2015. 

Sohbati, R., Liu, J., Jain, M., Murray, A., Egholm, D., Paris, R., and Guralnik, B.: Centennial- to millennial-scale hard rock erosion rates deduced from luminescence-depth profiles, Earth Planet. Sc. Lett., 493, 218–230, 2018. 

Solomina, O. N., Bradley, R. S., Hodgson, D. A., Ivy-Ochs, S., Jomelli, V., Mackintosh, A. N., Nesje, A., Owen, L. A., Wanner, H., Wiles, G. C., and Young, N. E.: Holocene glacier fluctuations, Quaternary Sci. Rev., 111, 9–34, 2015. 

Steinemann, O., Ivy-Ochs, S., Hippe, K., Christl, M., Haghipour, N., and Synal, H.-A.: Glacial erosion by the Trift glacier (Switzerland): Deciphering the development of riegels, rock basins and gorges, Geomorphology, 375, 107533, 2021. 

Stroeven, A. P., Fabel, D., Hättestrand, C., and Harbor, J.: A relict landscape in the centre of the Fennoscandian glaciation: cosmogenic radionuclide evidence of tors preserved through multiple glacial cycles, Geomorphology, 44, 145–154, 2002. 

Turowski, J. M. and Cook, K. L.: Field techniques for measuring bedrock erosion and denudation, Earth Surf. Process., 42, 109–127, 2017. 

Twidale, C. R.: Granite Landforms, Elsevier, Amsterdam, ISBN 978-0-444-42116-6, 372 pp., 1982. 

Uppala, S. M., Kallberg, P. W., Simmons, A. J., Andrae, U., Bechtold, V. D. C., Fiorino, M., Gibson, J. K., Haseler, J., Hernandez, A., Kelly, G. A., Li, X., Onogi, K., Saarinen, S., Sokka, N., Allan, R. P., Andersson, E., Arpe, K., Balmaseda, M. A., Beljaars, A. C. M., Van De Berg, L., Bidlot, J., Bormann, N., Caires, S., Chevallier, F., Dethof, A., Dragosavac, M., Fisher, M., Fuentes, M., Hagemann, S., Holm, E., Hoskins, B. J., Isaksen, L., Janssen, P. A. E. M., Jenne, R., Mcnally, A. P., Mahfouf, J.-F., Morcrette, J.-J., Rayner, N. A., Saunders, R. W., Simon, P., Sterl, A., Trenberth, K. E., Untch, A., Vasiljevic, D., Viterbo, P., and Woollen, J.: The ERA-40 re-analysis, Q. J. Roy. Meteor. Soc., 131, 2961–3012, 2005.  

Vafiadou, A., Murray, A., and Liritzis, I.: Optically stimulated luminescence (OSL) dating investigations of rock and underlying soil from three case studies, J. Archaeol. Sci., 34, 1659–1669, 2007. 

Valla, P. G., Van Der Beek, P. V., and Carcaillet, J.: Dating bedrock gorge incision in the French Western Alps (Ecrins-Pelvoux massif) using cosmogenic 10Be, Terra Nova, 22, 18–25, 2010. 

Willenbring, J. K. and Jerolmack, D. J.: The null hypothesis: Globally steady rates of erosion, weathering fluxes and shelf sediment accumulation during Late Cenozoic mountain uplift and glaciation, Terra Nova, 28, 11–18, 2016. 

Willenbring, J. K. and von Blanckenburg, F.: Long-term stability of global erosion rates and weathering during late-Cenozoic cooling, Nature, 465, 211–214, 2010. 

Wirsig, C.: Constraining the timing of deglaciation of the High Alps and rates of subglacial erosion with cosmogenic nuclides, PhD dissertation No. 22978, ETH, Zurich, 2016. 

Wirsig, C., Zasadni, J., Christl, M., Akçar, N., and Ivy-Ochs, S.: Dating the onset of LGM ice surface lowering in the High Alps, Quaternary Sci. Rev., 143, 37–50, 2016. 

Wirsig, C., Ivy-Ochs, S., Reitner, J. M., Christl, M., Vockenhuber, C., Bichler, M., and Reindl, M.: Subglacial abrasion rates at Goldbergkees, Hohe Tauern, Austria, determined from cosmogenic 10Be and 36Cl concentrations, Earth Surf. Process., 42, 1119–1131, 2017. 

Young, N. E., Briner, J. P., Maurer, J., and Schaefer, J. M.: 10Be measurements in bedrock constrain erosion beneath the Greenland Ice Sheet margin, Geophys. Res. Lett., 43, 11708–11719, 2016. 

Zhang, P., Molnar, P., and Downs, W. R.: Increased sedimentation rates and grain sizes 2–4 Myr ago due to the influence of climate change on erosion rates, Nature, 410, 891–897, 2001. 

Short summary
Glacial and non-glacial processes have left a strong imprint on the landscape of the European Alps, but further research is needed to better understand their long-term effects. We apply a new technique combining two methods for bedrock surface dating to calculate post-glacier erosion rates next to a Swiss glacier. Interestingly, the results suggest non-glacial erosion rates are higher than previously thought, but glacial erosion remains the most influential on landscape evolution.