Two sets of bubble-density foliation, one vertical and transverse to flow and one vertical and parallel to flow, developed in the ice of the terminus lobe of Variegated Glacier, Alaska, during the passage of the final stages of the 1982–83 surge in June and July 1983. The transverse foliation is ubiquitous in the lower part of the terminus lobe affected by the surge, but disappears abruptly over a 5 m wide zone at the boundary of the surge, as delineated by measurable strain rates during the surge. The ice at the outermost limit of the transverse foliation experienced compressive total strains no greater than −0.10 (logarithmic), but did briefly experience high strain rates, corresponding to compressive stresses in excess of 0.6 MPa. It is concluded that the transverse foliation was developed as a consequence of high compressive stress rather than high total strains. Experimental evidence of Kamb indicates that bubble migration occurs under similar tectonic conditions to those in the terminus lobe during the surge, and that bubble migration is related to stress rather than finite deformation or rotation of pre-existing structures. Bubble migration may be the mechanism by which the transverse foliation developed during the passage of the surge front.

In the summer of 1989 a glacio-meteorological experiment was carried out on Hintereisferner, Austria. Here, we discuss the interpretation of detailed ablation measurements on the tongue. Two profiles are considered: one along the glacier axis and one across the glacier tongue. At the crossing, an energy-balance station was operated, providing the necessary meteorological data. At all stakes the albedo was measured several times with a simple light-weight instrument (described separately in the Appendix). The absorption of solar radiation was calculated for all stake positions, taking into account local albedo, slope, aspect and shading effects. It appears that the differences in ablation can be explained almost entirely by the differences in surface albedo, both along and across the glacier tongue.

This paper describes fine-resolution measurements of glacier surface strain rates very close to the margin of Russell Glacier, West Greenland. Measurements at a small scale make possible detailed analysis of strain patterns close to the glacier margin, and suggest that strain rates vary over small areas. The strain pattern is determined by ice flexure over subglacial obstacles as well as by seasonally variable marginal retardation and by the orientation of the ice margin relative to the flow direction.

Three glaciers in southern Norway, with very different mass-balance characteristics, are studied with an energy-balance model of the ice/snow surface. The model simulates the observed mass-balance profiles in a satisfactory way, and can thus be used with some confidence in a study of climate sensitivity. Calculated changes in equilibrium-line altitude for a 1 K temperature increase are 110, 108 and 135 m for Nigardsbreen, Hellstugubreen and Alfotbreen, respectively. The corresponding changes in mass balance, averaged over the entire glacier area, are −0.88, −0.715 and −1.11 m year−1 (water equivalent).

Runs with an ice-flow model for Nigardsbreen, to which calculated mass-balance profiles arc imposed, predict that the front will advance by 3 km for a 1 K cooling, and will retreat by as much as 6.5 km for a 1 K warming. The response to a 10% increase in precipitation would be a 2 km advance of the snout, whereas a 4 km retreat is predicted for a 10% decrease. This large sensitivity (as compared to many other glaciers) is to a large extent due to the geometry of Nigardsbreen.

Although quantitative interpretation of the low-frequency electrical conductivity of ice cores from central Greenland is complicated by temperature variations of the measured core, annual layers can be recognized in sections of the core that are not impacted by non-seasonal features. Ambiguities in counting of annual layers can be minimized by comparing the electrical conductivity measurements to measurements of dust concentration and visual stratigraphy. A non-linear relationship between applied voltage and the current measured across two electrodes complicates the quantitative comparison of measurements made with different equipment, but does not affect the overall shape of the observed features.

Temperature compensation built into many conductivity meters becomes inaccurate in low-temperature waters typical of glacial melt streams. Experiments using simulated glacial waters show a linear dependence of electrical conductivity on temperature from 0.3° to 25°C. The slope of the temperature–conductivity relation is linearly dependent on the conductivity of the solution, allowing a numerical or analytical temperature correction to be made. To minimize error introduced by temperature compensation, measurements of electrical conductivity in glacial streams should be corrected to a low standard temperature, and 0°C is suggested.

Laws of snow metamorphism have been introduced in a numerical model which simulates the evolution of temperature, density and liquid-water profiles of snow cover as a function of weather conditions.

To establish these laws, the authors have summarized previous studies on temperature gradient and on wet-snow metamorphism and they have also conducted metamorphism experiments on dry or wet fresh-snow samples. An original formalism was developed to allow a description of snow with parameters evolving continuously throughout time.

