Local erosion rates versus active tectonics: cosmic ray exposure modelling in Provence (south-east France)

Over the past decade, in situ-produced cosmogenic nuclides have revolutionised the study of landscape evolution. In particular, numerous studies have demonstrated that, in active tectonic settings, cosmic ray exposure dating of deformed or displaced geomorphic features makes it possible to quantify long-term deformation rates. In western European countries, erosion due to climatically driven processes and human activities is probably the factor that most limits the accuracy of exposure ages and landscape modification rates. In this study, we present the results of a depth-profiling technique applied to alluvial terraces located along the Rhône and the Moyenne Durance rivers. The expected decrease with depth of the measured ¹⁰Be concentrations has been modelled using a χ² inversion method in order to constrain the exposure history of the alluvial sediments. The results suggest that: (1) over the Quaternary, the local surface erosion rates including both regional uplift and climatically driven processes acting on landforms are on the order of 30 m/Myr in southeastern France, and (2) providing a fairly good bracketing of the exposure age, the modelled abandonment age of alluvial terraces affected by the Moyenne Durance Fault allows estimating incision rates, comparing the alluvial terrace elevations with topographic river profiles, and a minimum vertical slip rate value of roughly 0.02 mm/yr for the southern segment of the Moyenne Durance Fault.     

Complex exposure histories of chert clasts in the late Pleistocene shorelines of Lake Lisan,southern Israel

Activities of ²⁶Al and ¹⁰Be in five chert clasts sampled from two beach ridges of late Pleistocene Lake Lisan, precursor of the Dead Sea in southern Israel, indicate low rates of chert bedrock erosion and complex exposure, burial, and by inference, transport histories. The chert clasts were derived from the Senonian Mishash Formation, a chert‐bearing chalk, which is widely exposed in the Nahal Zin drainage basin, the drainage system that supplied most of the material to the beach ridges. Simple exposure ages, assuming only exposure at the beach ridge sampling sites, range from 35 to 354 ky; using the ratio ²⁶Al/¹⁰Be, total clast histories range from 0·46 to 4·3 My, unrelated to the clasts' current position and exposure period on the late Pleistocene beach ridges, 160–177 m below sea level. Optically stimulated luminescence dating of fine sediments from the same and nearby beach ridges yielded ages of 20·0 ± 1·4 ka and 36·1 ± 3·3 ka. These ages are supported by the degree of soil development on the beach ridges and correspond well with previously determined ages of Lake Lisan, which suggest that the lake reached its highest stand around 27 000 cal. years BP . If the clasts were exposed only once and than buried beyond the range of significant cosmogenic nuclide production, then the minimum initial exposure and the total burial times before delivery to the beach ridge are in the ranges 50–1300 ky and 390–3130 ky respectively. Alternatively, the initial cosmogenic dosing could have occurred during steady erosion of the source bedrock. Back calculating such rates of rock erosion suggests values between 0·4 and 12 m My^?1. The relatively long burial periods indicate extended sediment storage as colluvium on slopes and/or as alluvial deposits in river terraces. Some clasts may have been stored for long periods in abandoned Pliocene and early Pleistocene routes of Nahal Zin to the Mediterranean before being transported again back into the Nahal Zin drainage system and washed on to the shores of Lake Lisan during the late Pleistocene. Copyright © 2003 John Wiley & Sons, Ltd.     

Late Pleistocene Glaciations in the Northwestern Sierra Nevada, California

Pleistocene fluvial landforms and riparian ecosystems in central California responded to climate changes in the Sierra Nevada, yet the glacial history of the western Sierra remains largely unknown. Three glacial stages in the northwestern Sierra Nevada are documented by field mapping and cosmogenic radionuclide surface-exposure (CRSE) ages. Two CRSE ages of erratic boulders on an isolated till above Bear Valley provide a limiting minimum age of 76,400±3800 ¹⁰Be yr. Another boulder age provides a limiting minimum age of 48,800±3200 ¹⁰Be yr for a broad-crested moraine ridge within Bear Valley. Three CRSE ages producing an average age of 18,600±1180 yr were drawn from two boulders near a sharp-crested bouldery lateral moraine that represents an extensive Tioga glaciation in Bear Valley. Nine CRSE ages from striated bedrock along a steep valley transect average 14,100±1500 yr and suggest rapid late-glacial ice retreat from lower Fordyce Canyon with no subsequent extensive glaciations. These ages are generally consistent with glacial and pluvial records in east-central California and Nevada.     

