Deformational structures at the surge-type glacier Kongsvegen, Svalbard, are displayed at the glacier surface and on a grounded cliff section at the terminus. A 300 m × 65 m grid of 200 MHz ground-penetrating radar (GPR) profiles has been collected adjacent to the cliff section in order to identify englacial structure.Two sub-horizontal reflectors have been imaged; the upper is interpreted as the glacier bed, and represents a transition between glacier ice and frozen subglacial sediments; while the lower is interpreted as a transition between frozen and unfrozen subglacial sediment. Dipping reflectors, corresponding to sediment-filled features on the cliff and glacier surface, do not cross the glacier bed. A small number of reflectors, interpreted as thrust faults, are visible below the bed reflector. A model is developed for structural development, which suggests that ice built up in a reservoir zone during quiescence. During the surge, ice propagated rapidly from this reservoir, creating a zone of compression which resulted in thrusting. Subsequently an extensional flow regime resulted in extensive fracture of the ice. We suggest dilated sediment was evacuated into these extensional crevasses from the glacier bed, accelerating surge termination.

Observations of surface elevation (s) and horizontal velocity components (u and v) are inverted to infer the topography (b) and lubrication (c) at the bed of an ice stream, based on a linearized perturbation theory of the transmission of flow disturbances through the ice thickness. Synthetic data are used to illustrate non-uniqueness in the inversion, but also demonstrate that effects of b and c can be separated when s, u and v are specified, even with added noise to simulate measurement errors. We have analyzed prominent short-horizontal-scale (∼2 km) features in topography and velocity pattern in a local 64 km by 32 km area of the surface of Ice Stream E,West Antarctica. Our preferred interpretation of bed conditions beneath the most prominent features on the surface identifies a deep trough in the basal topography with low lubrication in the base of the trough.

This investigation uses Landsat images from 12 days in 1999 and 2000 to study the spatial and temporal variation in surface albedo of a glacier with a rugged topography: Morteratschgletscher, Switzerland. Our retrieval method considers all processes that substantially influence the relationship between the satellite signal and the surface albedo. The correction for the anisotropy of the reflected radiation field of ice and snow ranges up to 0.10, depending on wavelength band, solar zenith angle and surface type. We analyzed the uncertainties in the retrieval method and mainly expect errors in satellite-derived albedos for areas with large variation in topography and high albedos. The latter is due to application of bidirectional reflectance distribution function (BRDF) parameterizations for the anisotropic correction to albedos outside the parameterizations’ validity. On average, satellite-derived albedos exceed the measured surface albedo by 0.03. The glacier tongue is characterized by bands of low and high albedo, relating to ice with higher and lower concentrations of debris. The ice albedo shows no dependence on altitude, except at 2000–2200 m a.s.l. It increases during summer, which is likely associated with rainfall, as concluded from a comparison between summer rainfall and measured albedos.

Automatic feature tracking on two Landsat images (acquired inJanuary 2000 and December 2001) generates a complete and accurate velocity field of Mertz Glacier, East Antarctica. This velocity field shows two main tributaries to the ice stream. Between the tributaries, a likely obstruction feature in the bedrock results in a slow-down of the flow. A third Landsat image, acquired in 1989 and combined with the 2000 image, permits the determination of the glacier mean velocity during the 1990s. Although some parts of the Mertz Glacier system show evidence of slight speed increase, we conclude that the Mertz flow speed is constant within our uncertainty (35 m a−1). Using this complete velocity field, new estimates of the ice discharge flux, 17.8 km3 a−1 (16.4 Gt a−1), and of the basal melting of the tongue, 11 m a−1 of ice, are given. Our results lead to an apparent imbalance of the drainage basin (ice discharge 3.5 km3 a−1 lower than the accumulation). Considering previous studies in the Mertz Glacier area, we then discuss the uncertainty of this imbalance and the problems with accumulation mapping for this region.

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.

The inverse barometer effect (IBE) is the isostatic response of ocean surface height to changes in atmospheric pressure (P air) at a rate of about 1 cm hPa−1. The IBE is a significant contributor to variability of ice-shelf surface elevation (η ice), as we demonstrate with simultaneous global positioning system measurements of η ice and local measurements of P air from the Amery, Brunt and Ross Ice Shelves, Antarctica. We find that an IBE correction is justified for frequencies (ω) covering the “weather band”, 0.03 < ω < 0.5 cpd (cycles per day). The IBE correction reduces the standard deviation of the weather-band signal of η ice from ∼9 cm to ∼3 cm. With this correction, the largest remaining high-frequency error signal in η ice is the inaccuracy of the present generation of Antarctic tide models, estimated to be of order 10 cm for most of Antarctica.

