High-resolution radar images of Glaciar Moreno, Argentina, acquired by the Shuttle Imaging Radar C (SIR-C) on 9 and 10 October 1994 at 24 cm wavelength (L-band), are utilized to map the glacier velocity both interferometrically and using the phase correlation method. The precision of the interferometric ice velocities is 1.8 cm d-1 (6 m a-1) (1σ). The phase correlation method measures ice velocity with a precision of 14 cm d-1 (50 m a-1) with image data at a 6 m sample spacing acquired 1 day apart. Averaged strain rates are measured with a precision of 10-4 d-1 at a 240 m sample spacing with the phase correlation method, and 10-5 d-l with radar interferometry. The phase correlation method is less precise than radar interferometry, but it performs better in areas of rapid flow, is more robust to temporal changes in glacier scattering and measures the glacier velocity in two dimensions with only one image pair. Using this technique, we find that Glaciar Moreno flows at 400 m a-1 in the terminal valley and 800 m a-1 at the calving front, in agreement with velocities recorded a decade ago. Assuming steady-state flow conditions, the vertical strain rates measured by SIR-C are combined with prior data on mass ablation to estimate the glacier thickness and ice discharge. The calculated discharge is 0.6 ± 0.2 km3 ice a-1 at 300 m elevation, and 1.1 ± 0.2 km3 ice a-1 at the equilibrium-line elevation (1150 m), which yields a balance accumulation of 6 ± 1 m ice a-1.

We present results from satellite imagery, ice-motion surveys and ice-penetrating radar studies of part of the north margin of Ice Stream C, one of the ice streams draining the West Antarctic ice sheet to the “Ross Embayment”. Our studies suggest that the shutdown of Ice Stream C about 150 years ago was not a single event, but a sequence involving stagnation of ice and migration of the ice-stream boundary. Ground-based studies confirm the inference from imagery that a series of former shear zones exist, decreasing in age towards the ice-stream center. A region of ice-stream trunk, including a former margin, lies sheared and folded between the (recent) inner and (older) outer margins of the area. ice-motion and topographic surveys give some constraint on the time of shutdown of the outer margin. The results provide a forum for discussing shutdown mechanisms. Possible causes for the stepwise migration of the north margin of Ice Stream C include a gradual decrease in ice flux, a reduction in the available water or hydrostatic pressure in the basal till, or a freezing of the till layers on the northern side.

Laboratory creep deformation experiments have been conducted on initially isotropic laboratory-made samples of polycrystalline ice. Steady-state tertiary creep rates, , were determined at strains exceeding 10% in either uniaxial-compression or simple-shear experiments. Isotropic minimum strain rates, , determined at ˜1 % strain, provide a reference for comparing the relative magnitude of tertiary creep rates in shear and compression through the use of strain-rate enhancement factors, E, defined as the ratio of corresponding tertiary and isotropic minimum creep rates, i.e. . The magnitude of strain-rate enhancement in simple shear was found to exceed that in uniaxial compression by a constant factor of 2.3. Results of experiments conducted at octahedral shear stresses of to = 0.040.80 MPa indicate a creep power-law stress exponent of n = 3 for isotropic minimum creep rates and n = 3.5 for tertiary creep rates. The difference in stress exponents for minimum and tertiary creep regimes can be interpreted as a t 0 stress-dependent level of strain-rate enhancement, i.e. .The implications of these results for deformation in complex multicomponent stress configurations and at stresses below those used in the current experiments are discussed.

The depth of reflecting layers in Arctic ice sheets has been determined by electromagnetic echo sounding, using a varying distance between transmitter and receiver to determine the radar wave velocity. The depth of the radar reflecting layers is compared with a profile of electrical conductivity measurements (ECMs) from the Greenland Ice Core Project (GRIP) ice core, in order to determine the velocity of the radar waves in the ice cap. By using several reflecting layers, it is possible to isolate the firn correction of the wave velocity and to estimate the accuracy of the calculated electromagnetic wave velocity. The measured firn correction is compared with the correction calculated from the density profile, and a comparison between the depth profiles of ECM and radar based on the corrected electromagnetic wave velocity is presented. This profile shows that acid layers, which originate from major volcanic eruptions, show up as reflecting radar horizons.

