Annual net balance eight North Cascade glaciers during the 1984-94 period has been determined by measurement of total mass loss firn and ice melt and ice melt and, residual snow depth at the end of the Summer season. Overall spatial density of measurment points is 200 points km−2. Mean annual balance of North Clascade glaciers from 1984 to 1994 has been −0.38 ma−1. The resulting 4.2 m loss in water-equivalent thickness is significant, since North Cascade glaciers have an average thickness of 30–50 m.

Cross-correlation of annual net balance Ior eight glaciers ranges from 0.83 to 0.97. This indicates the mass balances of the eight glaciers have been responding similarly to elimate conditions despite their range of topographic and geographic characteristics. Annual net balance of individual glaciers was correlated with climate records. The highest ablation-season correlation coefficient is mean May–August temperature, ranging from 0.63 to 0.84. The highest accumulation-season correlation coefficient is total accumulation-season precipitation, ranging from 0.35 to 0.59.

Detailed digital elevation models (DEMs) do not exist for much of the Greenland and Antartic ice sheets. Radar altimetry is at present the primary, in many cases the only, source of topographic data over the ice sheets, but the horizontal resolution of such data is coarse. Satellite-radar interferometry uses the phase difference between pairs of synthetic aperture radar (SAR) images to measure both ice-sheet topography and surface displacement. We have applied this technique using ERS-1 SAR data to make detailed (i.e. 80 m horizontal resolution) maps of surface topography in a 100 km by 300 km strip in West Greenland, extending northward from just above Jakobshavns Isbræ. Comparison with а 76 km long line of airborne laser-altimeter data shows that We have achieved a relative accuracy of 2.5 m along the profile. These observations provide a detailed view of dynamically Supported topography near the margin of an ice sheet. In the final section We compare our estimate of topography with phase contours due to motion, and confirm our earlier analysis concerning vertical ice-sheet motion and complexity in ERS-1 SAR interferograms.

Radio-echo soundings provide an effective tool for mapping the thermal regimes of polythermal glaciers on a regional scale. Radar signals of 320–370 MHz penetrate ice at sub-freezing temperatures but are reflected from the top of layers of ice which are at the melting point and contain water. Radar signals of 5–20 MHz, on the other hand, see through both the cold and the temperate ice down to the glacier bed. Radio-echo soundings at these frequencies have been used to investigate the thermal regimes of four polythermal glaciers in Svalbard: Kongsvegen, Uvérsbreen, Midre Lovénbreen and Austre Brøggerbreen. In the ablation area of Kongsvegen, a cold surface layer (50–160 m thick) was underlain by a warm basal layer which is advected from the temperate accumulation area. The surface ablation of this cold layer may be compensated by freezing at its lower cold-temperate interface. This requires that the free water content in the ice at the freezing interface is about 1 % of the volume. The cold surface layer is thicker beneath medial moraines and where cold-based hanging glaciers enter the main ice stream. On Uvérsbreen the thermal regime was similar to that of Kongsvegen. A temperate hole was found in the otherwise cold surface layer of the ablation area in a surface depression between Kongsvegen and Uvérsbreen where meltwater accumulates during the summer (near the subglacial lake Setevatnet, 250 m a.s.l.). Lovénbreen w as frozen to the bed at the snout and along all the mountain slopes but beneath the central part of the glacier a warm basal layer (up to 50 m thick) was fed by temperate ice from two cirques. On Austre Brøggerbreen, a temperate basal layer was not detected by radio-echo soundings but the basal ice was observed to be at the melting point in two boreholes.

The effects on very low-frequency surface-impedence measurements of lateral variations commonly found in ice environments have been measured and modelled numerically using die quasi-static two-dimensional boundary-element method. Results indicate that surface-impedance measurements made in the vicinity of crevasses oriented perpendicular to the plane Of incidence, and those made in the vicinity of moraines and melt streams, can all show significant changes to the measured apparent resistivity. It is, therefore, misleading to use such measurements in the interpretation of ice depth.

A finite-difference numerical model is used to simulate the temperature profile at the center of an ice sheet throughout the course of a glaciation. The ice sheet is gradually built to a thickness of 3000 m over about 10 000 years, starting on permafrost. A geothermal heat flux is applied at large depth. For an initial surface temperature of –12.5°C, our model shows that basal melting occurs 72000 years after the onset of the glaciation. The important parameters determining the basal temperatures are the initial temperature of the ice and substrate, the rate of downward advection of cold ice and, to a lesser extent, the thickness of the ice sheet. The growth history of the ice sheet does not significantly influence the time at which basal melting occurs. Our results show the possibility that the central parts of the continental ice sheets were cold-based for a significant part of their existence. Heating due to the geothermal heat flux cannot account for basal melting during most or all of a glacial cycle. These results may help to explain the existence of preserved land forms under the ice sheets.

