The North Greenland Icecore Project (NorthGRIP) was initiated in 1995 as a joint international programme involving Denmark, Germany, Japan, Belgium, Sweden, Iceland, the U.S.A., France and Switzerland. the main goal was to obtain undisturbed high-resolution information about the Eemian climatic period (115–130 kyr BP). the records from the Greenland Icecore Project (GRIP) and Greenland Ice Sheet Project 2 (GISP2) in central Greenland are different and disturbed down in the ice covering this period. Internal radio-echo sounding layers show that NorthGRIP, placed 325 km north-northwest of GRIP at the Summit of the Greenland ice sheet, is located on a gently sloping ice ridge with very flat bedrock and internal layers found so high that an undisturbed Eemian record is possible. Internal layers much farther above bedrock than their apparent counter parts at GRIP suggest that conditions are favourable for recovery of an undisturbed Eemian record. So far, a 1351 mdeep ice core (NorthGRIP1) and a 3001 mdeep ice core (NorthGRIP 2) have been recovered. the ice thickness is expected to be 3080 m, and the ice temperature at 3001 m is –5.6°C, so we expect basal melting at the bedrock. Most of the Eemian ice will be melted away, leaving only the last part and the transition between the Eem and the Last Glacial Period. At 3001 m the age of the ice is 110 kyr BP and the annual layers are of the order 1 cm.With modern methods the annual layers can be resolved, resulting in detailed information on the decline of the warm Eemian period into the Last Glacial Period.

The ice-core logging procedure used to log the North Greenland Icecore Project (NorthGRIP) ice cores is described. the existence of two deep ice cores, NorthGRIP 1and 2, drilled 25 mapart, allows an independent evaluation of the procedure. the logged depths of the NorthGRIP 1 and 2 cores deviate from the length of the unwound drill cable corrected for elongation of the cable when hanging in the hole, by 1.5‰ and 50.5‰ at depths of 1371 and 2931 m, respectively. Differences between logged depths of identified layers found in both cores are studied in the depth interval where they overlap. Layers are identified by electrical conductivity measurements and dielectric profiling. the difference between the logged depths of layers identified in both cores increases to 0.5 m close to the bottom of the NorthGRIP 1 core, which is <0.5 mm m–1 ice core. the comparison between the two cores shows that the NorthGRIP logging procedure is accurate and reproducible. Further, our results show that the temperature conditions and handling of the core during logging are important for obtaining a precise depth.

Near-surface firn temperatures were measured in 22 steam-drilled boreholes in the summit region of Mont Blanc (France and Italy) at 3800 –4800ma.s.l. in June 1998 and in 31 boreholes in the Monte Rosa area (Italy and Switzerland) at 3900– 4500ma.s.l. in May–July 1999. Borehole temperatures were logged to 22 m depth. the temperatures at 18 m depth ranged between temperate conditions and approximately –15˚C. In a small altitude band, the observed distribution pattern suggests a strong influence of shortwave radiation and turbulent heat exchange (being generally more effective at wind-exposed sites). These two energy fluxes mainly determine the melt-energy input into the snow and firn during summer and, thereby, the measured near-surface temperatures. A statistical analysis of the measured firn temperatures revealed altitude-dependent firn temperature gradients of –1.48 and of –2.36˚C (100m)–1 for the Mont Blanc and Monte Rosa areas, respectively. the high lapse rates, as compared to the air-temperature lapse rate, are the result of englacial latent-heat contribution. the parameters elevation, potential direct solar radiation, slope and accumulation explain 480% of the variation of the mean annual firn temperatures. Aspect-dependent lower boundaries for the cold-firn occurrence in the two areas ranged between 3500 and 4100 ma.s.l.