Snow accumulation is an important component of the mass balance of alpine glaciers. To improve our understanding of the processes related to accumulation and their representation in state-of-the-art mass-balance models, extensive field measurements are required. We present measurements of snow accumulation distribution on Findelengletscher, Switzerland, for April 2010 using (1) in situ snow probings, (2) airborne ground-penetrating radar (GPR) and (3) differencing of two airborne light detection and ranging (lidar) digital elevation models (DEMs). Calculating high-resolution snow depth from DEM-differencing requires careful correction for vertical ice-flow velocity and densification in the accumulation area. All three methods reveal a general increase in snow depth with elevation, but also a significant small-scale spatial variability. Lidar-differencing and in situ snow probings show good agreement for the mean specific winter balance (0.72 and 0.78 m w.e., respectively). The lidar-derived distributed snow depth reveals significant zonal correlations with elevation, slope and curvature in a multiple linear regression model. Unlike lidar-differencing, GPR-derived snow depth is not affected by glacier dynamics or firn compaction, but to a smaller degree by snow density and liquid water content. It is thus a valuable independent data source for validation. The simultaneous availability of the three datasets facilitates the comparison of the methods and contributes to a better understanding of processes that govern winter accumulation distribution on alpine glaciers.

We present a detailed, complete glacier inventory for Alaska and neighboring Canada using multi-sensor satellite data from 2000 to 2011. For each glacier, we derive outlines and 51 variables, including center-line lengths, outline types and debris cover. We find 86 723 km2 of glacier area (27 109 glaciers >0.025 km2), ∼12% of the global glacierized area outside ice sheets. Of this area 12.0% is drained by 39 marine-terminating glaciers (74 km of tidewater margin), and 19.3% by 148 lake- and river-terminating glaciers (420 km of lake-/river margin). The overall debris cover is 11%, with considerable differences among regions, ranging from 1.4% in the Kenai Mountains to 28% in the Central Alaska Range. Comparison of outlines from different sources on >2500 km2 of glacierized area yields a total area difference of ∼10%, emphasizing the difficulties in accurately delineating debris-covered glaciers. Assuming fully correlated (systematic) errors, uncertainties in area reach 6% for all Alaska glaciers, but further analysis is needed to explore adequate error correlation scales. Preliminary analysis of the glacier database yields a new set of well-constrained area/length scaling parameters and shows good agreement between our area–altitude distributions and previously established synthetic hypsometries. The new glacier database will be valuable to further explore relations between glacier variables and glacier behavior.

We combine European Remote-sensing Satellite (ERS-1) radar altimeter surface elevations (Fricker and others, 2000) with six different accumulation distributions to compute balance fluxes for the Lambert Glacier–Amery Ice Shelf drainage system. These interpolated balance fluxes are compared with fluxes derived from in situ measurements of ice thickness and velocity at 73 stations of the Lambert Glacier basin traverse and at 11 stations further downstream, to assess the system’s state of balance. For the upstream line we obtain a range of imbalance estimates, from −23.8% to +19.9% of the observed flux, reflecting the sensitivity to the accumulation distributions. For some of the accumulation distributions the imbalance estimates vary significantly between different parts of the line. Imbalance estimates for the downstream line range from −17.7% to +70.2%, with four of the estimates exceeding +30%, again reflecting the sensitivity of the result to input accumulation, and strongly suggesting that the mass balance of the region between the two lines is positive. Our results confirm the importance of accurate estimates of accumulation in ice-sheet mass-balance studies. Furthermore, they suggest that it is not possible to accurately determine the state of balance of large Antarctic drainage basins on the basis of currently available accumulation distributions.

In summer 2005, two pilot snow/firn cores were obtained at 5365 and 5206 m a.s.l. on Fedchenko glacier, Pamirs, Tajikistan, the world’s longest and deepest alpine glacier. The well-defined seasonal layering appearing in stable-isotope and trace element distribution identified the physical links controlling the climate and aerosol concentration signals. Air temperature and humidity/precipitation were the primary determinants of stable-isotope ratios. Most precipitation over the Pamirs originated in the Atlantic. In summer, water vapor was re-evaporated from semi-arid regions in central Eurasia. The semi-arid regions contribute to non-soluble aerosol loading in snow accumulated on Fedchenko glacier. In the Pamir core, concentrations of rare earth elements, major and other elements were less than those in the Tien Shan but greater than those in Antarctica, Greenland, the Alps and the Altai. The content of heavy metals in the Fedchenko cores is 2–14 times lower than in the Altai glaciers. Loess from Afghan–Tajik deposits is the predominant lithogenic material transported to the Pamirs. Trace elements generally showed that aerosol concentration tended to increase on the windward slopes during dust storms but tended to decrease with altitude under clear conditions. The trace element profile documented one of the most severe droughts in the 20th century.

