The fiordlands south of Merchants Bay contain an extensive, well-preserved moraine record of a late Foxe advance of local valley glaciers. This has allowed accurate reconstruction of former glacier margins and computation of former equilibrium-line altitudes (ELAs) by a variety of methods. Statistical comparison of three methods (maximum lateral-moraine elevation, median elevation, and accumulation area ratio (AAR)) shows that different techniques can give different results for the same glaciers. Lateral moraines gave estimates that were too low, probably due to post-glacial erosion or to non-deposition. Median elevations and the AAR method produced statistically similar results but only for glaciers of simple geometry. The median-elevation method fails to take into account variations in valley morphology and glaciological parameters, and so is not reliable in all situations. The AAR method is supported by empirical evidence and is the best of the three methods for estimating former ELAs.

Analysis of trend surfaces of present and late Foxe ELAs shows changes in elevation and orientation through time due to changing environmental factors. Present ELAs are strongly influenced by local factors, southerly storm tracks, and warm maritime conditions. Paleo-ELAs do not show this influence, suggesting that Davis Strait may have been ice-covered during the late Foxe stade and that storm tracks were from the north.

Various field experiments have shown that microwave radars can be used to distinguish multi-year from first-year ice, although optimum radar parameters are not yet fully defined.

This paper presents the results from two theoretical models that, using selected physical parameters of sea ice, are able to predict the backscattering from multi-year and first-year ice under cold conditions. The possible ranges of the backscattering coefficient under various conditions (surface roughness, salinity, temperature, density, and air-bubble size) are calculated for multi-year and first-year ice by adjusting the parameters within the reported range of values.

Although the calculations show no specific resonance that would favor any particular frequency or incidence angles, the results confirm the experimental findings that Ku- and X-band frequencies, and incidence angles greater than 30°, are better for distinguishing sea-ice types than lower frequencies.

Deep gravel-pit exposures reveal the distribution and structure of till and underlying sand and gravel in drumlins near Waukesha, Wisconsin. The subglacial sediment is interpreted to have moved laterally into the drumlin sites because the till thickens from the margin to the core of the drumlins, the stone orientation in the till is perpendicular and oblique to ice flow on the drumlin margins, and recumbent isoclinal folds occur in sand on the drumlin margins with axes parallel to the drumlin axes. The resulting accumulations of sediment presented obstacles to ice flow and were streamlined into the minimum-drag drumlin shape by erosion on the margins and by remolding of material in the core of the drumlins. These drumlin nuclei may have formed at spots where there was low effective stress on the bed. The subglacial sediment became mobile as a result of high pore pressure that may have developed as ground water and subglacial melt water were trapped behind a frozen bed at the ice margin. Under certain conditions, however, lateral sediment flow might also have occurred when the sediment was frozen.

A freezing slope, distinct from that of precipitation, is displayed on a δD–δ18O diagram by basal ice in different circumstances. However, if the subglacial reservoir allowed to freeze is mixed in the course of time with an input having a lighter isotopic composition, basal ice cannot be distinguished from glacier ice in terms of slope. Such a situation is encountered at the base of Grubengletscher and is indicated by a computer simulation using the open-system model of Souchez and Jouzel (1984). Suggested implications for the paleoclimatic interpretation of deep ice cores recovered from the bottom part of polar ice sheets are given.

Analyses of fluoride, chloride, sodium, sulfate, bromide, nitrate, and iron from a 3 m snow pit on Heard Island, collected at an elevation of 2450, m are used to assess the potential of glaciochemical studies on. Heard Island glaciers. Sources for the chemical species are identified and, in particular, chloride, sodium, and sulfate are found to be useful seasonal indicators. The total record measured is believed to be less than one mass-balance year.

Drift-ice abrasion marks are common along present-day rocky shorelines in cold regions. They include polished surfaces, scratches, striations, small grooves, and minor friction cracks. Most are found on relatively soft rocks like shale, sandstone, limestone, dolomite, and basalt, and occasionally on harder rocks like granite and gneiss. They were made by rock fragments frozen at the base of ice cover or by ice floes pushed onshore by wind or dragged along the bottom by waves, tides, and currents. They are found both in the modern and Pleistocene marine, lacustrine, and fluvial environments. Along the Hudson Bay eastern shoreline, these abrasion marks are superimposed on glacially polished and striated surfaces. Along the St. Lawrence Estuary, they are more common on boulders in the tidal zone. Characteristics of drift-ice abrasion marks are described and their significance is pointed out.

