Deep drilling of glaciers: Russian projects in the Arctic (1975-1995)

The data collection "Deep Drilling of Glaciers: Soviet-Russian projects in Arctic, 1975-1995" was collected by the following basic considerations: - compilation of deep (>100 m) drilling projects on Arctic glaciers, using data of (a) publications; (b) archives of IGRAN; (c) personal communication of project participants; - documentation of parameters, references. Accuracy of data and techniques applied to determine different parameters are not evaluated. The accuracy of some geochemical parameters (up to 1984 and heavy metalls) is uncertain. Most reconstructions of ice core age and of annual layer thickness are discussed; - digitizing of published diagrams (in case, when original numerical data were lost) and subsequent data conversion to equal range series and adjustment to the common units. Therefore, the equal-range series were calculated from original data or converted from digitized chart values as indicated in the metadata. For the methodological purpose, the equal-range series obtained from original and reconstructed data were compared repeatedly; the systematic difference was less then 5-7%. Special attention should be given to the fact, that the data for individual ice core parameters varies, because some parameters were originally measured or registered. Parameters were converted in equal-range series using 2 m steps; - two or more parameter values were determined, then the mean-weighted (i.e. accounting the sample length) value is assigned to the entire interval; - one parameter value was determined, measured or registered independently from the parameter values in depth intervals which over- and underlie it, then the value is assigned to the entire interval; - one parameter value was determined, measured or registered for two adjoining depth intervals, then the specific value is assigned to the depth interval, which represents >75% of sample length ; if each of adjoining depth intervals represents <75% of sample length, then the correspondent parameter value is assigned to both intervals of depth. This collection of ice core data (version 2000) was made available through the EU funded QUEEN project by S.M. Arkhipov, Moscow.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Stable isotope ratios of foraminifera from last glacial-interglacial sediments of the southern high latitudes

A set of numerical equations is developed to estimate past sea surface temperatures (SST) from fossil Antarctic diatoms. These equations take into account both the biogeographic distribution and experimentally derived silica dissolution. The data represent a revision and expansion of a floral data base used previously and includes samples resulting from progressive opal dissolution experiments. Factor analysis of 166 samples (124 Holocene core top and 42 artificial samples) resolved four factors. Three of these factors depend on the water mass distribution (one Subantarctic and two Antarctic assemblages); factor 4 corresponds to a 'dissolution assemblage'. Inclusion of this factor in the data analysis minimizes the effect of opal dissolution on the assemblages and gives accurate estimates of SST over a wide range of biosiliceous dissolution. A transfer function (DTF 166/34/4) is derived from the distribution of these factors versus summer SST. Its standard error is +/- 1°C in the -1 to +10 °C summer temperature range. This transfer function is used to estimate SST changes in two southern ocean cores (43°S and 55°S) which cover the last climatic cycle. The time scale is derived from the changes in foraminiferal oxygen and carbon isotopic ratios. The reconstructed SST records present strong analogies with the air temperature record over Antarctica at the Vostok site, derived from changes in the isotopic ratio of the ice. This similarity may be used to compare the oceanic isotope stratigraphy and the Vostok time scale derived from ice flow model. The oceanic time scale, if taken at face value, would indicate that large changes in ice accumulation rates occurred between warm and cold periods.

The influence of carbon surface oxygen groups on Dubinin–Astakhov equation parameters calculated from CO2 adsorption isotherm

We present the results of a systematic study of the influence of carbon surface oxidation on Dubinin–Astakhov isotherm parameters obtained from the fitting of CO2 adsorption data. Using GCMC simulations of adsorption on realistic VPC models differing in porosity and containing the most frequently occurring carbon surface functionalities (carboxyls, hydroxyls and carbonyls) and their mixtures, it is concluded that the maximum adsorption calculated from the DA model is not strongly affected by the presence of oxygen groups. Unfortunately, the same cannot be said of the remaining two parameters of this model i.e. the heterogeneity parameter (n) and the characteristic energy of adsorption (E0). Since from the latter the pore diameters of carbons are usually calculated, by inverse-type relationships, it is concluded that they are questionable for carbons containing surface oxides, especially carboxyls.

Sea surface temperature and age analysis from sediment core KT94-15_PC-9, XP98-PC-4 and XP98-PC-2