A theory of erosion, transport and deposition of unlithified sediments by glaciers is presented. It predicts the large-scale areal distribution of zones and rates of erosion and deposition in time and space through a complete glacial cycle, together with the resultant intensity of large-scale lineations (drumlins) which will be incised in the landscape. The theory also predicts the dispersal patterns of subglacial lithologies, together with the form of dispersal trains derived from distinctive sources and the vertical and horizontal distribution of lithologies within a till. It predicts major erosional discontinuities within tills and the formation of boulder pavements. It suggests that the dominant proportion of the lowland tills produced by Pleistocene mid-latitude ice sheets was generated by subglacial deformation and explains why they are predominantly fine-grained.

The theory is based on an analysis of glacier-dynamic processes and therefore can be used to infer the dynamic behaviour of former ice sheets from the distribution of tills and their lithologic composition.

Previous observations have shown spatial covariances between microwave emission from Antarctic firn at 6 cm wavelength, physical firn temperature and firn-density stratification. Such observations motivate us to understand the physics underlying such covariances and, based on that understanding, to develop estimation methods for firn in which density, and therefore dielectric permittivity, varies randomly in discrete layers with mean thicknesses on the order of centimeters. The model accounts for depth profiles of the physical temperature, mean density and variance of random density fluctuations from layer to layer. We also present a procedure to estimate emission-model input parameters objectively from in situ density-profile observations, as well as uncertainties in the input parameters and corresponding uncertainties in theoretical brightness-temperature predictions. We compare emission-model predictions with ground-based observations at four diverse sites in Antarctica which span a range of accumulation rates and other parameters. We use coincident characterization data to estimate model inputs. At two sites, layered-medium emission-model predictions based on the most probable input parameters (i.e. with no model tuning) agree with observations to within 3.5% for incidence angles≤50°. Corresponding figures for the other two sites are 7.5% and 10%. However, uncertainties in the input parameters are substantial due to the limited length and depth resolution of the characterization data. Uncertainties in brightness-temperature predictions are correspondingly substantial. Thus brightness-temperature predictions for the last-mentioned sites based on only slightly less probable input parameters are also in close agreement with observations. The significance of agreements and discrepancies could be clarified using characterization measurements with finer depth resolution.

This paper describes the development and testing of a distributed surface energy-balance model used to calculate rates of surface melting at Haut Glacier d’Arolla, Valais, Switzerland. The model uses a digital elevation model (DEM) of the glacier surface and surrounding topography together with meterological data collected at a site in front of the glacier to determine hourly or daily totals of the energy-balance components and hence of melting over the entire surface of the glacier with a spatial resolution of 20 m. The model can also be used to determine temporal and spatial variations in snow depth, snow-line position and glacier surface albedo. Calculations from the model are compared with observations made along the glacier centre line 1990, and in general the model performs very well. The correlation coefficients between calculated and measured snow-line elevation, albedo and ablation are 0.99, 0.85 and 0.81, respectively. The main source of error between modelled and measured values of these variables is probably inadequacies in the parameterization of albedo used in the model.

Analytical and numerical techniques are used to examine the flow response of a sloped slab of power-law fluid (power n) subjected to basal boundary conditions that vary spatially across the flow direction, as for example near an ice-stream margin with planar basal topography. The primary assumption is that basal shear stress is proportional to the basal speed times a spatially variable slip resistance. The ratio of mean basal speed to the speed originating from shearing through the thickness. denoted as r, gives a measure of how slippery the bed is. The principal conclusion is that a localized disturbance in slip resistance affects the basal stress and speed in a zone spread over a greater width of the flow. In units of ice thickness H, the spatial scale of spreading is proportional to a single dimensionless number Rn ≡ (r/n+ 1)1/n+1 derived from n and r. The consequence for a shear zone above a sharp jump in slip resistance is that the shearing is spread out over a boundary layer with a width proportional to Rn. For an ice stream caused by a band of low slip resistance with a half-width of w H, the margins influence velocity and stress in the central part of the band depending on Rn in comparison to w. Three regimes can be identified, which for n = 3 are quantified as follows: low r defined as R3 < 0.1w, for which the central flow is essentially unaffected by the margins and the driving stress is supported entire by by basal drag; high r defined as R3 > 1w, for which the boundary layers from both sides bridge across the full flow width and the driving stress in the center is supported almost entirely by side drag; intermediate r, for which the driving stress in the center is supported by a combination of basal and side drag. Shear zones that are narrower than predicted on the basis of this theory (≈ R3) would require localized softening of the ice to explain the concentration of deformation at a shorter scale.

