In the described study, we derive an uncertainty range for global mean sea level rise (SLR) contribution from Thwaites Glacier in a 200-year period under an extreme ocean warming scenario. We derive the spatial and vertical resolutions needed for bedrock data acquisition missions in order to limit global mean SLR contribution from Thwaites Glacier to ±2 cm in a 200-year period. We conduct sensitivity experiments in order to present the locations of critical regions in need of accurate mapping.

The presence of diatoms in Antarctic ice cores has been scarcely documented and poorly understood. Here we present a detailed analysis of the spatial and temporal distribution of the diatom record preserved in a set of Antarctic ice cores. Our results reveal that the timing and amount of diatoms deposited present a strong geographical division. This study highlights the potential of the diatom record preserved in Antarctic ice cores to provide useful information about past environmental changes.

Monitoring the snow water equivalent (SWE) in high mountain regions is highly important and a challenge. We explore the use of muon counts to infer SWE temporally continuously. We deployed muonic cosmic ray snow gauges (µ-CRSG) on a Swiss glacier over the winter 2020/21. Evaluated with manual SWE measurements and SWE estimates inferred from neutron counts, we conclude that the µ-CRSG is a highly promising method for remote high mountain regions with several advantages over other current methods.

We used satellite images to create a comprehensive record of annual glacier change in northwest Greenland from 1972 through 2021. We found that nearly all glaciers in our study area have retreated and glacier retreat accelerated from around 1996. Comparing these results with climate data, we found that glacier retreat is most sensitive to water runoff and moderately sensitive to ocean temperatures. These can affect glacier fronts in several ways, so no process clearly dominates glacier retreat.

Convective heat transfer (CHT) is important in affecting thermal regimes in permafrost regions. We quantified its thermal impacts by contrasting the simulation results from three scenarios in which the Simultaneous Heat and Water model includes full, partial, and no consideration of CHT. The results show the CHT commonly happens in shallow and middle soil depths during thawing periods and has greater impacts in spring than summer. The CHT has both heating and cooling effects on the active layer.

Recent research indicates the importance of lateral transport of dissolved carbon in the polygonal tundra, suggesting that the freeze-up period could further promote lateral carbon transport. We conducted subsurface tracer simulations on high-, flat-, and low-centered polygons to test the importance of the freeze–thaw cycle and freeze-up time for tracer mobility. Our findings illustrate the impact of hydraulic and thermal gradients on tracer mobility, as well as of the freeze-up time.

This study investigates the past snow reliability (1961–2019) of 16 ski resorts in the French Alps using state-of-the-art snowpack modelling. We used snowmaking investment figures to infer the evolution of snowmaking coverage at the individual ski resort level. Snowmaking improved snow reliability for the core of the winter season for the highest-elevation ski resorts. However it did not counterbalance the decreasing trend in snow cover reliability for lower-elevation ski resorts and in spring.

We modelled the response of the Larsen C Ice Shelf (LCIS) and its tributary glaciers to the calving of the A68 iceberg and validated our results with observations. We found that the impact was limited, confirming that mostly passive ice was calved. Through further calving experiments we quantified the total buttressing provided by the LCIS and found that over 80 % of the buttressing capacity is generated in the first 5 km of the ice shelf downstream of the grounding line.

This companion paper deals with numerical particularities of partial differential equations underlying 1D snow models. In this first part we neglect mechanical settling and demonstrate that the nonlinear coupling between diffusive transport (heat and vapor), phase changes and ice mass conservation contains a wave instability that may be relevant for weak layer formation. Numerical requirements are discussed in view of the underlying homogenization scheme.

Ice motion can be used to better understand how weather and climate change affect the ice. Antarctic sea ice extent has shown large variability over the observed period, and dynamical features may also have changed. Our method allows for the quantification of rotational motion caused by wind and how this may have changed with time. Cyclonic motion dominates the Atlantic sector, particularly from 2015 onwards, while anticyclonic motion has remained comparatively small and unchanged.

Proglacial lakes were pervasive along the retreating continental ice margins after the Last Glacial Maximum. Similarly to the marine ice boundary, interactions at the ice-lake interface impact ice sheet dynamics and mass balance. Previous numerical ice sheet modeling studies did not include a dynamical lake boundary. We describe the implementation of an adaptive lake boundary condition in PISM and apply the model to the glacial retreat of the Laurentide Ice Sheet.

