Deciphering the driving forces of erosion rates on millennial to million-year timescales in glacially impacted landscapes: An example from the Western Alps

In many regions, tectonic uplift is the main driver of erosion over million-year (Myr) timescales, but climate changes can markedly affect the link between tectonics and erosion, causing transient variations in erosion rates. Here we study the driving forces of millennial to Myr-scale erosion rates in the French Western Alps, as estimated from in situ produced cosmogenic 10Be and a newly developed approach integrating detrital and bedrock apatite fission-track thermochronology. Millennial erosion rates from 10Be analyses vary between ~0.27 and ~1.33 m/kyr, similar to rates measured in adjacent areas of the Alps. Significant positive correlations of millennial erosion rates with geomorphic measures, in particular with the LGM ice thickness, reveal a strong transient morphological and erosional perturbation caused by repeated Quaternary glaciations. The perturbation appears independent of Myr-scale uplift and erosion gradients, with the effect that millennial erosion rates exceed Myr-scale erosion rates only in the internal Alps where the latter are low (<0.4 km/Myr). These areas, moreover, exhibit channels that clearly plot above a general linear positive relation between Myr-scale erosion rates and normalized steepness index. Glacial erosion acts irrespective of rock uplift and thus not only leads to an overall increase in erosion rates but also regulates landscape morphology and erosion rates in regions with considerable spatial gradients in Myr-scale tectonic uplift. Our study demonstrates that climate change, e.g., through occurrence of major glaciations, can markedly perturb landscape morphology and related millennial erosion rate patterns, even in regions where Myr-scale erosion rates are dominantly controlled by tectonics.

Glacier response to the change in atmospheric circulation in the eastern Mediterranean during the Last Glacial Maximum

In this study, we document glacial deposits and reconstruct the glacial history in the Karagöl valley system in the eastern Uludağ in northwestern Turkey based on 42 cosmogenic 10Be exposure ages from boulders and bedrock. Our results suggest the Last Glacial Maximum (LGM) advance prior to 20.4 ± 1.2 ka and at least three re-advances until 18.6 ± 1.2 ka during the global LGM within Marine Isotope Stage-2. In addition, two older advances of unknown age are geomorphologically well constrained, but not dated due to the absence of suitable boulders. Glaciers advanced again two times during the Lateglacial. The older is exposure dated to not later than 15.9 ± 1.1 ka and the younger is attributed to the Younger Dryas (YD) based on field evidence. The timing of the glaciations in the Karagöl valley correlates well with documented archives in the Anatolian and Mediterranean mountains and the Alps. These glacier fluctuations may be explained by the change in the atmospheric circulation pattern during the different phases of North Atlantic Oscillation (NAO) winter indices.

Late Weichselian local ice dome configuration and chronology in Northwestern Svalbard: early thinning, late retreat

The chronology and configuration of the Svalbard Barents Sea Ice Sheet (SBSIS) during the Late Weichselian (LW) are based on few and geographically scattered data. Thus, the timing and configuration of the SBSIS has been a subject of extensive debate. We present provenance data of erratic boulders and cosmogenic 10Be ages of bedrock and boulders from Northwest Spitsbergen (NWS), Svalbard to determine the thickness, configuration and chronology during the LW. We sampled bedrock and boulders of mountain summits and summit slopes, along with erratic boulders from coastal locations around NWS. We suggest that a local ice dome over central NWS during LW drained radially in all directions. Provenance data from erratic boulders from northern coastal lowland Reinsdyrflya suggest northeastward ice flow through Liefdefjorden. 10Be ages of high-elevation erratic boulders in central NWS (687–836 m above sea level) ranging from 18.3 ± 1.3 ka to 21.7 ± 1.4 ka, indicate that the centre of a local ice dome was at least 300 m thicker than at present. 10Be ages of all high-elevation erratics (>400 m above sea level, central and coastal locations) indicate the onset of ice dome thinning at 25–20 ka. 10Be ages from erratic boulders on Reinsdyrflya ranging from 11.1 ± 0.8 ka to 21.4 ± 1.7 ka, indicate an ice cover over the entire Reinsdyrflya during LW and a complete deglaciation prior to the Holocene, but apparently later than the thinning in the mountains. Lack of moraine deposits, but the preservation of beach terraces, suggest that the ice covering this peninsula possibly was cold-based and that Reinsdyrflya was part of an inter ice-stream area covered by slow-flowing ice, as opposed to the adjacent fjord, which possibly was filled by a fast-flowing ice stream. Despite the early thinning of the ice sheet (25–20 ka) we find a later timing of deglaciation of the fjords and the distal lowlands. Several bedrock samples (10Be) from vertical transects in the central mountains of NWS pre-date the LW, and suggest either ice free or pervasive cold-based ice conditions. Our reconstruction is aligned with the previously suggested hypothesis that a complex multi-dome ice-sheet-configuration occupied Svalbard and the Barents Sea during LW, with numerous drainage basins feeding fast ice streams, separated by slow flowing, possibly cold-based, inter ice-stream areas.

