High-resolution mineral dust and sea ice proxy records from the Talos Dome ice core

In this study we report on new non-sea salt calcium (nssCa2+, mineral dust proxy) and sea salt sodium (ssNa+, sea ice proxy) records along the East Antarctic Talos Dome deep ice core in centennial resolution reaching back 150 thousand years (ka) before present. During glacial conditions nssCa2+ fluxes in Talos Dome are strongly related to temperature as has been observed before in other deep Antarctic ice core records, and has been associated with synchronous changes in the main source region (southern South America) during climate variations in the last glacial. However, during warmer climate conditions Talos Dome mineral dust input is clearly elevated compared to other records mainly due to the contribution of additional local dust sources in the Ross Sea area. Based on a simple transport model, we compare nssCa2+ fluxes of different East Antarctic ice cores. From this multi-site comparison we conclude that changes in transport efficiency or atmospheric lifetime of dust particles do have a minor effect compared to source strength changes on the large-scale concentration changes observed in Antarctic ice cores during climate variations of the past 150 ka. Our transport model applied on ice core data is further validated by climate model data. The availability of multiple East Antarctic nssCa2+ records also allows for a revision of a former estimate on the atmospheric CO2 sensitivity to reduced dust induced iron fertilisation in the Southern Ocean during the transition from the Last Glacial Maximum to the Holocene (T1). While a former estimate based on the EPICA Dome C (EDC) record only suggested 20 ppm, we find that reduced dust induced iron fertilisation in the Southern Ocean may be responsible for up to 40 ppm of the total atmospheric CO2 increase during T1. During the last interglacial, ssNa+ levels of EDC and EPICA Dronning Maud Land (EDML) are only half of the Holocene levels, in line with higher temperatures during that period, indicating much reduced sea ice extent in the Atlantic as well as the Indian Ocean sector of the Southern Ocean. In contrast, Holocene ssNa+ flux in Talos Dome is about the same as during the last interglacial, indicating that there was similar ice cover present in the Ross Sea area during MIS 5.5 as during the Holocene.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

Drilling into debris-rich basal ice at the bottom of the NEEM (Greenland) borehole

After the NEEM (Greenland) deep ice-core drilling was declared terminated with respect to developing stratigraphic climate reconstructions, efforts were turned toward collecting basal ice-sheet debris and, if possible, drilling into the bedrock itself. In 2010, several meters of banded debris-rich ice were obtained under normal ice-drilling operations with the NEEM version of the Hans Tausen (HT) drill, but further penetration was obstructed by a rock in the path of the drill head at 2537.36 m. During short campaigns in 2011 and 2012, attempts were made to penetrate further using various reinforced ice cutters mounted on the HT drill head, tailored to cut through rock. These had some success in penetrating coarse material, but produced severely damaged cutters. Additionally a 51 mm diameter diamond cutting tipped rock drill was adapted to fit the NEEM drill. With this device, several additional meters of core containing subglacial sediments, rocks and rock fragments were collected. With these tools 1.39 m of additional material were obtained during the 2011 field season, and 7.1 m during 2012. Subglacial water refreezing into the newly formed borehole hindered further penetration, and the bedrock interface was not reached before final closure of the NEEM Camp.

Retrieving the paleoclimatic signal from the deeper part of the EPICA Dome C ice core

An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ18Oice, δ18Oatm, total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica.

Deglaciation records of 17O-excess in East Antarctica: reliable reconstruction of oceanic 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., 2008) revealed differences up to 35 ppm in 17O-excess mean level and evolution for the three sites. Our data showed 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 13 ppm, whereas the TD data show no significant variation from the Last Glacial Maximum (LGM) to the Early Holocene (EH). 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 relative humidity (RH) at the oceanic source region (OSR) are 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 RH at the OSR.

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.

Melting of Northern Greenland during the last interglaciation

