Engabreen is an outlet glacier of the Svartisen Ice Cap located in Northern Norway. It is a unique glacier due to the Svartisen Subglacial Laboratory which allows direct access to the glacier bed. In this study, we combine both sub- and supraglacial observations with ice-flow modelling in order to investigate conditions at the bed of Engabreen both spatially and temporally. We use the full-Stokes model Elmer/Ice and satellite-based surface-velocity maps from 2010 and 2014 to infer patterns of basal friction. Direct measurements of basal sliding and deformation of lower layers of the ice are used to adjust the ice viscosity and provide essential input to the setup of our model and influence the interpretation of the results. We find a clear seasonal cycle in the subglacial conditions at the higher elevation region of the study area and discuss this in relation to the subglacial hydrological system. Our results also reveal an area with an overdeepening where basal friction is significantly lower than elsewhere on the glacier all year round. We attribute this to either water pooling at the base, or saturated sediments and increased strain heating at this location which softens the ice further.

The detection and monitoring of meltwater within firn presents a significant monitoring challenge. We explore the potential of small wireless sensors (ETracer+, ET+) to measure temperature, pressure, electrical conductivity and thus the presence or absence of meltwater within firn, through tests in the dry snow zone at the East Greenland Ice Core Project site. The tested sensor platforms are small, robust and low cost, and communicate data via a VHF radio link to surface receivers. The sensors were deployed in low-temperature firn at the centre and shear margins of an ice stream for 4 weeks, and a ‘bucket experiment’ was used to test the detection of water within otherwise dry firn. The tests showed the ET+ could log subsurface temperatures and transmit the recorded data through up to 150 m dry firn. Two VHF receivers were tested: an autonomous phase-sensitive radio-echo sounder (ApRES) and a WinRadio. The ApRES can combine high-resolution imaging of the firn layers (by radio-echo sounding) with in situ measurements from the sensors, to build up a high spatial and temporal resolution picture of the subsurface. These results indicate that wireless sensors have great potential for long-term monitoring of firn processes.

Ice shelves play an important role in buttressing land ice from reaching the sea, thus restraining the rate of grounded ice loss. Long-period gravity-wave impacts excite vibrations in ice shelves that can expand pre-existing fractures and trigger iceberg calving. To investigate the spatial amplitude variability and propagation characteristics of these vibrations, a 34-station broadband seismic array was deployed on the Ross Ice Shelf (RIS) from November 2014 to November 2016. Two types of ice-shelf plate waves were identified with beamforming: flexural-gravity waves and extensional Lamb waves. Below 20 mHz, flexural-gravity waves dominate coherent signals across the array and propagate landward from the ice front at close to shallow-water gravity-wave speeds (~70 m s−1). In the 20–100 mHz band, extensional Lamb waves dominate and propagate at phase speeds ~3 km s−1. Flexural-gravity and extensional Lamb waves were also observed by a 5-station broadband seismic array deployed on the Pine Island Glacier (PIG) ice shelf from January 2012 to December 2013, with flexural wave energy, also detected at the PIG in the 20–100 mHz band. Considering the ubiquitous presence of storm activity in the Southern Ocean and the similar observations at both the RIS and the PIG ice shelves, it is likely that most, if not all, West Antarctic ice shelves are subjected to similar gravity-wave excitation.

Increasing ship traffic and human activity in the Arctic has led to a growing demand for accurate Arctic weather forecast. High-quality forecasts obtained by models are dependent on accurate initial states achieved by assimilation of observations. In this study, a multi-variate nudging (MVN) method for assimilation of sea-ice variables is introduced. The MVN assimilation method includes procedures for multivariate update of sea-ice volume and concentration, and for extrapolation of observational information spatially. The MVN assimilation scheme is compared with the Ensemble Kalman Filter (EnKF) using the Los Alamos Sea Ice Model. Two multi-variate experiments are conducted: in the first experiment, sea-ice thickness from the European Space Agency's Soil Moisture and Ocean Salinity mission is assimilated, and in the second experiment, sea-ice concentration from the ocean and Sea Ice Satellite Application Facility is assimilated. The multivariate effects are cross-validated by comparing the model with non-assimilated observations. It is found that the simple and computationally cheap MVN method shows comparable skills to the more complicated and expensive EnKF method for multivariate update. In addition, we show that when few observations are available, the MVN method is a significant model improvement compared to the version based on one-dimensional sea-ice concentration assimilation.

