Hobson’s Choice Ice Island is a tabular iceberg that calved in 1982-83 from East Ward Hunt Ice Shelf, N.W.T., Canada. Four ice cores have been analyzed for ice-crystal size, structure and fabric, bulk density, liquid electrical conductivity, δ18O and tritium. This has enabled a complete characterization of the physical properties and the structural characteristics of the ice-shelf component of Hobson’s Choice Ice Island and the first ever study of the stratigraphy and growth history of East Ward Hunt Ice Shelf. The δ18O values range from -34.6 to and indicate that all the ice is derived directly and/or indirectly from precipitation. High tritium values occur only in the lowermost 5 m of the ice shelf in a layer named stratum B. The tritium is anthropogenic and indicates bottom accretion of fresh-water ice since 1952, most likely from fresh water flowing out of Disraeli Fiord below the eastern ice shelf. Above this deepest and youngest ice layer is a 35-38 m thick, unconformable layer (stratum A) comprising three ice types: iced-firn, slush ice and lake ice. This depositional-superimposed ice represents past surface accumulation, which, according to δ18O and ice-crystal structure and size variations, occurred in three major periods, each interrupted by major ablation periods. Fresh water flowing out of Disraeli Fiord below the ice shelf during those warm intervals was the most likely agent responsible for the ablation and eventual complete loss of the original sea-ice platform on which stratum A initially accumulated. The three sub-strata of stratum A vary in thickness from core to core and suggest that there has been an inversion of relief during each ablation period. The different properties and occurrences of the three superimposed ice types are due primarily to past topographic variability.

Results from a detailed profile in a 5.54 m multi-year sea-ice core from the rift area in the southern part of George VI Ice Shelf are presented. Stratigraphy, stable isotopes and Na content are used to investigate the growth processes of the ice cover and to relate them to melting processes at the bottom of the ice shelf.

The thickest multi-year sea ice in the sampling area appears to be second-year sea ice that has survived one melt season. Combined salinity/stable-isotope analyses show large-scale sympathetic fluctuations that can be related to the origin of the parent water. Winter accretion represents half of the core length and mainly consists of frazil ice of normal sea-water origin. However, five major dilution events of sea water, with fresh-water input from the melting base of the ice shelf reaching 20% on two occasions, punctuate this winter accretion. Two of them correspond to platelet-ice production, which is often related to the freezing of ascending supercooled water from the bottom of the ice shelf.

Brackish ice occurs between 450 and 530 cm in the core. It is demonstrated that this results from the freezing of brackish water (Jeffries and others, 1989) formed by mixing of normal sea water with melted basal shelf ice, with dilution percentages of maximum 80% fresh water.

A new method of hydrograph separation for bulk meltwaters draining Alpine glaciers is proposed. It is based on the two-component (subglacial and englacial) mixing model of Collins (1978), but allows the composition of the subglacial component to vary between ascending and descending lines of the hydrograph. The mean englacial component can be derived from linear relationships between sulphate concentrations and other ions in bulk meltwaters. On certain occasions during the ablation season, the maximum concentration of ions in the subglacial component can be determined from the linear relationship between bulk meltwater sulphate concentrations and discharge. The bulk discharge is then a direct measure of the mass fraction of the englacial component. At maximum discharge, the contribution of the subglacial component approaches zero, which has implications for the storage and mixing of waters in subglacial reservoirs. Further, the subglacial component is not of constant composition, and may itself be a mixture of dilute supraglacial and concentrated subglacial water.

Englacial temperature measurements in Arctic valley glaciers suggest in the ablation zone the existence of a basal layer of temperate ice lying below the bulk of cold ice. For such a polythermal glacier, a mathematical model is presented that calculates the temperature in the cold part and the position of the cold-temperate transition surface (CTS). The model is based on the continuum hypothesis for ice and the ice-water mixture, and on the conservation laws for moisture and energy. Temperate ice is treated as a binary mixture of ice and water at the melting point of pure ice. Boundary and transition conditions are formulated for the free surface, the base and the intraglacial cold-temperate transition surface. The model is reduced to two dimensions (plane flow) and the shallow-ice approximation is invoked. The limit of very small moisture diffusivity is analysed by using a stationary model with further reduction to one dimension (parallel-sided slab), thus providing the means of a consistent formulation of the transition conditions for moisture and heat flux through the CTS at the limit of negligibly small moisture diffusion.

