To quantify recent thinning of marine-terminating outlet glaciers in northwestern Greenland, we carried out field and satellite observations near the terminus of Bowdoin Glacier. These data were used to compute the change in surface elevation from 2007 to 2013 and this rate of thinning was then compared with that of the adjacent land-terminating Tugto Glacier. Comparing DEMs of 2007 and 2010 shows that Bowdoin Glacier is thinning more rapidly (4.1 ± 0.3 m a−1) than Tugto Glacier (2.8 ± 0.3 m a−1). The observed negative surface mass-balance accounts for <40% of the elevation change of Bowdoin Glacier, meaning that the thinning of Bowdoin Glacier cannot be attributable to surface melting alone. The ice speed of Bowdoin Glacier increases down-glacier, reaching 457 m a−1 near the calving front. This flow regime causes longitudinal stretching and vertical compression at a rate of −0.04 a−1. It is likely that this dynamically-controlled thinning has been enhanced by the acceleration of the glacier since 2000. Our measurements indicate that ice dynamics indeed play a predominant role in the rapid thinning of Bowdoin Glacier.

Donjek Glacier has an unusually short and regular surge cycle, with eight surges identified since 1935 from aerial photographs and satellite imagery with a ~12 year repeat interval and ~2 year active phase. Recent surges occurred during a period of long-term negative mass balance and cumulative terminus retreat of 2.5 km since 1874. In contrast to previous work, we find that the constriction where the valley narrows and bedrock lithology changes, 21 km from the terminus, represents the upper limit of surging, with negligible surface speed or elevation change up-glacier from this location. This positions the entire surge-type portion of the glacier in the ablation zone. The constriction geometry does not act as the dynamic balance line, which we consistently find at 8 km from the glacier terminus. During the 2012–2014 surge event, the average lowering rate in the lowest 21 km of the glacier was 9.6 m a−1, while during quiescence it was 1.0 m a−1. Due to reservoir zone refilling, the ablation zone has a positive geodetic balance in years immediately following a surge event. An active surge phase can result in a strongly negative geodetic mass balance over the surge-type portion of the glacier.

We have determined accumulation histories by identifying annual-layer horizons in records obtained by three independent methods: (1) glaciochemical analysis on a core, (2) density profiling in the borehole from which the core was taken, using the neutron-probe (NP) technique, and (3) borehole optical stratigraphy (BOS), again in the same borehole. We also used three different techniques for determining density to convert annual-layer thickness to accumulation: (1) gravimetric measurements on core samples, (2) measurement of density using NP and (3) a simple empirical model based on regional climatology. The result is nine different accumulation time series, three of which are completely independent. The chemical-analysis- and NP-derived accumulation time series are correlated, and the ∼70 year means are in agreement. The BOS-derived accumulation ∼70 year mean is slightly lower, probably due to a combination of the empirical density model’s underestimate of the density profile and the misidentification of sub-annual events in the shallow part of the borehole as annual horizons.

The longitudinal pattern of surface velocity of a large, predominantly cold, polythermal glacier (John Evans Glacier, Ellesmere Island, Canada) was measured over summer and winter periods. In the accumulation and upper ablation areas, where ice is predominantly cold-based, summer velocities were slightly higher than overwinter velocities. Predicted velocities due to ice deformation alone in these areas closely matched these observations in the winter, with limited basal motion likely in the summer. In the lower ablation area, where ice is likely warm-based, measured summer velocities were up to double overwinter velocities. Predicted ice deformation could not account for all of these measured velocities in either summer or winter. This suggests that basal motion occurs throughout the year over at least part of the lower ablation area. This finding is supported by radio-echo sounding, subglacial drainage reconstructions and analyses of early-summer meltwater chemistry, which suggest that subglacial water is present throughout the year in this region. In summer, basal motion may account for up to 75% of the total surface velocity throughout the lower ablation area. The inferred rate of basal motion increases sharply directly below a set of moulins by which most surface meltwater reaches the glacier bed.

