Optical satellite imagery was used to estimate glacier surface velocities and iceberg calving rates from Agassiz and western Grant Ice Caps, Nunavut, Canada, between 1999 and 2003. The largest mean annual surface velocities ranged from ∼400 to 700 m a−1, but velocities in the ∼100–200 m a−1 range were common. Summer velocities were up to an order of magnitude larger than annually averaged velocities. The highest velocity (∼1530 m a−1) was measured on the floating tongue of Lake Tuborg Glacier between 19 July and 19 August 2001. Calving rates from individual glaciers varied by up to a factor of two between successive years. Summer calving rates were ∼2–8 times larger than annual average rates. The average ratio of the calving flux due to terminus-volume change to that due to ice flow through the glacier terminus was ∼0.81 for the annual rates and ∼1.71 for summer rates. The estimated mean annual calving rate from Agassiz Ice Cap in the period 1999–2002 was 0.67 ± 0.15 km3 a−1, of which ∼54% emanated from Eugenie Glacier alone. This total rate is similar to a previously estimated calving rate from Devon Ice Cap.

Satellite imagery has been used to acquire seasonal terminus positions of tidewater Hubbard Glacier, Alaska, USA, from 1992 to 2006. During this 15 year time period, the width-averaged advance of the entire terminus has been ∼620 m at a mean rate of 35 m a−1. Seasonal fluctuation of the terminus ranges from 150 to 200 m on average and varies spatially. A section of the terminus, near a narrow gap where the glacier has now twice closed off 40 km long Russell Fiord, exhibited little to no mean advance during this time period but displayed seasonal fluctuations of 300–500 m. Seasonal variability in surface ice speeds and surface sea-water temperatures was also observed; both are potential forcing mechanisms for terminus fluctuations. Seasonal changes in sea-water temperature of 10–12°C, as well as seasonal changes in subglacial freshwater discharge, are inferred to influence calving and submarine melting at the terminus, driving seasonal variations. Displacements of the medial moraine separating Hubbard and Valerie Glaciers at the terminus suggest surge-like pulses of the latter, with a periodicity of several years. The timing of these pulses suggests they may influence the Hubbard terminus near Gilbert Point and have implications for future closures of Russell Fiord.

We present a re-analysis of the results obtained from a series of measurements on freshwater and saline ice beams under various centrifugal accelerations. The data show a strong influence of beam size, brine volume and centrifugal acceleration on the elastic modulus of ice. The data suggest a transition brine volume at around 9%, which might occur close to the melting point, at which the elastic modulus of ice drops rapidly due to a possible change of brine-pocket structure. Furthermore, for brine volumes less than 9%, there is a negligible increase in the elastic modulus measured under high centrifugal acceleration, but for brine volumes more than 9% the increase is considerable, approaching that measured with freshwater ice. This may be due to necking of brine drainage channels just above the ice/water interface at high centrifugal acceleration. A model of sea ice was constructed based on existing theories of brine inclusions in sea ice, which satisfactorily predicts the observed trends.

Flow stripes seen in satellite imagery of ice streams and ice shelves are caused by surface undulations with kilometer-scale spacing and meter-scale relief and generally indicate current or recent fast ice flow. On a similar scale, folding of internal ice stratigraphy depicted in cross-flow ice-penetrating radar profiles is also a common occurrence in ice streams, suggesting a possible relationship between the two sets of features. We have traced surface flow stripes in RADARSAT and MODIS imagery on Kamb Ice Stream, West Antarctica, from the onset of streaming flow into the near-stagnant trunk. We compare the morphology and evolution of the surface flow stripes to the folds seen in the internal stratigraphy in cross-ice-stream radar profiles. We find essentially no correspondence in the observed locations or spacings between the radar internal layer folds at depths greater than 100 m and the flow stripes on the surface. The gap in the radar data and the surface mappings in the top 100 m of firn prevents a precise depiction of how the flow stripes and fold patterns at depth diverge. We explore hypotheses about how flow stripes and internal stratigraphic folds can originate and evolve differently as ice flows downstream. We suggest that flow stripes are subject to surface processes that can modify their morphology independently of the internal stratigraphy, leading to changes in the pattern of flow stripes relative to the internal layers below.

Devon Ice Cap, Nunavut, Canada, has been losing mass since at least 1960. Laser-altimetry surveys, however, suggest that the high-elevation region (>1200 m) of the ice cap thickened between 1995 and 2000, perhaps because of anomalously high accumulation rates during this period. We derive an independent estimate of thickness change in this region by comparing ∼40 year mean annual net accumulation rates to mean specific outflow rates for 11 drainage basins. The area-averaged rate of thickness change across the whole region is within error of zero (0.01 ± 0.12 m w.e. a−1), but two drainage basins in the northwest are thickening significantly, and two basins in the south are thinning significantly. The laser-altimetry observations are biased towards the drainage basins where we find thickening. Recent changes in the rate of accumulation or the rate of firnification cannot explain the observed thickening, but decreased ice outflow, due to the penetration of Neoglacial cooling to, and subsequent stiffening of, the basal ice, may provide an explanation. Thinning in the south may result from increased ice outflow from basins in which fast flow and basal sliding extend above 1200 m.

