We evaluate the ability of the Community Earth System Model (CESM2) to simulate cryospheric–hydrological variables, such as glacier surface mass balance (SMB), over High Mountain Asia (HMA) by using a global grid (~111 km) with regional refinement (~7 km) over HMA. Evaluations of two different simulations show that climatological biases are reduced, and glacier SMB is improved (but still too negative) by modifying the snow and glacier model and using an updated glacier cover dataset.

We established a novel model that simulates the temporal changes in cryoconite hole (CH) depth using heat budgets calculated independently at the ice surface and CH bottom based on hole shape geometry. The simulations suggest that CH depth is governed by the balance between the intensity of the diffuse component of downward shortwave radiation and the wind speed. The meteorological conditions may be important factors contributing to the recent ice surface darkening via the redistribution of CHs.

We examine the thermal regime of the Grigoriev ice cap and the Sary-Tor glacier, both located in the inner Tien Shan in Kyrgyzstan. Our findings are important as the ice dynamics can only be understood and modelled precisely if ice temperature is considered correctly in ice flow models. The calibrated parameters of this study can be used in applications with ice flow models for individual ice masses as well as to optimise more general models for large-scale regional simulations.

We present a new approach for modelling debris area and thickness evolution. We implement the module into a combined mass-balance ice-flow model, and we apply it using different climate scenarios to project the future evolution of all glaciers in High Mountain Asia. We show that glacier geometry, volume, and flow velocity evolve differently when modelling explicitly debris cover compared to glacier evolution without the debris-cover module, demonstrating the importance of accounting for debris.

We simulate the current state and future evolution of the Mocho-Choshuenco ice cap in southern Chile (40°S, 72°W) with the ice-sheet model SICOPOLIS. Under different global warming scenarios, we project ice mass losses between 56 % and 97 % by the end of the 21st century. We quantify the uncertainties based on an ensemble of climate models and on the temperature dependence of the equilibrium line altitude. Our results suggest a considerable deglaciation in southern Chile in the next 80 years.

Debris-covered glaciers have a greater extent than their debris-free counterparts due to insulation from the debris cover. However, the transient response to climate change remains poorly understood. We use a numerical model that couples ice dynamics and debris transport and varies the climate signal. We find that debris cover delays the transient response, especially for the extent. However, adding cryokarst features near the terminus greatly enhances the response.

To better estimate the uncertainty in glacier mass change modeling during the 20th century we ran an established model with an ensemble of meteorological data sets. We find that the total ensemble uncertainty, especially in the early 20th century, when glaciological and meteorological observations at glacier locations were sparse, increases considerably compared to individual ensemble runs. This stems from regions with a lot of ice mass but few observations (e.g., Greenland periphery).

We show that plutonium is an effective tracer to identify ice originating from the early 1960s at the surface of a mountain glacier after a long time within the ice flow, giving unique information on the long-term former ice motion. Combined with ice flow modelling, the dating can be extended to the entire glacier, and we show that an airplane which crash-landed on the Gauligletscher in 1946 will likely soon be released from the ice close to the place where pieces have emerged in recent years.

We use a numerical flow model to simulate the behaviour of the Djankuat Glacier, a WGMS reference glacier situated in the North Caucasus (Republic of Kabardino-Balkaria, Russian Federation), in response to past, present and future climate conditions (1752–2100 CE). In particular, we adapt a more sophisticated and physically based debris model, which has not been previously applied in time-dependent numerical flow line models, to look at the impact of a debris cover on the glacier’s evolution.

Direct records of glacier changes rarely go back more than the last 100 years and are few and far between. We used a sophisticated glacier model to simulate glacier length changes over the last 1000 years for those glaciers that we do have long-term records of, to determine whether the model can run in a stable, realistic way over a long timescale, reproducing recent observed trends. We find that post-industrial changes are larger than other changes in this time period driven by recent warming.

We use simple, physics-based models to compare how marine-terminating glaciers respond to changes at their marine margin vs. inland surface melt. Initial glacier retreat is more rapid for ocean changes than for inland changes, but in both cases, glaciers will continue responding for millennia. We analyze several implications of these differing pathways of change. In particular, natural ocean variability must be better understood to correctly identify the anthropogenic role in glacier retreat.

We introduce a novel approach for simulating glacier mass balances using a deep artificial neural network (i.e. deep learning) from climate and topographical data. This has been added as a component of a new open-source parameterized glacier evolution model. Deep learning is found to outperform linear machine learning methods, mainly due to its nonlinearity. Potential applications range from regional mass balance reconstructions from observations to simulations for past and future climates.

To provide estimates of past glacier mass changes, an adequate initial state is required. However, information about past glacier states at regional or global scales is largely incomplete. Our study presents a new way to initialize the Open Global Glacier Model from past climate information and present-day geometries. We show that even with perfectly known but incomplete boundary conditions, the problem of model initialization leads to nonunique solutions, and we propose an ensemble approach.

We investigate thickness change of Bhutanese glaciers during 2004–2011 using repeat GPS surveys and satellite-based observations. The thinning rate of Lugge Glacier (LG) is > 3 times that of Thorthormi Glacier (TG). Numerical simulations of ice dynamics and surface mass balance (SMB) demonstrate that the rapid thinning of LG is driven by both negative SMB and dynamic thinning, while the thinning of TG is minimised by a longitudinally compressive flow regime.

We have implemented a frontal ablation parameterization into the Open Global Glacier Model and have shown that inversion methods based on mass conservation systematically underestimate the mass turnover (and therefore the thickness) of tidewater glaciers when neglecting frontal ablation. This underestimation can rise up to 19 % on a regional scale. Not accounting for frontal ablation will have an impact on the estimate of the glaciers’ potential contribution to sea level rise.

Submarine melting (SM) has been discussed as potentially triggering the recently observed retreat at outlet glaciers in Greenland. How much it may contribute in terms of future sea level rise (SLR) has not been quantified yet. When accounting for SM in our experiments, SLR contribution of 12 outlet glaciers increases by over 3-fold until the year 2100 under RCP8.5. Scaling up from 12 to all of Greenland's outlet glaciers increases future SLR contribution of Greenland by 50 %.

Iceberg calving is a major factor in the retreat of outlet glaciers of the Greenland Ice Sheet. Massive block overturning calving events occur at major outlet glaciers. A major calving event in 2009 was triggered by the release of a smaller block of ice from above the waterline. Using a numerical model, we investigate the feasibility of this mechanism to drive large calving events. We find that relatively small perturbations induce forces large enough to open cracks in ice at the glacier bed.

We present global glacier volume projections for the end of this century, under a range of high-end climate change scenarios, defined as exceeding 2 °C global average warming. The ice loss contribution to sea level rise for all glaciers excluding those on the peripheral of the Antarctic ice sheet is 215.2 ± 21.3 mm. Such large ice losses will have consequences for sea level rise and for water supply in glacier-fed river systems.