Much of the surface of an ice shelf is covered with series of undulations. These undulations or “pressure rollers” are particularly noticeable in the neighbourhood of ice rises or ice streams. To date, there is no satisfactory theoretical model explaining the formation of these waves. As a contribution to understanding this phenomenon, this paper investigates the stability of ice shelves to perturbations in the background stress and strain-rate distributions.

The perturbation analysis is based on Glen’s creep law and leads to a continuous eigenvalue problem for the wavelength of the disturbance as a function of growth-rate. It is shown that, provided these strain-rates are sufficiently compressive, waves of the type observed can be expected to form. It is shown that lateral extensional strain-rates have a destabilizing effect and pressure rollers are more likely to form when these are present. Comparison of predicted wavelengths is made with available field data.

Dielectric experiments have been undertaken at temperatures between −2° and −70° C in the frequency range 10Hz to 100 kHz on 14 firn and ice samples retrieved from the Antarctic Peninsula. This investigation shows that the dielectric behaviour of polar samples from the Antarctic Peninsula is very similar to that of polar firn and ice from Greenland and from elsewhere in Antarctica. In contrast, temperate samples from the Antarctic Peninsula have relaxation times up to ten times shorter for a given temperature between –20° and –70°C, and have higher values of high-frequency conductivity than those of polar samples. Consequently, the thermal regime (temperate or polar) can be distinguished by the dielectric behaviour of the samples.

High-frequency conductivities of polycrystalline samples from the Antarctic Peninsula match the trends of published conductivity data for HF- and HCl-doped laboratory ice; higher conductivities are associated with coastal sites where greater concentrations of marine ions occur in snow.

Annealing polar firn above −10°C results in elevated conductivities across all frequencies measured and shortened relaxation times. Thus, samples for dielectric analysis should not be warmed to above –10°C for risk of irreversibly altering their dielectric behaviour.

Laboratory data on supercooling and frazil-ice formation in sea-water are analysed using a boundary-layer model. The model is based on a turbulent channel-flow boundary-layer theory, in which buoyancy effects become important because of vertical gradients in temperature, salinity, and suspended frazil-ice crystals. The frazil-ice crystals are treated as thin uniform plates. By assuming a mean face diameter, a mean thickness, and a mean Nusselt number of 10−3 m, 10−4 m, and 4, respectively, the general experimental findings are well reproduced by the model.

Jökulhlaups, also known as outburst floods, have occurred every 1 to 5 years from Strandline Lake, one of the largest glacier-dammed lakes in North America. The development of a distinct calving embayment in the lobe of Triumvirate Glacier which dams the lake, as well as the filling of a number of supraglacier pools, appear to be reliable precursors to a jökulhlaup. Analysis of contour maps made from photographs taken immediately before and after the jökulhlaup of 17 September 1982 indicate that over 95% of the lake drains, a volume of approximately 7 × 108 m3 of water. The glacier lobe which dams the lake fractures and subsides during a jökulhlaup, indicating that the release mechanism is hydrostatic lifting of the ice off of a sub-glacial spillway. Evidence from the ice-free margins of the glacier suggests that the spillway may be controlled by bedrock. Large variation in the refilling period of Strandline Lake indicates that the subglacial drainage tunnels can remain open for as much as a few years after a jökulhlaup, before they become sealed by sediments and/or glacier ice.

An evolutionary mathematical model of ice sheets is presented. The model takes into account the basic climatic and geophysical parameters, with temperature parameterization. Some numerical data derived from experiments on the Greenland ice sheet are received. At present the Greenland ice sheet is found to be in a state essentially different from a stationary one corresponding to modern climatic conditions.