The Galactic Cosmic Ray Intensity over the Past 106–109 Years as Recorded by Cosmogenic Nuclides in Meteorites and Terrestrial Samples

Concentrations of stable and radioactive nuclides produced by cosmic ray particles in meteorites allow us to track the long term average of the primary flux of galactic cosmic rays (GCR). During the past ∼10 Ma, the average GCR flux remained constant over timescales of hundreds of thousands to millions of years, and, if corrected for known variations in solar modulation, also during the past several years to hundreds of years. Because the cosmic ray concentrations in meteorites represent integral signals, it is difficult to assess the limits of uncertainty of this statement, but they are larger than the often quoted analytical and model uncertainties of some 30%. Time series of concentrations of the radionuclide 10Be in terrestrial samples strengthen the conclusions drawn from meteorite studies, indicating that the GCR intensity on a ∼0.5 million year scale has remained constant within some ±10% during the past ∼10 million years. The very long-lived radioactive nuclide 40K allows to assess the GCR flux over about the past one billion years. The flux over the past few million years has been the same as the longer-term average in the past 0.5–1 billion years within a factor of ∼1.5. However, newer data do not confirm a long-held belief that the flux in the past few million years has been higher by some 30–50% than the very long term average. Neither does our analysis confirm a hypothesis that the iron meteorite data indicate a ∼150 million year periodicity in the cosmic ray flux, possibly related to variations in the long-term terrestrial climate.

Calibration of cosmogenic noble gas production in ordinary chondrites based on 36 Cl- 36 Ar ages. Part 1: Refined produced rates for cosmogenic 21 Ne and 38 Ar

We measured the concentrations and isotopic compositions of He, Ne, and Ar in bulk samples and metal separates of 14 ordinary chondrite falls with long exposure ages and high metamorphic grades. In addition, we measured concentrations of the cosmogenic radionuclides 10Be, 26Al, and 36Cl in metal separates and in the nonmagnetic fractions of the selected meteorites. Using cosmogenic 36Cl and 36Ar measured in the metal separates, we determined 36Cl-36Ar cosmic-ray exposure (CRE) ages, which are shielding-independent and therefore particularly reliable. Using the cosmogenic noble gases and radionuclides, we are able to decipher the CRE history for the studied objects. Based on the correlation 3He/21Ne versus 22Ne/21Ne, we demonstrate that, among the meteorites studied, only one suffered significant diffusive losses (about 35%). The data confirm that the linear correlation 3He/21Ne versus 22Ne/21Ne breaks down at high shielding. Using 36Cl-36Ar exposure ages and measured noble gas concentrations, we determine 21Ne and 38Ar production rates as a function of 22Ne/21Ne. The new data agree with recent model calculations for the relationship between 21Ne and 38Ar production rates and the 22Ne/21Ne ratio, which does not always provide unique shielding information. Based on the model calculations, we determine a new correlation line for 21Ne and 38Ar production rates as a function of the shielding indicator 22Ne/21Ne for H, L, and LL chondrites with preatmospheric radii less than about 65 cm. We also calculated the 10Be/21Ne and 26Al/21Ne production rate ratios for the investigated samples, which show good agreement with recent model calculations.

Chronostratigraphy of the 600,000 year old continental record of Lake Van

