(Table 1) Gas tie points used to correlate the ice cores NGRIP and EDML

Water isotope records from the EPICA Dronning Maud Land (EDML) and the NorthGRIP ice cores have revealed a one to one coupling between Antarctic Isotope Maxima (AIM) and Greenland Dansgaard-Oeschger (DO) events back to 50 kyr. In order to explore if this north-south coupling is persistent over Marine Isotopic Stage 5 (MIS 5), a common timescale must first be constructed. Here, we present new records of d18O of O2 (d18Oatm) and methane (CH4) measured in the air trapped in ice from the EDML (68-147 kyr) and NorthGRIP (70-123 kyr) ice cores. We demonstrate that, through the period of interest, CH4 records alone are not sufficient to construct a common gas timescale between the two cores. Millennial-scale variations of d18Oatm are evidenced over MIS 5 both on the Antarctic and Greenland ice cores and are coupled to CH4 profiles to synchronise the NorthGRIP and EDML records. They are shown to be a precious tool for ice core synchronisation. With this new dating strategy, we produce the first continuous and accurate sequence of the north-south climatic dynamics on a common ice timescale for the last glacial inception and the first DO events of MIS 5, reducing relative dating uncertainties to an accuracy of a few centuries at the onset of DO events 24 to 20. This EDML-NorthGRIP synchronisation provides new firm evidence that the bipolar seesaw is a pervasive pattern from the beginning of the glacial period. The relationship between Antarctic warming amplitudes and their concurrent Greenland stadial duration highlights the particularity of DO event 21 and its Antarctic counterpart. Our results suggest a smaller Southern Ocean warming rate for this long DO event compared to DO events of MIS 3.

The Antarctic ice core chronology (AICC2012)

Abstract of Bazin et al. (2013): An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120-800 ka. In this framework, new measurements of d18Oatm over Marine Isotope Stage (MIS) 11-12 on EDC and a complete d18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from d18Oatm, dO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one. Abstract of Veres et al. (2013): The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on d15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores.

Water stable oxygen isotope records from six Antarctic ice core sites for the present and last interglacial periods

We compare the present and last interglacial periods as recorded in Antarctic water stable isotope records now available at various temporal resolutions from six East Antarctic ice cores: Vostok, Taylor Dome, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML), Dome Fuji and the recent TALDICE ice core from Talos Dome. We first review the different modern site characteristics in terms of ice flow, meteorological conditions, precipitation intermittency and moisture origin, as depicted by meteorological data, atmospheric reanalyses and Lagrangian moisture source diagnostics. These different factors can indeed alter the relationships between temperature and water stable isotopes. Using five records with sufficient resolution on the EDC3 age scale, common features are quantified through principal component analyses. Consistent with instrumental records and atmospheric model results, the ice core data depict rather coherent and homogenous patterns in East Antarctica during the last two interglacials. Across the East Antarctic plateau, regional differences, with respect to the common East Antarctic signal, appear to have similar patterns during the current and last interglacials. We identify two abrupt shifts in isotopic records during the glacial inception at TALDICE and EDML, likely caused by regional sea ice expansion. These regional differences are discussed in terms of moisture origin and in terms of past changes in local elevation histories, which are compared to ice sheet model results. Our results suggest that elevation changes may contribute significantly to inter-site differences. These elevation changes may be underestimated by current ice sheet models

Aeolian dust and chemistry records from the EDML and EDC ice cores

To improve quantitative interpretation of ice core aeolian dust records a systematic methodical comparison has been made involving methods of water-insoluble particle counting (Coulter Counter and laser-sensing particle detector), soluble ions (ion chromatography, IC, and continuous flow analysis, CFA), elemental analysis (inductively coupled plasma mass spectroscopy, ICP-MS, at pH 1 and after full acid digestion), and water-insoluble elemental analysis (proton induced X-ray emission, PIXE). Ice core samples covering the last deglaciation have been used from the EPICA Dome C (EDC) and the EPICA Dronning Maud Land (EDML) ice cores. All methods correlate very well amongst each other. The ratios of glacial age concentrations to Holocene concentrations, which are typically a factor ~100, differ significantly between the methods, but differences are limited to a factor < 2 for most methods with insoluble particles showing the largest change. The recovery of ICP-MS measurements depends on the digestion method and is different for different elements and during different climatic periods. EDC and EDML samples have similar dust composition, which suggests a common dust source or a common mixture of sources for the two sites. The analysed samples further reveal a change of dust composition during the last deglaciation.

