The volume of glaciers in Iceland (∼3,400 km3in 2019) corresponds to about 9 mm of potential global sea level rise. In this study, observations from 98.7% of glacier covered areas in Iceland (in 2019) are used to construct a record of mass change of Icelandic glaciers since the end of the 19th century i.e. the end of the Little Ice Age (LIA) in Iceland. Glaciological (in situ) mass-balance measurements have been conducted on Vatnajökull, Langjökull, and Hofsjökull since the glaciological years 1991/92, 1996/97, and 1987/88, respectively. Geodetic mass balance for multiple glaciers and many periods has been estimated from reconstructed surface maps, published maps, aerial photographs, declassified spy satellite images, modern satellite stereo imagery, and airborne lidar. To estimate the maximum glacier volume at the end of the LIA, a volume–area scaling method is used based on the observed area and volume from the three largest ice caps (over 90% of total ice mass) at 5–7 different times each, in total 19 points. The combined record shows a total mass change of −540 ± 130 Gt (−4.2 ± 1.0 Gt a−1on average) during the study period (1890/91 to 2018/19). This mass loss corresponds to 1.50 ± 0.36 mm sea level equivalent or 16 ± 4% of mass stored in Icelandic glaciers around 1890. Almost half of the total mass change occurred in 1994/95 to 2018/19, or −240 ± 20 Gt (−9.6 ± 0.8 Gt a−1on average), with most rapid loss in 1994/95 to 2009/10 (mass change rate −11.6 ± 0.8 Gt a−1). During the relatively warm period 1930/31–1949/50, mass loss rates were probably close to those observed since 1994, and in the colder period 1980/81–1993/94, the glaciers gained mass at a rate of 1.5 ± 1.0 Gt a−1. For other periods of this study, the glaciers were either close to equilibrium or experienced mild loss rates. For the periods of AR6 IPCC, the mass change rates are −3.1 ± 1.1 Gt a−1for 1900/01–1989/90, −4.3 ± 1.0 Gt a−1for 1970/71–2017/18, −8.3 ± 0.8 Gt a−1for 1992/93–2017/18, and −7.6 ± 0.8 Gt a−1for 2005/06–2017/18.

Area and volume changes and the average geodetic mass balance of the non-surging outlet glaciers of the southeast Vatnajökull ice cap, Iceland, during different time periods between ~ 1890 and 2010, are derived from a multi-temporal glacier inventory. A series of digital elevation models (DEMs) (~ 1890, 1904, 1936, 1945, 1989, 2002, 2010) are compiled from glacial geomorphological features, historical photographs, maps, aerial images, DGPS measurements and a lidar survey. Given the mapped basal topography, we estimate volume changes since the end of the Little Ice Age (LIA) ~ 1890. The variable volume loss of the outlets to similar climate forcing is related to their different hypsometry, basal topography, and the presence of proglacial lakes. In the post-LIA period, the glacierized area decreased by 164 km2 (or from 1014 to 851 km2) and the glaciers had lost 10–30 % of their ~ 1890 area by 2010 (anywhere from 3 to 36 km2). The glacier surface lowered by 150–270 m near the terminus and the outlet glaciers collectively lost 60 ± 8 km3 of ice, which is equivalent to 0.15 ± 0.02 mm of sea-level rise. The volume loss of individual glaciers was in the range of 15–50%, corresponding to a geodetic mass balance between −0.70 and −0.32 m w.e. a−1. The annual rate of mass change during the post-LIA period was most negative in 2002–2010, on average −1.34 ± 0.12 m w.e. a−1, which is among the most negative mass balance values recorded worldwide in the early 21st century.

Sub-meter resolution, stereoscopic satellite images allow for the generation of accurate and high-resolution digital elevation models (DEMs) over glaciers and ice caps. Here, repeated stereo images of Drangajökull ice cap (NW Iceland) from Pléiades and WorldView2 (WV2) are combined with in situ estimates of snow density and densification of firn and fresh snow to provide the first estimates of the glacier-wide geodetic winter mass balance obtained from satellite imagery. Statistics in snow- and ice-free areas reveal similar vertical relative accuracy (< 0.5 m) with and without ground control points (GCPs), demonstrating the capability for measuring seasonal snow accumulation. The calculated winter (14 October 2014 to 22 May 2015) mass balance of Drangajökull was 3.33 ± 0.23 m w.e. (meter water equivalent), with ∼ 60 % of the accumulation occurring by February, which is in good agreement with nearby ground observations. On average, the repeated DEMs yield 22 % less elevation change than the length of eight winter snow cores due to (1) the time difference between in situ and satellite observations, (2) firn densification and (3) elevation changes due to ice dynamics. The contributions of these three factors were of similar magnitude. This study demonstrates that seasonal geodetic mass balance can, in many areas, be estimated from sub-meter resolution satellite stereo images.

Icelandic glaciers have been losing mass since the Little Ice Age in the mid‐to‐late 1800s, with higher mass loss rates in the early 21st century, followed by a slowdown since 2011. As of yet, it remains unclear whether this mass loss slowdown will persist in the future. By reconstructing the contemporary (1958–2019) surface mass balance of Icelandic glaciers, we show that the post‐2011 mass loss slowdown coincides with the development of the Blue Blob, an area of regional cooling in the North Atlantic Ocean to the south of Greenland. This regional cooling signal mitigates atmospheric warming in Iceland since 2011, in turn decreasing glacier mass loss through reduced meltwater runoff. In a future high‐end warming scenario, North Atlantic cooling is projected to mitigate mass loss of Icelandic glaciers until the mid‐2050s. High mass loss rates resume thereafter as the regional cooling signal weakens.

A simulation of the surface climate of Vatnajökull ice cap, Iceland, carried out with the regional climate model HIRHAM5 for the period 1980–2014, is used to estimate the evolution of the glacier surface mass balance (SMB). This simulation uses a new snow albedo parameterization that allows albedo to exponentially decay with time and is surface temperature dependent. The albedo scheme utilizes a new background map of the ice albedo created from observed MODIS data. The simulation is evaluated against observed daily values of weather parameters from five automatic weather stations (AWSs) from the period 2001–2014, as well as in situ SMB measurements from the period 1995–2014. The model agrees well with observations at the AWS sites, albeit with a general underestimation of the net radiation. This is due to an underestimation of the incoming radiation and a general overestimation of the albedo. The average modelled albedo is overestimated in the ablation zone, which we attribute to an overestimation of the thickness of the snow layer and not taking the surface darkening from dirt and volcanic ash deposition during dust storms and volcanic eruptions into account. A comparison with the specific summer, winter, and net mass balance for the whole of Vatnajökull (1995–2014) shows a good overall fit during the summer, with a small mass balance underestimation of 0.04 m w.e. on average, whereas the winter mass balance is overestimated by on average 0.5 m w.e. due to too large precipitation at the highest areas of the ice cap. A simple correction of the accumulation at the highest points of the glacier reduces this to 0.15 m w.e. Here, we use HIRHAM5 to simulate the evolution of the SMB of Vatnajökull for the period 1981–2014 and show that the model provides a reasonable representation of the SMB for this period. However, a major source of uncertainty in the representation of the SMB is the representation of the albedo, and processes currently not accounted for in RCMs, such as dust storms, are an important source of uncertainty in estimates of snow melt rate.