Postglacial to Holocene landscape evolution and process rates in steep alpine catchments

Climate change and high magnitude mass wasting events pose adverse societal effects and hazards, especially in alpine regions. Quantification of such geomorphic processes and their rates is therefore critical but is often hampered by the lack of appropriate techniques and the various spatiotemporal scales involved in these studies. Here we exploit both in situ cosmogenic beryllium-10 (¹⁰Be) and carbon-14 (¹⁴C) nuclide concentrations for deducing exposure ages and tracing of sediment through small alpine debris flow catchments in central Switzerland. The sediment cascade and modern processes we track from the source areas, through debris flow torrents to their final export out into sink regions with cosmogenic nuclides over an unprecedented five-year time series with seasonal resolution. Data from a seismic survey and a 90 m core revealed a glacially overdeepened basin, filled with glacial and paraglacial sediments. Surface exposure dating of fan boulders and radiocarbon ages constrain the valley fill from the last deglaciation until the Holocene and show that most of the fan existed in early Holocene times already. Current fan processes are controlled by episodic debris flow activity, snow (firn) and rock avalanches. Field investigations, digital elevation models (DEMs) of difference and geomorphic analysis agree with sediment fingerprinting with cosmogenic nuclides, highlighting that the bulk of material exported today at the outlet of the subcatchments derives from the lower fans. Cosmogenic nuclide concentrations steadily decrease from headwater sources to distal fan channels due to the incorporation of material with lower nuclide concentrations. Further downstream the admixture of sediment from catchments with less frequent debris flow activity can dilute the cosmogenic nuclide signals from debris flow dominated catchments but may also reach thresholds where buffering is limited. Consequently, careful assessment of boundary conditions and driving forces is required when apparent denudation rates derived from cosmogenic nuclide analysis are upscaled to larger regions. © 2018 John Wiley & Sons, Ltd.     

Basinga: A cell-by-cell GIS toolbox for computing basin average scaling factors,cosmogenic production rates and denudation rates

The calculation of denudation rates from the measured cosmogenic nuclide concentrations in river sediments requires assumptions and approximations. Several different approaches and numerical tools are available in the literature. A widely used analytical approach represents the muogenic production with one or two exponentials, assumes the attenuation length of muons to be constant and also neglects temporal variations in the Earth's magnetic field. The denudation rates are then calculated directly and analytically from the measured concentrations. A second numerical and iterative approach was more recently proposed and considers a more rigorous muogenic production law based on pre-calculated variable attenuation length of muons and accounts for temporal changes of the magnetic field. It also assumes a specific distribution of denudation rates throughout the basin and uses an iterative approach to calculate the basin average denudation rates. We tested the two approaches across several natural basins and found that both approaches provide similar denudation results. Hence, assuming exponential muogenic production and constant attenuation length of muons in the rock has little impact on the derived denudation rates. Therefore, unless a priori known distributions of denudation rates are to be tested, there does not appear to be any particular gain from using the second iterative method which is computationally less effective. Based on these findings, we developed and describe here Basinga , a new ArcGIS® and QGIS toolbox which computes the basin average scaling factors, cosmogenic production rates and denudation rates for several tens of drainage basins together. Basinga follows either the Lal/Stone or the Lifton/Sato/Dunai scaling schemes and includes several optional tools for correcting for topographic shielding, ice cover and lithology. We have also developed an original method for correcting the cosmogenic production rates for past variations in the Earth's magnetic field. © 2019 John Wiley & Sons, Ltd.     

Effects of long soil surface residence times on apparent cosmogenic nuclide denudation rates and burial ages in the Cradle of Humankind,South Africa

