Dry and wet air-extraction systems and precise analysis systems of the CO2 and CH4 concentrations for a polar ice core were developed to reconstruct their ancient levels. A dry-extraction system was capable of crushing an ice sample of 1000 g into fine powder within 2 min, and its air-extraction efficiency was found to be 98%. The CO2 and CH4 concentrations of extracted air were determined using gas chromatography with a flame-ionized detector. The overall precision of our measurements, including air extraction, was estimated to be better than ± 1 ppmv for CO2 and + 10 ppbv for CH4. Preliminary analysis of the ice core drilled at Mizuho Station, Antarctica, showed that the CO2 and CH4 concentrations at 3340–3700 year BP were about 280 ppmv and 700ppbv, respectively. The Yamato core drilled at the terminus of the glacial flow near the Yamato Mountains, Antarctica, yielded concentrations of 230–240 ppmv for CO2 and 520–550 ppbv for CH4, suggesting that the core had formed during the glacial period.

A mass-balance approach using hourly discharge and electrical conductivity values measured over a 10 d period during the ablation season was used to separate englacial and subglacial components of the total meltwater discharge from a small alpine glacier in the Sierra Nevada, California, U.S.A. Symmetrical diurnal hydrographs indicate that little delay occurred as water was tranferred through the drainage system. Electrical conductivity (which varied inversely with proglacial discharge) increased abruptly at each daily conductivity maximum, and cross-correlation analysis indicated that subglacial discharge peaked on the rising limb of the englacial hydrograph (about 2 h prior to the englacial peak). This suggests that a translatory flow process operates in which increased water pressure in the englacial system on the rising limb of the diurnal-discharge cycle forced subglacial water from beneath the glacier in advance of short residence-time meltwater. Net radiation dominated the energy balance at the glacier surface, explaining 86% of the variance in proglacial discharge, which was dominated by the englacial flow component.

Ice cores were obtained in January 1990 from the land-fast ice in McMurdo Sound for a study of variations in texture, fabric, sub-structure, composition and development. Two primary ice types were observed, congelation and platelet, with a minor amount of frazil ice. Congelation ice growth precedes platelet-ice accretion. Congelation-ice fabrics show frequent moderate to strong alignments, a phenomenon believed to be due to water-current control of selective ice-crystal growth. Platelet ice originates at the base of the congelation ice, initially as a porous latticework of tabular ice crystals which subsequently consolidate by congelation of the interstitial water. Interstitial congelation-ice fabrics generally have little or no alignment, indicating the reduced effect of currents within the platelet latticework prior to solidification. Platelet-crystal textures range from small, wavy-edged forms to large, blade-like forms. Platelet-crystal fabrics indicate that, in addition to being randomly oriented, the platelet latticeworks commonly include many crystals with their flat (0001) faces oriented both parallel and normal to the base of the overlying ice. Plate-width data suggest that the interstitial congelation ice-growth rates remain similar to those of the overlying congelation ice. This effective increase in growth rates probably happens because the latticework of accumulating platelets ahead of the freezing interface ensures that the water within the platelet layer is at the freezing point and less heat has to be removed from platelet-rich water than from platelet-free water for a given thickness of congelation ice to grow. The negative oceanic heat flux associated with platelet-ice formation in McMurdo Sound explains why McMurdo Sound fast ice is thicker than Ross Sea pack ice, and also why it reaches a greater thickness than Arctic fast ice grown in a similar polar marine climate. Plate widths in the McMurdo Sound congelation ice suggest, however, that it grows no faster than Arctic congelation ice.

In a pilot experiment, airborne particles were shown to exist in snow. In newly deposited snow they could be traced down to 17 cm below the surface. With our particle sensor, the snow was ventilated on the level of expected natural ventilation velocities. We show with a simple deposition model that air/snow exchange of airborne particles must be considered in the interpretation of impurities in snow and glacier ice. However, the relative magnitude of ventilation velocity compared to dry-deposition velocities at the surface is of crucial importance for determining total deposition. In particular, in ice sheets with high ventilation velocities, seasonal variations in atmospheric concentrations can be dampened and age distributions of deposited particles need to be considered similarly to the occlusion of gases.

A new method for high-resolution (3 mm) dielectric profiling of ice cores is described. The method measures the capacitance and conductance of the ice at a single frequency (chosen to be 50 kHz in this case). The technique requires only a few minutes per metre of core, and no sophisticated data analysis. The system is designed to operate on ice cores that have been cut longitudinally, providing one flat and one curved surface. The 50 kHz conductance is dependent on the acid and neutral-salt concentrations in the ice. The new method was tested successfully on the GRIP core from Summit, Greenland, in 1991. It is useful for detecting seasonal variations in impurities in both acidic Holocene ice and in normally alkaline Wisconsin ice.