Detailed ground-penetrating radar (GPR) surveys at 50 and 200 MHz on Hansbreen, a polythermal glacier in southern Svalbard, are presented and interpreted. Comparison of the variations in character of the radar reflections with borehole thermometry and water levels in moulins suggests that GPR can be used to study the hydrothermal properties of the glacier. The high resolution of the GPR data shows that the hydrothermal structure of the glacier is highly variable both along the centre line and on transverse profiles. Water contents for many places and depths within the glacier were calculated by estimating radar-wave velocities to point reflectors. We find typical water contents of 1-2% for the temperate ice, but wetter ice associated with surface crevassing and moulins (typically 4% water content). There is evidence that wet ice sometimes overlays drier ice. The hydrothermal structure is thus shown to be very complex. Temperature gradients in the cold ice indicate freezing rates of temperate ice below cold ice of 0.1-0.5 ma-1, while isolated point reflectors within the cold ice indicate large water-filled bodies that are probably related to the regular drainage structure of the glacier.

Seismic monitoring on three Cascade volcanoes in Washington State, U.S.A., has shown that small seismic events are generated by the active glaciers found on each mountain. Detailed seismic experiments have been conducted to investigate the sources for these icequakes. Considerable evidence indicates that the events are the result of a stick–slip type of motion taking place at the bed of the glaciers. The few events we have been able to locate had depths comparable with the glaciers’ thickness. The similarity of wave form from an explosion at the bottom of a glacier and natural icequakes suggests that the complexity of the seismic wave form is due to the path and not the source. The events exhibit an annual trend with more events being recorded in the summer than during the winter. A ten-fold increase in the number of events preceded a large ice avalanche that involved the entire glacier thickness, suggesting that seismic monitoring may be useful in predicting catastrophic ice movements.

Herein, we present results obtained from time-lapse imagery acquired by a digital single-lens reflex camera during 2014–18 used to monitor the Planpincieux Glacier on the Italian side of the Grandes Jorasses (Mont Blanc massif). We processed the images using image cross-correlation to measure the surface kinematics of the most active lobe of the glacier that presents a high velocity and frequent ice detachments. During the monitoring, we observed two or three periods of sharp acceleration per year that culminated with large break-offs followed by analogous decelerations. Overall, we registered more than 350 failures with a volume >100 m3, of which, 14 events had volumes larger than 10 000 m3. The study identified a monotonic positive relationship between the velocity and failure volume that may be used to estimate the volume of the collapses before an event. We identified the thresholds of velocity and acceleration that characterise the activation of the speed-up periods. The study allowed the characterisation of three different instability processes that lead to the break-off of ice chunks from the glacier terminus: (i) disaggregation, (ii) slab fracture and (iii) water tunnelling failure which can be differentiated based on the rheology, the volume involved and the trigger process.

Recent surface elevation changes of Kaskawulsh Glacier, Yukon, Canada, are quantified by comparing an air-photo derived DEM from 1977 and airborne lidar measurements from 1995, 2000 and 2007. Surface-area changes are assessed using historical aerial photography from 1956 and satellite imagery from 1977 to 2007. Combined, these measurements provide some of the first detailed records of volume change of a large Yukon glacier. Between 1977 and 2007, Kaskawulsh Glacier underwent a decrease in area of 1.53% and a decrease in volume of 3.27–5.94 km3 w.e.). The terminus also retreated by 655 m over the period 1956–2007. There was relatively minor volume change over the period 1977–95 (<+0.01 km3 w.e.a−1), while over the periods 1995–2000 and 2000–07 volume losses occurred at a relatively constant rate of −0.51 and −0.50 km3 a−1 w.e., respectively. Since 1995, thinning has been prominent throughout the ablation zone, while relative stability and even slight thickening has occurred in the accumulation zone. These findings are similar to those recently observed at other nearby Alaskan glaciers.

