Subglacial hydrologic systems in Antarctica and Greenland play a fundamental role in ice-sheet dynamics, yet critical aspects of these systems remain poorly understood due to a lack of observations. Ground-based electromagnetic (EM) geophysical methods are established for mapping groundwater in many environments, but have never been applied to imaging lakes beneath ice sheets. Here, we study the feasibility of passive- and active-source EM imaging for quantifying the nature of subglacial water systems beneath ice streams, with an emphasis on the interfaces between ice and basal meltwater, as well as deeper groundwater in the underlying sediments. We describe a suite of model studies that exam the data sensitivity as a function of ice thickness, water conductivity and hydrologic system geometry for models representative of a subglacial lake and a grounding zone estuary. We show that EM data are directly sensitive to groundwater and can image its lateral and depth extent. By combining the conductivity obtained from EM data with ice thickness and geological structure from conventional geophysical techniques, such as ground-penetrating radar and active seismic surveying, EM data have the potential to provide new insights on the interaction between ice, rock and water at critical ice-sheet boundaries.

The extensive dataset of Variegated Glacier's 1982/83 surge and antecedent quiescent phase is fragmentary and predominantly of flow line nature. Applying the raw centre line data in a flow line model to conduct research into the mechanism behind Variegated Glacier's surge behaviour inevitably entails problems. In this study, the incomplete dataset is extrapolated into the upper and lower parts of the glacier. Furthermore, the centre line data are adapted to account for differences between width-averaged and centre line surface mass balance and ice thickness change, inflow from tributaries, and changes in surface width over time. The modifications to the dataset are backed-up by observations and clear, plausible physical explanations. Moreover, the revised dataset meets the imposed constraints regarding ice volume flux, specific mass balance and net volume change. Hence, the final dataset is considered a satisfying revision that makes the dataset more valuable for future research. Subsequent application of the revised dataset corroborates the idea that glacier evolution during quiescence is basically the growth back towards steady state after the glacier was brought out of balance by the preceding surge.

Tidal flexure in ice shelf grounding zones has been used extensively in the past to determine grounding line position and ice properties. Although the rheology of ice is viscoelastic at tidal loading frequencies, most modelling studies have assumed some form of linear elastic beam approximation to match observed flexure profiles. Here we use density, radar and DInSAR measurements in combination with full-Stokes viscoelastic modelling to investigate a range of additional controls on the flexure of the Southern McMurdo Ice Shelf. We find that inclusion of observed basal crevasses and density dependent ice stiffness can greatly alter the flexure profile and yet fitting a simple elastic beam model to that profile will still produce an excellent fit. Estimates of the effective Young's modulus derived by fitting flexure profiles are shown to vary by over 200% depending on whether these factors are included, even when the local thickness is well constrained. Conversely, estimates of the grounding line position are relatively insensitive to these considerations for the case of a steep bed slope in our study region. By fitting tidal amplitudes only, and ignoring phase information, elastic beam theory can provide a good fit to observations in a wide variety of situations. This should, however, not be taken as an indication that the underlying rheological assumptions are correct.

We describe the use of a terrestrial laser scanner (TLS) to monitor the net mass balance of Urumqi Glacier No. 1, eastern Tien Shan. We used an ultra-long-range Riegl VZ®-6000 TLS, which is specially designed for surveying snow- and ice-covered terrain, to create repeated high spatiotemporal resolution DEMs, focusing on the monthly-scale (25 April–28 May 2015) net mass balance. According to the TLS-derived DEMs, the area of Urumqi Glacier No. 1 was 1.558 km2 on 25 April 2015 and the average surface elevation change was 0.225 m. By comparing the results from the use of TLS with the conventional glaciological mass-balance method, the correlation coefficient (R2) between glaciological elevation changes of individual stakes and the TLS-derived geodetic elevation change of corresponding points was 0.85. Considering the uncertainty of both methods, this is a promising result. Using the in situ measured snow densities (snow pits) of the glacier surface, the geodetic net mass balance was 0.074 m w.e., which is slightly positive. The mean uncertainty in the TLS-derived monthly net mass balance was 0.018 m w.e., showing that the TLS surveying system presented accurate and relevant results and is therefore suitable to monitor mass-balance evolution of mountain glaciers.

