Submarine melting of tidewater glaciers is proposed as a trigger for their recent thinning, acceleration and retreat. We estimate spring submarine melt rates (SMRs) of Kangiata Nunaata Sermia in southwest Greenland, from 2012 to 2014, by examining changes in along-fjord freeboard and velocity of the seasonal floating ice tongue. Estimated SMRs vary spatially and temporally near the grounding line, with mean rates of 1.3 ± 0.6, 0.8 ± 0.3 and 1.0 ± 0.4 m d−1 across the tongue in 2012, 2013 and 2014, respectively. Higher melt rates correspond with locations of emerging subglacial plumes and terminus calving activity observed during the melt season using time-lapse camera imagery. Modelling of subglacial flow paths suggests a dynamic system capable of rapid re-routing of subglacial discharge both within and between melt seasons. Our results provide an empirically-derived link between the presence of subglacial discharge plumes and areas of high spring submarine melting and calving along glacier termini.

We investigated the effects of microparticles and grain size on the microstructural evolutions and mechanical properties of polycrystalline ice. Uniaxial compression tests were conducted using fine-grained pure ice and silica-dispersed ice under various conditions. Deformation behavior of fine-grained ice was found to be characterized by stress exponent n ≈ 2 and activation energy Q ≈ 60 kJ mol−1. The derived strain rates of fine-grained ice were ≈ 1 order of magnitude larger than those of coarse-grained ice obtained in previous studies, and they were found to be independent of particle dispersion and dependent on the mean grain size of ice, with grain size exponent p ≈ 1.4. Work hardening was observed in dislocation creep, while the strain rate continued to decrease. These results indicate that the deformation mechanism of fine-grained ice is different from typical dislocation creep, often associated with n = 3. Although microparticles restricted grain growth, there was little direct effect on the deformation of fine-grained ice. Microstructural observations of the ice samples indicated that the grain boundaries were straight and that the subgrain boundary densities increased after deformation. Our experiments suggest that grain size and boundaries play important roles in the deformation processes of polycrystalline ice.

Glaciers on active volcanoes are subject to changes in both climate fluctuations and volcanic activity. Whereas many studies analysed changes on individual volcanoes, this study presents for the first time a comparison of glacier changes on active volcanoes on a continental scale. Glacier areas were mapped for 59 volcanoes across Latin America around 1986, 1999 and 2015 using a semi-automated band ratio method combined with manual editing using satellite images from Landsat 4/5/7/8 and Sentinel-2. Area changes were compared with the Smithsonian volcano database to analyse possible glacier–volcano interactions. Over the full period, the mapped area changed from 1399.3 ± 80 km2 to 1016.1 ± 34 km2 (−383.2 km2) or −27.4% (−0.92% a−1) in relative terms. Small glaciers, especially in tropical regions lost more of their area compared to large and extra–tropical glaciers. Interestingly, 46 out of 59 analysed glaciers (78%) showed a decelerating shrinkage rate in the second period (−1.20% a−1 before 1999 and −0.70% a−1 after 1999). We found a slightly higher (but statistically not significant) area loss rate (−1.03% a−1) for glaciers on volcanoes with eruptions than without (−0.86% a−1).

Experimental data from creep tests on polycrystalline ice samples highlight not only the non-Newtonian behavior of ice but also suggest a critical dependence of the various rheological parameters upon the applied hydrostatic pressure. We propose a new modeling framework, based on implicit theories of continuum mechanics, that generalizes two well-known constitutive models by taking into account the effect of the pressure in the description of ice in creep. To ascertain the validity of the proposed models, we fit the physical parameters with experimental data for the elongational flow of ice samples. The results show good agreement with the experimental creep curves. In particular, the proposed generalized models reproduce the increase of the creep rate due to the presence of hydrostatic pressure.

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.

Climate variability can complicate efforts to interpret any long-term glacier mass-balance trends due to anthropogenic warming. Here we examine the impact of climate variability on the seasonal mass-balance records of 14 glaciers throughout Norway, Sweden and Svalbard using dynamical adjustment, a statistical method that removes orthogonal patterns of variability shared between each mass-balance record and sea-level pressure or sea-surface temperature predictor fields. For each glacier, the two leading predictor patterns explain 27–81% of the winter mass-balance variability and 24–69% of the summer mass-balance variability. The spatial and temporal structure of these patterns indicates that accumulation variability for all of the glaciers is strongly related to the North Atlantic Oscillation (NAO), with the Atlantic Multidecadal Oscillation (AMO) also modulating accumulation variability for the northernmost glaciers. Given this result, predicting glacier change in the region may depend on NAO and AMO predictability. In the raw mass-balance records, the glaciers throughout southern Norway have significantly negative summer trends, whereas the glaciers located closer to the Arctic have negative winter trends. Removing the effects of climate variability suggests it can bias trends in mass-balance records that span a few decades, but its effects on most of the longer-term mass-balance trends are minimal.

