A glacierized terrain comprises different land covers, and their mapping using satellite data is challenged by their spectral similarity. We propose a hierarchical knowledge-based classification (HKBC) approach for differentiation of glacier terrain classes and mapping of glacier boundaries, using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery and Global Digital Elevation Model (GDEM). The methodology was tested over Kolahoi Glacier, Kashmir Himalaya. For the sequential extraction of various glacier terrain classes, several input layers were generated from the primary datasets by applying image-processing techniques. Noticeable differences in temperature and spectral response between supraglacial debris and periglacial debris facilitated the development of a thermal glacier mask and normalized-difference debris index, which together with slope enabled their differentiation. These and the other layers were then used in several discrete tests in HKBC, to map various glacier terrain classes. An ASTER visible near-infrared image and 42 field points were used to validate results. The proposed approach satisfactorily classified all the glacier terrain classes with an overall accuracy of 89%. The Z-test reveals that results obtained from HKBC are significantly (at 95% confidence level) better than those from a maximum likelihood classifier (MLC). Glacier boundaries obtained from HKBC were found to be plausibly better than those obtained from MLC and visual interpretation.

We present in situ firn temperatures from the extreme 2012 melt season in the southwestern lower accumulation area of the Greenland ice sheet. The upper 2.5 m of snow and firn was temperate during the melt season, when vertical meltwater percolation was inefficient due to a ~5.5 m thick ice layer underlying the temperate firn. Meltwater percolation and refreezing beneath 2.5 m depth only occurred after the melt season. Deviations from temperatures predicted by pure conductivity suggest that meltwater refroze in discrete bands at depths of 2.0–2.5, 5.0–6.0 and 8.0–9.0 m. While we find no indication of meltwater percolation below 9 m depth or complete filling of pore volume above, firn at 10 and 15 m depth was respectively 4.2–4.5°C and 1.7°C higher than in a conductivity-only simulation. Even though meltwater percolation in 2012 was inefficient, firn between 2 and 15 m depth the following winter was on average 4.7°C warmer due to meltwater refreezing. Our observations also suggest that the 2012 firn conditions were preconditioned by two warm summers and ice layer formation in 2010 and 2011. Overall, firn temperatures during the years 2009–13 increased by 0.6°C.

At least in conventional hydrostatic ice-sheet models, the numerical error associated with grounding line dynamics can be reduced by modifications to the discretization scheme. These involve altering the integration formulae for the basal traction and/or driving stress close to the grounding line and exhibit lower – if still first-order – error in the MISMIP3d experiments. MISMIP3d may not represent the variety of real ice streams, in that it lacks strong lateral stresses, and imposes a large basal traction at the grounding line. We study resolution sensitivity in the context of extreme forcing simulations of the entire Antarctic ice sheet, using the BISICLES adaptive mesh ice-sheet model with two schemes: the original treatment, and a scheme, which modifies the discretization of the basal traction. The second scheme does indeed improve accuracy – by around a factor of two – for a given mesh spacing, but $\lesssim 1$ km resolution is still necessary. For example, in coarser resolution simulations Thwaites Glacier retreats so slowly that other ice streams divert its trunk. In contrast, with $\lesssim 1$ km meshes, the same glacier retreats far more quickly and triggers the final phase of West Antarctic collapse a century before any such diversion can take place.

We investigated glacier changes in the Karatal river basin, the largest basin in Zhetysu (Dzhungar) Alatau, Kazakhstan, for the periods 1956–89, 1989–2001 and 2001–12, based on Landsat TM/ETM+ data analysis. In 1989, we found 243 glaciers with a total area of 142.8 km2; by 2012 these had shrunk to 214 glaciers with a total area of 109.3 km2, a decrease of 33.5 km2 over 23 years (1.02%a-1). This very high shrinkage rate is likely connected with a general trend of increasing temperatures, and small glaciers being situated at the relatively low altitude of the outer Zhetysu Alatau ranges. We also analyzed the shrinkage rate of glaciers based on their differences in size, altitude and aspect of slopes, as well as other topographic parameters, in four sub-basins where glacier shrinkage varied between 18% and 39%. Weather-station climate data showed a significant temperature increase and stable precipitation trends over the study period. We conclude that glacierized areas of the Karatal river basin are located in the most unfavorable conditions for glaciation, and as a result showed a higher shrinkage rate than other glacierized areas of the Tien Shan from 1956 to 2012.

