Glacier and lake variations in the Yamzhog Yumco basin in southern Tibet were studied by integrating series of spatial data from topographic maps and Landsat images at three different times: 1980, 1988/90 and 2000. The results indicate that the total glacier area has decreased from 218 km2 in 1980 to 215 km2 in 2000, a total reduction of 3 km2 (i.e. a 1.5% decrease). Glacier recession rates were clearly larger in the 1990s than the 1980s due to the warmer climate. The total lake area decreased by about 67 km2 during 1980–90 and increased by 32 km2 during 1990–2000. It is suggested that change of lake area in the basin was rapid and most likely caused primarily by the change in precipitation and evaporation in the basin, and secondarily by the increased water supply from melting glaciers.

Field measurements of surface deformation across the northern shear margin of Whillans Ice Stream, West Antarctica, were analyzed to better understand the processes controlling the position and migration of this margin. Four lines of poles extending ~6km from the interstream ridge into the shear margin were installed in November 1997 and resurveyed in January 2000, using GPS methods. A band of arcuate crevasses form where the surficial shear stress reaches ~130 kPa, and a chaotic zone of crevassing occurs where this stress reaches a maximum of ~270 kPa. Inspection of the RAMP mosaic indicates that the band of arcuate crevasses is separated from the chaotic zone by a narrow zone (~250m) with a relatively undisturbed surface morphology. A force-budget analysis suggests the transition from no basal sliding to full basal sliding must be restricted to a relatively narrow band, ~1.5 km wide beneath the chaotic zone. On the ice stream, resistance from the bed is near zero and the driving stress is balanced by gradients in lateral shearing. Basal drag reaches a maximum close to the shear margin where basal conditions transition from well lubricated (beneath the ice stream) to no basal slip (beneath the ridge). The zone of elevated basal drag is ~4 km wide. Within this zone, lateral drag is reduced and becomes increasingly restricted to upper layers. Localized meltwater production under the region of elevated basal drag is similar in magnitude to melting under the ice stream. The uniform bed morphology, inferred from radar sounding, excludes the presence of abundant basal water either in a water film at the glacier bed or in a thick saturated till layer. This suggests that most of the meltwater produced under the shear margin is diverted towards the ice stream, possibly playing a crucial role in maintaining lubricated conditions allowing streaming flow.

We use a global coupled ocean-sea ice model to test the hypothesis that the disintegration of the West Antarctic ice sheet (WAIS), or just its ice shelves, may modify ocean circulation and sea-ice conditions in the Southern Ocean. We compare the results of three model runs: (1) a control run with a standard (modern) configuration of landmask in West Antarctica, (2) a no-shelves run with West Antarctic ice shelves removed and (3) a no-WAIS run. In the latter two runs, up to a few million square kilometres of new sea surface area opens to sea-ice formation, causing the volume and extent of Antarctic sea-ice cover to increase compared with the control run. In general, near-surface waters are cooler around Antarctica in the no-shelves and no-WAIS model runs than in the control run, while warm intermediate and deep waters penetrate further south, increasing poleward heat transport. Varying regional responses to the imposed changes in landmask configuration are determined by the fact that Antarctic polynyas and fast ice develop in different parts of the model domain in each run. Model results suggest that changes in the extent of WAIS may modify oceanographic conditions in the Southern Ocean.

A measure of snow density is required to estimate water equivalent ice-sheet surface mass balance (SMB) from stake measurements. Previous studies have utilized the snow density at different depths within the snow. By considering the snow densification process in the time interval between stake height measurements, we find that use of the snow density at the base of the stake is more appropriate. We assume the stakes are firmly anchored at the bottom and that Sorge’s law holds, i.e. the density–depth profile does not change with time. Applying this method to the data for 36 snow stakes on Dome Fuji, the SMB in 2003 was 36.5 kg m-2 a−1, 27% larger than the previous estimate, which used surface snow density. Correct selection of the snow density for SMB estimations is important, especially for Antarctic inland areas where accumulation is low (e.g. Dome Fuji, Vostok, Dome C and South Pole) and where the snow density near the surface varies markedly.