Surface mass balance (SMB) is the net input of mass on a glacier's upper surface, composed of snow deposition, melt and erosion processes, and is a major contributor to the overall mass balance. Pine Island Glacier (PIG) in West Antarctica has been dynamically imbalanced since the early 1990s, indicating that discharge of solid ice into the oceans exceeds snow deposition. However, observations of the SMB pattern on the fast flowing regions are scarce, and are potentially affected by the firn's strain history. Here, we present new observations from radar-derived stratigraphy and a relatively dense network of firn cores, collected along a ~900 km traverse of PIG. Between 1986 and 2014, the SMB along the traverse was 0.505 m w.e. a−1 on average with a gradient of higher snow deposition in the South-West compared with the North-East of the catchment. We show that along ~80% of the traverse the strain history amounts to a misestimation of SMB below the nominal uncertainty, but can exceed it by a factor 5 in places, making it a significant correction to the SMB estimate locally. We find that the strain correction changes the basin-wide SMB by ~0.7 Gt a−1 and thus forms a negligible (1%) correction to the glacier's total SMB.

We conducted a mass-balance study of debris-free Trambau Glacier in the Rolwaling region, Nepal Himalaya, which is accessible to 6000 m a.s.l., to better understand mass-balance processes and the effect of precipitation on these processes on high-elevation Himalayan glaciers. Continuous in situ meteorological and mass-balance observations that spanned the three melt seasons from May 2016 are reported. An energy- and mass-balance model is also applied to evaluate its performance and sensitivity to various climatic conditions. Glacier-wide mass balances ranging from −0.34 ± 0.38 m w.e. in 2016 to −0.82 ± 0.53 m w.e. in 2017/18 are obtained by combining the observations with model results for the areas above the highest stake. The estimated long-term glacier mass balance, which is reconstructed using the ERA-Interim data calibrated with in situ data, is −0.65 ± 0.39 m w.e. a−1 for the 1980–2018 period. A significant correlation with annual precipitation (r = 0.77, p < 0.001) is observed, whereas there is no discernible correlation with summer mean air temperature. The results indicate the continuous mass loss of Trambau Glacier over the last four decades, which contrasts with the neighbouring Mera Glacier in balance.

Ice cliffs and ponds on debris-covered glaciers have received increased attention due to their role in amplifying local melt. However, very few studies have looked at these features on the catchment scale to determine their patterns and changes in space and time. We have compiled a detailed inventory of cliffs and ponds in the Langtang catchment, central Himalaya, from six high-resolution satellite orthoimages and DEMs between 2006 and 2015, and a historic orthophoto from 1974. Cliffs cover between 1.4% (± 0.4%) in the dry and 3.4% (± 0.9%) in the wet seasons and ponds between 0.6% (± 0.1%) and 1.6% (± 0.3%) of the total debris-covered tongues. We find large variations between seasons, as cliffs and ponds tend to grow in the wetter monsoon period, but there is no obvious trend in total area over the study period. The inventory further shows that cliffs are predominately north-facing irrespective of the glacier flow direction. Both cliffs and ponds appear in higher densities several hundred metres from the terminus in areas where tributaries reach the main glacier tongue. On the largest glacier in the catchment ~10% of all cliffs and ponds persisted over nearly a decade.

Cryoconite is a dark-coloured granular sediment that contains biological and mineralogical components, and it plays a pivotal role in geochemistry, carbon cycling and glacier mass balance. In this work, we collected cryoconite samples from Laohugou Glacier No. 12 (LHG) on the north-eastern Tibetan Plateau during the summer of 2015 and measured the spectral albedo. To explore the impacts of this sediment on surface ablation, the ice melting differences between the cryoconite-free (removed) ice and the intact layers were compared. The results showed that the mean concentrations of black carbon (BC), organic carbon (OC) and total iron (Fe) in the LHG cryoconite were 1.28, 11.18 and 39.94 mg g−1, respectively. BC was found to play a stronger role in solar light adsorption than OC and free Fe. In addition, ice covered by cryoconite exhibited the lowest mean reflectance (i.e., <0.1). Compared with the cryoconite-free ice surface, cryoconite effectively absorbed solar energy and enhanced glacial melting at a rate of 2.27–3.28 cm d−1, and free Fe, BC and OC were estimated to contribute 1.01, 0.99 and 0.76 cm d−1, respectively. This study provides important insights for understanding the role of cryoconite in the glacier mass balance of the northern Tibetan Plateau.