Using simulated climate data from the comprehensive coupled climate model IPSL CM4, we simulate the Greenland ice sheet (GrIS) during the Eemian interglaciation with the three-dimensional ice sheet model SICOPOLIS. The Eemian is a period 126 000 yr before present (126 ka) with Arctic temperatures comparable to projections for the end of this century. In our simulation, the northeastern part of the GrIS is unstable and retreats significantly, despite moderate melt rates. This result is found to be robust to perturbations within a wide parameter space of key parameters of the ice sheet model, the choice of initial ice temperature, and has been reproduced with climate forcing from a second coupled climate model, the CCSM3. It is shown that the northeast GrIS is the most vulnerable. Even a small increase in melt removes many years of ice accumulation, giving a large mass imbalance and triggering the strong ice-elevation feedback. Unlike the south and west, melting in the northeast is not compensated by high accumulation. The analogy with modern warming suggests that in coming decades, positive feedbacks could increase the rate of mass loss of the northeastern GrIS, exceeding the recent observed thinning rates in the south.

The impact of different glacial boundary conditions on atmospheric dynamics and precipitation in the North Atlantic region

Using a highly resolved atmospheric general circulation model, the impact of different glacial boundary conditions on precipitation and atmospheric dynamics in the North Atlantic region is investigated. Six 30-yr time slice experiments of the Last Glacial Maximum at 21 thousand years before the present (ka BP) and of a less pronounced glacial state – the Middle Weichselian (65 ka BP) – are compared to analyse the sensitivity to changes in the ice sheet distribution, in the radiative forcing and in the prescribed time-varying sea surface temperature and sea ice, which are taken from a lower-resolved, but fully coupled atmosphere-ocean general circulation model. The strongest differences are found for simulations with different heights of the Laurentide ice sheet. A high surface elevation of the Laurentide ice sheet leads to a southward displacement of the jet stream and the storm track in the North Atlantic region. These changes in the atmospheric dynamics generate a band of increased precipitation in the mid-latitudes across the Atlantic to southern Europe in winter, while the precipitation pattern in summer is only marginally affected. The impact of the radiative forcing differences between the two glacial periods and of the prescribed time-varying sea surface temperatures and sea ice are of second order importance compared to the one of the Laurentide ice sheet. They affect the atmospheric dynamics and precipitation in a similar but less pronounced manner compared with the topographic changes.

Simulated winter circulation types in the North Atlantic and European region for preindustrial and glacial conditions

[1] Winter circulation types under preindustrial and glacial conditions are investigated and used to quantify their impact on precipitation. The analysis is based on daily mean sea level pressure fields of a highly resolved atmospheric general circulation model and focuses on the North Atlantic and European region. We find that glacial circulation types are dominated by patterns with an east-west pressure gradient, which clearly differs from the predominantly zonal patterns for the recent past. This is also evident in the frequency of occurrence of circulation types when projecting preindustrial circulation types onto the glacial simulations. The elevation of the Laurentide ice sheet is identified as a major cause for these differences. In areas of strong precipitation signals in glacial times, the changes in the frequencies of occurrence of the circulation types explain up to 60% of the total difference between preindustrial and glacial simulations.

North-south palaeohydrological contrasts in the central Mediterranean during the Holocene: tentative synthesis and working hypotheses

