Two models that estimate glacier ablation from simple meteorological observations are described. Both the energy-balance and bulk-aerodynamic approaches are used to determine the bulk-exchange coefficient for transfer of sensible- and latent-heat energy to the melting ice surface from 17 summer days of measurement on Ivory Glacier, New Zealand. The behaviour of the coefficient is examined over progressively increasing time intervals and as a function of wind speed and gradient, and of atmospheric stability. The models are verified with 36 d of independent data from a different summer on the same glacier. The more accurate short-term estimates of ablation are obtained using the energy-balance approach with empirically but independently derived values of the exchange coefficient.

Calculation procedures are developed and results shown for the exact calculation of extinction and meteorological visual range using the blowing-snow mass in the atmosphere and particle radius. Results of the calculations show: (1) For monochromatic radiation, geometrical optics approximations of the extinction efficiency are found to provide results of only moderate accuracy in calculating the extinction of radiation by a single particle. (2) For broad-band radiation, the geometrical optics approximation is sufficiently accurate for many single-particle measurement instruments and applications, except in the infra-red band where Mie theory should be used. (3) For typical blowing-snow particle-size distributions, the shape parameter of the distribution of particle radii and the mean particle radius are very important in broad-band extinction and visual-range modelling. Estimates of blowing-snow quantities from broad-band extinction measurements or visual range from blowing-snow quantities should address the shape and mean value of the snow-particle radius distribution.

Two models are presented for the formulation of abrasion and basal drag due to rock–rock friction (debris drag) for the case of sparse debris entrained in the basal layers of a temperate glacier resting on a bedrock bed. The first model is formulated in terms of average basal melting rate, va, and the concentration, C, of basal debris fragments which make intermittent bed contact. The second model is formulated in terms of vn the component of ice velocity normal to the bed flowing around rock fragments contacting the bed, and Cc, the concentration of debris actually in contact with the bed. The relationship between the two models is given for the case of a sinusoidal bed. Generalizations are discussed as well as potentially important physical processes which remain to be investigated.

During the decade prior to its 1982–83 surge, Variegated Glacier experienced progressive changes in geometry and velocity. It thickened in the upper 60% and thinned in the lower 40% of its 20 km length. Thickness changes were up to 20%. Annual velocity increased by up to 500%, reaching a maximum of 0.7 m d−1 in the year before surge onset. Amplitude of seasonal variation in velocity increased up to 0.3 m d−1 by 1978, but did not increase markedly after that. The changes in velocity were larger than predicted from changes in deformation rate caused by changes in shear stress and depth. This anomalous velocity was especially large after 1978 in the zone of thickening on the upper glacier. If it is assumed to arise from basal sliding, the inferred pattern of sliding shows qualitative features consistent with a direct effect from basal shear stress and an inverse effect from effective normal stress. A drop in effective normal stress in a zone of decreasing surface slope up-glacier from the largest thickness increase may have been significant in the initiation of surge motion in 1982.

The mass balance, summer balance, winter balance, and equilibrium-line altitude of three Canadian glaciers (Peyto, Place, and Sentinel Glaciers) are compared with the meteorological records of neighbouring stations for the period 1966—84. While Peyto Glacier’s mass balance is almost entirely related to summer temperature, Sentinel Glacier’s mass balance is mostly controlled by winter precipitation. Place Glacier is influenced by both elements. Statistical reconstructions are presented for the three glaciers, using the best regression equations with the meteorological records since 1938.

When plotted on a δD-δ18O diagram, the isotope values of basal ice formed by refreezing of melt water beneath the cold-based margin of a corrie glacier lie on a straight line with a slope of 5.5. This freezing slope differs from a slope of 7.2 for the basal ice layer of the cold-based ice cap, Balgesvarri, which has not undergone refreezing. The latter slope is similar to that of 7.7 for the local precipitation. Co-isotope measurements from the two glaciers support the association of subglacial debris entrainment with refreezing processes.

