Decadal Variability in the Impact of Atmospheric Circulation Patterns on the Winter Climate of Northern Russia

: The Arctic continues to warm at a much faster rate than the global average. One process contributing to ‘‘Arcticampliﬁcation’’involveschangesinlow-frequencymacroscaleatmosphericcirculationpatternsandtheirconsequentinﬂuenceonregionalclimate.Here,usingERA5data,weexaminedecadalchangesintheimpactofsevensuchpatternson winter near-surface temperature (SAT) and precipitation (PPN) in northern Russia and calculate the temporal consistency of any statistically signiﬁcant relationships. We demonstrate that the 40-yr climatology hides considerable decadal variability in the spatial extent of such circulation pattern–climate relationships across the region, with few areas where their temporal consistency exceeds 60%. This is primarily a response to the pronounced decadal expansion/contraction and/or mobility of the circulation patterns’ centers of action. The North Atlantic Oscillation (NAO) is the dominant pattern (having the highest temporal consistency) affecting SAT west of the Urals. Farther east, the Scandinavian (SCA), Polar/ Eurasian (POL), and West Paciﬁc patterns are successively the dominant pattern inﬂuencing SAT across the West Siberian Plains, Central Siberian Plateau, and mountains of Far East Siberia, respectively. From west to east, the SCA, POL, and Paciﬁc–North American patterns exert the most consistent decadal inﬂuence on PPN. The only temporally invariant sig- niﬁcant decadal relationships occur between the NAO and SAT and the SCA and PPN in small areas of the North European Plain.


Introduction
The climate of the northern high latitudes is changing faster than elsewhere. Recent warming in near-surface Arctic temperatures has been more than twice that at lower latitudes (e.g., Overland et al. 2016), a phenomenon termed Arctic amplification. Moreover, regional precipitation has also increased (e.g., Box et al. 2019). Both changes have been at least partially ascribed to anthropogenic forcing (e.g., Chylek et al. 2014;Min et al. 2008). A significant proportion of the terrestrial Arctic comprises northern Russia, which stretches more than 7000 km from ;308E eastward to ;1708W (Fig. 1). This region has the potential to exert a huge influence on future greenhouse gas emissions because of the substantial amounts of carbon stored in its wetlands, soil, boreal forest, and terrestrial and shelf seas permafrost (e.g., Groisman et al. 2017, and references therein;Vasiliev et al. 2020). Many of these carbon sinks can only exist as such within narrow temperature ranges and are especially vulnerable to the ongoing warming, which was particularly marked across much of Siberia during the first half of 2020 (World Weather Attribution 2020). Therefore, it is imperative to understand the multiple processes that contribute to Arctic climate change, not least because of the divergent conclusions of studies assessing the impact of Arctic amplification on severe winter weather in midlatitudes (e.g., Cohen et al. 2020, and references therein). One such process concerns changes in the atmospheric circulation at northern high latitudes, which alters the spatial distribution of wind-induced heat advection and precipitationbearing clouds with associated regional changes in the radiation budget (Ye et al. 2015). Much of the low-frequency, macroscale atmospheric variability can be expressed as a series of atmospheric circulation (or teleconnection) patterns, which describe anomalous circulation conditions that occur simultaneously over remote regions. The locations of the maximum anomalies, usually defined in a pressure or geopotential height field, are termed ''centers of action,'' with their sign reversing between the opposite phases (or polarities) of the pattern. The majority of these circulation patterns are associated with large-scale changes in the atmospheric wave and jet stream patterns, and the centers of action appear as ''standing'' waves or north-south-oriented dipoles (e.g., Handorf and Dethloff 2012). They are driven primarily by either internal atmospheric dynamics or atmosphere-ocean-ice interactions, both of which can be modified by external anthropogenic forcing. These circulation patterns are generally more prevalent and have centers of action of greater magnitude during the Denotes content that is immediately available upon publication as open access.
Supplemental information related to this paper is available at the Journals Online website: https://doi.org/10.1175/JCLI-D-20-0566.s1. boreal winter, and thus engender a larger impact on climate variability in this season.
While these patterns may be considered as a continuum rather than discrete regimes (Franzke and Feldstein 2005), they are more usually described as separate entities following their definition in the pioneering work of Wallace and Gutzler (1981) and Barnston and Livezey (1987). The vast expanse of northern Russia means that its climate is shaped by circulation patterns originating in several oceanic basins: the North Atlantic Oscillation (NAO), East Atlantic pattern (EA), East Atlantic/Western Russia pattern (EAWR), and Scandinavian pattern (SCA) in the North Atlantic, the West Pacific pattern (WP) and Pacific-North American pattern (PNA) in the North Pacific, and the Polar/Eurasian pattern (POL) over the Arctic basin. Summary details of these seven circulation patterns are provided in Table 1. There is a large body of scientific research examining the impact of the North Atlantic and North Pacific patterns on the climate of western Europe and North America, respectively: however, there are far fewer examples in the western literature that investigate their influence on Eurasia including northern Russia. We note that while an alternative atmospheric catalogue (the Vangengeim-Girs classification) is sometimes employed in Russia and eastern Europe-see Degirmend zić and Ko_ zuchowski (2019) for a description and recent use of this classificationthere are marked similarities between this and many of the atmospheric circulation patterns listed above (Hoy et al. 2013;Gao et al. 2016).