An early study of the various components of the Greenland, Antarctic and Canadian Arctic ice-cap cores (Koerner, 1989) suggested that during the last interglacial period, the Greenland ice sheet suffered massive retreat and Canadian ice caps melted completely. Since then, modeling has helped support this interpretation (Cuffey and Marshall, 2000). Ice-core records of stable isotopes, melt layering and chemistry from the same Canadian ice cores, and others from the Russian Arctic islands, Svalbard and Greenland are presented as evidence for a more modest, but still substantial, retreat in the early Holocene. the sections representing the first half of the Holocene in many cores have less negative δ18O values (d values) and a higher percentage of melt layers than recently deposited ice, suggesting that temperatures were 1.3–3.5˚C warmer than today. Given that glacier balances are slightly negative today, they must have been substantially more negative during the early-Holocene thermal maximum, leading to retreat of the circumpolar ice caps. Evidence is presented to suggest that, with the exception of Academii Nauk ice cap, the ice in the Russian Arctic islands and Svalbard must have almost disappeared. In the Canadian Arctic, the larger Canadian ice caps retreated but survived. the cooling trend that followed this thermal maximum promoted re-expansion and new growth of most of the ice caps in the Russian Arctic islands and Svalbard.

The paper presents first results from the upper 54m of a 723.91m ice core drilled on Akademii Nauk ice cap, Severnaya Zemlya, Eurasian Arctic, in 1999– 2001, supplemented by data from shallow ice cores. the glacier’s peculiarity is the infiltration and refreezing of meltwater, which changes the original isotopic and chemical signals. Therefore, stratigraphical observations in these ice cores are more difficult than in those from central Greenland or Antarctica. However, the 1963 maximum of artificial radioactivity from atmospheric nuclear tests is clearly detectable in the deep ice core, and the δ18O profile of a 12.82 m shallow core shows annual variations. Consequently, at least for the upper part of the main core, an almost seasonal time resolution of palaeoclimate record could be expected. the Chernobyl layer is detected by increased 137Cs activity at depths of 11.81–12.51m related to the AD 2000 surface. the resulting mean annual net mass balance is 53±2 g cm–2 a–1. Data from dielectric profiling of the main core show considerable peaks in conductivity; one of them is interpreted as a volcano event. According to the resulting chronology, this part of the core represents approximately the last 100 years.

In order to reconstruct climatic and environmental changes in the Canadian Arctic, an 85 m deep ice core drilled in 1995 on Penny Ice Cap, Baffin Island, was analyzed for ions and δ18O. In addition to the core, snow-pit samples collected in 1994 and 1995 were also analyzed. Elution of ions caused by summer melting was observed in the pits. Due to the heavy summer melting on this ice cap, seasonal variations of ion chemistry and δ18O were not always present in the core. Comparisons of this core with a previously reported core drilled 2.5 maway show that the noise contained in single annual time series is 40–50% for ions and 25% for δ18O. the ice-core data, however, provide us with a reasonable proxy record of climatic and environmental changes during the last two centuries on better than a decadal basis. Sulfate and nitrate concentrations started to increase around 1900 and 1960, respectively, due to anthropogenic influx transported from the industrialized regions in North America. Sea-salt concentrations began to increase around the mid-19th century and were elevated throughout the 20th century. This trend of sea-salt concentrations is similar to that of melt percentage, which is a measure of summer temperature. Warming after the Little Ice Age would have reduced the sea-ice extent and led to the elevated sea-salt concentrations on Penny Ice Cap.

The U.S. Science Plan for Deep Ice Coring in West Antarctica calls for two ice cores to be collected. the first of these cores, from Siple Dome, was completed during the 1997/98 field season. the second core is to be collected from a site near the divide that separates ice flowing to the Ross Sea and to the Amundsen Sea.Using high-resolution, grid-based aerogeophysical surveys of the Ross/Amundsen ice-divide region, we identify seven candidate sites and assess their suitability for deep coring. We apply ice-flow and temperature calculations to predict time-scales and annual-layer resolution, and to assess the potential for basal melting for several selected sites. We conclude that basal melting is likely for sites with very thick ice, as was observed at the Byrd core site. Nevertheless, these sites are most attractive for coring since they promise recovery of a long climate record with comparatively high time resolution during the last glacial period.