A shallow firn core, recovered from the southern Patagonia ice cap (elevation 2680 m a.s.I.), has been analyzed for H stable isotope composition and for major soluble chemical compounds. The temperature measurement at 13.17 m depth (bottom of the core) shows that the ice cap is temperate. The chemical profiles indicate that some soluble impurities have been partly washed out by percolation, but seasonal deuterium content variations are relatively well preserved down to the bottom of the core, which allows the net accumulation rate of the site to be calculated (1.2 m water equivalent). The mean concentrations measured in the most recent year recorded (Cl−:4.85, nssSO42−: 0.65, NO3−:0.32, Na+:4.30, K+:0.64 and NH4+: 0.77, in μeq. 1−1) serve as a reference for the background chemical composition of precipitation at mid-southern latitudes under South Pacific meteorological conditions.

The presence of crevasses on the surface of ice masses indicates that a fracture criterion has been met. Understanding how crevasses form will provide information about the stress and strain-rate fields in the ice. This study derives a relationship between measurements of strain rate and observations of crevassing on the surface of ice masses. A literature search yielded 17 polar and alpine locations where strain rates had been measured and crevassing recorded. By plotting strain rates (converted to stresses using a creep law) using axes representing the surface-parallel principal stresses, failure envelopes were derived by enclosing measurements where surface crevassing was absent. The derived failure envelopes were found to conform well to theoretical ones predicted by the Coulomb and the maximum octahedral shear stress (von Mises) theories of failure. The derived failure envelopes were scaled by the tensile strength, which was found to vary from 90 to 320 kPa. There was no systematic variation of tensile strength with either temperature at 10 m depth or the method used to locate the crevasses. The observed variation in tensile strength could result from variations in ice properties (e.g. crystal size, impurity content or density) or could be related to uncertainty in the constitutive relation. Creep flow and fracture share a very similar temperature dependence, suggesting similar crystal-scale processes are responsible for both. The observed relationship will provide a supplementary tool with which to verify and test models of ice dynamics against remotely sensed imagery. The study also indicates that a temperature rise of a few degrees throughout the ice column will not result directly in any increase in calving rates from the large Antarctic ice shelves such as the Filchner–Ronne or Ross Ice Shelves.

Monthly ice ablation was measured at the margin of the Greenland ice sheet for June, July and August over 7 years (1980–86). The total winter ablation (September-May) has also been measured, and a simple ablation-temperature model used to assign ablation values to individual months. Under the present climate, the most ablation occurred in June-August (on average 81% of annual ablation), moderate ablation took place in May and September (17%) and very little ablation occurred in October-April (2%). The effect of climate change on ice ablation is simulated using the ablation model to recalculate ablation for higher temperatures. Summer ice ablation increases with temperature in the model, but there is proportionally greater increase for May and September, whereas the period from October to April is presently so cold that even a temperature rise of +5 °C will hardly increase ablation. The difference in annual ice ablation caused by future climate change will therefore depend upon the seasonal distribution of the temperature change. Changes in precipitation and accumulation will further modify the seasonal variation of ablation.

A simple model is developed to simulate the isotopic fractionation which accompanies Weertman regelation at the bed of temperate-based glaciers. The fractionation equations of Jouzel and Souchez (1982) are applied to multiple refreezing events over measured glacier-bed profiles, and mass balance is maintained as the basal ice and meltwater produced at one bedrock hummock enter the next. Simulation results indicate that undeformed regelation ice layers are on the order of millimetres to centimetres thick, often being completely melted at the stoss face of certain hummocks and exceptionally reaching a thickness in excess of 10 cm. Neither the internal morphology nor the isotopic composition of these layers is constant, but both vary down-glacier in accordance with bedrock configuration. A glacier-wide fractionation process is identified whereby heavy isotopes are preferentially removed from the basal meltwater film and incorporated into the basal ice. This process might go some way to explaining the anomalously “light” isotopic composition measured in base-flow waters leaving some glaciers. Vertical isotope profiles through undeformed basal ice layers are reconstructed and show that significant isotopic excursions can occur at a scale of millimeteres, while the range of isotopic compositions within such multi-layered regelation ice is greater than that which would occur in ice produced by a single refreezing event. In circumstances where the regelation system is disrupted by removal of film waters into a network of linked cavities, it is found that the remaining basal ice may be significantly enriched in heavy isotopes relative to the composition of the initial mass inputs to the system. Heavy isotope enrichment of this magnitude and consideration of the thickness of the basal ice layers concerned may explain the absence of recorded basal ice samples heavy enough to have been formed in equilibrium with subglacial precipitates sampled at one of the sites (Glacier de Tsanfleuron) and reported in an earlier paper.