Circulation pattern Acronym
Boreal winter centers of action relevant for this study (signs for positive polarity) References North Atlantic Oscillation NAO North-south dipole: Greenland/Iceland (2) and the Azores (1) Barnston and Livezey (1987) Hurrell (1995) East Atlantic EA North-south dipole: across the North Atlantic, generally displaced southeast from NAO pattern Barnston and Livezey (1987) West Pacific WP North-south dipole: Kamchatka Peninsula (2) and Southeast Asia/ western subtropical North Pacific (1) Wallace and Gutzler (1981) Barnston and Livezey (1987) Pacific-North American PNA Standing wave: intermountain region of North America (1) and south of Aleutian Islands (2) Barnston and Livezey (1987) East Atlantic/Western Russia EAWR (Eurasia-2) Standing wave: western Europe (1), north of Caspian Sea (2), northern China (1) Barnston and Livezey (1987) Lim (2015) Scandinavian SCA (Eurasia-1) Standing wave: western Europe (2), Scandinavia (1), eastern Russia (2) Barnston and Livezey (1987) Bueh and Nakamura (2007) (1) A positive NAO generally favors warmer and wetter conditions over much of northern Russia (e.g., Hurrell 1995;Hurrell et al. 2003;Shmakin and Popova 2006;Bader et al. 2011). For example, a positive winter NAO is usually associated with southerly flow over the far northwest of the region (e.g., Kislov and Matveeva 2020), leading to a warmer and wetter Kola Peninsula (Marshall et al. 2016), and is a major impact on winter climate extremes in this area (Casanueva et al. 2014;Marshall et al. 2020). Popova (2007) documented a negative relationship between the NAO and snow depth in northwest Russia because in mild winters precipitation is more likely to fall as rain. Other studies have analyzed the effects of the Arctic Oscillation (AO), which is broadly analogous to the NAO, and shown similar climatic relationships (e.g., Frey and Smith 2003;Yu et al. 2017). Previous work has also revealed that variability in autumn Eurasian snow cover, Arctic sea ice extent, North Atlantic sea surface temperature, and regional atmospheric circulation anomalies can feed back onto the subsequent winter NAO/AO (e.g., Chen et al. 2020;Cohen et al. 2020;Henderson et al. 2018;Kryjov 2015;Wu et al. 2013). The EA and EAWR patterns affect the intensity of the winter westerly wind belt over Eurasia (Gao et al. 2016). The positive polarity of the EAWR is associated with warming over eastern Siberia and has a wavelike appearance consistent with its wave train structure (Liu et al. 2014). It has also been linked to above average precipitation in the West Siberian Plains and the Kola Peninsula (Ionita 2014;Lim 2015). During the positive phase of the winter SCA most of Eurasia north of 408N is colder than normal (Liu et al. 2014), with an increase in the frequency of extreme cold events (Yu et al. 2017). Decreased precipitation occurs over northwestern Russia, the West Siberian Plains, and some of the Arctic coastal region, associated with reduced storm-track activity (Bueh and Nakamura 2007;Popova 2007;Liu et al. 2014). Crasemann et al. (2017) established that a positive winter SCA was more likely to occur following reduced Arctic sea ice in the preceding autumn although Peings (2019) found that blocking in this region was associated with internal variability rather than regional sea ice anomalies. There have been especially few studies regarding the effects of the North Pacific patterns-WP and PNA-and POL on Russian climate. Popova (2007) demonstrated an increase (decrease) in snow depth in the West Siberian Plains (the Far East region of Siberia) related to both positive PNA and POL. Interestingly, she also noted significant negative snow depth anomalies in the Kola Peninsula linked to the WP, despite this pattern having a center of action thousands of kilometers away in the North Pacific.
The majority of these previous studies have examined the climatological (long-term average) influence of circulation patterns on northern Russian climate over several decades. However, it is known that in addition to the relative frequency of the positive and negative polarities of a pattern varying through time (e.g., Hanna et al. 2015), the locations of the centers of action of these patterns are not temporally invariant and therefore neither is the spatial distribution of their climatic impacts (e.g., Goosse and Holland 2005; Vicente-Serrano and López-Moreno 2008; Zhang et al. 2008). The strength of the relationship between the circulation patterns and their proposed drivers (or statistical predictors) within the climate system may also be nonstationary, perhaps in response to changes in the background state (e.g., Kolstad and Screen 2019). Popova (2018) described decadal-scale variability in the climatic impact of the patterns in the north of Eurasia linked to switches in the rate of global warming. Moreover, modeling studies suggest that as the Arctic warms further the links between the dominant circulation patterns and Arctic surface climate are likely to vary considerably (van der Linden et al. 2017). Thus, in order to improve estimates of uncertainty in future climate projections it is important to understand how the climatic impacts of these atmospheric circulation patterns vary at decadal time scales and longer during the recent past. In this study we provide the first detailed analysis of the decadal variability in the influence of seven circulation patterns on the winter climate of northern Russia: such decadal-scale variability may be hidden in correlations calculated over longer time periods. Specifically, we examine seasonal gridded nearsurface air temperature (SAT) and precipitation (PPN) reanalysis data over the past four decades. Changes in these two parameters can have significant consequences: for example, disrupting fragile natural Arctic terrestrial ecosystems (Cooper 2014) and causing devastating socioeconomic impacts for indigenous peoples dependent on tundra reindeer nomadism (Forbes et al. 2016).
In section 2 we describe the data and statistical methodologies employed. In presenting the results (section 3) we 1) define the average circulation pattern-climate relationships for the past 40 years; 2) as examples of decadal variability, examine how these vary across four successive decades in response to changing locations of the centers of action; 3) summarize the temporal consistency of statistically significant decadal relationships and determine which circulation patterns are most temporally consistent with SAT and/or PPN across different regions of northern Russia; and 4) extend the length of the analysis using an alternative gridded climate dataset. Finally, in section 5, we summarize our principal conclusions, and discuss them in the context of previous and proposed future work.