Two medium-depth ice cores were retrieved from Berkner Island by a joint project between the Alfred-Wegener-Institut and the British Antarctic Survey in the 1994/95 field season. A 151m deep core from the northern dome (Reinwarthhöhe) of Berkner Island spans 700 years, while a 181 m deep core from the southern dome (Thyssenhöhe) spans approximately 1200 years. Both cores display clear seasonal cycles in electrical conductivity measurements, allowing dating by annual-layer counting and the calculation of accumulation profiles. Stable-isotope measurements (both δ18O and δD), together with the accumulation data, allow us to estimate changes in climate for most of the past millennium: the data show multi-decadal variability around a generally stable long-termmean. In addition, a full suite of major chemistry measurements is available to define the history of aerosol deposition at these sites: again, there is little evidence that the chemistry of the sites has changed over the past six centuries. Finally, we suggest that the southern dome, with an ice thickness of 950 m, is an ideal site from which to gain a climate history of the late stages of the last glacial and the deglaciation for comparison with the records from the deep Antarctic ice cores, and with other intermediate-depth cores such as Taylor Dome and Siple Dome.

Vapor transport in snow and firn plays a key role in post-depositional changes of many reactive species, in mass and energy balance of large snow-covered areas, and in metamorphic changes of the crystal structure of the snow itself. While conventional estimates of vapor transport and sublimation rate are based on diffusion of water vapor from the snow surface to the atmosphere, ventilation (airflow through interstitial pore spaces) in snow and firn can affect the top several meters of firn, yet the effects of ventilation on sublimation rates in firn have not been previously investigated. In this paper we present the first calculations of sublimation rates due to advective/ diffusive airflow and vapor transport in the snow and firn. Airflow velocities, vapor transport and sublimation rates are calculated using a two-dimensional finite-element model. the airflow patterns induce regions of sublimation and regions of condensation within the snow and firn. Because small surface sastrugi move in time as the snow surface profile changes due to wind redeposition, the subsurface flow and regions of condensation and sublimation will also likely change on a time-scale of days to weeks. However, if the roughness features are such that they move very little over time, it is likely that the regions of condensation and sublimation will have a noticeable effect on the microstructure of the firn. While the highest mass-transfer rates occur near the surface, the depth of vapor transport and phase change depends on firn properties, wavelength of the snow surface roughness, and temperature.

We use the upper 81 mof the record of stable isotopes of water from a 122m long ice core from Lomonosovfonna, central Spitsbergen, Svalbard, to construct an ice-core chronology and the annual accumulation rates over the icefield. the isotope cycles are counted in the ice-core record using a model that neglects short-wavelength and low-amplitude cycles. We find approximately the same number of δ18O cycles as years between known reference horizons, and assume these cycles represent annual cycles. Testing the validity of this assumption using cycles in δD shows that both records give similar numbers of cycles. Using the δ18O chronology, and decompressing the accumulation records using the Nye flow model, we calculate the annual accumulation for the ice-core site back to AD 1715. We find that the average accumulation rate from 1715 to 1950 was on average 0.30 mw.e. Accumulation rates increased about 25% during the later part of the 20th century to an average of 0.41 mw.e. for the period 1950–97. the accumulation rates show highly significant 2.1 and 21 year periodicities, which gives credibility to our time-scale.

Ice cores provide a valuable archive of climate history. for a complete understanding of this archive, it is important to understand air–snow exchange processes through the snow and firn in order to fully decode the ice-core record. Transport processes through the snow and firn are dependent upon their physical properties. In this paper, bidirectional permeabilities from selected sections of a 13 mcore from Summit, Greenland, are presented. Differences between lateral and vertical permeabilities are evident throughout the core in permeameter data and in microstructure statistics. Both lateral and vertical permeabilities are consistent with overall patterns of previous polar permeability data with depth. the differences between lateral and vertical permeability measurements for some samples can be attributed to equivalent sphere radius. Further studies examining mean free-path length may be helpful in chemical modeling and in deriving an equation relating permeability to microstructure.

