Transient sediment supply in a high-altitude Alpine environment evidenced through a 10Be budget of the Etages catchment (French Western Alps)

Although beryllium-10 (10Be) concentrations in stream sediments provide useful synoptic views of catchment-wide erosion rates, little is known on the relative contributions of different sediment supply mechanisms to the acquisition of their initial signature in the headwaters. Here we address this issue by conducting a 10Be-budget of detrital materials that characterize the morphogenetic domains representative of high-altitude environments of the European Alps. We focus on the Etages catchment, located in the Ecrins-Pelvoux massif (southeast France), and illustrate how in situ 10Be concentrations can be used for tracing the origin of the sand fraction from the bedload in the trunk stream. The landscape of the Etages catchment is characterized by a geomorphic transient state, high topographic gradients, and a large variety of modern geomorphic domains ranging from glacial environments to scarcely vegetated alluvial plains. Beryllium-10 concentrations measured in the Etages catchment vary from similar to 1 x 104 to 4.5 x 105 atoms per gram quartz, while displaying consistent 10Be signatures within each representative morphogenetic unit. We show that the basic requirements for inferring catchment-wide denudation from 10Be concentration measurements are not satisfied in this small, dynamic catchment. However, the distinct 10Be signature observed for the geomorphic domains can be used as a tracer. We suggest that a terrestrial cosmogenic nuclide (TCN) budget approach provides a valuable tool for the tracing of material origin in basins where the let nature do the averaging' principles may be violated.

The deglaciation history of the Simplon region (southern Swiss Alps) constrained by 10Be exposure dating of ice-molded bedrock surfaces

The deglaciation history of the Swiss Alps after the Last Glacial Maximum involved the decay of several ice domes and the subsequent disintegration of valley glaciers at high altitude. Here we use bedrock exposure dating to reconstruct the temporal and spatial pattern of ice retreat at the Simplon Pass (altitude: ∼2000 m) located 40 km southwest of the ‘Rhône ice dome’. Eleven 10Be exposure ages from glacially polished quartz veins and ice-molded bedrock surfaces cluster tightly between 13.5 ± 0.6 ka and 15.4 ± 0.6 ka (internal errors) indicating that the Simplon Pass depression became ice-free at 14.1 ± 0.4 ka (external error of mean age). This age constraint is interpreted to record the melting of the high valley glaciers in the Simplon Pass region during the warm Bølling–Allerød interstadial shortly after the Oldest Dryas stadial. Two bedrock samples collected a few hundred meters above the pass depression yield older 10Be ages of 17.8 ± 0.6 ka and 18.0 ± 0.6 ka. These ages likely reflect the initial downwasting of the Rhône ice dome and the termination of the ice transfluence from the ice dome across the Simplon Pass toward the southern foreland. There, the retreat of the piedmont glacier in Val d’Ossola was roughly synchronous with the decay of the Rhône ice dome in the interior of the mountain belt, as shown by 10Be ages of 17.7 ± 0.9 ka and 16.1 ± 0.6 ka for a whaleback at ∼500 m elevation near Montecrestese in northern Italy. In combination with well-dated paleoclimate records derived from lake sediments, our new age data suggest that during the deglaciation of the European Alps the decay of ice domes was approximately synchronous with the retreat of piedmont glaciers in the foreland and was followed by the melting of high-altitude valley glaciers after the transition from the Oldest Dryas to the Bølling–Allerød, when mean annual temperatures rose rapidly by ∼3 °C.

Minor inheritance inhibits the calibration of the 10Be production rate from the AD 1717 Val Ferret rock avalanche, European Alps

To calibrate the in situ 10Be production rate, we collected surface samples from nine large granitic boulders within the deposits of a rock avalanche that occurred in AD 1717 in the upper Ferret Valley, Mont Blanc Massif, Italy. The 10Be concentrations were extremely low and successfully measured within 10% analytical uncertainty or less. The concentrations vary from 4829 ± 448 to 5917 ± 476 at g−1. Using the historical age exposure time, we calculated the local and sea level-high latitude (i.e. ≥60°) cosmogenic 10Be spallogenic production rates. Depending on the scaling schemes, these vary between 4.60 ± 0.38 and 5.26 ± 0.43 at g−1 a−1. Although they correlate well with global values, our production rates are clearly higher than those from more recent calibration sites. We conclude that our 10Be production rate is a mean and an upper bound for production rates in the Massif region over the past 300 years. This rate is probably influenced by inheritance and will yield inaccurate (e.g. too young) exposure ages when applied to surface-exposure studies in the area. Other independently dated rock-avalanche deposits in the region that are approximately 103 years old could be considered as possible calibration sites.