We simulate the ice dynamics of the San Rafael Glacier (SRG) in the Northern Patagonia Icefield (46.7°S, 73.5°W), using glacier geometry obtained by airborne gravity measurements. The full-Stokes ice flow model (Elmer/Ice) is initialized using an inverse method to infer the basal friction coefficient from a satellite-derived surface velocity mosaic. The high surface velocities (7.6 km a−1) near the glacier front are explained by low basal shear stresses (<25 kPa). The modelling results suggest that 98% of the surface velocities are due to basal sliding in the fast-flowing glacier tongue (>1 km a−1). We force the model using different surface mass-balance scenarios taken or adapted from previous studies and geodetic elevation changes between 2000 and 2012. Our results suggest that previous estimates of average surface mass balance over the entire glacier (Ḃ) were likely too high, mainly due to an overestimation in the accumulation area. We propose that most of SRG imbalance is due to the large ice discharge (−0.83 ± 0.08 Gt a−1) and a slightly positive Ḃ (0.08 ± 0.06 Gt a−1). The committed mass-loss estimate over the next century is −0.34 ± 0.03 Gt a−1. This study demonstrates that surface mass-balance estimates and glacier wastage projections can be improved using a physically based ice flow model.

The positive degree-day (PDD) model provides a particularly simple approach to estimate surface melt from land ice based solely on air temperature. Here, we use a climate and snow pack simulation of the Greenland ice sheet (Modèle Atmosphérique Régional, MAR) as a reference, to analyze this scheme in three realizations that incorporate the sub-monthly temperature variability differently: (i) by local values, (ii) by local values that systematically overestimate the dampened variability associated with intense melting or (iii) by one constant value. Local calibrations reveal that incorporating local temperature variability, particularly resolving the dampened variability of melt areas, renders model parameters more temperature-dependent. This indicates that the negative feedback between surface melt and temperature variability introduces a non-linearity into the temperature – melt relation. To assess the skill of the individual realizations, we hindcast melt rates from MAR temperatures for each realization. For this purpose, we globally calibrate Greenland-wide, constant parameters. Realization (i) exhibits shortcomings in the spatial representation of surface melt unless temperature-dependent instead of constant parameters are calibrated. The other realizations perform comparatively well with constant parametrizations. The skill of the PDD model primarily depends, however, on the consistent calibration rather than on the specific representation of variability.

Tropical glaciers constitute an important source of water for downstream populations. However, our understanding of glacial melt processes is still limited. One observed process that has not yet been quantified for tropical glaciers is the enhanced melt caused by the longwave emission transfer. Here, we use high-resolution surface temperatures obtained from the thermal infrared imagery of the Cuchillacocha Glacier, in the Cordillera Blanca, Peru in June 2014 to calculate a margin longwave flux. This longwave flux, reaching the glacier margin from the adjacent exposed rock, varies between 81 and 120 W m−2 daily. This flux is incorporated into a physically-based melt model to assess the net radiation budget at the modeled glacier margin. The simulation results show an increase in the energy available for melt by an average of 106 W m−2 during the day when compared with the simulation where the LWmargin flux is not accounted for. This value represents an increase in ablation of ~1.7 m at the glacier margin for the duration of the dry season. This study suggests that including the quantification of the glacier margin longwave flux in physically-based melt models results in an improved assessment of tropical glacier energy budget and meltwater generation.