Lake Van sediment cores from the Ahlat Ridge and Northern Basin drill sites of the ICDP project PALEOVAN contain a wealth of information about past environmental processes. The sedimentary sequence was dated using climatostratigraphic alignment, varve chronology, tephrostratigraphy, argon-argon single-crystal dating, radiocarbon dating, magnetostratigraphy, and cosmogenic nuclides. Based on the lithostratigraphic framework, the different age constraints are compiled and a robust and precise chronology of the 600,000 year-old Lake Van record is constructed. Proxy records of total organic carbon content and sediment color, together with the calcium/potassium-ratios and arboreal pollen percentages of the 174-meter-long Ahlat Ridge record, mimic the Greenland isotope stratotype (NGRIP). Therefore, the proxy records are systematically aligned to the onsets of interstadials reflected in the NGRIP or synthesized Greenland ice-core stratigraphy. The chronology is constructed using 27 age control points derived from visual synchronization with the GICC05 timescale, an absolutely-dated speleothem record (e.g., Hulu, Sanbao, Linzhu cave) and the Epica Dome C timescale. In addition, the uppermost part of the sequence is complemented with four ages from Holocene varve chronology and two calibrated radiocarbon ages. Furthermore, nine argon-argon ages and a comparison of the relative paleointensity record of the magnetic field with reference curve PISO-1500 confirm the accuracy of the age model. Also the identification of the Laschamp event via measurements of 10Be in the sediment confirms the presented age model. The chronology of the Ahlat Ridge record is transferred to the 79-meter-long event-corrected composite record from the Northern Basin and supplemented by additional radiocarbon dating on organic marco-remains. The basal age of the Northern Basin record is estimated at ~90 ka. The variations of the time series of total organic carbon content, the Ca/K ratio, and the arboreal pollen percentages illustrate that the presented chronology and paleoclimate data are suited for reconstructions and modeling of the Quaternary and Pleistocene climate evolution in the Near East at millennial timescales. Furthermore, the chronology of the last 250 kyr can be used to test other dating techniques.

Bedrock bedding, landsliding and erosional budgets in the Central European Alps

We explore the controls of the litho-tectonic architecture on the erosional flux in the 370-km2 Glogn basin (European Alps). In this basin, the bedding and schistosity of the bedrock dip parallel to the topographic slope on the NW valley flank, leading to a non-dip slope situation on the opposite SE valley side. While the dip slope condition has promoted the occurrence of landslides (e.g. the c. 30-km2 deep-seated Lumnezia landslide), the opposite non-dip slope side of the valley hosts >100-m-deeply incised tributary streams. 10Be concentrations of stream sediments yield catchment-averaged denudation rates that vary between 0.27 ± 0.03 and 2.19 ± 0.37 mm a−1, while the spatially averaged denudation rate of the entire basin is 1.99 ± 0.34 mm a−1. Our 10Be-based approach reveals that the Lumnezia landslide front contributes c. 30–65% of the entire sediment budget, although it covers <5% of the Glogn basin. This suggests a primary control of the bedrock bedding on erosion rates and processes.

Efficiency of frost-cracking processes through space and time: An example from the eastern Italian Alps

It is widely accepted that climate has a strong impact and exerts important feedbacks on erosional processes and sediment transport mechanisms. However, the extent at which climate influences erosion is still a matter of debate. In this paper we test whether frost-cracking processes and related temperature variations can influence the sediment production and surface erosion in a small catchment situated in the eastern Italian Alps. To this extent, we first present a geomorphic map of the region that we complement with published 10Be-based denudation rates. We then apply a preexisting heat-flow model in order to analyze the variations of the frost-cracking intensity (FCI) in the study area, which could have controlled the sediment production in the basin. Finally, we compare the model results with the pattern of denudation rates and Quaternary deposits in the geomorphic map. The model results, combined with field observations, mapping, and quantitative geomorphic analyses, reveal that frost-cracking processes have had a primary role in the production of sediment where the intensity of sediment supply has been dictated and limited by the combined effect of temperature variations and conditions of bedrock preservation. These results highlight the importance of a yet poorly understood process for the production of sediment in mountain areas.

Timing of retreat of the Reuss Glacier (Switzerland) at the end of the Last Glacial Maximum

We used cosmogenic 10Be and 36Cl to establish the timing of the onset of deglaciation after the Last Glacial Maximum of the Reuss Glacier, one of the piedmont lobes of the Alpine ice cap that reached the northern Alpine foreland in Switzerland. In this study, we sampled erratic boulders both at the frontal position in the foreland (Lenzburg and Wohlen, canton Aargau) and at the lateral Alpine border position (Seeboden moraine, Rigi, canton Schwyz). The minimum age for the beginning of retreat is 22.2 ± 1.0 ka at the frontal (terminal) position and 20.4 ± 1.0 ka at the lateral position. These ages are directly comparable with exposure ages from the other piedmont lobes in the northern Alpine foreland. Our data from the mountain called Rigi, do not support the hypothesis that boulders located external to the Seeboden moraine were deposited prior to the last glacial cycle. We present a first exposure age from an erratic boulder in a retreat position in the Alpine foreland. The Reuss Glacier was approximately 12 km behind the maximal extent no later than at 18.6 ± 0.9 ka.