The Tibetan Plateau (TP) lakes are sensitive to climate change due to ice-albedo feedback, but almost no study has paid attention to the ice albedo of TP lakes and its potential impacts. Here we present a recent field experiment for observing the lake ice albedo in the TP, and evaluate the applicability of the Moderate Resolution Imaging Spectroradiometer (MODIS) products as well as ice-albedo parameterizations. Most of the observed lake ice albedos on TP are <0.12, and the clear blue ice albedo is only 0.075, much lower than reported in the previous studies. Even that of ice covered with snow patches is only 0.212. MOD10A1 albedo product has the best agreement with observations, followed by those of MYD10A1. MCD43A3 product is consistently higher than the observations. Due to an error of snow flag and inconsistent time windows in MCD43A2 and MCD43A3, at certain times, the albedo of the ice without snow is even higher than that covered with snow. When the solar zenith angle is not considered, there is no significant correlation between the albedo and the ice surface temperature. None of the existing ice-albedo parameterizations can reproduce well the observed relationship of the albedo and surface temperature.

It is generally assumed that the mass-balance gradient on glaciers is more or less conserved under climatic change. In studies of the dynamic response of glaciers to climatic change, one of the following assumptions is normally made: (i) the mass-balance perturbation is independent of altitude or (ii) the mass-balance profile does not change — it simply shifts up and down. Observational evidence for such an approach is not convincing; on some glaciers the inter-annual changes in mass balance seem to be independent of altitude, on others not at all. Moreover, it is questionable whether inter-annual variation can be “projected“ on different climatic states.

To see what a physical approach might contribute, we developed an altitude-dependent mass-balance model. It is based on the energy balance of the ice/snow surface, where precipitation is included in a parameterized form and numerical integrations are done through an entire balance year (with a 30 min time step). Atmospheric temperature, snowfall, and atmospheric transmissivity for solar radiation are all dependent on altitude, so a mass-balance profile can be calculated. Slope and exposure of the ice/snow surface are taken into account (and the effects of these parameters studied). In general, the calculations were done for 100m elevation intervals.

Climatological data from the Sonnblick Observatory (Austria; 3106 m a.s.l.) and from Vent (2000 m a.s.l.; Oetztal Alps, Austria) served as input for a number of runs. Simulation of the mass-balance profiles for Hinterseisferner (north-easterly exposure) and Kesselwandferner (south-easterly exposure) yields reasonable results. The larger balance gradient on Kesselwandferner is produced by the model, so exposure appears to be an important factor here.

Sensitivity of mass-balance profiles to shading effects, different slope, and exposure are systematically studied. Another section deals with the sensitivity to climatic change. Perturbations of air temperature, cloudiness, albedo, and precipitation are imposed to see their effects on the mass-balance profiles. The results clearly show that, in general, mass-balance perturbations depend strongly on altitude. They generally increase down-glacier, and are not always symmetric about the reference state.

For typical climatic conditions in the Alps, we found that a 1 K temperature change leads to a change in equilibrium-line altitude of 130 m. Three factors contribute to this large value; turbulent heat flux, longwave radiation from the atmosphere, and fraction of precipitation falling as snow. Here, the albedo feed-back increases the sensitivity in a significant way.

A new classification of snow on the ground is based on the major physical processes involved in the metamorphism of a snow cover. The major divisions are based on (I) the mechanical damage to snow crystals during precipitation, (II) the transport of water vapor at constant temperature because of surface-energy differences, (III) the transport of water vapor along a thermal gradient, and (IV) firnification because of melting and refreezing, and pressure consolidation.

A review of available information on six self-dumping glacier-dammed lakes indicates that instantaneous discharge during a jökulhlaup is proportional to some power of the cumulative volume of water released rather than being a direct function of time. Information on these and four other self-dumping lakes, moreover, indicates that the peak discharge at the end of each jökulhlaup is approximately proportional to the two-thirds power of the total volume of water released during the flood.