Previous estimates of melt and surface sublimation on glaciers of the subtropical semiarid Andes (29–34°S) have been obtained at few specific locations, but it is not clear how ablation components vary across the entire extent of a glacier in this dry environment. Here, we simulate the distributed energy and mass balance of Juncal Norte Glacier (33°S) during a 2-month summer period. Forcing fields of near-surface air temperature and wind speed are generated using two methods accounting for the main physical processes that shape their spatial variations. Simulated meteorological variables and ablation agree well with observations on the glacier tongue and reveal complex patterns of energy and mass fluxes. Ablation decreases from 70 mm w.e. d−1 at the low-albedo glacier terminus (~3000 m), where almost 100% of total ablation corresponds to melt, to <5 mm w.e. d−1 at wind-exposed, strong-radiated sites above 5500 m, where surface sublimation represents >75% of total ablation. Our simulations provide the first glacier-scale estimates of ablation components on a glacier in the study region and better reproduce the observed and expected spatial variations of melt and surface sublimation, in comparison with more simple assumptions, such as linear gradients and uniform wind speeds.

Measurements of glacier ice cliff evolution are sparse, but where they do exist, they indicate that such areas of exposed ice contribute a disproportionate amount of melt to the glacier ablation budget. We used Structure from Motion photogrammetry with Multi-View Stereo to derive 3-D point clouds for nine ice cliffs on Khumbu Glacier, Nepal (in November 2015, May 2016 and October 2016). By differencing these clouds, we could quantify the magnitude, seasonality and spatial variability of ice cliff retreat. Mean retreat rates of 0.30–1.49 cm d−1 were observed during the winter interval (November 2015–May 2016) and 0.74–5.18 cm d−1 were observed during the summer (May 2016–October 2016). Four ice cliffs, which all featured supraglacial ponds, persisted over the full study period. In contrast, ice cliffs without a pond or with a steep back-slope degraded over the same period. The rate of thermo-erosional undercutting was over double that of subaerial retreat. Overall, 3-D topographic differencing allowed an improved process-based understanding of cliff evolution and cliff-pond coupling, which will become increasingly important for monitoring and modelling the evolution of thinning debris-covered glaciers.

We performed an Empirical Orthogonal Function (EOF) analysis to assess the intraseasonal variability of 5–60 day band-pass filtered Antarctic sea-ice concentration in austral winter using a 20-year daily dataset from 1995 to 2014. Zonal wave number 3 dominated in the Antarctic, especially so across the west Antarctic. Results showed the coexistence of stationary and propagating wave components. A spectral analysis of the first two principal components (PCs) showed a similar structure for periods up to 15 days but generally more power in PC1 at longer periods. Regression analysis upon atmospheric fields using the first two PCs of sea-ice concentration showed a coherent wave number 3 pattern. The spatial phase delay between the sea-ice and mean sea-level pressure patterns suggests that meridional flow and associated temperature advection are important for modulating the sea-ice field. EOF analyses carried out separately for El Niño, La Niña and neutral years, and for Southern Annular Mode positive, negative and neutral periods, suggest that the spatial patterns of wave number 3 shift between subsets. The results also indicate that El Niño-Southern Oscillation and Southern Annular Mode affect stationary wave interactions between sea-ice and atmospheric fields on intraseasonal timescales.

We examine repeat surface altimetry and radio echo observations of two supraglacial lakes in the percolation zone of the Greenland ice sheet to investigate the changes in firn conditions leading to lake formation and implications for meltwater storage within firn. Both lakes formed in 2011, when an anomalously high melt season was followed by low winter accumulation, resulting in reduced infiltration and storage in the near surface. The lakes expanded during the 2012 record melt season and retained liquid meltwater through the following winter. The lakes then contracted, with one lake slowly draining and refreezing and another rapidly draining to the subsurface. The lack of observable change in firn conditions surrounding the lakes indicates increased run-off in the near surface firn, likely along low-permeability ice layers formed during the previous melt seasons. This implies a reduced ability of the firn to absorb increased meltwater.