Calving mechanisms are still poorly understood and stress states in the vicinity of ice-shelf fronts are insufficiently known for the development of physically motivated calving laws that match observations. A calving model requires the knowledge of maximum tensile stresses. These stresses depend on different simulation approaches and material models. Therefore, this study compares results of a two-dimensional (2-D) continuum approach using finite elements with results of a one-dimensional (1-D) beam model elaborated in Reeh (1968). A purely viscous model, as well as a viscoelastic Maxwell model, is applied for the 2-D case. The maximum tensile stress usually appears at the top surface of an ice shelf. Its location and magnitude are predominantly influenced by the thickness of the ice shelf and the height of the freeboard, the traction-free part at the ice front. More precisely, doubling the thickness leads to twice the stress maximum, while doubling the freeboard, based on changes of the ice density, results in an increase of the stress maximum by 61%. Poisson's ratio controls the evolution of the maximum stress with time. The viscosity and Young's modulus define the characteristic time of the Maxwell model and thus the time to reach the maximum principal stress.

The retention and release of liquid water in glacierized basins was modelled with a conceptual, semi-distributed model of the water and ice balance designed for long-term averages with monthly resolution for 100 m elevation bands. Here we present the components of the liquid water balance of 86 mostly glacierized basins on either side of the main Alpine divide between 10 and 13°E in the period 1998–2006 and compare them with the records of 30 basins monitored from 1970 to 1997. Basin average of liquid water retention has maxima in excess of 100 mm per month in May, often followed by maximum release when the retaining snow matrix melts. Glacier storage peaks in August partly due to ice melt and the ensuing filling of the englacial reservoirs and partly on account of a precipitation maximum. These two components combined to a common maximum of storage in summer in the first period 1970–97 and developed two distinct maxima in the warmer period 1998–2006. A further maximum of liquid water storage that was often found in October is most likely due to a peak in precipitation in the southern part of the study region.

The Institute Ice Stream (IIS) rests on a reverse-sloping bed, extending >150 km upstream into the ~1.8 km deep Robin Subglacial Basin, placing it at the threshold of marine ice-sheet instability. Understanding IIS vulnerability has focused on the effect of grounding-line melting, which is forecast to increase significantly this century. Changes to ice-flow dynamics are also important to IIS stability, yet little is known about them. Here we reveal that the trunk of the IIS occurs downstream of the intersection of three discrete subglacial features; a large ‘active’ subglacial lake, a newly-discovered sharp transition to a zone of weak basal sediments and a major tectonic rift. The border of IIS trunk flow is confined by the sediment on one side, and by a transition between basal melting and freezing at the border with the Bungenstock Ice Rise. By showing how basal sediment and water dictate present-day flow of IIS, we reveal that ice-sheet stability here is dependent on this unusual arrangement.

Glacier naledi are extrusive ice masses that appear in front of glaciers as a consequence of refreezing of meltwater seepage during the accumulation season. These structures provide a unique opportunity to understand subglacial drainage activity during the accumulation season; however, only few detailed studies have previously focused on their characteristics. Here, we investigated glacier-derived naled assemblages in the proglacial zone of the polythermal glacier Werenskioldbreen (27.4 km2) in SW Svalbard. We determined the spatial distribution of naledi using ground penetrating radar surveys. The main subglacial drainage pattern was related to a channel under the medial moraine, and three sources are linked to a distributed subglacial drainage network. The relation between atmospherically-corrected (Ca2+ + Mg2+) and (SO4 2−) in sub-naled waters was closely related to sulphide oxidation coupled with carbonate dissolution (r = 0.99; slope = 1.6). This is consistent with the local lithology, which is dominated by schist containing carbonates. We also found high carbonate saturation indices in pale white ice layers within the naled. We conclude that sulphide oxidation coupled with carbonate dissolution is the dominant chemical weathering process in the subglacial drainage system of Werenskioldbreen during the accumulation season.

