Radar-sounding data taken at Uvérsbreen glacier, Svalbard, are used to estimate the relative water content of the temperate part of the glacier. The relative change in water content of the temperate part of the glacier. The relative change in water content is estimated from the strength of the backscttered radar energy without knowing the size of the scatterers. The main findings are that the water content is on average changing by a factor of five in the ice, with the exception of some bright spots believed to stem fro larger water moulins. The general trend is also that the water content is increasing with lower altitude. The applied technique may be a useful method for large-scale mapping of the relative water-content variation in temperate and polythermal glaciers.

Models for calculating glacier mass balance are sensitive to surface reflectivity variation. Fieldwork carried out on Peyto Glacier, Alberta, Canada, contributes to the data set available for ice-reflectivity parameterization in such models. Hemispherical reflectivity in the visible and near-infrared parts of the solar spectrum was obtained for rock, snow and three contrasting glacier surfaces to examine temporal and spatial variations. Glacier-ice near-infrared reflectivity displays only minor spatial variation (0.12–0.17) in comparison with the visible range (0.23–0.40), the latter being influenced primarily by surface impurity content. Surface roughness is of minor importance compared with impurities. Temporal variation of reflectivity was weak at all glacier-ice and rock locations; slight variations observed were due to changes in either solar zenith angle or cloud amount. Snow reflectivity displayed pronounced diurnal asymmetry and a larger response to cloud cover. The minimal temporal variation in glacier-ice reflectivity simplifies its parameterization. This behaviour is additionally useful for satellite-based measurements of the reflectivity field on larger glaciers, as images obtained within a 6 h window centred on solar noon are likely to yield values which are within 2–3% of daily mean values.

A simple geometric model of glacier volume is derived. The model is based on the assumption that the net mass balance averaged over the glacier surface is related to the summer temperature and winter accumulation at a representative height on the glacier. This height varies with time as climate changes, in ways that are determined by the geometry of the glacier surface. Expressions are derived for the equilibrium glacier volume (Veq) as a function of summer temperature and winter accumulation and for the glacier-response time (τ) as a function of volume.

The model is used to reconstruct the volume of Storglaciären over the period AD 500–1992. Measured net mass-balance data for Storglaciren (1946–92) are used to estimate the model parameters. For the summer temperature forcing, the long tree-ring-based temperature reconstruction for northern Fennoscandia is used to extend a temperature record near the glacier back to AD 500. For the past accumulation forcing, a range of assumptions is tested.

The results show that the prime cause of the decrease in volume of Storglaciären between 1946 and 1992 was relatively low accumulation over 1946–89 as opposed to warm temperatures as previously supposed. Reconstructed volume changes agree well with geomorphological evidence. Where differences occur, deductions can be made about past accumulation. For example, the maximum reconstructed glacier volume is associated with the cold period from 1580 to 1740 but relatively low accumulation probably limited the volume to values not much larger than that achieved around 1916.

A relatively lightweight and simple airborne system for surface elevation profiling of glaciers in narrow mountain valleys has been developed and tested. The aircraft position is determined by kinematic global positioning system (GPS) methods. The distance to the glacier surface is determined with a laser ranger. The accuracy is about 0.3 m, sufficient to permit future changes to be observed over short time intervals. Long-term changes can be estimated by comparison of profiles with existing maps. Elevation profiles obtained in 1993–94 from three glaciers in central and south-central Alaska are compared with maps made about 1950. The resulting area-averaged, seasonally corrected thickness changes during the interval are: Gulkana Glacier (central Alaska Range)–11 m, Worthington Glacier (central Chugach Mountains) +7 m, and Bear Lake Glacier (Kenai Mountains) −12 m. All three glaciers retreated during the interval of comparison. The estimated uncertainty in the average thickness change is ±5 m. which is mainly due to errors in the existing maps. Constraints on the accuracy of the maps are obtained by profiling in proglacial areas.

The geographical distribution of surge-type glaciers worldwide displays a remarkably non-random pattern. Surge-type glaciers tend to be concentrated in certain glacierized areas and to be completely absent in others. This observation suggests that special conditions are required for surges to occur. However, the factors controlling the spatial occurrence of surge-type behaviour are not known. To investigate this problem we performed probability statistical analysis on a sample population of 615 glaciers in Svalbard. The probability that a glacier in the sample population is surge-type is 36.4%. Within the sampled area there is a spatial variation in the concentration of surge-type glaciers. Several geometric and environmental factors associated with glaciers in the sample population were measured and tested to determine if they are related to the probability of surging. Of the geometric factors tested (length, slope, elevation, orientation and presence or absence of tributaries), only glacier length is related to surging, with surge probability increasing with increasing length. Elevated probabilities of surging were also found for glaciers associated with sedimentary subglacial rocks and sub-polar thermal regimes. The distributions of related factors were used to predict the spatial distribution of surge-type glaciers. However, in each case the individual factors were unable to reproduce the observed pattern of surge-type glacier distribution.