The main perturbation in the mass balance of Storglaciären during the twentieth century was caused by a sudden 1°C increase in the summer mean temperature around 1910. Later perturbations of the climate have been of minor importance in relation to the mass balance. Annual field surveys suggest that the mass budget on Storglaciären has been in near balance for the last 15 years. Because of this major step-like change, we can establish the validity of theoretical models giving response times for Storglaciären of the order of 50 years. According to these models, Storglaciären could by now have reached a profile in balance with the present climate. To study this problem, the emergence velocity was calculated and compared with the net balance. The result shows that the emergence velocity either balances or exceeds the net balance for the entire tongue except for the lowermost part, where it decreases to about half of the net balance. A slight further recession of the front position would thus be expected with today’s climate.

Calculated balance velocities also suggest that most of the present profile is close to a steady-state profile, if the mean annual sliding velocity is about 50% of the surface velocity. Lower sliding velocities would imply a thickening of the tongue and a thinning of the accumulation area during years of balanced mass budget.

Measurements of snow properties across and down snow slopes have been used to calculate a safety margin — the difference between the basal shear strength and the applied static stress. Areas of basal deficit exist when the applied shear stress exceeds the basal shear strength (the safety margin is negative), and basal areas are pinned when the safety margin is positive. As the size of deficit increases, stresses within the overlying slab also increase, and these may be sufficient to cause an avalanche.

Measurements made on five slopes (four of which had avalanched) were characterized by considerable spatial variability, and the safety margin has been treated as a random function which varies over the slope. Statistical models of Vanmarcke (1977[a], 1983) have been applied to determine the most likely size of deficit required for avalanching (95% confidence). In one case, an avalanche occurred when the length of deficit was only 2.9 m, and in the other cases the length was always less than 7 m. This size of deficit is small compared with the total area of many avalanche slopes which suggests that avalanches initiate from small zones of deficit, and makes it difficult to locate a deficit with just a few tests.

The optimum sampling interval and number of tests required to yield an adequate estimate of the statistical parameters of the safety margin are also discussed.

As part of a larger experiment, detailed albedo measurements were carried out during the austral summer of 1985-86 in the dry-snow zone (1560 m) of Terre Adélie, eastern Antarctica. The following results were found:

• (1) Mean albedo values were high (around 82.6%). On clear days, the albedo showed some dependency on the solar elevation. The dependency was slight for solar elevations above 12° but became larger with low Sun angles.

• (2) The albedo was found to be a function of cloud amount and type, increasing with the amount and thickness of clouds. In white-out conditions, very high albedos were found (>90%).

• (3) The albedo showed a dependency on the type of snow. New snow displayed higher values than older snow, whose crystals had been destroyed by mechanical action.

• (4) A simple model was developed to assess the influence of sastrugi on the albedo. This model could explain the asymmetric diurnal variation about solar noon of the measured albedo above a sastrugi field.

The above four dependencies might explain the considerable discrepancies which can be found in the literature concerning the snow albedo of Antarctica.

The accumulation of rime ice on structures, due to the impact and freezing of small water droplets, has been modelled as a stochastic process. Individual droplets are introduced into the flow field about the structure at a random position. Their trajectories are then calculated to determine the position of impact on the structure, or on previously impacted droplets. By assuming that the droplets maintain their shape on impact, the modelled accretion is gradually built up, one droplet at a time.

In the present paper, attention has been limited to a circular cylinder as the collecting structure, and it has been assumed that the flow field and the ice accumulation are strictly two-dimensional. With these assumptions, the influence of the droplet/cylinder diameter ratio and of the air speed upon the resulting predictions has been investigated. The main feature of interest in the model prediction is the development, near the edges of the accumulation, of discrete structures called “rime feathers”. The mechanism for the growth of these rime feathers is described, and a comparison is made between the characteristics of the predicted structures and of some natural rime feathers grown in an icing wind tunnel.