The mechanical properties of ice are strongly affected by crystal texture and c-axis alignment. In this paper, we develop a general quantitative method for analysis of uniaxial crystal-orientation data. These data are represented as unit vectors from the origin with end points on the surface of a unit sphere. An orthogonal least-squares error measure is used to develop equations that define the closest plane and line through the data. The resulting eigenvalue problem is identical to that obtained by other investigators using different methods. However, here we identify an implicit assumption in the method, and observe that the error measure represents physical distance and quantifies the goodness of fit to the data of idealized structures. Also, a method is developed to transform the data and the results for viewing on Schmidt nets drawn in the best plane and the predominant basal plane of a sample, in addition to the standard xy-plane. Applications of the analysis to sea-ice samples include both numerical and Schmidt-net presentations of results.

The debris-rich ice from the bottom 6 m of the 82 m deep CAROLINE (Coastal Antarctic Record of Last Interglacial Natural Environment) ice core reaching bedrock, and from five 2 m long surface cores at Moraine Prudhomme in Terre Adélie (Antarctica) is described and compared to debris-laden ice from the core-drilling site DIO. Isotopic, total-gas content, CO2 concentration and SEM investigations of embedded particles, together with ice textures and fabrics, rule out “pressure-melting” regelation around bed obstacles or “freezing-on” as possible mechanisms for the debris entrainment at the ice-bedrock interface. It is suggested that the debris entrapment by purely mechanical means (e.g. shearing) is an efficient process in forming basal ice layers (BIL) at sub-freezing temperatures. This process might be dominant at the margin of the Antarctic ice sheet where no ice shelf exists and where a ramp terminus or a buttressing coastal relief induces compressive flow.

Individual variables found to be significant from a correlation analysis are analyzed as a function of probability of avalanching for data from Kootenay Pass, British Columbia. The analysis is compared with a similar study for data from Alta, Utah, U.S.A. The results show that the variable significance is very similar for the two areas. Primary variables include: snowfall rate, weight of new snow, water equivalent of new precipitation, total storm snow and new snow depth. Secondary variables include wind speed and direction, and new-snow density.

The relationship of glacier response due to mass-balance changes is of fundamental importance when climate variations are to be understood. In this paper, two aspects of the problem are analyzed from field data: (1) advance/retreat of the glacier terminus due to changes in mass balance, and (2) cross-correlation of mass-balance data from two glaciers in the same climate zone. The results show: (1) the terminus can respond quickly in accordance with expected minimum time-scale, and (2) two glaciers in the same general climate zone may have very different yearly mass balance and advance/retreat behaviour. This latter result indicates the importance of local climate variations.

In this paper we develop a theoretical model describing water motion in a coupled borehole-subglacial flow system. The theory applies to basal drainage systems having multiple and extensive interconnected flow paths. Within this domain it encompasses a broad range of flow regimes, from laminar Darcian flow in a thick permeable unit to turbulent sheet flow in a very thin layer. Important terms in the model are highlighted by recasting the problem in dimensionless form. The non-dimensional formulation indicates that there are four free parameters in the coupled system. These parameters characterize skin friction in the borehole, and diffusion, transmissivity and turbulent transport in the subglacial flow layer. Dimensionless results show that, under most circumstances, the effects of skin friction in the borehole are negligible. Diffusion, transmissivity and especially turbulent transport in the basal layer are found to influence subglacial water flow strongly. We use our model to predict fluctuations of borehole-water levels that result from different types of disturbances. We show how this framework can be used to estimate subglacial hydraulic properties by comparing model results with data collected during field experiments on Trapridge Glacier, Yukon Territory, Canada in 1989 and 1990.

Image-analysis procedures have been developed for a simple PC-based system to evaluate the textural features of ice thin sections automatically. The method yields parameters that describe both ice crystals and pores in an objective, reproducible manner. Specimens are recorded in linearly and circularly polarized light (ice grains) as well as in transmitted and incident plain light (pores) under standardized conditions. After preparative filtering, images are segmented through union of two Sobel-filtered images (ice crystals/grain boundaries) or through thresholding (ice crystals/pores). The quantification of texture is based on the evaluation of image contrast as well as on linear analysis and non-linear transforms to obtain information on grain-sizes and their distribution. In the study of ice pores, stereological effects of finite section thickness have to be corrected. The discussion focuses on errors involved in automated textural analysis and the glaciological yield of measured and derived parameters.

