We applied a newly developed three-dimensional, time-dependent, thermomechanical icesheet model including higher-order stress gradients, to simulate the ice flow across subglacial Vostok lake, East Antarctica. Simulations of both ‘lake’ and ‘no lake’ conditions (by treating the ice/lake interface as a stress-free surface, similar to an ice shelf) demonstrate the effects of the subglacial lake, such as pronounced surface flattening and ice-flow turning, on the overall ice dynamics in the vicinity, although subglacial lake dynamics are not treated explicitly. When buoyancy forces and hydrostatic equilibrium of the ice above the lake are taken into account, the along-lake surface slope is preserved and the ice-flow pattern is in accord with sparse observations. Model experiments point to a local ice speed-up in the northern part of the lake, which can be associated with the onset of an enhanced ice-flow feature, more precisely the onset of the Totten Glacier catchment.

Recent investigations of crevassing on alpine glaciers and ice shelves have been based on linear elastic fracture mechanics (LEFM). However, LEFM is unable to explain some aspects of crevasse formation such as the initiation of crevasse propagation from crystal-scale (mm) microcracks, the slow propagation of large fractures in ice shelves, and the acceleration of crevasse opening before breaking of the ice terminus. Here another mechanism to account for these observations is proposed: subcritical crevassing. Subcritical crack growth, documented in many materials though not yet explored in ice, is characterized by a crack velocity that scales as a power of the tensile stress intensity factor, but is much less than that associated with critical crack propagation. This mechanism allows crevasse propagation from mm-scale microcracks at velocities much lower than body wave speeds, and explains crevasse-opening accelerations in a natural way. Subcritical crevassing is theoretically explored for several simplified situations but is limited by a lack of available data on crevasse evolution.

The upper 20—30 m of ice-rich permafrost at three sites overridden by the northwest margin of the Laurentide ice sheet in the Tuktoyaktuk Coastlands, western Arctic Canada, comprise massive ice beneath ice-rich diamicton or sandy silt. The diamicton and silt contain (1) truncated ice blocks up to 15 m long, (2) sand lenses and layers, (3) ice veins dipping at 20—30°, (4) ice lenses adjacent and parallel to sedimentary contacts, and (5) ice wedges. The massive ice is interpreted as intrasedimental or buried basal glacier ice, and the diamicton and silt as glacitectonite that has never thawed. Deformation of frozen ground was mainly ductile in character. Deformation was accompanied by sub-marginal erosion of permafrost, which formed an angular unconformity along the top of the massive ice and supplied ice clasts and sand bodies to the overlying glacitectonite. After deformation and erosion ceased, postglacial segregated ice and ice- wedge ice developed within the deformed permafrost.

Over the period 1972–98 the height of the snow surface at eight Antarctic sites in Palmer Land and on Alexander Island has been measured with respect to fixed points on local nunataks. From these data an empirical relation between height changes over a given period and three key variables has been derived. These variables are (i) the local mean annual surface air temperature, (ii) a regional estimate of energy available for melt over the period (derived from the nearby Rothera air-temperature record) and (iii) a regional estimate of accumulation over the period (derived from the nearby Gomez Nunatak ice-core accumulation record). Using this relation, the contribution of the Antarctic Peninsula to sea-level rise for warming from climatic conditions (averaged over the last 30 years) is estimated to be −0.006 ± 0.002 mm a−1 K−1. If recent warm conditions persist, however, and meltwater can run off to the sea, the contribution to sea-level rise from ablation is calculated to be 0.07 ± 0.02 mm a−1 K−1.

It has been an underlying assumption in many studies that near-surface layers imaged by ground-penetrating radar (GPR) can be interpreted as depositional markers or isochrones. It has been shown that GPR layers can be approximately reproduced from the measured electrical properties of ice, but these material layers are generally narrower and more closely spaced than can be resolved by typical GPR systems operating in the range 50−400 MHz. Thus GPR layers should be interpreted as interference patterns produced from closely spaced and potentially discontinuous material layers, and should not be assumed to be interpretable as precise markers of isochrones. We present 100 MHz GPR data from Lyddan Ice Rise, Antarctica, in which near-surface (<50 m deep) layers are clearly imaged. The growth of the undulations in these layers with depth is approximately linear, implying that, rather than resulting from a pattern of vertical strain rate, they do correspond to some pattern of snowfall variation. Furthermore, comparison of the GPR layers with snow-stake measurements suggests that around 80% of the rms variability in mean annual accumulation is present in the GPR layers. The observations suggest that, at least in this case, the GPR layers do approximate isochrones, and that patterns of snow accumulation over Lyddan Ice Rise are dominated by extremely persistent spatial variations with only a small residual spatial variability. If this condition is shown to be widely applicable it may reduce the period required for measurements of surface elevation change to be taken as significant indications of mass imbalance.