On the basis of a multi-proxy approach and a strategy combining lacustrine and marine records along a north–south transect, data collected in the central Mediterranean within the framework of a collaborative project have led to reconstruction of high-resolution and well-dated palaeohydrological records and to assessment of their spatial and temporal coherency. Contrasting patterns of palaeohydrological changes have been evidenced in the central Mediterranean: south (north) of around 40° N of latitude, the middle part of the Holocene was characterised by lake-level maxima (minima), during an interval dated to ca. 10 300–4500 cal BP to the south and 9000–4500 cal BP to the north. Available data suggest that these contrasting palaeohydrological patterns operated throughout the Holocene, both on millennial and centennial scales. Regarding precipitation seasonality, maximum humidity in the central Mediterranean during the middle part of the Holocene was characterised by humid winters and dry summers north of ca. 40° N, and humid winters and summers south of ca. 40° N. This may explain an apparent conflict between palaeoclimatic records depending on the proxies used for reconstruction as well as the synchronous expansion of tree species taxa with contrasting climatic requirements. In addition, south of ca. 40° N, the first millennium of the Holocene was characterised by very dry climatic conditions not only in the eastern, but also in the central- and the western Mediterranean zones as reflected by low lake levels and delayed reforestation. These results suggest that, in addition to the influence of the Nile discharge reinforced by the African monsoon, the deposition of Sapropel 1 has been favoured (1) by an increase in winter precipitation in the northern Mediterranean borderlands, and (2) by an increase in winter and summer precipitation in the southern Mediterranean area. The climate reversal following the Holocene climate optimum appears to have been punctuated by two major climate changes around 7500 and 4500 cal BP. In the central Mediterranean, the Holocene palaeohydrological changes developed in response to a combination of orbital, ice-sheet and solar forcing factors. The maximum humidity interval in the south-central Mediterranean started ca. 10 300 cal BP, in correlation with the decline (1) of the possible blocking effects of the North Atlantic anticyclone linked to maximum insolation, and/or (2) of the influence of the remnant ice sheets and fresh water forcing in the North Atlantic Ocean. In the north-central Mediterranean, the lake-level minimum interval began only around 9000 cal BP when the Fennoscandian ice sheet disappeared and a prevailing positive NAO-(North Atlantic Oscillation) type circulation developed in the North Atlantic area. The major palaeohydrological oscillation around 4500–4000 cal BP may be a non-linear response to the gradual decrease in insolation, with additional key seasonal and interhemispheric changes. On a centennial scale, the successive climatic events which punctuated the entire Holocene in the central Mediterranean coincided with cooling events associated with deglacial outbursts in the North Atlantic area and decreases in solar activity during the interval 11 700–7000 cal BP, and to a possible combination of NAO-type circulation and solar forcing since ca. 7000 cal BP onwards. Thus, regarding the centennial-scale climatic oscillations, the Mediterranean Basin appears to have been strongly linked to the North Atlantic area and affected by solar activity over the entire Holocene. In addition to model experiments, a better understanding of forcing factors and past atmospheric circulation patterns behind the Holocene palaeohydrological changes in the Mediterranean area will require further investigation to establish additional high-resolution and well-dated records in selected locations around the Mediterranean Basin and in adjacent regions. Special attention should be paid to greater precision in the reconstruction, on millennial and centennial timescales, of changes in the latitudinal location of the limit between the northern and southern palaeohydrological Mediterranean sectors, depending on (1) the intensity and/or characteristics of climatic periods/oscillations (e.g. Holocene thermal maximum versus Neoglacial, as well as, for instance, the 8.2 ka event versus the 4 ka event or the Little Ice Age); and (2) on varying geographical conditions from the western to the eastern Mediterranean areas (longitudinal gradients). Finally, on the basis of projects using strategically located study sites, there is a need to explore possible influences of other general atmospheric circulation patterns than NAO, such as the East Atlantic–West Russian or North Sea–Caspian patterns, in explaining the apparent complexity of palaeoclimatic (palaeohydrological) Holocene records from the Mediterranean area.

Holocene climatic change and the development of the lake-effect snowbelt in Michigan, USA