White Glacier is a valley glacier at 79.5°N with an area of 38.7 km2. Its mass balance has been measured, over 32 years with a 3 year gap, by standard techniques using the stratigraphic system with a stake density of the order of one stake per km2. Errors in stake mass balance are about ±(200–250) mm, due largely to the local unrepresentativeness of measurements. Errors in the whole-glacier mass balance B are of the same order as single-slake errors. However, the lag-1 autocorrelation in the time series of B is effectively zero, so it consists of independent random samples, and the error in the long-term “balance normal” 〈B〉 is noticeably less. 〈B〉 is −100 ± 48 mm. The equilibrium-line altitude (ELA) averages 970 m. with a range of 470–1400 m. Mass balance is well correlated with ELA, but detailed modelling shows that the equilibrium line is undetectable on visible-band satellite images. A reduced network of a few stakes could give acceptable but less accurate estimates of the mass balance, as could estimates based on data from a weather station 120 km away. There is no evidence of a trend in the mass balance of White Glacier. To detect a climatologically plausible trend will require a ten-fold reduction of measurement error, a conclusion which may well apply to most estimates of mass balance based on similar stake densities.

A system for automatic, reliable, semi-instantaneous estimation of avalanche activity is presented in an attempt to check deterministic models and improve the surveillance of risk zones. The principle is the seismic detection of avalanches. Two experimental sites equipped with standard seismological equipment are described. Avalanche seismic signals are recorded, as well as many extraneous signals of natural or human origin. After several years of a posteriori identification of the signals, we are now able to set up criteria for the automatic recognition of non-avalanche signals. We have therefore developed an automatic analysis/decision system to discriminate between avalanche signals and others. This system worked satisfactorily in pre-operational conditions during the winter of 1995-96. The results are presented and compared with other data related to avalanche activity. Although there is still room for improvement, our system seems to be able to estimate avalanche activity better than human visual observations.

Experiments performed on ice single crystals oriented for basal slip indicate that the steady-state creep rate is only marginally affected by confining pressure up to 19 MPa, at a constant absolute temperature of 263 K. The observations Contradict earlier work at similar pressures and the disparity is examined in terms of experimental errors.

Satellite radar interferometry provides glaciologists with an important new tool for determining the motion and topography of large ice sheets. We examine the sources of error in interferometrically derived ice-motion measurements, including those errors due to inaccurate estimates of the interfero-metric baseline. Several simulations are used to assess baseline accuracy in terms of tie-point error and the number and distribution of tie points. These results give insight into how best to select tie points, and also demonstrate the level of accuracy that can be achieved. Examination of two representative cases likely to occur in mapping ice-sheet motion leads to the conclusion that with adequate tie-point information ice velocity can be measured accurately to within a few meters per year. A method to correct horizontal velocity estimates for the effect of vertical displacement using surface slopes is also developed. Finally, we estimate the single-component velocity field for an area on Humboldt Glacier, northern Greenland, using interferograms formed from ERS-1 SAR image. We estimatе that these velocity measurements are accurate to within 2.3 m year1.

Submersible remotely operated vehicles (ROVs) are valuable research tools for data collection in dangerous or inaccessible environments associated with glaciers terminating in the sea. At tide-water ice cliffs, iceberg calving makes close approaches for extended time periods In manned vessels dangerous. ROVs can be operated from relatively safe distances (hundreds of metres); they can also descend to considerably greater depths (hundreds rather than tens of metres) than scuba diving permits. They can provide data on glacier grounding-line and sea-floor morphology and water-column characteristics (e.g. salinity, turbidity, current velocity). They are also used for. diving under floating glacier tongues and ice shelves where no other access is possible. They can be fitted with a variety of oceanographic sensors, imaging sensors, tracking devices and water and sediment samplers, making them versatile research instruments that can supply qualitative and quantitative data for process studies in logistically difficult environments.

The temporal and spatial variation in the surface albedo of the Greenland ice sheet during the ablation season of 1991 is investigated. The study focuses on an area east of Søndre Strømfjord measuring 200 km by 200 km and centred at 67°5′ N, 48° 13′W. The analysis is based on satellite radiance measurements carried out by the Advanced Very High Resolution Radiometer (AVHRR). The broad-band albedo is estimated from the albedos in channel 1 (visible) and channel 2 (near-infrared). The results are calibrated with the surface albedo of sea and dry snow.

Satellite-derived albedos are compared with GIMEX ground measurements at three stations. There is a high degree of consistency in temporal variation at two of the three stations. Large systematic differences are attributed to albedo variations on sub-pixel scale.