The introduction of laws of metamorphism has improved significantly the derivation of the settlement of internal layers and of snow-covered albedo, which depend on the simulated stratigraphy, i.e. the type and size of snow grains of different layers of the snow cover.

The model was tested during a whole winter season without any re-initialization. Comparison between the simulated characteristics of the snow cover and the observations made in the field are described in detail. The model proved itself to be very efficient in simulating accurately the evolution of the snow-cover stratigraphy throughout the whole winter season.

The geometry of the vein system in ice has been investigated using photographs of enlarged veins in ice samples that were grown in the laboratory. The veins, which are non-uniform, act as tiny triangular-shaped, water-filled prisms that refract the light passing through them.

The three vein widths in the cross-section of a vein can be deduced from two photographs taken from different directions. The dihedral angle along a given vein edge can be observed directly by viewing it at a node, where four veins meet, from a particular direction. The dihedral angles range from 25° ± 1° to 105° ± 1°. It is shown that the vein cross-section can be constructed, given the three widths of a vein and one of the dihedral angles, providing that the radius of curvature around the vein wallsrv is a constant. This assumption can be checked if the values of at least two of the dihedral angles associated with the vein cross-section are known. Ifrv is a constant, then the solid-liquid interfacial energy ϒsl must be isotropic for the veins in question and any deviations from uniform equilibrium geometry must derive primarily from anisotropy in the grain-boundary energy ϒss. The cross-sections of three veins that meet in a particular node are constructed. The assumption of isotropic ϒsl is found to hold for this node.

In this study, an attempt is made to simulate the historical length variations of Hintereisferner (Austria). Key elements in the model are: (1) reconstruction of the local climate (Vent); (2) multiple linear regression analysis between measured climate and mass-balance variations; (3) reconstruction of the mass balance with (1) and (2); (4) a one-dimensional numerical ice-flow model forced by the mass-balance history. Progress relative to the previous, similar studies can be summarized as two points: (1) the mass balance is reconstructed more carefully, and (2) sensitivity experiments show that uncertainties and assumptions in the flow model hardly affect the response of the glacier to climate fluctuations with a time-scale of the order of (or larger than) the response time of the glacier (around 100 year).

In the regression equation, mass-balance fluctuations are determined by the temperature during the ablation season and the annual precipitation. Cloudiness has no effect. The local temperature can be reconstructed as far back as 1755 with the temperature series from Basel, but the reconstruction of precipitation is unreliable prior to 1900. Mainly due to this point, a simulation starting in 1755 failed completely. However, in a run starting in 1894, the total retreat since then (2.2 km) was reproduced within 10%.

In 1985, 1987 and 1988, we measured deformation of boreholes in the ablation zone of Storglaciären. A different hole was used each year. Deformation profiles were determined for four consecutive time periods between mid-July and early September. All three holes were in an area where the glacier lies in an overdeepened part of the bed, up-glacier from a riegel. The first hole was approximately on the glacier center line while the others were midway between the center line and one or the other margin. Despite large standard errors, variations in the deformation profiles during the melt season are broadly consistent with other data. We thus believe that they are real, and that they reflect changes in the stress field resulting from changes in drag at the bed. A decrease in drag early in the melt season, coupled with convergence of the bed contours towards the riegel, is believed to be responsible for a previously documented increase in vertical velocity just up glacier from the riegal at this time of year. Also noteworthy, but perhaps problematical, was an abrupt transverse shear of at a depth of ~25 m near the center of the glacier in July 1985. The shear took place in the zone in which cold ice near the surface gives way to temperate ice deeper in the glacier.

Ice-flow modeling requires a flow law relating strain rates to stresses in situ, but a flow law cannot be measured directly in ice sheets. Microscopic processes such as dislocation glide and boundary diffusion control both the flow law for ice and the development of physical properties such as grain-size and c-axis fabric. These microscopic processes can be inferred from observations of the physical properties, and the flow law can then be estimated from the microscopic processes.

A review of available literature shows that this approach can be imperfectly successful. Interior regions of large ice sheets probably have depth-varying flow-law “constants”, with the stress exponent, n, for power-law creep less than 3 in upper regions and equal to 3 only in deep ice; n probably equals 3 through most of the thickness of ice shelves and ice streams.