The evolution of Northern Hemisphere ice sheets through the last glacial cycle is simulated with the glacial index method by using the climate forcing from one General Circulation Model, COSMOS. By comparing the simulated results to geological reconstructions, we first show that the modelled climate is capable of capturing the main features of the ice-sheet evolution. However, large deviations exist, likely due to the absence of nonlinear interactions between ice sheet and other climate components. The model uncertainties of the climate forcing are examined using the output from nine climate models from the Paleoclimate Modelling Intercomparison Project Phase III. The results show a large variability in simulated ice sheets between the different models. We find that the ice-sheet extent pattern resembles summer surface air temperature pattern at the Last Glacial Maximum, confirming the dominant role of surface ablation process for high-latitude Northern Hemisphere ice sheets. This study shows the importance of the upper boundary condition for ice-sheet modelling, and implies that careful constraints on climate output is essential for simulating realistic glacial Northern Hemisphere ice sheets.

Increasing ocean and air temperatures have contributed to the removal of floating ice shelves from several Greenland outlet glaciers; however, the specific contribution of these external forcings remains poorly understood. Here we use atmospheric, oceanographic and glaciological time series data from the ice shelf of Petermann Gletscher, NW Greenland to quantify the forcing of the ocean and atmosphere on the ice shelf at a site ~16 km from the grounding line within a large sub-ice-shelf channel. Basal melt rates here indicate a strong seasonality, rising from a winter mean of 2 m a−1 to a maximum of 80 m a−1 during the summer melt season. This increase in basal melt rates confirms the direct link between summer atmospheric warming around Greenland and enhanced ocean-forced melting of its remaining ice shelves. We attribute this enhanced melting to increased discharge of subglacial runoff into the ocean at the grounding line, which strengthens under-ice currents and drives a greater ocean heat flux toward the ice base.

We derive the surface and basal radar reflectance and backscatter coefficients of the southern McMurdo Ice Shelf (SMIS) and part of the nearby Ross Ice Shelf (RIS), Antarctica, from radar statistical reconnaissance using a 60-MHZ airborne survey. The surface coefficients are further inverted in terms of snow density and roughness, providing a spatial distribution of the processes contributing to the surface boundary conditions. We disentangle the basal coefficients from surface transmission losses, and we provide the basal coherent content, an indicator of the boundary geometric disorder that is also self-corrected from englacial attenuation. The basal radar properties exhibit sharp gradients along specific iso-depths, suggesting an abrupt modification of the ice composition and geometric structure. We interpret this behavior as locations where the pressure-melting point is reached, outlining fields of freezing and melting ice. Basal steps are observed at both SMIS and RIS, suggesting a common geometric expression of widespread basal processes. This technique offers a simultaneous view of both the surface and basal boundary conditions to help investigate the ice-shelf stability, while its application to airborne data significantly improves coverage of the difficult-to-observe ice–ocean boundary. It also provides constraints on thermohaline circulation in ice shelves cavities, which are analogs for ice-covered ocean worlds.

To study subglacial hydrological condition and its influence on the glacier dynamics, we drilled Johnsons Glacier on Livingston Island in the Antarctic Peninsula region. Subglacial water pressure was recorded in boreholes at two locations over 2 years, accompanied by high-frequency ice-speed measurements during two summer melt seasons. Water pressure showed two different regimes, namely high frequency and large amplitude variations during the melt season (January–April) and small fluctuations near the overburden pressure the rest of the year. Speed-up events were observed several times in each summer measurement period. Ice motion during these events substantially contributed to total glacier motion, for example, fast ice flow over 1 week accounted for ~70% of the total displacement over a 25-day long measurement period. We did not find a clear relationship between subglacial water pressure and ice speed. This was probably because subglacial hydraulic conditions were spatially inhomogeneous and thus our borehole data did not always represent a large-scale subglacial condition. Ice temperature measurements in the boreholes confirmed the existence of a cold ice layer near the glacier surface. Our data provide a basis to better understand the dynamic and hydrological conditions of relatively unstudied glaciers in the Antarctic Peninsula region.