The response of glaciers to climate change has major implications for sea-level change and water resources around the globe. Large-scale glacier evolution models are used to project glacier runoff and mass loss, but are constrained by limited observations, which result in models being over-parameterized. Recent systematic geodetic mass-balance observations provide an opportunity to improve the calibration of glacier evolution models. In this study, we develop a calibration scheme for a glacier evolution model using a Bayesian inverse model and geodetic mass-balance observations, which enable us to quantify model parameter uncertainty. The Bayesian model is applied to each glacier in High Mountain Asia using Markov chain Monte Carlo methods. After 10,000 steps, the chains generate a sufficient number of independent samples to estimate the properties of the model parameters from the joint posterior distribution. Their spatial distribution shows a clear orographic effect indicating the resolution of climate data is too coarse to resolve temperature and precipitation at high altitudes. Given the glacier evolution model is over-parameterized, particular attention is given to identifiability and the need for future work to integrate additional observations in order to better constrain the plausible sets of model parameters.

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.

On 29 April 2015, four beacons were deployed onto an ice island in the Strait of Belle Isle to record positional data. The ice island later broke up into many fragments, four of which were tracked by the beacons. The relative influences of wind drag, current drag, Coriolis force, sea surface height gradient and sea-ice force on the drift of the tracked ice island fragments were analyzed. Using atmospheric and oceanic model outputs, the sea-ice force was calculated as the residual of the fragments' net forces and the sum of all other forces. This was compared against the force obtained through ice concentration-dependent relationships when sea ice was present. The sea-ice forces calculated from the residual approach and concentration-dependent relationships were significant only when sea ice was present at medium-high concentrations in the vicinity of the ice island fragments. The forces from ocean currents and sea surface tilt contributed the most to the drift of the ice island fragments. Wind, however, played a minimal role in the total force governing the drift of the four ice island fragments, and Coriolis force was significant when the fragments were drifting at higher speeds.

The accumulation region of Fedchenko Glacier represents an extensive snow reservoir in the Pamir Mountains feeding the longest glacier in Central Asia. Observed elevation changes indicate a continuous ice loss in the ablation region of Fedchenko Glacier since 1928, while the mass balance of the accumulation region is largely unknown. In this study, we show that accumulation varies considerably in the main accumulation basin, with accumulation rates up to 2400 mm w.e. a−1 in the West, decreasing to <1000 mm w.e. a−1 in the center, although the elevation difference is <200 m. The combination of snow/firn samples and ground-penetrating radar profiles suggests that this accumulation pattern is persistent during the recent past. The recent accumulation history is reconstructed from internal radar reflectors using a firn densification model and shows strong interannual variations, but near constant mean values since 2002. Modeling of trajectories, based on accumulation and glacier geometry, results in an estimate of the depth/age relation close to the main divide. This region provides one of the most suitable locations for retrieving climate information with temporal high resolution for the last millennium, with a potential to cover most of the Holocene in less detail.

Diurnal depth cycles of decimeter scale are observed in a supraglacial lake on the McMurdo Ice Shelf, Antarctica. We evaluate two possible causes: (1) tidal tilt of the ice shelf in response to the underlying ocean tide, and (2) meltwater input variation. We find the latter to be the most likely explanation of our observations. However, we do not rule out tidal tilt as a source of centimeter scale variations, and point to the possibility that other, larger supraglacial lake systems, particularly those on ice shelves that experience higher amplitude tidal tilts, such as in the Weddell Sea, may have depth cycles driven by ocean tide. The broader significance of diurnal cycles in meltwater depth is that, under circumstances where the ice shelf is thin, tidal-tilt amplitudes are high, and meltwater runoff rates are large, there may be associated flexure stresses that can contribute to ice-shelf fracture and destabilization. For the McMurdo Ice Shelf (~20–50 m thickness, ~ 1 m tidal amplitude and ~10 cm water-depth variations), these stresses amount to several 10's of kPa.

Sources and implications of black carbon (BC) and mineral dust (MD) on two glaciers on the central Tibetan Plateau were estimated based on in situ measurements and modeling. The results indicated that BC and MD accounted for ~11 ± 1% and 4 ± 0% of the albedo reduction relative to clean snow, while the radiative forcing varied between 11 and 196 and 1–89 W m−2, respectively. Assessment of BC and MD contributions to the glacier melt can reach up 88 to 434 and 35 to 187 mm w.e., respectively, contributing ~9–23 and 4–10% of the total glacier melt. A footprint analysis indicated that BC and MD deposited on the glaciers originated mainly from the Middle East, Central Asia, North China and South Asia during the study period. Moreover, a potentially large fraction of BC may have originated from local and regional fossil fuel combustion. This study suggests that BC and MD will enhance glacier melt and provides a scientific basis for regional mitigation efforts.

Glaciers depicted on old maps reveal their historical extents, before the advent of aerial and satellite remote sensing. Digital glacier inventories produced from these maps can be employed in assessments of centennial-scale glacier change. This study reconstructs the ~1899 (covering the period 1882–1916) glacier extent in Nordland, northern Norway, from historical gradteigskart maps, with an emphasis on examining the accuracy of the mapped glaciers. Glacier outlines were digitised from georectified scans of the analogue maps in a raster graphics editor and were subsequently inventoried in a GIS. The accuracy of the historical glacier extent was established from written descriptions and landscape photographs created during the original field surveys, and further validated against independent glacier outlines of (1) the maximum Little Ice Age extent derived from geomorphological evidence, and (2) the 1945 extent derived from vertical aerial photographs. An overall uncertainty of ±17% is associated with our inventory. Nordland's glaciers covered an area of 1712 ± 291 km2 in 1899. By 2000, total ice cover had decreased by 47% (807 ± 137 km2) at a rate of 6% 10 a−1 (80 ± 14 km2 10 a−1). The approach presented here may serve as a blueprint for future studies intending to derive glacier inventories from historical maps.