Twin glaciers collapsed in 2016 near Aru Co, western Tibet and caused extreme loss to human beings. In this study, we attempted to track the dynamics of glaciers in the region, for example the glacier area and mass changes in Aru Co for the period 1971–2016, which were determined using topographic maps and Landsat images and ASTER-derived DEMs (2011–16), the Shuttle Radar Terrain Mission DEM (2000) and topographic maps (1971). Our results showed that the glacier area of Aru Co decreased by −0.4 ± 4.1% during 1971–2016. The geodetic mass-balance results showed that the glaciers in Aru Co lost mass at a rate of −0.15 ± 0.30 m w.e. a−1 during 1971–99, while they gained mass at a rate of 0.33 ± 0.61 m w.e. a−1 for the period 1999–2016. The twin glaciers experienced a larger negative mass budget than the others in the region before 1999. This process produced large amounts of meltwater, followed by a sustained increase in the meltwater on the pressure melting point, possibly in response to a period of positive mass balance (1999–2016) and then, transferred to the glacier bed until the glaciers collapsed.

We have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

Subglacial hydrology plays a key role in many glaciological processes, including ice dynamics via the modulation of basal sliding. Owing to the lack of an overarching theory, however, a variety of model approximations exist to represent the subglacial drainage system. The Subglacial Hydrology Model Intercomparison Project (SHMIP) provides a set of synthetic experiments to compare existing and future models. We present the results from 13 participating models with a focus on effective pressure and discharge. For many applications (e.g. steady states and annual variations, low input scenarios) a simple model, such as an inefficient-system-only model, a flowline or lumped model, or a porous-layer model provides results comparable to those of more complex models. However, when studying short term (e.g. diurnal) variations of the water pressure, the use of a two-dimensional model incorporating physical representations of both efficient and inefficient drainage systems yields results that are significantly different from those of simpler models and should be preferentially applied. The results also emphasise the role of water storage in the response of water pressure to transient recharge. Finally, we find that the localisation of moulins has a limited impact except in regions of sparse moulin density.

Ice rheology governs how glaciers flow and respond to environmental change. The rheology of glacier ice evolves in response to a variety of mechanisms, including damage, heating, melting and the development of crystalline fabric. The relative contributions of these rheological mechanisms are not well understood. Using remotely sensed data and physical models, we decouple the influence of each of the aforementioned mechanisms along the margins of Rutford Ice Stream, a laterally confined outlet glacier in West Antarctica. We show that fabric is an important control on ice rheology in the shear margins, with an inferred softening effect consistent with a single-maximum fabric. Fabric evolves to steady state near the onset of streaming flow, and ice progressively softens downstream almost exclusively due to shear heating. The rate of heating is sensitive to local shear strain rates, which respond to local changes in bed topography as ice is squeezed through the basal trough. The impact of shear heating on the downstream evolution of ice rheology in a laterally confined glacier suggests that the thermoviscous feedback – wherein faster ice flow leads to higher rates of shear heating, further softening the ice – is a fundamental control on glacier dynamics.

Thermal and hydrological dynamics and their impacts on the stability of a moraine dam were analyzed and simulated for the Longbasaba Lake in the Himalaya, based on soil temperature, moisture and heat flux data observed at different depths in the dam from 2012 to 2016. Annual average heat income is greater than heat expenditure on the dam surface. The mean annual temperature at observed the depths of 0–150 cm is >0°C, although the average annual air temperature was −3.6°C over the dam, indicating a relatively larger temperature difference between moraine dam and air. The volumetric soil moisture content is relatively low with an annual average of 5%, peaking after the snow cover melting and active layer thawing. Simulation results indicate that the average yearly maximum thawing depth has been ~0.3 m deeper than the average yearly maximum freezing depth during the observation period. In the past 55 years, the yearly maximum thawing depth has increased, while yearly maximum freezing depth has decreased, implying that the permafrost in the dam has been deteriorating. The annual surplus heat and increasing permafrost thawing depth will result in further deterioration of permafrost and melting of buried ice in the dam, thereby decreasing its stability.

