Duration of Greenland Stadial 22 and ice-gas Δage from counting of annual layers in Greenland NGRIP ice core

High-resolution measurements of chemical impurities and methane concentrations in Greenland ice core samples from the early glacial period allow the extension of annual-layer counted chronologies and the improvement of gas age-ice age difference (Δage) essential to the synchronization of ice core records. We report high-resolution measurements of a 50 m section of the NorthGRIP ice core and corresponding annual layer thicknesses in order to constrain the duration of the Greenland Stadial 22 (GS-22) between Greenland Interstadials (GIs) 21 and 22, for which inconsistent durations and ages have been reported from Greenland and Antarctic ice core records as well as European speleothems. Depending on the chronology used, GS-22 occurred between approximately 89 (end of GI-22) and 83 kyr b2k (onset of GI-21). From annual layer counting, we find that GS-22 lasted between 2696 and 3092 years and was followed by a GI-21 pre-cursor event lasting between 331 and 369 yr. Our layer-based counting agrees with the duration of stadial 22 as determined from the NALPS speleothem record (3250 ± 526 yr) but not with that of the GICC05modelext chronology (2620 yr) or an alternative chronology based on gas-marker synchronization to EPICA Dronning Maud Land ice core. These results show that GICC05modelext overestimates accumulation and/or underestimates thinning in this early part of the last glacial period. We also revise the possible ranges of NorthGRIP Δdepth (5.49 to 5.85 m) and Δage (498 to 601 yr) at the warming onset of GI-21 as well as the Δage range at the onset of the GI-21 precursor warming (523 to 654 yr), observing that temperature (represented by the δ15N proxy) increases before CH4 concentration by no more than a few decades.

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.

Glacier mass balance reconstruction by sublimation induced enrichment of chemical species on Cerro Tapado (Chilean Andes)

A 36 m long ice core down to bedrock from the Cerro Tapado glacier (5536 m a.s.l, 30°08' S, 69°55' W) was analyzed to reconstruct past climatic conditions for Northern Chile. Because of the marked seasonality in the precipitation (short wet winter and extended dry summer periods) in this region, major snow ablation and related post-depositional processes occur on the glacier surface during summer periods. They include predominantly sublimation and dry deposition. Assuming that, like measured during the field campaign, the enrichment of chloride was always related to sublimation, the chemical record along the ice core may be applied to reconstruct the history of such secondary processes linked to the past climatic conditions over northern Chile. For the time period 1962–1999, a mean annual net accumulation of 316 mm water equivalent (weq) and 327 mm weq loss by sublimation was deduced by this method. This corresponds to an initial total annual accumulation of 539 mm weq. The annual variability of the accumulation and sublimation is related with the Southern Oscillation Index (SOI): higher net-accumulation during El-Niño years and more sublimation during La Niña years. The deepest part of the ice record shows a time discontinuity; with an ice body deposited under different climatic conditions: 290 mm higher precipitation but with reduced seasonal distribution (+470 mm in winter and –180 mm in summer) and –3°C lower mean annual temperature. Unfortunately, its age is unknown. The comparison with regional proxy data however let us conclude that the glacier buildup did most likely occur after the dry mid-Holocene.

Radiocarbon analysis in an Alpine ice core: Record of anthropogenic and biogenic contributions to carbonaceous Radiocarbon analysis in an Alpine ice core: Record of anthropogenic and biogenic contributions to carbonaceous aerosols in the past (1650-1940)