a. Climate data
We utilize monthly SAT and PPN data from ERA5, the latest European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis, which is currently available from 1979 to the present. Seasonal winter values were simply calculated as the mean of the December-February monthly values for 1980-2019 where the year refers to the January. ERA5 has both higher spatial and temporal resolution than its predecessors: 31 km and hourly, respectively. In combination with improved modeling of surface conditions and more consistent sea surface temperature and sea ice analyses, this helps give improved SAT in the Arctic (Hersbach et al. 2020). Moreover, precipitation over Arctic sea ice was found to be more representative of the limited observations than earlier reanalyses (C. Wang et al. 2019).
To extend the length of our analysis and test the robustness of the results we use the Japanese Meteorological Agency (JMA) 55-year Reanalysis (JRA-55) (Kobayashi et al. 2015), for which data are available from 1958 to the present. The spatial resolution of JRA-55 is approximately half that of ERA5, at ;55 km, and the monthly SAT (PPN) data are the mean of four daily 6-h analyses (forecasts). Importantly, Marshall et al. (2018) found that JRA-55 reproduced the climate of Arctic Fennoscandia with a similarly high accuracy across its whole time series.
Time series of the seven atmospheric circulation patternsthe NAO, EA, WP, PNA, EAWR, SCA, and POL-were obtained from the National Weather Service Climate Prediction Center (CPC), where the methodology for calculating them, as employed by Barnston and Livezey (1987), is also described. In short, rotated principal component analysis was used to identify the patterns based on monthly mean standardized 500-hPa height anomalies derived from the NCEP-NCAR reanalysis (Kalnay et al. 1996) in the region 208-908N for 1950-2000. An examination of the 12 calendar monthly sets of rotated modes of variability revealed 10 dominant circulation patterns, of which seven are analyzed here. Although some of these patterns only explain a small proportion of the total variance in the height fields, previous work has indicated they are all physically meaningful at seasonal frequencies (Barnston and Livezey 1987). The sign conventions of some of the patterns vary between different past studies: here we use those of the CPC. We note that the North Pacific Oscillation (NPO) may be considered to be the surface manifestation of the WP (Linkin and Nigam 2008) but here we use the CPC WP index as it is compatible with the time series of the other circulation patterns. We also note that the exact impact a circulation pattern has on regional climate may be dependent on how it is defined (e.g., Pokorná and Huth 2015). Of particular relevance to this study, Chen and Song (2018) demonstrated that the effect of the NPO (similar to the WP) on winter SAT in parts of Far East Siberia varied by several degrees Celsius across six definitions.
Seasonal winter indices were calculated in the same way as the seasonal reanalysis data. We do not include the Arctic Oscillation (Thompson and Wallace 1998) because it is highly correlated to the NAO (r 5 0.75 for the 40 winters examined here) and has a broadly similar impact on climate across northern Russia, although with slightly weaker (stronger) relationships in the west (east) of the region. To clarify, in this study we investigate how the climatic impacts of the atmospheric circulation patterns, rather than the patterns themselves, vary at decadal time scales.