Glaciological characteristics and carbon dioxide record of ice cores from Antarctica

We report a decadally resolved record of atmospheric CO2 concentration for the last 1000 years, obtained from the West Antarctic Ice Sheet (WAIS) Divide shallow ice core. The most prominent feature of the pre-industrial period is a rapid ~7 ppm decrease of CO2 in a span of ~20-50 years at ~1600 A.D. This observation confirms the timing of an abrupt atmospheric CO2 decrease of ~10 ppm observed for that time period in the Law Dome ice core CO2 records, but the true magnitude of the decrease remains unclear. Atmospheric CO2 variations over the time period 1000-1800 A.D. are statistically correlated with northern hemispheric climate and tropical Indo-Pacific sea surface temperature. However, the exact relationship between CO2 and climate remains elusive due to regional climate variations and/or uneven geographical data density of paleoclimate records. We observe small differences of 0 ~2% (0 ~ 6 ppm) among the high-precision CO2 records from the Law Dome, EPICA Dronning Maud Land and WAIS Divide Antarctic ice cores. However, those records share common trends of CO2 change on centennial to multicentennial time scales, and clearly show that atmospheric CO2 has been increasing above preindustrial levels since ~1850 A.D.

Abundance of plant macrofossils, mollusca and ostracoda in sediments obtained near Tessin, north Germany

A depression filled with Late Glacial and Holocene sediments was excavated during the geological exploration and recovery of a dump area near Tessin close to Rostock, and initiated the studies of the present paper. Pebble analysis of three exposed or respectively drilled till horizons as well as pollenanalytical, carpological and faunistical studies carried out allow the stratigraphical subdivision of the Quaternary sequence of the dump area. The basal till was probably the result of dead ice decay, and was lithostratigraphically assigned to the Pomerian Stage (qw2). The palynological results of boreholes RKS 19/93 and A/92 reveal pre-Allerod and other sediments instead of the expected interweichselian deposits. Based on the palynological and carpological findings, we correlated the beginning of the late glacial development in the locality with the end of the Meiendorf-lnterstadial sensu Menke in Bock et al. (1985, doi:10.3285/eg.35.1.18). The limnic-telmatic sedimentation could be observed pollen floristically probably starting with the Meiendorf-lnterstadial (Hippophae-Betula nana-phase) followed by the Bolling-(Betula nana-B. alba s.l.-Artemisia-Helianthemum-Poaceae-phase) and the Allerad-lnterstadial [Betula alba s.l.-(Pinus)-Cyperaceae-phase] lasting up to the Younger Dryas (Juniperus-Artemisia-Poaceae-phase). Sedimentation closed during the Younger Dryas with the accumulation of fine sands. It was reactivated later during the Holocene due to the anthropogene influence (Older and Younger Subatlantic, dampness of the depression by clearing).

Physico-chemical characteristics and ikaite content of brine and ice samples taken during SIPEX and DDU campaigns in East Antarctica, 2007

Calcium carbonate precipitation in sea ice is thought to potentially drive significant CO2 uptake by the ocean. However, little is known about the quantitative spatial and temporal distribution of CaCO3 within sea ice, although it is hypothesized that high quantities of dissolved organic matter and/or phosphate (common in sea ice) may inhibit its formation. In this quantitative study of hydrous calcium carbonate as ikaite, sea ice cores and brine samples were collected from pack and land fast sea ice between September and December 2007 during two expeditions, one in the East Antarctic sector and the other off Terre Adélie. Samples were analysed for CaCO3, salinity, dissolved organic carbon/nitrogen, inorganic phosphate, and total alkalinity. No relationship between these parameters and CaCO3 precipitation was evident. Ikaite was found mostly in the uppermost layers of sea ice with maximum concentrations of up to 126 mg ikaite per litre melted sea ice being measured, although both the temporal and horizontal spatial distributions of ikaite were highly heterogeneous. The precipitate was also found in the snow on top of the sea ice at some of the sampling locations.