Long-term concentration records of carbonaceous particles (CP) are of increasing interest in climate research due to their not yet completely understood effects on climate. Nevertheless, only poor data on their concentrations and sources before the 20th century are available. We present a first long-term record of organic carbon (OC) and elemental carbon (EC) concentrations – the two main fractions of CP – along with the corresponding fraction of modern carbon (fM) derived from radiocarbon (14C) analysis in ice. This allows a distinction and quantification of natural (biogenic) and anthropogenic (fossil) sources in the past. CP were extracted from an ice archive, with resulting carbon quantities in the microgram range. Analysis of 14C by accelerator mass spectrometry (AMS) was therefore highly demanding. We analysed 33 samples of 0.4 to 1 kg ice from a 150.5 m long ice core retrieved at Fiescherhorn glacier in December 2002 (46°33'3.2" N, 08°04'0.4" E; 3900 m a.s.l.). Samples were taken from bedrock up to the firn/ice transition, covering the time period 1650–1940 and thus the transition from the pre-industrial to the industrial era. Before ~1850, OC was approaching a purely biogenic origin with a mean concentration of 24 μg kg−1 and a standard deviation of 7 μg kg−1. In 1940, OC concentration was about a factor of 3 higher than this biogenic background, almost half of it originating from anthropogenic sources, i.e. from combustion of fossil fuels. The biogenic EC concentration was nearly constant over the examined time period with 6 μg kg−1 and a standard deviation of 1 μg kg−1. In 1940, the additional anthropogenic input of atmospheric EC was about 50 μg kg−1.

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.

Synchronisation of the EDML and EDC ice cores for the last 52 kyr by volcanic signature matching

A common time scale for the EPICA ice cores from Dome C (EDC) and Dronning Maud Land (EDML) has been established. Since the EDML core was not drilled on a dome, the development of the EDML1 time scale for the EPICA ice core drilled in Dronning Maud Land was based on the creation of a detailed stratigraphic link between EDML and EDC, which was dated by a simpler 1D ice-flow model. The synchronisation between the two EPICA ice cores was done through the identification of several common volcanic signatures. This paper describes the rigorous method, using the signature of volcanic sulfate, which was employed for the last 52 kyr of the record. We estimated the discrepancies between the modelled EDC and EDML glaciological age scales during the studied period, by evaluating the ratio R of the apparent duration of temporal intervals between pairs of isochrones. On average R ranges between 0.8 and 1.2 corresponding to an uncertainty of up to 20% in the estimate of the time duration in at least one of the two ice cores. Significant deviations of R up to 1.4–1.5 are observed between 18 and 28 kyr before present (BP), where present is defined as 1950. At this stage our approach does not allow us unequivocally to find out which of the models is affected by errors, but assuming that the thinning function at both sites and accumulation history at Dome C (which was drilled on a dome) are correct, this anomaly can be ascribed to a complex spatial accumulation variability (which may be different in the past compared to the present day) upstream of the EDML core.

The Southern Hemisphere at glacial terminations: insights from the Dome C ice core

The many different proxy records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core allow for the first time a comparison of nine glacial terminations in great detail. Despite the fact that all terminations cover the transition from a glacial maximum into an interglacial, there are large differences between single terminations. For some terminations, Antarctic temperature increased only moderately, while for others, the amplitude of change at the termination was much larger. For the different terminations, the rate of change in temperature is more similar than the magnitude or duration of change. These temperature changes were accompanied by vast changes in dust and sea salt deposition all over Antarctica. Here we investigate the phasing between a South American dust proxy (non-sea-salt calcium flux, nssCa2+), a sea ice proxy (sea salt sodium flux, ssNa+) and a proxy for Antarctic temperature (deuterium, δD). In particular, we look into whether a similar sequence of events applies to all terminations, despite their different characteristics. All proxies are derived from the EPICA Dome C ice core, resulting in a relative dating uncertainty between the proxies of less than 20 years. At the start of the terminations, the temperature (δD) increase and dust (nssCa2+ flux) decrease start synchronously. The sea ice proxy (ssNa+ flux), however, only changes once the temperature has reached a particular threshold, approximately 5°C below present day temperatures (corresponding to a δD value of −420‰). This reflects to a large extent the limited sensitivity of the sea ice proxy during very cold periods with large sea ice extent. At terminations where this threshold is not reached (TVI, TVIII), ssNa+ flux shows no changes. Above this threshold, the sea ice proxy is closely coupled to the Antarctic temperature, and interglacial levels are reached at the same time for both ssNa+ and δD. On the other hand, once another threshold at approximately 2°C below present day temperature is passed (corresponding to a δD value of −402‰), nssCa2+ flux has reached interglacial levels and does not change any more, despite further warming. This threshold behaviour most likely results from a combination of changes to the threshold friction velocity for dust entrainment and to the distribution of surface wind speeds in the dust source region.

Climatic and human impacts on mountain vegetation at Lauenensee (Bernese Alps, Switzerland) during the last 14,000 years

Lake sediments from Lauenensee (1381 m a.s.l.), a small lake in the Bernese Alps, were analysed to reconstruct the vegetation and fire history. The chronology is based on 11 calibrated radiocarbon dates on terrestrial plant macrofossils suggesting a basal age of 14,200 cal. BP. Pollen and macrofossil data imply that treeline never reached the lake catchment during the Bølling–Allerød interstadial. Treeline north of the Alps was depressed by c. 300 altitudinal meters, if compared with southern locations. We attribute this difference to colder temperatures and to unbuffered cold air excursions from the ice masses in northern Europe. Afforestation started after the Younger Dryas at 11,600 cal. BP. Early-Holocene tree-Betula and Pinus sylvestris forests were replaced by Abies alba forests around 7500 cal. BP. Continuous high-resolution pollen and macrofossil series allow quantitative assessments of vegetation dynamics at 5900–5200 cal. BP (first expansion of Picea abies, decline of Abies alba) and 4100–2900 cal. BP (first collapse of Abies alba). The first signs of human activity became noticeable during the late Neolithic c. 5700–5200 cal. BP. Cross-correlation analysis shows that the expansion of Alnus viridis and the replacement of Abies alba by Picea abies after c. 5500 cal. BP was most likely a consequence of human disturbance. Abies alba responded very sensitively to a combination of fire and grazing disturbance. Our results imply that the current dominance of Picea abies in the upper montane and subalpine belts is a consequence of anthropogenic activities through the millennia.

Simulating the temperature and precipitation signal in an Alpine ice core

Accumulation and delta O-18 data from Alpine ice cores provide information on past temperature and precipitation. However, their correlation with seasonal or annual mean temperature and precipitation at nearby sites is often low. This is partly due to the irregular sampling of the atmosphere by the ice core (i.e. ice cores almost only record precipitation events and not dry periods) and the possible incongruity between annual layers and calendar years. Using daily meteorological data from a nearby station and reanalyses, we replicate the ice core from the Grenzgletscher (Switzerland, 4200m a.s.l.) on a sample-by-sample basis by calculating precipitation-weighted temperature (PWT) over short intervals. Over the last 15 yr of the ice core record, accumulation and delta O-18 variations can be well reproduced on a sub-seasonal scale. This allows a wiggle-matching approach for defining quasi-annual layers, resulting in high correlations between measured quasi-annual delta O-18 and PWT. Further back in time, the agreement deteriorates. Nevertheless, we find significant correlations over the entire length of the record (1938-1993) of ice core delta O-18 with PWT, but not with annual mean temperature. This is due to the low correlations between PWT and annual mean temperature, a characteristic which in ERA-Interim reanalysis is also found for many other continental mid-to-high-latitude regions. The fact that meteorologically very different years can lead to similar combinations of PWT and accumulation poses limitations to the use of delta O-18 from Alpine ice cores for temperature reconstructions. Rather than for reconstructing annual mean temperature, delta O-18 from Alpine ice cores should be used to reconstruct PWT over quasi-annual periods. This variable is reproducible in reanalysis or climate model data and could thus be assimilated into conventional climate models.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

The onset of Neoglaciation 6000 years ago in western Mongolia revealed by an ice core from the Tsambagarav mountain range

Glacier highstands since the Last Glacial Maximum are well documented for many regions, but little is known about glacier fluctuations and lowstands during the Holocene. This is because the traces of minimum extents are difficult to identify and at many places are still ice covered, limiting the access to sample material. Here we report a new approach to assess minimal glacier extent, using a 72-m long surface-to-bedrock ice core drilled on Khukh Nuru Uul, a glacier in the Tsambagarav mountain range of the Mongolian Altai (4130 m asl, 48°39.338′N, 90°50.826′E). The small ice cap has low ice temperatures and flat bedrock topography at the drill site. This indicates minimal lateral glacier flow and thereby preserved climate signals. The upper two-thirds of the ice core contain 200 years of climate information with annual resolution, whereas the lower third is subject to strong thinning of the annual layers with a basal ice age of approximately 6000 years before present (BP). We interpret the basal ice age as indicative of ice-free conditions in the Tsambagarav mountain range at 4100 m asl prior to 6000 years BP. This age marks the onset of the Neoglaciation and the end of the Holocene Climate Optimum. The ice-free conditions allow for adjusting the Equilibrium Line Altitude (ELA) and derive the glacier extent in the Mongolian Altai during the Holocene Climate Optimum. Based on the ELA-shift, we conclude that most of the glaciers are not remnants of the Last Glacial Maximum but were formed during the second part of the Holocene. The ice core derived accumulation reconstruction suggests important changes in the precipitation pattern over the last 6000 years. During formation of the glacier, more humid conditions than presently prevailed followed by a long dry period from 5000 years BP until 250 years ago. Present conditions are more humid than during the past millennia. This is consistent with precipitation evolution derived from lake sediment studies in the Altai.

Observations of the northern seasonal polar cap on Mars: I. Spring sublimation activity and processes

Spring sublimation of the seasonal CO2 northern polar cap is a dynamic process in the current Mars climate. Phenomena include dark fans of dune material propelled out onto the seasonal ice layer, polygonal cracks in the seasonal ice, sand flow down slipfaces, and outbreaks of gas and sand around the dune margins. These phenomena are concentrated on the north polar erg that encircles the northern residual polar cap. The Mars Reconnaissance Orbiter has been in orbit for three Mars years, allowing us to observe three northern spring seasons. Activity is consistent with and well described by the Kieffer model of basal sublimation of the seasonal layer of ice applied originally in the southern hemisphere. Three typical weak spots have been identified on the dunes for escape of gas sublimed from the bottom of the seasonal ice layer: the crest of the dune, the interface of the dune with the interdune substrate, and through polygonal cracks in the ice. Pressurized gas flows through these vents and carries out material entrained from the dune. Furrows in the dunes channel gas to outbreak points and may be the northern equivalent of southern radially-organized channels ("araneiform" terrain), albeit not permanent. Properties of the seasonal CO2 ice layer are derived from timing of seasonal events such as when final sublimation occurs. Modification of dune morphology shows that landscape evolution is occurring on Mars today, driven by seasonal activity associated with sublimation of the seasonal CO2 polar cap.