Understanding the effects of climate on glaciers requires precise estimates of ice volume change over several decades. This is achieved by the geodetic mass balance computed by two means: (1) the digital elevation model (DEM) comparison (SeqDEM) allows measurements over the entire glacier, however the low contrast over glacierized areas is an issue for the DEM generation through the photogrammetric techniques and (2) the profiling method (SePM) is a faster alternative but fails to capture the spatial variability of elevation changes. We present a new framework (SSD) that relies upon the spatial variability of the elevation change to densify a sampling network to optimize the surface-elevation change quantification. Our method was tested in two small glaciers over different periods. We conclude that the SePM overestimates the elevation change by ~20% with a mean difference of ~1.00 m (root mean square error (RMSE) = ~3.00 m) compared with results from the SeqDEM method. A variogram analysis of the elevation changes showed a mean difference of <0.10 m (RMSE = ~2.40 m) with SSD approach. A final assessment on the largest glacier in the French Alps confirms the high potential of our method to compute the geodetic mass balance, without going through the generation of a full-density DEM, but with a similar accuracy than the SeqDEM approach.

We investigate the appropriateness of calving or crevasse models from the literature using linear elastic fracture mechanics (LEFM). To this end, we compare LEFM model-predicted stress intensity factors (SIFs) against numerically computed SIFs using the displacement correlation method in conjunction with the finite element method. We present several benchmark simulations wherein we calculate the SIF at the tips of water-filled surface and basal crevasses penetrating through rectangular ice slabs under different boundary conditions, including grounded and floating conditions. Our simulation results indicate that the basal boundary condition significantly influences the SIF at the crevasse tips. We find that the existing calving models using LEFM are not generally accurate for evaluating SIFs in grounded glaciers or floating ice shelves. We also illustrate that using the ‘single edge crack’ weight function in the LEFM formulations may be appropriate for predicting calving from floating ice shelves, owing to the low fracture toughness of ice; whereas, using the ‘double edge crack’ or ‘central through crack’ weight functions is more appropriate for predicting calving from grounded glaciers. To conclude, we recommend using the displacement correlation method for SIF evaluation in real glaciers and ice shelves with complex geometries and boundary conditions.

Band ogives are a striking and enigmatic feature of Mer de Glace glacier flow. The surface mass balances (SMBs) of these ogives have been thoroughly investigated over a period of 12 years. We find similar cumulative SMBs over this period, ranging between −64.1 and −66.2 m w.e., on the dark and light ogives even though the dark ogive albedo is ~40% lower than that of the light ogives. We, therefore, looked for another process that could compensate for the large difference of absorbed short-wave radiation between dark and light ogives. Based on in situ roughness measurements, our numerical modeling experiments demonstrate that a significant difference in turbulent flux over the dark and light ogives due to different surface roughnesses could compensate for the difference in radiative forcing. Our results discard theories for the genesis of band ogives that are based on the assumption of a strong ice ablation contrast between dark and light ogives. More generally, our study demonstrates that future roughness changes are as important to analyze as the radiative impacts of a potential increase of aerosols or debris at the surface of glaciers.

We present a comprehensive review of the status and changes in glacier length (since the 1850s), area and mass (since the 1960s) along the Himalayan-Karakoram (HK) region and their climate-change context. A quantitative reliability classification of the field-based mass-balance series is developed. Glaciological mass balances agree better with remotely sensed balances when we make an objective, systematic exclusion of likely flawed mass-balance series. The Himalayan mean glaciological mass budget was similar to the global average until 2000, and likely less negative after 2000. Mass wastage in the Himalaya resulted in increasing debris cover, the growth of glacial lakes and possibly decreasing ice velocities. Geodetic measurements indicate nearly balanced mass budgets for Karakoram glaciers since the 1970s, consistent with the unchanged extent of supraglacial debris-cover. Himalayan glaciers seem to be sensitive to precipitation partly through the albedo feedback on the short-wave radiation balance. Melt contributions from HK glaciers should increase until 2050 and then decrease, though a wide range of present-day area and volume estimates propagates large uncertainties in the future runoff. This review reflects an increasing understanding of HK glaciers and highlights the remaining challenges.