Surface-to-bedrock cores obtained with a CRREL thermal drill were taken in 1972 and 1973 from the top of the Devon Island ice cap. There are very pronounced variations in oxygen isotope, micro-particle concentration, and ice texture in the lowermost 5 m of the core. There is a section of isotopically cold, very fine bubbly ice with high micro-particle concentrations between 2.6 and 4.4 m above the bed, considered to represent the Last Ice Age. There is coarse, isotopically warm, clean ice above and below this. For 1.2 m above the bed, the ice is finer again with high micro-particle concentrations but it shows very low bubble concentration and is isotopically the warmest in the core. While the broad variations are common to both cores, in detail there are significant variations despite the fact that the cores were taken only 27 m apart. The variations, when analysed statistically, show that at least 25–30% of the originally continuous profile is missing from each core. Faulting within the near-bedrock ice may be responsible for some of the effect but bubble fabric also gives evidence for irregular non-laminar flow. Because of the strong relationship between crystal size and micro-particle concentrations in the Devon Island cores, it is suggested that the fine-grained nature of dirty layers in the Antarctic and Greenland ice sheets is due to the effect of the dirt inclusions and not of shearing. Steep isotopic gradients in the Devon Island cores are shown to be evidence for possible shearing, which does not effect any change in the crystal texture. Clear ice near the bed is considered a tectonic feature, but the lack of effect on its bed by the ice cap confirms the non-erosional nature of an ice cap frozen to its bed.

In terms of paleoclimatic history, it means that, because of bedrock effects, ice caps of intermediate depth (i.e. <400 m) can give continuous information only over the last approximate 5 000 years. Between 5 000 and 10 000 B.P. the time series becomes slightly discontinuous and beyond 10 000 B.P. so discontinuous as to allow only broad climatic inferences to be drawn.

A 141m ice core was recovered from Combatant Col (51.385° N, 125.258° W; 3000ma.s.l.), Mount Waddington, Coast Mountains, British Columbia, Canada. Records of black carbon, dust, lead and water stable isotopes demonstrate that unambiguous seasonality is preserved throughout the core, despite summer surface snowmelt and temperate ice. High accumulation rates at the site (>4 m ice eq. a-1) limit modification of annual stratigraphy by percolation of surface meltwater. The ice-core record spans the period 1973–2010. An annually averaged time series of lead concentrations from the core correlates well with historical records of lead emission from North America, and with ice-core records of lead from the Greenland ice sheet. The depth-age scale for the ice core provides sufficient constraint on the vertical strain to allow estimation of the age of the ice at bedrock. Total ice thickness at Combatant Col is ~250 m; an ice core to bedrock would likely contain ice in excess of 200 years in age. Accumulation at Combatant Col is significantly correlated with both regional precipitation and large-scale geopotential height anomalies.

Salinity profiles of young sea ice are presented which show salinity differences of the order of 2‰ to 3‰ within a ~3 cm. distance, horizontally as well as vertically. This ice is also shown, by use of a penetrating dye, to be very permeable which encumbers any attempt to explain such pronounced salinity variations. Two generations of brine-drainage features are described; the first is inherited from the liquid and the later-generation channels develop from brine-pocket and Tyndall-figure migration.

It is proposed that at least part of the brine contained in the ice mass migrates downward in the following way: cold dense brine cascades from higher levels, causing the more dilute, warmer brine at lower levels to partially freeze. This would also account for some of the local variations in salinity.

Extensive observations on Nivlisen, an ice shelf on Antarctica’s Atlantic coast, are analyzed and combined to obtain a new description of its complex glaciological regime. We generate models of ice thickness (primarily from ground-penetrating radar), ellipsoidal ice surface height (primarily from ERS-1 satellite altimetry), freeboard height (by utilizing precise sea surface information) and ice-flow velocity (from ERS-1/-2 SAR interferometry and GPS measurements). Accuracy assessments are included. Exploiting the hydrostatic equilibrium relation, we infer the ‘apparent air layer thickness’ as a useful measure for a glacier’s density deviation from a pure ice body. This parameter exhibits a distinct spatial variation (ranging from ≈2 to ≈16m) which we attribute to the transition from an ablation area to an accumulation area. We compute mass-flux and mass-balance parameters on a local and areally integrated scale. The combined effect of bottom mass balance and temporal change averaged over an essential part of Nivlisen is –654 ± 170 kg m–2 a–1, which suggests bottom melting processes dominate. We discuss our results in view of temporal ice-mass changes (including remarks on historical observations), basal processes, near-surface processes and ice-flow dynamical features. The question of temporal changes remains open from the data at hand, and we recommend further observations and analyses for its solution.

