This paper presents a glacier ice-flow model that simulates changes to alpine glaciers of various sizes and their runoff response to climate change in the Yili river basin in the Tien Shan mountains, northwestern China. It is suggested that the sensitivity of glaciers to climatic change is determined by glacier size. The change in glacial runoff does not keep pace with climatic change. As climate warms and glaciers retreat, the glacier runoff tends to increase and then decrease. The runoff peak and its timing depend not only on glacier size but also on the rate of air-temperature rise.

A laboratory experiment was conducted in which new ice crystals were nucleated from the vapor phase onto large existing ice crystals obtained from Antarctic lake ice. Flat, smooth ice-crystal surfaces were prepared, with c axes oriented either vertically or horizontally. When these were subjected to a supersaturated vapor environment, multiple individual crystals nucleated onto the substrates adopting the same crystallographic orientation as the parent. A dominant grain-growth scenario for kinetic-growth metamorphism in snow, which in some ways is analogous to the oriented morphologies in lake ice, is hypothesized. In the lake-ice-growth scenario, optimally oriented crystals will grow at the expense of those less preferentially positioned.The proposed dominant grain-growth theory for snow is in agreement with the observed decrease in the number of grains and the proximal similarity of crystal habit in kinetic-growth metamorphism in snow. Similarly, kinetic crystal growth on the interior of gas inclusions in Antarctic lake ice will also acquire the crystallographic orientation of the substrate ice. These small-faceted interior crystals significantly influence light scattering and penetration in the lake-ice cover.

A simple approach to glacier dynamics is explored in which there is postulated to be a relationship between area and volume with three parameters: the time for area to respond to changes in volume, a thickness scale, and an area characterizing the condition of the initial state.This approach gives a good fit to the measurements of cumulative balance and area on South Cascade Glacier from 1970–97; the area time-scale is roughly 8 years, the thickness scale about 123 m, and the 1970 area roughly 4% larger than required for adjustment with volume. Combining this relationship with a version of mass continuity expressed in terms of area and volume produces a theory of glacier area and volume response to climate in which another time constant, the volume time-scale, appears. Area and volume both respond like a damped spring and mass system. The damping of the South Cascade response is approximately critical, and the volume time-scale is roughly 48 years, six times the area time-scale. The critically damped spring and mass analogy reproduces the time dependence predicted by the more complicated traditional theory of Nye.

Two disintegration events leading to the loss of Larsen A and B ice shelves, Antarctic Peninsula, in 1995 and 2002, respectively, proceeded with extreme rapidity (order of several days) and reduced an extensive, seemingly integrated ice shelf to ajumble of small fragments. These events strongly correlate with warming regional climate and accumulation of surface meltwater, supporting a hypothesis that meltwater-induced propagation of pre-existing surface crevasses may have initiated ice-shelf fragmentation.We address here an additional, subsequent mechanism that may sustain and accelerate the ice-shelf break-up once it begins. The proposed mechanism involves the coherent capsize of narrow (less than thickness) ice-shelf fragments by rolling 90° in a direction toward, or away from, the ice front. Fragment capsize liberates gravitational potential energy, forces open ice-shelf rifts and contributes to further fragmentation of the surrounding ice shelf.

A statistical analysis of the lengths of grain-boundary and transgranular cracks induced during the initial straining of columnar-grain ice by a compressive load applied perpendicular to the long direction of the columns is presented. The analysis shows that the crack lengths are randomly distributed and form distinct but correlated populations.The lognormal distribution function is shown to be a good descriptor of the populations for 5–90% of their range. Statistical models are presented for the lognormal behaviour of the crack-length distribution and for the strain dependence of the crack density. The models assume that a change in the value of the random variable of the respective population depends on the population value of the variable at the time of the change. It is shown that the model for the strain dependence of the crack density is suitable for the strain dependence of the acoustic emission, measured in both columnar-grain and granular ice subject to constant compressive loads. Evidence is also presented for a lognormal dependence of the dislocation density on strain. The analysis demonstrates that the cracks that form during the initial straining of polycrystalline ice are independent, random events and that the resulting crack populations are precursors to failure by fracture.

Taku Glacier is one of the few glaciers in Alaska, U.S.A., that has advanced over the last century: 7 km since 1890. This advance slowed substantially during the past decade, but in summer 2001 the glacier terminus began to readvance at a rate of 30 cm d−1. The advance produced dramatic proglacial sediment deformation up to 200 m in front of the terminus. Two to three large bulges and several secondary bulges developed in the proglacial sediments as a result of glacial compression along a 1 km wide portion of the terminus. The bulge nearest the terminus was 10 m high and 65 m wide. The middle bulge (7 m high) advanced at 15 cm d−1 and the distal bulge (3 m high and 50 m wide) at 9 cm d−1. Crenulations and prominent fractures developed in the overlying vegetation layer. The frontal lobes of the bulges were steep and overlaid a shear zone, where sediments were being thrust up and over the ground surface. Ice-proximal push moraines, 1–10 m high, formed along much of the 9 km wide terminus, although deformation was minimal at some locations.

