We use the apparent horizontal shift of an englacial tephra layer outcrop to calculate local glacier mass balance on Sléttjökull, a lobe of Mýrdalsjökull in Southern Iceland. For this approach, the dipping angle of the englacial tephra layer in the glacier upstream of the outcrop and the flow velocity of the ice need to be known. An earlier investigation was expanded by the application of ground-penetrating radar, detecting the depth of the tephra along tracks with a total length of 10 km. Interpolation between the tracks enables us to derive the dipping angle of the layer along several flow lines. Together with glacier surface velocities, determined from feature tracking, we are able to estimate the local surface mass balance from the horizontal displacement of the tephra outcrop using freely available satellite imagery without additional fieldwork. The earlier local balance series was extended to the period 2014/15 to 2019/20. Although the results for the individual profiles differ slightly from each other, they show the same temporal pattern and clear variations from year to year. The results are compared to traditional mass-balance data from Hofsjökull. The two series show a good agreement in their interannual variability.

Due to adverse snow and cloud conditions, only a few inventories are available for the maritime glaciers in New Zealand. These are difficult to compare as different approaches and baseline data have been used to create them. In consequence, glacier fluctuations in New Zealand over the past two decades are only known for a few glaciers based on field observations. Here we present the results of a new inventory for the ‘year 2000’ (some scenes are from 2001 and 2002) that is based on glacier outlines from a recently published inventory for the year 2016 and allowed consistent change assessment for nearly 3000 glaciers over this period. The year 2000 inventory was created by manual on-screen digitizing using Landsat ETM+ satellite imagery (15 m panchromatic band) in the background and the year 2016 outlines as a starting point. Major challenges faced were late and early seasonal snow, clouds and shadow, the geo-location mismatch between Landsat and Sentinel-2 as well as the correct interpretation of ice patches and ice under debris cover. In total, we re-mapped 2967 glaciers covering an area of 885.5 km2 in 2000, which is 91.7 km2 (or 10.4%) more than the 793.8 km2 mapped in 2016. Area change rates (mean rate −0.65% a−1) increase towards smaller glaciers. Strongest area loss from 2000 to 2016 occurred at elevations ~1900 m but the highest relative loss was found below 800 m a.s.l. In total, 109 glaciers split into two or more entities and 264 had wasted away by 2016.

At any one time 130 000 icebergs are afloat in the Southern Ocean; 97% of these are too small to be registered in current satellite-based databases, yet the melting of these small icebergs provides a major input to the Southern Ocean. We use a unique set of visual size observations of 53 000 icebergs in the South Atlantic Ocean, the SCAR International Iceberg Database, to derive average iceberg dissolution rates. Fracture into two parts is the dominant dissolution process for tabular icebergs, with an average half-life of 30 days for icebergs <4 km length and 60 days for larger icebergs. Complete shatter producing many icebergs <1 km length is rare. A side attrition rate of 0.23 m d−1 combined with drift speed of 6 km d−1, or any proportional change in both numbers fits the observed changes in iceberg distribution. The largest injection into the Southern Ocean of fresh water and any iceberg-transported material takes place in a ~2.3 × 10⁶ km2 zone extending east-northeast from the Antarctic Peninsula to the Greenwich meridian. The iceberg contribution to salinities and temperatures, with maximum contribution north of the Weddell Sea, differs in some regions, from those indicated by tracking large icebergs.

Over the last two decades, anomalous warming events have been observed in coastal Antarctic regions. While these events have been documented in the Ross Sea sector, the Antarctic interior is believed to have been buffered from warming. In this work, we present data from lakes located near Mt. Heekin and Thanksgiving Valley (~85° S) along the Shackleton Glacier, which are believed to be the southern-most Antarctic dry valley lakes. In 2018, the lakes were characterized, repeat satellite images were examined, and lake water chemistry was measured. Our analysis shows that lake areas recently increased, and the water-soluble ion chemistry indicates a flushing of salts from periglacial soils, likely from increased glacial melt as illustrated by water isotope data. Our results show that high southern latitude ice-free areas have likely been affected by warm pulses over the past 60 years and these pulses may be quasi-synchronous throughout the Transantarctic Mountains.