The Japan Sea experienced bottom water anoxia at the last glacial maximum (LGM) since it is surrounded by four shallow straits, the sill depths of which are close to, or shallower than, the drop in sea level (~120 m) that occurred then. A distinctive negative d18O excursion of planktonic foraminifera also took place during the LGM. This excursion has been interpreted from foraminiferal data as recording a drop in the paleosalinity of surface waters on the assumption of a constant low sea surface temperatures between 34 and 11 ka. We present here a profile of alkenone-based sea surface temperatures (alkenone-SSTs) over the past 36 kyr. Our results suggest that SSTs during the LGM were much higher than those previously assumed. After considering the factors that might affect estimation of alkenone-SSTs and comparisons of core-top alkenone-SSTs values with values for modern seawater we conclude that the higher alkenone-SSTs during the LGM are reliable and reasonable. These warm SSTs were probably caused by radiative equilibrium associated with the development of stable water stratification in the Japan Sea during the LGM.

Stratigraphic Variation Within Polar Firn Caused by Differential Accumulation and Ice Flow: Interpretation of a 400 Mhz Short-Pulse Radar Profile from West Antarctica

We investigate causes of the stratigraphic variation revealed in a 177 km, 400 MHz short-pulse radar profile of firn from West Antarctica. The profile covers 56 m depth, and its direction was close to those of the ice flow and mean wind. The average, near-surface accumulation rates calculated from the time delays of one radar horizon consistently show minima on leeward slopes and maxima on windward slopes, confirming an earlier study based on stake observations. The stratigraphic variation includes up to 30 m depth variation in individual horizons over tens of km, fold limbs that become progressively steeper with depth, and fold-hinge loci that change direction or propagate down-ice with depth over distances far less than predicted by the ice speeds. We use an accumulation rate model to show how local rate anomalies and the effect of ice speed upon a periodic variation in accumulation rate cause these phenomena, and we reproduce two key features seen in the stratigraphic variations. We conclude that the model provides an explanation of changes in spatial stratigraphy and local measures of accumulation history given the constraints of surface topography, ice and wind velocities, and a general accumulation rate for an area.

Chronology for the EDML 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.

Density, d18O and accumulation rates from snow pits on the Filchner-Ronne Ice Shelf

Accumulation rates in the eastern part of Ronne Ice Shelf were determined by isotopic stratigraphy (18O). The samples were taken from snow-pits dug during the Filchner I and II operations in 1984 and 1986. In general, the accumulation rate decreases towards the south; the greatest decrease, from 21.3 to 13.3 g/cm**2/a, was observed between Filchner Station and measuring point 341, sited 270 km up-stream of the ice edge. The d18O values of the near-surface layers vary between -25 and -29 per mil. The 18O content in the more southerly part is progressively depleted in the direction of Möllereisstrom, paralleling a decrease in the accumulation rate. Near the ice edge the 18O content decreases to the west. A 100 m ice core drilled in 1984 at point 340, 22 km from the ice edge, probably goes back to A.D. 1460; it has been dated by isotopic stratigraphy. The accumulation rate up-stream of the drilling site was deduced from the sequence of annual layers, using a simple ice-flow model. The accumulation shows strong variations over the last 200 years, which may be caused in part by local variations in the accumulation on Ronne Ice shelf.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

Chemical concentrations and fluxes from EPICA ice cores EDML and EDC

Continuous sea salt and mineral dust aerosol records have been studied on the two EPICA (European Project for Ice Coring in Antarctica) deep ice cores. The joint use of these records from opposite sides of the East Antarctic plateau allows for an estimate of changes in dust transport and emission intensity as well as for the identification of regional differences in the sea salt aerosol source. The mineral dust flux records at both sites show a strong coherency over the last 150 kyr related to dust emission changes in the glacial Patagonian dust source with three times higher dust fluxes in the Atlantic compared to the Indian Ocean sector of the Southern Ocean (SO). Using a simple conceptual transport model this indicates that transport can explain only 40% of the atmospheric dust concentration changes in Antarctica, while factor 5-10 changes occurred. Accordingly, the main cause for the strong glacial dust flux changes in Antarctica must lie in environmental changes in Patagonia. Dust emissions, hence environmental conditions in Patagonia, were very similar during the last two glacials and interglacials, respectively, despite 2-4 °C warmer temperatures recorded in Antarctica during the penultimate interglacial than today. 2-3 times higher sea salt fluxes found in both ice cores in the glacial compared to the Holocene are difficult to reconcile with a largely unchanged transport intensity and the distant open ocean source. The substantial glacial enhancements in sea salt aerosol fluxes can be readily explained assuming sea ice formation as the main sea salt aerosol source with a significantly larger expansion of (summer) sea ice in the Weddell Sea than in the Indian Ocean sector. During the penultimate interglacial, our sea salt records point to a 50% reduction of winter sea ice coverage compared to the Holocene both in the Indian and Atlantic Ocean sector of the SO. However, from 20 to 80 ka before present sea salt fluxes show only very subdued millennial changes despite pronounced temperature fluctuations, likely due to the large distance of the sea ice salt source to our drill sites.