Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Age and chemical composition of apatites from five DSDP and ODP Sites in the Pacific, Atlantic and Indian Oceans

Since studies on deep-sea cores were carried out in the early 1990s it has been known that ambient temperature may have a marked affect on apatite fission track annealing. Due to sluggish annealing kinetics, this effect cannot be quantified by laboratory annealing experiments. The unknown amount of low-temperature annealing remains one of the main uncertainties for extracting thermal histories from fission track data, particularly for samples which experienced slow cooling in shallow crustal levels. To further elucidate these uncertainties, we studied volcanogenic sediments from five deep-sea drill cores, that were exposed to maximum temperatures between ~10° and 70°C over geological time scales of ~15-120 Ma. Mean track lengths (MTL) and etch pit diameters (Dpar) of all samples were measured, and the chemical composition of each grain analyzed for age and track length measurements was determined by electron microprobe analysis. Thermal histories of the sampled sites were independently reconstructed, based on vitrinite reflectance measurements and/or 1D numerical modelling. These reconstructions were used to test the most widely used annealing models for their ability to predict low-temperature annealing. Our results show that long-term exposure to temperatures below the temperature range of the nominal apatite fission track partial annealing zone results in track shortening ranging between 4 and 11%. Both chlorine content and Dpar values explain the downhole annealing patterns equally well. Low chlorine apatite from one drill core revealed a systematic relation between Si-content and Dpar value. The question whether Si-substitution in apatite has direct and systematic effects on annealing properties however, cannot be addressed by our data. For samples, which remained at temperatures <30°C, and which are low in chlorine, the Laslett et al. [Laslett G., Green P., Duddy I. and Gleadow A. (1987) Thermal annealing of fission tracks in apatite. Chem. Geol. 65, 1-13] annealing model predicts MTL up to 0.6 µm longer than those actually measured, whereas for apatites with intermediate to high chlorine content, which experienced temperatures >30°C, the predictions of the Laslett et al. (1987) model agree with the measured MTL data within error levels. With few exceptions, predictions by the Ketcham et al. [Ketcham R., Donelick R. and Carlson W. (1999) Variability of apatite fission-track annealing kinetics. III: Extrapolation to geological time scales. Am. Mineral. 84/9, 1235-1255] annealing model are consistent with the measured data for samples which remained at temperatures below ~30°C. For samples which experienced maximum temperatures between ~30 and 70°C, and which are rich in chlorine, the Ketcham et al. (1999) model overestimates track annealing.

Fly-through movie for the reconstructed folds in the onset region of the Petermann Glacier, North Greenland

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.

(Table 1) Moisture content of buried snowpack deposit samples from Pearse Valley, McMurdo Dry Valleys

Buried snowpack deposits are found within the McMurdo Dry Valleys of Antarctica, which offers the opportunity to study these layered structures of sand and ice within a polar desert environment. Four discrete buried snowpacks are studied within Pearse Valley, Antarctica, through in situ observations, sample analyses, O-H isotope measurements and numerical modelling of snowpack stability and evolution. The buried snowpack deposits evolve throughout the year and undergo deposition, melt, refreeze, and sublimation. We demonstrate how the deposition and subsequent burial of snow can preserve the snowpacks in the Dry Valleys. The modelled lifetimes of the buried snowpacks are dependent upon subsurface stratigraphy but are typically less than one year if the lag thickness is less than c. 7 cm and snow thickness is less than c. 10 cm, indicating that some of the Antarctic buried snowpacks form annually. Buried snowpacks in the Antarctic polar desert may serve as analogues for similar deposits on Mars and may be applicable to observations of the north polar erg, buried ice at the Mars Phoenix landing site, and observations of buried ice throughout the martian Arctic. Numerical modelling suggests that seasonal snows and subsequent burial are not required to preserve the snow and ice on Mars.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.

(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ~34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved under¬standing of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Age and alkenone-derived Holocene sea-surface temperature records of sediment core M40-4-SL78

This Special Issue of The Holocene contains 16 research papers based on a symposium at the 11th International Meeting of the European Union of Geosciences held in Strasbourg in April 2001. The aim of the symposium was a state-of-the-art assessment of empirical studies of postglacial marine and terrestrial climatic archives and their integration with numerical climate models. This editorial places the individual papers in the broader context of natural climate variability and anthropogenic impacts on the global climate system, regional differences in climate between maritime and continental areas, and the need for an improved theoretical basis for understanding the underlying causes of environmental change. The focus of the Special Issue is the dynamic and relatively well-understood climate of the North Atlantic and the European realm, where, in relation to the steepest offshore temperature gradient on Earth, observational data are abundant and many recent advances have been made in climate reconstruction from proxy archives. The editorial also contains a summary and overview of the papers included in the four main sections of the Special Issue, which emphasize: (1) numerical modelling experiments; (2) models of glacier buildup and equilibrium-line altitude; (3) marine and terrestrial proxy records of climatic change; and (4) multiproxy palaeoenvironmental reconstruction of a Portuguese lagoonal system.