In situ cosmogenic nuclides are an important tool for quantifying landscape evolution and dating fossil-bearing deposits in the Cradle of Humankind (CoH), South Africa. This technique mainly employs cosmogenic 10-Beryllium (¹⁰Be) in river sediments to estimate denudation rates and the ratio of 26-Aluminium (²⁶Al) to ¹⁰Be (²⁶Al/¹⁰Be), to constrain ages of sediment burial. Here, we use ¹⁰Be and ²⁶Al concentrations in bedrock and soil above the Rising Star Cave (the discovery site of Homo naledi) to constrain the denudation rate and the exposure history of soil on the surface. Apparent ¹⁰Be-derived denudation rates obtained from pebble- to cobble-sized clasts and coarse-sand in soil (on average 3.59 ± 0.27 m/Ma and 3.05 ± 0.25 m/Ma, respectively) are 2-3 times lower than the bedrock denudation rates (on average 9.46 ± 0.68 m/Ma). In addition, soil samples yield an average ²⁶Al/¹⁰Be ratio (5.12 ± 0.27) that is significantly lower than the surface production ratio of 6.75, which suggests complex exposure histories. These results are consistent with prolonged surface residence of up to 1.5 Ma in vertically mixed soils that are up to 3 m thick. We conclude that the ¹⁰Be concentrations accumulated in soils during the long near-surface residence times can potentially cause underestimation of single-nuclide (¹⁰Be) catchment-wide denudation rates in the CoH. Further, burial ages of cave sediment samples that consist of an amalgamation of sand-size quartz grains could be overestimated if a pre-burial ²⁶Al/¹⁰Be ratio calculated from the surface production is assumed. © 2019 John Wiley & Sons, Ltd.     

Climate and relief-induced controls on the temporal variability of denudation rates in a granitic upland

How soil erosion rates evolved over the last about 100 ka and how they relate to environmental and climate variability is largely unknown. This is due to a lack of suitable archives that help to trace this evolution. We determined in situ cosmogenic beryllium-10 (¹⁰Be) along vertical landforms (tors, boulders and scarps) on the Sila Massif to unravel their local exhumation patterns to develop a surface denudation model over millennia. Due to the physical resistance of tors, their rate of exhumation may be used to derive surface and, thus, soil denudation rates over time. We derived soil denudation rates that varied in the range 0–0.40 mm yr^-1. The investigated boulders, however, appear to have experienced repositioning processes about ~20–25 ka bp and were therefore a less reliable archive. The scarps of the Sila upland showed a rapid bedrock exposure within the last 8–15 ka. Overall, the denudation rates increased steadily after 75 ka bp but remained low until about 17 ka bp . The exhumation rates indicate a denudation pulse that occurred about 17–5 ka bp . Since then the rates have continuously decreased. We identify three key factors for these developments – climate, topography and vegetation. Between 75 and 17 ka bp , climate was colder and drier than today. The rapid changes towards warmer and humid conditions at the Pleistocene–Holocene transition apparently increased denudation rates. A denser vegetation cover with time counteracted denudation. Topography also determined the extent of denudation rates in the upland regime. On slopes, denudation rates were generally higher than on planar surfaces. By determining the exhumation rates of tors and scarps, soil erosion rates could be determined over long timescales and be related to topography and particularly to climate. This is key for understanding geomorphic dynamics under current environmental settings and future climate change. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Resistant rock layers amplify cosmogenically-determined erosion rates