b. Statistical methods
All regression and correlation analyses were undertaken on detrended time series, thus assuming no a priori link between the two datasets. Regression coefficients were calculated using standard least squares methodology and the statistical significance obtained using Student's t test. In the latter calculation, we assumed there is no correlation between successive winter values of the circulation patterns or meteorological data (e.g., Feldstein 2002). In addition, we employed the Brown-Forsythe test for the equality of variances (Brown and Forsythe 1974) and the false discovery rate test, as described in Wilks (2006).

a. Climatological winter relationships
To set the decadal variability of the influence of each of the different atmospheric circulation patterns on the winter climate of northern Russia in context, first we describe the mean climatological relationships for the 40 years from 1980 to 2019. In Fig. 2 we show the spatial pattern of the regression coefficients between each of the positive polarities of the circulation patterns and sea level pressure (SLP)-the contour spacing of which is analogous to the associated near-surface wind anomalies-and the regions where a significant relationship (p , 0.10) exists between the pattern and SAT and PPN. Note that circumpolar plots to 508N of the regression coefficients are provided in Fig. S1 in the online supplemental material, which enable that fraction of the spatial pattern over northern Russia to be set in the context of the full pattern. The relatively high spatial resolution of ERA5 allows us to detect the impact that major orographic barriers, such as the Urals, have on these relationships.

1) NORTH ATLANTIC OSCILLATION
The NAO center of action in the region is the northeast extension of the Icelandic low that stretches north of Fennoscandia (Fig. 2i). The resultant cyclonic circulation for positive NAO (NAO1) leads to anomalous southwesterly flow over northwestern Russia, bringing warmer, moist air from the south. The widespread impact of the NAO on SAT in northern Russia is clearly indicated in Fig. 2ii, with a significant positive relationship between the two across almost all land north of 608N and west of 1308E. The positive relationship between PPN and NAO across northwest Russia is less spatially coherent (Fig. 2iii): nevertheless, there is still a significant relationship between the NAO and PPN across the majority of land west of 908E. The rainshadow effect of the Urals at ;608E can be perceived by the lack of a significant relationship between the NAO and PPN on the lee (east) side.
2) EAST ATLANTIC EA1 has a center of action over the Arctic Ocean giving anticyclonic flow over much of northern Russia in combination with negative centers of action farther south (Fig. 2iv). The easterly circulation anomaly over the Arctic coast leads to significantly colder SAT in this region between 458-1208E, including Novaya Zemlya and Severnaya Zemlya (Fig. 2v). EA1 is also associated with significant decreases in PPN over parts of the same region together with smaller areas of increased PPN farther south, immediately east of the Urals-in contrast to NAO1 the locally anomalous easterly flow associated with EA1 means that the rain-shadow effect of the Urals is observed on the western side-and in parts of the Central Siberian Plateau (;1008E) ( Fig. 2vi).

3) WEST PACIFIC
The negative northern center of action of WP1 is located over the Bering Sea (Fig. 2vii), which gives weak northeasterly circulation anomalies in Far East Siberia. While WP1 is linked to significantly colder SAT over the mountains of eastern Siberia (Fig. 2viii), there are no large coherent areas of anomalous PPN in northern Russia associated with this circulation pattern (Fig. 2ix).

4) PACIFIC-NORTH AMERICA
Similar to WP1, PNA1 also has a negative center of action in the North Pacific. However, Fig. 2x reveals that the northerly wind anomalies are much greater than for WP1 (contours closer together). There is also a positive center of action over the Arctic Ocean, north of the Laptev Sea, which means that the anomalous circulation associated with PNA1 over northern Russia is confined to the far east of the region. In contrast to WP1, PNA1 has very little impact on regional SAT ( Fig. 2xi) but is concomitant with significantly reduced PPN over parts of Far East Siberia (Fig. 2xii).

5) EAST ATLANTIC/WESTERN RUSSIA
The winter EAWR1 pattern comprises four centers of action, one of which can be observed in Fig. 2xiii, with two of the others in Fig. S1v in the online supplemental material. This is a broad negative SLP anomaly, centered north of the Caspian Sea (;608E), which gives anomalous northerly flow over the Barents Sea and weaker southwesterly flow across much of northern Russia east of 908E. The associated winter SAT anomalies are warmer (colder) in parts of the Central Siberian Plateau (the East Siberian Sea) (Fig. 2xiv). The latter observation is noteworthy because one might expect the anomalously southerly flow associated with EAWR1 to give warmer SAT: it arises because the climatological winter temperature over the East Siberian Sea is warmer than the region of Far East Siberia to the south. The stronger northerly flow over the Barents Sea during EAWR1 gives greater PPN over the northern Kola Peninsula (Fig. 2xv). There are also areas of significant changes in PPN of both signs in the Central Siberian Plateau associated with EAWR1, with increases generally found farther south, matching the regions of anomalous warming.

