Deep gravel-pit exposures reveal the distribution and structure of till and underlying sand and gravel in drumlins near Waukesha, Wisconsin. The subglacial sediment is interpreted to have moved laterally into the drumlin sites because the till thickens from the margin to the core of the drumlins, the stone orientation in the till is perpendicular and oblique to ice flow on the drumlin margins, and recumbent isoclinal folds occur in sand on the drumlin margins with axes parallel to the drumlin axes. The resulting accumulations of sediment presented obstacles to ice flow and were streamlined into the minimum-drag drumlin shape by erosion on the margins and by remolding of material in the core of the drumlins. These drumlin nuclei may have formed at spots where there was low effective stress on the bed. The subglacial sediment became mobile as a result of high pore pressure that may have developed as ground water and subglacial melt water were trapped behind a frozen bed at the ice margin. Under certain conditions, however, lateral sediment flow might also have occurred when the sediment was frozen.

Glaciological results of the continuing investigations of the Quelccaya ice cap located at lat. 13° 56’ S., long. 70° 50’ W., in the Cordillera Oriental of southern Peru are presented. Ice cores to a depth of 15 m have been retrieved from the summit dome (5650 m), middle dome (5543 m), and south dome (5480 m) and sampled in detail for microparticle, oxygen-isotope, and total-β-activity measurements. Results of these core analyses indicate that although the summit of this ice cap is only 300 m above the annual snow line and the firn is temperate, an interpretable stratigraphic record is preserved. The marked seasonal ice stratigraphy is produced by the marked seasonal variation in regional precipitation. High concentrations of microparticles and β- radioactive material occur during the dry season (May-August). Microparticles deposited during the rainy season are larger than those deposited during the dry season. On the Quelccaya ice cap the most negative δ18O values occur during the warmer rainy season (the opposite occurs in polar regions). The near-surface mean δ value of – 21‰ is remarkably low for this tropical site where the measured mean annual air temperature is – 3°C The seasonality of the microparticles, total β activity, and isotope ratios offers the prospect of a climatic ice-core record from this tropical ice cap.

In 1983 three directional surveys were made in the bore hole from which a deep ice core was obtained in the summers of 1979–81. The inclination and azimuth of the bore hole were measured on three surveys, temperature was included on two surveys, fluid pressure and hole diameter on one of the surveys. Fluid-pressure measurements show that the ice-overburden pressure was undercompensated in the upper few hundred meters and overcompensated at the bottom of the hole. Diameter measurements show closure in the upper portion and expansion near the bottom beginning at the transition from the Holocene to Wisconsin ice at 1784 m. The hole expansion and increase in inclination correlate with dust and silt content in the Wisconsin ice. Changes in azimuth are due to flow of the ice and are consistent with the direction of flow at the surface. Temperature measurements show that the hole is at or near equilibrium. The gradient of 0.012 K/m below 1400 m is less than the 0.018 K/m at Camp Century. There is a slight reduction in gradient near the bottom from internal friction in the silty ice.

The Tian Shan mountain range, known as the water towers of Central Asia, plays a key role in local water supply, yet large uncertainties remain about the amount of water that is stored in its glaciers. In this study, we assess the impact of the boundary conditions on ice thickness estimates using two inversion models: a mass conservation (MC) model and a basal shear stress (BS) model. We compare the widely used Randolph Glacier Inventory version 6 with the updated Glacier Area Mapping for Discharge from the Asian Mountains glacier inventory, as well as two digital elevation models (SRTM DEM and Copernicus DEM). The results show that the ice volume (in ~2000 CE) in the Tian Shan range is 661.0 ± 163.5 km3 for the MC model and 552.8 ± 85.3 km3 for the BS model. There are strong regional differences due to inventory, especially for glaciers in China (17–25%). However, the effect of different DEM sources on ice volume estimation is limited. By the end of the 21st century, the projected mass loss differences between inventories are higher than between adjacent emission scenarios, illustrating the vital importance of high-quality inventories. These differences should be carefully considered during water resource planning.

A general survey of the preliminary results of a three-year program of measurements on the Amery Ice Shelf by A.N.A.R.E. are presented, together with theoretical considerations of the velocity and stress distributions and the mass and energy regimes of the ice shelf.

In order to explain the observed velocity distribution it has been found necessary to extend Weertman’s theory of ice-shelf creep to an ice shelf bounded at its sides. The resulting theoretical velocity profile applied to the results of the Amery Ice Shelf provides estimates of the average values of the power flow-law parameters for the ice shelf.