Age determination, sedimentological, biogeochemical, and plant macrofossils analyses of outcrop PG2038-1

Thermokarst lakes are a widespread feature of the Arctic tundra, in which highly dynamic processes are closely connected with current and past climate changes. We investigated late Quaternary sediment dynamics, basin and shoreline evolution, and environmental interrelations of Lake El'gene-Kyuele in the NE Siberian Arctic (latitude 71°17'N, longitude 125°34'E). The water-body displays thaw-lake characteristics cutting into both Pleistocene Ice Complex and Holocene alas sediments. Our methods are based on grain size distribution, mineralogical composition, TOC/N ratio, stable carbon isotopes and the analysis of plant macrofossils from a 3.5-m sediment profile at the modern eastern lake shore. Our results show two main sources for sediments in the lake basin: terrigenous diamicton supplied from thermokarst slopes and the lake shore, and lacustrine detritus that has mainly settled in the deep lake basin. The lake and its adjacent thermokarst basin rapidly expanded during the early Holocene. This climatically warmer than today period was characterized by forest or forest tundra vegetation composed of larches, birch trees and shrubs. Woodlands of both the HTM and the Late Pleistocene were affected by fire, which potentially triggered the initiation of thermokarst processes resulting later in lake formation and expansion. The maximum lake depth at the study site and the lowest limnic bioproductivity occurred during the longest time interval of ~7 ka starting in the Holocene Thermal Maximum and lasting throughout the progressively cooler Neoglacial, whereas partial drainage and an extensive shift of the lake shoreline occurred ~0.9 cal. ka BP. Correspondingly, this study discusses different climatic and environmental drivers for the dynamics of a thermokarst basin.

Isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) of sediment core SO188_342 from the Northern Bay of Bengal for the last 18 ka

The Indian Summer Monsoon (ISM) is a major global climatic phenomenon. Long-term precipitation proxy records of the ISM, however, are often fragmented and discontinuous, impeding an estimation of the magnitude of precipitation variability from the Last Glacial to the present. To improve our understanding of past ISM variability, we provide a continuous reconstructed record of precipitation and continental vegetation changes from the lower Ganges-Brahmaputra-Meghna catchment and the Indo-Burman ranges over the last 18,000 years (18 ka). The records derive from a marine sediment core from the northern Bay of Bengal (NBoB), and are complemented by numerical model results of spatial moisture transport and precipitation distribution over the Bengal region. The isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) are compared to results from an isotope-enabled general atmospheric circulation model (IsoCAM) for selected time slices (pre-industrial, mid-Holocene and Heinrich Stadial 1). Comparison of proxy and model results indicate that past changes in the dD of precipitation and plant waxes were mainly driven by the amount effect, and strongly influenced by ISM rainfall. Maximum precipitation is detected for the Early Holocene Climatic Optimum (EHCO; 10.5-6 ka BP), whereas minimum precipitation occurred during the Heinrich Stadial 1 (HS1; 16.9-15.4 ka BP). The IsoCAM model results support the hypothesis of a constant moisture source (i.e. the NBoB) throughout the study period. Relative to the pre-industrial period the model reconstructions show 20% more rain during the mid-Holocene (6 ka BP) and 20% less rain during the Heinrich Stadial 1 (HS1), respectively. A shift from C4-plant dominated ecosystems during the glacial to subsequent C3/C4-mixed ones during the interglacial took place. Vegetation changes were predominantly driven by precipitation variability, as evidenced by the significant correlation between the dD and d13C alkane records. When compared to other records across the ISM domain, precipitation and vegetation changes inferred from our records and the numerical model results provide evidence for a coherent regional variability of the ISM from the Last Glacial to the present.

Plant-wax, alkenone and major element analyses of ODP Site 659 for the Last Glacial cycle and associated coretop data