The Qinghai-Tibetan Plateau (QTP) is characterized by a cold climate and a large number of lakes. The long ice season necessitates study of the widespread ice covers in the region. An unprecedented multidisciplinary field campaign was conducted on lake ice processes in the central QTP during the period 2019–13. The study lake generally froze up in late October or early November, and broke up in mid or late April, with a maximum ice thickness of 50–70 cm. The mass balances at both ice surface and bottom were measured continuously. Significant ice surface sublimation/ablation was detected and accounted for up to 40% of the whole ice thickness over the ice season. A simple heat-transfer model was developed for the surface ice loss. The calculated values were in good agreement with the observations. They also indicated that atmospheric conditions, including low air humidity and prevailing strong winds, are the primary drivers of the ice surface sublimation.

Supraglacial ponds on debris-covered glaciers present a mechanism of atmosphere/glacier energy transfer that is poorly studied, and only conceptually included in mass-balance studies of debris-covered glaciers. This research advances previous efforts to develop a model of mass and energy balance for supraglacial ponds by applying a free-convection approach to account for energy exchanges at the subaqueous bare-ice surfaces. We develop the model using field data from a pond on Lirung Glacier, Nepal, that was monitored during the 2013 and 2014 monsoon periods. Sensitivity testing is performed for several key parameters, and alternative melt algorithms are compared with the model. The pond acts as a significant recipient of energy for the glacier system, and actively participates in the glacier’s hydrologic system during the monsoon. Melt rates are 2-4 cm d-1 (total of 98.5 m3 over the study period) for bare ice in contact with the pond, and <1 mmd-1 (total of 10.6m3) for the saturated debris zone. The majority of absorbed atmospheric energy leaves the pond system through englacial conduits, delivering sufficient energy to melt 2612 m3 additional ice over the study period (38.4 m3 d-1). Such melting might be expected to lead to subsidence of the glacier surface. Supraglacial ponds efficiently convey atmospheric energy to the glacier’s interior and rapidly promote the downwasting process.

Neighboring tidewater glaciers often exhibit asynchronous dynamic behavior, despite relatively uniform regional atmospheric and oceanic forcings. This variability may be controlled by a combination of local factors, including glacier and fjord geometry, fjord heat content and circulation, and glacier surface melt. In order to characterize and understand contrasts in adjacent tidewater glacier and fjord dynamics, we made coincident ice-ocean-atmosphere observations at high temporal resolution (minutes to weeks) within a 10 000 km2 area near Uummannaq, Greenland. Water column velocity, temperature and salinity measurements reveal systematic differences in neighboring fjords that imply contrasting circulation patterns. The observed ocean velocity and hydrography, combined with numerical modeling, suggest that subglacial discharge plays a major role in setting fjord conditions. In addition, satellite remote sensing of seasonal ice flow speed and terminus position reveal both speedup and slow-down in response to melt, as well as differences in calving style among the neighboring glaciers. Glacier force budgets and modeling also point toward subglacial discharge as a key factor in glacier behavior. For the studied region, individual glacier and fjord geometry modulate subglacial discharge, which leads to contrasts in both fjord and glacier dynamics.

Supraglacial streams are a significant part of the glacial hydrological system and important for understanding the connection between glacial hydrology and glacier dynamics. Here we determine the factors that influence the development of step-pool formation and pulsating flow in a supraglacial stream on Bylot Island, Nunavut. Results show that during the second week of a 2-week study, multiple successive rainfall events occurred, stream temperature increased and ablation decreased; which also caused stream discharge to decrease. In addition, the stream, which flowed over a 13 m high waterfall off the front of Fountain Glacier, rapidly formed 21 step-pools and began to pulsate. The pulsating phenomenon involved the complete stoppage of flow over the waterfall and the subsequent restart between 8 and 20 s later. Pulsating flow resulted from rapid changes in the streambed morphology. In particular, the formation of the step-pool sequence was caused by helical flow around meander bends and hydrologically induced slippage along transverse shear planes, evidenced by observations of high-pressure artesian flow from transverse fractures. Contrary to previous literature, this study shows that high discharge is not necessarily the cause of step-pool formation and pulsating flow within supraglacial streams.