In the course of the ablation season four zones appear, each parallel to the ice edge. It is proposed that these are, in order of increasing altitude: (I) clean and dry ice, (II) ice with surface water, (III) superimposed ice, and (IV) snow. An extensive description of these zones is given on the basis of the situation on 25 July 1991. Zones I, III and IV reveal fairly constant albedos (0.46, 0.65 and 0.75 on average), whereas zone II is characterised by an albedo minimum (0.34). Survey of the western margin of the Greenland ice sheet (up to 71° N) shows that the zonation occurs between 66° and 70° N.

Dye-tracer techniques are widely used in infer the character of subglacial drainage systems. Quantitative analysis of dye breakthrough curves focuses on the determination of the water through flow velocity (u), the dispersion coefficient (D) and the dispersivity parameter (d = D/u). Together, these parameters describe the rate of passage of tracer through the drainage system and the extent to which the dye cloud becomes spread out during passage. They have been used to infer the nature of flow conditions within a drainage system and temporal changes in system morphology. Estimates of all three parameters, however, are dependent upon the sampling interval at which measurements of dye concentration used to define breakthrough curves are made. For a given breakthrough curve, the dispersion coefficient increases with the sampling interval, while the through flow velocity shows no systematic variation. As a result, the dispersivity also tends to increase with the sampling interval. Investigations of the sensitivity of parameter estimates to the sampling interval reveal that reliable estimates can be obtained only if the sampling interval is less than 1/16 of the time from dye injection to peak dye concentration. As a general guide, we Suggest that, ideally, quantitative analyses of dye breakthrough curves should therefore be conducted only when this criterion can be met.

We describe a high-precision (0.1–1.0 mK) borehole-temperature (BT) logging system developed at the United States Geological Survey (USGS) for use in remote polar regions. We discuss calibration, operational and data-processing procedures, and present an analysis of the measurement errors. The system is modular to facilitate calibration procedures and field repairs. By interchanging logging cables and temperature sensors, measurements can be made in either shallow air-filled boreholes or liquid-filled holes up to 7 km deep. Data can be acquired in either incremental or continuous-logging modes. The precision of data collected by the new logging system is high enough to detect and quantify various thermal effects at the milli-Kelvin level. To illustrate this capability, we present sample data from the 3 km deep borehole at GISP2, Greenland, and from a 130 m deep air-filled hole at Taylor Dome, Antarctica. The precision of the processed GISP2 continuous temperature logs is 0.25–0.34 mK, while the accuracy is estimated to be 4.5 mK. The effects of fluid convection and the dissipation of the thermal disturbance caused by drilling the borehole are clearly visible in the data. The precision of the incremental Taylor Dome measurements varies from 0.11 to 0.32 mK. depending on the wind strength during the experiments. With this precision, we found that temperature fluctuations and multi-hour trends in the BT measurements correlate well with atmospheric-pressure changes.

Data from Columbia Glacier are used to identify processes that control calving from a temperate tidewater glacier and to re-evaluate models that have been proposed to describe iceberg calving. Since 1981, Columbia Glacier has been retreating rapidly, with an almost seven-fold increase in calving rate from the mid-1970s to 1993. At the same time, the speed of the glacier increased almost as much, so that the actual rate of retreat increased more slowly. According to the commonly accepted model, the calving rate is linearly related to the water depth at the terminus, with retreat of the glacier snout into deeper water, leading to larger calving rates and accelerated retreat. The Columbia Glacier data show that the calving rate is not simply linked to observed quantities such as water depth or stretching rate near the terminus. During the retreat, the thickness at the terminus appears to be linearly correlated with the water depth; at the terminus, the thickness in excess of flotation remained at about 50 m. This suggests that retreat may be initiated when the terminal thickness becomes too small, with the rate of retreat controlled by the rate at which the snout is thinning and by the basal slope. The implication is that the rapid retreat of Columbia Glacier (and other comparable tidewater glaciers) is not the result of an increase in calving as the glacier retreated into deeper water. Instead, the retreat was initiated and maintained by thinning of the glacier. For Columbia Glacier, the continued thinning is probably associated with the increase in glacier speed and retreat may be expected to continue as long as these large speeds are maintained. It is not clear what mechanism may be responsible for the speed-up but the most likely candidate is a change in basal conditions or subglacial drainage. Consequently, the behavior of tidewater glaciers may be controlled by processes acting at the glacier bed rather than by what happens at the glacier terminus.

Large, single ice crystals containing nо air bubbles and free of both small-angle grain boundaries and visible stress birefringence can be grown using a very simple growth chamber within a temperature-controlled, outer enclosure. The method relies upon the spontaneous formation of an ice crystal with its с axis accurately normal to a free, slightly supercooled, water surface.