Remote-sensing methods are the primary ones used for ice mapping in the Baltic Sea. A major methodological improvement is now being introduced by satellite radars due to their weather independency and high resolution. To learn how to use ERS-1 synthetic aperture radar (SAR) data, an extensive field programme BEPERS (Bothnian Experiment in Preparation for ERS-1) with airborne SARs has been arranged. The BEPERS pilot study was undertaken in 1987 using the French VARAN-S X-band SAR. The SAR was flown on 1 day over four study areas of size approximately 10 km x 50 km, and intensive validation observations were made. The data were most useful for the education they provided on how to work with SAR in sea-ice mapping. They have been used for developing SAR image-analysis methods, back-scatter modelling investigations and geophysical validation of SAR imagery. Cleaning-up of images consisted of speckle reduction and segmentation. Back-scatter characteristics of undeformed ice and ridges were examined. Ice-type classification was based on the box-classification method. Eight ice types were defined but basically only two types, undeformed ice/open water and deformed ice, could be discriminated. Two basic problems of high practical importance remained: how to discriminate between (1) open water and undeformed ice, and (2) ridged ice and brash ice. The data further showed illustrative examples of SAR imagery over sea ice.

Experimental studies are reported concerning the thermal behaviour of the water-vein system in ice grown in the laboratory from dilute solutions. The temperature versus vein-size behaviour of these samples is determined. The measurements show that the solute, which is concentrated in the veins, remains in the liquid phase during a temperature change. The mass of the solute per unit length of vein M is found to be of the order of M ≈ 10−8 mol m−1 for samples grown from singly-distilled water. M is seen to vary with temperature only because of the volume expansion (contraction) on freezing (melting) which causes the liquid to flow along the veins. The effect of these flows on the sample is studied. They are found to provide a mechanism for the transport of impurities along the veins and to and from the sample surface. Samples grown from doubly-distilled water doped with small amounts of NaCl or H2SO4 are studied and are found to display the same general behaviour. However, M is an order of magnitude higher in the H2SO4-doped samples than in either the NaCl-doped samples or the samples grown from singly-distilled water. The approach to equilibrium of distortions in the vein-system geometry is studied. It is suggested that these distortions are due to variations in M along the length of the veins and that equilibration is therefore governed by diffusion of the solutes.

Two practical problems in the use of time-lapse photography for the measurement of speed were encountered during the recent surge of West Fork Glacier in the central Alaska Range, Alaska, U.S.A. The first is severe rotational camera instability; we show how natural, unsurveyed features on the valley wall can be used to make the necessary corrections. The second problem is the computation of absolute speed when many different, unsurveyed glacier-surface features are used as targets. We give a method for connecting the data obtained from different targets, and for determining the scale using limited information obtained by surveying. Severe systematic errors can occur unless the angle between the axis of the lens and the direction of horizontal motion is determined.

The long-term monitoring of Lewis Glacier on Mount Kenya serves as a basis for the study of glacier evolution in response to climatic forcing through modeling of its ice flow and mass budget. Following up on an earlier modeling and prediction study to 1990, this paper examines the ice-mass and flow changes in relation to the net-balance conditions over 1986–90. A model experiment using as climatic forcing the observed 1978–86 vertical net-balance profile yielded a volume loss and slow down of ice flow more drastic than observed during 1986–90. The causes of this discrepancy were examined in successive model experiments. Realistic simulations of mass-budget and thickness changes over 1986–90 are obtained using as input the net-balance forcing for the same period rather than for the preceding 1978–86 interval, and approximate flow velocities. With the same net-balance forcing and a completely stagnant Lewis Glacier, the elimination of mass redistribution by ice flow acts to mitigate the loss of volume and thickness in the upper glacier, and to accentuate it in the lower glacier. Accordingly, the observed 1986–90 net-balance profile along with the 1990 ice-flow velocities provide suitable input for the modeling of Lewis Glacier changes to 1994. Under continuation of the 1986–90 climatic forcing, ice thinning ranging from less than 1 m in the upper glacier to more than 7 m in the lower glacier, and a total volume loss of order 57 × 104 m3, are anticipated over the 1990–94 time interval.

Systematic isotopic studies based on natural and artificial radio-isotopes (32Si, 137Cs), stable isotopes (δ18O) and total β activity measurements have been carried out on Chhota Shigri glacier, Himachal Pradesh, central Himalaya, to study the dynamics of the ice, meltwater composition and to identify the deposition of the Chernobyl fall-out in the Himalayan region.