Basal slip accounts for a large part of the flow of ice streams draining ice from Antarctica and Greenland into the ocean. Therefore, an appropriate representation of basal slip in ice flow models is a prerequisite for accurate sea level rise projections. Various friction laws have been proposed to describe basal slip in models. Here, we compare the influence on grounding line (GL) dynamics of four friction laws: the traditional Weertman law and three effective pressure-dependent laws, namely the Schoof, Tsai and Budd laws. It turns out that, even when they are tuned to a common initial reference state, the Weertman, Budd and Schoof laws lead to thoroughly different steady-state positions, although the Schoof and Tsai laws lead to much the same result. In particular, under certain circumstances, it is possible to obtain a steady GL located on a reverse slope area using the Weertman law. Furthermore, the predicted transient evolution of the GL as well as the projected contributions to sea level rise over a 100-year time horizon vary significantly depending on the friction law. We conclude on the importance of choosing an appropriate law for reliable sea level rise projections and emphasise the need for a coupling between ice flow models and physically based subglacial hydrological models.

Glacial earthquakes are slow earthquakes of magnitude M~5 associated with major calving events at near-grounded marine-terminating glaciers. These globally detectable earthquakes provide information on the grounding state of outlet glaciers and the timing of large calving events. Seismic source modeling of glacial earthquakes provides information on the size and orientation of forces associated with calving events. We compare force orientations estimated using a centroid-single-force technique with the calving-front orientations of the source glaciers at or near the time of earthquake occurrence. We consider earthquakes recorded at four glaciers in Greenland – Kangerdlugssuaq Glacier, Helheim Glacier, Kong Oscar Glacier, and Jakobshavn Isbræ – between 1999 and 2010. We find that the estimated earthquake force orientations accurately represent the orientation of the calving front at the time of the earthquake, and that seismogenic calving events are produced by a preferred section of the calving front, which may change with time. We also find that estimated earthquake locations vary in a manner consistent with changes in calving-front position, though with large scatter. We conclude that changes in glacial-earthquake source parameters reflect true changes in the geometry of the source glaciers, providing a means for identifying changes in glacier geometry and dynamics that complements traditional remote-sensing techniques.

The albedo of seasonal snow cover plays an important role in the global climate system due to its influence on Earth's radiation budget and energy balance. Volunteer Community Collaborative Rain, Hail, and Snow-Albedo (CoCoRaHS-Albedo) observers collected 3249 individual daily albedo, snow depth and density measurements using standardized techniques at dozens of sites across New Hampshire, USA over four winter seasons. The data show that albedo increases rapidly with snow depth up to ~0.14 m. Multiple linear regression models using snowpack age, snow depth or density, and air temperature provide reasonable approximations of surface snow albedo during times of albedo decay. However, the linear models also reveal systematic biases that highlight an important non-linearity in snow albedo decay. Modeled albedo values are reasonably accurate within the range of 0.6–0.9, but exhibit a tendency to overestimate lower albedo values and underestimate higher albedo values. We hypothesize that rapid reduction in high albedo fresh snow results from a decrease in snow specific surface area, while during melt-events the presence of liquid water in the snowpack accelerates metamorphism and grain growth. We conclude that the CoCoRaHS-Albedo volunteer observer network provides useful snow albedo, depth and density measurements and serves as an effective model for future measurement campaigns.

Time series ice-draft data were obtained from moored ice-profiling sonar (IPS), in the coastal northeastern Chukchi Sea during 2009/10. Time series data show seasonal growth of sea-ice draft, occasionally interrupted by coastal polynya. The sea-ice draft distribution indicates a slightly lower abundance of thick, deformed ice compared with the eastern Beaufort Sea. In January, a rapid increase in the abundance of thick ice coincided with a period of minimal drift indicating compaction again the coast and dynamical thickening. The overall mean draft and corresponding derived thickness are 1.27 and 1.38 m, respectively. The evolution of modal ice thickness observed can be explained mostly by thermodynamic growth. The derived ice thicknesses are used to estimate heat losses based on ERA-interim data. Heat losses from the raw, 1 s IPS data are ~50 and 100% greater than those calculated using IPS data averaged over spatial scales of ~20 and 100 km, respectively. This finding demonstrates the importance of subgrid-scale ice-thickness distribution for heat-loss calculation. The heat-loss estimate based on thin ice data derived from AMSR-E data corresponds well with that from the 1 s observed ice-thickness data, validating heat-loss estimates from the AMSR-E thin ice-thickness algorithm.