Thermodynamics of a seasonal supraglacial lake were examined based on field data from three summers. At maximum, the lake body consisted of an upper layer with thin ice on top, and a lower layer with slush, hard ice and sediment at the bottom. Sublimation from the upper ice surface averaged to 0.7 mm d−1, and melting in the ice interior averaged to 9.1 mm d−1 during summer. Albedo was on average 0.6 and light attenuation coefficient was ~1 m−1. Averaged over December and January, and over 3 different years, we found that the net solar heating was 137 W m−2, while the losses averaged to 62 W m−2 for the longwave radiation, 16 W m−2 for the sensible heat flux, 24 W m−2 for the latent heat flux and 3 W m−2 for the bottom flux. The depth scale is determined by the light attenuation distance and thermal diffusion coefficient, and the net liquid water volume ranged from 0.5 to 1.0 m in different years. The potential winter growth is more than summer melting, and thus the lake freezes up completely in winter in the present climate.

An assessment of glacier shrinkage (reduction of area) for all of High Mountain Asia requires a complete compilation of measured rates of change and also a methodology for objective comparison of rates. I present a compilation from 155 publications reporting glacier area changes, and also a methodology that overcomes the main obstacles hindering comparison. Glacier areas are not always assigned uncertainties, and this problem is addressed with an error model derived from published estimates. The problem of discordant survey dates is addressed by interpolating measured areas to fixed dates at pentadal intervals. Interpolation error depends only incoherently on the time span between measurements, but strongly on glacier size: smaller glaciers, in addition to changing more rapidly on average, exhibit more variable rates of change. The overlapping boundaries of study regions are reconciled by mapping all of the information to a 0.5° geographical grid. When coupled with glacier area information from the Randolph Glacier Inventory, the widely observed inverse dependence of shrinkage rates on glacier size shows promise as a tool for treating incomplete spatial coverage. Over High Mountain Asia as a whole from 1960 to 2010, the unweighted average shrinkage rate is –0.57% a–1, but corrections for variable glacier size raise the average to –0.34% a–1, and filling unmeasured gridcells with rates based on size dependence alters the latter estimate to –0.40% a–1. The uncertainties in these rates are large. The Karakoram anomaly is found to be a zonal feature extending well to the east of the Karakoram proper.

High-latitude atmospheric warming is impacting freshwater cycling, requiring techniques for monitoring the hydrology of sparsely-gauged regions. The submarine runoff of tidewater glaciers presents a particular challenge. We evaluate the utility of Moderate Resolution Imaging Spectroradiometer (MODIS) imagery for monitoring turbid meltwater plume variability in the glacier lagoon Jökulsárlón, Iceland, for a short interval before the onset of the main melt season. Total Suspended Solids concentrations (TSS) of surface waters are related to remotely-sensed reflectance via empirical calibration between in-situ-sampled TSS and reflectance in a MODIS band 1-equivalent wavelength window. This study differs from previous ones in its application to an overturning tidewater glacier plume, rather than one derived from river runoff. The linear calibration improves on previous studies by facilitating a wider range of plume metrics than areal extent, notably pixel-by-pixel TSS values. Increasing values of minimum plume TSS over the study interval credibly represent rising overall turbidity in the lagoon as melting accumulates. Plume extent responds principally to consistently-strong offshore winds. Further work is required to determine the temporal persistence of the calibration, but remote plume observation holds promise for monitoring hydrological outputs from ungauged or ungaugeable systems.

In Northeast Greenland, the Norske Øer Ice Barrier (NØIB) abuts Nioghalvfjerdsfjorden (79N) and Zachariae Isstrøm (ZI), two floating outlets of the Northeast Greenland Ice Stream. NØIB is an extensive region of perennially fast sea ice, which varies in size from year to year, but with complete breakup a rare event in the 20th century. It reportedly broke up in the 1950s and was seen to break up in 1997. Since 2000 the NØIB has broken up during 11 of the last 14 summers. The forcings driving the increased frequency of ice barrier breakup are poorly understood, and it is not clear if the breakup is a purely local phenomenon or an indicator of regional changes in the ocean and atmosphere. Here we use a logistic regression statistical model to show that the odds of breakup are linked to June positive degree days and July wind speeds at a nearby weather station. It is too soon to know if subtle changes detected on 79N and ZI in the last decade are connected to breakups of the NØIB but, if they are, it suggests a complex interaction between the atmosphere, ocean and outlet glaciers in this part of Greenland.