The application of the model to the steady-state Laika Glacier, using present-day conditions, results in a wholly cold glacier with a cold sole, in sharp contrast to observations. The present polythermal state of this glacier is suspected to be a remnant of the varying climatic conditions and glacier geometry during the past few centuries. Steady-state solutions representing a polythermal structure can indeed be found within a range of prescribed conditions which are judged to be realistic for Laika Glacier at the last maximum extent of the glacier.

Jakobshavns Isbræ is a large, fast-moving ice stream/outlet glacier in West Greenland which ends at a floating, calving front. It drains about 6.5% of the area of the Greenland ice sheet. Studies of its surface morphology are described in this paper. The surface is relatively steep (0.01–0.03) and the thickness is large (up to 2600 m along the center line (Clarke and Echelmeyer, 1989)), indicating very high driving stresses (200–300 kPa). The ice stream is about 6 km wide and 85–90 km long, all of which is in an area of surface melting. The base of the ice stream, and of much of the drainage area, is below sea level. Marginal crevasse zones have a width on the order of the width of the ice stream itself. Unique surficial features are ice blisters and lakes; the latter have a sequence of ogive-like features on their floating ice cover which can be used to determine velocity. There is a pinning point near the terminus which may act as a stabilizing influence, possibly playing a role in halting, at least temporarily, a recent retreat of the terminus. Ice-thickness estimates at the terminus lead to a flux which is less than previously assumed by others (e.g. Bindschadler, 1984; Pelto and others, 1989) when estimating Jakobshavns Isbræ’s drainage basin to be nearly in balance.

The driving stresses on Jakobshavns Isbræ are an order of magnitude higher than those of the ice streams of West Antarctica. Its crevasse patterns are much less localized. Its relatively unconfined terminus is more comparable to that of relatively unbuttressed ice streams such as Pine Island and Thwaites Glaciers than it is to other West Antarctic ice streams which terminate in large, confined ice shelves.

The widespread, uniform till sheets of the southern margin of the Laurentide ice sheet were deposited by fast-moving, wet-bedded ice over relatively short periods of time. Characteristics of these till sheets are entirely consistent with deposition from deforming subglacial sediment layers; it is difficult (although not impossible) to explain their origin through debris transport in basal ice. Careful estimation of debris budgets and studies of sorted sediments contained in the till sheets may clarify their origin further.

The effects of notch acuity (crack-tip sharpness) on the fracture toughness of S2 ice were investigated using six groups of single-edge-notched-bend (SENB) specimens with different crack (or notch) root radii fabricated by six different methods. The mean value and standard deviations of the apparent fracture-toughness values KQ of the specimens with blunt notches were significantly higher than those of the specimens with sharp cracks. The results presented in a plot of fracture toughness versus , where p is the crack-tip radius, provide an estimate of the required notch acuity for fracture-toughness testing. The microstructural features in the immediate vicinity of the crack tip and the crack-tip geometry govern the magnitude of the apparent fracture toughness of the ice. The crack-tip-opening-displacement (CTOD) measured in this study indicates the existence of microplastic deformation in the vicinity of the crack tip at the initiation of unstable fracture.

We report on laboratory experiments on the motion of powder-snow avalanches along a bent chute. The avalanches are simulated as turbulent boundary-layer flows of polystyrene particles in still water along a chute consisting of a straight inclined part, a curved part and a second, possibly inclined, run-out zone. An ultrasonic measuring technique is used to determine mean particle speeds (via the Doppler shift of the reflected signal) and the particle concentration (via the attenuation of the echoes). By ensemble averaging, individual turbulence effects are eliminated. As measuring procedures, profiles were determined for particle velocity and density across the boundary layer; these were taken (i) for the avalanche tail along the entire track, i.e. in the steep part and the run-out zone, and (ii) for the avalanche head in the run-out zone below the kink in the terrain. Moreover, time sections of velocity and density (i.e. time series at fixed positions along a line through the boundary layer) were recorded and the particle mass deposited in the various zones of the track was measured.