In the standard solution for the stresses in a glacier or ice sheet obeying Glen’s law, the down-slope component of the weight is supported by the basal shear stress, and the longitudinal deviatoric stress is second order. However, it has been found necessary to account for the longitudinal stress gradient when relating surface to bed topography with empirical data. In addition, during rapid stretching of the glacier, perhaps during a surge, the longitudinal stress gradient becomes comparable to or larger than the shear stress, and the standard solution is not entirely valid. In this paper, we consider the analysis of the stresses and strain-rates in a glacier when the longitudinal deviatoric stresses are comparable to the basal shear stresses. In some circumstances the down-slope component of weight is not borne completely by basal shear stress to leading order and some of the weight is shifted to the longitudinal deviatoric stress gradient. This case has also been examined. The results are used to obtain expressions for basal shear stress in terms of glacier thickness, slope, surface strain-rate gradient, and ice properties.

The Haut Glacier d’Arolla displays three medial moraines: an ablation-dominant moraine east of the glacier centre line, an ice-stream interaction moraine west of the glacier centre line, and a medial-lateral supraglacial moraine complex along the western margin of the glacier. The eastern moraine is formed by the emergence of coarse angular clasts from a series of short longitudinal debris bands 1 km from the glacier terminus; this debris is derived from rock falls at the glacier head-wall, and follows a low-level transport path through the glacier. The western moraine is formed partly by the direct fall of weathered debris on to the ice surface at a point 2.6 km from the snout, but it receives further increments of sediment from extensive longitudinal englacial debris bands developed from the convergence of two ice streams with lateral englacial debris septa. The medial-lateral moraine complex is formed largely by rock falls from the western valley wall, but is also nourished by a debris layer developed at the base of the small tributary Glacier de la Mitre. This layer crops out at the eastern margin of the re- sultant inset ice stream, and releases more rounded, abraded and comminuted rock fragments to supplement the surface debris cover of supraglacial origin.

The transfer of basal velocity anomalies to the surface of a glacier is investigated using a model of a planar parallel-sided slab (thickness H) of linear viscous rheology. Surface velocity parallel (us) and normal (vs) to the surface is calculated for various spatial distributions of basal velocity anomalies with components parallel (ub) and normal (vb) to the surface. Four scales of differing behavior can be identified depending on the spatial length L of the basal anomalies. At very short scales (L ≤ 1H) there is essentially no response at the surface. At short scales (1H ≤ L ≤ 5H), a basal anomaly ub induces a response in both us and vs. The spatial pattern of us is such that velocity peaks in us can be shifted from peaks in ub, and may differ in number. The amplitude of us is up to about 0.3|ub|. The amplitude of the cross-component effect vs may be greater than the amplitude of us. A basal anomaly vb induces a response in both vs and us. The pattern of vs is the same as the pattern of vb, and the amplitude of vs is up to about 0.7 |vb|. The amplitude of the cross-component effect us is less than the amplitude of vs. At intermediate scales (5H ≤ L ≤ 10H), results differ from the short scale in two respects; velocity peaks in us correspond with peaks in ub; and surface amplitudes are increased, except for cross-component effects for which surface amplitudes are of the same order as at the short scale. These cross-component effects at the short and intermediate scales show in particular that substantial anomalous surface-normal motions can be induced by deformation, even though the basal velocity anomaly is parallel to the surface. At long scales (10H ≤ L), the velocity anomaly at the surface is essentially the same as the anomaly at the bed. For all scales, the longitudinal strain-rate averaged over depth is larger in magnitude than the longitudinal strain-rate at the surface and, at the short scale, it may differ in sign, so that vs cannot be easily estimated from surface strain-rate. Although the simplifications of the model do not allow its rigorous quantitative application to field measurements, the results indicate the need for caution in interpreting surface-velocity variations in terms of basal velocity anomalies. It is important to establish the spatial pattern of surface motions for any chance of a confident interpretation in terms of basal motions.

Eight stakes situated in the ablation and the accumulation areas of a sub-polar glacier in West Greenland were surveyed at intervals of 10 days during the summers of 1982 and 1983, The horizontal velocity in both the ablation and the accumulation areas increased distinctly during the short summer season. This indicates that melt water reaches the bed and that the glacier is sliding. It is proposed that melt water produced in the ablation area is forced up-glacier through a subglacial water system. The vertical displacement of stakes showed variations indicating an apparent uplift of the glacier during the summer. However, this is interpreted as the result of seasonal variation in atmospheric refraction.