Measurements of changes in surface height carried out with sonic altimeters mounted on automatic weather stations in Dronning Maud Land (DML) and on Berkner Island, Antarctica, are used to derive the surface mass balance. The surface mass balance is positive at all sites, i.e. accumulation outweighs ablation. The spatial and temporal variability in accumulation is high. Accumulation occurs in numerous small events and a few large events per year. The larger events contribute more to the annual accumulation than the small events; ∼50% of all accumulation is contributed by 10–25% of all events. The accumulation generally decreases with increasing distance from the coast and elevation. Annual averaged values range from ∼375 ± 59 mm w.e. a−1 near the coast to ∼33 ± 25 mm w.e. a−1 on theAntarctic plateau and are in good agreement with long-term averaged annual accumulation rates obtained from snow pits and firn cores. The records show seasonal dependency of the amount of accumulation, with a maximum in winter in the coastal and escarpment region of DML and in summer on Berkner Island and on the plateau. The seasonal cycles are significant on Berkner Island, and in the coastal area and part of the escarpment region.

A simple equation is derived relating the net mass-balance and hypsometric curves of a steady-state valley glacier. It is used to examine how valley shape is linked to disparate extents and responses of glaciers subjected to similar climatic conditions. Examples are given which show that area-based indices (e.g. AAR) for estimating the equilibrium line altitude (ELA) may be subject to a substantial built-in variance because they implicitly rely upon similarity of glacier shape and regimen over a region. If accurate topographic maps are available, the equation may be used to infer the regimen of modern glaciers in the form of a dimensionless ratio of net mass-balance gradients. Alternatively, if similar information is available concerning regional glacier regimen, disparate extents and responses may be collectively utilized to estimate values of ELA or to infer climatic influence, taking glacier hypsometry into account.

The discussion on global change has led to increased interest in glacier mass balance since glaciers can be used as climatic indicators. To meet the need for high-quality mass-balance data requires critical examination of traditional mass-balance methods and their possible errors. One issue regarding mass-balance measurements that has received little attention is internal accumulation. Our study shows that internal accumulation in the firn layer of Storglaciären, Sweden, significantly affects the mass balance of the glacier. This occurs because the winter cold wave penetrates below the previous year’s summer surface and into underlying firn. We estimated internal accumulation from measurements of temperature and water content in firn. The depth of the 0°C isotherm correlated with snow depth and air temperature, so that low snow depth and low air temperature separately cause a deeper 0°C isotherm. We determined irreducible gravimetric water content in firn to 2–3%, which corresponds to an irreducible water saturation of 6–8%. Our value for firn is relatively high compared with that for snow, probably due to trapped water in isolated firn pores. Refreezing of percolating meltwater in spring accounted for ~30% of annual internal accumulation. The remaining 70% was due to re-freezing of retained capillary water in firn pores during winter. Disregarding internal accumulation would lead to underestimation of annual net mass balance by 0.04–0.06 m w.e., corresponding to 3–5% of annual accumulation of the entire glacier in an average year. Hence, internal accumulation potentially becomes a source for systematic error if not accounted in mass-balance measurements.

Drilling bore holes in deep, cold ice masses by hot-water methods and maintaining these holes with sufficient diameter to allow down-hole experimentation poses a major obstacle to the investigation of conditions beneath ice sheets and ice streams. Closure of the water-filled holes by refreezing is the dominant difficulty. In this paper, we describe calculations of heat transfer from the drilling system to the ice and the subsequent time-dependent motion of the phase boundary defining the bore-hole wall. Results are presented with the view of optimizing the bore-hole radius at depth for a fixed drill performance and a variable rate of drilling.

Calculation of melting/refreezing rates at the bore-hole wall requires the use of a one-dimensional, time-dependent numerical heat-flow model with a distorting mesh which follows the changing hole size. The delay of hole closure is discussed with a view to keeping holes open long enough to allow instruments to be lowered to the glacier bed, while realizing that drilling-system performance may be marginal because of logisitical and/or expenditure constraints. The relative merits of drilling a large hole, which is very time consuming with a small drill, and the use of water-soluble antifreezes, which have a history of creating plugs of ice slush, are discussed. A method of creating a stable hole filled with antifreeze in which ice slush does not occur is described.

