The climate in Svalbard is undergoing amplified change compared to the global mean. This has major implications for runoff from glaciers and seasonal snow on land. We use a coupled energy balance–subsurface model, forced with downscaled regional climate model fields, and apply it to both glacier-covered and land areas in Svalbard. This generates a long-term (1957–2018) distributed dataset of climatic mass balance (CMB) for the glaciers, snow conditions, and runoff with a 1 km×1 km spatial and 3-hourly temporal resolution. Observational data including stake measurements, automatic weather station data, and subsurface data across Svalbard are used for model calibration and validation. We find a weakly positive mean net CMB (+0.09 m w.e. a⁻¹) over the simulation period, which only fractionally compensates for mass loss through calving. Pronounced warming and a small precipitation increase lead to a spatial-mean negative net CMB trend (−0.06 m w.e. a⁻¹ decade⁻¹), and an increase in the equilibrium line altitude (ELA) by 17 m decade⁻¹, with the largest changes in southern and central Svalbard. The retreating ELA in turn causes firn air volume to decrease by 4 % decade⁻¹, which in combination with winter warming induces a substantial reduction of refreezing in both glacier-covered and land areas (average −4 % decade⁻¹). A combination of increased melt and reduced refreezing causes glacier runoff (average 34.3 Gt a⁻¹) to double over the simulation period, while discharge from land (average 10.6 Gt a⁻¹) remains nearly unchanged. As a result, the relative contribution of land runoff to total runoff drops from 30 % to 20 % during 1957–2018. Seasonal snow on land and in glacier ablation zones is found to arrive later in autumn (+1.4 d decade⁻¹), while no significant changes occurred on the date of snow disappearance in spring–summer. Altogether, the output of the simulation provides an extensive dataset that may be of use in a wide range of applications ranging from runoff modelling to ecosystem studies.

Since the first estimates of Svalbard-wide glacier mass balance were made in the early 2000s, there has been great progress in remote sensing and modelling of mass balance, existing field records have been extended, field records at new locations have been added, and there has been considerable environmental change. There is a wide spread in the available estimates of both total mass balance and climatic mass balance, but there is overall agreement that the glaciers on Svalbard have been losing mass since the 1960s, with a tendency towards more negative mass balance since 2000. We define criteria to select data that are representative and of high credibility; this subset shows a more coherent evolution and reduced spread. In addition, we combine individual field mass balance records collected by different groups into a single dataset that samples glaciers across Svalbard and a range of different size classes. We find a close relationship between measured glacier mass balance and size of the glacier, in such a way that smaller glaciers experience more negative climatic mass balances. A qualitatively similar relationship between the accumulation area ratio and glacier area is found for all glaciers in the Svalbard, suggesting that the relation derived from glaciological records is not only an artefact caused by the limited number of samples (n=12). We apply this relation to upscale measured climatic mass balance for a new estimate for all glaciers of Svalbard. Our reconciled estimates are for the climatic mass balance -7 ± 2 Gt a-1 (2000-2019) and a total mass balance of -8 ± 6 Gt a-1. From the difference and the related uncertainty, we derive an estimate of frontal ablation amounting to ca. 5-10 Gt a-1. While this is similar to a previous estimate of Svalbard-wide frontal ablation, the uncertainties are large. Furthermore, several large and long-lasting surges have had considerable and multi-year impact on the total mass balance, and in particular on calving rates, emphasizing the need for better-resolved and more frequently updated estimates of frontal ablation.