GPS campaigns on glaciers during jökulhlaups show how subglacial floods affect glacier motion and shed light on the dynamics of such floods. Three such campaigns have been carried out on southern and western Vatnajökull, southeast Iceland, over known jökulhlaup paths. Two slowly rising jökulhlaups from Grímsvötn and two rapidly rising jökulhlaups from the western and eastern Skaftá cauldrons were captured in these campaigns, with maximum discharge ranging from 240 to 3300 m3 s−1. Glacier surface movements measured in these campaigns are presented along with the corresponding discharge curves. The measurements are interpreted as indicating: (1) initiation of rapidly rising jökulhlaups with a propagating subglacial pressure wave, (2) decreased glacier basal friction during jökulhlaups, (3) subglacial accumulation of water in slowly rising jökulhlaups and (4) lifting of the glacier caused by subglacial water pressure exceeding overburden in both rapidly and slowly rising jökulhlaups. The latter two observations are inconsistent with assumptions that are typically made in theoretical and numerical modelling of jökulhlaups. Both viscous and elastic deformation of the glacier as well as turbulent hydraulic fracture at the ice/bedrock interface are important in the dynamics of the subglacial pressure wave at the front of rapidly rising jökulhlaups.

This study outlines a consistent methodology for identifying glacier surfaces from Landsat 5, 7 and 8 imagery that is applied to map all mainland North Asian glaciers, providing the first methodologically consistent and complete glacier inventory for the region ~2010. We identify 5065 glaciers covering a planimetric area of 2326 ± 186 km2, most of which is located in the Altai mountain subregion. The total glacier count is 15% higher, but the total glacier area is 32 ±11.6% lower, than the estimated glacier coverage provided in version 4.0 of the Randolph Glacier Inventory. We investigate the distribution of glacier size within North Asia and find that the majority of glaciers (82%) are smaller than 0.5 km2 but only account for a third of the total glacier area, with the largest 1 % (60 glaciers ≥ 5 km2 ) accounting for 28% of the total area. We present hypsometric characterizations of North Asian glaciers, largely substantiating existing findings that glaciers in this region are dominated by cold, relatively dry conditions. We provide a detailed assessment of errors and determine the uncertainty in our area estimate to be ±8.0%, with snow-cover uncertainty the largest contributing factor. Based on this assessment, the new glacier inventory presented here is more complete and of higher quality than other currently available data sources.

Glaciers in the Bugyai Kangri are located in a transition zone from southeast Tibet, where monsoonal temperate glaciers dominate, to inner Tibet, where continental glaciers dominate. Here we analyze glacier and glacial lake changes in this region using multi-year inventories based on Landsat images from 1981–2013. Results show that the total area of 141 glaciers in the region decreased by 30.44 ±0.89 km2 from 198.35 ±9.54 km2 (1980s) to 167.93 ±4.52 km2 (2010s). The annual area shrinkage rate (–0.48% a–1) is lower than that reported for southeastern Tibet but higher than that of inner Tibet. Both the number and total area of glacial lakes increased between 1981 and 2013. Among all lakes, proglacial lakes contribute most (~81 %) to the expansion. The total area of ten proglacial lakes increased by 150.3 ± 13.17% and of these ten lakes the four that expanded most sharply showed increased calving at their upper margins, resulting in more rapid retreat of lake-terminating glaciers than land-terminating glaciers. Owing to rapid calving, several lakes may undergo further growth in the near future, increasing the potential risk of glacial lake outburst floods.

A suite of surface and basal measurements during and after borehole drilling is used to perform in situ investigation of the local basal drainage system and pressure forcing in western Greenland. Drill and borehole water temperature were monitored during borehole drilling, which was performed with dyed hot water. After drilling, borehole water pressure and basal dye concentration were measured concurrently with positions in a GPS strain diamond at the surface. Water pressure exhibited diurnal changes in antiphase with velocity. Dye monitoring in the borehole revealed stagnant basal water for nearly 2 weeks. The interpretation of initial connection to an isolated basal cavity is corroborated by the thermal signature of borehole water during hot water drilling. Measurement-based estimates of cavity size are on the order of cubic meters, and analysis indicates that small changes in its volume could induce the observed pressure variations. It is found that longitudinal coupling effects are unable to force necessary volume changes at the site. Sliding-driven basal cavity opening and elastic uplift from load transfer are plausible mechanisms controlling pressure variations. Elastic uplift requires forcing from a hydraulically connected reach, which observations suggest must be relatively small and in close proximity to the isolated cavity.