The recent application of these theoretical ideas to the planning and implementation of successful hot-water drilling programs in Antarctica and Greenland is also presented.

In 1997 a 121 m ice core was retrieved from Lomonosovfonna, the highest ice field in Spitsbergen, Svalbard (1250 m a.s.l.). Radar measurements indicate an ice depth of 126.5 m, and borehole temperature measurements show that the ice is below the melting point. High-resolution sampling of major ions, oxygen isotopes and deuterium has been performed on the core, and the results from the uppermost 36 m suggest that quasi-annual signals are preserved. The 1963 radioactive layer is situated at 18.5–18.95 m, giving a mean annual accumulation of 0.36 m w.e. for the period 1963–96. The upper 36 m of the ice core was dated back to 1920 by counting layers provided by the seasonal variations of the ions in addition to using a constant accumulation rate, with thinning by pure shear according to Nye (1963). The stratigraphy does not seem to have been obliterated by meltwater percolation, in contrast to most previous core sites on Svalbard. The anthropogenic influence on the Svalbard environment is illustrated by increased levels of sulphate, nitrate and acidity. Both nitrate and sulphate levels started to increase in the late 1940s, remained high until the late 1980s and have decreased during the last 15 years. The records of δ18O, MSA (methane-sulphonic acid), and melt features along the core agree with the temperature record from Longyearbyen and the sea-ice record from the Barents Sea at a multi-year resolution, suggesting that this ice core reflects local climatic conditions.

Modelling the response of Himalayan glaciers to rapid climate change is an important problem. The poorly understood effects of debris cover and the variable response of the glaciers have made it difficult to understand their dynamics. We propose a simple model for debris-covered glaciers and validate it against data from Dokriani Glacier, India. Numerical investigations of the model show that the response of debris-covered glaciers to a warming climate has two timescales. There is a period when the glacier loses ice by thinning but the front is almost stationary and it develops a long, slow-flowing tongue. This stationary period, which can be >100 years for glaciers with a large extent of debris cover, is negligible for bare glaciers. The quasi-stagnant tongue does not develop in response to cooling. An analysis of remote-sensing data in the light of these results indicates that the variable response of the glaciers in the Himalaya is consistent with a climate that is warming on average, but has considerable spatial variability in the warming rates. We estimate the average warming rate to be about the same as the global average.

Both lake-calving Yakutat Glacier (337 km2), Alaska, USA, and its parent icefield (810 km2) are experiencing strong thinning, and under current climate conditions will eventually disappear. Comparison of digital elevation models shows that Yakutat Glacier thinned at area-averaged rates of 4.76 ± 0.06 m w.e.a−1 (2000–07) and 3.66 ± 0.03 m w.e.a−1 (2007–10). Simultaneously, adjacent Yakutat Icefield land-terminating glaciers thinned at lower but still substantial rates (3.79 and 2.94 m w.e.a−1 respectively for the same time periods), indicating lake-calving dynamics helps drive increased mass loss. Yakutat Glacier terminates into Harlequin Lake and for over a decade sustained a ∼3 km long floating tongue, which started to disintegrate into large tabular icebergs in 2010. Such floating tongues are rarely seen on temperate tidewater glaciers. We hypothesize that this difference is likely due to the lack of submarine melting in the case of lake-calving glaciers. Floating-tongue ice losses were evaluated in terms of overall mass balance and contribution to sea-level rise. The post-Little Ice Age collapse of Yakutat Icefield was driven in part by tidewater calving retreats of adjacent glaciers, the lake-calving retreat of Yakutat Glacier, a warming climate and by the positive feedback mechanisms through surface lowering.

This work presents a method of mapping deformation-related sublimation patterns, formed on the surface of ice specimens, at microscopic resolution (3–4 μm pixel−1). The method is based on the systematic sublimation of a microtomed piece of ice, prepared either as a thick or a thin section. The mapping system consists of an optical microscope, a CCD video camera and a computer-controlled xy-stage. About 1500 images are needed to build a high-resolution mosaic map of a 4.5 × 9 cm section. Mosaics and single images are used to derive a variety of statistical data about air inclusions (air bubbles and air clathrate hydrates), texture (grain size, shape and orientation) and deformation-related features (subgrain boundaries, slip bands, subgrain islands and loops, pinned and bulged grain boundaries). The most common sublimation patterns are described, and their relevance for the deformation of polar ice is briefly discussed.

