We examine the validity of two methods for estimating glacier equilibrium-line altitudes (ELAs) from topographic maps. The ELA determined by contour inflection (the kinematic ELA) and the mean elevation of the glacier correlate extremely well with the ELA determined from mass-balance data (observed ELA). However, the range in glacier elevations above sea level is much larger than the variation in ELA, making this correlation unhelpful. The data were normalized and a reasonable correlation (r2 = 0.59) was found between observed and kinematic ELA.The average of the normalized kinematic ELAs was consistently located down-glacier from the observed ELA, consistent with theory. The normalized mean elevation of the glacier exhibited no correlation and suggests that the toe–headwall altitude ratio is not a good approximation for the ELA. Kinematic waves had no effect on the position of the kinematic ELA. Therefore, topographic maps of glacier surfaces can be used to infer the position of the ELA and provide a method for estimating past ELAs from historic topographic maps.

The surface velocity of a predominantly cold polythermal glacier (John Evans Glacier, Ellesmere Island, Canada) varies significantly on both seasonal and shorter time-scales. Seasonal variations reflect the penetration of supraglacial water to the glacier bed through significant thicknesses of cold ice. Shorter-term events are associated with periods of rapidly increasing water inputs to the subglacial drainage system. Early-season short-term events immediately follow the establishment of a drainage connection between glacier surface and glacier bed, and coincide with the onset of subglacial outflow at the terminus. A mid-season short-term event occurred as surface melting resumed following cold weather, and may have been facilitated by partial closure of subglacial channels during this cold period. There is a close association between the timing and spatial distribution of horizontal and vertical velocity anomalies, the temporal pattern of surface water input to the glacier, and the formation, seasonal evolution and distribution of subglacial drainage pathways. These factors presumably control the occurrence of highwater-pressure events and water storage at the glacier bed. The observed coupling between surface water inputs and glacier velocity may allow predominantly cold polythermal glaciers to respond rapidly to climate-induced changes in surface melting.

Permeability and meltwater flow have been studied in sea ice in the Siberian and central Arctic during the summers of 1995 and 1996. A bail-test technique has been adapted to allow for measurements of in situ permeability, found to range between 10−11 and 10−8 m2. Permeability varied by about a factor of 2 between 1995 (above-normal melt rates) and 1996 (below-normal melt rates). Release of fluorescent tracers (fluorescein, rhodamine) furthermore allowed the derivation of flow velocities and assessment of the relevant driving forces. Hydraulic gradients in rough ice and wind stress in ponded ice were found to be particularly important, driving meltwater over distances of several meters per day. The mid- to late-summer ice was found to be permeable enough to completely divert meltwater from the surface into the ice interior. It is shown, however, that lower permeabilities of the upper ice layers as well as refreezing of meltwater, particularly during the early melt season, foster the development of surface melt ponds.

The effective viscosity of ice depends upon many factors, including temperature, deviatoric stress, crystal orientation and impurities. A flow law that includes these factors and is simple to implement is a requirement for numerically efficient ice-flow models. The dominant microscale flow mechanism changes as temperature, deviatoric stress or grain-size changes. For both anisotropic and isotropic constitutive relations, this shift in dominant flow mechanism is expressed as a change in the stress exponent. We study the effects of this shift in stress exponent on ice flow using a two-term flow law for isotropic ice. Our stress–strain-rate relationship does not explicitly describe the microscale processes of ice deformation; however, it encompasses a range of deformation behaviors with a simple law. In terrestrial ice, a flow-mechanism shift may occur in low-deviatoric-stress regions near ice divides, resulting in a near-linear constitutive relationship for ice flow. Compared to a non-linear (Glen) divide, a divide dominated by a near-linear flow mechanism has vertical-velocity profiles that are similar at divide and flank sites, internal layers that do not develop a Raymond bump, and a steady-state surface profile that is more rounded near the divide.

Shallow lakes cover >25% of Alaska’s Arctic Coastal Plain. These remain frozen and snow-covered from October to June. The lake snow is thinner, denser, harder and has less water equivalent than snow on the surrounding tundra. Itcontains less depth hoar than land snow, yet paradoxically is subject to stronger temperature gradients. It also has fewer layers and these have been more strongly affected by wind. Dunes and drifts are better developedon lakes; they have wavelengths of 5–20 m, compared to <5 m on land. Because of these differences, lake snow has roughly half the thermal insulating capacity of land snow. The winter mass balance on lakes is also different because (1) some snow falls into the water before the lakes freeze, (2) some snow accumulates in drifts surrounding the lakes, and (3) prevailing winds lead to increased erosion and thinner snow on the eastern lake sides. Physical models that extrapolate land snow over lakes without appropriate adjustments for depth, density, distribution and thermal properties will under-predict ice thickness and winter heat losses.