A study intended to be mainly of practical use in the prediction of ice-avalanche hazards was carried out. About 100 ice avalanches, mostly from the Alps, were documented. Starting zones of these ice avalanches could be classified by using simple terrain characteristics. Ice avalanches from some starting zones at relatively low altitudes and with large, homogeneously inclined bedrock planes occur predominantly in summer and autumn. No such seasonal variation in frequency was found in connection with other types of starting zones occurring either at higher altitudes or involving abrupt changes of the bedrock’s gradient. One- and two-parameter models for the estimation of run-out distances of landslides and snow avalanches were tested for their use with ice avalanches. Introduction of the second (v2 proportional) frictional parameter leads only to moderately improved accuracy in the prediction of run-out distances. For relatively short run-out distances (several hundred meters), an alternative method of estimation, based on terrain characteristics, is proposed.

A comprehensive set of ice-velocity and thickness data from traverses within the IAGP study area (bounded by long. 90°E. and 135°E., and north of lat. 80°S.) is compared with steady-state mass-flux calculations based on Scott Polar Research Institute (SPRI) map compilations.

The results of previous regional mass-budget estimates are reviewed and followed by a description of the new field measurements and the basis upon which a computer “grid–point” program is used to calculate balance fluxes.

A comparison of measured and balance fluxes indicates that the ice sheet in this region of East Antarctica is unlikely to be significantly out of balance.

The ratio of average column to surface velocity is discussed and calculated to be 0.89.

An analysis of the mean shear strain-rate (VS/Z), versus down-slope basal shear stress (τb = ρgᾱZ), suggests that power flow-law parameters of n = 3.21 and k = 0.023 bar−n m−1 are appropriate for the effective basal shear zone in this region of the Antarctic ice sheet.

Studies on cresceatic fractures and crescentic gouges were carried out in the central part of the Alpine piedmont glaciation and the Quaternary Alpine mountain ice sheet. The test sites were situated in the pro-glacial area of the Hornkees south Berliner Hütte (Oberer Zemmgrund, Zillertal, Tyrol, Austria). This is also the area of the former confluence with another glacier which flowed at about 90° to Hornkees during the late-glacial period. It was therefore possible to make a comparison between different older glacier systems, different flow directions, and different ice thicknesses. About 12 000 minor features were measured by close-range photogrammetry; the result was 240 pairs of stereophotographs which were interpreted on first-order equipment. The single forms were registered point by point by coordinates. It was therefore possible to obtain three-dimensionally registered forms. One important result was that the late-glacial and the post-glacial glaciers of the Zemmgrund produced minor features with a large variation in size. A further result was the definition of possible glacier-flow directions by interpretation of the different directions of crescentic fractures and crescentic gouges. By this means it was possible to obtain several working reconstructions of the late- and post-glacial stream lines of the merged glaciers.

The two-dimensional, time-dependent flow of an arbitrarily shaped ice mass can be successfully modeled with the finite-element technique on a small computer. Methods developed for automatically generating the mesh data greatly simplify the data preparation and optimize the numerical simulations. Using quadratic basis functions permits the flow to be approximated quite adequately by only two element rows (five nodes vertically). Mixed-order basis functions, however, must be used so that numerical oscillations do not set in, and the ends of the ice mass, where the thickness tends to zero, must be treated carefully. Time simulations to a steady-state condition are necessary to test such numerical models adequately.

South Cascade Glacier, Washington, is currently close to equilibrium. A bedrock sill dominates the bed topography in the lower half of the glacier, rising to a height of about 20% of the ice thickness. This sill produces a maximum increase in the overall thickness of about 6–7% compared to what the thickness would have been if the sill were not present. Finally, this glacier does not appear to be sliding much, if at all, despite its maritime alpine environment. This could help explain the difficulties encountered when trying to measure sliding and basal water pressures on the same glacier (Hodge, 1979), or it could imply that drag exerted by the valley walls has a significantly greater effect than conventional shape-factor concepts imply.

During the summer of 1983 a 2 kW array of solar voltaic panels was used to power a drill which collected two cores of 154 m and 163 m length to bedrock on the Quelccaya ice cap in Peru. The operation and performance of this array are discussed, considering the need for future requirements in remote pristine areas of the world.