A scheme for katabatic turbulent heat transfer proposed by Oerlemans and Grisogono (2002), here referred to as the OG parameterization, is compared with bulk heat-transfer estimates on Peyto Glacier, Alberta, Canada. Automatic weather stations (AWSs) provide off-glacier data to drive the parameterization and glacier data for bulk estimates. Micrometeorological datasets are used to assess two schemes that employ the Monin-Obukhov stability parameter, z/L, to modify logarithmic, or neutral, bulk heat-transfer equations to allow for stability. Both schemes fail at >1 m above the surface, where the AWS sensors are located, unless a modified approach is used in which the stability correction is constant for z/L ≥1/3. Then the bulk sensible-heat-flux density falls to ≈0.93 of its neutral estimate at all measurement levels, thus providing a basis for comparison with the parameterization. The results of the comparison are very good, indicating that a one-to-one relationship between bulk and parameterized values can be achieved by optimizing the fit with a background exchange coefficient and, because there is only one off-glacier AWS, using a sinusoidal function to model the diurnal variation of the potential temperature lapse rate.

Investigating snow avalanches using a purely statistical approach raises several issues. First, even in the heavily populated areas of the Alps, there are few data on avalanche motion or extension. Second, most of the field data are related to the point of furthest reach in the avalanche path (run-out distance or altitude). As data of this kind are tightly dependent on the avalanche path profile, it is a priori not permissible to extrapolate the cumulative distribution function fitted to these data without severe restrictions or further assumptions. Using deterministic models is also problematic, as these are not really physically based models. For instance, they do not include all the phenomena occurring in the avalanche movement, and the rheological behaviour of the snow is not known. Consequently, it is not easy to predetermine extreme-event extensions. Here, in order to overcome this problem, we propose to use a conceptual approach. First, using an avalanche-dynamics numerical model, we fitted the model parameters (friction coefficients and the volume of snow involved in the avalanches) to the field data. Then, using these parameters as random variables, we adjusted appropriate statistical distributions. The last steps involved simulating a large number of (fictitious) avalanches using the Monte Carlo approach. Thus, the cumulative distribution function of the run-out distance can be computed over a much broader range than was initially possible with the historical data. In this paper, we develop the proposed method through a complete case study, comparing two different models.

This paper presents two methods for glacier monitoring on Svalbard using synthetic aperture radar (SAR) satellite images. Both methods were developed on glaciers in the Kongsfjorden area. The first method monitors the firn area extent and the firn line over time by thresholding and filtering the SAR image. Manual detection of the threshold is preferable, but using a constant threshold for all images also gives adequate results. A retreat of the firn-line position is visible, especially on Kongsvegen, corresponding to consecutive years of negative mass balance. The second method applies a k-means classification to three clusters on the glacier surface. The areal extent of the resulting class on the upper part of the glacier correlates remarkably well with the independently measured mass balance of Kongsvegen, having a correlation coefficient of around 0.89 for the various glaciers. This is because the snow from the accumulation area influences the k-means classification. Thus, on glaciers where mass-balance values are available, new mass-balance values can be predicted solely from SAR images. For glaciers where no mass balance is available, the area change cannot be calibrated to absolute mass-balance values, but relative changes can be predicted.

Long-term ice-core records of Alpine glaciers are usually taken from cold-firn areas at high altitudes, as on Mont Blanc and Monte Rosa. Perennial ice bodies at lower altitudes, however, also bear information about the past. Recent findings from the remains of such ice (the Oetztal iceman found in Austria; wooden bows at Lötschen Pass, Switzerland) clearly indicate the hitherto little-recognized fact that small, more-or-less static perennial ice bodies which are cold and frozen to the underlying bedrock may contain very old ice and, hence, important palaeoclimatic information about warm periods with minimum ice extent in the Alps. Since autumn 1998, investigations have been initiated on a crest-type location or “miniature ice cap” at Piz Murtèl, Engadine, Swiss Alps. First results from shallow drilling, temperature data-logging, geodetic surveying, visual observation, finite-element modelling of simplified basic two-dimensional configurations and comparison to earlier measurements at similar sites provide promising perspectives concerning a little-studied phenomenon with considerable scientific-environmental research potential. Specific characteristics of the investigated site, and probably of many other comparable mountain sites, are: cold ice (about –4°C at 10 m depth), no basal sliding, small mass turnover, striking lack of a firn zone, accumulation mainly by superimposed ice, and direct access to old layers (centuries, millennia?) at the ice/bedrock interface.

