Ground-penetrating radar of supraglacial debris thickness is used to study local thickness variability. Freshly emergent debris cover appears to have higher skewness and kurtosis than more mature debris covers. Accounting for debris thickness variability in ablation models can result in markedly different ice ablation than is calculated using the mean debris thickness. Slope stability modelling reveals likely locations for locally thin debris with high ablation.

Calving is an important mechanism that controls the dynamics of Greenland outlet glaciers. We test and compare four calving laws and assess which calving law has better predictive abilities. Overall, the calving law based on von Mises stress is more satisfactory than other laws, but new parameterizations should be derived to better capture the detailed processes involved in calving.

Snow avalanches show complicated flow behaviour, characterized by several flow regimes which coexist in one avalanche. In this work, we analyse flow regime transitions where a powder snow avalanche transforms into a plug flow avalanche by incorporating warm snow due to entrainment. Prediction of such a transition is very important for hazard mitigation, as the efficiency of protection dams are strongly dependent on the flow regime, and our results should be incorporated into avalanche models.

Snowfall observations over Antarctica are scarce and currently limited to information from the CloudSat satellite. Here, a first evaluation of the CloudSat snowfall record is performed using observations of ground-based precipitation radars. Results indicate an accurate representation of the snowfall climatology over Antarctica, despite the low overpass frequency of the satellite, outperforming state-of-the-art model estimates. Individual snowfall events are however not well represented.

A high-resolution spatially lumped hydrologic surface routing model is proposed to simulate meltwater transport over bare ice surfaces. In an ice-covered catchment, meltwater is routed by slow interfluve flow (~10⁻³–10⁻⁴ m s⁻¹) followed by fast open-channel flow (~10⁻¹ m s⁻¹). Seasonal evolution of supraglacial stream-river networks substantially alters the magnitude and timing of moulin discharge with implications for subglacial hydrology and ice dynamics.

We derive the first continuous record of total ice discharged from all large Greenland outlet glaciers over the 2000–2016 period, resolving a distinct pattern of seasonal variability. We compare these results to glacier retreat and meltwater runoff and find that while runoff has a limited impact on ice discharge in summer, long-term changes in discharge are highly correlated to retreat. These results help to better understand Greenland outlet glacier sensitivity over a range of timescales.

Regional climate models (RCMs) used to estimate the surface mass balance (SMB) of Antarctica depend on boundary forcing fields including sea surface conditions. Here, we assess the sensitivity of the Antarctic SMB to perturbations in sea surface conditions with the RCM MAR using unchanged atmospheric conditions. Significant SMB anomalies are found for SSC perturbations in the range of CMIP5 global climate model biases.

Drifting snowstorms that are hundreds of meters in depth are reproduced using a large-eddy simulation model combined with a Lagrangian particle tracking method, which also exhibits obvious spatial structures following large-scale turbulent vortexes. The horizontal snow transport flux at high altitude, previously not observed, actually occupies a significant proportion of the total flux. Thus, previous models may largely underestimate the total mass flux and consequently snow sublimation.

When the porous compacted snow layers on an ice shelf extend below sea level, seawater is able to infiltrate onto the shelf. Here it can affect measurements of ice shelf thickness by changing the average density and affect iceberg calving if the seawater enters fractures. Seawater infiltration has only been directly observed in a few locations around Antarctica. Using continent-wide geometry and snow density data we show that it may be more widespread than previously realised.

Thwaites Glacier, West Antarctica, has experienced rapid grounding line retreat and mass loss in the past decades. In this study, we simulate the evolution of Thwaites Glacier over the next century using different model configurations. Overall, we estimate a 5 mm contribution to global sea level rise from Thwaites Glacier in the next 30 years. However, a 300 % uncertainty is found over the next 100 years, ranging from 14 to 42 mm, depending on the model setup.