A comprehensive review of historical changes in snow metrics in HA regions is discussed in the following sections. Different observation data sources (in situ and remote sensing) have been used to identify and interpret historical changes in snow metrics due to changing climate. Historical changes in snow cover were examined using the snow metrics outlined in
Table 1, which include precipitation (P), snow cover extent (SCE), snow cover duration (SCD), snow depth (SD) and snow water equivalent (SWE).
Most previous studies that evaluated and discussed changes in snow cover examined them over the whole Northern Hemisphere or pan-Arctic region. In this review, only published studies focusing on the historical and projected changes of snow cover over the HA region have been considered and included. Studies conducted in specific HA sectors were first reviewed, followed by studies made at the local scale. The trends in each snow metric are presented in tables and maps according to their location. Studies conducted at the regional scale were distinguished by italic format and in different colours inside the tables from the studies conducted on specific HA sectors. The direction of trends is indicated by arrows (upward: strongly positive; oblique upward: slightly positive; downward: strongly negative; oblique downward: slightly negative; horizontal arrow: no trend), and trend magnitude by a red (decreasing) to blue (increasing) colour scale. For studies that only reported the direction of the trend without its magnitude, a single plain colour was used to indicate the sign of the trend (red: decreasing; blue: increasing). Moreover, as studies used different units for reporting the trend in the snow metrics, trends were reported in common units to allow comparison between different studies (P: mm/decade, SD: cm/decade, SWE: mm water equivalent/decade, SCD: days/decade and SCE: percent/decade).
2.1. Precipitation
In the HA, snowfall is the main type of precipitation, mostly occurring in the wintertime and remaining on the ground between early September and early June [
1,
2]. Snowfall thus has a critical influence on the amount and duration of the snow cover [
2,
3]. Several studies of precipitation over the Arctic have reported a rising trend in precipitation during the last four decades [
4,
6] (
Table 2). This increase in precipitation was linked with an increase in water vapor in the troposphere [
7].
Winter precipitation has been increasing at high latitudes over most of the HA region (
Table 2). On the long time (1930–2010), the Atlantic sector has experienced a significant increase in winter precipitation. For example, Callaghan et al. [
22] analysed the trends in seasonal precipitation totals between October and May, corresponding to the snowfall season, for long-term meteorological stations and drifting stations located north of 70° N, for the 1936–2009 period. Their results indicated that cold season precipitation increased between 1936 and 1980 by as much as 9.1 mm per decade and 10.1 mm per decade between 1980 and 2010 (
Table 2). However, there has been spatial variability in precipitation trends at the local scale. Heterogeneous and weakly positive precipitation trends were notably reported in Svalbard, in the Atlantic sector; while the trend was positive in the north of the archipelago (1.5 mm per decade), it was negative in the south (−0.5 mm per decade), for the period of 1961–2012 [
17,
23].
As for the Atlantic sector, the North European sector also experienced a significant positive trend of winter precipitation, with an increase of 6.2 mm per decade between 1936 and 1980, and 9.9 mm per decade between 1980 and 2010 (
Table 2). Over the West Siberian sector, negative trends were reported for the period 1936–1980, while for the more recent 1980–2009 period, this sector encountered the greatest increase in winter precipitation (17.1 mm per decade) of all the HA sectors. In the west-Siberian sector of the HA, Frey and Smith [
24] reported a reduction of winter precipitation of greater magnitude (−13 mm per decade) during 1958–1993 (
Table 2).
The East Siberian, Chukchi and Alaskan sectors experienced an overall decrease in precipitation of, respectively, −3.5, −8.4 and −0.4 mm per decade between 1930 and 1980, while precipitation increased in the other HA sectors (
Table 2). However, after 1980, Chukchi still had a negative, albeit smaller trend with (−3.3 per decade), while the Siberian and Alaskan sectors showed a positive trend of 1 and 10.1 mm per decade, respectively.
The Canadian sector experienced an increase in cold season precipitation in the long-term, i.e., over 1930–2010 (
Table 2). The greatest increase happened during the 1930–1980 period (8.5 mm per decade), while an increase with smaller magnitude occurred afterward (2.9 mm per decade during 1980–2010). Another study by Zhang et al. [
25] reported a similar trend but with smaller magnitude of 5.2 mm per decade for longer period (1961–1998). Some studies investigated precipitation trends at the regional and local scales across the Canadian HA (
Table 2). As for Callaghan et al. [
22], they found a positive trend in snowfall and a significant increase in annual (35%) and winter (45%) precipitation (indicating snowfall) over the Canadian HA over that period. In general, the ratio of snowfall to total precipitation has been rising in the Canadian HA, due mainly to an increase in winter precipitation that generally falls as snow (an average of 25%). Another study at the regional scale found a positive trend in cold season snowfall with a magnitude of 54% between 1948 and 2012 [
8]. The same trend albeit of slightly smaller magnitude was also reported at the local scale: Young et al. [
26] analysed precipitation trends in Resolute (Nunavut, Canada) for the period 1950–2009 and found that snowfall has increased by around 20%. A study by Bjorkman et al. [
27] reported an increase in total winter snowfall of 10.9 cm per decade at Eureka (Nunavut) over the period 1995–2010.
While snowfall has overall increased during the cold season in the HA, Screen and Simmonds [
28] revealed a pronounced decline in summer snowfall across the Arctic Ocean and Canadian Arctic Archipelago from 1989 to 2009, based on in situ and an atmospheric reanalysis data (ERA-Interim). They concluded that the snowfall decrease is almost entirely driven by changes in precipitation phase (snow turning to rain), which they attributed to the Arctic warming trend over that period.
2.2. Snow Cover Duration
Trends in snow cover duration (SCD) observed at different scales over the last century show that SCD has declined across most of the HA sectors. For example, Derksen et al. [
20] reported that annual snow cover duration is decreasing by up to 5 days per decade over the Arctic since 1967. Callaghan et al. [
22] reported even higher declines over the HA sectors, with the most extreme reduction noted in the Canadian and north European sectors of the HA, by as much as 10 days per decade over 37 years (1972–2009).
There has been evidence of regional variability in SCD reduction over the different HA sectors. For example, Derksen et al. [
29] revealed that during the period 1980–2015 snow cover onset occurred earlier than normal over the eastern Canadian HA and later than normal over the Eurasian and Alaskan HA sectors. Snow-off dates were near normal over most of the Eurasian HA over the period 1981–2010. Bulygina et al. [
30] selected different snow depth thresholds to derive SCD and quantify the impact of climate change on SCD over the Eurasian sector for the period 1966–2007. When using a 1 cm depth threshold for snow presence, they found that there was a significant reduction in SCD across all the HA sectors, by 2–6 days per decade. The main reason for the decline was attributed to decreasing snow depth over most of Eurasia. Zhang and Ma [
31] also studied the trends in snow cover onset and disappearance dates, and snow cover persistence across Eurasia over the period of 1976–2006. Their results showed that both annual and seasonal SCD have decreased across the entire Eurasian HA sector, which is consistent with the earlier findings by Bulygina et al. [
30] and Brown and Robinson [
32]. Brown and Dreksen [
29] used daily snow depth observations from Canada and concluded that SCD has reduced by 4.4 days per decade annually over the period 1955–2013. However, SCD declined by 1.8 days per decade from August to January and by 2.3 days from February to July. The northern part of Europe and Eurasia demonstrated the same trend but with lower values. However, the most northwestern part of the Siberian sector showed a positive trend of 2–4 days/decade.
Several studies have attributed the negative trends in SCD to be mostly caused by earlier snow disappearance (
Figure 2). Euskirchen et al. [
33] studied seasonal snow cover variability over the 70–75° N latitude band, i.e., the southern region of the HA, over two warming periods of the 20th century (1910–1940 and 1970–2000) and found a greater sensitivity of snow cover to warming temperature during the more recent period. They also showed that spring SCD declined faster (1.1 days per decade) than the fall SCD (0.7 days per decade). Brown et al. [
34] performed an extensive study on SCD based on ten different sources of data including SD observations, reanalysis, and satellite data.
Figure 2 shows the trend in fall (left) and spring (right) SCD (days per decade) over the 1972–2008 period obtained from the NOAA weekly SCE data set maintained at Rutgers University [
32]. Fall SCD mostly decreased over the Canadian, the east Eurasian, and Chukchi HA sectors between 2 to 10 days per decades, while it increased over the west and east Siberia between 2 to 6 days per decades (
Figure 2).
Unlike the fall SCD, spring SCD showed a reduction between 4 to 10 days per decade within the HA during the 1972–2009 period (
Figure 2). The pan-Arctic study of Callaghan et al. [
22] showed that fall SCD decreased over all the HA except the East Siberian sector, with the largest decline in SCD observed over the North European and Canadian HA sectors. Spring SCD, on the other hand, decreased over all HA sectors except over the East Siberian and Atlantic sectors. This was consistent with the pattern of positive surface air temperature anomalies in the spring, observed over all Arctic land areas except for eastern North America [
12]. However, in some regions (western Eurasia), declining SCD trends were attributed to a later winter onset [
13]. Another study by Derksen et al. [
20] confirmed this trend of declining SCD mostly driven by earlier spring melt over the whole HA. Using snow cover map produced from the NOAA daily Interactive Multisensor Snow and Ice Mapping System (IMS) for the period 1998–2010, they reported an earlier snow cover onset in the fall over much of the HA; however, the northern part of the North European sector and east part of the Atlantic sector showed an opposite trend (earlier onset by 10–20 days/decades).
At the local scale, Kankaanpää et al. [
14] monitored the spatiotemporal variability of snowmelt patterns within the Zackenberg valley, northeast Greenland, over 2006–2014 [
14]. They revealed that despite large interannual variations in snowmelt timing, the spatial pattern of snowmelt remained itself similar among sites. On Bylot Island in the Canadian HA sector, Bilodeau et al. [
16] reported a significant reduction of 2 days per decade in SCD during the last two decades (2000s and 2010s). In Svalbard, Van Pelt et al. [
17] found that the snow onset date in the fall increased by around 2 days per decade due to warming, whereas there was no significant trend in the snow disappearance date in response to spring/summer warming the period 1961–2012. Also, a recent study in archipelago (Svalbard) found a clear patterns of earlier snow melt date in the southern and central parts of the archipelago, while the northernmost parts exhibit little change or trend toward later snow melt date for the period 2000–2019 [
18]. Young et al. [
26] studied the variability of SCD in Bathurst and Cornwallis Island in the Canadian HA using a combination of in situ (point-scale) snow survey and snowmelt measurements, as well landscape-scale snow cover observations from remote sensing. They found no reduction in SCD at the point scale due to insufficient data, but they observed a shallower snowpack around ponds and plateau sites while they found deeper snow in incised valleys.
2.3. Snow Cover Extent
The Arctic land surface is covered by snow during a large part of the year, and approximately snow free during the short (July–August) summer, which makes the fall and spring transition seasons very important in terms of characterizing the temporal variability in SCE [
20,
22]. Over the last century, SCE has shown a significant decrease across the HA, with the sharpest decline observed since the early 1970s [
32,
36,
37].
Several studies have used the weekly NOAA snow chart climate data record [
38] to analyse changes in SCE over different time periods since 1966 [
20,
30,
34,
39]. Mudryk et al. [
39] used the NOAA snow chart climate data record (CDR, [
32]) to analyse spatial trends in snow cover over the Northern Hemisphere during the period 1981–2010. A large-scale decrease in SCE was found during spring over all HA sectors, with the most drastic reduction (20%) occurred over the northern European and Alaskan sectors. However, a significant increase in SCE was observed during winter and fall over the northern European and eastern Atlantic sectors. We recomputed linear trends in snow cover extent using the same dataset (the NOAA snow chart climate data record) but for the extended 1972–2021 period (
Figure 3). Despite some differences among these studies and our updated trend maps regarding the seasonality of SCE trends, there is strong agreement that (i) annual SCE has declined overall; (ii) there is no evidence of significant changes in winter and early spring (JFM) SCE in the HA region, while (iii) extensive negative trends are seen during late spring (AMJ). Studies by Estilow et al. [
38] and Thackeray et al. [
40] reported that SCE declines were largest during summer, particularly during June and July over Arctic region. The updated trend map in
Figure 3 shows negative trends in both spring (AMJ) and summer (JAS) in the HA, however significant trends are less extensive in summer than spring. No significant trends were found in the fall (ND).
In the Eurasian HA, Zhang and Ma [
31] found no changes in SCE in winter during the 1975–2005 period, while declining SCE was observed in other seasons, particularly in summer with maximum reductions of 1.8% per decade. In Svalbard, North Atlantic HA, Van Pelt et al. [
17] showed that end-of-summer SCE declined strongly (from 48 to 36%) for the period 1961–2012. Despite differences in trend magnitude due to the use of different datasets and time periods, all these studies consistently indicated reductions in spring and summer SCE over the 1981–2015 in HA regions except in the Atlantic sector (mainly Greenland), while no trends were reported in the fall across the Canadian and Atlantic sectors [
38].
To compare and visualize the difference in SCE trends among cold regions (northern hemisphere, Arctic and High Artic), we analysed the seasonal trends in SCE using the NOAA weekly snow cover data [
38] for each domain (
Figure 4).
Despite the difference in the magnitude of the SCE reduction, SCE has declined overall across all the regions. However, both the Northern Hemisphere (NH) and the Arctic have experienced a significant SCE reduction since 1970 (1.03 and 0.3 million square kilometres per decade, respectively) while the HA has encountered a comparatively smaller reduction in spring SCE (0.03 million square kilometre per decade), due to the small area of the region. When looking at relative trends, the SCE reduction in spring has been much slower in the HA (−1.7%/decade) than the whole Arctic (−7.5%/decade) and NH (−5.5%/decade). SCE in the HA has reduced faster in summer (−6%/decade) than spring (−1.7%/decade), as is the case for the whole Arctic, due to the longer lasting snow cover. Overall, trends are slower in the HA compared to the Arctic-wide and NH-wide regions, as expected for the colder climate (
Figure 4).
2.4. Snow Depth
Recent studies have indicated that (SD) has increased overall in the last few decades across the HA, except over the Canadian sector (
Table 3). Snow depth was reported to increase over the Atlantic sector, for example by Van pelt et al. [
17] who found a positive trend in mean snow depth in Svalbard. However, a study by Pedersen et al. [
41] for Zackenberg showed no changes in snow depth. The greatest increase across the HA sectors was reported in the North European and Siberian sectors (
Table 3). For instance, Bulygina et al. [
30] analysed data for 1966 to 2007 over all the Eurasian sectors and found that the winter-averaged SD has increased in the western part of Russia, with maximum increases in Northern West Siberia (4 cm per decade). A similar positive trend was reported by Park et al. [
42] across all the HA Siberian sectors (2012) with a magnitude of 2 cm per decade. Another study by Kitaev et al. [
43] showed a similar trend but with a smaller magnitude of 1 cm per decade over all Siberian sectors between 1940 and 2002. Analysing snow depth in February from the same period, Popova [
44] concluded that the most significant changes in snow depth were strongly connected to the strengthening of the zonal circulation at high latitudes that have continued since the beginning of the 1970s. The anomaly of the planetary wave pattern, induced by the declining sea ice in the Barents and Kara Seas, reinforced moisture transport and caused greater snow depth across western Siberia [
45]. This is in accordance with earlier modelling results by Park et al. [
17], who found a significant negative correlation between sea-ice extent in September and snow depths across north-eastern Siberia during autumn and winter, i.e., reduced sea-ice leading to higher snow depths.
Various studies for the Canadian HA reported significant decreases in snow depth (
Table 3). For instance, Olsen et al., [
12] found that the most pronounced reductions in snow cover and SD have occurred in the Canadian HA over the last century, while SD was increasing in a large portion of Eurasia. The analysis of winter and early spring snow depths from meteorological stations for the 1945–1995 period over the Canadian HA showed that the greatest decreases occurred in February and March. Also, maximum snow depth was reported to decrease by as much as 1 to 3 cm per decade in the Canadian HA during the 1960–2015 period [
46]. Using observations from climate stations in Canadian HA regions, Vincent et al. [
47] reported that the annual maximum snow depth decreased over the longer 1950–2012 period. These decreased SD were followed by significant decreases in spring and summer snow cover duration over most of the Arctic [
46]. Snow depths trends were not linear, but rather were marked by a quick reduction from the mid-1970s onward that overlapped with a widely recorded change in atmospheric circulation in the Pacific-North American sector [
48]. Park et al. [
42] reported that the Alaskan HA exhibited a decrease in maximum snow depth of −2 cm per decade during the 1960–2010 period.
Local scale studies also revealed negative trends in SD within the Canadian and Alaskan HA sectors. A study in Alert (Nunavut) by Smith et al. [
49] revealed that maximum winter snow depth obtained from Environment and Climate Change Canada meteorological stations, has decreased over 1950 to 2000 which is consistent with the larger-scale studies by Brown et al. [
34] and Mudryk et al. [
8]. Also, Stieglitz et al. [
50] found that the mean snow depth over Prudhoe Bay (Alaskan sector) decreased significantly from 1940 to 2000 with the most dramatic reduction occurring in the 1970s. However, Bilodeau et al. [
16] using in situ data from a meteorological station and manual measurements found that both the variability and the average winter snow depth have decreased over the last three decades in the western part of Bylot Island.
Table 3.
Historical trends in SD (cm per decade). Positive trends are delineated by the blue colour scale and negative trends by the red colour scale (upward arrow: strongly positive; oblique upward arrow: slightly positive; downward arrow: strongly negative; oblique downward arrow: slightly negative; horizontal arrow: no trend). References
2.5. Snow Water Equivalent
Compared to other snow metrics, monitoring snow water equivalent (SWE) is more complicated due to limited surface observations and lingering challenges in retrieving SWE from remote-sensing data [
51]. Despite this, recent studies have demonstrated significant reductions in SWE over much of the last century. The results from these studies are summarized in
Table 4. A combination of in situ, remote sensing data and climate model outputs have been used to study the trends and variability of SWE across the Arctic. Despite the uncertainties associated with SWE retrieval algorithms, satellite data are best suited to investigate trends in SWE at large spatial scales.
Muskett [
52] used multi-satellite sensors to derive the variations and trends in mean SWE at high latitudes of the Northern Hemisphere from 1978 to 2010. Their study revealed increasing SWE trends with regional variations in magnitude across the Siberian HA, while significant decreasing SWE trends were found across the western Siberia HA. Jeong et al. [
53] used the GlobSnow observation dataset to characterize the temporal changes in SWE in the spring (February–April) between 1980–2012 for the Northern Hemisphere. Based on their results, the trend of observed SWE for February, March, and April showed tremendous spatial variability. In the HA, the Chukchi sector was the only region with increasing SWE during the spring. This could be due to significant increases in the cumulative precipitation from the previous winter to spring months. Other sectors showed overall decreases over the 1980–2012 period, which could be largely attributable to rising temperatures for the same months and/or decreasing winter precipitation [
53].
Table 4.
Historical trends in SWE (mm per decade). Positive trends are delineated by the blue colour scale and negative trends by the red colour scale (upward arrow: strongly positive; oblique upward arrow: slightly positive; downward arrow: strongly negative; oblique downward arrow: slightly negative; horizontal arrow: no trend). References
The southern Canadian HA showed significant spatiotemporal variability in SWE trends [
54]. Decreasing trends were more pronounced in spring (10 mm/year) and in the southeast and south, while an increase (10 mm/year) was found in the southwest [
8]. With advancing in the spring season (March to April) the larger area of snow cover in the west demonstrated a negative trend of SWE. There were two main reasons for these trends: mostly because of increasing temperatures for the same months and/or decreasing cumulative precipitation from the December to spring months [
8,
53]. Li et al. [
55] used three different satellite products to estimate SWE trends at the scale of the Northern Hemisphere over a 32-year period (1979/80–2010/11). Their results showed a significant reduction in the total mass of snow for the HA regions. The most remarkable changes in total snow mass, −50 mm per year, occurred during January in northern Europe. Zhang and Ma [
31] investigated the trends in seasonal and annual snow water equivalent (SWE) over Eurasia from 1979 to 2005. They found that annual SWE has decreased all over the region, with North Europe and west-Siberia showing the maximum reduction (−5.8 mm/decade).
Trends in SWE have shown pronounced seasonality, with the most dramatic reduction (−15.7 mm/decade) occurring in the spring. Brown et al. [
10] analysed the trends in SWE using daily snow depth observations from the Global Historical Climatology Network (GHCN) and snow densities from snow surveys over the Arctic except Greenland. Their results indicated that SWE has decreased by −0.7 cm/decade across the Canadian HA during the period 1955–2013. It was also shown that SWE has decreased overall by −0.2 mm/decade over the Russian HA during that period. However, SWE increased during the earlier 1967–1993 period by 6.6 mm/decade, and then decreased afterwards by −4.9 mm/decade from 1994–2011 over the Russian Arctic. Mudryk et al. [
8] showed considerable spatial heterogeneity in maximum SWE (SWEmax) trends for 1981–2015 across the Canadian HA. They discovered that SWEmax has decreased by 5–7.5% per decade in the north and by 2–5% in the southern part, while the west part has increased by as much as 2–5% per decade. There were no changes in the eastern part of the HA. This is consistent with the results of Jeong et al. [
53] over the same period. A study on Svalbard snow cover changes published by Van Pelt et al. [
17] showed that SWEmax in winter/spring has increased slightly (1 cm/decade) over the period 1961–2012.