Using operational sea-ice maps, we provide first insight into the seasonal evolution of fast ice in the East Siberian Sea for the period between 1999 and 2021. The fast ice season tends to start later by 4.7 d per decade and to end earlier by 9.7 d per decade. As a result, there is a trend towards a shorter length of fast ice season by 2 weeks per decade. The analysis of air temperatures indicates that onset and end of the fast ice season are largely driven by thermodynamic processes. Two spatial modes (large, L-mode and small, S-mode) of East Siberian fast ice cover which have significant areal differences were distinguished. The occurrence of L- and S-modes was linked to the polarity of the Arctic Oscillation (AO) index. Negative AO phase leads to increased sea-ice convergence in the region, which in turn favours sea-ice grounding and promotes the development of large fast ice extent (L-mode). Lower deformation rates in the region during positive AO phase does not allow the formation of grounded features which results in small fast ice extent (S-mode). An analysis of sea-ice divergence confirms that L-mode seasons are characterised by higher on-shore convergence compared with S-mode seasons.

Arctic rain-on-snow (ROS) events can have significant impacts on Arctic wildlife and socio-economic systems. This study addresses the meteorology of two different Arctic ROS events. The first, occurring near Nuuk, Greenland, generated significant impacts, including slush avalanches. The second, less severe, event occurred within the community of Iqaluit, Nunavut, Canada. This research utilizes atmospheric reanalysis, automated surface observation station data and atmospheric soundings to determine the meteorological conditions driving these events and the differences between each case. In both cases, atmospheric blocking played a leading role in ROS initiation, with atmospheric rivers – narrow bands of high water vapor transport, typically originating from the tropics and subtropics – having both direct and indirect effects. Cyclone-induced low-level jets and resultant ‘warm noses’ of higher air temperatures and moisture transport were other key features in ROS generation. To our knowledge, our study is the first to visualize how the varying strength and manifestation of these coupled features contribute to differences in the severity of Arctic ROS events. The meteorological drivers identified here find support from other studies on Arctic ROS events and are similar to weather features associated with Arctic precipitation events of extreme magnitude.

Extending the record of glacier area changes into the past improves our understanding of climate change impacts. Although analogue maps showing historic glacier extents are abundant, digital outlines from before the satellite era are sparse as the digitisation of moraines and trimlines on freely available satellite images is challenging. With the now available very high-resolution images provided by Web Map Services (WMS), new doors are open for the precise digitisation. Here, we used the ESRI WMS to digitise Little Ice Age (LIA) glacier extents and present area changes since the LIA in four selected regions along with a detailed uncertainty analysis. We used modern glacier outlines as a starting point and additionally consulted Sentinel-2 images, the ArcticDEM and historic maps for interpretation. Dating records from the literature allowed calculating area change rates. In total, 493 LIA glaciers (4640 km2, now 891 ice bodies with 3590 km2) were digitised, yielding relative area changes of −20% (−0.14% a−1), −15% (−0.10% a−1), −26% (−0.16% a−1) and −61% (−0.19% a−1) for Alaska, Baffin Island, Novaya Zemlya and the tropics, respectively. The ESRI WMS images are a great asset to precisely map moraines and trimlines, but information about the timing of the related extents requires further sources.

Aldegondabreen is a relatively small (5.3 km2) land-terminating glacier located in Nordeskiöld Land of Svalbard, ~10 km southwest of Barentsburg. Cumulative mass balance during 2002–20 equalled −21.79 m w.e., which corresponds to 37% of the total mass loss. The annual mass balance (Ba) varied from −0.24 to −2.19 m w.e., while the winter mass balance (Bw) ranged between 0.36 and 0.85 m w.e. Ba and Bw were strongly correlated with the positive degree-day sums and solid precipitation amounts, respectively, measured at the Barentsburg weather station. There was also a strong correlation (r = 0.76) between Ba and Bw, which shows that winter accumulation amplifies the consecutive summer ice melt by changing the surface albedo. The trend in both observational series is not detectable because the period from 2005 to 2013 was characterized by relatively high Bw and not very negative Ba values. This was also observed on the other Svalbard glaciers, and was related to prevailing north-westerly flows over Svalbard during the summer. Therefore, the decadal periodicity of the Aldegondabreen mass balance follows general archipelago patterns that are determined by regional-scale factors. Thus, the surface mass-balance time series, which is now the longest one in the central part of the Spitsbergen Island, is representative for the archipelago.