For our work in the Filchner-Ronne Ice Shelf Programme (FRISP), we have developed a new technique for measuring the bottom-melting rate with high reliability. The method is based on time-domain reflectometry (TDR) measurements of transmission lines inserted into melt holes. The TDR-data are digitally recorded on magnetic tape. System resolution has been estimated at 0.2 m. Hence, re-measuring after 1 year gives an accuracy of 10% for melting rates of 2m a-1. Two transmission lines for TDR measurements were installed during the German FRIS Expedition 1989–90. This paper describes the design of the system. Examples of recorded wave forms are given.

This paper is concerned with the motion of an unconfined finite mass of granular material down an inclined plane when released from a rest position in the shape of a circular or elliptical paraboloid. The granular mass is treated as a frictional Coulomb-like continuum with a constant angle of internal friction. The basal friction force is assumed to be composed of a Coulomb-type component with a bed-friction angle that is position-dependent and a viscous Voellmy-type resistive stress that is proportional to the velocity squared. The model equations are those of Hutter and others (in press b) and form a spatially two-dimensional set for the evolution of the avalanche height and the depth averaged in-plane velocity components; they hold for a motion of a granular mass along a plane surface.

Similarity solutions, i.e. solutions which preserve the shape and the structure of the velocity field, are constructed by decomposing the motion into that of the centre of mass and the deformation relative to it. This decomposition is possible provided the effect of the Voellmy drag on the deformation is ignored. With it, the depth and velocities relative to those of the centre of mass of the moving pile can be determined analytically. It is shown that the pile has a parabolic cap shape and contour lines are elliptical. The semi-axes and the position and velocity of the centre of mass are calculated numerically. We explicitly show that

• (i) For two-dimensional spreading, a rigid-body motion does not exist, no matter what be the values of the bed-friction angle and the coefficient of viscous drag.

• (ii) A steady final velocity of the centre of the mass cannot be assumed, but the motion of the centre of mass depends strongly on the value of the Voellmy coefficient.

• (iii) The geometry of the moving pile depends on the variation of the bed-friction angle with position, as well as on the value of the coefficient of viscous drag.

The results of measurements of radio-wave velocities (RWV) by wide-angle reflection (WAR) methods in the temperate Abramov Glacier in the Alai Mountain Ridge and the “two-layered” sub-polar Fridtjovbreen and Hansbreen on Svalbard using a low-frequency (2–13 MHz) radar are considered and discussed. The experimental data obtained and the data from the literature show that the values of RWV could be a good indicator of the hydrothermal state of glaciers. As such, these data enable the identification of cold, temperate and transitional (two-layered) glaciers, and can be used for estimation of the water content in glaciers and changes in the hydrothermal state.

The Seasat A satellite scatterometer radar, initially designed to measure ocean-wind intensity and direction, also provided observations on the Antarctic ice sheet. The signal of the back-scatter coefficient decreases strongly from 10 to −20 dB when the incidence angle of the observations increases from 0° to 65°. An additional 5 dB signal is found, which is correlated with the direction and intensity of katabatic winds, independent of the incidence angle and polarization of the signal. By using simplified models of the volume-scattering within the snowpack (which is mostly sensitive to snow grain-size) and surface-scattering from the air-snow interface (which depends on roughness), it is evident that the signal of the scatterometer could result from the effects of snow dunes at low incidence angle, and of micro-roughness and volume back-scatter at incidence angles greater than 25°. The instrument therefore provides a means of measuring the direction and intensity of katabatic winds.

Numerical models to assess the principal response of large ice caps to climatic changes are used as a guide to the interpretation of field evidence of changes in the glaciological regime in the coastal areas of Victoria Land and north Victoria Land, Antarctica. Based on the theoretical work, the following scenario is proposed: areas within about 300 km from the coast of Victoria Land experienced (i) significantly shallower ice slopes and a lesser degree of glaciation during most of the late Tertiary, (ii) steep slopes and thicker ice than today during glacial stages, and (iii) moderate thinning of the ice in intervening interglacial stages.

The model predicts, for central regions of Antarctica, a slightly thinner ice cap (lower snow-accumulation rate) during glacial stages, but an estimated 200 m thicker ice cap in warmer Tertiary climates than today. The calculated “Tertiary ice caps” indicate a probable tendency of periodic surges due to basal melting at the outer fringes.

Only modest changes of the ice thickness in reasonably good agreement with the model predictions for the current interglacial stage have been observed on four blue-ice fields, all located within 250 km of the coastline. Investigated ice fields include two meteorite traps — the Allan Hills Icefield and the Frontier Mountain meteorite fields. Antarctic meteorite traps are sustained by very specific glaciological conditions — in particular, only moderate changes in ice thickness over time. The sub-ice topography of these fields was determined by radar measurements and reveals a former, very different glaciological regime, which is interpreted as being associated with glacial processes, operative during the late Tertiary.

Field evidence for a hypsithermal event during the Holocene is presented.