The first major incision of the Swiss Deckenschotter landscape

The Swiss Deckenschotter (“cover gravels”) is the oldest Quaternary units in the northern Swiss Alpine Foreland. They are a succession of glaciofluvial gravel layers intercalated with glacial and/or overbank deposits. This lithostratigraphic sequence is called Deckenschotter because it “covers” Molasse or Mesozoic bedrock and forms mesa-type hill-tops. Deckenschotter occurs both within and beyond the extent of the Last Glacial Maximum glaciers. The Swiss Deckenschotter consist of two sub-units: Höhere (Higher) and Tiefere (Lower) Deckenschotter. Although the Höhere Deckenschotter sub-unit (HDS) is topographically higher than the Tiefere Deckenschotter, it is older. The only available age for the Swiss Deckenschotter is 2.5–1.8 Ma based on mammal remains found in HDS at the Irchel site. In this study, we present an exposure age for the topographically lowest HDS, calculated from a cosmogenic 10Be depth-profile. Our results show that the first phase of the Deckenschotter glaciations in the Swiss Alps terminated at least 1,020+80−120 ka ago, which is indicated by a significant fluvial incision. This line of evidence seems to be close to synchronous with the beginning of the Mid-Pleistocene Revolution, when the frequency of the glacial-interglacial cyclicity changed from 41 to 100 ka and the amplitude from low to high, between marine isotope stages 23 and 22.

36Cl production rate from K-spallation in the European Alps (Chironico landslide, Switzerland)

The abundant production of in situ cosmogenic 36Cl from potassium renders 36Cl measurements in K-rich rocks or minerals, such as K-feldspars, potentially useful for precisely dating rock surfaces, either in single-nuclide or in multi-nuclide studies, for example combined with 10Be measurements in quartz. However, significant discrepancies in experimentally calibrated 36Cl production rates from spallation of potassium (36PK-sp), referenced to sea-level/high-latitude (SLHL), limit the accuracy of 36Cl dating from K-rich lithologies. We present a new 36Cl calibration using K-feldspars, in which K-spallation is the most dominant 36Cl production pathway (>92% of total 36Cl), thus minimizing uncertainties from the complex multi-pathway 36Cl production systematics. The samples are derived from boulders of an ∼13.4 ka-old landslide in the Swiss Alps (∼820 m, 46.43°N, 8.85°E). We obtain a local 36PK-sp of 306 ± 16 atoms 36Cl (g K)−1 a−1 and an SLHL 36PK-sp of 145.5 ± 7.7 atoms 36Cl (g K)−1 a−1, when scaled with a standard scaling protocol (‘Lm’). Applying this SLHL 36PK-sp to determine 36Cl exposure ages of K-feldspars from 10Be-dated moraine boulders yields excellent agreement, confirming the validity of the new SLHL 36PK-sp for surface exposure studies, involving 36Cl in K-feldspars, in the Alps.

Quaternary uplift rates of the Central Anatolian Plateau, Turkey: insights from cosmogenic isochron-burial nuclide dating of the Kızılırmak River terraces

The Central Anatolian Plateau (CAP) in Turkey is a relatively small plateau (300 × 400 km) with moderate average elevations of ∼1 km situated between the Pontide and Tauride orogenic mountain belts. Kızılırmak, which is the longest river (1355 km) within the borders of Turkey, flows within the CAP and slowly incises into lacustrine and volcaniclastic units before finally reaching the Black Sea. We dated the Cappadocia section of the Kızılırmak terraces in the CAP by using cosmogenic burial and isochron-burial dating methods with 10Be and 26Al as their absolute dating can provide insight into long-term incision rates, uplift and climatic changes. Terraces at 13, 20, 75 and 100 m above the current river indicate an average incision rate of 0.051 ± 0.01 mm/yr (51 ± 1 m/Ma) since ∼1.9 Ma. Using the base of a basalt fill above the modern course of the Kızılırmak, we also calculated 0.05–0.06 mm/yr mean incision and hence rock uplift rate for the last 2 Ma. Although this rate might be underestimated due to normal faulting along the valley sides, it perfectly matches our results obtained from the Kızılırmak terraces. Although up to 5–10 times slower, the Quaternary uplift of the CAP is closely related to the uplift of the northern and southern plateau margins respectively.