Using in-situ measured data from Qiyi Glacier, in combination with meteorological and run-off data from stations, a distributed degree-day model was developed for 631 investigated glaciers in the Beida River catchment to explore glacier mass change and its effect on streamflow. The results showed that the average mass balance was −272 ± 67 mm w.e. a−1, with an ice loss of 3.99 Gt during 1957–2013. Assuming a continuous linear trend, equilibrium line altitude rose by 242 m. Compared with morpho-topographic variables, climatic control is a more important factor affecting glacier change. Mass-balance sensitivity to air temperature was −239 mm w.e.°C−1 a−1, while to precipitation it was +1.1 mm w.e. mm−1 a−1. That is, a 210 mm increase in precipitation would be needed to compensate for the net mass loss induced by an air temperature increase of 1°C. Average annual glacier meltwater runoff was 1.51 × 108 m3 from 1957 to 2013, accounting for 15.2% of surface runoff. The time series of meltwater runoff changed abruptly in 2000, and its contribution to surface runoff increased from 13.9 to 20.4%.

We investigate how temperature gradient and initial density influence depth-hoar growth in snow and seek to better define the range of conditions under which cohesive, or hard, depth hoar forms. Samples of 400 kg m−3 sieved snow were exposed to temperature gradients of 20–80°C m−1, and samples of four different densities were exposed to a 40°C m temperature gradient. Following exposure to temperature gradients, penetrometer tests were made on samples to determine the presence of solid and/or hard depth hoar. Grain bond orientation was analyzed in section planes by two-dimensional stereological techniques where hard depth hoar developed. Results indicate that hard cohesive depth hoar forms from rounded-grain snow having a density of 400 kg m−3 or greater, following exposure to a temperature gradient of 20°C m−1 or greater. Hard depth hoar appears to consist of solid-type depth-hoar grains connected by necks, with vertically preferred directions of grain elongation and organization of grain-to-grain chains. This work corroborates Atikaya’s (1974) results, but extends his observation of formation of hard depth hoar to weaker temperature gradients for high-density snow. Our results also indicate that hard depth hoar is composed of faceted solid-type (anhedral) grains.

An experimental investigation into the mechanical behaviour of polycrystalline ice in triaxial compression has been conducted using conditions generally favourable to brittle fracture and microcracking. Under triaxial stresses at high strain rate, ice failure occurs by abrupt shear fracturing, generally inclined at about 45° to the maximum principal stress. At −20°C, such failure is suppressed by the imposition of a small confining pressure, allowing a transition to ductile-type flow accompanied by distributed microcracking, but at —40°C shear fracture persists under confinement of up to at least 50 MPa. For low confining pressures (< 10 MPa), brittle strength is strongly pressure-dependent; above this it is pressure-independent. Evidence is presented that suggests this may reflect a change from a fracture process influenced by friction to fracture initiated by localized yielding. Ductile yield strength is found to be little influenced by confining pressure despite the inhibition of cracking that leads to greatly contrasting observed crack densities. Flow conforms to the well-known power law for ice with Q = 69 J mol−1 and n = 4.2 over the temperature range −20° to −4-5° C Under these conditions, microcracking in ice appears to remain remarkably stable and non-interacting.

The effects of silt-sized particles (average diameter of 50 μm) on the compressive creep of polycrystalline ice have been studied at stress levels from 0.1 to 1.45 MPa and temperatures of –12ºC and –10°C. Dislocation densities during creep have been estimated using a dislocation-based model of anelasticity. The results indicate that at low concentrations (up to 4wt.%), particles increase the minimum creep rate. Power-law behavior with an exponent of 3 was observed for both particle-free ice and ice with 1 wt.% particles when the stress was >0.3 MPa. In contrast, linear behavior was observed when the stress was <0.3 MPa. Calculations show that the linear behavior is associated with a constant dislocation density, and the power-law behavior is associated with increasing dislocation densities with increasing stress.