S. L. Buchardt

A synchronized dating of three Greenland ice cores throughout the Holocene

As part of the effort to create the new Greenland Ice Core Chronology 2005 (GICC05) a synchronized stratigraphical timescale for the Holocene parts of the DYE- 3, GRIP and NGRIP ice cores is made by using volcanic reference horizons in electri- cal conductivity measurements to match the cores. The main annual layer counting is carried out on the most suited records only, exploit- ing that the three ice cores have been drilled at locations with different climatic con- ditions and differences in ice flow. However, supplemental counting on data from all cores has been performed between each set of reference horizons in order to verify the valid- ity of the match. After the verification, the main dating is transferred to all records us- ing the volcanic reference horizons as tie points. An assessment of the mean annual layer thickness in each core section confirms that the new synchronized dating is consistent for all three cores. The data used for the main annual layer counting of the past 7900 years are the DYE- 3, GRIP and NGRIP stable isotope records. As the high accumulation rate at the DYE- 3 drill site makes the seasonal cycle in the DYE-3 stable isotopes very resistant to firn diffusion, an effort has been made to extend the DYE-3 Holocene record. The new syn- chronized dating relies heavily on this record of �75,000 stable isotope samples. The dat- ing of the early Holocene consists of an already established part of GICC05 for GRIP and NGRIP which has now been transferred to the DYE-3 core. GICC05 dates the Younger Dryas termination, as defined from deuterium excess, to 11,703 b2k; 130 years earlier than the previous GRIP dating.

Holocene thinning of the Greenland ice sheet

On entering an era of global warming, the stability of the Greenland ice sheet (GIS) is an important concern, especially in the light of new evidence of rapidly changing flow and melt conditions at the GIS margins. Studying the response of the GIS to past climatic change may help to advance our understanding of GIS dynamics. The previous interpretation of evidence from stable isotopes (δ18O) in water from GIS ice cores was that Holocene climate variability on the GIS differed spatially and that a consistent Holocene climate optimum—the unusually warm period from about 9,000 to 6,000 years ago found in many northern-latitude palaeoclimate records—did not exist. Here we extract both the Greenland Holocene temperature history and the evolution of GIS surface elevation at four GIS locations. We achieve this by comparing δ18O from GIS ice cores with δ18O from ice cores from small marginal icecaps. Contrary to the earlier interpretation of δ18O evidence from ice cores, our new temperature history reveals a pronounced Holocene climatic optimum in Greenland coinciding with maximum thinning near the GIS margins. Our δ18O-based results are corroborated by the air content of ice cores, a proxy for surface elevation. State-of-the-art ice sheet models are generally found to be underestimating the extent and changes in GIS elevation and area; our findings may help to improve the ability of models to reproduce the GIS response to Holocene climate.

Estimating the basal melt rate at NorthGRIP using a Monte Carlo technique

Abstract From radio-echo sounding (RES) surveys and ice core data it can be seen that the ice sheet is melting at the base in a large area in Northern Greenland. The RES images reveal internal layers in the ice. The layers are former deposition surfaces and are thus isochrones. Undulations of the isochrones in regions where the base is smooth suggest that the basal melt rate changes over short distances. This indicates that the geothermal heat flux is very high and has large spatial variability in Northern Greenland. In this study, the basal melt rate at the NorthGRIP drill site in North-Central Greenland is calculated by inverse modelling. We use simple one- and two-dimensional flow models to simulate the ice flow along the NNW-trending ice ridge leading to NorthGRIP. The accumulation is calculated from a dynamical model. Several ice flow parameters are unknown and must be estimated along with the basal melt rate using a Monte Carlo method. The Monte Carlo inversion is constrained by the observed isochrones, dated from the timescale established for the NorthGRIP ice core. The estimates of the basal melt rates around NorthGRIP are obtained from both the one- and two-dimensional models. Combining the estimated basal melt rates with the observed borehole temperatures allows us to convert the basal melt rates to geothermal heat flow values. From the two-dimensional model we find the basal melt rate and geothermal heat flux at NorthGRIP to be 6.1 mma–1 and 129 mWm–2, respectively.

Water isotope diffusion rates from the NorthGRIP ice core for the last 16,000 years – Glaciological and paleoclimatic implications

Abstract A high resolution (0.05 m) water isotopic record ( δ O 18 ) is available from the NorthGRIP ice core. In this study we look into the water isotope diffusion history as estimated by the spectral characteristics of the δ O 18 time series covering the last 16,000 years. The diffusion of water vapor in the porous medium of the firn pack attenuates the initial isotopic signal, predominantly having an impact on the high frequency components of the power spectrum. Higher temperatures induce higher rates of smoothing and thus the signal can be used as a firn paleothermometer. We use a water isotope diffusion model coupled to a steady-state densification model in order to infer the temperature signal from the site, assuming the accumulation and strain rate history as estimated using the GICC05 layer counted chronology and a Dansgaard–Johnsen ice flow model. The temperature reconstruction accurately captures the timing and magnitude of the Bolling–Allerod and Younger Dryas transitions. A Holocene climatic optimum is seen between 7 and 9 ky b2k with an apparent cooling trend thereafter. Our temperature estimate for the Holocene climatic optimum, points to a necessary adjustment of the ice thinning function indicating that the ice flow model overestimates past accumulation rates by about 10% at 8 ky b2k. This result, is supported by recent gas isotopic fractionation studies proposing a similar reduction for glacial conditions. Finally, the record presents a climatic variability over the Holocene spanning millennial and centennial scales with a profound cooling occurring at approximately 4000 years b2k. The new reconstruction technique is able to provide past temperature estimates by overcoming the issues apparent in the use of the classical δ O 18 slope method. It can at the same time resolve temperature signals at low and high frequencies.

At what depth is the Eemian layer expected to be found at NEEM

Abstract No continuous record from Greenland of the Eemian interglacial period (130–115 ka BP) currently exists. However, a new ice-core drill site has been suggested at 77.449 ˚ N, 51.056 ˚ Win north-west Greenland (North Eemian or NEEM). Radio-echo sounding images and flow model investigations indicate that an undisturbed Eemian record may be obtained at NEEM. In this work, a two-dimensional ice flow model with time-dependent accumulation rate and ice thickness is used to estimate the location of the Eemian layer at the new drill site. The model is used to simulate the ice flow along the ice ridge leading to the drill site. Unknown flow parameters are found through a Monte Carlo analysis of the flow model constrained by observed isochrones in the ice. The results indicate that the Eemian layer is approximately 60m thick and that its base is located approximately 100m above bedrock.