Fast ice flow on the Antarctic continent constitutes much of the mass loss from the ice sheet. However, geophysical methods struggle to constrain ice flow history at depth, or separate the signatures of topography, ice dynamics and basal conditions on layer structure. We develop and demonstrate a methodology to compare layer signatures in multiple airborne radar transects in order to characterize ice flow at depth, or improve coverage of existing radar surveys. We apply this technique to generate synthetic, along-flow radargrams and compare different deformation regimes to observed radargram structure. Specifically, we investigate flow around the central sticky spot of Whillans Ice Stream, West Antarctica. Our study suggests that present-day velocity flowlines are insufficient to characterize flow at depth as expressed in layer geometry, and streaklines provide a better characterization of flow around a basal sticky spot. For Whillans Ice Stream, this suggests that ice flow wraps around the central sticky spot, supported by idealized flow simulations. While tracking isochrone translation and rotation across survey lines is complex, we demonstrate that our approach to combine radargram interpretation and modeling can reveal critical details of past ice flow.

To exploit the high accuracy potential of laser ranging, suitable laser-altimeter calibration procedures combining laboratory and in-flight calibration steps are needed. Previous calibration schemes have limited the calibration process to flying over planar surfaces, and systematic errors remained in the system. In this paper we present a laser calibration algorithm that utilizes natural surfaces as calibration sites. This is especially important for calibrating systems in the interior of ice sheets. The proposed approach is demonstrated in calibrating the U.S. National Science Foundation Support Office for Aerogeophysical Research laser altimetry system. The results show a significant improvement in accuracy.

In this work, we study the size of the particles involved in a powder-snow avalanche. To determine this value, we study all the phenomena encountered by the particles before they arrive in the “body” of the avalanche. We study the boundary layer which is at the bottom of the avalanche. We determine, with the help of experimental data, the range of sizes of particles that can be entrained by the avalanche. We then examine the possibility of these particles reaching the top of the boundary layer and thus taking part in the avalanche. Our final result is that the typical size of particles suspended in a powder-snow avalanche is < 200 μm.

The Woore Moraine, situated within the Cheshire Plain, is one of the largest glacial landforms in Britain. However, its key characteristics, internal structure and glaciological significance are poorly understood due to limited exposure of its key geological features. This paper focuses on the results of a suite of complementary geophysical techniques (ground-penetrating radar, electrical resistivity tomography and seismic refraction) that were used to image the interior structure of the moraine and help determine its origin. Geophysical imaging reveals a two-layer moraine structure comprising a thin and heterogeneous carapace of intercalated sediments overlying a core of subglacial clay-rich diamicton and proglacial sands deformed by a series of fore- and back-thrusts. These features indicate that the moraine is glaciotectonic in origin, involving initial compression of proglacial sediments against a bedrock obstruction and subsequent modification due to overriding.

The three-parameter hydraulic model of snow avalanche dynamics including the coefficients of dry and turbulent friction and the coefficient of new-snow-mass entrainment was investigated. The ‘Domestic’ avalanche site in Elbrus region, Caucasus, Russia, was chosen as the model avalanche range. According to the model, the fixed avalanche run-out can be achieved with various combinations of model parameters. At the fixed value of the coefficient of entrainment me , we have a curve on a plane of the coefficients of dry and turbulent friction. It was found that the family of curves ( me is a parameter) are crossed at the single point. The value of the coefficient of turbulent friction at the cross-point remained practically constant for the maximum and average avalanche run-outs. The conclusions obtained are confirmed by the results of modelling for six arbitrarily chosen avalanche sites: three in the Khibiny mountains, Kola Peninsula, Russia, two in the Elbrus region and one idealized site with an exponential longitudinal profile. The dependences of run-out on the coefficient of dry friction are constructed for all the investigated avalanche sites. The results are important for the statistical simulation of avalanche dynamics since they suggest the possibility of using only one random model parameter, namely, the coefficient of dry friction, in the model. The histograms and distribution functions of the coefficient of dry friction are constructed and presented for avalanche sites Nos 22 and 43 (Khibiny mountains) and ‘Domestic’, with the available series of field data.

Relatively little data on the distribution of Antarctic icebergs were available prior to 1980. The published literature included size data of about 5000 icebergs, and position data of 12 000 icebergs. There were indications that the size data were biased in favour of larger icebergs.

A programme of systematic iceberg observations was therefore initiated by Norsk Polarinstitutt in 198! through the SCAR Working Group on Glaciology. This programme is based on standard “blue” forms distributed to all ships going to Antarctica. The icebergs are recorded every 6 h and in Five length groups: 10–50, 50–200, 200–500, and 500–1000 m, and those over 1000 m are described individually.

