The ice-cored Galena Creek Rock Glacier, Wyoming, USA, has been the subject of a number of studies that sought to determine the origin of its ice. We present new observations of the rock glacier's internal structure from ground-penetrating radar to constrain ice and debris distribution and accumulation. We imaged dipping reflectors in the center of the glacier that are weak and discontinuous, in contrast to strong reflectors toward the edge of the cirque beneath large debris-avalanche chutes. These reflectors form a network of concave-up, up-glacier dipping layers. We interpret these as englacial debris bands formed by large debris falls buried by subsequent ice and snow accumulation. They are discontinuous where ice outpaces debris accumulation, but with sufficient debris accumulation an interleaved pattern of ice and debris layers can form. We propose a model in which the ice in these interleaved layers is snowfall preserved by debris-facilitated accumulation. Large debris falls that occur in early spring bury sections of the snowpack, which are then preserved through summer and incorporated into the rock glacier body over time. This study highlights the importance of sequential accumulation of ice and debris for understanding the dynamics of rock glaciers and debris-covered glaciers.

Antarctic lakes with perennial ice covers provide the opportunity to investigate in-lake processes without direct atmospheric interaction, and to study their ice-cover sensitivity to climate conditions. In this study, a numerical model – driven by radiative, atmospheric and turbulent heat fluxes from the water body beneath the ice cover – was implemented to investigate the impact of climate change on the ice covers from two Antarctic lakes: west lobe of Lake Bonney (WLB) and Crooked Lake. Model results agreed well with measured ice thicknesses of both lakes (WLB – RMSE= 0.11 m over 16 years of data; Crooked Lake – RMSE= 0.07 m over 1 year of data), and had acceptable results with measured ablation data at WLB (RMSE= 0.28 m over 6 years). The differences between measured and modeled ablation occurred because the model does not consider interannual variability of the ice optical properties and seasonal changes of the lake's thermal structure. Results indicate that projected summer air temperatures will increase the ice-cover annual melting in WLB by 2050, but that the ice cover will remain perennial through the end of this century. Contrarily, at Crooked Lake the ice cover becomes ephemeral most likely due to the increase in air temperatures.

We present a 3-year record of seasonal variations in ice speed and frontal ablation of 10 marine-terminating outlet glaciers along the coast of Prudhoe Land in northwestern Greenland. The glaciers showed seasonal speedup initiated between late May and early June, and terminated between late June and early July. Ice speed subsequently decreased from July to September. The timing of the speedup coincided with the onset of the air temperature rise to above freezing, suggesting an influence of meltwater availability on the glacier dynamics. No clear relationship was found between the speedup and the terminus position or the sea-ice/ice-mélange conditions. These results suggest that the meltwater input to the glacier bed triggered the summer speedup. The excess of summer speed (June–August) over the mean for the rest of the year accounted for 0.5–13% of the annual ice motion. Several glaciers showed seasonal frontal variations, i.e. retreat in summer and advance in winter. This was not due to ice-speed variations, but was driven by seasonal variations in frontal ablation. The results demonstrate the dominant effect of glacier surface melting on the seasonal speedup, and the importance of seasonal speed patterns on longer-term ice motion of marine-terminating outlet glaciers in Greenland.

The high-Alpine ice-core drilling site Colle Gnifetti (CG), Monte Rosa, Swiss/Italian Alps, provides climate records over the last millennium and beyond. However, the full exploitation of the oldest part of the existing ice cores requires complementary knowledge of the intricate glacio-meteorological settings, including glacier dynamics. Here, we present new ice-flow modeling studies of CG, focused on characterizing the flow at two neighboring drill sites in the eastern part of the glacier. The3-D full Stokes ice-flow model is thermo-mechanically coupled and includes firn rheology, firn densification and enthalpy transport, and is implemented using the finite element software Elmer/Ice. Measurements of surface velocities, accumulation, borehole inclination, density and englacial temperatures are used to validate the model output. We calculate backward trajectories and map the catchment areas. This constrains, for the first time at this site, the so-called upstream effects for the stable water isotope time series of the two ice cores drilled in 2005 and 2013. The model also provides a 3-D age field of the glacier and independent ice-core chronologies for five ice-core sites. Model results are a valuable addition to the existing glaciological and ice-core datasets. This especially concerns the quantitative estimate of upstream conditions affecting the interpretation of the deep ice-core layers.