6) SCANDINAVIAN
The SCA pattern has a primary center of action over the Kola Peninsula, which, as the CPC defines the SCA, is a positive SLP anomaly for SCA1 (Fig. 2xvi). Therefore, the resultant anticyclonic circulation causes anomalous northerly flow into northern Russia immediately to the east, over the Kara Sea. Given that the climatological SLP regression patterns of NAO1 and SCA1 are almost opposite in spatial extent (cf. Figs. 2i and 2xvi), it is not surprising that their associated climatic anomalies also generally ''mirror'' each other. Thus, there are significant negative SAT anomalies across most of northern Russia west of 1308E (Fig. 2xvii). Interestingly, in contrast to the NAO, the impact of the SCA on SAT does not extend as far north as the Arctic coast in parts of this region, such as the Kola, Yamal, and Taymyr Peninsulas. However, the region of significant decreases in PPN accompanying SCA1 is larger and more spatially coherent than the equivalent increases associated with NAO1 (cf. Figs. 2iii and 2xviii). 7) POLAR/EURASIA The principal regional center of action of the POL1 pattern comprises a broad negative anomaly centered over the Kara and Laptev Seas (Fig. 2xix), analogous to an intensified polar vortex in the Eastern Hemisphere. This circulation pattern leads to a more zonally directed flow at ;658N that inhibits meridional exchanges of heat and moisture between the Arctic and midlatitudes. Significant SAT anomalies linked to POL1 therefore form a north-south dipole, with colder (warmer) temperatures over the Arctic Ocean, including Zemlya Frantsa Iosifa (Central Siberian Plateau) (Fig. 2xx). Almost all the associated PPN anomalies are positive and although such areas are scattered across much of northern Russia, including the Barents Sea coast and parts of western and central Siberia, they are sufficiently large in extent to signify they are not simply noise (Fig. 2xxi).