Results of laboratory uniaxial compression tests over the stress range 0.18–0.52 MPa and the strain range 0.5–8.6% at approximately –5 and –20°C are presented. Grain-size analysis and comparisons with annealing tests confirm that grain-growth reducing processes are active during deformation. Microstructural observations reveal that subgrain-rotation recrystallization and grain-shape changes due to strain-induced grain-boundary migration are the causes of the grain-growth deceleration. Further results from microstructural observations show that obstacle formation by dislocation walls and subgrain boundaries is the reason for isotropic hardening during creep. Subgrainboundary types that are likely to be relevant for studies on the activity of different dislocation types are described.

A major difficulty in assimilating geomorphological information with ice-sheet models is the lack of a consistent methodology to systematically compare model output and field data. As an initial step in establishing a quantitative comparison methodology, automated proximity and conformity analysis (APCA) and automated flow direction analysis (AFDA) have been developed to assess the level of correspondence between modelled ice extent and ice-marginal features such as end moraines, as well as between modelled basal flow directions and palaeo-flow direction indicators, such as glacial lineations. To illustrate the potential of such an approach, an ensemble suite of 40 numerical simulations of the Fennoscandian ice sheet were compared to end moraines of the Last Glacial Maximum and the Younger Dryas and to glacial lineations in northern Sweden using APCA and AFDA. Model experiments evaluated in this manner were ranked according to level of correspondence. Such an approach holds considerable promise for optimizing the parameter space and coherence of ice-flow models by automated, quantitative assessment of multiple ensemble experiments against a database of geological or glaciological evidence.

We show that it is possible to extract a high-resolution (annual) paleoclimate record from the surface of a blue-ice area (BIA). The variability of the surface stable-isotope values suggests that almost all the surface ice in Scharffenbergbotnen BIA, East Antarctica, is of Holocene age. The isotopic changes across the BIA show that the modern climate there is warmer than the climate in the early-Holocene optimum (11 kyr BP). A volume-conserving ice flow model for the BIA constrained by isotopic variability and layer thicknesses, and a series of 14C ages indicate both that the BIA has been smaller than now, and that the surface velocities were considerably smaller during the Last Glacial Maximum. Changes in ice-sheet thickness drive the BIA towards present-day conditions. The relatively young age of the majority of the BIA also explains the lack of meteorite finds in this area, and may be typical for many BIAs in low-elevation nunatak areas.

Samples of snow and firn from accumulation zones on Clark, Commonwealth, Blue and Victoria Upper Glaciers in the McMurdo Dry Valleys (∼77–78º S, 161–164º E), Antarctica, are evaluated chemically and isotopically to determine the relative importance of local (site-specific) factors vs regional-scale influences in defining glaciochemistry. Spatial variation in snow and firn chemistry confirms documented trends within individual valleys regarding major-ion deposition relative to elevation and to distance from the coast. Sodium and methylsulfonate (MS–), for example, follow a decreasing gradient with distance from the coast along the axis of Victoria Valley (350–119 µg L−1 for Na+; 33–14 µg L−1 for MS–); a similar pattern exists between Commonwealth and Newall Glaciers in the Asgaard Range. When comparing major-ion concentrations (e.g. Na+, MS−, Ca2+) or trace metals (e.g. Al, Fe) among different valleys, however, site-specific exposures to marine and local terrestrial chemical sources play a dominant role. Because chemical signals at all sites respond to particulates with varying mixtures of marine and terrestrial sources, each of these influences on site glaciochemistry must be considered when drawing temporal climate inferences on regional scales.