This paper presents an overview of firn accumulation in Dronning Maud Land (DML), Antarctica, over the past 1000 years. It is based on a chronology established with dated volcanogenic horizons detected by dielectric profiling of six medium-length firn cores. In 1998 the British Antarctic Survey retrieved a medium-length firn core from western DML. During the Nordic EPICA (European Project for Ice Coring in Antarctica) traverse of 2000/01, a 160 m long firn core was drilled in eastern DML. Together with previously published data from four other medium-length ice cores from the area, these cores yield 50 possible volcanogenic horizons. All six firn cores cover a mutual time record until the 29th eruption. This overlapping period represents a period of approximately 1000 years, with mean values ranging between 43 and 71 mm w.e. The cores revealed no significant trend in snow accumulation. Running averages over 50 years, averaged over the six cores, indicate temporal variations of5%. All cores display evidence of a minimum in the mean annual firn accumulation rate around AD 1500 and maxima around AD 1400 and 1800. The mean increase over the early 20th century was the strongest increase, but the absolute accumulation rate was not much higher than around AD 1400. In eastern DML a 13% increase is observed for the second half of the 20th century.

The melt cycle of snow is investigated by combining ground-based microwave radiometric measurements with conventional and meteorological data and by using a hydrological snow model. Measurements at 2000 m a.s.l in the basin of the Cor-devole river in the eastern Italian Alps confirm the high sensitivity of microwave emission at 19 and 37 GHz to the snow melt−freeze cycle, while the brightness at 6.8 GHz is mostly related to underlying soil. Simulations of snowpack changes performed by means of hydrological and electromagnetic models, driven with meteorological and snow data, provide additional insight into these processes and contribute to the interpretation of the experimental data.

Global positioning system (GPS) data are now routinely used for many glaciological applications. In some common cases, systematic errors are unmodelled at the data-processing stage, although they are often presumed insignificant. In this paper, I investigate these assumptions for three different scenarios: (1) measurements on a moving glacier; (2) measurements on a floating ice shelf; and (3) precise height determination over large elevation ranges, such as for aircraft positioning in lidar/laser altimeter missions. In each case, systematic errors are shown to be present in the coordinate solutions that have a far greater magnitude than the formal error estimates produced by the GPS processing software, under certain conditions. If these coordinate biases go undetected, short- and long-term measurements of horizontal ice velocity or rates of ice-thickness change may be erroneous and the coordinates could not be expected to match rigorously processed data or results from different processing techniques. More rigorous processing strategies are discussed that allow for bias-free parameter estimation.

In the framework of the European Project for Ice Coring in Antarctica (EPICA), a comprehensive glaciological pre-site survey has been carried out on Amundsenisen, Dronning Maud Land, East Antarctica, in the past decade. Within this survey, four intermediate-depth ice cores and 13 snow pits were analyzed for their ionic composition and interpreted with respect to the spatial and temporal variability of volcanic sulphate deposition. The comparison of the non-sea-salt (nss)-sulphate peaks that are related to the well-known eruptions of Pinatubo and Cerro Hudson in AD 1991 revealed sulphate depositions of comparable size (15.8±3.4 kg km–2) in 11 snow pits. There is a tendency to higher annual concentrations for smaller snow-accumulation rates. The combination of seasonal sodium and annually resolved nss-sulphate records allowed the establishment of a time-scale derived by annual-layer counting over the last 2000 years and thus a detailed chronology of annual volcanic sulphate deposition. Using a robust outlier detection algorithm, 49 volcanic eruptions were identified between AD 165 and 1997. The dating uncertainty is ±3 years between AD 1997 and 1601, around ±5 years between AD 1601 and 1257, and increasing to ±24 years at AD 165, improving the accuracy of the volcanic chronology during the penultimate millennium considerably.