The energy and mass budget considerations, together with the recorded change in form of the ice front, suggest that the ice-shelf regime is not in a continual state of balance but may fluctuate as the ice shelf changes in form over a period of about forty years.

Crystalline texture and c-axis orientation of the 2083 m ice core at Vostok Station, covering more than 150kyear, reveal the existence of strong anisotropics. Changes in crystal size with depth are compatible with the growth of grains driven by the free energy of grain boundaries. A smaller growth rate appears to be associated with cold periods. A gradual increase in the horizontal elongation of grains was observed between 350 and 680 m. But, the mean value of the coefficient of the linear dimensional orientation of grains does not change below 700 m.

The c-axis orientation of ice grains tends to orientate perpendicular to the direction of the elongation of grains, forming a vertical girdle pattern. This characteristic fabric has been interpreted as resulting from the gradual rotation of grains by basal glide under uniaxial longitudinal tension. The rotation of grains was calculated with respect to the total strain, simulating the formation of the girdle fabric pattern. The fabric-enhancement factor was calculated at various depths. It appears that Vostok ice hardens gradually with depth when considering the transverse convergent flow. No significant variation of the enhancement factor was observed with changes in climate and impurity content.

Avalanche accidents are increasing in Scotland and fatalities now average two each year. The victims are either mountaineers or cross-country skiers as the problem is confined to the steeper and remoter slopes. Studies in the Cairngorm Mountains have shown that Scottish snow differs from alpine snow in three respects: depth hoar does not occur in Scotland, higher values are found for density and rammsonde penetration, and ice layers are more abundant. Slab avalanches are found on slopes of about 35° and sluffs on slopes of 30° and steeper. Although most avalanches are small, crown fracture lines reach a maximum length of 500 m, and slabs may be up to eight metres thick. Preliminary correlations between avalanche release and weather conditions have shown the importance of cold periods following heavy snow falls and the significance of thaws. Much more needs to be done before an avalanche-forecasting service could be provided, not the least of which is persuading likely beneficiaries that such a service would be worthwhile.

A degree-day model developed for parameterizing melt rates on the Greenland ice sheet is adapted to the climatic conditions on glaciers in southern Norway. The model is calibrated by means of observed average mass-balance-elevation relationships (1963–90) for three glaciers in a west-east transect in southern Norway and 30 year normals (1961–90) of temperature and precipitation observed at nearby climate stations. The calibration gives a surprisingly small variation of the model parameters (degree-day factors for snow-and ice-melt, and precipitation-elevation gradient) from one glacier to another. The derived values of the parameters are used to estimate the change of the mass-balance-elevation relationship for different climatic scenarios. The study indicates that a low-lying glacier in the maritime, high-precipitation environment near the Atlantic coast is more sensitive to both temperature and precipitation changes than the high elevated glaciers in the dry, more continental climate farther away from the coast. However, all of the glaciers studied will lose mass in a warmer climate, unless the warming is accompanied by a dramatic increase in the precipitation of 25–40% deg−1 warming.

The spatial distribution of basal water critically impacts the evolution of ice sheets. Current estimates of basal water distribution beneath the Greenland Ice Sheet (GrIS) contain large uncertainties due to poorly constrained boundary conditions, primarily from geothermal heat flux (GHF). The existing GHF models often contradict each other and implementing them in numerical ice-sheet models cannot reproduce the measured temperatures at ice core locations. Here we utilize two datasets of radar-detected basal water in Greenland to constrain the GHF at regions with a thawed bed. Using the three-dimensional ice-sheet model SICOPOLIS, we iteratively adjust the GHF to find the minimum GHF required to reach the bed to the pressure melting point, GHFpmp, at locations of radar-detected basal water. We identify parts of the central-east, south and northwest Greenland with significantly high GHFpmp. Conversely, we find that the majority of low-elevation regions of west Greenland and parts of northeast have very low GHFpmp. We compare the estimated constraints with the available GHF models for Greenland and show that GHF models often do not honor the estimated constraints. Our results highlight the need for community effort to reconcile the discrepancies between radar data, GHF models, and ice core information.

