An algorithm is developed for adjusting glacier surface-velocity vectors, given on the nodes of a square grid, so that they obey a central-difference approximation of the continuity equation. Also required on the grid nodes are the glacier thickness, the ratio of the surface-velocity to the average velocity in the column, and the difference between the mass balance and the thickness change. All these other variables are assumed to be known exactly, and only the surface-velocity field is adjusted. The result is optimum in the sense that the magnitude of the adjustment is minimized. Either the relative or the absolute adjustment can be minimized, depending on how weights are specified. No restriction is placed on the shape of the solution region, and no boundary condition is required. The algorithm is not iterative. The algorithm first forms a parallel flow field that satisfies the continuity equation, and then uses a stream function to add a divergenceless field to it. The stream function that leads to the minimum velocity adjustment is obtained as four independent, interlacing solutions covering the solution region. For each of the four, a well-conditioned, sparse-matrix system of simultaneous linear equations is solved. A compact, sub-optimum, well-behaved iterative procedure is also developed for transforming part of the velocity adjustment into an adjustment of the thickness field.

Micrographs of the surface of evaporating polycrystalline ice containing impurities are presented. The micrographs show that impurities that enter the ice lattice substitutionally enhance the fine structure on the surface, but impurities rejected by the lattice suppress the fine structure. It is concluded that the observed fine structure is not an artifact of the machine but a real property of evaporating ice dependent upon the defects at the surface.

Mass-balance measurements have been renewed on two small ice caps on north-eastern Ellesmere Island. Original stake networks were established in 1972 and 1976. Since then, both ice caps have experienced significant mass losses averaging –70 to –140 kg m−2 a−1. They have also decreased in area. The equilibrium line in this area has averaged around 1150 m for the last decade or so. The ice caps are remnants of former climatic conditions and are out of equilibrium with contemporary climate.

The interaction of point detects in ice has been neglected for a long time. Experimental data obtained from dielectric measurements on HF-doped crystals stimulated a new evaluation of the possibility of an interaction between Bjerrum defects and ions. In a previous paper it has been shown that this leads us to assume the existence of aggregates of Bjerrum defects and ions. In this paper these aggregates and Bjerrum defects are used to explain the dielectric properties of ice, especially the temperature dependence of the product of the high and low frequency conductivity σ0σ∞.

The interaction of Bjerrum defects and impurity molecules leads to a dependence of the concentration of frenkel pairs on Bjerrum-defect concentration. At HF concentrations above the native Bjerrum-defect concentration the formation of a Frenkel pair is enhanced. This leads to the fast out-diffusion which has been studied in highly doped crystals by means of NMR techniques.

We recently reported on the radio-frequency attenuation length of cold polar ice at Summit Station, Greenland, based on bi-static radar measurements of radio-frequency bedrock echo strengths taken during the summer of 2021. Those data also allow studies of (a) the relative contributions of coherent (such as discrete internal conducting layers with sub-centimeter transverse scale) vs incoherent (e.g. bulk volumetric) scattering, (b) the magnitude of internal layer reflection coefficients, (c) limits on signal propagation velocity asymmetries (‘birefringence’) and (d) limits on signal dispersion in-ice over a bandwidth of ~100 MHz. We find that (1) attenuation lengths approach 1 km in our band, (2) after averaging 10 000 echo triggers, reflected signals observable over the thermal floor (to depths of ~1500 m) are consistent with being entirely coherent, (3) internal layer reflectivities are ≈–60$\to$–70 dB, (4) birefringent effects for vertically propagating signals are smaller by an order of magnitude relative to South Pole and (5) within our experimental limits, glacial ice is non-dispersive over the frequency band relevant for neutrino detection experiments.

Using historical and recent aerial photography and structure from motion (SfM) multiview stereo (MVS) techniques, we reconstruct the 1959 and 2018 ice surface topography and determine the geodetic mass balance of Bowman Glacier, a small mountain glacier on northern Ellesmere Island. This is combined with optical satellite imagery to reconstruct the evolution in extent of the glacier over six decades, and ground-penetrating radar measurements of ice thickness to estimate the remaining ice volume. Between 1959 and 2020, Bowman Glacier lost 78% of its extent (reducing from 2.75 to 0.61 km2), while average annual area loss rates have nearly tripled in the past two decades. Over the 1959–2018 period, glacier-wide ice-thickness change averaged −22.7 ± 4.7 m, corresponding to a mean specific annual mass balance of −347.0 ± 71.4 mm w.e. a−1. Projecting rates of area and volume change into the future indicates that the glacier will likely entirely disappear between 2030 and 2060. This study demonstrates the potential of SfM-MVS processing to generate elevation products from 1950/60s historical aerial photographs, and to extend observations of ice elevation and glacier volume change for the Canadian Arctic, prior to the satellite record.