Evidence of changing basal and internal ice properties near the grounding line was derived from airborne radio-echo-sounder observations of the ice sheet around the Sør Rondane Mountains, Antarctica. From the trailing figure of the bottom-echo signal, the roughness of the ice bottom near the grounding line was inferred. Results show that the specular components of scattering begin to appear on the ice-shelf side of the grounding line. Furthermore, double-trip echoes were observed with a strong scattering in the shelf area, and their boundary of occurrence was very close to the grounding line. This is evidence of interaction between ice and sea water at the bottom of the ice shelf. We also examined the occurrence of internal layered echoes. In most of the area around the mountains, internal echoes were observed continously, but they were not found at or close to the ice shelf. The boundary between the appearance and disappearance of internal-layer echoes is distinct, and occurs 20–30 km inland from the grounding line. These results suggest that some major change may occur in the internal ice on the inland side of the grounding line.

Black Rapids Glacier, a surge-type glacier in the Alaska Range, most recently surged in 1936–37 and is currently in its quiescent phase. Mass balance, ice velocity and thickness change have been measured at three to ten sites from 1972 to 1994. The annual speed has undergone cyclical fluctuations of as much as 45% about the mean speed. Ice thickness and surface slope did not change enough to cause the speed fluctuations through changes in ice deformation, which indicates that they are being drinven by changes in basal motion. The behavior of Black Rapids Glacier during this quiescent phase is significantly different from that of Variegated Glacier, another well-studied surge-type glacier in Alaska. The present medial-moraine configuration of Black Rapids Glacier indicates that a surge could occur at any time. However, ice velocity data indicate that the next surge may not be imminent. We belive that there is little chance that the next surge will cross and dam the Delta River.

The stability of a tidewater terminus is controlled by glacial dynamics, calving processes and sedimentary processes at the grounding line. An investigation of grounding-line sediment dynamics and morainal-bank sediment budgets in Glacier Bay, Alaska, U.S.A., has yielded data that enable us to determine the debris fluxes of Grand pacific, Margerie and Muir Glaciers. Debris flux ranges from 105 to 106 m3 a−1, one to two orders of magnitude lower than the glacifluvial sediment fluxes (106−107 m3 a−1). Combined, these fluxes represent the highest yields known for glacierized basins. Large debris fluxes reflect the combined effects of rapid glacier flow, driven by the maritime climate of southeast Alaska, and highly erodible bedrock. Englacial-debris distribution is affected by valley width and relief, both of which control the availability of sediment. The number of tributaries controls the distribution and volume of debris in englacial and supraglacial moraines. At the terminus, iceberg-rafting removes up to two orders of magnitude more sediment from the ice-proximal environment than is deposited by melt-out or is dumped during calving events. Rough estimates of the sediment flux by deforming beds suggests that soft-bed deformation may deliver up to an order of magnitude more sediment to the terminus than is released from within the glacier ice.

Reduced amounts of snow in the eight winters from 1986-87 to 1993-94 at Nagaoka, Japan, seem to be due to a winter air-temprature rise. The winter air temprature has shown cyclic varition gradual increase in the past 100years. The linear rate of the temperature rise in the past century was calculated as 1.35°C per 100 years. Both the maximum Snow depth and winter precipitation showed an inversely positive correlation with winter mean air temperature, The square of the statistical correlation coefficient r2 was calculated as 0.321 and 0.107. respectively. Statistically smoothed curves or the maximum snow depth and winter precipitation showed maximum values in 1940, Fluctuations in deviation of the maximum Snow depth showed smaller values than in precipitation. The minimum winter mean air temperature obtained from a 10 year moving average curve was found in 1942, and the deviation fom the climatic mean changed from negative to positive in 1949. The change in sign or the temperature deviation and the increase of the deviation may be attributable to global warming.

The triangular shape of an ice-ridge sail in the Baltic Sea was studied statistically. The dimensions and orientation of individual ice blocks were measured in several ridges. All measurements were carried out from the standpoint of backscattering research. The results confirm that given geometrical properties of the two sides of ice ridges are in general not equal. A slight negative correlation was found between the slope angle and width of the sail. All three orientations and all three size parameters of ice blocks on both sides of the ridges studied can be considered as normally distributed. A log-normal distribution, however, fits slightly better the dimensions of the orientationally best rectangular approximations of the polygonal main facets of flat ice blocks. Three-dimensional modelling of ice ridges is essential, since the total side-facet area visible is typically at least as large as the total main-facet area. Calculated incidence-angle distrubutions show that the broad distributions of orientations and dimensions make different redges appear very similar on radar images.

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.

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.

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.

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.