Thermal convection in the Antarctic and Greenland ice sheets has been dismissed on the grounds that radio-echo stratigraphy is undisturbed for long distances. However, the undisturbed stratigraphy lies, for the most part, above the density inversion in polar ice sheets and therefore does not disprove convection. An echo-free zone is widespread below the density inversion, yet nobody has cited this as a strong indication that convection is indeed present at depth. A generalized Rayleigh criterion for thermal convection in elastic-viscoplastic polycrystalline solids heated from below is developed and applied to ice-sheet convection. An infinite Rayleigh number at the onset of primary creep decreases with time and becomes constant when secondary creep dominates, suggesting that any thermal buoyancy stress can initiate convection but convection cannot be sustained below a buoyancy stress of about 3 kPa. An analysis of the temperature profile down the Byrd Station core hole suggests that about 1000 m of ice below the density inversion will sustain convection. Creep along the Byrd Station strain network, radar sounding in East Antarctica, and seismic sounding in West Antarctica are examined for evidence of convective creep superimposed on advective creep. It is concluded that the evidence for convection is there, if we look for it with the intention of finding it.

A number of measurements of ice flow in the coastal regions of Antarctica are given. Observations show that the general outward movement of the continental ice, termed “sheet flow”, is locally accelerated where “ice streams” are formed. Estimates indicate that ice streams, which only occupy a small fraction of the total length of coast, are responsible for the removal of more ice from the continent than the “sheet flow” over the remaining length of coast. Further estimates suggest that the great bulk of Antarctic icebergs are produced by ice shelves, but that data on ice shelf movement are inadequate at the time of writing.

A method is described for measuring the thermal conductivity of ice in the temperature range from 0.5 to 5.0 K using a 3He apparatus. The results from our first experiments are not too far from the theoretical law for the low-temperature thermal conductivity of ice . Measurements at still lower temperatures are necessary to confirm our results.