b. Variability across successive decades
In this section we use four successive decades-1980-89, 1990-99, 2000-09, and 2010-19-as independent examples of the decadal variability of the impact of the seven atmospheric circulation patterns on the winter climate of northern Russia. Plots of the significant relationships between the positive polarities of the patterns and SAT and PPN are shown in Figs. 3 and 4, respectively, with equivalent plots of the regression coefficient provided in Figs. S2 and S3, respectively. while the seasonal regression between the patterns and SLP is shown in Fig. 5. The temporal variability of significant (p , 0.10) decadal pattern-climate relationships in different regions between the four decades (simply defined as when the majority of the region has such a relationship) is summarized in Fig. 6. Figures 3 and 4 reveal that there is indeed very marked decadal variability in the spatial locations and extent over which an individual circulation pattern exerts a statistically significant influence on the climate of northern Russia. This is despite there often being concomitant long-term ''climatological'' relationships over the same 40-yr period (cf. Fig. 2). Figure 6 demonstrates that there was not a single statistically significant ''regional-scale'' relationship between a circulation pattern and either SAT or PPN that occurred across all four decades. Some decadal differences are especially striking: for example, the change in spatial extent of the significant impact of the NAO between the 1980s and 1990s (cf. Figs. 3i and 4i and Figs. 3ii and 4ii). There are also some instances of reversals in the sign of significant regional-scale relationships, such as between the PNA and SAT in Far East Siberia in the 2000s and 2010s (cf. Figs. 3xv and 3xvi). Figures 3 and 4 indicate that the areas where there is a significant relationship at p , 0.10 but not at p , 0.05 are generally very small. Therefore, the results in Fig. 6 are not substantially dependent on using the lower significance level. However, there are one or two exceptions, such as the relationship between WP and SAT over the Northern European Plain during the 2000s (Fig. 4xi) and that between EAWR and PPN over the West Siberian Plains in the same decade (Fig. 5xix).
Examination of Fig. 5 suggests that a high proportion of the decadal variability is likely to be a direct response to the changing location and/or extent of a pattern's center(s) of action, which modifies the spatial distribution of advection of heat and moisture into and across the region. In the 1980s the NAO center of action north of Fennoscandia, the Icelandic low, expanded farther east than average, as far as the Laptev Sea (Fig. 5i), whereas during the 1990s the center of action retracted back westward (Fig. 5ii); broadly similar eastward extents of the Icelandic low were observed in the 2000s and 2010s, respectively (Figs. 5iii and iv). The SCA and PNA also display comparable decadal variability in the extent of their principal centers of action, situated in the west (Figs. 5xxi-xxiv) and east (Figs. 5xiii-xvi) of northern Russia, respectively. The principal regional center of action of both the EAWR and POL patterns was present in all four decades, but its location exhibited considerable variability (Figs. 5xvii-xx and 5xxv-xxviii, respectively). Moreover, in the case of the EA and WP patterns, there were changes in both the sign and location of the associated centers of action over northern Russia (Figs. 5v-viii and 5ix-xii).
Another potential source for the differences between the impact of the atmospheric circulation patterns in the four successive decades is that the pattern itself is dominated by decadal variability. For example, if in a given decade the polarity remained positive then the shape of its distribution would be different from a decade when the polarity varied across both signs, and thus we could expect differences in the regression coefficient for a given significance level. To address this issue, we compare the variance of the four successive decades against the climatology using the Brown-Forsythe test (Brown and Forsythe 1974). This reveals that none of the circulation patterns has a variability in any of the four decades that is significantly different from the climatology (p , 0.10). Moreover, if   Fig. S2ii). Thus, we conclude that decadal changes in climate-circulation pattern relationships across northern Russia result predominantly from decadal changes in the anomalous circulation associated with that pattern. c. Decadal variability of the spatial impact of the circulation patterns As a measure of this variability, we next examine the spatial patterns of temporal consistency. This is simply a measure of the probability that a statistically significant correlation between an atmospheric circulation pattern and SAT or PPN occurs in an arbitrary decade. Maps of the temporal consistency, derived from the 31 decades within the 1980-2019 period and using a significance level of p , 0.10, are shown in Fig. 7. Thus, if a significant correlation occurs in 20 of these 31 decades at a location, then the temporal consistency of the relationship there is ;65%. Note that when significant circulations of both signs occur the temporal consistency refers only to the most frequent pattern. Figure 2ii indicates that there is a significant 40-yr climatological relationship between the NAO and SAT that extends across northern Russia as far east as the Verkhoyanskiy Mountains. However, Fig. 7i reveals that the proportion of individual decades for which a significant relationship exists is considerably less than the maximum possible across the majority of this region. The most consistent significant decadal impact of the NAO on SAT is observed in the Northern European Plain, where the temporal consistency is ;75%. Indeed, there is a small area where the NAO significantly impacts SAT in more than 90% of decades (i.e., where the relationship is essentially temporally invariant). The temporal consistency generally decreases rapidly to the north and especially east, although there are regions, such as the Kola Peninsula and West Siberian Plateau, where it exceeds 60%. The equivalent plot for PPN (Fig. 7ii) indicates that the spatial extent where the NAO-PPN temporal consistency is at a given level is markedly smaller than for NAO and SAT: the only regions where the temporal consistency is more than 60% are the southern Kola Peninsula and a small part of Karelia. With a few exceptions, such as the Taymyr Peninsula, areas where the NAO-PPN relationship has a temporal consistency greater than 30% only extend as far east as ;458E.

2) EAST ATLANTIC
Despite the significant EA-SAT climatological relationship between 458 and 1208E (Fig. 2v), there are very few terrestrial areas where the temporal consistency of the relationship exceeds 30% (Fig. 7iii). Interestingly, there is not a single decade in 1980-2019 when there is a significant relationship between EA and SAT across the majority of the North European Plain and West Siberian Plateau. There are similarly few areas of northern Russia where a significant EA-PPN relationship occurs in more than 30% of decades (Fig. 7iv). They include a spatially coherent region in the West Siberian Plateau south of ;628N, while smaller areas encompass Zemlya Frantsa Iosifa and Severnaya Zemlya.

3) WEST PACIFIC
The WP has the most consistent significant impact on SAT in the mountainous regions of Far East Siberia. Here, the FIG. 6. Summary of the spatial distribution of significant winter decadal relationships (p , 0.10) in four consecutive decades-1) 1980s, 2) 1990s, 3) 2000s, and 4) 2010s-between the circulation patterns and (a) SAT and (b) PPN across selected regions of northern Russia (cf. Fig. 1). Red and blue indicate a significant positive and negative relationship, respectively, while white indicates no significant relationship in that decade. temporal consistency generally exceeds 30%, and increases to more than 60% in some areas (Fig. 7v). Interestingly, a significant decadal relationship between the WP and SAT is also observed in the Barents Sea region during more than 30% of the period examined. In contrast, there are very few decades when WP has any significant impact on PPN in Far East Siberia or indeed across the whole of northern Russia (Fig. 7vi).

4) PACIFIC-NORTH AMERICA
There are no areas of northern Russia where the relationship between the PNA pattern and SAT has a temporal consistency of more than 30% (Fig. 7vii), not unexpected given the climatology (Fig. 2xi). Unsurprisingly, the areas with the most consistent PNA impact on SAT are found in Far East Siberia, the region closest to the PNA center of action in the North Pacific. The land areas where the PNA-PPN relationship is significant in more than 30% of decades are limited to parts of Far East Siberia south of ;618N and the Taymyr Peninsula (Fig. 7viii).