The orientations of individual crystals within a polycrystalline aggregate subjected to stress have a strong influence on its bulk strain rate and flow behavior. The ability to include the effect of crystal fabric and recrystallization processes in an ice flow law, especially at the bottom of glaciers and ice sheets where temperature is close to the pressure-melting point, is important because the stratigraphy of the ice body may be affected and the paleoclimate reconstruction hampered. We present herein three newly developed deformation apparatuses offering the possibility, from single experiments, of investigating different finite strain stages and their corresponding c-axis fabric and grain texture patterns in various confined, shear flow configurations (simple shear, pure shear and compression/extension bending). The technical set-ups and major advantages compared to classical methods are explained, and results from experiments are discussed in order to illustrate the functioning and purposes of the methods. In all experiments, significant variations in the microstructural development have been observed that reflect the varying orientations of the anisotropy and its relationship to the stress pattern. In monocrystalline ice-bending experiments, the pre-existing c-axis fabric is shown to have a profound influence on the response to stress and possibly to the type of slip system activated.

Stable isotope ratios and salinities of ice samples obtained from a submarine ice cliff at Explorers Cove demonstrate that the upper parts of the ice cliff have frozen directly from sea-water and are an underwater expression of permafrost, whereas the lower parts appear to be partially glacial in origin. These results indicate that there may be ice cores in the moraines of Explorers Cove, in which case the coastline of McMurdo Sound is more extensively ice-cored than previously known.

At lat. 7° S. the fern line was at about 3,500 m. during the Riss glaciation, at 3,650 m. during the Würm maximum and at 3,700 m. during a later stage (Würm II?). The moraines were deposited by equatorial-type glaciers from which there was relatively little melt water.

At lat. 14° S., during the Riss glaciation, an important ice cap covered a 3,900–4,000 m. high plateau, the firn line being at about 3,950 m. The role of melt waters was important in the formation of the moraines, suggesting their deposition from tropical-type glaciers. From the evidence of periglacial gravels it appears that this southern region has not been glaciated since Riss times.

Holes drilled into thin areas of the Brunt Ice Shelf encounter a layer of liquid brine less than 1 m thick approximately at sea-level. Assuming the brine to be moving horizontally, analysis of its effects on thermal equilibrium gives an estimate of steady-state annual brine flow that is in good agreement with the value deduced from a percolation model. The effect of firn density on percolation rates is such that the slope of an active brine layer increases rapidly as ice thickness increases. However, the heat transport model predicts that brine layers are unlikely to be active in both very thick and very thin ice shelves.

The authors describe a new method that allows an accurate, cheap and quick characterization of snow grains from natural snowpacks. The method is based on stopping all metamorphism of snow grains by immersing them in cold isooctane. Snow grains are then carried from the field, stored in a cold room and easily observed through a microscope. It then makes possible recording of snow-grain pictures in a state similar to their natural state, even after storage of several months.

The 2037 m deep bore hole at Dye 3 in south Greenland was surveyed in 1981, 1983, 1985, and 1986. The directional surveys show the ice flow is planar with a surface velocity of 12.2m/year at an azimuth of 060°, which agrees with surface velocity measured by navigation satellites. Measurements of hole diameter and inclination are highly correlated with dust content in the ice. The temperature measurements show strong convection in the hole fluid with a cell height of about 20 m and an amplitude of 0.1 K. The calculated mean in-situ ice density is 921.3 ± 1.5kg/m3. Due to ice deformation, the lowest 4 m of the hole were not accessible in 1985 and the lowest 180 m were not accessible in 1986.

Mountain glaciers have response times that govern retreat due to anthropogenic climate change. We use geometric attributes to estimate individual response times for 383 glaciers in the Cascade mountain range of Washington State, USA. Approximately 90% of estimated response times are between 10 and 60 years, with many large glaciers on the short end of this distribution. A simple model of glacier dynamics shows that this range of response times entails consequential differences in recent and ongoing glacier changes: glaciers with decadal response times have nearly kept pace with anthropogenic warming, but those with multi-decadal response times are far from equilibrium, and their additional committed retreat stands well beyond natural variability. These differences have implications for changes in glacier runoff. A simple calculation highlights that transient peaks in area-integrated melt, either at the onset of forcing or due to variations in forcing, depend on the glacier's response time and degree of disequilibrium. We conclude that differences in individual response times should be considered when assessing the state of a population of glaciers and modeling their future response. These differences in response can arise simply from a range of different glacier geometries, and the same basic principles can be expected in other regions as well.