We recorded the burial times of temperature sensors mounted on a specially constructed tower to determine snow accumulation during individual storms in the summit caldera of Mount Wrangell, Alaska, USA, (62° N, 144° W; 4100 m a.s.l.) during the accumulation year June 2005 to June 2006. The experiment showed most of the accumulation occurred in episodic large storms, and half of the total accumulation was delivered in late summer. The timing of individual events correlated well with storms recorded upwind, at Cordova, the closest Pacific coastal weather station (200 km south-southeast), although the magnitude of events showed only poor correlation. Hence, snow accumulation at Mount Wrangell appears to be a reflection of synoptic-scale regional weather systems. The accumulation at Mount Wrangell’s summit (>2.5 m w.e.) exceeded the precipitation at Cordova. Although the direct relationship between accumulation of individual storms at the summit of Mount Wrangell and precipitation events at Cordova may be unique in the region, it is useful for interpreting ice cores obtained on Mount Wrangell. This is especially the case here because the high rate of accumulation allows high time resolution within the core.

Processes governing meltwater penetration into cold firn remain poorly constrained. Here, in situ experiments are used to develop a grain-scale model to investigate physical limitations on meltwater infiltration in firn. At two sites in Greenland, drilling pumped water into cold firn to >75 m depth, and the thermo-hydrologic evolution of the firn column was measured. Rather than filling all available pore space, the water formed perched aquifers with downward penetration halted by thermal and density conditions. The two sites formed deep aquifers at ~40 m depth and at densities considerably less than the air pore close-off density (~725 kg m−3 at −18°C, and ~750 kg m−3 at −14°C), demonstrating that some pore space at depth remains inaccessible. A geometric grain-scale model of firn is constructed to quantify the limits of a descending fully saturated wetting front in cold firn. Agreement between the model and field data implies the model includes the first-order effects of water and heat flow in a firn lattice. The model constrains the relative importance of firn density, temperature and grain/pore size in inhibiting wetting front migration. Results imply that deep infiltration, including that which leads to firn aquifer formation, does not have access to all available firn pore space.

Outlet glaciers undergo rapid spatial and temporal changes in flow velocity during calving events. Observing such changes requires both high temporal and high spatial resolution methods, something now possible with terrestrial radar interferometry. While a single such radar provides line-of-sight velocity, two radars define both components of the horizontal flow field. To assess the feasibility of obtaining the two-dimensional (2-D) flow field, we deployed two terrestrial radar interferometers at Jakobshavn Isbrae, a major outlet glacier on Greenland's west coast, in the summer of 2012. Here, we develop and demonstrate a method to combine the line-of-sight velocity data from two synchronized radars to produce a 2-D velocity field from a single (3 min) interferogram. Results are compared with the more traditional feature-tracking data obtained from the same radar, averaged over a longer period. We demonstrate the potential and limitations of this new dual-radar approach for obtaining high spatial and temporal resolution 2-D velocity fields at outlet glaciers.

The growth of bonds between snow grains or ice spheres has been variously attributed to surface diffusion, volume diffusion, and evaporation–condensation. To distinguish among these possible mechanisms the unconfined compressive strength of two groups of snow samples was determined as a function of time. One group was allowed to sinter under atmospheric conditions while the other group was kept immersed in silicone oil. The much lower rate of strengthening of the latter group suggests that evaporation–condensation must be the major mechanism of mass transport in snow under atmospheric conditions. The possible magnitudes of the various mass transfer coefficients are discussed.

The inaugural data from the first systematic program of sea-ice observations in Kotzebue Sound, Alaska, in 2018 coincided with the first winter in living memory when the Sound was not choked with ice. The following winter of 2018–19 was even warmer and characterized by even less ice. Here we discuss the mass balance of landfast ice near Kotzebue (Qikiqtaġruk) during these two anomalously warm winters. We use in situ observations and a 1-D thermodynamic model to address three research questions developed in partnership with an Indigenous Advisory Council. In doing so, we improve our understanding of connections between landfast ice mass balance, marine mammals and subsistence hunting. Specifically, we show: (i) ice growth stopped unusually early due to strong vertical ocean heat flux, which also likely contributed to early start to bearded seal hunting; (ii) unusually thin ice contributed to widespread surface flooding. The associated snow ice formation partly offset the reduced ice growth, but the flooding likely had a negative impact on ringed seal habitat; (iii) sea ice near Kotzebue during the winters of 2017–18 and 2018–19 was likely the thinnest since at least 1945, driven by a combination of warm air temperatures and a persistent ocean heat flux.

