A new continuum approach to snowdrift modelling is introduced. In addition, numerical studies are carried out to identify the influence of time-varying wind conditions on snowdrift simulations. We compare the snowdrift patterns at Grimming mountain, Austria, derived using a time-averaged wind field and a time-varying wind field obtained from the numerical weather prediction model INCA. The results show significant differences in the deposition patterns and snow depth even after a 6 hour drift period. Using time-averaged boundary conditions leads to an underprediction of the resulting snow depth caused by averaging the wind speed, which lets gusts of wind disappear while snow transport is a non-linear function of the wind speed. Using numerical weather prediction models for snowdrift simulation therefore provides enhanced knowledge of the snow depth for local avalanche warning services.

Interpretation of ice mass elevation changes observed by satellite altimetry demands quantification of the proportion of elevation change which is attributable to variations in firn densification. Detailed stratigraphic logging of snowpack structure and density was carried out at ~1km intervals along a 47 km transect on Devon Ice Cap, Canada, in spring (pre-melt) and autumn (during/ after melt) 2004 and 2006 to characterize seasonal snowpack variability across the full range of snow facies. Simultaneous meteorological measurements were gathered. Spring (pre-melt) snowpacks show low variability over large spatial scales, with low-magnitude changes in density. The end-of-summer/ autumn density profiles show high variability in both 2004 and 2006, with vastly different melt regimes generating dissimilar patterns of ice-layer formation over the two melt seasons. Dye-tracing experiments from spring to autumn 2006 reveal that vertical and horizontal distribution of meltwater flow within and below the annual snowpack is strongly affected by the pre-existing, often subtle stratigraphic interfaces in the snowpack, rather than its bulk properties. Strong interannual variability suggests that using a simple relationship between air temperature, elevation and snowpack densification to derive mass change from measurements of elevation change across High Arctic ice caps may be misguided. Melt timing and duration are important extrinsic factors governing snowpack densification and ice-layer formation in summer, rather than averaged air temperatures.

Liligo glacier, in the central eastern Karakoram, Pakistan, is a small, south-to-north-flowing glacier situated in a transverse valley on the left (south) side of Baltoro glacier. New processing of satellite imagery enables a better quantification of terminus oscillations over the past 30 years. From the beginning of the 1970s to the beginning of the 21st century, Liligo glacier advanced about 2 km (60ma–1). The progress was characterized by a significant evolution of terminus morphology, similar to that observed on the same glacier during the advance event near the beginning of the 20th century, and to those of many other Karakoram glaciers. This suggests indications of a surge-type mechanism. Field observations performed in 2004 indicated there was probably no confluence at that time between Liligo and Baltoro glaciers and that a quiescent phase had started in the early years of the 21st century.

A simple combined heat and ice-sheet model has been used to calculate temperatures at the base of the Laurentide ice sheet. We let the ice sheet surge when the basal temperature reaches the pressure-melting temperature. Driving the system with the observed accumulation and temperature records from the GRIP ice core, Greenland, produces surges corresponding to the observed Heinrich events. This suggests that the mechanism of basal sliding, initiated when the basal temperature reaches the melting point, can explain the surges of the Laurentide ice sheet. This study highlights the importance of the surface temperature and accumulation rate as a means of forcing the timing and strength of the Heinrich events, thus implying important ice-sheet climate feedbacks.

Two-or three-dimensional avalanche-simulation models offer a wide range of applications; however, a challenging model-verification process is demanded, accompanied by a reliable determination of model-input parameters. We show that a verification process can be arranged with remote-monitoring data from an artificially triggered avalanche, leading to the calculation of avalanche mass balance. Two numerical methods are applied to increase the quality of the parameter fit and to reduce the number of simulations. The quality of the parameter fit is verified by comparing measured and simulated run-out lengths. In addition, a cross-check is performed using velocities derived from Doppler radar measurements.

We used remote-sensing and in situ measurements of surface accumulation rate, ice surface velocity, thickness and elevation to evaluate the mass budgets of grounded ice-flow regimes that form the Lambert Glacier–Amery Ice Shelf system. Three distinct drainage regimes are considered: the western and eastern margins of the ice shelf, and the southern grounding line at the major outlet glacier confluence, which can be identified with drainage zones 9, 11 and 10 respectively of Giovinetto and Zwally (2000). Our findings show the entire grounded portion of the basin is approximately in balance, with a mass budget of –4.2±9.8 Gt a–1. Drainages 9, 10 and 11 are within balance to the level of our measurement uncertainty, with mass budgets of –2.5±2.8 Gt a–1, –2.6±7.8 Gt a–1 and 0.9±2.3 Gt a–1, respectively. The region upstream of the Australian Lambert Glacier basin (LGB) traverse has a net mass budget of 4.4±6.3 Gt a–1, while the downstream region has –8.9±9.9 Gt a–1. These results indicate that glacier drainages 9, 10 and 11, upstream and downstream of the Australian LGB traverse, are in balance to within our measurement error.