The analysis of the data reveals the following results: a concave change in the terrain topography acts as a very efficient mechanism for particle sedimentation, thus affecting particle concentration and velocity, and considerable reduction close to the ground, whereas the density is reduced throughout the depth, thus leading to a considerable reduction of dynamic pressure close to the ground and leaving it somewhat greater at higher altitudes. We conclude with practical considerations for the field glaciologist.

Annual accumulation records were derived from two firn cores drilled at GD03 (69° 00’S, 115° 30’E, 1835m) and GD15 (69° 00’S, 130° 48’E, 2155m), which are separated by 750 km along the 2000 m contour in eastern Wilkes Land. The accumulation records, spanning the period 1930–85, show that annual accumulation has increased by 25% from 1960 to 1985, resulting in the highest accumulation rates in the 55 year records. Annual accumulation layers were identified in the firn cores using a firn-stratigraphic model (based on the observed characteristics and variability of the seasonal snow-surface layer) together with seasonal oxygen-isotope ratios. The accumulation records support other evidence for a similar recent increase across Antarctica. The greater increase over the long-term mean for 1930–85 in eastern Wilkes Land was observed at GD15 which is located on a northeast aspect slope near Terre Adèlie. It is suggested that the increase is the result of changes in the general atmospheric circulation pattern which have produced a higher annual frequency of precipitation events.

Surface and basal long profiles are reconstructed for 13 outlet glaciers that drained ice from a large ice field (80 km by 120 km) that formed in the western Grampians of Scotland during part of the Late-glacial period (c, 14000–10000years BP). Basal shear stresses are calculated at 5 km intervals along the central flowlines of the reconstructed outlet glaciers. Individual basal shear stresses for the outlet glaciers range from 10 to 204 kPa. Variations in calculated basal shear stresses within and between the glaciers are mainly explained by differences in bedrock topography, extending and compressional flow, and by differences in basal boundary conditions. Low basal shear stresses (<53kPa) calculated for the terminal zones of Creran, Menteith and Lomond glaciers are partly explained by the overriding of glaciomarine clays with inferred high pore-water pressures and a low yield strength that may have led to rapid basal sliding and thinning of the ice lobes.

Quantitative understanding of the processes that couple the lower atmosphere to the upper surface of ice sheets is necessary for interpreting ice-core records. Of special interest are those processes that involve the exchange of energy or atmospheric constituents. One such process, wind pumping, entails both possibilities and provides a possible mechanism for converting atmospheric kinetic energy into a near-surface heat source within the firn layer. The essential idea is that temporal and spatial variations in surface air pressure, resulting from air motion, can diffuse into permeable firn by conventional Darcy flow. Viscous friction between moving air and the solid firn matrix leads to energy dissipation in the firn that is equivalent to a volumetric heat source.

Initial theoretical work on wind pumping was aimed at explaining anomalous near-surface temperatures measured at sites on Agassiz Ice Cap, Arctic Canada. A conclusion of this preliminary work was that, under highly favourable conditions, anomalous warming of as much as 2°C was possible. Subsequent efforts to confirm wind-pumping predictions suggest that our initial estimates of the penetration depth for pressure fluctuations were optimistic. These observations point to a deficiency of the initial theoretical formulation — the surface-pressure forcing was assumed to vary temporally, but not spatially. Thus, within the firn there was only a surface-normal component of air flow. The purpose of the present contribution is to advance a three-dimensional theory of wind pumping in which air flow is driven by both spatial and temporal fluctuations in surface pressure. Conclusions of the three-dimensional analysis are that the penetration of pressure fluctuations, and hence the thickness of the zone of frictional interaction between air and permeable firn, is related to both the frequency of the pressure fluctuations and to the spatial coherence length of turbulence cells near the firn surface.