The accuracy of sea-ice motion products provided by the National Snow and Ice Data Center (NSIDC) and the Ocean and Sea Ice Satellite Application Facility (OSI-SAF) was validated with data collected by ice drifters that were deployed in the western Arctic Ocean in 2014 and 2016. Data from both NSIDC and OSI-SAF products exhibited statistically significant (p < 0.001) correlation with drifter data. The OSI-SAF product tended to overestimate ice speed, while underestimation was demonstrated for the NSIDC product, especially for the melt season and the marginal ice zone. Monthly Lagrangian trajectories of ice floes were reconstructed using the products. Larger spatial variability in the deviation between NSIDC and drifter trajectories was observed than that of OSI-SAF, and seasonal variability in the deviation for NSIDC was observed. Furthermore, trajectories reconstructed using the NSIDC product were sensitive to variations in sea-ice concentration. The feasibility of using remote-sensing products to characterize sea-ice deformation was assessed by evaluating the distance between two arbitrary positions as estimated by the products. Compared with the OSI-SAF product, relative errors are lower (<11.6%), and spatial-temporal resolutions are higher in the NSIDC product, which makes it more suitable for estimating sea-ice deformation.

Cylindrical samples containing 0.59 mm to 0.84 mm diameter silica sand at about 97% and 55% ice saturation (the ratio of ice volume to sand pore volume) were tested at a temperature of −12° C in triaxial compression. Both constant axial strain-rate tests and step-stress creep tests provide information on the influence of confining pressure on the shear strength and creep behavior of the sand–ice material. Changes in the degree of ice saturation help show the influence of the ice matrix versus the sand material on the mechanical behavior. Data are discussed in terms of the Mohr–Coulomb failure law and creep theories. It is shown that the cohesive component of strength depends on response of the ice matrix, whereas the frictional component of strength responds in a manner very similar to unfrozen sand tested at high confining pressures. Experimental data show that creep rates decrease exponentially and creep strength increases with an increase in confining pressure.

Speke Glacier, the best studied and largest single glacier in the Ruwenzori Range (Uganda, East Africa), was revisited in January 1990. A comparison of its present state with former observations reveals that its terminus receded 35–45 m between 1958 and 1977; since 1977 it has retreated more than 150 m. Nearly half the width of the snout has been lost, mainly on its western side, and the entire western section of the south-facing glacier has shrunk substantially. The observed fluctuation conforms with changes occurring on other tropical glaciers. Explanations concerning the reasons for the glacier recession are presented.

A hybrid analytical and random-walk model has been developed to predict the shape of ice accreted on a horizontal cylindrical insulator due to freezing rain. The freezing rain occurs with the temperature of the vertically falling raindrops above the freezing point and the air temperature below freezing. The analytical model calculates the angular distribution of the water-film temperature and the location where freezing begins. The random-walk model predicts the accretion shape. The two random-walk model parameters, the freezing probability and the shedding parameter, are expressed as functions of the atmospheric conditions. The model predicts a variety of realistic accretion shapes from cylindrical to icicle-like. Model verification based on comparisons with other models and with experimental results demonstrates quantitatively and qualitatively the credibility of this new modelling approach.

Using changes observed in daily seismic reflections, we have investigated the basal morphology of Black Rapids Glacier, Alaska, U.S.A. The englacial drainage of ice-marginal lakes caused significant changes in the daily reflections, as well as dramatic increases in basal motion. Changes in reflection arrival times and amplitudes indicate that there is a basal till layer at least 5 m thick at some locations beneath this surge-type glacier. Rapid changes in the observed reflection coefficients during the drainage events indicate that changes in till properties must occur throughout the entire 5 m thick layer, they must last for several days following the lake drainages and they must be completely reversible over as little as 36 min. Our seismic analysis shows that changes in effective pressure of the till are unlikely to cause the required changes in the reflection coefficients, but that a decrease in till saturation is likely. We therefore interpret the cause of the seismic anomalies as being a temporary decrease in saturation as water is input to the subglacial hydraulic system, and propose that such a change may occur quickly and reversibly by a redistribution of overburden pressure. Higher water pressures within the hydraulic system cause that region to support more of the glacier’s weight, leaving the remaining areas to support less. Any till within these areas of decreased normal stress would experience a consequent decrease in pore-water pressure, causing gas to exolve, thus decreasing saturation. This decrease in saturation would cause a change in the strength of the basal layer and may affect basal dynamics.