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.

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.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Diffusive equilibration of N2, O2 and CO2 mixing ratios in a 1.5-million-years-old ice core

In the framework of the International Partnerships in Ice Core Sciences, one of the most important targets is to retrieve an Antarctic ice core that extends over the last 1.5 million years (i.e. an ice core that enters the climate era when glacial–interglacial cycles followed the obliquity cycles of the earth). In such an ice core the annual layers of the oldest ice would be thinned by a factor of about 100 and the climatic information of a 10 000 yr interval would be contained in less than 1 m of ice. The gas record in such an Antarctic ice core can potentially reveal the role of greenhouse gas forcing on these 40 000 yr cycles. However, besides the extreme thinning of the annual layers, also the long residence time of the trapped air in the ice and the relatively high ice temperatures near the bedrock favour diffusive exchanges. To investigate the changes in the O2 / N2 ratio, as well as the trapped CO2 concentrations, we modelled the diffusive exchange of the trapped gases O2, N2 and CO2 along the vertical axis. However, the boundary conditions of a potential drilling site are not yet well constrained and the uncertainties in the permeation coefficients of the air constituents in the ice are large. In our simulations, we have set the drill site ice thickness at 2700 m and the bedrock ice temperature at 5–10 K below the ice pressure melting point. Using these conditions and including all further uncertainties associated with the drill site and the permeation coefficients, the results suggest that in the oldest ice the precessional variations in the O2 / N2 ratio will be damped by 50–100%, whereas CO2 concentration changes associated with glacial–interglacial variations will likely be conserved (simulated damping 5%). If the precessional O2 / N2 signal will have disappeared completely in this future ice core, orbital tuning of the ice-core age scale will be limited.

Direct linking of Greenland and Antarctic ice cores at the Toba eruption (74 ka BP)

The Toba eruption that occurred some 74 ka ago in Sumatra, Indonesia, is among the largest volcanic events on Earth over the last 2 million years. Tephra from this eruption has been spread over vast areas in Asia, where it constitutes a major time marker close to the Marine Isotope Stage 4/5 boundary. As yet, no tephra associated with Toba has been identified in Greenland or Antarctic ice cores. Based on new accurate dating of Toba tephra and on accurately dated European stalagmites, the Toba event is known to occur between the onsets of Greenland interstadials (GI) 19 and 20. Furthermore, the existing linking of Greenland and Antarctic ice cores by gas records and by the bipolar seesaw hypothesis suggests that the Antarctic counterpart is situated between Antarctic Isotope Maxima (AIM) 19 and 20. In this work we suggest a direct synchronization of Greenland (NGRIP) and Antarctic (EDML) ice cores at the Toba eruption based on matching of a pattern of bipolar volcanic spikes. Annual layer counting between volcanic spikes in both cores allows for a unique match. We first demonstrate this bipolar matching technique at the already synchronized Laschamp geomagnetic excursion (41 ka BP) before we apply it to the suggested Toba interval. The Toba synchronization pattern covers some 2000 yr in GI-20 and AIM-19/20 and includes nine acidity peaks that are recognized in both ice cores. The suggested bipolar Toba synchronization has decadal precision. It thus allows a determination of the exact phasing of inter-hemispheric climate in a time interval of poorly constrained ice core records, and it allows for a discussion of the climatic impact of the Toba eruption in a global perspective. The bipolar linking gives no support for a long-term global cooling caused by the Toba eruption as Antarctica experiences a major warming shortly after the event. Furthermore, our bipolar match provides a way to place palaeo-environmental records other than ice cores into a precise climatic context.

Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.