Energy exchange between the atmosphere and a melting glacier surface is mediated by the presence of a water layer. Under conditions of rapid melt and/or heavy rainfall, the possibility exists that a supraglacial run-off layer can advect sensible heat and influence the spatial variations of melt. The potential magnitude of such advection was investigated by numerically solving differential equations expressing the mass and energy balances of a two-dimensional run-off layer. Solutions were obtained for conditions typical of rainfall events, in which the potential for supraglacial heat advection should be maximal. The solutions indicate that advection cannot influence macro-scale melt patterns and surface morphology, except perhaps under heavy rainfall and/or rapid melt conditions, but can possibly cause micro-scale variations in ice melt. One-dimensional energy-balance models, which have normally been applied over glacier surfaces, should remain valid for most conditions.

Many glaciers in Svalbard and in other glacierized areas of the world are known to surge. However, the time series of observations required to assess the duration of fast motion is very restricted. Data on active-phase duration in Svalbard come from aerial photographs, satellite imagery, field surveys and airborne reconnaissance. Evidence on surge duration is available for eight Svalbard ice masses varying from 3 to 1250 km2. Worldwide, active-phase duration is recorded for less than 50 glaciers. Few observations are available on high polar ice masses. The duration of the active phase is significantly longer for Svalbard glaciers than for surge-type glaciers in other areas from which data are available. In Svalbard, the active phase may last from 3 to 10 years. By contrast, a surge duration of 1–2 years is more typical of ice masses in northwest North America, Iceland and the Pamirs. Ice velocities during the protracted active phase on Svalbard glaciers are considerably lower than those for many surge-type glaciers in these other regions. Mass is transferred down-glacier more slowly but over a considerably longer period. Svalbard surge-type glaciers do not exhibit the very abrupt termination of the active phase, over periods of a few days, observed for several Alaskan glaciers. The duration of the active phase in Svalbard is not dependent on parameters related to glacier size. The quiescent phase is also relatively long (50–500 years) for Svalbard ice masses. Detailed field monitoring of changing basal conditions through the surge cycle is required from surge-type glaciers in Svalbard in order to explain the significantly longer length of the active phase for glaciers in the archipelago, which may also typify other high polar ice masses. The finding that surge behaviour, in the form of active-phase duration, shows systematic differences between different regions and their environments has important implications for understanding the processes responsible for glacier surges.

The dynamics of the glaciers which form the headwaters of the Susitna River in central Alaska exhibit several interesting features, including a spectrum of surge-type behavior. A difference between balance flux and actual down-glacier transport, which is taken to be an indicator of surge behavior, shows West Fork Glacier and two tributaries of Susitna Glacier to be surge-type of varying strengths, while East Fork Glacier and one tributary of Susitna Glacier to be non-surge-type. The main trunk of Susitna Glacier and its two unstable tributaries surge simultaneously with a period estimated to be 50–60 years. Having last surged in 1951–52, its next surge should be expected sometime in the first decade of the next century. West Fork Glacier last surged in about 1935 and again in 1987–88, indicating a similar surge period of about 50 years. Significant seasonal velocity variations were observed during the glacier’s quiescent phases, with increases of 30–100% over background occurring early in the melt season. In some cases, the annual minimum occurred during late summer, implying that basal motion contributes measurably to winter velocity on at least some glaciers in the area. Events of rapid motion lasting a day or so also occur occasionally during the melt season and may account for brief velocity increases of up to 300%.

Water is present in glaciers in the form of veins at the three-grain junctions. This water remains unfrozen even many degrees below the normal freezing point, mainly because it contains much of the soluble impurity content of a glacier, but also because of the microscopic curvature of the ice–water interfaces. As the temperature is lowered and the veins shrink, the concentration of impurities in them increases, and the curvature effect also increases. The predicted relation between vein size and temperature has now been verified by laboratory experiments.

Because of the latent heat of the vein water, the ice behaves macroscopically as a continuum with an anomalous specific heat capacity that depends strongly on temperature. From this point of view, a polythermal glacier is a single medium with continuously varying properties, rather than consisting of distinct cold and temperate phases with sharp boundaries between them. The paper sets up differential equations for heat diffusion in such a continuum. To explain the local uniformity of the vein system seen under the microscope, it is found necessary to include the effect of diffusion of solutes along the veins.