Extensive albedo data from Storglaciären, Sweden, during nine summers are analyzed, focusing on the effect of surface slope on measurements and on the influence of clouds on albedo of both snow and ice surfaces. On clear-sky days, albedo continuously dropped throughout the day by > 0.3 when derived from measurements in a horizontal plane over the slightly sloping surface. When we correct for the tilt effect, over frozen surfaces the apparent decrease in albedo largely disappeared, while over melting surfaces it became significantly less pronounced. This emphasizes the need to correct horizontally measured data for slope when they are used, for example, as ground truth for satellite-derived albedo, or in high-resolution glacier melt runoff models. Albedo varied by > 0.1 due to cloud fluctuations. To quantify cloud effects, albedo changes between successive (half-)hourly and daily data were regressed against corresponding changes in cloudiness expressed by the ratio of global radiation to top-of-atmosphere solar radiation. The resulting relationship over snow explains 66% of oberved albedo changes using hourly, and 42% using daily, time-steps. Over ice, the corresponding values are 28% and 0%. Hence, cloudiness is a dominating factor in determining short-term variability of albedo over snow, whereas over ice the observed large day-to-day variability ( > 20%) is mainly attributed to changes in the physical characteristics of the weathering crust. Results suggest that cloud effects should be incorporated into models treating snow and ice surfaces separately.

The genesis of some types of patterned ground, including hummocks, frost boils and sorted stone circles, has been attributed to differential frost heave (DFH). However, a theoretical model that adequately describes DFH has yet to be developed and validated. In this paper, we present a mathematical model for the initiation of DFH, and discuss how variations in physical (i.e. soil/vegetation properties) and environmental (i.e. ground/air temperatures) properties affect its occurrence and length scale. Using the Fowler and Krantz multidimensional frost-heave equations, a linear stability analysis anda quasi-steady-state real-time analysis are performed. Results indicate that the following conditions positively affect the spontaneous initiation of DFH: silty soil, small Young’s modulus, small non-uniform surface heat transfer or cold uniform surface temperatures, and small freezing depths. The initiating mechanism for DFH is multidimensional heat transfer within the freezing soil. Numerical integration of the linear growth rates indicates that expression of surface patterns can become evident on the 10–100 year time-scale.

A deformation-controlled triaxial apparatus has been developed in order to investigate the mechanical behaviour of alpine snow. Cylindrical specimens (58 mm in diameter, 126 mm in height) can be axially deformed by applying tensile or compressive strain rates between 2.4 × 10−7 s−1 and 2.4 × 10−2 s−1. Compressed air is used to laterally load the snow specimens with pressures of up to 40 k Pa.The volumetric deformation is found by carefully encapsulating the snow samples in an airtight Mylar foil and measuring the expelled pore air during sample deformation. The multidimensional stress–strain behaviour can therefore be determined, and constitutive theories of snow advanced. The device is located in the cold rooms of the Weissfluhjoch research station of the Swiss Federal Institute for Snow and Avalanche Research. Some 180 tests have been performed over the last four winters. The density of the snow specimens varied between 190 and 435 kg m−3. Part I of this two-part paper discusses the development of the apparatus and the applied testing methodology. Qualitative results are provided. In Part II, a detailed analysis of the results is presented and the viscoelastic properties of snow quantified.

Fine-grained, dry snow with a density range of 190–435 kg m−3 was tested in triaxial compression at −12°C with confining pressures varying between 0 and 40 kPa. The tests were strain-rate controlled, with strain rates ranging between 7.4 × 10−7 s−1 and 6.6 × 10−5 s−1. The analysis of the test results revealed that the relationship between yield stress and viscous strain rate is best given by a power law, similar to polycrystalline ice. However, the power-law exponent n is a function of density and varies between 1.8 (low-density snow, ρ < 200 kg m−3) and 3.6 (high-density snow, ρ > 320 kg m−3). The tests also showed that lower-density snow displays a significant non-linear stress–strain response before yielding. Two further aspects of the constitutive behaviour of snow were identified: (1) the strainrate independence of the post-yield work-hardening behaviour in compression and (2) the independence of the axial yield stress in relation to the confining pressure. The experimental observations are discussed with respect to the mechanical properties of polycrystalline ice, which is the constituent material of the load-bearing ice skeleton.

