Deglaciation records of 17O-excess in East Antarctica: reliable reconstruction of oceanic normalized relative humidity from coastal sites

We measured δ17O and δ18O in two Antarctic ice cores at EPICA Dome C (EDC) and TALDICE (TD), respectively, and computed 17O-excess with respect to VSMOW. The comparison of our 17O-excess data with the previous record obtained at Vostok (Landais et al., 2008a) revealed differences up to 35 ppm in 17O-excess mean level and evolution for the three sites. Our data show that the large increase depicted at Vostok (20 ppm) during the last deglaciation is a regional and not a general pattern in the temporal distribution of 17O-excess in East Antarctica. The EDC data display an increase of 12 ppm, whereas the TD data show no significant variation from the Last Glacial Maximum (LGM) to the Early Holocene (EH). A Lagrangian moisture source diagnostic revealed very different source regions for Vostok and EDC compared to TD. These findings combined with the results of a sensitivity analysis, using a Rayleigh-type isotopic model, suggest that normalized relative humidity (RHn) at the oceanic source region (OSR) is a determining factor for the spatial differences of 17O-excess in East Antarctica. However, 17O-excess in remote sites of continental Antarctica (e.g. Vostok) may be highly sensitive to local effects. Hence, we consider 17O-excess in coastal East Antarctic ice cores (TD) to be more reliable as a proxy for RHn at the OSR.

Late Pleistocene glaciation in the Central Andes: Temperature versus humidity control A case study from the eastern Bolivian Andes (17°S) and regional synthesis

A glacier–climate model was used to calculate climatic conditions in a test site on the east Andean slope around Cochabamba (17°S, Bolivia) for the time of the maximum Late Pleistocene glaciation. Results suggest a massive temperature reduction of about − 6.4 °C (+ 1.4/− 1.3 °C), combined with annual precipitation rates of about 1100 mm (+ 570 mm/− 280 mm). This implies no major change in annual precipitation compared with today. Summer precipitation was the source for the humidity in the past, as is the case today. This climate scenario argues for a maximum advance of the paleo-glaciers in the eastern cordillera during the global Last Glacial Maximum (LGM, 20 ka BP), which is confirmed by exposure age dates. In a synthesized view over the central Andes, the results point to an increased summer precipitation-driven Late Glacial (15–10 ka BP) maximum advance in the western part of the Altiplano (18°S–23°S), a temperature-driven maximum advance during full glacial times (LGM) in the eastern cordillera, and a pre- and post-LGM (32 ka BP/14 ka BP) maximum advance around 30°S related to increased precipitation and reduced temperature on the western slope of the Andes. The results indicate the importance of understanding the seasonality and details of the mass balance–climate interaction in order to disentangle drivers for the observed regionally asynchronous past glaciations in the central Andes.

Colonization history of the Swiss Rhine basin by the bullhead (Cottus gobio): inference under a Bayesian spatially explicit framework

The present distribution of freshwater fish in the Alpine region has been strongly affected by colonization events occurring after the last glacial maximum (LGM), some 20,000 years ago. We use here a spatially explicit simulation framework to model and better understand their colonization dynamics in the Swiss Rhine basin. This approach is applied to the European bullhead (Cottus gobio), which is an ideal model organism to study fish past demographic processes since it has not been managed by humans. The molecular diversity of eight sampled populations is simulated and compared to observed data at six microsatellite loci under an approximate Bayesian computation framework to estimate the parameters of the colonization process. Our demographic estimates fit well with current knowledge about the biology of this species, but they suggest that the Swiss Rhine basin was colonized very recently, after the Younger Dryas some 6600 years ago. We discuss the implication of this result, as well as the strengths and limits of the spatially explicit approach coupled to the approximate Bayesian computation framework.

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.

Paleoglaciations in Anatolia: A Schematic Review and First Results