Paleohydraulic reconstruction of a 40 kyr-old terrace sequence implies that water discharge was larger than today

The evolution of landscapes crucially depends on the climate history. This is particularly evident in South America where landscape responses to orbital climate shifts have been well documented. However, while most studies have focused on inferring temperature variations from paleoclimate proxy data, estimates of water budget changes have been complicated because of a lack of adequate physical information. Here, we present a methodology and related results, which allowed us to extract water discharge values from the sedimentary record of the 40 Ka-old fluvial terrace deposits in the Pisco valley, western Peru. In particular, this valley hosts a Quaternary cut-and-fill succession that we used, in combination with beryllium-10 (10Be)-based sediment flux, gauging records, channel geometries and grain size measurements, to quantitatively assess sediment and water discharge values c. 40 Ka ago in relation to present-day conditions. We compare these discharge estimates to the discharge regime of the modern Pisco River and find that the water discharge of the paleo-Pisco River, during the Minchin pluvial period c. 40 Ka ago, was c. 7–8 times greater than the modern Pisco River if considering the mean and the maximum water discharge. In addition, the calculations show that inferred water discharge estimates are mainly dependent on channel gradients and grain size values, and to a lesser extent on channel width measures. Finally, we found that the c. 40 Ka-old Minchin terrace material was poorer sorted than the modern deposits, which might reflect that sediment transport during the past period was characterized by a larger divergence from equal mobility compared to the modern situation. In summary, the differences in grain size distribution and inferred water discharge estimates between the modern and the paleo-Pisco River suggests that the 40 Ka-old Minchin period was characterized by a wetter climate and more powerful flood events.

Nd, Pb, and Be isotopes of ferromanaganese crusts of the Atlantic Ocean

The isotopic composition of Nd in present-day deep waters of the central and northeastern Atlantic Ocean is thought to fingerprint mixing of North Atlantic Deep Water with Antarctic Bottom Water. The central Atlantic Romanche and Vema Fracture Zones are considered the most important pathways for deep water exchange between the western and eastern Atlantic basins today. We present new Nd isotope records of the deepwater evolution in the fracture zones obtained from ferromanganese crusts, which are inconsistent with simple water mass mixing alone prior to 3 Ma and require additional inputs from other sources. The new Pb isotope time series from the fracture zones are inexplicable by simple mixing of North Atlantic Deep Water and Antarctic Bottom Water for the entire past 33 Myr. The distinct and relatively invariable Nd and Pb isotope records of deep waters in the fracture zones appear instead to have been controlled to a large extent by contributions from Saharan dust and the Orinoco/Amazon Rivers. Thus the previously observed similarity of Nd and Pb isotope time series from the western and eastern North Atlantic basins is better explainable by direct supply of Labrador Seawater to the eastern basin via a northern pathway rather than by advection of North Atlantic Deep Water via the Romanche and Vema Fracture Zones.

Chronology for the EDML ice core from Dronning Maud Land, Antarctica, over the last 150 000 years

A chronology called EDML1 has been developed for the EPICA ice core from Dronning Maud Land (EDML). EDML1 is closely interlinked with EDC3, the new chronology for the EPICA ice core from Dome-C (EDC) through a stratigraphic match between EDML and EDC that consists of 322 volcanic match points over the last 128 ka. The EDC3 chronology comprises a glaciological model at EDC, which is constrained and later selectively tuned using primary dating information from EDC as well as from EDML, the latter being transferred using the tight stratigraphic link between the two cores. Finally, EDML1 was built by exporting EDC3 to EDML. For ages younger than 41 ka BP the new synchronized time scale EDML1/EDC3 is based on dated volcanic events and on a match to the Greenlandic ice core chronology GICC05 via 10Be and methane. The internal consistency between EDML1 and EDC3 is estimated to be typically ~6 years and always less than 450 years over the last 128 ka (always less than 130 years over the last 60 ka), which reflects an unprecedented synchrony of time scales. EDML1 ends at 150 ka BP (2417 m depth) because the match between EDML and EDC becomes ambiguous further down. This hints at a complex ice flow history for the deepest 350 m of the EDML ice core.