Recent speleological surveys of meltwater drainage systems in cold and polythermal glaciers have documented dynamic englacial and in some cases subglacial conduits formed by the ‘cut-and-closure’ mechanism. Investigations of the spatial distribution of such conduits often require a combination of different methods. Here, we studied the englacial drainage system in the cold glacier Longyearbreen, Svalbard by combining speleological exploration of a 478 m long meltwater conduit with a high-resolution ground penetrating radar (GPR) survey with two different centre-frequencies (25 and 100 MHz). The results yielded a 3-D documentation of the present englacial drainage system. The study shows that the overall form of englacial conduits can be detected from velocity−depth converted GPR data, and that the 3-D model can facilitate a method to pinpoint the reflections in a radargram corresponding with the englacial drainage system, although fine detail cannot be resolved. Visible reflections approximately parallel to the mapped englacial water drainage system likely result from sediment incorporated in the ice or from abandoned parts of the englacial drainage system.

Surface melt on the coastal Antarctic ice sheet (AIS) determines the viability of its ice shelves and the stability of the grounded ice sheet, but very few in situ melt rate estimates exist to date. Here we present a benchmark dataset of in situ surface melt rates and energy balance from nine sites in the eastern Antarctic Peninsula (AP) and coastal Dronning Maud Land (DML), East Antarctica, seven of which are located on AIS ice shelves. Meteorological time series from eight automatic and one staffed weather station (Neumayer), ranging in length from 15 months to almost 24 years, serve as input for an energy-balance model to obtain consistent surface melt rates and energy-balance results. We find that surface melt rates exhibit large temporal, spatial and process variability. Intermittent summer melt in coastal DML is primarily driven by absorption of shortwave radiation, while non-summer melt events in the eastern AP occur during föhn events that force a large downward directed turbulent flux of sensible heat. We use the in situ surface melt rate dataset to evaluate melt rates from the regional atmospheric climate model RACMO2 and validate a melt product from the QuikSCAT satellite.

The spatial distribution of geothermal heat flux (GHF) under ice sheets is largely unknown. Nonetheless, it is an important boundary condition in ice-sheet models, and suggested to control part of the complex surface velocity patterns observed in some regions. Here we investigate the effect of including subglacial hydrology when modelling ice streams with elevated GHF. We use an idealised ice stream geometry and a thermomechanical ice flow model coupled to subglacial hydrology in the Ice Sheet System Model (ISSM). Our results show that the dynamic response of the ice stream to elevated GHF is greatly enhanced when including the interactive subglacial hydrology. On the other hand, the impact of GHF on ice temperature is reduced when subglacial hydrology is included. In conclusion, the sensitivity of ice stream dynamics to GHF is likely to be underestimated in studies neglecting subglacial hydrology.

In this paper, we give an overview of changes in area, length, surface elevation and mass balance of glaciers in mainland Norway since the 1960s. Frontal advances have been recorded in all regions except the northernmost glaciers in Troms and Finnmark (Storsteinsfjellbreen, Lyngen and Langfjordjøkelen). More than half of the observed glaciers, 27 of 49, had marked advances in the 1990s. The glaciological mass-balance values for the period 1962–2018, where 43 glaciers have been measured, show great inter-annual variability. The results reveal accelerated deficit since 2000, the most negative decade being 2001–2010. Some years with a positive mass balance (or less negative) after 2010s can be attributed to variations in large-scale atmospheric circulation. A surface elevation change and geodetic mass balance were calculated for a sample of 131 glaciers covering 817 km2 in the ‘1960s’ and 734 km2 in the ‘2010s’, giving an area reduction of 84 km2, or 10%. The sample covers many of the largest glaciers in Norway, and they had an overall change in surface elevation of −15.5 m for the ~50 year period. Converted to a geodetic mass balance this gives a mean mass balance of −0.27 ± 0.05 m w.e. a−1.

Rapidly-rising jökulhlaups, or glacial outburst floods, are a phenomenon with a high potential for damage. The initiation and propagation processes of a rapidly-rising jökulhlaup are still not fully understood. Seismic monitoring can contribute to an improved process understanding, but comprehensive long-term seismic monitoring campaigns capturing the dynamics of a rapidly-rising jökulhlaup have not been reported so far. To fill this gap, we installed a seismic network at the marginal, ice-dammed lake of the A.P. Olsen Ice Cap (APO) in NE-Greenland. Episodic outbursts from the lake cause flood waves in the Zackenberg river, characterized by a rapid discharge increase within a few hours. Our 6 months long seismic dataset comprises the whole fill-and-drain cycle of the ice-dammed lake in 2012 and includes one of the most destructive floods recorded so far for the Zackenberg river. Seismic event detection and localization reveals abundant surface crevassing and correlates with changes of the river discharge. Seismic interferometry suggests the existence of a thin basal sedimentary layer. We show that the ballistic part of the first surface waves can potentially be used to infer medium changes in both the ice body and the basal layer. Interpretation of time-lapse interferograms is challenged by a varying ambient noise source distribution.