The amount of data has increased greatly from the start in 1981–82. The position of 70 000 icebergs, including 50 000 that had been size classified, were on file at Norsk Polarinstitutt by December 1985, and the data set is growing rapidly. Most ships travelling to and from Antarctica now participate in collection of the data. (Fig.1 shows the locations of the icebergs sighted.) Fig. 1.

Location of iceberg observations under the programme initiated in 1981. Main ship tracks are clearly reflected. The average observation represents 14 icebergs.

The size distribution of the classified icebergs observed under this programme up to December 1985 is given in Table I: Table I

The “standard size” (length, width, and thickness) is based on our observations from three Antarctic expeditions which carried out dedicated iceberg studies. Many icebergs are of course not right-angled parallelepipedal in shape, but this is a good approximation for most of the larger icebergs.

The data are based both on visual sightings and on radar observations. Duplicate observations from a ship moving at slow or zero speed are as far as possible eliminated, both during observation, and by critical appraisal before the data are filed. The data editing also includes evaluation of data quality, especially in connection with radar observations, and comparison of positions and dimensions of the large icebergs in order to reduce to a minimum repeated observations from different vessels of icebergs >1000 m. These account for most of the iceberg mass (see Table I).

Consideration of iceberg-distribution patterns and the observed area of the Southern Ocean, and of duplicate observations, indicates more than 300 000 icebergs south of the Antarctic Convergence, with a total ice mass of about 1016 kg. Consideration of mean residence times indicates an annual iceberg production from the continent of 23–1015 kg, which is considerably higher than most other recent estimates. This also suggests that the Antarctic ice sheet is in balance.

The data indicate large regional differences in iceberg sizes, the most noticeable being between the two sides of the Antarctic Peninsula, and between the Amery Ice Shelf/ Prydz Bay area and the remainder of East Antarctica. These differences are probably mainly related to different calving sites.

About one-third of the observed icebergs are over the continental shelf of Antarctica. The total under-water area of these icebergs is two orders of magnitude less than the under-water area of the Antarctic ice shelves. The annual total iceberg melting and its effect on the water masses over the continental shelf has been calculated from ocean-water temperature variations at 200 m depth and estimated melt rates. This turns out to be an order of magnitude less than the annual effect of melting sea ice. The iceberg data considered here are probably under-represented with respect to the smallest sizes, and they do not include icebergs that have become <10 m. Inclusion of these ice bodies would increase the total melt.

GPS measurements of tidal modulation of ice flow and seismicity within the grounding zone of Beardmore Glacier show that tidally induced fluctuations of horizontal flow are largest near the grounding line and decrease downstream. Seismic activity is continuous, but peaks occur on falling and rising tides. Beamforming methods reveal that most seismic events originate from two distinct locations, one on the grid-north side of the grounding zone, and one on the grid-south side. The broad pattern of deformation generated as Beardmore Glacier merges with the Ross Ice Shelf results in net extension along the grid-north side of the grounding zone and net compression along the grid-south side. During falling tides, seismic activity peaks on both sides because of increased vertical flexure across the grounding line. During rising tides, seismic activity in the region of extension on the grid-north side is relatively low because the tidal influence on both horizontal strain rate and vertical flexure is small. On the grid-south side during rising tides, however, tidally induced horizontal strain rates promote increased seismicity in regions of long-term compressional flow paths. Our study highlights how concurrent geodetic and seismic measurements provide insight into grounding-zone mechanics and their influence on ice-shelf buttressing.

Radio-echo sounding techniques are very useful for fast profiling of seasonal snowpack. Ground-penetrating radar (GPR) is used widely for various cryospheric applications, such as snow/glacier depth estimation, snow layer identification and snow water equivalent assessment. The dielectric constant of snow is an important input parameter for the acquisition and interpretation of GPR data from the snowpack. In this study, snow dielectric constant was measured along with physical properties of snow using a snow fork operating at 1 GHz frequency. Experiments were conducted at field observatories of the Snow and Avalanche Study Establishment located in different Himalayan ranges: Patseo (Greater Himalayan range), Dhundhi and Solang (Pir Panjal range). Interseasonal spatial and temporal variations in snow dielectric constant and associated snowpack properties were analysed for five winter seasons (2010-14). The mean seasonal snow dielectric constant is higher at Dhundhi (1.82 ±0.02) than at Patseo (1.69 ±0.02). The measured snow dielectric constant was used to derive snow density and liquid-water content (LWC). A better correlation between snow dielectric constant and LWC is observed for high-density snow (>300kgm-3; R2 = 0.95) than for low-density snow (<200kgm-3; R2 = 0.73). Snow-fork-derived snow density was in good agreement with manually measured values. The snow dielectric constant database generated during this study can be used as a reference for various field applications of GPR in snow-related studies.