The seasonal ice-free period in the Hudson Bay Complex (HBC) has grown longer in recent decades in response to warming, both from progressively earlier sea-ice retreat in summer and later sea-ice advance in fall. Such changes disrupt the HBC ecosystem and ice-based human activities. In this study, we compare 102 simulations from 37 models participating in phase 6 of the Coupled Model Intercomparison Project to the satellite passive microwave record and atmospheric reanalyses. We show that, throughout the HBC, models simulate an ice-free period that averages 30 d longer than in satellite observations. This occurs because seasonal sea-ice advance is unrealistically late and seasonal sea-ice retreat is unrealistically early. We find that much of the ice-season bias can be linked to a warm bias in the atmosphere that is associated with a southerly wind bias, especially in summer. Many models also exhibit an easterly wind bias during winter and spring, which reduces sea-ice convergence on the east side of Hudson Bay and impacts the spatial patterns of summer sea-ice retreat. These results suggest that, for many models, more realistic simulation of atmospheric circulation would improve their simulation of HBC sea ice.

This study describes the morphological and dynamic changes of Parkachik Glacier, Suru River valley, Ladakh Himalaya, India. We used medium-resolution satellite images; CORONA KH-4, Landsat and Sentinel-2A from 1971–2021, and field surveys between 2015 and 2021. In addition, we used the laminar flow-based Himalayan Glacier Thickness Mapper and provide results for recent margin fluctuations, surface ice velocity, ice thickness, and identified glacier-bed overdeepenings. The results revealed that overall the glacier retreated by −210.5 ± 80 m with an average rate of 4 ± 1 m a−1 between 1971 and 2021. Whereas a field study suggested that the glacier retreat increased to −123 ± 72 m at an average rate of −20 ± 12 m a−1 between 2015 and 2021. Surface ice velocity was estimated using COSI-Corr on the Landsat data. Surface ice velocity in the lower ablation zone was 45 ± 2 m a−1 in 1999–2000 and 32 ± 1 m a−1 in 2020–2021, thus reduced by 28%. Further, the maximum thickness of the glacier is estimated to be ~441 m in the accumulation zone, while for glacier tongue it is ~44 m. The simulation results suggest that if the glacier continues to retreat at a similar rate, three lakes of different dimensions may form in subglacial overdeepenings.

The Greenland Ice Sheet has become an increasingly larger contributor to sea level rise in the past two decades and is projected to continue to lose mass. Climate variability is expected to increase under future warming, but the effect of climate variability on the Greenland Ice Sheet volume is poorly understood and is adding to the uncertainty of the projected mass loss. Here we quantify the influence of inter-annual temperature variability on mass loss from the Greenland Ice Sheet using the PISM model. We construct an ensemble of temperature-forcing fields that accounts for inter-annual variability in temperature using reanalysis data from NOAA-CIRES over the period 1851–2014. We investigate the steady-state and transient response of the Greenland Ice Sheet. We find that the simulated steady-state ice-sheet volume decreases by 1.9 ± 0.4 cm of sea level equivalent when forced with a varying temperature forcing compared to a constant temperature forcing, and by 11.5 ± 1.4 cm when the variability is doubled. The northern basins are particularly sensitive with a change in volume of 0.9–1.1%. Our results emphasize the importance of including temperature variability in projections of future mass loss.