We used surface climate fields from high-resolution (~0.5660.56˚) European Centre for Medium-RangeWeather Forecasts (ECMWF) operational analyses (1992–98), together with meteorological and glaciological models of snow accumulation and surface meltwater runoff/retention, to produce novel maps of Greenland ice sheet (GIS) net accumulation, net runoff and surface mass balance (SMB). We compared our runoff maps with similar-scaled runoff (melt minus refreezing) maps based on passive-microwave satellite data. Our gross spatial/temporal patterns of runoff compared well with those from the satellite data, although amounts of modelled runoff are likely too low. Mean accumulation was 0.287 (0.307)ma–1, and mean runoff was 0.128 (0.151)ma–1, averaged across the W. Abdalati (T. L. Mote) GIS mask. Corresponding mean SMB was 0.159 (0.156)ma–1, with considerable interannual variability (standard deviation ~0.11ma–1) primarily due to variations in runoff. Considering best estimates of current iceberg calving, overall the GIS is probably currently losing mass. Our study shows great promise for meaningfully modelling SMB based on forthcoming ``second-generation’’ ECMWF re-analysis (ERA-40) data, and comparing the results with ongoing laser/radarmeasurements of surface elevation. This should help elucidate to what extent surface elevation changes are caused by short-term SMB variations or other factors (e.g. ice dynamics).

Extensive aircraft-based radar ice-thickness measurements over the interior and outlet-glacier regions of the Greenland ice sheet have been obtained by the University of Kansas since 1993, with the latest airborne surveys conducted in May 2001. the radar has evolved during this period to a highly versatile system capable of characterizing ice thickness over a wide variety of ice-sheet conditions. Before 1997, the digital system was limited, only capable of storing incoherent data or coherent data with a very large number of presumed signals at a low pulse-repetition frequency. In 1998, the radar was upgraded with modern components allowing coherent data to be stored with a small number of presumed returns for 1024 range cells at a high pulse-repetition frequency. the new data on ice thickness of Greenland outlet glaciers are archived and made available to the scientific community in the form of radar echograms and derived ice thickness at http://tornado.rsl.ukans.edu/Greenlanddata.htm. the U.S. National Snow and Ice Data Center (NSIDC) also provides a link to these data, and NSIDC will eventually serve as the permanent archive of these data. Improvements in radar sensitivity in outlet-glacier regions have been achieved by collecting coherent radar data and applying various signal-processing techniques. Deep outlet-glacier channels that were previously unresolved with incoherent data can now be mapped using a coherent signal, signal conditioning and synthetic aperture radar (SAR) processing.

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. the results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. the sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. the results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

A 45.97 m long ice core was recovered in the accumulation area of Glaciar Tyndall (50˚59’05’’ S, 73˚31’12’’W; 1756ma.s.l.), Campo de Hielo Patagόnico Sur (southern Patagonia icefield), during December 1999. the firn core was subjected to visual stratigraphic observation and bulk density measurements in the field, and later to analyses of water isotopes (δ18O, δD), major dissolved ions and snow algal biomass. the drillhole remained dry down to about 43 m depth, where a water-soaked layer appeared. Seasonal cycles were found for δ18O, δD and the D-excess, although the amplitudes of the cycles decreased with depth. Major dissolved ions (Na+, K+, Mg2+, Ca2+, Cl–, SO42–) and algal biomass exhibit rapid decreases in the upper 3 m, probably due to meltwater elution. Annual increments defined by the δ18O and D-excess peaks suggest that the minimum net accumulation rates at this location were 17.8ma–1 in 1997/98–1998/99 and 411.0 ma–1 in 1998/99–1999/2000. These are much higher values than those previously obtained from past ice-core studies in Patagonia, but are of the same order of magnitude as those predicted from various observations in ablation areas of Patagonian glaciers.

Since 1979/80, glaciological studies have been carried out at Ekströmisen, Antarctica, including accumulation-stake measurements, snow-pit and shallow-firn-core studies. Snow stratigraphy, chemical properties and stable-isotope ratios (δ18O) were investigated. This study focuses on three cores taken between 1982 and 1998. the 1998 core was dated using dielectric profiling, δ18O profiles and stake measurements. Accumulation rates show high interannual and spatial variability due to the extreme wind influence. No significant trend was found for the last 50 years; during the first half of the 20th century, accumulation decreased. the high spatial and interannual variability, however, means that trends must be interpreted with care. In spite of the highly irregular accumulation distribution, stable-isotope ratios show little spatial variability. the mean annual δ18O values of cores B04 and FB0198 agree fairly well for the time period 1955–82 covered by both cores. δ18O values have increased during most of the 20th century; since the late 1980s a decrease is observed. This change is not related to air temperature, since mean annual air temperatures at Neumayer show no significant trend over the last two decades.