Dust particles from high latitudes have a potentially large local, regional, and global significance to climate and the environment as short-lived climate forcers, air pollutants, and nutrient sources. Identifying the locations of local dust sources and their emission, transport, and deposition processes is important for understanding the multiple impacts of high-latitude dust (HLD) on the Earth's systems. Here, we identify, describe, and quantify the source intensity (SI) values, which show the potential of soil surfaces for dust emission scaled to values 0 to 1 concerning globally best productive sources, using the Global Sand and Dust Storms Source Base Map (G-SDS-SBM). This includes 64 HLD sources in our collection for the northern (Alaska, Canada, Denmark, Greenland, Iceland, Svalbard, Sweden, and Russia) and southern (Antarctica and Patagonia) high latitudes. Activity from most of these HLD sources shows seasonal character. It is estimated that high-latitude land areas with higher (SI ≥0.5), very high (SI ≥0.7), and the highest potential (SI ≥0.9) for dust emission cover >1 670 000 km², >560 000 km², and >240 000 km², respectively. In the Arctic HLD region (≥60 ^∘ N), land area with SI ≥0.5 is 5.5 % (1 035 059 km²), area with SI ≥0.7 is 2.3 % (440 804 km²), and area with SI ≥0.9 is 1.1 % (208 701 km²). Minimum SI values in the northern HLD region are about 3 orders of magnitude smaller, indicating that the dust sources of this region greatly depend on weather conditions. Our spatial dust source distribution analysis modeling results showed evidence supporting a northern HLD belt, defined as the area north of 50^∘ N, with a “transitional HLD-source area” extending at latitudes 50–58^∘ N in Eurasia and 50–55^∘ N in Canada and a “cold HLD-source area” including areas north of 60^∘ N in Eurasia and north of 58^∘ N in Canada, with currently “no dust source” area between the HLD and low-latitude dust (LLD) dust belt, except for British Columbia. Using the global atmospheric transport model SILAM, we estimated that 1.0 % of the global dust emission originated from the high-latitude regions. About 57 % of the dust deposition in snow- and ice-covered Arctic regions was from HLD sources. In the southern HLD region, soil surface conditions are favorable for dust emission during the whole year. Climate change can cause a decrease in the duration of snow cover, retreat of glaciers, and an increase in drought, heatwave intensity, and frequency, leading to the increasing frequency of topsoil conditions favorable for dust emission, which increases the probability of dust storms. Our study provides a step forward to improve the representation of HLD in models and to monitor, quantify, and assess the environmental and climate significance of HLD.

In recent years extreme winter warming events have been reported in arctic areas. These events are characterized as extraordinarily warm weather episodes, occasionally combined with intense rainfall, causing ecological disturbance and challenges for arctic societies and infrastructure. Ground-ice formation due to winter rain or melting prevents ungulates from grazing, leads to vegetation browning, and impacts soil temperatures. The authors analyze changes in frequency and intensity of winter warming events in the Nordic arctic region—northern Norway, Sweden, and Finland, including the arctic islands Svalbard and Jan Mayen. This study identifies events in the longest available records of daily temperature and precipitation, as well as in future climate scenarios, and performs analyses of long-term trends for climate indices aimed to capture these individual events. Results show high frequencies of warm weather events during the 1920s–30s and the past 15 years (2000–14), causing weak positive trends over the past 90 years (1924–2014). In contrast, strong positive trends in occurrence and intensity for all climate indices are found for the past 50 years with, for example, increased rates for number of melt days of up to 9.2 days decade−1 for the arctic islands and 3–7 days decade−1 for the arctic mainland. Regional projections for the twenty-first century indicate a significant enhancement of the frequency and intensity of winter warming events. For northern Scandinavia, the simulations indicate a doubling in the number of warming events, compared to 1985–2014, while the projected frequencies for the arctic islands are up to 3 times higher.

Compared to other Arctic ice masses, Svalbard glaciers are low-elevated with flat interior accumulation areas, resulting in a marked peak in their current hypsometry (area-elevation distribution) at ~450 m above sea level. Since summer melt consistently exceeds winter snowfall, these low-lying glaciers can only survive by refreezing a considerable fraction of surface melt and rain in the porous firn layer covering their accumulation zones. We use a high-resolution climate model to show that modest atmospheric warming in the mid-1980s forced the firn zone to retreat upward by ~100 m to coincide with the hypsometry peak. This led to a rapid areal reduction of firn cover available for refreezing, and strongly increased runoff from dark, bare ice areas, amplifying mass loss from all elevations. As the firn line fluctuates around the hypsometry peak in the current climate, Svalbard glaciers will continue to lose mass and show high sensitivity to temperature perturbations.