Lake-effect snow is an important constraint on ecological and socio-economic systems near the North American Great Lakes. Little is known about the Holocene history of lake-effect snowbelts, and it is difficult to decipher how lake-effect snowfall abundance affected ecosystem development. We conducted oxygen-isotope analysis of calcite in lake-sediment cores from northern Lower Michigan to infer Holocene climatic variation and assess snowbelt development. The two lakes experience the same synoptic-scale climatic systems, but only one of them (Huffman Lake) receives a significant amount of lake-effect snow. A 177-cm difference in annual snowfall causes groundwater inflow at Huffman Lake to be 18O-depleted by 2.3‰ relative to O'Brien Lake. To assess when the lake-effect snowbelt became established, we compared calcite-δ18O profiles of the last 11,500 years from these two sites. The chronologies are based on accelerator-mass-spectrometry 14C ages of 11 and 17 terrestrial-plant samples from Huffman and O'Brien lakes, respectively. The values of δ18O are low at both sites from 11,500 to 9500 cal yr BP when the Laurentide Ice Sheet (LIS) exerted a dominant control over the regional climate and provided periodic pulses of meltwater to the Great Lakes basin. Carbonate δ18O increases by 2.6‰ at O'Brien Lake and by 1.4‰ at Huffman Lake between 9500 and 7000 cal yr BP, suggesting a regional decline in the proportion of runoff derived from winter precipitation. The Great Lakes snowbelt probably developed between 9500 and 5500 cal yr BP as inferred from the progressive 18O-depletion at Huffman Lake relative to O'Brien Lake, with the largest increase of lake-effect snow around 7000 cal yr BP. Lake-effect snow became possible at this time because of increasing contact between the Great Lakes and frigid arctic air. These changes resulted from enhanced westerly flow over the Great Lakes as the LIS collapsed, and from rapidly rising Great Lakes levels during the Nipissing Transgression. The δ18O difference between Huffman and O'Brien lakes declines after 5500 cal yr BP, probably because of a northward shift of the polar vortex that brought increasing winter precipitation to the entire region. However, δ18O remains depleted at Huffman Lake relative to O'Brien Lake because of the continued production of lake-effect snow.

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.

Detailed insight into Arctic climatic variability during MIS 11c at Lake El'gygytgyn, NE Russia

Here we present a detailed multi-proxy record of the climate and environmental evolution at Lake El'gygytgyn, Far East Russian Arctic during the period 430–395 ka covering the marine isotope stage (MIS) 12/11 transition and the thermal maximum of super interglacial MIS 11c. The MIS 12/11 transition at Lake El'gygytgyn is characterized by initial warming followed by a cold reversal implying similarities to the last deglaciation. The thermal maximum of MIS 11c is characterized by full and remarkably stable interglacial conditions with mean temperatures of the warmest month (MTWM) ranging between ca. 10–15 °C; annual precipitation (PANN) ranging between ca. 300–600 mm; strong in-lake productivity coinciding with dark coniferous forests in the catchment; annual disintegration of the lake ice cover; and full mixis of the water column. Such conditions persisted, according to our age model, for ca. 27 ± 8 kyr between ca. 425–398 ka. The Lake El'gygytgyn record closely resembles the climate pattern recorded in Lake Baikal (SE Siberia) sediments and Antarctic ice cores, implying interhemispheric climate connectivity during MIS 11c.

Lake Floor Morphology and Sediment Architecture of Lake Torneträsk, Northern Sweden

ABSTRACT. Here we present datasets from a hydroacoustic survey in July 2011 at Lake Torneträsk, northern Sweden. Our hydroacoustic data exhibit lake floor morphologies formed by glacial erosion and accumulation processes, insights into lacustrine sediment accumulation since the beginning of deglaciation, and information on seismic activity along the Pärvie Fault. Features of glacial scouring with a high-energy relief, steep slopes, and relative reliefs of more than 50 m are observed in the large W-basin. The remainder of the lacustrine subsurface appears to host a broad variety of well preserved formations from glacial accumulation related to the last retreat of the Fennoscandian ice sheet. Deposition of glaciolacustrine and lacustrine sediments is focused in areas situated in proximity to major inlets. Sediment accumulation in distal areas of the lake seldom exceeds 2 m or is not observable. We assume that lack of sediment deposition in the lake is a result of different factors, including low rates of erosion in the catchment, a previously high lake level leading to deposition of sediments in higher elevated paleodeltas, tributaries carrying low suspension loads as a result of sedimentation in upstream lakes, and an overall low productivity in the lake. A clear off-shore trace of the Pärvie Fault could not be detected from our hydroacoustic data. However, an absence of sediment disturbance in close proximity to the presumed fault trace implies minimal seismic activity since deposition of the glaciolacustrine and lacustrine sediments.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.