The spatial distribution of surface mass balance on the Greenland ice sheet is mapped on a 50 km grid using a combination of methods depending on a zonal characterization of the diagenetic snow fades. In the zones of dry snow and upper percolation fades, the accumulation rate is calculated from microwave emissivities derived from satellite measurements using a model that is calibrated with field accumulation data. In the lower percolation zone, accumulation rates are obtained from visual interpolation of previously compiled field data, with some modification so the balance is zero at the equilibrium line In the ablation zone, ablation rates are calculated as a function of ice-surface elevation and latitude. Average values of the surface balance are 263 kg m–2a–1 in the accumulation zone, (-) 1259 kg m–2 a–1 in the ablation zone and 1286 kg m–2 a–1 overall. Compared to the findings of a previous study using practically the same approach but different models, our bulk estimate of balance (216 Gt a–1) is 57% smaller, but the differences in the estimates of net accumulation and net ablation are, respectively, 30% and 172% larger. In this and other comparisons, there is evidence that the differences in estimates are primarily due to differences in the delineation of the equilibrium line and the estimate of ablation, and secondarily to the estimate of accumulation and interpolation of field data. The differences noted with six other estimates reported in the last two decades are all of a size close to the composite variation of the difference (±50Gtσ–1) . Our surface balance is smaller than three estimates, larger than one and in agreement with two. If substituted in the latest mass-budget estimate that indicates equilibrium, our surface balance estimate would suggest a negative budget of 55 Gt a–1 and thus a positive contribution to sea-level change of 0.15 mm a–1.

Retention of meltwater runoff by percolation and/or refreezing in the snowpack cannot be neglected when studying the surface mass balance of the Greenland ice sheet. In this paper, we make a detailed comparison of several treatments proposed in the literature to account for this process in large-scale mass-balance parameterizations. The melt on the Greenland ice sheet is calculated with a revised degree-day model using updated datasets of surface elevation and precipitation rate on a 5 km grid. Crucial model parameters are recalibrated by comparing mass-balance characteristics with available observations on a regional basis. We discuss the role of meltwater retention in the light of the overall mass balance of the Greenland ice sheet and its sensitivity to climatic change, and display patterns of effective-retention fractions for the various methods. As a main conclusion it appears that overall results are quite similar for the various models, but that meltwater retention has a large spatial variation not described by the simple treatments. Using the most comprehensive retention model, the sensitivity of the runoff is found to be +0.35 mm ˚C–1 of sea-level change per year. We also present a new map of the different zones (facies) that characterize the accumulation area of the Greenland ice sheet, which is useful for interpreting field data and calibrating satellite observations.

We have measured ice flow and detailed topography in northeastern Greenland using satellite-based synthetic-aperture radar (SAR) interferometry. The subject of this study is the large ice stream that drains this quadrant of the ice sheet. A single SAR interferogram allows the measurement of one component of motion over a several-day-long interval. We have used a set of such measurements from multiple look directions to produce a mosaic of ice-flow velocity. The resulting flow field is tied to an estimated balance-velocity distribution in slow-moving areas and assumes flow to be locally surface-parallel. The velocity field is the most detailed, consistent data-set available over a flow feature of this size. It compares with global positioning system surveyed velocity measurements at the 5 m a–1 level In the process of mapping ice-flow velocity an enhanced elevation model of the ice stream was produced. The elevation model is based on a blend of interferometrically measured short-wavelength topography and radar-altimetry-determined longer-wavelength topography. This enhanced model has improved information on local surface slope, which is useful for estimating the horizontal components of the velocity field.

Finsterwalderbreen (77°26’ N, 15° 15’ E) is a 35 km2 polythermal, surge-type glacier in southern Spitsbergen, Svalbard. Extensive field investigations have been carried out on the glacier, including mapping of bed topography by radio-echo soundings, mass-balance measurements and hydrological and meteorological registrations. In this paper we have used surface and bed topography and observations of summer ablation to estimate total runoff and subglacial drainage pattern by means of spatial map analysis in standard geographic information system software. The location and relative importance of the water-channel outlets from the glacier were estimated, as well as the change in drainage pattern over the period 1970–90. The results could be compared to direct observations of drainage pattern and discharge recordings. The analyses indicate that the location of the subglacial channels is dominated by the surface topography and thus that the channels are partly pressurized.