Four different algorithms for retrieving ice concentration from passive microwave imagery are applied to SSM/I data collected over the Weddell Sea in September 1989. Comparison of the results along a typical transect from the ice edge to the coast shows point-wise differences of up to 45% in ice concentration. The observed differences can largely be explained by the different combinations of the 19 and 37 GHz, horizontal and vertical polarization, data channels utilized by each algorithm. Through their frequency/polarization signatures, and with the aid of coincident surface observations, the differences are further related to ice type, specifically the presence of grease ice or thin ice, and surface conditions.

A reconnaissance program has been carried out to identify problems caused by glaciers in a large proposed hydroelectric development in the Susitna River basin of Alaska. Balance measurements on the major glaciers have been initiated, and long-term balance between 1949 and 1980 has been estimated from existing photo sets. From the latter it appears that shrinking of the glac!iers, which comprise 4% of the basin area, may have contributed appreciably to the measured basin runoff. A potential instability in the drainage of Eureka Glacier, on the edge of the basin, has been identified. The glaciers of the basin seem to be largely temperate, and most of them are surging or pulsing types. Velocity measurements show seasonal variations that suggest appreciable contribution to the motion from basal sliding. A study of the moraines of Susitna Glacier, which is a surging type, indicates that no surge is imminent. Glacier-dammed lakes exist in the basin; they are small but could be enlarged by surging or other mechanisms. Some general problems in the estimation of the transport of suspended sediment are noted.

The repetitive synoptic ice data obtainable by satellite sensing provide a means of studying the time-dependent behavior of both sea ice and ice sheets on climatic time scales. Examples of sea-ice parameters which may be measured are extent, concentration, and multiyear fraction; and examples of ice-sheet/ice-shelf parameters are surface elevation, ice-front position, extent and duration of summer melting, and ice accumulation rates. Desired snow-cover parameters include extent and snow depth or water-equivalent depth. The unique ability of satellites to measure such ice parameters and the characteristics of the consequent data sets significantly influence the structure of ice models that can be successfully used with the data. Ice data sets recently acquired by satellite sensing are described. The past decade of sea-ice data provides a detailed description of the interannual variability of sea ice on a regional and seasonal basis. Because of the longer time scales involved in ice-sheet variations, a comparable record of ongoing ice-sheet variations has not yet been established, but important baseline data sets are being developed.

New categories for the climatic division of snowy areas according to their snow-cover character in mid-winter are proposed. They are a wet-snow region, a dry-snow region, an intermediate snow region and a depth-hoar region. The wet-snow region is defined as the region in which every layer of deposited snow is wet due to percolation of snowmelt water throughout the winter. In contrast, areas in which the snow cover is dry, at least in the coldest period of the winter season, are classified into two categories, that is the dry-snow region and the depth-hoar region. In the latter region, the small snow depth and low air temperature induce development of depth hoar. The intermediate snow region was introduced to indicate an intermediate character between the dry-snow and wet-snow regions. From the climatic dataset calculated by the Japanese Meteorological Agency and from snow surveys, it has been found that in snowy areas, which have a climatic monthly mean temperature in January (Tjan) higher than 0.3°C, snow would be expected to be wet throughout the winter and, in areas that have Tjan, lower than −1.1°C, to be dry at least in the coldest period. Snow covers, where Tjan is between these two values, are expected to have intermediate characters. Therefore, these temperatures are supposed to be critical values among the wet, dry and intermediate snow regions. The criterion that separates the depth-hoar region from the dry-snow areas was found to be given by a climatic mean temperature gradient. This value lies between 10 and 12°Cm−1, which is derived by dividing the absolute value of the average of the climatic monthly mean air temperature, which is always below 0°C, by the average of the monthly maximum snow depth during January and February.

Accurate co-location and geo-registration of AVHRR (Advanced Very High Resolution Radiometer) is necessary for most applications involving these data. To obtain the necessary accuracy, most of the current systems available for AVHRR image rectification require corrections to errors arising from several different sources. This paper presents a software package which is able to correct registration errors from all sources using effective roll, pitch, and yaw values determined from ground control point offsets. These effective values are used to establish a baseline registration against which other models are compared. Spacecraft and sampling geometry are revisited in order to systemize and better define the procedure of AVHRR geo-registration. This systematic approach has resulted in a modular software package which has a registration accuracy of 1 satellite pixel, and is easily modified and improved.