The 2037 m deep bore hole at Dye 3 in south Greenland was surveyed in 1981, 1983, 1985, and 1986. The directional surveys show the ice flow is planar with a surface velocity of 12.2m/year at an azimuth of 060°, which agrees with surface velocity measured by navigation satellites. Measurements of hole diameter and inclination are highly correlated with dust content in the ice. The temperature measurements show strong convection in the hole fluid with a cell height of about 20 m and an amplitude of 0.1 K. The calculated mean in-situ ice density is 921.3 ± 1.5kg/m3. Due to ice deformation, the lowest 4 m of the hole were not accessible in 1985 and the lowest 180 m were not accessible in 1986.

Observations of the physical properties of the snow cover and underlying young fast ice in Resolute Passage, N.W.T., Canada, were made during the winter of 1982. Detailed measurements of snow density, and ice and snow temperatures, salinities, and brine volumes were made over a period of 46 d, beginning when the ice was 0.4 m thick and about 8 d old. The recorded values are used in a theoretical mixture model to predict the dielectric properties of the snow cover over the microwave frequency range. The results of this analysis are then used to investigate the effects of the snow properties on the radar back-scatter signatures of young sea ice. The results show that back-scatter is a function of the incidence angle and can change significantly over short periods of time during the early evolutionary phase of ice and snow-cover development. This has important consequences for the identification of young ice forms from Synthetic Aperture or Side Looking Airborne Radar images.

Among the various artificial radioactive markers (mainly from atmospheric nuclear tests), the contamination from the Chernobyl accident has been found at five sampling sites on Northern Hemisphere glaciers. Total beta-activity measurements reveal a very high radio-active level. However, due to the short time of its occurrence, the temporal resolution of this event in the snow layers can be generally quite low and positive results require careful sampling. The size of the signal also depends on the trajectory of the contaminant cloud, the amount of precipitation, and the surface conditions during deposition.

Depth hoar in polar firn forms when large temperature gradients act on low-density firn, but high-density firn does not develop into depth hoar. Low densities in firn may be depositional (burial of surface hoar or still-air snowfall) or diagenetic (mass loss to the free atmosphere, typically in autumn); however, diagenesis is sufficiently strong to cause significant mass loss only in the top 50–100 mm of firn. Between about 50–100 mm and 2 m depth, grain growth and densification are accelerated strongly by temperature-gradient effects; from 2 to 5 m, temperature gradients have a small but measurable effect on rates of processes in firn, and below 5 m rates essentially have isothermal values.

Diagenetic depth-hoar layers typically develop in the autumn, are relatively thick, and have smooth bases. Depositional depth-hoar layers may develop at any season, are relatively thin, and have irregular bases. In low-accumulation regions, visual stratification may preserve only an annual signal, but in high-accumulation regions individual storms or other features may be recognizable.

The mass balance is computed along the Ohio State University (OSU) transect near the Arctic Circle and along the Expédition Glaciologique Internationale au Groenland (EGIG) line. Measured surface velocities are compared with velocities calculated from up-glacial accumulation rate, flow-line spreading, ice thickness, and the depth variation in horizontal velocity. The depth variation in velocity is calculated using the constitutive relation for ice, calculated temperatures within the glacier, computed shear and longitudinal stresses, and allowance for impurity content and ice-crystal orientation. The resulting mass balance is +0.6 ± 0.14 m a−1 for the OSU transect and 0 ± 0.07 m a−1 along the EGIG line. The errors arise mainly from uncertainties in measured accumulation rate and flow-line spreading, and perhaps in flow-enhancement factors due to ice anisotropy or impurities. The results for the EGIG line differ from prior estimates mainly because earlier works placed greater emphasis on short-term accumulation rates.

The annual mass of snow moved by avalanches was observed at 45 avalanche paths over a period of 19 years by measuring the volume and density of each individual avalanche. Several other methods to estimate the total annual avalanche mass were applied but they were found less efficient owing to inaccuracy and difficulties in timing.