Formulae for calculating temperature profiles in steady-state ice sheets are derived in terms of the horizontal and vertical velocities of the surface ice, and of the variation in surface temperature along the flow lines. Two models are considered: parallel flow and radial flow. In each case the result consists of a “stationary” term (equivalent to Robin's equation), a term dependent on the horizontal velocity, and a term arising from the temperature variations. This last term produces reversals in the temperature field such as are measured in practice.

Stable-oxygen-isotope ratios (δ18O) of precipitation and ice-core samples collected from the headwaters of the Urümqi river, Tien Shan, China, were used to test the relationship between δ18O and contemporaneous surface air temperature (Ta). A strong temporal relationship is found between δ18O in precipitation and Ta, particularly for the monthly averages which remove synoptic-scale influences such as changes in condensation level, condensation temperature and moisture sources ( Yao and others, 1996). Linear fits as high as 0.95‰° C-1 for precipitation events and 1.23‰° C-1 for monthly averages are found. Although the δ18O amplitude in ice cores drilled at the nearby Ürümqi glacier No.1(~2 km from the precipitation sampling site) decreased dramatically compared to the precipitation samples, the ice-core records of annually averaged δ18 are still positively correlated with contemporaneous air temperature, especially summer air temperature, at the nearby Daxigou meteorological station. Nevertheless, the relationship between the ice-core δ18O records and contemporaneous air temperature is less significant than that for the precipitation samples due to depositional and post-depositional modification processes, which are highlighted by the successive snow-pit δ18O profiles from the Ürümqi glacier No.1. Our results might extend the application of high-altitude and subtropical ice-core δ18O —Ta records for paleoclimate reconstruction.

An analysis is made of the rate of bubble coalescence in an ice mass that is deforming. A total strain of at least 8 is required before appreciable coalescence occurs. The analysis has been applied to deforming ice shelves and ice sheets. No appreciable coalescence is expected in ice shelves but coalescence should occur in ice sheets (or glaciers) if the shear strain-rate at the bottom surface is of the order of 0·075/year or larger. Measurements of bubble concentration are capable of setting limits on paleo-strain-rates of the present ice sheets. Bubble migration down temperature gradients presents complications to the study of bubble coalescence.

The design and operating characteristics are presented for a visible and near-infrared scanning photometer which can measure incident and reflected spectral irradiance from 400 nm to 2 450 nm. The instrument is designed to record over 98% of the solar radiation incident upon and reflected from the Earth’s surface using a turret-type cosine collector and a circular variable interference filter. The apparatus has been tested successfully on Arctic sea ice at temperatures down to –20°C. It is self-contained and easily portable, and its mass is less than 16 kg. Some preliminary results are presented.

The movement of the Irish Sea glacier is discussed for the St George’s Channel sector during the last major glaciation. It is considered that the ice was already diverging when it crossed Pembrokeshire, after passing through a constriction farther north. This constriction was probably caused by local Irish ice in the west and local Welsh ice in the east. In view of this, it is suggested that the Cardigan Bay coasts could not have been ice-free during the last glaciation.

The ice movement on the Byrd Glacier, Antarctica, was measured as function of time, on the basis of aerial photography. The method of aero-levelling strip triangulation with a minimum of geodetic control was applied. The methods and procedures were verified over Laufen-Bauma test area.

A study of the movement of Rusty Glacier was undertaken and continued through four summers because it is believed to be a surging glacier in the last stages of the inactive phase preceding a surge. The entire glacier is very slow moving, essentially motionless in the lower third and most rapid in an area well above the firn line. Unusually steep flow-line emergence angles and higher than average longitudinal compression rates in the lower–middle part of the glacier indicate gradual thickening of the ice above the stagnant lower tongue. There is no clear correlation between local variations in flow rates and surface or bottom topography. The glacier is mostly colder than 0° C to the bottom, and in the one known area of 0° C bottom temperature, flow rates are not greater than elsewhere. Although the glacier is everywhere very thin, maximum flow rates seem clearly related only to variations in ice thickness. The nearby Trapridge Glacier is also a surging glacier and exhibits a strikingly similar flow pattern.

In the summer of 1972, the amount of ice melt obtained by stake measurements on the McCall Glacier, Brooks Range, was compared with the discharge of a small controlled run-off site, and with heat balance calculations. For daily values the agreement was not always good; however, for the 11 d period for which the study was carried out, the agreement can be considered as excellent. The amount of ice melt calculated from run-off and heat balance gave 3.5 and 1.2% higher values, respectively, than the direct measurements. This is encouraging, insofar as it suggests that it may be possible to calculate the combined ice, water, and heat balance for the glacier as a whole.

The results of 66 Schlumberger geoelectric soundings on first-year sea ice are analysed to obtain average values of d.c. electrical resistivity both parallel ρ1 and perpendicular ρp to the long brine cell direction. This analysis yields ρ1 = 104 Ω m, 80 Ω m, 45 Ω m, and ρp = 298 Ω m, 230 Ω m, 129 Ω m for brine volumes νb of 40‰, 67‰, and 82‰ respectively. The observed dependence of ρ1 on νb is used in conjunction with a mixing formula for the electrical resistivity of a two-phase heterogeneous system to arrive at values of “average” brine cell lengths l for relatively warm sea ice. This treatment yields l ≈ 1.7 cm for 40‰ ≤νb≤ 70‰ with a rapid extension of brine cell lengths for brine volumes above c. 70‰.

Glacial ripping involves glaciotectonic disintegration of rock hills and extensive removal of rock at the ice-sheet bed, triggered by hydraulic jacking caused by fluctuating water pressures. Evidence from eastern Sweden shows that glacial ripping caused significant subglacial erosion during the final deglaciation of the Fennoscandian ice sheet, distinct from abrasion and plucking (quarrying). Here we analyse the ice drag forces exerted onto rock obstacles at the base of an ice sheet, and the resisting forces of such rock obstacles: glaciotectonic disintegration requires that ice drag forces exceed the resisting forces of the rock obstacle. We consider rock obstacles of different sizes, shapes and fracture patterns, informed by natural examples from eastern Sweden. Our analysis shows that limited overpressure events, unfavourable fracture patterns, low-transmissivity fractures, slow ice and streamlined rock hamper rock hill disintegration. Conversely, under fast ice flow and fluctuating water pressures, disintegration is possible if the rock hill contains subhorizontal, transmissive fractures. Rock steps on previously smooth, abraded surfaces, caused by hydraulic jacking, also enhance drag forces and can cause disintegration of a rock hill. Glacial ripping is a physically plausible erosion mechanism, under realistic glaciological conditions prevalent near ice margins.

A two-dimensional finite element computer program has been used to compute the elastic stress distribution in realistic multi-layered snow packs. Computations have been done on three-layered and five-layered snow packs intended to simulate conditions on the Lift Gully at Berthoud Pass, Colorado. Calculations have been performed to determine the effect of a layer of new snow and the effect of a weak sub-layer. Stress levels were obtained which are reasonable compared with available snow strength data.

Inclinometry with embedded probes is analyzed with a Stokes model of a solid body floating in a fluid with much smaller viscosity for a two-dimensional flow field. The assumption that such a probe behaves like a Lagrangian unit vector is only justified for probes embedded in a Newtonian fluid with lengths at least four times their width. A fluid with Glen-type rheology results in a slightly smaller rotation rate of the probe compared to Newtonian fluids.

Since its inception in 1975, the Northwest Avalanche Center (NWAC) has developed and produced micro- and mesoscale weather forecasts to support avalanche forecast and control needs for the Olympic and Cascades Mountains of Washington and Oregon, U.S.A. This paper describes NWAC’s array of data, observational results, and analytical techniques that make “avalanche weather” forecasting possible. In addition, NWAC’s operational program and the general terrain and climate of the area are described. A sample forecast is also included.