A two-dimensional ice-flow model based on the shallow-ice approximation (SIA) is used to investigate the dynamics of Glacier de Saint-Sorlin, France. This glacier is well suited for this kind of study. First, the particular geometry of the glacier itself as well as that of the bedrock surface allows for correct applicability of the SIA (zeroth-order equations), provided that thickness changes and termini positions rather than short-scale dynamics are considered. Secondly, the wealth of available data for the glacier including mass-balance series and records of glacier changes provides a reliable forcing and a powerful constraining set for the model. Steady-state simulations show realistic results and the capabilities of the model in reproducing the glacier extent at the beginning of the 20th century. An extensive parameter study of ice rheology and sliding intensity is also carried out and the results are checked against the thickness changes as well as the glacier termini positions since 1905. It is possible to find a parameter combination that best matches these two types of data with an ice-flow rate factor of 2 × 10−24 Pa−3 s−1 and a Weertman-type sliding factor of 5 × 10−14 m8 N−3 a−1 which both appear to be in agreement with similar inferences from recent modelling attempts.

Glaciers are an important element of the environment in the source regions of the Yangtze and Yellow rivers on the Tibetan Plateau. Using Geographical Information System techniques, we have studied changes in the location of glacier margins in two areas: the Geladandong area in the headwaters of the Yangtze, and the A’nyêmaqên Shan mountains in the headwaters of the Yellow River. Marginal positions during the Little Ice Age (LIA) maximum, in 1969 in the Geladandong area, in 1966 in the A’nyêmaqên Shan, and in 2000 in both areas, were determined using aerial photographs, satellite images, topographical maps and digital elevation models. Extrapolating the results to the entire source regions of the Yangtze and Yellow rivers, we estimate that the total glacierized area decreased about 1.7% between 1969 and 2000 in the Geladandong area and about 17% between 1966 and 2000 in the A’nyêmaqên Shan. Glaciers were stable or advanced slightly between 1969 and 1995 in the Geladandong area, and between 1966 and 1981 in the A’nyêmaqên Shan, but have retreated since the mid-1990s in the former and since the 1980s in the latter. Significant increases in summer air temperature and decreases in annual precipitation are the causes of the present retreat. As a consequence of the retreat, water storage, as ice, in the Yangtze and Yellow river drainages is decreasing by 65–70 × 106 m3 a−1.

Catchment-wide information on glacier snow-cover depth, surface albedo and surface roughness is important input data for distributed models of glacier energy balance. In this study, we investigate the small-scale (mm to 100 m) spatial variability in these properties, with a view to better simulating this variability in such models. Data were collected on midre Lovénbreen, a 6 km2 valley glacier in northwest Svalbard. The spatial variability of all three properties was found to be self-similar over the range of scales under investigation. Snow depth and albedo exhibit a correlation length within which measurements were spatially autocorrelated. Late-winter and summer properties of snow depth differed, with smaller depths in summer due to melt, and shorter correlation lengths. Similar correlation lengths for snow depth and surface albedo may suggest that snow-depth variation is an important control on the small-scale spatial variability of glacier surface albedo. For surface roughness, the data highlight a possible problem in energy-balance studies which use microtopographic surveys to calculate aerodynamic roughness, in that the scale of the measurements made affects the calculated roughness value. This suggests that further investigations of the relationships between surface form and aerodynamic roughness of glacier surfaces are needed.

Three early-melt-season high-velocity events (or “spring events”) occurred on Haut Glacier d’Arolla, Switzerland, during the melt seasons of 1998 and 1999. The events involve enhanced glacier velocity during periods of rapidly increasing bulk discharge in the proglacial stream and high subglacial water pressures. However, differences in spatial patterns of surface velocity, internal ice deformation rates, the spatial extent of high subglacial water pressures and in rates of subglacial sediment deformation suggest different hydrological and mechanical controls. The data from two of the events suggest widespread ice–bed decoupling, particularly along a subglacial drainage axis creating the highest rates of basal motion and “plug flow” in the overlying ice. The other event showed evidence of less extensive ice–bed decoupling and sliding along the drainage axis with more mechanical support for ice overburden transferred to areas adjacent to decoupled areas. We suggest that: (1) plug flow may be a common feature on glaciers experiencing locally induced reductions in basal drag; (2) under certain circumstances, enhanced surface motion may be due in part to non-locally forced enhanced bed deformation; and (3) subglacial sediment deformation is confined to a depth of the order of centimetres to decimetres.