A major surge of the glacier complex of Sortebræ, East Greenland, occurred between 1992 and 1995. The impact of this surge on the topography of Sortebræ was examined through the production of a pre-surge (1981) digital elevation model (DEM) from interpolated digital map data, and a post-surge DEM from differential synthetic aperture radar interferometry using European Remote-sensing Satellite (ERS) imagery. The combined vertical error is typically 30–40 m; however, downdraw of up to 270 m in the reservoir zone, and uplift of up to 145 m in the receiving zone were measured. The upper glacier reservoir area discharged in excess of 24.3 ± 9.5 km3 volume over the surge, of which ∼12.5 km3 was stored in the advanced lower glacier, the balance being lost, predominantly to calving.

As a complement to earlier measurements on the friction of both granular fresh-water ice and S2 columnar salt-water ice, new experiments were performed on the friction of S2 columnar fresh-water ice sliding against itself at low velocities (5 × 10−7 to 5 × 10−1 m s−1) and at −10°C, using the same double-shear device as was used earlier. The results showed that under a given set of experimental conditions the kinetic coefficient of friction of S2 fresh-water ice compares favorably with that of the other two variants.The experiments also revealed friction-induced surface cracks and recrystallized grains.These deformation features are explained, respectively, in terms of fracture mechanics and an earlier model of dynamic recrystallization in ice.

Quantitative ice-core paleoclimatology must account for post-depositional processes, such as vapor-phase diffusion in the firn. After pore close-off, diffusion continues to smooth the stable-isotope records δ18O and δD that are eventually recovered from the ice, leading to the loss of high-frequency information. Johnsen and others (1997) found much higher rates of diffusive smoothing in the Greenland Icecore Project (GRIP) Holocene ice than would be predicted by diffusion through solid ice alone, and Nye (1998) argued that transport through liquid veins might explain this apparent excess diffusion. However, the analysis of Johnsen and others (2000) indicates that the required vein dimensions may be unrealistically large. Here, we model the diffusion of stable isotopes in polycrystalline ice and show that the predictions of Nye (1998) and those of Johnsen and others (2000) actually represent two end-members in a range of potential behavior. Our model determines which of these asymptotic regimes more closely resembles the prevailing conditions and quantifies the role of pre-melted liquid in the smoothing of isotopic signals. The procedure thereby ties together the two approaches and provides a rostrum for accurate analysis of isotope records and paleotemperature reconstructions.

The longitudinal pattern of surface velocity of a large, predominantly cold, polythermal glacier (John Evans Glacier, Ellesmere Island, Canada) was measured over summer and winter periods. In the accumulation and upper ablation areas, where ice is predominantly cold-based, summer velocities were slightly higher than overwinter velocities. Predicted velocities due to ice deformation alone in these areas closely matched these observations in the winter, with limited basal motion likely in the summer. In the lower ablation area, where ice is likely warm-based, measured summer velocities were up to double overwinter velocities. Predicted ice deformation could not account for all of these measured velocities in either summer or winter. This suggests that basal motion occurs throughout the year over at least part of the lower ablation area. This finding is supported by radio-echo sounding, subglacial drainage reconstructions and analyses of early-summer meltwater chemistry, which suggest that subglacial water is present throughout the year in this region. In summer, basal motion may account for up to 75% of the total surface velocity throughout the lower ablation area. The inferred rate of basal motion increases sharply directly below a set of moulins by which most surface meltwater reaches the glacier bed.

We present some new results describing energy exchange processes of drifting sea ice in the marginal ice zone (MIZ) in the Barents Sea, Arctic Ocean. All measurements and observations of meteorological parameters and ice conditions were taken on board the Norwegian research vessel Lance from 3 to 22 May 1999. Components of surface heat balance were measured and correlated with ice conditions and synoptic observations. These results can be used in atmospheric boundary layer modelling as lower boundary conditions. A relationship was found between modelled turbulent heat fluxes and observed sea-ice concentrations.

Glacier responses to tephra deposition are shown to be highly variable where wind-transported eruption plumes produce narrow distal fallout zones with steep lateral thickness gradients. Significant but short-lived advances of faster-flowing glaciers can be triggered by deposition from modest eruptions. The 1947 eruption of Hekla, south Iceland, covered nearby glaciers with variable thicknesses of tephra, causing dramatic spatial variations in ablation rate. Relative snow ablation increased by an estimated >80% at the eastern margin of Eyjafjallajökull, but decreased by >54% at the western margin. Relative ice ablation increased by only 4% in the east but decreased by >75% at the western margin, only 15 km distant. The effect on mass balance therefore depends on tephra distribution as well as on the nature of the glacier surface. On Gígjökull, retardation of ablationwas greatest at the terminus, associated with an anomalous advance of ≥328 m between 1947 and 1954. Other glaciers in the fallout zone show either no recognizable response, or accelerated retreat. Advances will be short-lived due to the rapid redistribution of supraglacial tephra in a maritime climate, and the potential for discriminating between volcanically and climatically forced advances in the glacio-geomorphological record is low. The long-term glaciological effect of volcanism is to create thin, low-albedo covers after reworking, which make conditions less favourable for glaciation.