For the data set, the yield ratio, i.e. the percentage of snow removed annually by avalanches, had a mean value of 11.2% and a 30 year maximum of 30.9%. The yield ratio varied strongly among avalanche paths and yearly, and could not be explained satisfactorily. The exposure to wind, inclination of the avalanche track, control by artillery, and the winter weather can all have an influence on the value of the parameter.

The temporal variability of avalanche yield is erratic but weather patterns frequently associated with high or low yield ratios could be identified. High values are usually associated with large avalanches triggered by major weather events in the late part of the winter.

Creep curves obtained from the polycrystalline ice samples containing low concentrations of antifreeze glycopeptides (AFGPs) do not show the re-accelerating or tertiary creep that is found in pure ice samples. Previous work has shown that AFGPs inhibit grain-boundary migration in ice and the present results are consistent with this, as tertiary creep is associated with nucleation and growth of new ice grains.

A comprehensive multi-mechanism theory of pressure sintering has been applied to the densification of two polar ice sheets. The comparison, which is made using pressure-sintering mechanism maps, indicates that power-law creep is the controlling mechanism between 50% and 98% theoretical density. Lattice diffusion becomes dominant at low porosities. The densification rates predicted by the theory are in good agreement with the data, and suggest that a reasonable estimation of the densification behaviour of a polar ice sheet can be made using the theory, based on information obtained from a relatively shallow core.

The annual balance has been measured for ten North Cascade glaciers in 1983–84, 1984–85, 1985–86, and 1986–87 (1984, 1985, 1986, and 1987). Based on these data, an annual balance prediction method was designed and tested. Comparison of measured versus predicted annual balances indicates an accuracy of ±0.22–0.30 m. The method is based on annual measurement of the accumulation area ratio (AAR), and determination of the perennially constant activity index and area-altitude distribution on each glacier. The accumulation area ratio is determined from aerial and ground photographs at the end of the ablation season. The activity index is identified from observation of the rise of the snow line with time, compared to measured snow depths above the snow line. The AAR-activity index method was used to calculate the annual balance of 47 North Cascade glaciers in 1984, 1985, 1986, and 1987. The mean balance during the 4 year period was —0.33 m.

From the mass-balance records, it is apparent that North Cascade glaciers can be divided into six climatic sensitivity groups. Each glacier type responds differently to specific climatic conditions. The mass-balance variation for glaciers of the same type is small.

Since 1977, warmer, drier climatic conditions have prevailed in the North Cascades, resulting in the retreat of 42 of the 47 glaciers examined.

A flow-line model is presented for calculating the surface profile and the velocity, strain-rate, and stress fields in an ice sheet with given base-elevation profile, ice thickness at the dome (divide), flow-law parameters, mass-balance distribution, and convergence/divergence conditions along the flow line. The model, which is based on a “quasi-similarity” hypothesis as regards the horizontal velocity-depth profiles, accounts for changes along the flow line in the depth distributions of temperature, normal stress deviators, and possible enhanced flow of deep ice of Wisconsin origin. A curvilinear coordinate system is applied with horizontal axes along flow lines and surface-elevation contours, respectively. The flow equations are reduced to two differential equations, one for the surface-elevation profile, and the other for a profile function that determines the depth distributions of velocities and strain-rates. The two equations are coupled through a profile parameter that communicates the influence of velocity-profile changes to the surface-profile equation. It is shown that the variation along the flow line of this parameter should also be considered when deriving flow-law parameters from ice-sheet flow-line data. For a symmetric dome, explicit expressions are derived for the depth distributions of the vertical velocity, strain-rates, and stresses. The strain-rate profiles display an inflection about half-way down the ice sheet, and, in the case of isothermal ice, have surface values 2.2 times their depth-averaged values. The depth distribution of the vertical velocity indicates that a relatively thick layer of almost stagnant ice is present at the ice-sheet base below a dome.