We improved the DOE's E3SM land model (ELMv1-ECA) simulations of soil temperature, zero-curtain period durations, cold-season CH₄, and CO₂ emissions at several Alaskan Arctic tundra sites. We demonstrated that simulated CH₄ emissions during zero-curtain periods accounted for more than 50 % of total emissions throughout the entire cold season (Sep to May). We also found that cold-season CO₂ emissions largely offset warm-season net uptake currently and showed increasing trends from 1950 to 2017.

Surface meltwater lakes have been observed on several Antarctic ice shelves in field studies and optical images. Meltwater lakes can drain and refreeze, increasing the fragility of the ice shelves. The combination of synthetic aperture radar (SAR) backscatter and interferometric information (InSAR) can provide the cryosphere community with the possibility to continuously assess the dynamics of the meltwater lakes, potentially helping to facilitate the study of ice shelves in a changing climate.

This study shows an elevation difference between the radar interferometric measurements and the optical measurements from a coordinated campaign over the snow-covered deformed sea ice in the western Weddell Sea, Antarctica. The objective is to correct the penetration bias of microwaves and to generate a precise sea ice topographic map, including the snow depth on top. Excellent performance for sea ice topographic retrieval is achieved with the proposed model and the developed retrieval scheme.

We investigate the long-term (1953–2017) morphodynamic changes in rock glaciers in Kaunertal valley, Austria. Using a combination of historical aerial photographs and laser scanning data, we derive information on flow velocities and surface elevation changes. We observe a loss of volume and an acceleration from the late 1990s onwards. We explain this by changes in the meteorological forcing. Individual rock glaciers react to these changes to varying degrees.

Cold content is the energy required to attain an isothermal (0 °C) state and resulting in the snow surface melt. This study focuses on determining the multi-layer cold content (30 min time steps) relying on field measurements, snow temperature profile, and empirical formulation in four distinct forest sites of Montmorency Forest, eastern Canada. We present novel research where the effect of forest structure, local topography, and meteorological conditions on cold content variability is explored.

The paper investigates the performance of altimetric satellite instruments to detect river ice onset and melting dates and to retrieve ice thickness of the Ob River. This is a first attempt to use satellite altimetry for monitoring ice in the challenging conditions restrained by the object size. A novel approach permitted elaboration of the spatiotemporal ice thickness product for the 400 km river reach. The potential of the product for prediction of ice road operation was demonstrated.

We use surface and borehole measurements to investigate the generation and fate of basal meltwater in the ablation zone of western Greenland. The rate of basal meltwater generation at borehole study sites increases by up to 20 % over the winter period. Accommodation of all basal meltwater by expansion of isolated subglacial cavities is implausible. Other sinks for water do not likely balance basal meltwater generation, implying water evacuation through a connected drainage system in winter.

This companion paper deals with numerical particularities of partial differential equations underlying one-dimensional snow models. In this second part we include mechanical settling and develop a new hybrid (Eulerian–Lagrangian) method for solving the advection-dominated ice mass conservation on a moving mesh alongside Eulerian diffusion (heat and vapor) and phase changes. The scheme facilitates a modular and extendable solver strategy while retaining controls on numerical accuracy.

Understanding the long-term behaviour of ice sheets is essential to projecting future changes due to climate change. In this study, we use rocks deposited along the margin of the David Glacier, one of the largest glacier systems in the world, to reveal a rapid thinning event initiated over 7000 years ago and endured for ~ 2000 years. Using physical models, we show that subglacial topography and ocean heat are important drivers for change along this sector of the Antarctic Ice Sheet.

Sea ice is retreating with rapid warming in the Arctic. It will continue and approach the worst predicted pathway released by the IPCC. The irreversible tipping point might show around 2060 when the oldest ice will have completely disappeared. It has a huge impact on human production. Ordinary merchant ships will be able to pass the Northeast Passage and Northwest Passage by the midcentury, and the opening time will advance to the next 10 years for icebreakers with moderate ice strengthening.

Snow depth data are essential to monitor the impacts of climate change on sea ice volume variations and their impacts on the climate system. For that purpose, we present and assess the altimetric snow depth product, computed in both hemispheres from CryoSat-2 and SARAL satellite data. The use of these data instead of the common climatology reduces the sea ice thickness by about 30 cm over the 2013–2019 period. These data are also crucial to argue for the launch of the CRISTAL satellite mission.

This study investigated the aerodynamic roughness length (z₀) – an important parameter to determine the surface roughness – of crevassed tidewater glaciers on Svalbard using drone data. The results point out that the range of z₀ values across a crevassed glacier is large but in general significantly higher compared to non-crevassed glacier surfaces. The UAV approach proved to be an ideal tool to provide distributed z₀ estimates of crevassed glaciers which can be used to model turbulent fluxes.

We present O₂/N₂ and Ar/N₂ records from the Dome Fuji ice core through the bubbly ice, bubble–clathrate transition, and clathrate ice zones without gas-loss fractionation. The insolation signal is preserved through the clathrate formation. The relationship between Ar/Ν₂ and Ο₂/Ν₂ suggests that the fractionation for the bubble–clathrate transition is mass independent, while the bubble close-off process involves a combination of mass-independent and mass-dependent fractionation for O₂ and Ar.