Several authors relate accumulation (or precipitation) at the glacier equilibrium-line altitude (ELA) to summer mean temperature using exponential or power-law functions. I analyze the accumulation–temperature relation at the ELA with a degree-day model using data from the 1992 paper by A. Ohmura and others. The dataset includes estimates at the ELA of winter balance and of ‘winter balance plus summer precipitation’ which represent respectively low and high estimates of annual accumulation, which is seldom measured. The Ohmura dataset only lists summer mean temperature, but I recover monthly temperatures for the whole year for 66 of the glaciers by assuming sinusoidal temperature variation through the year and using annual temperature range from a gridded climatology. Monthly degree-day sums are then estimated from monthly mean temperature and summed to give annual totals so degree-day factors for melting snow at the ELA are obtained. The degree-day factors fall close to those reported in the literature for glacier snowmelt, with averages of 3.5 ± 1.4 and 4.6 ± 1.4 mm d−1 K−1 for low- and high-accumulation estimates on the 66 glaciers. The degree-day model gives a family of accumulation–temperature curves that depend upon the annual temperature range, representing the contrast between maritime and continental climates.

We describe field measurements (ground-penetrating radar (GPR), geodetic survey and ice-core drilling) to provide new information on the movement mechanism and internal structure of a polar rock glacier on James Ross Island, Antarctic Peninsula. We collected GPR data along longitudinal and transverse profiles. The longitudinal GPR profiles identify inter-bedded debris-rich layers that dip up-glacier, similar to the thrust structures in the compression zone of a valley glacier. The transverse GPR profiles indicate a syncline structure inclined towards the central part of the rock glacier, resembling the transverse foliation of a valley glacier. The stratigraphy of two boreholes shows that the rock glacier consists primarily of bubbly ice with thin debris-rich layers, an internal structure similar to the ‘nested spoons’ structure common in the interior of valley glaciers. These results indicate that the glacier motion is controlled by shear movement, common in valley glaciers. The geodetic survey confirms that flow velocities decrease towards the lower part of the rock glacier. Such heterogeneous movement causes longitudinal compression and forms thrusts which then create the debris-rich layer by uplifting basal ice and debris. Pushing of the upstream ice against the downstream ice bends the surface layers, forming transverse ridges on the rock glacier surface.

The first analyses of biological components in glaciers of the Chilean lake district are presented based on microalgae biovolume, pollen and other microorganisms detected in shallow (∼10 m) firn/ice cores. Three cores were retrieved, two at Volcán Mocho-Choshuenco (39°55′ S, 72°02′ W; summit at 2422 m a.s.l; east glacier at 2000 m a.s.l.), and one at the summit of Volcán Osorno (41°06′ S, 72°30′ W; 2652 m a.s.l.). Microalgae, protozoa and pollen quantified in the samples obtained from the two summit cores show clear fluctuations interpreted as seasonal signs. In contrast, δD and many chemical species from the summit cores show strong dampening at depth, probably due to water percolation. The limited information provided by isotopic and chemical analyses is used to support the seasonal interpretation of biological parameters from the summit cores, with microorganism maxima inferred to occur in summer and pollen maxima in spring. A good comparison is found between mass-balance estimations from the Volcán Mocho-Choshuenco summit core and values obtained near that site by means of the stake method. It is concluded that biological analyses of firn/ice cores provide reliable estimations of annual and seasonal markers from these temperate glaciers.

On an 11 m firn/ice core from Kongsvegen, Svalbard, we have used dielectric profiling (DEP) to measure electrical properties, and digital photography to measure a core optical stratigraphy (COS) profile. We also used a neutron-scattering probe (NP) to measure a density profile in the borehole from which the core was extracted. The NP- and DEP-derived density profiles were similar, showing large-scale (>30 cm) variation in the gravimetric densities of each core section. Fine-scale features (<10 cm) are well characterized by the COS record and are seen at a slightly lower resolution in both the DEP and NP records, which show increasing smoothing. A combination of the density accuracy of NP and the spatial resolution of COS provides a useful method of evaluating the shallow-density profile of a glacier, improving paleoclimate interpretation, mass-balance measurement and interpretation of radar returns.