Twice-daily satellite observations from the Special Sensor Microwave Imager (SSM/I) indicate melt onset and refreeze on southeast-Alaskan icefields. Melt and refreeze are based on 37 GHz vertically polarized brightness temperatures (Tb) and diurnal-amplitude variations (DAV). Two types of melt regime have different summer characteristics. Onset is characterized by increasing average daily Tb and a switch from low- to high-amplitude DAV. Melt timing, calibrated using Juneau Icefield temperatures, correlates well with nearby stream hydrographs. Some pixels maintain high Tb throughout the melt season and return to low-amplitude DAV after melt onset. Refreeze on these pixels is identified by decrease in Tb and accompanying high-amplitude DAV. Other pixels maintain high DAV throughout the summer, indicating nocturnal refreeze. Fall refreeze is determined by the end of high-amplitude DAV. Interannual variability in melt timing and ablation-season length is high. Melt onset and refreeze timing show a regional tendency toward earlier glacier-melt onset and longer ablation seasons from 1988–98.

Based on aerial photographs, topographical maps and the Landsat-5 image data, we have analyzed fluctuations of glaciers in the western Qilian Shan, north-west China, from the Little Ice Age (LIA) to 1990. The areas and volumes of glaciers in the whole considered region decreased 15% and 18%, respectively, from the LIA maximum to 1956. This trend of glacier shrinkage continued and accelerated between 1956 and 1990. These latest decreases in area and volume were about 10% in 34 years. The recent shrinkage may be due either to a combination of higher temperatures and lower precipitation during the period 1956–66, or to continuous warming in the high glacierized mountains from 1956 to 1990. As a consequence, glacier runoff from ice wastage between 1956 and 1990 has increased river runoff by 6.2 km3 in the four river basins under consideration. Besides, the equilibrium-line altitude (ELA) rise estimated from the mean terminus retreat of small glaciers <1 km long is 46 m, which corresponds to a 0.3°C increase of mean temperatures in warm seasons from the LIA to the 1950s.

Marine-geological and -geophysical data collected from the continental shelf in Pine Island Bay, Antarctica, reveal a complex paleo-subglacial drainage system controlled by bedrock topography and subglacial meltwater discharge. Significant amounts of freely flowing meltwater existed beneath former ice sheets in Pine Island Bay. Subglacial drainage is characterized by descriptions of glacial landforms imaged on the sea floor and sedimentary deposits collected in piston cores. Bedrock geology is characterized using seismic data. Large-scale landforms on the shelf include channels and cavities incised into impermeable crystalline bedrock. There is a transition from randomly oriented channels on the inner shelf to a dendritic pattern of elongate channels on the middle shelf. On the outer shelf, a change in basal conditions occurs where sedimentary deposits bury crystalline bedrock. No evidence for flowing meltwater exists on sedimentary substrates. Instead, meltwater formed at the ice–sediment contact was incorporated into the sediments, contributing to development of a deforming bed, which was sampled in piston cores. Characterization of subglacial meltwater processes that occurred in the past may aid in understanding the role meltwater plays in stability of the West Antarctic ice sheet today.

The structural and stable-oxygen-isotopic properties of sea-ice samples collected on land-fast sea ice in the northern Baltic Sea, at 60.2–65.7° N, 21.3–26.9° E, were used to estimate the contribution of snow ice to the total sea-ice thickness, and the fraction of snow in the snow-ice layers and in the total sea-ice thickness. The textural analysis of the sea-ice cores revealed a general structure of a granular surface layer and a columnar basal layer. Using different isotopic criteria to distinguish snow-ice layers, we estimate that snow-ice layers contributed on average 24–32% of the total ice thickness in the study region. Transition ice (intermediate granular–columnar ice) also contributed significantly to the total ice thickness. The snow in the snow-ice layers contributed on average 18–21% of the total sea-ice thickness (by mass), a significant contribution when compared to other regions with sea-ice cover. The methods used could be improved to give more accurate estimates, especially by improving knowledge of the temporal variability of the isotopic properties of the snow incorporated into the sea-ice cover. However, our results can be regarded as plausible estimates for the contribution of snow ice and snow to sea-ice growth in the study area.

The response of Arctic ice masses to climate change is studied using ice cores containing information on past climatic and environmental features. Interpretation of this information requires accurate chronological data. Absolute dating of ice cores from sub-polar Arctic glaciers is possible using well-known radioactive layers deposited by atmospheric nuclear tests (maximum fallout in 1963) and the Chernobyl accident (1986). Analysis of several isotopes (3H, 137Cs) shows that 3H provides the most accurate dating of the 1963 maximum, as indicated also in comparison with results from total-beta measurements (90Sr and 137Cs). Mean annual net mass balances are derived from the dated ice cores from 1963 up to the date of the drillings. The 137Cs and 3H deposited by nuclear tests, after decay correction, are used to define a melt index for all 13 ice cores studied. The relative strength of melting and percolation post-depositional processes is studied on the basis of these 137Cs and 3H deposits.