Glacier albedo controls the surface energy budget, the variability of which affects the glacier surface melt rate and, in turn, impacts the mass balance of the glacier. During 2013 and 2014, spatial and temporal variations of albedo were investigated using 18 Landsat images of Urumqi Glacier No. 1. Factors influencing these spatiotemporal profiles were analyzed. An established retrieval process, including geolocation, radiometric calibration, atmospheric, topographic, and anisotropic correction and narrow- to broadband conversion, was applied for the first time to Landsat-8 images. Differences between Landsat image derived albedo values and albedo values measured using a handheld spectroradiometer ranged from −0.024 to 0.049. Spatial and temporal variations of surface albedo were significant, especially in the ablation area. The variability of the values of ice albedo ranged from 0.06 to 0.44 due to topographic effects and light-absorbing impurities. The results suggest that this retrieval method can be used to investigate the spatial and temporal variability of surface albedo from Landsat-8 images on mountain glaciers. Moreover, as constant albedo values for ice and snow cannot be assumed, the distribution of albedo was not completely dependent on altitude under conditions of more intense ablation, and by reason of light-absorbing impurities during the melt season.

We analyse ice flow and structural glaciology of Shackleton Glacier, a dendritic glacier with multiple icefalls in the Canadian Rockies. A major tributary-trunk junction allows us to investigate the potential of tributaries to alter trunk flow and structure, and the formation of bedrock steps at confluences. Multi-year velocity-stake data and structural glaciology up-glacier from the junction were assimilated with glacier-wide velocity derived from Radarsat-2 speckle tracking. Maximum flow speeds are 65 m a−1 in the trunk and 175 m a−1 in icefalls. Field and remote-sensing velocities are in good agreement, except where velocity gradients are high. Although compression occurs in the trunk up-glacier of the tributary entrance, glacier flux is steady state because flow speed increases at the junction due to the funnelling of trunk ice towards an icefall related to a bedrock step. Drawing on a published erosion model, we relate the heights of the step and the hanging valley to the relative fluxes of the tributary and trunk. It is the first time that an extant glacier is used to test and support such model. Our study elucidates the inherent complexity of tributary/trunk interactions and provides a conceptual model for trunk flow restriction by a tributary in surge-type glaciers.

Radioisotopes (239Pu, 240Pu, 236U) formed during atmospheric nuclear weapons testing (NWT) can be used for dating and therefore be applied as markers of the beginning of the Anthropocene Epoch. Moreover, 240Pu/239Pu ratios enable source identification (general fallout, local emission sources). Ice core segments from a 57.6 m ice core to bedrock from eastern Tien Shan, China were selected for 239,240Pu and 236U analyses by using compact low-energy accelerator mass spectrometry. The observed 240Pu/239Pu atom ratios were 0.18 ± 0.02, with one exception, indicating global fallout. No evidence for emissions from the nearby local sources Lop Nor was observed. The total deposition rates during NWT for 239Pu and 240Pu amount to 1.55 × 109 atoms·cm−2 and 3.1 × 108 atoms·cm−2, respectively. With the higher undisturbed fallout of 239Pu compared with values reported for glaciers from European Alps at Col du Dome and Colle Gnifetti as well as from the polar region, the eastern Tien Shan glacier would be an ideal site for defining the start of the Anthropocene. 236U is a rather new environmental tracer, while little is known about its global fallout from NWT. The observed deposition flux of 236U during NWT amounts to 3.5 × 108 atoms·cm−2 at the drill site. The average 236U/239Pu ratio of 0.27 ± 0.09 is in good agreement with literature values from global fallout.