Solutions are presented for a model in which two semi-infinite slabs, initially having different temperatures, impurity concentrations and vein sizes, are instantaneously brought into contact. In this way, transition thicknesses between cold and temperate ice are estimated, and also the velocities of various kinds of waves that are generated from the original discontinuity at the interface.

Optical measurements have been made on the water lenses which form under pressure at grain boundaries in polycrystalline ice. Monochromatic light from a point source is focused by the lenses but, because the lenses are microscopic in size, the image is blurred by diffraction. The diffraction pattern observed under a microscope has been compared with the computed diffraction pattern to deduce the angle 2θ at the rim of each lens. This is the dihedral angle for water at a grain boundary in ice, and gives the ratio of the grain-boundary energy to that of an ice-water interface. The most sensitive measurements are those made on the rings of the virtual diffraction pattern formed on the object side of the lens. They give θ = 12.5 ± 0.5° for the grain boundary under observation, which is 26% lower than the previous value for θ found by ignoring diffraction.

A technique to estimate undisturbed ice temperature is discussed for sensors placed in boreholes that have been heated to the melting point during drilling, and for which only a limited time span of temperature record is available. A short temperature record after the hole refreezes commonly results when using hot-water or steam drills, where measurements are constrained by logistics, ice deformation, sensor drift or other problems, or where the refreezing time is long because of near-freezing ice temperatures or large hole sizes. Short data records are also typical in ongoing drilling programs where temperature information may be necessary for the program itself. Building on analyses by Lachenbruch and Brewer (1959) and a numerical model by Jarvis and Clarke (1974), it is shown that estimates of undisturbed temperatures can be made from records of temperature that extend only marginally beyond the initial refreezing. Complex effects of hole size, heating history, and the thermodynamic and geometrical effects of a moving boundary (the freezing borehole walls) are important to temperature decay immediately after freeze-up, so that the standard technique of comparing temperature decay to an inverse of time model is not applicable, and comparsion has to be made to a numerical model of heat flow to a refreezing borehole. Data from Ice Stream B, Antarctica, are compared to the numerical model to illustrate the technique. Data are also compared to simpler (inverse time) thermal models, and a potential for error is pointed out, since a short data record can be spuriously matched with the simpler, one or two free-parameter, models.

Thirty-nine samples of suspended sediment from the meltwater stream of Austre Okstindbreen, Norway, were analysed using both moments and parameters of the hyperbolic distribution. The samples turned out to be non-normally distributed. No direct correlation was found between discharge and grain-size distribution, but it is supposed that the hydraulic influence is hidden behind factors like composition of morainic or stream-bed material. Through comparison of the grain-size parameters from two periods in the field season, it was possible to distinguish between two sediment populations.

The authors describe a new method that allows an accurate, cheap and quick characterization of snow grains from natural snowpacks. The method is based on stopping all metamorphism of snow grains by immersing them in cold isooctane. Snow grains are then carried from the field, stored in a cold room and easily observed through a microscope. It then makes possible recording of snow-grain pictures in a state similar to their natural state, even after storage of several months.

The different snow and ice types on a glacier may be subdivided according to the glacier-facies concept. The surficial expression of some facies may be detected at the end of the balance year by the use of visible and near-infrared image data from the Landsat multispectral scanner (MSS) and thematic mapper (TM) sensors. Ice and snow can be distinguished by reflectivity differences in individual or ratioed TM bands on Brúarjökull, an outlet glacier on the northern margin of the Vatnajökull ice cap, Iceland. The Landsat scene shows the upper limit of wet snow on 24 August 1986. Landsat-derived reflectance is lowest for exposed ice and increases markedly at the transient snow line. Above the slush zone is a gradual increase in near-infrared reflectance as a result of decreasing grain-size of the snow, which characterizes drier snow. Landsat data are useful in measuring the areal extent of the ice facies, the slush zone within the wet-snow facies, the snow facies (combined wet-snow, percolation and dry-snow facies), and the respective positions of the transient snow line and the slush limit. In addition, fresh snowfall and/or airborne contaminants, such as soot and tcphra, can limit the utility of Landsat data for delineation of the glacier facies in some cases.