Glacier phototroph blooms on the surfaces of ice sheets and glaciers cause albedo reduction, leading to increased melting rates. We observed seasonal changes in the abundance of phototrophs on the Qaanaaq Ice Cap in northwestern Greenland from June to August 2014, and reproduced these changes using numerical and empirical models. The phototroph community on the ice surface mainly consisted of the glacier alga Ancylonema nordenskioldii and the cyanobacterium Phormidesmis priestleyi. The glacier alga appeared on the ice surface in late June, after which its abundance increased exponentially throughout the melting period. A logistic growth model designed for snow algal growth reproduced the measured exponential increases, suggesting that growth could be explained using the model as a function of the ice melting duration. Cyanobacteria appeared and their abundance increased in late July but did not change exponentially thereafter. The abundance of cyanobacteria was explained with an empirical model expressed as a function of the amount of mineral dust on the bare ice surface. Our numerical and empirical models for reproducing glacier algae and cyanobacteria could be useful for quantifying the albedo reduction caused by their growth and the melt rates of the Greenland ice sheet and glaciers in the future.

Airborne radio-echo sounding (RES) surveys are widely used to measure ice-sheet bed topography. Measuring bed topography as accurately and widely as possible is of critical importance to modelling ice dynamics and hence to constraining better future ice response to climate change. Measurement accuracy of RES surveys is influenced both by the geometry of bed topography and the survey design. Here we develop a novel approach for simulating RES surveys over glaciated terrain, to quantify the sensitivity of derived bed elevation to topographic geometry. Furthermore, we investigate how measurement errors influence the quantification of glacial valley geometry. We find a negative bias across RES measurements, where off-nadir return measurement error is typically −1.8 ± 11.6 m. Topographic highlands are under-measured an order of magnitude more than lowlands. Consequently, valley depth and cross-sectional area are largely under-estimated. While overall estimates of ice thickness are likely too high, we find large glacier valley cross-sectional area to be under-estimated by −2.8 ± 18.1%. Therefore, estimates of ice flux through large outlet glaciers are likely too low when this effect is not taken into account. Additionally, bed mismeasurements potentially impact our appreciation of outlet-glacier stability.

Current sea-level rise partly stems from increased surface melting and meltwater runoff from the Greenland ice sheet. Multi-year snow, also known as firn, covers about 80% of the ice sheet and retains part of the surface meltwater. Since the firn cold content integrates its physical and thermal characteristics, it is a valuable tool for determining the meltwater-retention potential of firn. We use gap-filled climatological data from nine automatic weather stations in the ice-sheet accumulation area to drive a surface-energy-budget and firn model, validated against firn density and temperature observations, over the 1998–2017 period. Our results show a stable top 20 m firn cold content (CC20) at most sites. Only at the lower-elevation Dye-2 site did CC20 decrease, by 24% in 2012, before recovering to its original value by 2017. Heat conduction towards the surface is the main process feeding CC20 at all nine sites, while CC20 reduction occurs through low-cold-content fresh-snow addition at the surface during snowfall and latent-heat release when meltwater refreezes. Our simulations suggest that firn densification, while reducing pore space for meltwater retention, increases the firn cold content, enhances near-surface meltwater refreezing and potentially sets favourable conditions for ice-slab formation.

We constrain the history of the Greenland Ice Sheet margin during the Holocene at Upernavik Isstrøm, a major ice stream in northwestern Greenland. Radiocarbon-dated sediment sequences from proglacial-threshold lakes adjacent to the present ice margin constrain deglaciation of the sites to older than 9.6 ± 0.1 ka. This age of deglaciation is confirmed with 10Be ages of 9.9 ± 0.1 ka from an island adjacent to the historical ice position. The lake sediment sequences also constrain the ice margin to have been less extensive than it is today for the remainder of the Holocene until ~ 1100 to ~ 700 yr ago, when it advanced into two lake catchments. The ice margin retreated back out of these lake catchments in the last decade. The early Holocene deglaciation in Melville Bugt, one of few locations around Greenland where a vast stretch of the current ice margin is marine-based, preceded deglaciation in most other parts of Greenland. Earlier deglaciation in this ice-sheet sector may have been caused by additional ablation mechanisms that apply to marine-based ice margins. Furthermore, despite ice-sheet models depicting this sector of Greenland as relatively stable throughout the Holocene, our data indicate a > 20 km advance-retreat cycle within the last millennium.

The full history of ice sheet and climate interactions is recorded in the vertical profiles of geochemical tracers in polar ice sheets. Numerical simulations of these archives promise great advances both in the interpretation of these reconstructions and the validation of the models themselves. However, fundamental mathematical shortcomings of existing models subject tracers to spurious diffusion, thwarting straightforward solutions. Here, I propose a new vertical discretization for ice-sheet models that eliminates numerical diffusion entirely. Vertical motion through the model mesh is avoided by mimicking the real-world flow of ice as a thinning of underlying layers. A new layer is added to the surface at equidistant time intervals, isochronally, thus identifying each layer uniquely by its time of deposition and age. This new approach is implemented for a two-dimensional section through the summit of the Greenland ice sheet. The ability to directly compare simulations of vertical ice cores with reconstructed data is used to find optimal model parameters from a large ensemble of simulations. It is shown that because this tuning method uses information from all times included in the ice core, it constrains ice-sheet sensitivity more robustly than a realistic reproduction of the modern ice-sheet surface.

Surface meltwater can refreeze within firn layers and crevasses to warm ice through latent-heat transfer on decadal to millennial timescales. Earlier work posited that the consequent softening of the ice might accelerate ice flow, potentially increasing ice-sheet mass loss. Here, we calculate the effect of meltwater refreezing on ice temperature and softness in the Pâkitsoq (near Swiss Camp) and Jakobshavn Isbræ regions of western Greenland using a numeric model and existing borehole measurements. We show that in the Jakobshavn catchment, meltwater percolation within the firn warms the ice at depth by 3–5°C. By contrast, meltwater refreezing in crevasses (cryo-hydrologic warming) at depths of ~300 m warms the ice in Pâkitsoq by up to 10°C, but this causes minimal increase in ice motion (<10 m a−1). Pâkitsoq is representative of western Greenland's land-terminating ice, where the slow movement of ice through a wide ablation zone provides ideal conditions for cryo-hydrologic warming to occur. We find that only ~37% of the western Greenland ice flux, however, travels through such areas. Overall, our findings suggest that cryo-hydrologic warming will likely have only a limited effect on the dynamic evolution of the Greenland ice sheet.

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.