Mount Achernar moraine is a terrestrial sediment archive that preserves a record of ice-sheet dynamics and climate over multiple glacial cycles. Similar records exist in other blue ice moraines elsewhere on the continent, but an understanding of how these moraines form is limited. We propose a model to explain the formation of extensive, coherent blue ice moraine sequences based on the integration of ground-penetrating radar (GPR) data with ice velocity and surface exposure ages. GPR transects (100 and 25 MHz) both perpendicular and parallel to moraine ridges at Mount Achernar reveal an internal structure defined by alternating relatively clean ice and steeply dipping debris bands extending to depth, and where visible, to the underlying bedrock surface. Sediment is carried to the surface from depth along these debris bands, and sublimates out of the ice, accumulating over time (>300 ka). The internal pattern of dipping reflectors, combined with increasing surface exposure ages, suggest sequential exposure of the sediment where ice and debris accretes laterally to form the moraine. Subsurface structure varies across the moraine and can be linked to changes in basal entrainment conditions. We speculate that higher concentrations of debris may have been entrained in the ice during colder glacial periods or entrained more proximal to the moraine sequence.

The retention of meltwater in the accumulation area of the Greenland ice sheet and other Arctic ice masses buffers their contribution to sea level change. However, sustained warming also results in impermeable ice layers or ‘ice slabs’ that seal the underlying pore space. Here, we use a 1-D, physically based, high-resolution model to simulate the surface mass balance (SMB), percolation, refreezing, ice layer formation and runoff from across the high-elevation area of Devon Ice Cap, Canada, from 2001 to 2016. We vary the thickness of the ‘impermeable’ ice layer at which underlying firn becomes inaccessible to meltwater. Thick near-surface ice layers are established by an initial deep percolation, the formation of decimetre ice layers and the infilling of interleaving pore space. The cumulative SMB increases by 48% by varying impermeable layer thickness between 0.01 and 5 m. Within this range we identify narrower range (0.25–1 m) that can simulate both the temporal variability in SMB and the observed near-surface density structure. Across this range, cumulative SMB variation is limited to 6% and 45–49% of mass retention takes place within the annually replenished snowpack. Our results indicate cooler summers after intense mid-2000s warming have led to a partial replenishment of pore space.

We developed a tilt sensor for studying ice deformation and installed our tilt sensor systems in two boreholes drilled close to the shear margin of Jarvis Glacier, Alaska to obtain kinematic measurements of streaming ice. We used the collected tilt data to calculate borehole deformation by tracking the orientation of the sensors over time. The sensors' tilts generally trended down-glacier, with an element of cross-glacier flow in the borehole closer to the shear margin. We also evaluated our results against flow dynamic parameters derived from Glen's exponential flow law and explored the parameter space of the stress exponent n and enhancement factor E. Comparison with values from ice deformation experiments shows that the ice on Jarvis is characterized by higher n values than that is expected in regions of low stress, particularly at the shear margin (~3.4). The higher n values could be attributed to the observed high total strains coupled with potential dynamic recrystallization, causing anisotropic development and consequently sped up ice flow. Jarvis' n values place the creep regime of the ice between basal slip and dislocation creep. Tuning E towards a theoretical upper limit of 10 for anisotropic ice with single-maximum fabric reduces the n values by 0.2.

The response of snowpack to a +2°C global warming relative to the present climate was estimated in Hokkaido, Japan, using a physical snowpack model driven by dynamically downscaled (DDS) data, after model evaluation. The evaluation revealed that the snowpack model successfully reproduced the height of snow cover (HS), snow water equivalent (SWE) and snow-covered days (SCDs), but had a moderate bias in the thickness ratios of melt form (MF) and hoar category (HC). The DDS-forced simulation predicted that the seasonal-maximum HS and SWE would decrease by 30–40% in the southwestern and eastern parts of Hokkaido due to a large decrease in snowfall during the accumulation period, and that the HS and SWE in the north would decrease, albeit not significantly due to uncertain atmospheric forcing. The number of SCDs in Hokkaido was predicted to decline by ~30 d. Additionally, ~50% of snowpack thickness during a season would be MF in most areas, whereas HC would be <50% all over Hokkaido.