As warming intensifies across the Greenland ice sheet, an increasing number of shallow coring and radar studies are targeting the melt-impacted firn column to investigate meltwater processes. Highly inhomogeneous infiltration and refreezing, however, redistributes mass, distorting age–depth relationships and confounding comparisons between different cores. Here, we utilize a dynamic time warping algorithm for time series alignment to investigate and quantify the heterogeneous impact of melt processes on nine closely spaced (within 50 m) firn core-density profiles. The 10 m cores were collected relatively high in Greenland's percolation zone, where melt alteration is minimal compared to lower elevation. Our analysis demonstrates the effectiveness of dynamic time warping as a tool for assessing heterogeneity between ice core records. We find that the optimal alignment of density profiles in the nine cores requires vertical stretching and compression of individual profiles, ranging from, on average, <1 to ~16% of the core lengths. We identified four depth zones of mass redistribution that appear to coincide with observed ice layers. Further, ~75% of density measurements from each core do not align with an age model-derived density profile that assumes no mass redistribution of meltwater, indicating the pervasive impact of melt processes.

Recent satellite measurements of glacier mass balance show mountain glaciers all over the world had generally negative mass balances in the first decades of the 21st century. We analyse archived data for surface mass balance and summer temperature for 37 Northern Hemisphere glaciers with data for 1961–2020. We compare mean annual balances for 1961–90 and 1991–2020, and for 25 glaciers explain the changes in annual balance by changes in winter and summer balances. Mean balances 1961–90 were already substantially negative for 19 out of the 25 glaciers. Changes in winter balances from 1961–90 to 1991–2020 average close to zero but changes in summer balance are strongly negative. Mean balance 1991–2020 is strongly correlated with change in summer balance, weakly correlated with winter balance change, and strongly correlated with mean balance 1961–90. We estimate 1991–2020 summer temperature anomalies for the 37 glaciers and confirm that summer temperature anomalies for 1991–2020 were higher for the Alps, by nearly 1.5°C, than for other areas. Substantial variations in the temperature – sensitivity of summer balances for individual glaciers of −0.2 to −1.0 m w.e. a−1 °C−1 deserve further study.

The uncertainty of glacier change projections is largely influenced by glacier models. In this study, we focus on temperature-index mass-balance (MB) models and their calibration. Using the Open Global Glacier Model (OGGM), we examine the influence of different surface-type dependent degree-day factors, temporal climate resolutions (daily, monthly) and downscaling options (temperature lapse rates, temperature and precipitation corrections) for 88 glaciers with in-situ observations. Our findings indicate that higher spatial and temporal resolution observations enhance MB gradient representation due to an improved calibration. The addition of surface-type distinction in the model also improves MB gradients, but the lack of independent observations limits our ability to demonstrate the added value of increased model complexity. Some model choices have systematic effects, for example weaker temperature lapse rates result in smaller projected glaciers. However, we often find counter balancing effects, such as the sensitivity to different degree-day factors for snow, firn and ice, which depends on how the glacier accumulation area ratio changes in the future. Similarly, using daily versus monthly climate data can affect glaciers differently depending on the shifting balance between melt and solid precipitation thresholds. Our study highlights the importance of considering minor model design differences to predict future glacier volumes and runoff accurately.

The mass balance of glaciers requires more detailed and continuous observations to understand their seasonal change in relation to climate. Here, we designed and installed an automated real-time monitoring platform at 4645 m a.s.l. on the Baishui River Glacier No.1 to collect continuous high-resolution observational data, and analyzed the seasonal dynamic from glacier movement and surface mass balance from glacier melting and snow accumulation. The results showed that the platform moved northeastward ~12.9 m at a rate of 0.06 ± 0.02 m d−1 between September 2021 and April 2022. The surface mass balance showed a varied temporal period. July and August were the main ablation periods, while ablation decreased and ceased in September. The glacier neither melted nor accumulated much between October and December, but began to have rapid snow accumulation in January. The glacier surface temperature varied with the air temperature and showed significant inter-seasonal differences among monsoon, post-monsoon and winter seasons. The surface mass balance also exhibited a strong response to the air temperature changes, with an average decrease of 1°C the point mass balance increased by 0.11 m w.e. from monsoon to post-monsoon and 0.22 m w.e. from post-monsoon to winter. Moreover, we found snowfall caused a decrease in the glacier surface temperature by increasing the surface albedo.