Spatial and temporal variations of the flow of Hansbreen, a tidewater glacier in southern Spitsbergen, Svalbard, are investigated. During summer 1999, surface flow velocities were measured in the ablation zone of Hansbreen with a temporal resolution of 1—2 hours. Short events with strongly increased surface velocities and a typical duration of 1—2 days were observed. These “speed-up events” are related to periods of strongly increased water input to the glacier, due to rainfall or enhanced surface melt. A close relation is found between the surface velocities and water pressure recorded in a moulin. However, there are indications from a short time lag between velocity and water-pressure peak as well as from observed vertical surface uplifts that basal motion is related to basal water storage rather than directly to basal water pressure. The observed short-term velocity variations and associated processes on Hansbreen are very similar to those observed on land-based valley glaciers and suggest that the relevant mechanisms and physical processes that control the flow and its temporal variations are similar. In contrast to the flow of land-based glaciers, sliding velocities on Hansbreen are observed to be high all year round and velocities increase towards the calving front.

Despite global warming and unlike their Himalayan neighbours, glaciers in the Karakoram mountains do not show signs of significant retreat. Here we report high velocity variations of Biafo Glacier, central Karakoram, which occurred between 2001 and 2009 and which indicate considerable dynamics in its flow behaviour. We have generated a dense time series of glacier surface velocities, based on cross-correlation of optical satellite images, which clearly shows seasonal and interannual velocity variations, reaching 50% in some places. The interannual velocity variations resemble the passing of a broad wave of high velocities, with peak velocities during 2005 and some diffusion down-glacier over a period of at least 4 years. High interannual velocity variations are also observed at other glaciers in the vicinity, suggesting a common cause, although these appear to partly comprise longer acceleration phases. Analysis of weather station data provides some indications of meteorological conditions that could have promoted sustained sliding events during this period, but this does not explain the wave-like nature of the acceleration at Biafo Glacier, and the regular, protracted velocity changes.

Broadband radiation schemes (parameterizations) are commonly used tools in glacier mass-balance modelling, but their performance at high altitude in the tropics has not been evaluated in detail. Here we take advantage of a high-quality 2 year record of global radiation (G ) and incoming longwave radiation (L ↓) measured on Kersten Glacier, Kilimanjaro, East Africa, at 5873 m a.s.l., to optimize parameterizations of G and L ↓. We show that the two radiation terms can be related by an effective cloud-cover fraction neff , so G or L ↓ can be modelled based on neff derived from measured L ↓ or G, respectively. At neff = 1, G is reduced to 35% of clear-sky G, and L ↓ increases by 45–65% (depending on altitude) relative to clear-sky L ↓. Validation for a 1 year dataset of G and L ↓ obtained at 4850 m on Glaciar Artesonraju, Peruvian Andes, yields a satisfactory performance of the radiation scheme. Whether this performance is acceptable for mass-balance studies of tropical glaciers is explored by applying the data from Glaciar Artesonraju to a physically based mass-balance model, which requires, among others, G and L ↓ as forcing variables. Uncertainties in modelled mass balance introduced by the radiation parameterizations do not exceed those that can be caused by errors in the radiation measurements. Hence, this paper provides a tool for inclusion in spatially distributed mass-balance modelling of tropical glaciers and/or extension of radiation data when only G or L ↓ is measured.

Snow stratigraphy and water percolation are key contributing factors to avalanche formation. So far, only destructive methods can provide this kind of information. Radar technology allows continuous, non-destructive scanning of the snowpack so that the temporal evolution of internal properties can be followed. We installed an upward-looking ground-penetrating radar system (upGPR) at the Weissfluhjoch study site (Davos, Switzerland). During two winter seasons (2010/11 and 2011/12) we recorded data with the aim of quantitatively determining snowpack properties and their temporal evolution. We automatically derived the snow height with an accuracy of about ± 5 cm, tracked the settlement of internal layers (± 7 cm) and measured the amount of new snow (± 10 cm). Using external snow height measurements, we determined the bulk density with a mean error of 4.3% compared to manual measurements. Radar-derived snow water equivalent deviated from manual measurements by 5%. Furthermore, we tracked the location of the dry-to-wet transition in the snowpack until water percolated to the ground. Based on the transition and an independent snow height measurement it was possible to estimate the volumetric liquid water content and its temporal evolution. Even though we need additional information to derive some of the snow properties, our results show that it is possible to quantitatively derive snow properties with upGPR.