We have shown through field experiments that the amount of wave energy dissipated in landfast ice, sea ice attached to land, is much larger than in broken ice. By comparing our measurements against predictions of contemporary wave–ice interaction models, we determined which models can explain our observations and which cannot. Our results will improve our understanding of how waves and ice interact and how we can model such interactions to better forecast waves and ice in the polar regions.

About 10 % of Himalayan glaciers flow directly into lakes. This study finds, using satellite imagery, that such glaciers show higher flow velocities than glaciers without ice–lake contact. In particular near the glacier tongue the impact of a lake on the glacier flow can be dramatic. The development of current and new meltwater bodies will influence the flow of an increasing number of Himalayan glaciers in the future, a scenario not currently considered in regional ice loss projections.

Using the reanalysis datasets and the Canadian Ice Island Drift, Deterioration and Detection database, a probabilistic model was developed to quantify ice island fracture probability under various atmospheric and oceanic conditions. The model identified water temperature as the most dominant variable behind ice island fracture events, while ocean currents played a minor role. The developed model offers a predictive capability and could be of particular interest to offshore and marine activities.

We propose a generalized form for the damage parameterization such that super-critical stresses can return to the yield with different final sub-critical stress states. In uniaxial compression simulations, the generalization improves the orientation of sea ice fractures and reduces the growth of numerical errors. Shear and convergence deformations however remain predominant along the fractures, contrary to observations, and this calls for modification of the post-fracture viscosity formulation.

Antarctica is shrinking, and part of the mass loss is caused by higher temperatures leading to more snowmelt. We use computer models to estimate the amount of melt, but this can be inaccurate – specifically in the areas with the most melt. This is because the model cannot account for small, darker areas like rocks or darker ice. Thus, we trained a computer using artificial intelligence and satellite images that showed these darker areas. The model computed an improved estimate of melt.

We developed a novel automated method for the detection and the quantification of ocean waves generated by glacier calving. This method was applied to data recorded with a terrestrial radar interferometer at Eqip Sermia, Greenland. Results show a high calving activity at the glacier front sector ending in deep water linked with more frequent meltwater plumes. This suggests that rising subglacial meltwater plumes strongly affect glacier calving in deep water, but weakly in shallow water.

Along the edges of the Greenland Ice Sheet surface melt lubricates the bed and causes large seasonal fluctuations in ice speeds during summer. Accurately understanding how these ice speed changes occur is difficult due to the inaccessibility of the glacier bed. We show that by using surface velocity maps with high temporal resolution and numerical modelling we can infer the basal conditions that control seasonal fluctuations in ice speed and gain insight into seasonal dynamics over large areas.

Estimating how much ice loss from Greenland and Antarctica will contribute to sea level rise is of critical societal importance. However, our analysis shows that recent efforts are not trustworthy because the models fail at reproducing contemporary ice melt. Here we present a roadmap towards making more credible estimates of ice sheet melt.

We did a systematic analysis of the location of inclusions in the EGRIP ice core, the first ice core from an ice stream. We combine this with crystal orientation and grain size data, enabling the first overview about the microstructure of this unique ice core. Micro-inclusions show a strong spatial variability and patterns (clusters or horizontal layers); roughly one-third is located at grain boundaries. More holistic approaches are needed to understand deformation processes in the ice better.

With the increasing melt of the Greenland Ice Sheet, which contributes to sea level rise, the surface of the ice darkens. The dark surfaces absorb more radiation and thus experience increased melt, resulting in the melt–albedo feedback. Using a simple surface melt model, we estimate that this positive feedback contributes to an additional 60 % ice loss in a high-warming scenario and additional 90 % ice loss for moderate warming. Albedo changes are important for Greenland’s future ice loss.

There is a widespread controversy on the existence of the elevation-dependent warming (EDW) phenomenon due to the limited observations in high mountains. This study provides new evidence of EDW from the Chinese Tian Shan based on a high-resolution (1 km, 6-hourly) air temperature dataset. The result reveals the significant EDW on a monthly scale. The warming rate of the minimum temperature in winter showed a significant elevation dependence (p < 0.01), especially above 3000 m.

Lakes form on the surface of the Antarctic Ice Sheet during the summer. These lakes can generate further melt, break up floating ice shelves and alter ice dynamics. Here, we describe a new automated method for mapping surface lakes and apply our technique to the Amery Ice Shelf between 2005 and 2020. Lake area is highly variable between years, driven by large-scale climate patterns. This technique will help us understand the role of Antarctic surface lakes in our warming world.

Ambient noise correlation is a broadly used method in seismology to monitor tiny changes in subsurface properties. Some environmental forcings may influence this method, including snow. During one winter season, we studied this snow effect on seismic velocity of the medium, recorded by a pair of seismic sensors. We detected and modeled a measurable effect during early snowfalls: the fresh new snow layer modifies rigidity and density of the medium, thus decreasing the recorded seismic velocity.