This study adjusts and compares digital elevation models (DEMs) created from photogrammetric and interferometric synthetic aperture radar techniques to determine volume and surface elevation changes of five icefields in a remote region of southwest British Columbia, Canada, between the mid-1980s and 1999. Preliminary differences between the DEMs in ice-free and vegetation-free areas indicated variable elevation offsets with increasing altitude (11 m km−1) and with increasing slope (2.7 m (10°)−1). Results indicate a surface elevation change of −6.0 ± 2.7 m (−0.5 ± 0.2 m a−1) and a total volume loss of −19.4 ± 8.8 km3 (−1.5 ± 0.7 km3 a−1), which represents a potential sea-level rise contribution of 0.004 ± 0.002 mm a−1. Temperature and snowfall data from four nearby meteorological stations indicate that increased temperatures and decreased snowfall throughout the late 1980s and 1990s are a likely cause of the thinning. Glacier terminus positions were compared between a historical map (1927) and satellite images (1974, 1990/91 and 2000/01). All observed glaciers were in retreat between 1927 and 1974, as well as between 1990/91 and 2000/01, but many glaciers advanced or significantly slowed in their retreat between 1974 and 1990/91.

Surface-elevation change and ice velocities have been measured over the debris-covered tongue of Chacraraju, Peru. Elevation change was measured by reflectorless survey at a 1 m horizontal resolution over three separate areas of the glacier between 2004 and 2005. Area-averaged change revealed general lowering, with two of the surveyed areas experiencing surface lowering of 0.58 and 0.77 m, and the third a rise of 0.07 m. Combining all three areas (43 216 m2) resulted in a mean net lowering of 0.43 m a−1, which is at the higher end of the range of long-term studies in the region. Velocity was measured over 7 days by the repeated optical survey of 12 prisms attached to stakes inserted directly into the glacier’s surface. Results indicate that velocity increases approximately with distance squared from the glacier’s terminus, from <10 mm d−1 near the terminus to approaching 100 mm d−1 at the base of the glacier’s icefall, located ∼1.7 km up-glacier. Velocity vectors also changed systematically along the glacier, from a consistent down-glacier orientation near the icefall to more variable orientations within ∼300 m of the terminus. No up-glacier motion component was measured.

A finite-volume model is used to simulate 9 years (1995–2003) of snow temperatures at the South Pole. The upper boundary condition is skin-surface temperature derived from routine upwelling longwave radiation measurements, while the lower boundary condition is set to the seasonal temperature gradient at 6.5 m depth, taken from prior measurements at the South Pole. We focus on statistics of temperature, heat fluxes, heating rates and vapour pressures in the top metre of snow, but present results from the full depth of the model (6.5 m). The monthly mean net heat flux into the snow agrees with results from previous studies performed at the South Pole. On shorter timescales, the heating rates and vapour pressures show large variability. The net heat flux into the snow, which is between ±5 W m−2 in the monthly mean, can be greater than ±20 W m−2 on hourly timescales. On sub-daily timescales, heating rates exceed 40 K d−1 in the top 10 cm of the snow. Subsurface temperatures, and therefore heating rates, are more variable during winter (April–September) due to increased synoptic activity and the presence of a strong, surface-based, atmospheric temperature inversion. The largest vapour pressures (60–70 Pa) and vertical gradients of vapour pressure are found in the top metre of snow during the short summer (December–January). In contrast, during the long winter, the low temperatures result in very small vapour pressures and insignificant vapour-pressure gradients. The high summertime vapour-pressure gradients may be important in altering the isotopic composition of snow and ice on the Antarctic plateau.

We present spatial patterns of glacier fluctuations from the Cordillera Blanca, Peru, (glacier area, terminus elevations, median elevations and hypsography) at decadal timescales derived from 1970 aerial photography, 2003 SPOT5 satellite data, Geographic Information Systems (GIS) and statistical analyses. We derived new glacier outlines from the 2003 SPOT images, and ingested them in the Global Land and Ice Measurements from Space (GLIMS) glacier database. We examined changes in glacier area on the eastern and western side of the Cordillera in relation to topographic and climate variables (temperature and precipitation). Results include (1) an estimated glacierized area of 569.6 ± 21 km2 in 2003, (2) an overall loss in glacierized area of 22.4% from 1970 to 2003, (3) an average rise in glacier terminus elevations by 113 m and an average rise in the median elevation of glaciers by 66 m, showing a shift of ice to higher elevations, especially on the eastern side of the Cordillera, and (4) an increase in the number of glaciers, which indicates disintegration of ice bodies. Annual air temperature showed a significant upward trend in the last 30 years, with larger temperature increases at lower elevations. There was a slight but not significant decrease in precipitation. Our results are consistent with glacier retreat and warming trends noted in the last three decades in the tropics.