Anatolia is situated in the Eastern Mediterranean region between 36 – 42N and 26 – 45E. The geological records of paleoglaciations in the high terrains of Anatolia are key archives to quantify paleoclimate change in the Eastern Mediterranean area. The climate of the Eastern Mediterranean region is influenced by three main atmospheric systems: the main middle to high latitude westerlies, the mid-latitude subtropical high-pressure systems, and the monsoon climate. Glacial geological studies in Turkey have started in the late 19th century. Glacial deposits are found mainly in the eastern, northeastern and southern part of the Anatolian Peninsula. Anatolia is the fundamental element to understand the interactions between paleoenvironment, climatic variations, and development of the human societies. As the Taurus and Black Sea Mountains are sensitively situated for the paleoclimatic reconstructions, a chronostratigraphic framework on the paleoglaciation should be elaborated. The timing of the Last Glacial Maximum (LGM) in Anatolia is still unknown. Our first results from Kavron Valley (Kaçkar Mountains, NE Turkey) are encouraging for the reconstruction of paleoglaciations in Turkey and related paleoclimatological interpretations although it is presently difficult to pinpoint the classical Last Glacial Maximum – Younger Dryas – Little Ice Age moraine sequences in the field.

Water versus ice: The competing roles of modern climate and Pleistocene glacial erosion in the Central Alps of Switzerland

Recent studies have identified relationships between landscape form, erosion and climate in regions of landscape rejuvenation, associated with increased denudation. Most of these landscapes are located in non-glaciated mountain ranges and are characterized by transient geomorphic features. The landscapes of the Swiss Alps are likewise in a transient geomorphic state as seen by multiple knickzones. In this mountain belt, the transient state has been related to erosional effects during the Late Glacial Maximum (LGM). Here, we focus on the catchment scale and categorize hillslopes based on erosional mechanisms, landscape form and landcover. We then explore relationships of these variables to precipitation and extent of LGM glaciers to disentangle modern versus palaeo controls on the modern shape of the Alpine landscape. We find that in grasslands, the downslope flux of material mainly involves unconsolidated material through hillslope creep, testifying a transport-limited erosional regime. Alternatively, strength-limited hillslopes, where erosion is driven by bedrock failure, are covered by forests and/or expose bedrock, and they display oversteepened hillslopes and channels. There, hillslope gradients and relief are more closely correlated with LGM ice occurrence than with precipitation or the erodibility of the underlying bedrock. We relate the spatial occurrence of the transport- and strength-limited process domains to the erosive effects of LGM glaciers. In particular, strength-limited, rock dominated basins are situated above the equilibrium line altitude (ELA) of the LGM, reflecting the ability of glaciers to scour the landscape beyond threshold slope conditions. In contrast, transport-limited, soil-mantled landscapes are common below the ELA. Hillslopes covered by forests occupy the elevations around the ELA and are constrained by the tree line. We conclude that the current erosional forces at work in the Central Alps are still responding to LGM glaciation, and that the modern climate has not yet impacted on the modern landscape.

Identification of erosional mechanisms during past glaciations based on a bedrock surface model of the central European Alps

The bedrock topography beneath the Quaternary cover provides an important archive for the identification of erosional processes during past glaciations. Here, we combined stratigraphic investigations of more than 40,000 boreholes with published data to generate a bedrock topography model for the entire plateau north of the Swiss Alps including the valleys within the mountain belt. We compared the bedrock map with data about the pattern of the erosional resistance of Alpine rocks to identify the controls of the lithologic architecture on the location of overdeepenings. We additionally used the bedrock topography map as a basis to calculate the erosional potential of the Alpine glaciers, which was related to the thickness of the LGM ice. We used these calculations to interpret how glaciers, with support by subglacial meltwater under pressure, might have shaped the bedrock topography of the Alps. We found that the erosional resistance of the bedrock lithology mainly explains where overdeepenings in the Alpine valleys and the plateau occur. In particular, in the Alpine valleys, the locations of overdeepenings largely overlap with areas where the underlying bedrock has a low erosional resistance, or where it was shattered by faults. We also found that the assignment of two end-member scenarios of erosion, related to glacial abrasion/plucking in the Alpine valleys, and dissection by subglacial meltwater in the plateau, may be adequate to explain the pattern of overdeepenings in the Alpine realm. This most likely points to the topographic controls on glacial scouring. In the Alps, the flow of LGM and previous glaciers were constrained by valley flanks, while ice flow was mostly divergent on the plateau where valley borders are absent. We suggest that these differences in landscape conditioning might have contributed to the contrasts in the formation of overdeepenings in the Alpine valleys and the plateau.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Greenland accumulation and its connection to the large-scale atmospheric circulation in ERA-Interim and paleoclimate simulations