Simple models to calculate melt and refreezing are reviewed. Both degree-day and energy-balance models can give distributed melt inputs to ice-dynamics models but have only been tested extensively in West Greenland, and more data are needed from the remoter parts of Greenland. The energy-balance model is more realistic but needs input data that are not generally available over the whole ice sheet. On the other hand, degree-day factors vary from situation to situation although a value of about 8 kg m−2 d−1 deg−1 for ice ablation is a reasonable first approximation as assumed in recent ice-dynamics models. Meltwater refreezing in the accumulation area can be modelled very simply but more detailed physical models are needed to describe the shifts in accumulation zones under different climates.

For the first time, a cold ice cover of a summit in the central Swiss Alps has been sampled from the surface to the bed for determining its isotopic composition in δD and δ18O. Results of the analyses show a striking decrease of δ-values with depth. The δ-value differences are greater than those explicable by a direct temperature effect, but a substantial increase in melt water percolation through the firn since the formation of the deepest ice layer may explain the situation of this high-altitude ice.

The time evolution of drifting snow under a steady wind is estimated using a new numerical model of drifting snow. In the model, Lagrangian stochastic theory is used to incorporate the effect of turbulence on the motion of drifting-snow particles. This method enables us to discuss both the saltation and the suspension process. Aerodynamic entrainment, grain/bed collision (splash process), wind modification and particle size distribution are also taken into account. The calculations show that the time needed by the total mass flux to reach a steady state appears to be 3–5 s. Vertical profiles of horizontal mass flux near the surface show a similar tendency. In contrast, it takes >50 s for the wind speed and the whole mass-flux profile to reach a steady state. This longer time depends on the time-scale of the turbulent diffusion, which is responsible for the mass flux extending to an order of a few meters in height. Applying Taylor’s hypothesis, the estimated length scale at which drifting snow reaches equilibrium is around 400 m. This result is comparable with previously reported field observations.

Recent observations suggest that high-frequency forcing events have profound influence on the air-sea-ice interactions in the polar region. Studying these events with sea-ice models requires close examination of the model sensitivity that may arise from the high-frequency variability of the forcing. We show that the maximum layer thickness is dictated by the time-scale of the forcing variability and that the computation of the surface temperature develops enhanced sensitivity at high-frequency forcing. We resolve these constraints by developing an "adjustable-layering" thermodynamic formulation for ice and snow that re-computes the number of layers required each time-step to satisfy this maximum thickness, which preserves the total enthalpy and general internal thermal gradients. The conservation equations form a tri-diagonal system ideal for a fast and efficient implicit solution. Furthermore, we resolve the issue of the high sensitivity of the surface flux balance by solving the linearized version of the flux boundary condition simultaneously with the overall conservation system. In this paper we develop the analyses specifying the model requirements, describe the model system and test its algorithmic implementation.

Very little data exist on the fatigue behavior of ice, yet fatigue may have an important role to play in a number of ice failure situations. Whilst reversed-bending tests and cyclic-compression tests have been performed on ice, no attempts have been made previously to measure fatigue-crack growth rates in ice. This paper presents the results of fatigue-crack growth experiments on S2 (i.e. columnar grained with the c axes randomly aligned within the horizontal plane) fresh-water ice. There is clear evidence that crack growth has occurred under cyclic loading but, as yet, the data are insufficient to determine the validity of da/dΝ versus δK type behavior. Preliminary results are presented and the implications of the results are discussed.

The densification of water-saturated firn, which had formed just above the firn-ice transition in the wet-snow zone of temperate glaciers, was investigated by compression tests under pressures ranging from 0.036 to 0.173 MPa, with special reference to the relationship between densification rate, time and pressure. At each test, the logarithm of the densification rate was proportional to the logarithm of the time, and its proportionality constant increased exponentially with increasing pressure. The time necessary for ice formation in the firn aquifer was calculated using the empirical formula obtained from the tests. Consequently, the necessary time decreased exponentially as the pressure increased, which shows that the transformation from firn in ice can be completed within the period when the firn aquifer exists, if the overburden pressure acting on the water-saturated firn is above 0.12–0.14 MPa. This critical value of pressure was in good agreement with the overburden pressure obtained from depth–density curves of temperate glaciers. It was concluded that the depth of firn–ice transition was self-balanced by the overburden pressure to result in the concentration between 20 and 30 m.

We summarize here data on in situ nitrate ion concentrations in snow pits and firn cores over the last ∼3 250 a. Nitrate fluctuations show seasonal, 11 and 22 a periodicities, and long-term changes both at South Pole station and Vostok. High nitrate levels conform to winter darkness and solar activity peaks. Long-term lows and highs conform to solar activity minima and maxima. The data available support the hypothesis that nitrate is fixed in the upper atmosphere by some solar-mediated phenomenon causing a periodicity in East Antarctica snow. Background levels and non-periodic spikes in nitrate come from other sources.