Using dense networks of three-component seismometers installed in direct contact with the ice, the seismic activity of Gornergletscher, Switzerland, was investigated during the summers of 2004 and 2006, as subglacial water pressures varied drastically. These pressure variations are due to the diurnal cycle of meltwater input as well as the subglacial drainage of Gornersee, a nearby marginal ice-dammed lake. Up to several thousand seismic signals per day were recorded. Whereas most icequakes are due to surface crevasse openings, about 200 events have been reliably located close to the glacier bed. These basal events tend to occur in clusters and have signals with impulsive first arrivals. At the same time, basal water pressures and ice-surface velocities were measured to capture the impact of the lake drainage on the subglacial hydrological system and the ice-flow dynamics. Contrary to our expectations, we did not observe an increase of basal icequake activity as the lake emptied, thereby raising the subglacial water pressures close to the flotation level for several days. In fact, the basal icequakes were usually recorded during the morning hours, when the basal water pressure was either low or decreasing. During the high-pressure period caused by the drainage of the lake, no basal icequakes were observed. Furthermore, GPS measurements showed that the glacier surface was lowering during the basal seismic activity. These observations lead us to conclude that such icequakes are connected to the diurnal variation in glacier sliding across the glacier bed.

The mass balance of Glacier de Sarennes, French Alps, has been measured since 1949, using the glaciological method based on core and ablation stake data, and area extrapolations, to find the overall glacier balance. The cumulative balance obtained in this way is very dependent on systematic errors that can increase linearly with the number, N, of measurement years, whereas random errors rise with . The volumetric-balance method based on aerial photogrammetry provides results whose errors do not depend on the number of years. This method was used to test field measurements for the period 1952–2003 and gives a mass balance of −32.30 ± 1.04 m w.e. compared to −34.89 ± 1.15 m w.e. based on field data. The discrepancy between the two methods is discussed on the basis of a careful error analysis. Moreover, the possibility of using the volumetric method to detect biases in field measurements is evaluated in terms of two types of errors. The number and locations of measurement sites required to account for all the spatial and temporal variabilities of the mass balance is discussed by variance analysis. Methodological implications and recommendations are presented to provide mass-balance measurements of the best possible accuracy.

The microstructures of particle-free granular freshwater ice and ice containing 1 wt.% of 50 ± 10 mm uniformly distributed particles were investigated before and after compressive creep to ∼10% strain with stresses of 1.45 MPa at −10°C and 0.4 MPa at −5°C. Creep rates of particle-containing ice were always higher than those of particle-free ice. For an initial stress of 1.45 MPa at −10°C, dynamic recrystallization occurred with new grains nucleating and growing along grain boundaries for both sets of specimens, and the ice with particles showed a higher nucleation rate. Under creep with an initial stress of 0.4 MPa at −5°C, dynamic recrystallization also occurred by the nucleation and growth of new grains along the grain boundaries for ice containing particles, but recrystallization in the particle-free ice occurred through grain boundary migration.

The altitude of the 0°C isotherm obtained from radiosonde data of the aerological Chilean stations Antofagasta, Quintero/Santo Domingo, Puerto Montt and Punta Arenas are analyzed, along with surface temperature and precipitation records from nearby stations. The strong effect of the 1976/77 climate shift due to a change in the Pacific Decadal Oscillation is evident in the temperature and precipitation data. The data are used as input for an empirical model which reconstructs annually the equilibrium-line altitude (ELA) for the last 49 years on the western side of the southern Andes. The model takes air temperature, precipitation and altitude as the main parameters, and was first developed by Fox (1993) and applied by Condom and others (2007). From the radiosonde data, a significant positive trend of the 0°C isotherm has occurred in the northern, central and southern regions, indicating an ELA rise due to regional warming. General glacier retreat, ice thinning and negative mass balance observed during the past few decades in virtually all the Chilean Andes concur with the observed ELA reconstruction. In the Punta Arenas radiosonde record there is slight evidence for precipitation increase but no evidence for significant warming in the past few decades. This results in a slight lowering of the ELA according to the model reconstruction, which does not agree with the strong and increased glacier retreat observed in recent decades in Patagonia.

We describe an automated method to generate an inventory of glaciers and glacier morphometry from a digital topographic database containing glacier boundaries and a digital elevation model for British Columbia, Canada. The inventory contains over 12 000 glaciers with a total cumulative area that exceeds 25 000 km2, based on mapping from aerial photographs circa the mid-1980s. We use the inventory to examine dimensional characteristics among glaciers, namely the scaling relations between glacier length, width and area. Glacier length is a good predictor of glacier area, and its predictive ability improves when glaciers are stratified by the number of up-valley accumulation basins. The spatial pattern of glacier mid-range altitude parallels glaciation limits previously mapped for British Columbia and similarly reflects large-scale controls of orographic precipitation and continentality. The inventory is also used to refine models that relate glacier mid-range and terminus altitudes to regional position, aspect and, in the case of terminus altitudes, an index of glacier shape. Relations between glacier altitude limits and controlling spatial and topographic factors are used to make further climatic and mass-balance inferences from the glacier inventory.