The geodetic mass balance of a glacier corresponds to glacier-wide volume changes, converted to mass changes using density assumptions. It is typically calculated by differencing multi-temporal digital elevation models. In this study, we show how the annual geodetic mass balance of a glacier can be derived from uncrewed aerial vehicle (UAV) data. The presented workflow is applied to two small- to medium-sized glaciers in the Kyrgyz Tien Shan (Central Asia): Bordu glacier and Sary-Tor glacier. The obtained geodetic mass balance is compared with the glaciological mass balance derived from a network of ablation stakes and snow pits. A previously calibrated mass-balance model is used to correct for the difference in acquisition dates. The results show that the determined geodetic mass balance matches closely with the glaciological mass balance. Besides, for both glaciers the geodetic mass balance does not seem to be particularly sensitive to the assumptions regarding volume-to-mass conversion. Therefore, our results demonstrate that UAVs can serve as a valuable instrument to quantify the annual geodetic mass balance and to validate the glaciological mass balance. The conventional glaciological mass-balance estimation often relies on interpolation and extrapolation methods, whereas UAVs offer the potential for direct data acquisition over the entire glacier surface.

Glaciers of Baffin Island and nearby islands of Arctic Canada have experienced rapid mass losses over recent decades. However, projections of loss rates into the 21st century have so far been limited by the availability of model calibration and validation data. In this study, we model the surface mass balance of the largest ice cap on Baffin Island, Penny Ice Cap, since 1959, using an enhanced temperature index model calibrated with in situ data from 2006–2014. Subsequently, we project changes to 2099 based on the RCP4.5 climate scenario. Since the mid-1990s, the surface mass balance over Penny Ice Cap has become increasingly negative, particularly after 2005. Using volume–area scaling to account for glacier retreat, peak net mass loss is projected to occur between ~2040 and 2080, and the ice cap is expected to lose 22% (377.4 Gt or 60 m w.e.) of its 2014 ice mass by 2099, contributing 1.0 mm to sea level rise. Our 2015–2099 projections are approximately nine times more sensitive to changes in temperature than precipitation, with an absolute cumulative difference of 566 Gt (90 m w.e.) between +2 and −2°C scenarios, and 63 Gt (10 m w.e.) between +20% and −20% precipitation scenarios.

We analyze the internal structure of two polythermal glaciers, Hurd and Johnsons, located on Livingston Island, Antarctica, using 200 and 750 MHz GPR data collected in 2003/04, 2008/09 and 2016/17 field campaigns. Based on the different permittivities of snow and ice, we determined the thickness distribution of the end-of winter snow cover and of the cold ice layer. Their knowledge is fundamental for mass balance and glacier dynamics studies due to the different densities and rheological properties of such media. The average measured thicknesses for the snow and cold ice layers (the latter including the snow layer) were of 1.44 ± 0.09 and 29.1 ± 1.5 m, and their corresponding maxima were of 2.45 ± 0.21 and 80.8 ± 2.5 m. GPR snow profiling allowed for extension of the coverage of the snow thickness survey, but added little information to that supplied by snow pits, stake readings and manual snow probing, because of the multiplicity of reflections within the seasonal snowpack caused by internal ice layers and lenses. The polythermal structure determined for Hurd Glacier fits into the so-called Scandinavian type, seldom reported for the Antarctic region.

Melting and calving of glaciers and ice caps in Antarctica and Greenland could potentially contribute significantly to global sea level rise. Updates to existing outlines that provide critical glacier baseline information in both regions could help in the analysis of particular changes in glacier parameters such as area and volume from time-series inventories. Here we synthesize previously established techniques and apply new multi-source datasets to update glacier outlines in selected test areas of Antarctica and Greenland, as well as to reduce uncertainties and errors during the mapping process. The workflow includes mapping glacier boundaries, subdividing glaciers by watersheds and assigning glacier attributes. Complicated glacier scenarios and updating challenges in polar regions are discussed and demonstrated by representative case studies. For the first time in Antarctica, we analyze the effect of terminus types on mapped glacier areas, and in Greenland we compare the differences with glacier mapping results using Landsat OLI and ETM+. With new data sources, the methods described in this study might help to create glacier outlines on a larger scale in Antarctica and Greenland. Although data sources can be substituted, the enormous amount of manual labor required to update glacier inventories remains a significant challenge.