Modeling gas transport in ice sheets from surface to close-off is key to interpreting climate archives. Estimates of the diffusion coefficient and permeability of snow and firn are required but remain a large source of uncertainty. We present a new dataset of diffusion coefficients and permeability from 20 to 120 m depth at two Antarctic sites. We suggest predictive formulas to estimate both properties over the entire 100–850 kg m³ density range, i.e., anywhere within the ice sheet column.

Sea ice thickness anomalies during the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) winter in January, February and March 2020 were simulated with the coupled Regional Arctic climate System Model (RASM) and compared with CryoSat-2/SMOS satellite data. Hindcast and ensemble simulations indicate that the sea ice anomalies are driven by nonlinear interactions between ice growth processes and wind-driven sea-ice transports, with dynamics playing a dominant role.

We analyze state-of-the-art climate models’ ability to describe snow mass and whether biases in modeled temperature or precipitation can explain the discrepancies in snow mass. In winter, biases in precipitation are the main factor affecting snow mass, while in spring, biases in temperature becomes more important, which is an expected result. However, temperature or precipitation cannot explain all snow mass discrepancies. Other factors, such as models’ structural errors, are also significant.

Climate changes require a dynamic description of glaciers in hydrological models. In this study we focus on the local modelling of snow and firn. We tested our model at the site of Colle Gnifetti, 4400–4550 m a.s.l. The model shows that wind erodes all the precipitation of the cold months, while snow is in part conserved between April and September since higher temperatures protect snow from erosion. We also compared modelled and observed firn density, obtaining a satisfying agreement.

Thawing permafrost releases carbon to the atmosphere, enhancing global warming. Part of the permafrost soils have been flooded by rising sea levels since the last ice age, becoming subsea permafrost (SSPF). The SSPF is less studied than the part on land. In this study we use a global model to obtain rates of thawing of SSPF under different future climate scenarios until the year 3000. After the year 2100 the scenarios strongly diverge, closely connected to the eventual disappearance of sea ice.

In this study, we improve the regional climate model RACMO2 and investigate the climate of Antarctica. We have implemented a new radiative transfer and snow albedo scheme and do several sensitivity experiments. When fully tuned, the results compare well with observations and snow temperature profiles improve. Moreover, small changes in the albedo and the investigated processes can lead to a strong overestimation of melt, locally leading to runoff and a reduced surface mass balance.

In recent years, seven glaciers are confirmed in the northern Japanese Alps. However, their mass balance has not been clarified. In this study, we calculated the seasonal and continuous annual mass balance of these glaciers during 2015–2019 by the geodetic method using aerial images and SfM–MVS technology. Our results showed that the mass balance of these glaciers was different from other glaciers in the world. The characteristics of Japanese glaciers provide new insights for earth science.

A record minimum July sea ice extent, since 1979, was observed in 2020. Our results reveal that an anomalously high advection of energy and water vapor prevailed during spring (April to June) 2020 over regions with noticeable sea ice retreat. The large-scale atmospheric circulation and cyclones act in concert to trigger the exceptionally warm and moist flow. The convergence of the transport changed the atmospheric characteristics and the surface energy budget, thus causing a severe sea ice melt.

We describe, apply, and validate the Environment and Climate Change Canada automated sea ice tracking system (ECCC-ASITS) that routinely generates large-scale sea ice motion (SIM) over the pan-Arctic domain using synthetic aperture radar (SAR) images. The ECCC-ASITS was applied to the incoming image streams of Sentinel-1AB and the RADARSAT Constellation Mission from March 2020 to October 2021 using a total of 135 471 SAR images and generated new SIM datasets (i.e., 7 d 25 km and 3 d 6.25 km).

We develop a regional linear Markov model consisting of four modules with seasonally dependent variables in the Pacific sector. The model retains skill for detrended sea ice extent predictions for up to 7-month lead times in the Bering Sea and the Sea of Okhotsk. The prediction skill, as measured by the percentage of grid points with significant correlations (PGS), increased by 75 % in the Bering Sea and 16 % in the Sea of Okhotsk relative to the earlier pan-Arctic model.