At the German wintering base Neumayer, an intensive glacio-meteorological programme has been carried out during the last two decades. A complete meteorological dataset and data from surface snow samples, snow pits, firn cores and weekly accumulation measurements from a stake array are available. We first investigated the attenuation of the seasonal δ18O signal due to water-vapour diffusion in the snowpack. A comparison of surface snow samples and firn cores of different ages shows that only one-third of the seasonal δ18O signal of the surface snow samples remains in the cores after the first year. No further significant change in the amplitude of the seasonal δ18O signal is found later. Changes in the seasonal distribution of accumulation can lead to a bias in ice-core properties. This is studied on a short time-scale, using high-time-resolution data of accumulation, stable-isotope ratios and air temperature. Mean annual δ18O values from firn cores are not well correlated to annual mean air temperatures. However, the correlation is improved considerably by calculating an annual mean air temperature using monthly mean temperatures weighted by monthly accumulation. At Neumayer, it is mainly the cyclonic activity in late winter/early spring that determines whether and how the core data are biased.

Snow-accumulation rates and rates of ice-thickness change (mass balance) are studied at several sites on Siple Dome,West Antarctica. Accumulation rates are derived from analyses of gross beta radioactivity in shallow firn cores located along a 60 km transect spanning both flanks and the crest of the dome. There is a north–south gradient in snow-accumulation rate across the dome that is consistent with earlier radar mapping of internal stratigraphy. Orographic processes probably control this distribution. Mass balance is inferred from the difference between global positioning system (GPS)-derived vertical velocities and snow-accumulation rates for sites close to the firn-core locations. Results indicate that there is virtually no net thickness change at four of the five sites. the exception is at the northernmost site where a small amount of thinning is detected, that appears to be inconsistent with other studies. A possible cause of this anomalous thinning is recent retreat of the grounding line of Ice Stream D.

Records of recent oxygen isotope ratios (δ18O) and accumulation rates are presented for the region of Wilhelm II Land, East Antarctica, between 78˚ and 93˚E and from the coast to 2100m elevation. These records were derived from analysis of 21 shallow firn cores collected during the 1997/98 and 1998/99 Australian National Antarctic Research Expeditions summer operations. the accumulation rates were determined using comparisons between detailed analyses of density, δ18O, hydrogen peroxide (H2O2) levels and electrical conductivity. the δ18O distribution follows an approximately linear relationship with snow surface elevation, with values from –22‰ near the coast to –32‰ towards 2000m elevation. Accumulation-rate distribution does not display this simple relationship with topography. South of the West Ice Shelf the contours run parallel to lines of latitude (oblique to the coast and topography), with 400 kg m–2 a–1 towards the coast and 2000m elevation, and a lower zone of 300 kg m–2 a–1 along an axis of 68.4˚ S. This pattern of accumulation is also evident along the Mirny–Vostok traverse route. Southwest of the West Ice Shelf the rate of accumulation drops gradually from 300 to 200 kg m–2 a–1 towards Lambert Glacier basin. Surface-snow redistribution and variations in accumulation rate cause variability in the clarity of core records, but several sites show sufficient stratigraphic preservation to suggest potential for extraction of extended palaeoenvironmental records through further drilling.

It is generally recognized at present that ice-core deuterium excess measurements are potentially useful for reconstructions of vapor source-region temperature and humidity history, and that such measurements provide a method for correcting isotopic paleothermometers for effects of source temperature variations. Here we use a zonally averaged vapor-transport and isotopic-distillation model to show that deuterium excess of precipitation on the ice sheets is sensitive to a wide variety of source-region climate changes in addition to those changes of temperature and humidity that affect the composition of evaporate. Moreover, it is demonstrated that this wide variety of source-region changes all cause anticorrelated changes in deuterium excess with δD and δ18O over the ice sheets, suggesting that deuteriu 14m excess is a generally more useful tool for correcting isotopic thermometers than is currently recognized.