5) EAST ATLANTIC/WESTERN RUSSIA
The primary region where significant relationships between the EAWR pattern and both SAT and PPN have a temporal consistency greater than 30% is the eastern side of the Central Siberian Plateau at ;1208E (cf. Figs. 7ix and 7x). In addition, the EAWR has a significant decadal relationship with PPN in the northern Kola Peninsula for more than 30% of the time.

6) SCANDINAVIAN
The extent of northern Russia where the temporal consistency between the SCA and SAT is greater than 30% aligns approximately with that where there is a significant climatological relationship (cf. Figs. 7xi and 2xvii), an observation that applies to other patterns too. However, across the majority of this region the temporal consistency is more than 60%, predominantly in the West Siberian Plains and Central Siberia Plateau. In the latter region the SCA-SAT relationship is more temporally consistent than NAO-SAT. Regarding PPN, Fig. 7xii indicates that the area where the temporal consistency exceeds 60% is shifted west relative to SAT and in parts of the North European Plain the SCA-PPN relationship is temporally invariant. Similar to the climatology, at a given latitude of the North European Plain the temporal consistency of the SCA-PPN relationship is generally greater than for SCA-SAT (Figs. 7xi and xii).

7) POLAR/EURASIA
The temporal consistency of winter decadal POL-SAT and POL-PPN relationships in northern Russia are shown in Figs. 7xiii and 7xiv, respectively. For SAT, areas with a temporal consistency of 30% comprise the Central Siberian Plateau south of ;658N, while for PPN they consist of the western and eastern flanks of the same region but with reduced temporal consistency in the center. For both parameters there are small areas where the proportion of decades with a significant relationship with POL is greater than 60%. Significant decadal relationships between POL and SAT and/or PPN that have a temporal consistency exceeding 30% are also observed over some Arctic island archipelagos.
Next, we determine the ''dominant'' circulation pattern, defined as having the highest temporal consistency in a region, and examine how that varies spatially across northern Russia. Figure 8i reveals that the NAO is the dominant pattern affecting winter SAT west of the Urals, and also in some Arctic coastal areas to the east, including the Yamal and Taymyr Peninsulas, and over the Barents and Kara Seas and Novaya Zemlya. Farther east, as far as ;1108E, the SCA is generally the dominant pattern, including over much of the Central Siberian Plateau. Immediately to the east of this area and south of ;648N, the POL becomes dominant, and is also the most temporally consistent pattern affecting SAT over the Laptev Sea and both Zemlya Frantsa Iosifa and Sevenaya Zemlya. In the mountainous regions of Far East Siberia the WP is dominant. While Fig. 7 indicates that no circulation pattern has a temporal consistency exceeding 30% in the farthest east region of northern Russia, it is actually the Atlantic-focused EA pattern that is the dominant pattern affecting SAT there. Moreover, Fig. 8i demonstrates that there are no spatially coherent areas where either the PNA or EAWR is the dominant pattern influencing SAT in northern Russia.
The spatial distribution of the dominant circulation pattern influencing winter PPN is shown in Fig. 8ii. In general, it is less clearly defined than for SAT, especially over the ocean and Arctic islands, but there are still distinct spatial patterns. The   Fig. 8ii).
The most extreme response of the observed marked decadal variability in the impact of atmospheric circulation patterns is when statistically significant relationships of both signs arise. While such areas do exist, being most prevalent in Far East Siberia, a false discovery rate test (e.g., Wilks 2006) reveals that the significant majority probably arise through chance.