Comparison of the terminus locations of Rongbuk Glacier, Mount Everest, measured in 1966 and 1997 shows that in the past 30 years the glacier has retreated 170–270 m equivalent to a retreat speed of 5.5–8.7 m a–1 . During summer 1997, a 15 m firn core was recovered from Dasuopu glacier (28°23’ N, 85°44’ E; 7000 m a.s.D on the northwest margin of Xixabangma Feng, Xizang (Tibet). The seasonal variations of δ18O values in the core indicate that monsoon signals are clearly recorded in the glacier. δ18O values are controlled by the amount effect in the monsoon season; more negative δ18O is representative of the monsoon season in snow layers. Analysis of the relationship between ice-core δ18O, sampled from 6500 m a.s.l. on the north side of Mount Everest, and instrumental series representing regional-scale precipitation, atmospheric circulation and temperature suggests a change in the relative influence of these parameters on δ18O since the 1940s. The results of the comparison add to and lengthen the sparse array of instrument data available for the Tibetan (Qinghai-Xizang) Plateau and demonstrate a recent decline in moisture flux for at least the southern part of the plateau. Glacier retreat, associated with a recent increase in temperature in the region, is coincident with this period of decreased moisture flux.

Data on the mass balance of Xiao Dongkemadi glacier in the Tanggula mountains, central Tibetan Plateau, were obtained over 5 5 years from 1989 to 1995. These are the first continuous mass-balance data for a continental-type glacier on the Tibetan Plateau, where the glacier accumulates during the summer monsoon (summer-accumulation-type glacier). Mass-balance vs altitude profiles were steeper in the negative than in the positive mass-balance years. This is considered to have resulted from the effect of summer accumulation. The annual mass balance is compared with air temperature, precipitation, and black-body temperature in the area including the glacier, which is calculated from infrared radiation observations by theJapanese Geostationary Meteorological Satellite. It was found that the interannual variation in the glacier mass balance was not closely related to maximum monthly mean air temperature, while it did have a relatively good correlation with maximum monthly mean black-body temperature.

A large number of multispectral and stereo-image data are expected to become available as part of the Global Land Ice Measurements from Space project. We investigate digital elevation model extraction, anisotropic reflectance correction and selected glacier analysis tasks that must be developed to achieve full utility of these new data. Results indicate that glaciers in the Karakoram and Nanga Parbat Himalaya, northern Pakistan, exhibit unique spectral, spatial and geomorphometric patterns that can be exploited by various models and algorithms to produce accurate information regarding glacier extent, supraglacial features and glacier geomorphology The integration of spectral, spatial and geomorphometric features, coupled with approaches for advanced pattern recognition, can help geoscientists study glacier mass balance, glacier erosion, sediment-transfer efficiency and landscape evolution.

An isopleth map showing the spatial distribution of net mass accumulation at the surface on the Antarctic ice sheet, excluding Graham Land, the Larsen Ice Shelf and eastern Palmer Land, is produced based on field data from approximately 2000 sites. A database of accumulation values for 5365 gridpoint locations with 50 km spacing is interpolated from the isopleth map, giving a bulk accumulation of 2151 Gt a–1 and a mean of 159 kg m–2 a–1 for an area of 13.53 × 106 km2. Following the implementation of deflation and ablation adjustments applicable to sectors of the coastal zone, the accumulation values are reduced to 2020 Gt a–1 and 149 kg m–2 a–1. The new accumulation distribution is compared with another recent distribution, which was based on essentially the same field data using different analysis and interpolation criteria. Differences between the distributions are assessed using residuals for the 50 km gridpoint locations and by comparing average accumulation values for 24 drainage systems. The assessment based on residuals indicates that the two distributions show patterns of accumulation that are coherent at the continental scale, a shared attribute underscored by a small mean residual value of 6 kg m–2 a–1 (a difference of <4%). However, the regional assessment based on average accumulation values for the drainage systems shows differences that are larger than the assessment error (>22%) for six systems that collectively comprise approximately 4/10 of the ice-sheet area and 3/10 of the accumulation.

The characteristics of summer energy budgets in the ablation zones during the summer at the surface of two glaciers in the Antarctic Peninsula are investigated and compared to the findings of previous studies. The study areas are located on King George Island (62° S) and in Marguerite Bay (68° S). The summer energy balance was computed from automatic weather station data. The results reveal that turbulent heat fluxes dominate over radiation balance in Marguerite Bay whereas on King George Island ablation is driven by net radiation. The summer energy balance on King George Island reflects the more maritime subpolar climate along the northwest tip of the peninsula in contrast to a more continentally toned polar climate in Marguerite Bay and areas further south. The terms of the energy balances are partitioned very differently but ad-vection from northerly directions causes the highest summer snowmelt rates at both study sites. It is concluded that sensitivity studies should consider not only the mean air-temperature increase, but also possible changes in other climate parameters.

A network of relatively fast-flowing tributaries in the catchment basins of the West Antarctic ice streams transport ice from the inland reservoir to the heads of the ice streams. Branches of the network follow valleys in basal topography but not all valleys contain tributaries. We investigate the circumstances favoring tributary flow upstream of Ice Streams D and E, using a combination of observation and numerical modelling. No consistent pattern emerges. The transition from tributary to ice-stream flow occurs smoothly along the main tributary feeding into the onset of Ice Stream D, with ice thickness being relatively more important upstream, and sliding being relatively more important downstream. Elsewhere, the downstream pattern of flow is more complicated, with local increases and decreases in the contribution of sliding to ice speed. Those changes may be due to variations in basal water storage, subglacial geologic properties or a combination of the two.

An approach to deriving the thickness of the Antarctic ice sheet and inferring the bedrock elevation to fill in gaps in the compilations of ice-thickness observations is presented. It combines assumptions about the general state of balance of the ice sheet and of the dynamics of ice flow, with information about ice accumulation and the topography of the ice-sheet surface elevation to infer the ice thickness and bedrock. A simplified version of the scheme already shows the potential of this approach, as is demonstrated by an application to the Lambert Glacier basin, East Antarctica.

The Chinese National Antarctic Research Expedition (GHINARE) carried out three traverses from Zhongshan station to Dome A, Princess Elizabeth Land and Inaccessible Area, East Antarctic ice sheet, during the 1996/97 to 1998/99 Antarctic field seasons. The expeditions are part of the Chinese International Trans-Antarctic Scientific Expedition program. In this project, glaciological investigations of mass balance, ice temperature, ice flow, stratigraphy in snow pits and snow/firn ice cores, as well as the glaciochemical study of surface snow and shallow ice cores, have been carried out. In the 1998/99 field season, CHINARE extended the traverse route to 1128 km inland from Zhongshan station. The density profiles show that firnification over Princess Elizabeth Land and Inaccessible Area (290–1100 km along the route) is fairly slow, and the accumulation rate recovered from snow pits along the initial 460 km of the route is 4.6–21 cm (46–210 kg m–2a–1 ) water equivalent. The initial 460 km of the route can be divided into four sections based on the differences of accumulation rate. This pattern approximately coincides with the study on the Lambert Glacier basin (LGB) by Australian scientists. During the past 50 years, the trends of both air temperature and accumulation rate show a slight increase in this area, in contrast to the west side of the LGB. Data on surface accumulation rates and their spatial and temporal variability over ice-drainage areas such as the LGB are essential for precise mass-balance calculation of the whole ice sheet, and are important for driving ice-sheet models and testing atmospheric models.

The motion field of the northern Larsen Ice Shelf, Antarctic Peninsula, was analyzed using radar interferometry in combination with field measurements and finite-element model calculations. The ice shelf between Jason Peninsula and Seal Nunataks has been in steady retreat since January 1995. Model calculations suggest that the ice shelf is in a stage of irreversible retreat since the last calving events in summer 1998/99. The interferometric analysis is based on synthetic-aperture radar data of the tandem mission of European remote-sensing satellites ERS-1 and –2 in austral spring 1995. The phase contributions due to tidal motion were estimated from the vertical displacement at those parts of the grounding zone where the horizontal motion is close to zero" in order to separate the vertical and horizontal motion components over the ice shelf. Satellite-derived velocities compare well with the long-term field measurements along a transverse and a longitudinal profile. The real interferograms and synthetic interferograms, calculated from model velocities, show reasonable agreement over the main parts of the ice shelf, but differ in the boundary zones where the details are not resolved by the model.

A high-resolution time-dependent two-dimensional flowline model was developed, which takes into account all relevant stresses (i.e. both shear and longitudinal stresses), and which is solved on a fixed finite-difference grid. Special attention was paid to the computational efficiency. A comparison is made with a similar model that solves the velocity field according to the shallow-ice approximation. Both models were applied to an East Antarctic flowline for which the data on surface and bedrock elevation were sampled at different spatial resolutions. The numerical analysis deals with the importance of longitudinal stresses at the ice-sheet/ice-shelf interface at different grid sizes" and the importance of basal drag vs driving stress in basal sliding laws. At a grid resolution of 5–10 km the force balance indicates that all stress components are of equal importance in an area that is larger than the nominal grid size. When sliding becomes dominant, this transition zone considerably widens. Furthermore, the longitudinal stress deviator is found to vary considerably with depth due to the thermomechanical coupling.

A numerical glacier-flow model (finite-element method) is used to suggest the processes that control the flow behind the calving front of a tidewater glacier. The model is developed for grounded calving glaciers and includes an effective-pressure-dependent sliding law The sliding law is implemented by adding a soft basal layer with a variable viscosity The model is applied on Hansbreen, a tidewater calving glacier in Svalbard. Comparison between modeled surface velocities and observed velocity data shows good agreement. We conclude that the flow of a grounded calving glacier can be modeled with an effective-pressure-dependent sliding law

Outbursts of subglacial water from numerous alpine glaciers have been observed and documented. Such events tend to occur in spring and are thus attributed to an inability of the winter subglacial drainage system (characterized by high water pressure and low capacity) to accommodate a sudden and profuse influx of surface meltwater. Prior to a release event, bursts of glacier motion are common, and the release then precipitates the restoration of summer plumbing that damps or terminates surface acceleration. The events bear witness to the importance of interactions between surface melt, runoff, en-glacial water storage and internal routing, in addition to subglacial drainage morphology. Using a distributed numerical model to simultaneously solve surficial, englacial and subglacial water-transport equations, we investigate the role of these components in a hydro-mechanical event observed at Trapridge Glacier, YukonTerritory, Canada, in July 1990.

An experimental program at Worthington Glacier, Alaska, U.S.A., has yielded the first in situ measurement of the full stress tensor in glacier ice. Measurements were made with an array of stiff (low-compliance) normal-force sensors frozen into a borehole at 120 m depth. Freezing in temperate ice was accomplished by a down-hole heat exchanger which extracted heat at a rate of 15 W. Under slowly varying stress conditions, relaxation of stress anomalies by viscous creep following drilling of the hole and installation of the sensors allows for equilibration of measured stresses with far-field stresses. Equilibration of local and far-field stresses was confirmed and pressure sensors calibrated in laboratory experiments prior to the field program. Results of the stress measurements show principal axes of the stress tensor oriented in directions consistent with the geometry of the glacier and broadly consistent with measured englacial strain rate. The magnitudes of stress-tensor components are more error-prone and more sensitive to uncertainty in sensor magnitude than uncertainty in sensor orientation. Mean stress determined by pressure measurements agrees with estimated lithostatic overburden to within approximately 15%. Unexpected results include a stress perturbation lasting about 5 days that caused a rotation of the orientations of the principal stress axes of approximately 5° about an axis pointing in the down-flow direction.

Dislocation/grain-boundary (GB) interactions have been studied in situ in polycrystalline ice using synchrotron X-ray topography in the temperature range 0° to –15°C GBs were observed to act both as sources of lattice dislocations and as strong obstacles to dislocation motion. Dislocations were observed to form pile-ups at GBs upon loading. Generally the basal slip system with the highest Schmid factor was found to be the most active, and dislocations were emitted from GB facets as semi-hexagonal loops in order to relieve the stress build-up from GB sliding. When the relative orientation of two adjacent grains and the orientation of the GB between them with respect to the loading direction discouraged GB sliding, thus suppressing dislocation nucleation at the GB, dislocations originating in one grain piled up at the GB and led to slip transmission through the GB The latter geometrical arrangement is rarely encountered, suggesting that slip transmission through grain boundaries in ice is a rare event. When basal slip was suppressed, i.e. when the loading direction lay in the basal plane, slip occurred by the glide of a fast edge segment on non-basal planes.

The fluctuation of mountain glaciers is recognized as a high-confidence indicator of air-temperature trends and as a valuable element of a strategy for early detection of possible Man-induced climate changes. The Terrestrial Observation Panel for Climate therefore recommended that glacier mass and area be monitored as part of the Global Climate Observing System (GCOS) established in 1992 by the World Meteorological Organization, the Intergovernmental Oceanographic Commission, the United Nations Environment Programme and the International Council of Scientific Unions. A tiered Global Hierarchical Observing Strategy was developed to be used for all GCOS terrestrial variables. According to this system of tiers, the regional to global representativeness in space and time of the records relating to glacier mass and area should be assessed by more numerous observations of glacier length changes as well as by compilation of regional glacier inventories repeated at time intervals of a few decades, the typical dynamic response time of smaller mountain glaciers.

During the 1970s, Fritz Müller directed the Permanent Service on the Fluctuations of Glaciers and the Temporary Technical Secretariat for the World Glacier Inventory. These two bodies were combined in 1986 to form the World Glacier Monitoring Service, which is now responsible for internationally coordinated glacier monitoring, working in close collaboration with the World Data Center for Glaciology, Boulder.

In order to investigate the temperature dependence of vibrational spectra of ice, we measured the Raman spectra of artificial ice Ih in the temperature range 198–270 K. The frequency of translational lattice vibrations decreases with increases in the measurement temperature, and the rate of decrease discontinuously changes at 237 K. The discontinuous change of the rate of decrease is consistent with the phase transition from a proton-disordered arrangement to a proton-ordered arrangement at the ice temperature Ti = 237 K in polar ice sheets, as proposed by Fukazawa and others (1998b). We then compared the Raman spectra of translational lattice vibrations in ice Ih and Antarctic ice and examined the possible effects of the phase transition on the geophysical properties of Antarctic ice. We also report preliminary results of the neutron powder diffraction of Dome Fuji (DF) Antarctic ice, and discuss arrangements of protons in DF ice.