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.

Versatile instruments assembled from off-the-shelf sensors and open-source electronics are used to record wave propagation and damping measured by Inertial Motion Units (IMUs) in a grease ice slick near the shore in Adventfjorden, Svalbard. Viscous attenuation of waves due to the grease ice slick is clearly visible by comparing the IMU data recorded by the different instruments. The frequency dependent spatial damping of the waves is computed by comparing the power spectral density obtained from the different IMUs. We model wave attenuation using the one-layer model of Weber from 1987. The best-fit value for the effective viscosity is ν = (0.95 ± 0.05 × 10−2)m2 s−1, and the coefficient of determination is R2 = 0.89. The mean absolute error and RMSE of the damping coefficient are 0.037 and 0.044m−1, respectively. These results provide continued support for improving instrument design for recording wave propagation in ice-covered regions, which is necessary to this area of research as many authors have underlined the need for more field data.

Ice-stream dynamics are strongly controlled by processes taking place at the ice/bed interface where subglacial water both lubricates the base and saturates any existing, underlying sediment. Large parts of the former Eurasian ice sheet were underlain by thick sequences of soft, marine sediments and many areas are imprinted with geomorphological features indicative of fast flow and wet basal conditions. Here, we study the effect of subglacial water on past Eurasian ice-sheet dynamics by incorporating a thin-film model of basal water flow into the ice-sheet model SICOPOLIS and use it to better represent flow in temperate areas. The adjunction of subglacial hydrology results in a smaller ice-sheet building up over time and generally faster ice velocities, which consequently reduces the total area fraction of temperate basal ice and ice streaming areas. Minima in the hydraulic pressure potential, governing water flow, are used as indicators for potential locations of past subglacial lakes and a probability distribution of lake existence is presented based on estimated lake depth and longevity.

We show results from a positive degree-day (PDD) model of Greenland ice sheet (GrIS) surface mass balance (SMB), 1870–2012, forced with reanalysis data. The model includes an improved daily temperature parameterization as compared with a previous version and is run at 1 km rather than 5 km resolution. The improvements lead overall to higher SMB with the same forcing data. We also compare our model with results from two regional climate models (RCMs). While there is good qualitative agreement between our PDD model and the RCMs, it usually results in lower precipitation and lower runoff but approximately equivalent SMB: mean 1979–2012 SMB (± standard deviation), in Gt a−1, is 382 ± 78 in the PDD model, compared with 379 ± 101 and 425 ± 90 for the RCMs. Comparison with in situ SMB observations suggests that the RCMs may be more accurate than PDD at local level, in some areas, although the latter generally compares well. Dividing the GrIS into seven drainage basins we show that SMB has decreased sharply in all regions since 2000. Finally we show correlation between runoff close to two calving glaciers and either calving front retreat or calving flux, this being most noticeable from the mid-1990s.

Ice-shelf basal melting is the largest contributor to the negative mass balance of the Antarctic ice sheet. However, current implementations of ice/ocean interactions in ice-sheet models disagree with the distribution of sub-shelf melt and freezing rates revealed by recent observational studies. Here we present a novel combination of a continental-scale ice flow model and a calibration technique to derive the spatial distribution of basal melting and freezing rates for the whole Antarctic ice-shelf system. The modelled ice-sheet equilibrium state is evaluated against topographic and velocity observations. Our high-resolution (10-km spacing) simulation predicts an equilibrium ice-shelf basal mass balance of −1648.7 Gt a−1 that increases to −1917.0 Gt a−1 when the observed ice-shelf thinning rates are taken into account. Our estimates reproduce the complexity of the basal mass balance of Antarctic ice shelves, providing a reference for parameterisations of sub-shelf ocean/ice interactions in continental ice-sheet models. We perform a sensitivity analysis to assess the effects of variations in the model set-up, showing that the retrieved estimates of basal melting and freezing rates are largely insensitive to changes in the internal model parameters, but respond strongly to a reduction of model resolution and the uncertainty in the input datasets.

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.