A heavy down-hole hammer actuated from the surface by a light composition rope was used to place instrumented probes into the active, 7m thick, clast-rich till underlying a site on Black Rapids Glacier, Alaska, USA, where the ice is 500m thick. A till penetration of about 2.5m was obtained, and greater depths seem possible. The probes measured pore-water pressure and two axes of tilt, which they broadcasted, without wires, to a receiver just above the ice–till interface.

The behavior of majorions, δ18O, dissolved organic carbon (DOC) and trace elements was studied during the initial freezing of low-saline water (3 practical salinity units) in a freezing experiment. Samples were also collected from first-year sea ice from pack ice in the Bothnian Bay, northern Baltic Sea. During initial ice formation, the major-ion ratios in sea ice indicated variable behavior, with some ions showing relative enrichment (sulfate, calcium and magnesium), conservative behavior (sodium) or relative depletion (potassium) compared to sea water at the same salinity. DOC, iron and aluminum showed enrichment in the ice, while zinc was depleted to salinity. Lead was detected in surface snow-ice layers only, implying atmospheric accumulation. First-year sea ice, with a variable growth and thermal history, showed behavior for major ions similar to that observed in new ice. However, for trace elements the picture was much more complicated, most likely due to active secondary processes such as atmospheric supply and biological activity. Ice growth has a potential impact on the chemical budgets and cycling of some elements, especially those which are selectively rejected/retained during sea-ice formation, particularly in the shallow parts of the Bothnian Bay covered with a land-fast ice cover.

Mertz Glacier, East Antarctica, is characterized by a 140 km long, 25 km wide floating ice tongue. In this paper, we combine a large number of remotely sensed datasets, including in situ global positioning system measurements, satellite radar altimetry, airborne radio-echo sounding and satellite synthetic aperture radar imagery and interferometry. These various datasets allow us to study the interaction of the ice tongue with the tides and currents. However, the inverse barometer effect needs to be applied to sea-level variations affecting the tongue. We find that the tide-induced currents exert a small lateral pressure on the tongue which, when integrated over the large surface of the tongue, induce a flexure of up to 2 m amplitude per day. Simple elastic modelling of the flexure confirms the observations and helps validate the boundary conditions necessary to explain different eastward and westward tongue deflections. In addition, the along-flow velocity of the tongue does vary daily from 1.9 to 6.8 md-1 depending on the tidal current. When the current pushes the tongue toward the eastern boundary of the fjord, the tongue is retarded by the drag and the velocity decreases. The accumulated stress is released, allowing the tongue to flow very rapidly when the current pushes the tongue westward. These forcing and boundary conditions on the floating ice flow are important and must be taken into account when studying glacier discharge and calving.

An instrument is described that automatically records ice ablation while eliminating the need for ablation stakes. A pressure transducer placed at the bottom of a hole drilled into the ice is connected by a hose to a bladder lying on the surface. Ice ablation is detected as a reduction in the hydrostatic pressure measured by the transducer.

We have developed a novel ice-penetrating radar system that can be carried on a backpack. Including batteries for a 3 hour continuous measurement, the total weight is 13 kg. In addition, it operates reliably down to –25°C, has a low power consumption of 24 W, and is semi-waterproof. The system has a built-in-one controller with a high-brightness display for reading data quickly, a receiver with 12-bit digitizing, and a 1 kV pulse transmitter in which the pulse amplitude varies by <0.2%. Optical communications between components provides low-noise data acquisition and allows synchronizing of the pulse transmission with sampling. Measurements with the system revealed the 300 m deep bed topography of a temperate valley glacier in the late ablation season.

During two cruises in 1998 and 1999, we examined drift and ridging characteristics of sea ice in the Ross Sea, Antarctica. Mean ice thickness in the western Ross Sea in autumn was 0.5 m, while higher level-ice thicknesses, greater areal coverages of ridges and higher sails were found in the central and eastern Ross Sea in summer. Near the continent, ice drifted westward near the coast and turned eastward further north. We use a regional sea-ice−mixed-layer−pycnocline model to initiate backward trajectories at the time and location of field observations and examine the dynamic and thermodynamic processes that determine ice thickness along these trajectories. Model results agree with previously published field data to indicate that thermodynamic and dynamic thickening and snow-ice formation each contribute significantly to the ice mass of the summer ice field in the central and eastern Ross Sea. For first-year ice in the western Ross Sea, model results and field data both indicate that thermodynamic thickening is the dominant process that determines ice thickness, with dynamic thickening also contributing 20% to the net ice-thickening rate. However, model results fail to reproduce the prevalence of snow- ice formation that was seen in field data.