We use ice flow modelling to simulate the englacial stratigraphy of Blåskimen Island, an ice rise in Dronning Maud Land and elucidate the evolution of this data-sparse region. We apply a thermo-mechanically coupled Elmer/Ice model to a profile along flowlines and through the ice-rise summit, where surface mass balance (SMB), flow velocity and ice stratigraphy were recently measured. We conclude that: (i) the ice rise is presently thickening at a rate of 0.5~0.6 m ice equivalent per year (mieq a−1), which is twice an earlier estimate using the field data and the input–output method; (ii) present thickening started 20–40 years in the past, before which the ice rise was in a steady state; (iii) SMB contrast between the upwind and downwind slopes was stronger than the present value by ~23% (or 0.15 mieq a−1) prior to ~1100 years ago. Since then, this contrast has been decreasing overall. We surmise that these SMB changes are likely a result of synoptic-scale atmospheric changes, rather than local atmospheric changes controlled by local ice topography. Our technique effectively assimilates geophysical data, avoiding the complexity of ice flow beneath the ice divide. Thus, it could be applied to other ice rises to elucidate the recent glacial retreat.

We have mapped the full crystallographic orientation of sea ice using electron backscatter diffraction (EBSD). This is the first time EBSD has been used to study sea ice. Platelet ice is a feature of sea ice near ice shelves. Ice crystals accumulate as an unconsolidated sub-ice platelet layer beneath the columnar ice (CI), where they are subsumed by the advancing sea–ice interface to form incorporated platelet ice (PI). As is well known, in CI the crystal preferred orientation comprises dominantly horizontal c-axes, while PI has c-axes varying between horizontal and vertical. For the first time, this study shows the a-axes of CI and PI are not random. Misorientation analysis has been used to illuminate the possible drivers of these alignments. In CI the misorientation angle distribution from random pairs and neighbour pairs of grains are indistinguishable, indicating the distributions are a consequence of crystal preferred orientation. Geometric selection during growth will develop the a-axis alignment in CI if ice growth in water is fastest parallel to the a-axis, as has previously been hypothesised. In contrast, in PI random-pair and neighbour-pair misorientation distributions are significantly different, suggesting mechanical rotation of crystals at grain boundaries as the most likely explanation.

Determining the position and stability of the grounding line of a marine ice sheet is a major challenge for ice-sheet models. Here, we investigate the role of lateral shear and ice-shelf buttressing in grounding line dynamics by extending an existing boundary layer theory to laterally confined marine ice sheets. We derive an analytic expression for the ice flux at the grounding line of confined marine ice sheets that depends on both local bed properties and non-local ice-shelf properties. Application of these results to a laterally confined version of the MISMIP 1a experiment shows that the boundary condition at the ice-shelf front (i.e. the calving law) is a major control on the location and stability of the grounding line in the presence of buttressing, allowing for both stable and unstable grounding line positions on downwards sloping beds. These results corroborate the findings of existing numerical studies that the stability of confined marine ice sheets is influenced by ice-shelf properties, in contrast to unconfined configurations where grounding line stability is solely determined by the local slope of the bed. Consequently, the marine ice-sheet instability hypothesis may not apply to buttressed marine ice sheets.

Simple models of glacier volume evolution are important in understanding features of glacier response to climate change, due to the scarcity of data adequate for running more complex models on a global scale. Two quantities of interest in a glacier's response to climate changes are its response time and its volume sensitivity to changes in the equilibrium line altitude (ELA). We derive a simplified, computationally inexpensive model of glacier volume evolution based on a block model with volume–area–length scaling. After analyzing its steady-state properties, we apply the model to each mountain glacier worldwide and estimate regionally differentiated response times and sensitivities to ELA changes. We use a statistical method from the family of global sensitivity analysis methods to determine the glacier quantities, geometric and climatic, that most influence the model output. The response time is dominated by the climatic setting reflected in the mass-balance gradient in the ablation zone, followed by the surface slope, while volume sensitivity is mainly affected by glacier size, followed by the surface slope.