d. Extending the length of the analysis
One shortcoming of using ERA5 for this analysis is that 40 years is a relatively short time period over which to examine decadal variability. Given that the CPC indices of the atmospheric circulation patterns begin in 1950, there is the potential to use the JRA-55, for which data are available from 1958 to the present, to extend the length of the analysis to 61 years and hence determine whether our findings remain valid over this longer time period. Figure S4 shows the temporal consistency of the significant decadal relationships between the seven circulation patterns and SAT and PPN based on JRA-55 data, for the 1980-2019 overlap period with ERA5. The two reanalyses are generally very similar, with no fundamental differences, indicating that the ERA5 results are robust. The regions of high temporal consistency (.60%) tend to be more spatially coherent in JRA-55, likely due to the poorer spatial resolution of this reanalysis. This means that some local details are missed, such as the clear change in temporal consistency between the POL pattern and SAT associated with the Verhoyansk Mountains seen in the ERA5 data (cf. Fig. 7xiii and Fig. S4xiii). However, overall the two reanalyses are sufficiently similar to justify using the JRA-55 as a direct comparison in order to increase the length of the analysis.
Plots of the winter decadal temporal consistency between the circulation patterns and SAT and PPN based on 61 years of the JRA-55 from 1959 to 2019 are shown in Fig. 9. There are no major disparities between this figure and that for the shorter ERA5based analysis (Fig. 7). Nonetheless, there are a few noteworthy differences. For example, the JRA-55 data reveal reduced temporal consistency between the NAO and SAT over the Northern European Plain and West Siberian Plains, although it increases over the Yamal Peninsula (cf. Figs. 7i and 9i). Conversely, the temporal consistency between the NAO and PPN in these regions is generally greater in the longer time series (cf. Figs. 7ii and 9ii). This is also true of the SCA relationships with climate: in particular there is a clear reduction (increase) in temporal consistency between the SCA and SAT (PPN) over the West Siberian Plains (cf. Figs. 7xi,xii and 9xi,xii). The region in the mountains of Far East Siberia where the temporal consistency between the WP and SAT is more than 30% is expanded in the JRA-55 data and now reaches .90%, but the higher consistency in the Barents Sea has disappeared (cf. Figs. 7v and 9v). The longer time series also demonstrates reductions in the extent of regions where the temporal consistency of both POL-SAT and POL-PPN is greater than 30%. In contrast, the temporal consistency of regional climate relationships with the EA, PNA, and EAWR is generally slightly higher in the JRA-55 data but there remain relatively few areas that exceed 30%.

Discussion and conclusions
In this paper we have investigated the decadal variability of the impact of seven atmospheric circulation patterns on the winter climate of northern Russia, as represented by SAT and PPN. Our results provide a measure against which climate model output for recent decades can be assessed (e.g., Gong et al. 2017), and future projections compared and estimates of their uncertainty improved.
The analysis, based on 40 years of ERA5 data, has demonstrated that there is indeed marked decadal variability in the spatial locations and extent of northern Russia over which these patterns exert a statistically significant influence on climate, often concomitant with long-term significant relationships. Examination of four consecutive decades established that some of this variability is a direct response to the changing location of a pattern's center(s) of action, which affects the spatial distribution of advection of heat and moisture into and across the region. These changes are likely a response to both natural variability within the climate system and external anthropogenic forcing (e.g., Dong et al. 2010;L. Wang et al. 2019). The temporal consistency of significant winter decadal relationships between each circulation pattern and SAT or PPN was determined by analyzing each of the 31 decades contained within the 40-yr study period. Regions of northern Russia with climatological, statistically significant relationships between atmospheric circulation patterns and climate approximately match those with a decadal temporal consistency of only 30%, and thus the climatology masks considerable decadal variability. Extending the period examined from 40 to 61 years using JRA-55 did not significantly alter any of our conclusions, although there are limited expansions and/or contractions of areas with higher temporal consistency for some circulation pattern-climate relationships. Our principal findings regarding each of the atmospheric circulation patterns are as follows: In this study we have examined the changes in decadal impact of each atmospheric circulation pattern in isolation. However, a small number of studies have investigated the changing locations of the north-south dipole of the NAO and shown that its interplay with the EA and SCA patterns exerts one of the primary controls governing the position of its centers of action (Moore et al. 2013;Comas-Bru and McDermott 2014;Mellado-Cano et al. 2019) and subsequent impacts on, for example, the carbon cycle (Bastos et al. 2016) and precipitation (Álvarez-García et al. 2019). Based on the CPC sign conventions used in the present study, when the NAO and EA are of the same sign there is a southwest movement in the centers of action while when the NAO and SCA are of the same sign there is a counterclockwise rotation of the centers (and vice versa). If we examine the 1980s as an example, during this decade the most frequent winter NAO-EA combination was positive NAO and negative EA and similarly the most frequent NAO-SCA combination was positive NAO and negative SCA. According to the two premises above, the northerly NAO center of action should have been negative and shifted both to the northeast movement and rotated clockwise relative to its climatological position. Comparing Figs. 2i and 5i indicates that this is indeed the case in the 1980s. The other three successive decades also follow the premises reasonably well, with the NAO-EA interplay appearing to have a greater influence than that between the NAO and SCA, possibly due to the inversely weighted contribution to regional SLP variability between the NAO and SCA, as noted by Kislov and Matveeva (2020).
Given that we have established that the magnitude of the eastward extension of the northern NAO center has a major impact on the winter climate of much of northern Russia, such interactions between these and different combinations of circulation patterns are likely to contribute to the overall climate variability of the region, including the persistence of SAT anomalies into the following spring Wu and Chen 2020). Therefore, future work will use machine-learning techniques, based on nonlinear regression analysis, to describe the interplay between the various different patterns and to quantify the extent to which it contributes to the marked winter decadal-scale variability in the impact of the individual patterns on northern Russian climate as described in this study.