The spatial pattern of accumulation rate can be inferred from internal layers in glaciers and ice sheets. Non-dimensional analysis determines where finite strain can be neglected (‘shallow-layer approximation’) or approximated with a local one-dimensional flow model (‘local-layer approximation’), and where gradients in strain rate along particle paths must be included (‘deep layers’). We develop a general geophysical inverse procedure to infer the spatial pattern of accumulation rate along a steady-state flowband, using measured topography of the ice-sheet surface, bed and a ‘deep layer’. A variety of thermomechanical ice-flow models can be used in the forward problem to calculate surface topography and ice velocity, which are used to calculate particle paths and internal-layer shapes. An objective tolerance criterion prevents over-fitting the data. After making site-specific simplifications in the thermomechanical flow algorithm, we find the accumulation rate along a flowband through Taylor Mouth, a flank site on Taylor Dome, Antarctica, using a layer at approximately 100 m depth, or 20% of the ice thickness. Accumulation rate correlates with ice-surface curvature. At this site, gradients along flow paths critically impact inference of both the accumulation pattern, and the depth-age relation in a 100 m core.

A number of improvements have been made to an existing two-dimensional ice-flow model applied to an alpine glacier. Analysis of the results of the existing model revealed several shortcomings. The first concerns the lack of mass conservation of the applied alternating-direction-implicit (ADI) scheme. A semi-implicit (SI) scheme is therefore proposed and the effects on mass conservation assessed by a comparison with the ADI scheme. The comparison is first carried out with a simple theoretical glacier for which the improvement is significant. Concerning the real case of Glacier de Saint-Sorlin, France, the initial deviation in mass conservation was much less pronounced such that the new scheme, although improving mass conservation, does not significantly change the modelled dynamics. However, other shortcomings that have a more profound impact on the modelling of glacier behaviour have been identified. The ice thickness may become negative over some gridpoints, leading to an inconsistency. The problem is partly resolved by incorporating extra checks on critical gridpoints at the glacier border. Finally, with the help of ice particle tracking, unrealistic ice settlement above the bergschrund has been identified as the main reason for spurious dynamic effects and has been corrected.

The Australian Antarctic Program’s iceberg dataset (from ship-based observations), including information from the austral summer seasons 1984/85 to 1999/2000, is examined and used to extend earlier studies. Using ‘snapshots’ of the iceberg population to provide an idea of the iceberg life cycle, the distribution of icebergs between 60 and 150° E is discussed in terms of calving regions and ocean currents. Temporal changes are also examined. The discussion leads us to the point where we can define an area, bounded to the north by the maximum sea-ice limit and to the south by the Antarctic Divergence, in which icebergs are confined as they drift eastward. This allows estimation of total dissolution, in terms of iceberg numbers and volume, within 10° longitudinal sectors and, with knowledge of drift speeds, iceberg movement rates and freshwater input across the sector. Iceberg dissolution rates are found to be ~0.03–0.05 m d–1 and the total mass contribution of fresh water to the ocean as the icebergs traverse our 30° of longitude study sector is ~32 Gt. This amounts to a contribution equivalent to precipitation of ~15.5 cm a–1, accounting for ~2% of the total iceberg discharge from the Antarctic ice sheet.

Variations in annual maximum and accumulated snow depths, snow-cover duration, precipitation and air temperature have been analyzed using daily snow depth, monthly air temperature and monthly precipitation data from 1960 to 1999 from six meteorological stations in the source regions of the Yangtze and Yellow Rivers in China. Annual maximum snow depth, snow-cover duration and precipitation increased by ~0.23, ~0.06 and ~0.05% a–1, respectively, during the study period, while annual accumulated snow depth increased by ~2.4% a–1. Annual mean air temperature increased by ~0.6°C over the study period. An unusually heavy snow cover in 1985 coincided with historically low air temperatures. Data from Tuotuohe and Qingshuihe meteorological stations are used to examine inter-station variability. The annual maximum and accumulated snow depths increased by ~0.35 and ~10.6% a–1 at Tuotuohe, and by ~0.42 and ~2.3%a–1 at Qingshuihe. However, from the late 1980s until 1999 the climate in the study region has become warmer and drier. The precipitation decrease in the 1990s (and not the rapid rise in measured temperature) is thought to be the primary cause of the decrease in snow depth in those years.