The transition of inland to complex ice flow is investigated in the downstream area of Shirase Glacier catchment, East Antarctica. The 250 km long Mizuho Triangulation Chain (MTC; Naruse, 1978), 230 km upstream of the grounding line of Shirase Glacier, is re-analyzed by a force-balance analysis, in combination with balance-velocity estimates. Numerical model experiments were carried out with a two-dimensional thermomechanical higher-order ice-sheet model along a central flowline extending from the ice divide (Dome Fuji) to the grounding line of Shirase Glacier, passing through the MTC. Results show that complex flow originates near the MTC, where longitudinal normal stresses and normal drag in particular play a decisive role in the force balance. The fast-flow area, however, is geologically controlled and confined to the “bottleneck” or ice-flow convergence area, a few tens of kilometres upstream of the grounding line. In between the MTC and the grounding zone, longitudinal normal stress and basal sliding are dominant features of the ice-flow regime.

We present an analysis of many European Remote-sensing Satellite (ERS) synthetic aperture radar (SAR) images of Vatnajökull, Iceland, by comparing them with Advanced Very High Resolution Radiometer (AVHRR) images, mass-balance observations and modelled firn stratigraphy. Summer SAR and AVHRR images both detect the surface firn line (or the transient snow-line when it lies below the firn line) as a distinct boundary.Winter SAR images of Vatnajökull generally display the late-summer surface firn line, and often also a transition at a higher elevation that possibly reflects a subsurface continuation of the firn line. On some winter images the firn line is not adequately detected, possibly due to a shallow firn pack. We found no differences in back-scatter between melting snow and firn, and for most years no differences in reflectance either. For Vatnajökull, equilibrium-line altitude is therefore not a particularly useful estimator of the mean specific mass balance (Bm). Unlike albedo images, SAR images do not display interannual variations of the signal within the accumulation area that are clearly related to Bm. Hence, SAR images of Vatnajökull contain less information about Bm than albedo images. However, for some drainage basins, mean transient snow-line/firn-line altitude during the melting season can be used to infer Bm.

Glaciar Coronas is one of the small cirque glaciers located in the Maladeta Range, central Spanish Pyrenees. The ice extends from 3220 to 3065 m a.s.l. Using early and modern photographs and aerial photographs, we study the evolution of the glacier during the 20th century, with particular focus on evaluating the variations in area, ice thickness and the withdrawal of the ice front. This information is compared with some basic climatic data (quarterly and annual evolution of temperatures and precipitation) recorded in the nearby Seira observatory. The results show that ice wastage, clearly associated with the main regional climatic evolution, has been dominant during the 20th century at Glaciar Coronas, but the rate has varied and has not been continuous. Three periods have been detected: (1) 1920s–early 1960s: ice depletion period; (2) early 1960s–mid-1980s: stabilization of the deglaciation process; (3) mid-1980s–1998: marked ice depletion period. The study shows a decrease in area of Glaciar Coronas from 19 ha in 1928 to 2.4 ha in 1998. The decrease in length from 1928 to 1998 has been 325 m.

A simple model uses once-daily meteorological values in the U.S. National Centers for Environmental Prediction and U.S. National Center for Atmospheric Research (NCEP–NCAR) re-analysis database to estimate summer balance of South Cascade Glacier, Washington, U.S.A., each year over 1959–99. The rms error, 0.30 m w.e. (r2 = 0.71), is comparable to measurement error. The model relates summer balance linearly to temperature T > 0°C at 2000 m andto snow flux at 1650 m, the altitudes in recent years of the equilibrium line and terminus. The snow flux is the product of the humidity and westerly wind component at 850 hPa when temperature T <+2°C at 1650 m. Temperatures are interpolated linearly between the 850 and 700 hPa levels. Both the positive 2000 m temperature and the snow flux are summed from 26 April to 4 October. When the summer estimates are combined with those from a winter balance model using the same database, the rms error in estimating net balance is 0.40 m (r2 = 0.81). The indicated sensitivities of balance to warming of 1°C are −0.51 m for summer and −0.24 m for winter. On the assumption that the total −0.75 m °C−1 sensitivity exists at all altitudes, a warming of only 0.7°C would be sufficient to overcome the 1986–98 average net balance +0.5 m at the top of the glacier.

Radar returns from valley side-walls may complicate the identification of the glacier bed when examining airborne radio-echo sounding (RES) records. We describe simple visualization techniques for aiding analysis of RES data with the aim of resolving these ambiguities, making use of high-resolution satellite imagery. An application of the technique is illustrated using RES data obtained from part of Stygge Glacier, Ellesmere Island, Nunavut, Canada.