Chronology of Lateglacial ice flow reorganization and deglaciation in the Gotthard Pass area, Central Swiss Alps, based on cosmogenic 10Be and in situ 14C

We reconstruct the timing of ice flow reconfiguration and deglaciation of the Central Alpine Gotthard Pass, Switzerland, using cosmogenic 10Be and in situ14C surface exposure dating. Combined with mapping of glacial erosional markers, exposure ages of bedrock surfaces reveal progressive glacier downwasting from the maximum LGM ice volume and a gradual reorganization of the paleoflow pattern with a southward migration of the ice divide. Exposure ages of ∼16–14 ka (snow corrected) give evidence for continuous early Lateglacial ice cover and indicate that the first deglaciation was contemporaneous with the decay of the large Gschnitz glacier system. In agreement with published ages from other Alpine passes, these data support the concept of large transection glaciers that persisted in the high Alps after the breakdown of the LGM ice masses in the foreland and possibly decayed as late as the onset of the Bølling warming. A younger group of ages around ∼12–13 ka records the timing of deglaciation following local glacier readvance during the Egesen stadial. Glacial erosional features and the distribution of exposure ages consistently imply that Egesen glaciers were of comparatively small volume and were following a topographically controlled paleoflow pattern. Dating of a boulder close to the pass elevation gives a minimum age of 11.1 ± 0.4 ka for final deglaciation by the end of the Younger Dryas. In situ14C data are overall in good agreement with the 10Be ages and confirm continuous exposure throughout the Holocene. However, in situ14C demonstrates that partial surface shielding, e.g. by snow, has to be incorporated in the exposure age calculations and the model of deglaciation.

New constraints on the gas age-ice age difference along the EPICA ice cores, 0-50 kyr

Gas is trapped in polar ice sheets at ~50–120 m below the surface and is therefore younger than the surrounding ice. Firn densification models are used to evaluate this ice age-gas age difference (Δage) in the past. However, such models need to be validated by data, in particular for periods colder than present day on the East Antarctic plateau. Here we bring new constraints to test a firn densification model applied to the EPICA Dome C (EDC) site for the last 50 kyr, by linking the EDC ice core to the EPICA Dronning Maud Land (EDML) ice core, both in the ice phase (using volcanic horizons) and in the gas phase (using rapid methane variations). We also use the structured 10Be peak, occurring 41 kyr before present (BP) and due to the low geomagnetic field associated with the Laschamp event, to experimentally estimate the Δage during this event. Our results seem to reveal an overestimate of the Δage by the firn densification model during the last glacial period at EDC. Tests with different accumulation rates and temperature scenarios do not entirely resolve this discrepancy. Although the exact reasons for the Δage overestimate at the two EPICA sites remain unknown at this stage, we conclude that current densification model simulations have deficits under glacial climatic conditions. Whatever the cause of the Δage overestimate, our finding suggests that the phase relationship between CO2 and EDC temperature previously inferred for the start of the last deglaciation (lag of CO2 by 800±600 yr) seems to be overestimated.

Direct north-south synchroniszation of abrupt climate change record in ice cores using Beryllium 10

A new, decadally resolved record of the 10Be peak at 41 kyr from the EPICA Dome C ice core (Antarctica) is used to match it with the same peak in the GRIP ice core (Greenland). This permits a direct synchronisation of the climatic variations around this time period, independent of uncertainties related to the ice age-gas age difference in ice cores. Dansgaard-Oeschger event 10 is in the period of best synchronisation and is found to be coeval with an Antarctic temperature maximum. Simulations using a thermal bipolar seesaw model agree reasonably well with the observed relative climate chronology in these two cores. They also reproduce three Antarctic warming events observed between A1 and A2.

EDML1: a chronology for the EPICA deep 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.

Exposure dating of Late Glacial and pre-LGM moraines in the Cordon de Doña Rosa, Northern/Central Chile (~ 31°S)