Retreat of Mendenhall Glacier nearJuneau, Alaska, U.S.A., has exposed a bedrock ridge spotted with “siltskins”, patchy coatings of calcite-cemented clay-to sandsized lithic grains. Coatings 0.5–20 mm thick occur in two distinct morphologies. Thin, striated siltskins coat mainly stoss faces. Thicker, corrugated siltskins on lee faces consist of parallel micro-ridges elongated downslope. Thin-section analysis shows that siltskins consist of a basal, calcite-rich layer overlain by microlaminated layers of calcite-cemented lithic grains. Scanning electron microscope (SEM) energy-dispersive spectrometer (EDS) analysis of laminae and surfaces shows laterally persistent Ca/Si differences. Isotopic values for δO18 and δO13 ranged from −19.52‰ to −12.74‰ and −6.18‰ to −3.44‰, respectively, consistent with deposition from subglacial waters of varying isotopic composition and with derivation of carbon from inorganic sources. Corrugated siltskins are complex depositional features modified by erosional processes. Parallel micro-ridges spaced 1–10 mm apparently formed as sediment-rich water dripped down lee-slope rock faces. Ice–rock separation, flow energy and the transported sediment controlled the layering and depositional forms. Siltskins probably formed when a subglacial cavity system was active on the rock ridge and provide clues about how microscale hydrologic processes interact with larger-scale subglacial systems.

Using the finite-element we have modeled the stress field near the calving face of an idealized tidewater glacier under a variety of assumptions about submarine calving-face height, subaerial calving-face height, and ice rheology. These simulations all suggest that a speed maximum should be present at the calving face near the waterline. In experiments without crevassing, the decrease in horizontal velocity above this maximum culminates in a zone of longitudinal compression at the surface somewhat up-glacier from the face. This zone of compression appears to be a consequence of the non-linear rheology of ice. It disappears when a linear rheology is assumed. Explorations of the near-surface stress field indicate that when pervasive crevassing of the surface ice is accounted for in the simulations (by rheological softening), the zone of compressive strain rates does not develop. Variations in the pattern of horizontal velocity with glacier thickness support the contention that calving rates should increase with water depth at the calving face. In addition, the height of the subaerial calving face may have an importance that is not visible in current field data owing to the lack of variation in height of such faces in nature. Glaciers with lower calving faces may not have sufficient tensile stress to calve actively, while tensile stresses in simulated higher faces are sufficiently high that such faces will be unlikely to build in nature.

Knowledge of iceberg calving is important for understanding instabilities of tidewater glaciers and ice sheets. Since 1995 the terminus of LeConte Glacier, Alaska, U.S.A., has retreated about 2 km and the glacier has thinned approximately 120 m at its 1999 terminus position. Our focus is short-term (hours to weeks) variability of the frequency and magnitude of calving events and calving flux. Both photogrammetric and visual observations are employed in a temporal analysis over a several-week period. We combined these data with measurements of ice speed, tide level, surface water input and water-storage estimates in an attempt to better understand the calving process. Contrary to results obtained over longer time-scales on other glaciers, our results show no correlation between ice speed and the frequency of calving. However, calving events do not appear to occur randomly; often they are a response to measurable changes in other parameters within the terminus region. Caclving can often be attributed to buoyancy perturbations and possibly flexure of the nearly floating terminus. Given the multiple possibilities for buoyancy perturbations, we have found no simple relationship between any specific forcing parameter and calving at short time-scales.

This paper describes the use of a neutron probe to measure detailed stratigraphy in ice and snow. The Wallingford neutron probe, developed for measurement of soil moisture, consists of an annular radioactive source of fast neutrons around the centre of a cylindrical detector for slow (thermal) neutrons. In snow and ice, the fast neutrons lose energy by scattering from hydrogen atoms, and the number of slow neutrons arriving at the detector (the count rate) is related to the density of the medium. Calibration equations for count rate as a function of snow density and borehole diameter have been derived. Snow-density profiles from boreholes obtained using the probe show that, despite the smoothing produced by the neutron-scattering process, annual variations in density can be resolved. The potential contribution of the neutron probe to improvements in mass-balance monitoring is discussed.