Using 32Si concentrations, the snout ice has been dated at ~ 250 years, based on which the past average surface ice-flow rate has been estimated as ~ 28 m year−1. Based on δ18O variations, in a shallow ice core, the accumulation rate of the ice has been estimated at ~520kgm−2 year−1. 32Si measurements of snout ice and englacial meltwaters indicate that at least 55% of the snow meltwater mixed with 45% of the old ice-melt water that emerged from englacial streams in the month of August 1987. Deposition of the artificial radionuclide (137Cs) and the very high total β activity observed in snow samples on Chhota Shigri glacier give the first evidence of Chernobyl fall-out deposition in the Indian Himalaya.

Annual ice ablation has been measured at three locations at the margin of the Greenland ice sheet: Nordbogletscher (for 6 years: 1977–78 to 1982– 83), Qamanârssûp sermia (8 years: 1979–80 to 1986–87) and at Paakitsup Akuliarusersua (7 years: 1982–83 to the present). As the data sets cover different periods, it is difficult to compare ablation directly between the three sites. However, measured series at each site can be extended over the last 30 years (1961–90) by simulations using climatic data, and the extended data series can be compared. The pattern of calculated ablation variations at the three locations is remarkably similar with no sign of any trend towards increased ablation in recent years. There was generally low ablation from the mid-1970s to the mid-1980s that may explain the recent thickening of the ablation area which has been detected by satellite radar altimetry. Ablation varies substantially from year to year, e.g. with a standard deviation of the order of ±0.5 m water a−1, and any future monitoring programme must detect trends of increasing ablation against this background of natural variations.

We calculate the present ice budget for Antarctica from measurements of accumulation minus iceberg calving, run-off and in situ melting beneath the floating ice shelves. The resulting negative mass balance of 469 Gt year−1 differs substantially from other recent estimates but some components are subject to high temporal variability and budget uncertainties of 20–50%. Annual accumulation from an earlier review is adjusted to include the Antarctic Peninsula for a total of 2144 Gt year−1. An iceberg production rate of 2016 Gt year−1 is obtained from the volume of large icebergs calculated from satellite images since 1978, and from the results of an international iceberg census project. Ice-shelf melting of 544 Gt year−1 is derived from physical and geochemical observations of meltwater outflow, glaciological field studies and modeling of the sub-ice ocean circulation. The highest melt rates occur near ice fronts and deep within sub-ice cavities. Run-off from the ice-sheet surface and from beneath the grounded ice is taken to be 53 Gt year−1. Less than half of the negative mass balance need come from the grounded ice to account for the unattributed 0.45 mm year−1 in the IPCC “best estimate” of the recent global sea-level rise.

The thermal effects of air flow forced through a snow sample are investigated numerically. A new method for calculating vapor transport in snow is presented which allows for the determination of the effects of sublimation. In this method, the snow is not assumed to be saturated with water vapor. Results of the model show very good agreement with analytical and experimental results. The effects of heat conduction, heat advected by the dry air and heat associated with sublimation are examined in a comprehensive theory, and then each effect is isolated to determine its overall contribution. It is demonstrated that the heat transfer associated with vapor transport is significant in the determination of the overall temperature profile of a ventilated snow sample but that the major effects are controlled by the heat carried by the dry air flow through the snow and heat conduction due to the temperatures imposed at the boundaries. The thermal effects of ventilation of snow are more likely to be observed when there is a smaller temperature gradient over the entire snowpack and a greater flow rate of air in the snow than would be observed when there is a greater overall temperature gradient and lesser air-flow rate.

The subglacial drainage system of South Cascade Glacier is inferred by examining how discharges of water, solutes and artificial tracers in streams draining the glacier respond to meltwater and precipitation on the glacier surface. Results indicate that the glacier is divided into three (two large and one small) drainage basins, each drained by a single stream. The positions of drainage divides on the glacier surface do not always lie directly over those at the base of the glacier. Comparisons of water discharge and cation load in each of the two main streams indicate that subglacial hydraulic processes differ between drainage basins. One stream drains from a conduit that is isolated in its lower reach from the surrounding subglacial region and receives water routed englacially from the surface. The upper reach of the conduit also receives water routed englacially from the surface as well as from a distributed subglacial flow system. The other main stream drains from a conduit coupled to a debris layer beneath the glacier. Observations of the layer in natural ice tunnels indicate that the water may flow within a thin layer of debris. A one-dimensional model of flow through the debris layer can explain both the base-flow and diurnal variations of the second main stream.