Despite the important role of the Central Andes (15–30° S) for climate reconstruction, knowledge about the Quaternary glaciation is very limited due to the scarcity of organic material for radiocarbon dating. We applied 10Be surface exposure dating (SED) on 22 boulders from moraines in the Cordon de Doña Rosa, Northern/Central Chile (~31° S). The results show that several glacial advances in the southern Central Andes occurred during the Late Glacial between ~14.7±1.5 and 11.6±1.2 ka. A much more extensive glaciation is dated to ~32±3 ka, predating the temperature minimum of the global LGM (Last Glacial Maximum: ~20 ka). Reviewing these results in the paleoclimatic context, we conclude that the Late Glacial advances were most likely caused by an intensification of the tropical circulation and a corresponding increase in summer precipitation. High-latitude temperatures minima, e.g. the Younger Dryas (YD) and the Antarctic Cold Reversal (ACR) may have triggered individual advances, but current systematic exposure age uncertainties limit precise correlations. The absence of LGM moraines indicates that moisture advection was too limited to allow significant glacial advances at ~20 ka. The tropical circulation was less intensive despite the maximum in austral summer insolation. Winter precipitation was apparently also insufficient, although pollen and marine studies indicate a northward shift of the westerlies at that time. The dominant pre-LGM glacial advances in Northern/Central Chile at ~32 ka required lower temperatures and increased precipitation than today. We conclude that the westerlies were more intense and/or shifted equatorward, possibly due to increased snow and ice cover at higher southern latitudes coinciding with a minimum of insolation.

Paleo erosion rates and climate shifts recorded by Quaternary cut-and-fill sequences in the Pisco valley, central Peru

Fluvial cut-and-fill sequences have frequently been reported from various sites on Earth. Nevertheless, the information about the past erosional regime and hydrological conditions have not yet been adequately deciphered from these archives. The Quaternary terrace sequences in the Pisco valley, located at ca. 13°S, offer a manifestation of an orbitally-driven cyclicity in terrace construction where phases of sediment accumulation have been related to the Minchin (48–36 ka) and Tauca (26–15 ka) lake level highstands on the Altiplano. Here, we present a 10Be-based sediment budget for the cut-and-fill terrace sequences in this valley to quantify the orbitally forced changes in precipitation and erosion. We find that the Minchin period was characterized by an erosional pulse along the Pacific coast where denudation rates reached values as high as 600±80 mm/ka600±80 mm/ka for a relatively short time span lasting a few thousands of years. This contrasts to the younger pluvial periods and the modern situation when 10Be-based sediment budgets register nearly zero erosion at the Pacific coast. We relate these contrasts to different erosional conditions between the modern and the Minchin time. First, the sediment budget infers a precipitation pattern that matches with the modern climate ca. 1000 km farther north, where highly erratic and extreme El Niño-related precipitation results in fast erosion and flooding along the coast. Second, the formation of a thick terrace sequence requires sufficient material on catchment hillslopes to be stripped off by erosion. This was most likely the case immediately before the start of the Minchin period, because this erosional epoch was preceded by a >50 ka-long time span with poorly erosive climate conditions, allowing for sufficient regolith to build up on the hillslopes. Finally, this study suggests a strong control of orbitally and ice sheet forced latitudinal shifts of the ITCZ on the erosional gradients and sediment production on the western escarpment of the Peruvian Andes at 13° during the Minchin period.

Denudation rates of small transient catchments controlled by former glaciation: The Hörnli nunatak in the northeastern Swiss Alpine Foreland

Extensive glaciers repeatedly occupied the northern Alpine Foreland during the Pleistocene and left a strongly glacially overprinted low slope landscape. Only few islands appeared as nunataks standing above the surface of the large piedmont glacier lobes. These nunatak areas kept their original shape, manifested in steep catchments with mean slopes up to 33 . Even though not glaciated, these catchments where significantly affected by base-level changes occurring as a consequence of phases of glacier advances and retreats. Both domains, the glacially eroded and non-eroded, are therefore prone to different mechanisms and time-scales of fluvial and colluvial re-adjustment. In this study we investigate these effects by exploring the spatial distribution and magnitude of denudation in the Hörnli region of the eastern Swiss Alpine Foreland in the present Interglacial. The area represents both domains in a relatively small area with largely uniform tectonic, lithologic and climatic conditions. The differences in Holocene andscape evolution are investigated using topographic analyses and catchment-averaged denudation rates derived from 10Be concentrations in fluvial quartz sand. We find that in formerly non-glaciated, fluvially dominated catchments close hillslope-channel coupling prevails and that these catchments yield high average denudation rates of 350 mm/ka. Glacially overprinted catchments yielded catchment-wide denudation rates an order of magnitude lower. These low denudation rates are hypothesized to be the consequence of both (i) a dominance of slow hillslope processes and (ii) admixture of high concentration, pre-LGM glacial sediment. This suggests that a) a careful field investigation must accompany the denudation rate studies and b) that the concept of area-weighted cosmogenic nuclide denudation rates must be considered in light of the predominant catchment processes.

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.