Prior numerical modeling work has suggested that incision into sub-horizontal layered stratigraphy with variable erodibility induces non-uniform erosion rates even if base-level fall is steady and sustained. Erosion rates of cliff bands formed in the stronger rocks in a stratigraphic sequence can greatly exceed the rate of base-level fall. Where quartz in downstream sediment is sourced primarily from the stronger, cliff-forming units, erosion rates estimated from concentrations of cosmogenic beryllium-10 (¹⁰Be) in detrital sediment will reflect the locally high erosion rates in retreating cliff bands. We derive theoretical relationships for threshold hillslopes and channels described by the stream-power incision model as a quantitative guide to the potential magnitude of this amplification of ¹⁰Be-derived erosion rates above the rate of base-level fall. Our analyses predict that the degree of erosion rate amplification is a function of bedding dip and either the ratio of rock erodibility in alternating strong and weak layers in the channel network, or the ratio of cliff to intervening-slope gradient on threshold hillslopes. We test our predictions in the cliff-and-bench landscape of the Grand Staircase in southern Utah, USA. We show that detrital cosmogenic erosion rates in this landscape are significantly higher (median 300 m/Ma) than the base-level fall rate (~75 m/Ma) determined from the incision rate of a trunk stream into a ~0.6 Ma basalt flow emplaced along a 16 km reach of the channel. We infer a 3–6-fold range in rock strength from near-surface P-wave velocity measurements. The approximately four-fold difference between the median ¹⁰Be-derived erosion rate and the long-term rate of base-level fall is consistent with our model and the observation that the stronger, cliff-forming lithologies in this landscape are the primary source of quartz in detrital sediments. © 2020 John Wiley & Sons, Ltd.     

Quantifying erosional equilibrium across a slowly eroding,soil mantled landscape

The textbook concept of an equilibrium landscape, which posits that soil production and erosion are balanced and equal channel incision, is rarely quantified for natural systems. In contrast to mountainous, rapidly eroding terrain, low relief and slow-eroding landscapes are poorly studied despite being widespread and densely inhabited. We use three field sites along a climosequence in South Africa to quantify very slow (2-5 m/My) soil production rates that do not vary across hillslopes or with climate. We show these rates to be indistinguishable from spatially invariant catchment-average erosion rates while soil depth and chemical weathering increase strongly with rainfall across our sites. Our analyses imply landscape-scale equilibrium although the dominant means of denudation varies from physical weathering in dry climates to chemical weathering in wet climates. In the two wetter sites, chemical weathering is so significant that clay translocates both vertically in soil columns and horizontally down hillslope catenas, resulting in particle size variation and the accumulation of buried stone lines at the clay-rich depth. We infer hundred-thousand-year residence times of these stone lines and suggest that bioturbation by termites plays a key role in exhuming sediment into the mobile soil layer from significant depths below the clay layer. Our results suggest how tradeoffs in physical and chemical weathering, potentially modulated by biological processes, shape slowly eroding, equilibrium landscapes. © 2019 John Wiley & Sons, Ltd.     

Quantifying glacial erosion on a limestone bed and the relevance for landscape development in the Alps

Glacial erosion is the basic process that has shaped the landscapes of the Alps. Despite intense research over centuries, and the use of various techniques, determination of glacial erosion rates remains challenging. This is not only because the location where the process occurs is almost inaccessible, but also because it is dependent on many different factors, including ice thickness and velocity, glacier thermal regime and lithology. Reported glacial erosion rates range over several orders of magnitude (0.01 to >10 mm a⁻¹). Most studies focus on crystalline bedrock, whereas few researchers have investigated glacial erosion on limestone. Here we analyse glacially polished bedrock surfaces at the recently deglaciated forefield of the Tsanfleuron glacier, Swiss Alps. The nearly horizontally bedded limestone hosts a well-developed karst system. Meltwater from the glacier drains into the subsurface within a few metres of the ice margin. By combining geomorphological mapping, measurement of cosmogenic ³⁶Cl concentrations of glacially eroded bedrock surfaces and a numerical model (MECED), we quantify at each sample location the amount of rock removed during glacier occupation. The glacial erosion rates calculated from these values range from 0 to 0.08 mm a⁻¹. These are orders of magnitude lower than values measured at comparable sites on crystalline bedrock. The high ³⁶Cl concentrations we measured show that the Tsanfleuron glacier was unable to effectively erode the gently dipping, strongly karstified limestone. We suggest that this effect may play a key role in formation and preservation over many glacial cycles of high-elevation, low-relief limestone plateaus in the Alps. © 2020 John Wiley & Sons Ltd.