Like most ice caps and glaciers worldwide, Icelandic glaciers are retreating in a warming climate. Here, the evolution of Vatnajökull ice cap, Iceland, from 1980 to 2300 is simulated by forcing the Parallel Ice Sheet Model (PISM) with output from Regional Climate Models (RCMs). For climate simulations of the recent past, HARMONIE-AROME reanalysis-forced simulations are used, while for future climate conditions, high-resolution (5.5 km) simulations from the RCM HIRHAM5 are used in addition to available CORDEX simulations (12 km). The glacier evolution is modelled using the RCP 4.5 and RCP 8.5 scenarios until 2100. To extend the time series, the 2081–2100 climate forcing is repeated until 2300. For RCP 4.5, the ice cap loses 31–64% of its volume and 13–37% of its area by 2300 depending on the used model forcing. For RCP 8.5, the volume decrease is 51–94% and the area decrease is 24–80% by 2300. In addition, the effect of elevation feedbacks is investigated by adding a precipitation and temperature lapse rate to the HIRHAM5 simulations. By 2300, the lapse rate runs have a 9–14% smaller volume and a 9–20% smaller area than the runs without a lapse rate correction.

Retreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.

Climate change has led to a significant shrinkage of glaciers in the Tropical Andes during the last decades. Recent multi-temporal quantifications of ice mass loss at mountain range to regional scale are missing. However, this is fundamental information for future water resource planning and glacier change projections. In this study, we measure temporally consistent glacier area changes and geodetic mass balances throughout the Bolivian Cordillera Real and Tres Cruces based on multi-sensor remote-sensing data. By analyzing multi-spectral satellite images and interferometric SAR data, a glacier recession of 81 ± 18 km2 (29%; 5.1 ± 1.1 km2 a−1), a geodetic mass balance of −403 ± 78 kg m−2 a−1 and a total ice mass loss of 1.8 ± 0.5 Gt is derived for 2000–2016. In the period 2013–2016, ice mass loss was 21% above the average rate. A retreat rate of 15 ± 5 km2 a−1 and a mass budget of −487 ± 349 kg m−2 a−1 are found in this more recent period. These higher change rates can be attributed to the strong El Niño event in 2015/16. The analyses of individual glacier changes and topographic variables confirmed the dependency of the mass budget and glacier recession on glacier aspect and median elevation.

Accurate estimations of ice thickness and volume are indispensable for ice flow modelling, hydrological forecasts and sea-level rise projections. We present a new ice thickness estimation model based on a mass-conserving forward model and a Bayesian inversion scheme. The forward model calculates flux in an elevation-band flow-line model, and translates this into ice thickness and surface ice speed using a shallow ice formulation. Both ice thickness and speed are then extrapolated to the map plane. The model assimilates observations of ice thickness and speed using a Bayesian scheme implemented with a Markov chain Monte Carlo method, which calculates estimates of ice thickness and their error. We illustrate the model's capabilities by applying it to a mountain glacier, validate the model using 733 glaciers from four regions with ice thickness measurements, and demonstrate that the model can be used for large-scale studies by fitting it to over 30 000 glaciers from five regions. The results show that the model performs best when a few thickness observations are available; that the proposed scheme by which parameter-knowledge from a set of glaciers is transferred to others works but has room for improvements; and that the inferred regional ice volumes are consistent with recent estimates.

Passive seismology is becoming increasingly popular for glacier/ice-sheet structure investigations in Polar regions. Single-station passive seismic methods including P-wave receiver functions (PRFs), horizontal-to-vertical spectral ratio (HVSR) and a recently proposed autocorrelation method have been used to retrieve glacier/ice-sheet structures. Despite their successful applications, analysis regarding their detection abilities in different glaciological environments has not been reported. In this study, we compare ice thicknesses and vp/vs ratios obtained from the three methods using data collected at GAMSEIS and POLENET/ANET seismic arrays in Antarctica. Ice thickness estimates derived from the three methods are found to be consistent. Comparisons conducted under various model setups, including those involving tiled layers and sedimentary layers, show that the effectiveness of the autocorrelation method is not superior to the PRF method for retrieving ice-sheet structures. The autocorrelation method however can complement other methods as it only requires a single component seismic record.

Ground-penetrating radar data acquired in the 2016/17 austral summer on Sørsdal Glacier, East Antarctica, provide evidence for meltwater lenses within porous surface ice that are conceptually similar to firn aquifers observed on the Greenland Ice Sheet and the Arctic and Alpine glaciers. These englacial water bodies are associated with a dry relict surface basin and consistent with perennial drainage into an interconnected englacial drainage system, which may explain a large englacial outburst flood observed in satellite imagery in the early 2016/17 melt season. Our observations indicate the rarely-documented presence of an englacial hydrological system in Antarctica, with implications for the storage and routing of surface meltwater. Future work should ascertain the spatial prevalence of such systems around the Antarctic coastline, and identify the degree of surface runoff redistribution and storage in the near surface, to quantify their impact on surface mass balance.