Changes in Greenland accumulation and the stability in the relationship between accumulation variability and large-scale circulation are assessed by performing time-slice simulations for the present day, the preindustrial era, the early Holocene, and the Last Glacial Maximum (LGM) with a comprehensive climate model. The stability issue is an important prerequisite for reconstructions of Northern Hemisphere atmospheric circulation variability based on accumulation or precipitation proxy records from Greenland ice cores. The analysis reveals that the relationship between accumulation variability and large-scale circulation undergoes a significant seasonal cycle. As the contributions of the individual seasons to the annual signal change, annual mean accumulation variability is not necessarily related to the same atmospheric circulation patterns during the different climate states. Interestingly, within a season, local Greenland accumulation variability is indeed linked to a consistent circulation pattern, which is observed for all studied climate periods, even for the LGM. Hence, it would be possible to deduce a reliable reconstruction of seasonal atmospheric variability (e.g., for North Atlantic winters) if an accumulation or precipitation proxy were available that resolves single seasons. We further show that the simulated impacts of orbital forcing and changes in the ice sheet topography on Greenland accumulation exhibit strong spatial differences, emphasizing that accumulation records from different ice core sites regarding both interannual and long-term (centennial to millennial) variability cannot be expected to look alike since they include a distinct local signature. The only uniform signal to external forcing is the strong decrease in Greenland accumulation during glacial (LGM) conditions and an increase associated with the recent rise in greenhouse gas concentrations.

Late Pleistocene Interactions of East and West Antarctic Ice-Flow Regimes: Evidence from the McMurdo Ice Shelf

We present new interpretations of deglaciation in McMurdo Sound and the western Ross Sea, with observationally based reconstructions of interactions between East and West Antarctic ice at the last glacial maximum (LGM), 16 000, 12 000, 8000 and 4000 sp. At the LGM? East Antarctic ice from Mulock Glacier split, one branch turned westward south of Ross Island but the other branch rounded Ross Island before flowing southwest into McMurdo Sound. This flow regime, constrained by an ice saddle north of Ross Island, is consistent with the reconstruction of Stuiver and others (1981a). After the LGM, grounding-line retreat was most rapid in areas with greatest water depth, especially along the Victoria Land coast. By 12 000 sp, the ice-now regime in McMurdo Sound changed to through-flowing Mulock Glacier ice, with lesser contributions from Koettlitz, Blue and Ferrar Glaciers, because the former ice saddle north of Ross Island was replaced by a dome. The modern flew regime was established similar to 4000 BP. Ice derived from high elevations on the Polar Plateau but now stranded on the McMurdo Ice Shelf, and the pattern of the Transantarctic Mountains erratics support our reconstructions of Mulock Glacier ice rounding Minna Bluff but with all ice from Skelton Glacier ablating south of the bluff. They are inconsistent with Drewry's (1979) LGM reconstruction that includes Skelton Glacier ice in the McMurdo-Sound through-flow. Drewry's (1979) model closely approximates our results for 12 000-4000 BP. Ice-sheet modeling holds promise for determining whether deglaciation proceeded by grounding-line retreat of an ice sheet that was largely stagnant, because it never approached equilibrium flowline profiles after the Ross Ice Shelf, grounded, or of a dynamic ice sheet with flowline profiles kept low by active ice streams that extended northward from present-day outlet glaciers after the Ross Ice Shelf grounded.

Antarctic Glacial History Since the Last Glacial Maximum: An Overview of the Record on Land

This overview examines available circum-Antarctic glacial history archives on land, related to developments after the Last Glacial Maximum (LGM). It considers the glacial-stratigraphic and morphologic records and also biostratigraphical information from moss banks, lake sediments and penguin rookeries, with some reference to relevant glacial marine records. It is concluded that Holocene environmental development in Antarctica differed from that in the Northern Hemisphere. The initial deglaciation of the shelf areas surrounding Antarctica took place before 10000 C-14 yrs before present(sp), and was controlled by rising global sea level. This was followed by the deglaciation of some presently ice-free inner shelf and land areas between 10000 and 8000 yr sp. Continued deglaciation occurred gradually between 8000 yr sp and 5000 yr sp. Mid-Holocene glacial readvances are recorded from various sites around Antarctica. There are strong indications of a circum-Antarctic climate warmer than today 4700-2000 yr sp. The best dated records from the Antarctic Peninsula and coastal Victoria Land suggest climatic optimums there from 4000-3000 yr sp and 3600-2600 yr sp, respectively. Thereafter Neoglacial readvances are recorded. Relatively limited glacial expansions in Antarctica during the past few hundred years correlate with the Little Ice Age in the Northern Hemisphere.

Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.