We present a hydrologic reconstruction of the Sahara-Sahel transition, covering the complete last glacial cycle (130 ka), based on a combination of plant-wax-specific hydrogen (dD) and carbon isotopes (d13C). The dD and d13C signatures of long-chain n-alkanes from ODP Site 659 off NW Africa reveal a significant anti-correlation. Complementary to published pollen data, we infer that this plant-wax signal reflects sensitive responses of the vegetation cover to precipitation changes in the Sahel region, as well as varying contributions from biomes north of the Sahara (C3 domain) by North-East Trade Winds (NETW). During arid phases, especially the northern parts of the Sahel likely experienced crucial water stress, which resulted in a pronounced contraction of the vegetation cover, thus reducing the amount of C4 plant waxes from the region. The increase in NETW strength during dry periods further promoted a more pronounced C3-plant-wax signal derived from the North African C3 plant domain. During humid periods, the C4-dominated Sahelian environments spread northward into the Saharan realm, in association with lower NETW inputs of C3 plant waxes. Arid-humid cycles deduced from plant-wax dD are in accordance with concomitant changes in weathering intensity reflected in varying major element distributions. Environmental shifts are generally linked to periods with large fluctuations in Northern Hemisphere summer insolation. During Marine Isotope Stages 2 and 3, when insolation variability was low, coupling of the hydrologic regime to alkenone-based estimates of NE Atlantic sea-surface temperatures becomes apparent.

Synchronization of the NGRIP, GRIP, and GISP2 ice cores

We here present a synchronization of the NGRIP, GRIP, and GISP2 ice cores based mainly on volcanic events over the period 14.9-32.45 ka b2k (before AD 2000), corresponding to Marine Isotope Stage 2 (MIS 2) and the end of MIS 3. The matching provides a basis for applying the recent NGRIP-based Greenland Ice Core Chronology 2005 (GICC05) time scale to the GRIP and GISP2 ice cores, thereby making it possible to compare the synchronized palaeoclimate profiles of the cores in detail and to identify relative accumulation differences between the cores. Based on the matching, a period of anomalous high accumulation rates in the GISP2 ice core is detected within the period 16.5-18.3 ka b2k. The d18O and [Ca2+] profiles of the three cores are presented on the common GICC05 time scale and generally show excellent agreement across the stadial-interstadial transitions and across the two characteristic dust events in Greenland Stadial 3. However, large differences between the d18O and [Ca2+] profiles of the three cores are seen in the same period as the 7-9% increase in the GISP2 accumulation rate. We conclude that changes of the atmospheric circulation are likely to have occurred in this period, altering the spatial gradients in Greenland and resulting in larger variations between the records.

Carbon/Nitrogen ratio of plant tissues and biomass of belowground consumers after 11 years of experimental treatment, Toolik lake

Theory and observation indicate that changes in the rate of primary production can alter the balance between the bottom-up influences of plants and resources and the top-down regulation of herbivores and predators on ecosystem structure and function. The Exploitation Ecosystem Hypothesis (EEH) posited that as aboveground net primary productivity (ANPP) increases, the additional biomass should support higher trophic levels. We developed an extension of EEH to include the impacts of increases in ANPP on belowground consumers in a similar manner as aboveground, but indirectly through changes in the allocation of photosynthate to roots. We tested our predictions for plants aboveground and for phytophagous nematodes and their predators belowground in two common arctic tundra plant communities subjected to 11 years of increased soil nutrient availability and/or exclusion of mammalian herbivores. The less productive dry heath (DH) community met the predictions of EEH aboveground, with the greatest ANPP and plant biomass in the fertilized plots protected from herbivory. A palatable grass increased in fertilized plots while dwarf evergreen shrubs and lichens declined. Belowground, phytophagous nematodes also responded as predicted, achieving greater biomass in the higher ANPP plots, whereas predator biomass tended to be lower in those same plots (although not significantly). In the higher productivity moist acidic tussock (MAT) community, aboveground responses were quite different. Herbivores stimulated ANPP and biomass in both ambient and enriched soil nutrient plots; maximum ANPP occurred in fertilized plots exposed to herbivory. Fertilized plots became dominated by dwarf birch (a deciduous shrub) and cloudberry (a perennial forb); under ambient conditions these two species coexist with sedges, evergreen dwarf shrubs, and Sphagnum mosses. Phytophagous nematodes did not respond significantly to changes in ANPP, although predator biomass was greatest in control plots. The contrasting results of these two arctic tundra plant communities suggest that the predictions of EEH may hold for very low ANPP communities, but that other factors, including competition and shifts in vegetation composition toward less palatable species, may confound predicted responses to changes in productivity in higher ANPP communities such as the MAT studied here.

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.