INTERNATIONAL JOURNAL OF CLIMATOLOGY
Int. J. Climatol. 24: 33–44 (2004)
Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/joc.984

CHANGES IN THE SUB-DECADAL COVARIABILITY BETWEENNORTHERN HEMISPHERE SNOW COVER AND THE GENERALCIRCULATION OF THE ATMOSPHEREKAZUYUKI SAITO,a, * TETSUZO YASUNARIa,b and JUDAH COHENc
a Frontier Research System for Global Change, 3173-25 Showa-machi Kanazawa-ku, Yokohama City, Kanagawa, 236-0001, Japanb Hydrospheric Atmospheric Research Center, Nagoya University, Furo-cho, Chigusa-ku, Nagoya City, Aichi, 464-8601, Japanc Atmospheric and Environmental Research, Inc., Lexington, MA 02421, USA

Received 23 June 2003
Revised 21 October 2003
Accepted 21 October 2003

ABSTRACT
Details of the sub-decadal covariability relationship between continental snow cover extent anomalies and the dominantmode of atmospheric variability, referred to as the Arctic oscillation (AO) or North Atlantic oscillation (NAO), for theperiod 1971–2001 are explored. On the seasonal time scale, the winter AO is found to be significantly correlated with thepreceding autumn Eurasian snow cover (SNCEUR ) throughout the period observed. Consistent with this finding, SNCEUR
variability led the AO variability on the sub-decadal time scale in the early half of the record. However, starting in
the mid 1980s, the AO and SNCEUR vary in phase. Analyses of the seasonal relationship and persistence of snow andatmospheric variables illustrate a phase shift in the sub-decadal variability between the AO and SNCEUR due to the lossof autumn–winter SNCEUR autocorrelation replaced by a significant winter–spring persistence and the emergence of aconcurrent SNCEUR –AO connection in winter and spring. Similar analysis shows that the sub-decadal NAO variation ismostly described by the fluctuation in summer North American snow cover. Copyright  2004 Royal MeteorologicalSociety.
KEY WORDS: Northern Hemisphere; snow–climate interaction; interannual variability; atmospheric general circulation; observational

data analysis

1. INTRODUCTION
The dominant mode of Northern Hemisphere (NH) atmospheric variability, the North Atlantic oscillation(NAO) or Arctic oscillation (AO), is known to have significant spectral peak(s) at interannual and/or
longer time scales (Marshall et al., 2001; Robertson, 2001). Feldstein (2002) demonstrated that the increased
interannual variance of the dominant mode of atmospheric variability after the 1960s is larger than expected
if it is solely driven by internal (intraseasonal) fluctuations of the atmosphere, suggesting coupling with othercomponents of the climate system (hydro- and/or cryo-sphere) and/or forcings external to the system. Links
and covariabilities of the interannual fluctuation of the dominant mode within the system have been found insea-surface temperature, sea ice, and snow cover (e.g. Czaja and Frankignoul, 1998; Cohen and Entekhabi,
1999; Wang and Ikeda, 2000).
Saito and Cohen (2003) show that the dominant mode and continental snow cover have significant coherency
at the sub-decadal (6–10 year) period, and that snow cover leads the atmosphere by several months. The
results in that previous paper are, however, derived from analysis in the frequency domain and, therefore, do
not provide a detailed profile in the temporal domain. In this article we report a more detailed description of* Correspondence to: Kazuyuki Saito, Frontier Research System for Global Change, 3173-25 Showa-machi Kanazawa-ku, Yokohama
City, Kanagawa, 236-0001, Japan; e-mail: ksaito@jamstec.go.jp

Copyright  2004 Royal Meteorological Society

34

K. SAITO, T. YASUNARI AND J. COHEN

changes in the sub-decadal covariability between the continental snow cover extent and the dominant mode
during the remotely sensed period. The identity of and the difference between the AO and NAO remain
controversial. We use ‘the dominant mode’ as a general term to cover both; however, we will refer to AO andNAO variability individually in association with fluctuations with continental snow cover. Previous studies
have shown (and we will also show) that Eurasian snow cover (SNCEUR ) is more closely related to the AO,
and that North American snow (SNCNA ) is related to the NAO (Cohen and Entekhabi, 1999; Bojariu andGimeno, 2003; Saito and Cohen, 2003; Saunders et al., 2003).

2. DATA AND RESULTS
We use monthly NAO and AO indices, as defined by Hurrell and van Loon (1997) and by Thompson
and Wallace (1998) respectively. Continental snow cover areas for Eurasia and North America (excluding
Greenland) are calculated from satellite observations and provided by the Snow Data Resource Center,
Rutgers University Climate Lab (Robinson et al., 1993). We use National Centers for Environmental
Prediction–National Center for Atmospheric Research reanalysis (Kalnay et al., 1996) for atmosphericvariables. Surface temperatures are provided by the Climate Research Unit, University of East Anglia (Jonesand Moberg, 2003), for land and by the Hadley Centre (Rayner et al., 2002) for ocean.
First, we look at 70–120 month period variation, bandpass-filtered by wavelet analysis with Mexican hat(Marr) wavelet (Torrence and Compo, 1998), of the monthly normalized series of the AO index and theEurasian snow cover extent anomaly (SNCEUR ) from 1971 through to 2001 in Figure 1(a). The e-folding
time for that period range corresponds to from 24.9 to 42.7 months, or about 2 to 3.5 years; both ends of therange are influenced by edge effects due to the finite length of the series analysed, and hence the amplitude
is presumably reduced. The filtered AO and SNCEUR show a distinct change in phase between the first andthe second halves of the remotely sensed period delineated by the mid 1980s (Figure 1(a)). For the first half,
SNCEUR tends to lead the AO variability by several or more months, whereas for the second half the twoseries vary almost in phase.
Next we examine the season most representative of the observed SNCEUR and AO sub-decadal variability.
The two filtered time series shown in Figure 1(a), are presented in Figure 1(b) and (c) respectively, including
a 3-year running mean (thin line) of detrended seasonal anomalies. Until the mid 1980s, for SNCEUR it is thewinter that best resembles the sub-decadal fluctuation, whereas in the later period the warm season (spring
to autumn) accounts for more of the interannual variability (Figure 1(b)). For the AO, it is primarily winter(December to February; DJF) that represents most of the variation in the sub-decadal fluctuations through theperiod. Spring (March to May; MAM) variability only exhibits a strong relation with winter variability after
the mid 1980s (Figure 1(c)).
A corresponding change between the two periods is found in the seasonal relation between SNCEUR and the
AO. Figure 2(a) shows the Spearman (rank) correlation of SNCEUR between successive seasons for running
12-year periods. Rank correlation is used in this study to minimize outlier influence, although Pearson’s linear
method produces substantially the same results. The statistically significant autocorrelation between autumn(September to November; SON) and the following DJF SNCEUR found in the earlier period is lost after the
emergence of a significant DJF–MAM autocorrelation in the mid 1980s. A similar emergence of a significant
DJF–MAM lagged cross-correlation is also observed between the DJF AO and the MAM SNCEUR , andconcurrent cross-correlation between the two variables in winter and spring is consistently seen throughout
the entire period (Figure 2(b)).

3. DISCUSSION
The changes in the seasonal autocorrelation of SNCEUR and its correlation with the AO, as indicated in
Figure 2, may be responsible for the changes in the sub-decadal covariability shown in Figure 1(a). Given
the above results, we present the following hypothesis to explain the difference in the covariability between
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

CHANGES IN SUB-DECADAL SNOW-ATMOSPHERE AVAILABILITY

35

0.2

-0.2

0.1

-0.1

0.0

0.0
0.1

-0.1
AO
SNCEUR

-0.2
1970

AO

SNCEUR(a) 70--120-mo filtered (Marr): SNCEUR and AO

1975

1980

1985

1990

1995

0.2

2000

Year
(b) SNCEUR
0.2
0.1
0.0
-0.1

JJA
DJF
Filtered

-0.2
1970

1975

1980

1985

1990

1995

2000

Year
(c) AO
-0.2
-0.1
0.0
0.1

MAM
DJF
Filtered

0.2
1970

1975

1980

1985

1990

1995

2000

Year

Figure 1. (a) 70–120 month bandpass-filtered monthly anomaly time series of Eurasian snow cover (solid; SNCEUR ) and the AO (dashed;scale at right). (b) 70–120 month bandpass-filtered monthly anomaly (thick) and 3-year running mean of detrended JJA (long-dashed)
and DJF (dash-dotted) seasonal anomaly of SNCEUR . (c) As (b), but for the AO (thick) and its MAM (dashed) and DJF (dash-dotted)
seasonal anomaly. This figure is available in colour online at http://www.interscience.wiley.com/ijoc

SNCEUR anomalies and the AO for the two periods. For the earlier period, the autumn SNCEUR persisted into
winter while also leading the winter AO variability (Figure 2(a) and (b)). Also, during the earlier period, thespring SNCEUR had only a weak relation with the winter AO and winter SNCEUR (Figure 2(a) and (b)). This
resulted in SNCEUR leading the AO by several or more months through their mutual oscillations (Figure 1(a)).
For the later period, although autumn SNCEUR still leads the winter AO, it is no longer correlated strongly
with the winter SNCEUR ; rather, winter SNCEUR exhibits strong persistence into the spring. Spring SNCEUR is
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

36

K. SAITO, T. YASUNARI AND J. COHEN

1.0(a) SNCEUR - SNCEUR

0.8

Correlation

0.6
0.4
0.2
0.0
-0.2

SON - DJF
DJF - MAM

-0.4

1985
1990
1980
Year (middle of 12-year period)

1995(b) SNCEUR - AO
0.4
0.2

SON - DJF
DJF - DJF
MAM - DJF
MAM - MAM

Correlation

0.0
-0.2
-0.4
-0.6
-0.8
-1.0

1985
1990
1980
Year (middle of 12-year period)

1995

Figure 2. (a) Rank correlation of SNCEUR between autumn and winter (solid) and winter and spring (dashed) for running 12-year period,
plotted at the middle of each period. Filled (open) diamonds denote significant at 5% (10%) confidence level. (b) As (a), but for the AO
and SNCEUR between different combination of seasons. This figure is available in colour online at http://www.interscience.wiley.com/ijoc

now strongly connected to the winter and spring AO (Figure 2(a) and (b)); therefore, sub-decadal variability
for both climate variables is now dominated by winter to spring variations (Figure 1(b) and (c)), resulting intheir variation in phase (Figure 1(a)).
Still an important question remains: Why did the persistence of the SNCEUR and winter AO–spring SNCEUR
relation change between the two periods whereas the autumn SNCEUR –winter AO relationship held throughout
the entire period despite small-amplitude fluctuations (Figure 2(b))?
We select two 12-year periods, i.e. 1974–85, and 1986–97, for further analysis and present selected environmental fields for each period in an effort to illustrate and explain the observed change and difference between
the periods delineated by the mid 1980s. Figure 3 shows composite differences between five extreme yearsof each sign with respect to SON SNCEUR to show the location and magnitude of maximum differences
between both extremes. The month shown in the figure denotes the centre of a 1–2–1 weight applied to threeadjacent months. The time series and environmental fields are detrended within a semi-period. Significanceis tested using a t-test. Field significance is tested against a null distribution of areas of significant local
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

37

CHANGES IN SUB-DECADAL SNOW-ATMOSPHERE AVAILABILITY(a) Cmpst: WAFz100 (Dec) by SNCEUR (Oct) [74-85](c) Cmpst: WAFz100 (Dec) by SNCEUR (Oct) [86-97]
0

-2

0
0

-2

2

0

-2

2

4

-4

-4
-2

2

0

-4

2
0

2 -6
--4

4
0

-2 02
4
6

2

-6

2

-4
-2

0

68

0

2

0

8
10

0

--42
-6

6

-2

4

0

2

4

4

2
2

2

-2

6
4

810

0

12

0

-2

-2

-4

2

0

0

-8
-6

0

0

2
4

0

-4

0
2

0
-2

2

-4
0

-.99

-.95

-.90

.90

.95

.99

(b) Cmpst: Snowcover (Jan) by SNCEUR (Oct) [74-85]

-.99

-.95

-.90

.90

.95

.99

(d) Cmpst: Snowcover (Jan) by SNCEUR (Oct) [86-97]

0

4
84

-4

-4
4

0

0

0

0

12

4

8

0

0

16

-4

8

20

0

4

16

0
0

2
116

4

0

0

8

4

16
812

0

0

0

0

0

0

4
-4

0

0

-8

12
8

-.99

-.95

-.90

.90

0

0

.95

.99

-.99

-.95

-.90

0

4

0

0

0

0

4
8

12

.90

.95

.99

Figure 3. (a) Composite of December vertical WAF (WAFz ) at 100 hPa between five extensive and sparse October SNCEUR for 1974–85.
Contour interval is 2 m2 s−2 , and positive (negative) values are contoured in thick (thin). Significance at 10, 5, 1% level in positive
and negative directions is shown by grey-scale gradation. WAFz field and SNCEUR series are 1–2–1 weighted with those of adjacentmonths. (b) As (a) but for January snow cover. Contour interval is 4%. (c) As (a) and (d) as (b), except for the period 1986–97

correlation at the 5% significance level, which is constructed by the time series permutated in a random
order. The domain is defined to include focused regions of each relationship. The results are summarized inTable I.
Previous studies have demonstrated that autumn SNCEUR influences the winter AO by modulating the
upward propagation of the stationary Rossby waves as diagnosed by the wave activity flux (Plumb, 1985;
WAFz ), both with observational (Saito et al., 2001) and numerical modelling studies (Gong et al., 2003).
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

38

K. SAITO, T. YASUNARI AND J. COHEN

Table I. Results of field significance test. Percentage area of local significance denotes total
areas of grid cells that show significant local correlation at 5% significance level within thedomainFigure no.

3

5

Figure part

DomainPercentage areaof local significancea

(a)
(c)
(b)
(d)

[80–180 ° E, 40–65 ° N]

7.38+
7.96+
6.07∗
2.46

(a)
(d)
(b)
(e)
(c)
(f)

[0–360 ° E, 30–85 ° N]

[0–180 ° E, 30–70 ° N]

[0–180 ° E, 30–70 ° N]
[0–180 ° E, 30–70 ° N]

17.3∗
21.1∗
4.53
19.2∗
3.72
6.56∗

a +: within upper 10%; ∗ : within upper 5%.

Changes to the atmospheric circulation, forced by anomalous SNCEUR , are consistent with existing theories
on troposphere–stratosphere coupling (Zhou et al., 2002). The upward WAFz anomaly, field significant overnortheastern Eurasia at 10% (Table I), is observed in late autumn–early winter for both periods. Although
the location of the anomalous upward WAFz does vary between the two semi-observational periods, bothregions of anomalous flux lie within the major Siberian wave train of upward flux (Plumb, 1985), consistent
with previously demonstrated results. Therefore, during both semi-periods, variations in autumn SNCEUR are
likely to initiate changes in the AO consistent with known troposphere–stratosphere coupling.
We next present snow-cover expansion from October through to January for both semi-periods by acomposite difference of January snow cover with respect to October SNCEUR (Figure 3(b) and (d)). SON–DJF
persistence in the earlier semi-period is due to an organized, almost zonally uniform advance of snow cover.
These areas cover most of the regions of ‘active’ snow-cover fluctuation, as shown in figure 3 of Clark et al.
(1999). In contrast, during the latter period, no organized expansion is seen except for a limited region ofcentral Europe. Clark et al. (1999) argue that the atmosphere controls the Eurasian snow cover extent during
winter (DJF) and demonstrate that three teleconnection patterns have distinct and different regional forcingof the extent of Eurasian snow cover.
We define HA as the area-weighted mean of the normalized geopotential height anomaly Z north of 60 ° N
at each pressure level (Cohen et al., 2002). HA serves as an excellent proxy for the AO (see Figure 4(a)), yet
it also captures the dominant regional variability of height anomaly at high latitudes, including the strength ofthe polar vortex at each pressure level. The monthly series of the AO, defined as the first principal componentof the geopotential height anomaly field, and HA at the same pressure level have substantial coherency onmonthly to sub-decadal time scales (not shown).
Differences in tropospheric persistence between the two periods are shown by HA500 (Figure 4(b)), i.e. HAdefined at 500 hPa, which exhibits similar seasonal characteristics as SNCEUR (Figure 2(a)). The 500 hPa levelis chosen for display, but similar results are obtained for 700 hPa through to 200 hPa. Through composite
analysis, HA500 variability in the earlier semi-period is found to relate to a pattern similar to the EurasianType 1 pattern (Barnston and Livezey, 1987; EU1) from autumn to winter in the pressure/height fields, andto anomalous surface temperature over northeastern Eurasia (not shown). Temperatures in this region are
known to be more sensitive to the EU1 pattern than southwestern Eurasia (Clark et al., 1999). The strong
autocorrelation in autumn–winter SNCEUR before the mid 1980s may be due to the persistence of the EU1
pattern and its influence on surface temperatures. In the later semi-period the composite difference of winter
HA500 is reminiscent of the NAO pattern instead, whose impact on snow cover is confined to western Eurasia(not shown).
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

CHANGES IN SUB-DECADAL SNOW-ATMOSPHERE AVAILABILITY

39(a) AO - HA1000

1.0

Correlation

0.8
0.6
0.4

DJF - DJF
MAM - MAM

0.2
0.0

1980
1985
1990
Year (middle of 12-year period)(b) HA500 - HA500

1.0
0.8

SON - DJF
DJF - MAM

0.6
Correlation

1995

0.4
0.2
0.0
-0.2
-0.4
1980
1985
1990
Year (middle of 12-year period)

1995(c) NAO - HA1000

1.0

Correlation

0.8
0.6
0.4
DJF - DJF
MAM - MAM

0.2
0.0

1980

1985

1990

1995

Year (middle of 12-year period)

Figure 4. (a) Concurrent rank correlation between the AO and area-weighted geopotential height anomaly (HA) at 1000 hPa, HA1000 ,in winter (solid) and spring (dashed). (b) As (a) except for HA at 500 hPa, HA500 . (c) As (a), except between the NAO and HA1000 .
This figure is available in colour online at http://www.interscience.wiley.com/ijoc

Next, we discuss the strong winter–spring relation that emerged after the mid 1980s. Figure 4(c) shows
the concurrent seasonal correlation between the NAO and HA1000 in winter and spring. The figure illustrates
that the AO–NAO connection is strongest in the cold season in the earlier semi-period whereas the relation
between the two in the spring is only significant after the mid 1980s. This difference in the warmer season
connection may be due to the snow–atmosphere relation, which is found only in the later period. Spatial
patterns of sea-level pressure associated with the March AO are shown in Figure 5(a) and (d). In the Atlantic
sector, before the 1980s, the March AO pattern is limited to the North Atlantic basin, having little influence
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

40

K. SAITO, T. YASUNARI AND J. COHEN(a) Cmpst: SLP (Mar) by AO (Mar) [74-85](d) Cmpst: SLP (Mar) by AO (Mar) [86-97]

0

0

0

0

0

-2

-2

-8

-6

-8
-6

-4

4

0-2
-4

2

2

2
0
-2

-6

4

0

2

2

0

0

24

6
8

0

2

0

4

4
2

0

-.95

-.90

.90

.95

.99

-.99

-.95

.90

.95

.99

.8

0.0

-0

-0.8

0.0

-.90

(e) Cmpst: Ts+SST (Mar) by AO (Mar) [86-97]

1.

0.
0
0.8

0.0

6
0.8

0.8

0.0

0.
8

-0.8

0.0

0.0

0.
0

0.

00.

0.0

-1.6

0.8

8

-0.8

0.0
0.0

2.3
4.2

-.99

(b) Cmpst: Ts+SST (Mar) by AO (Mar) [74-85]

-0.8
0.

0

0.0

8
0.

-.90

.90

.95

.99

(c) Cmpst: Snowcover (Mar) by AO (Mar) [74-85]

-.90

.90

.95

.99

(f) Cmpst: Snowcover (Mar) by AO (Mar) [86-97]

0

-8

0

0

-8

-8

0

8

0

-8

0

-8

0

0

0

0
0

0
0

0

0

0

-.95

0

-8

-.99

0

-.95

-8
-216
4

-.99

0

0.

0.0

0

0

0

8

-8

0

0

-.99

-.95

-.90

.90

.95

.99

-.99

-.95

-.90

.90

.95

.99

Figure 5. (a) Same as Figure 3(a) but for March sea-level pressure and the March AO. Contour interval is 2 hPa. (b) As (a), but for the
combined March surface temperature over land (Ts ) and ocean (SST). Contour interval is 0.8 K. (c) As (a), but for March snow cover.
Contour interval is 8%. (d) As (a), (e) as (b) and (f) as (c), except for the period 1986–97
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

CHANGES IN SUB-DECADAL SNOW-ATMOSPHERE AVAILABILITY

41

over Eurasia, whereas after the mid 1980s the March AO pattern extends into western Eurasia. This pattern ispersistent, beginning in January in the later semi-period (not shown), and implies possible interaction between
snow extent and the atmospheric circulation associated with the dominant mode variability from winter tospring.
Previous studies have established that the observed positive trend in the dominant mode is associated withthe surface temperature Ts increase over the northern part of Eurasia, especially northern Europe and theCentral Siberian Plateau (area around 100 ° E; Hurrell and van Loon, 1997; Clark et al., 1999; Trigo et al.,2002). As expected from the difference in the spring AO pattern between the two semi-periods (Figure 5(a)and (d)), only for the later semi-period are significant areas of warm Ts (associated with the AO) observed
over northern Eurasia, with the strongest signal over northeastern Europe (along 30 ° E longitude) and theCentral Siberian Plateau (Figure 5(b) and (e)). Consequently, a significant AO–snow-extent relation is only
observed in the later semi-period. During the later period, snow-cover extent is mostly influenced over the area
along 30 ° E where snow-cover extent is sensitive to intraseasonal temperature variance (Figure 5(c) and (f)).
The composite maps from January through to March produce substantially the same results, except for thegradual northeastward retreat of those areas of snow cover influenced by variability in the AO (not shown).
The duration of the relation between the AO, Ts and snow cover is likely responsible for the significant
winter AO–spring SNCEUR relation depicted in Figure 2(b) for the later semi-period. Figure 6 summarizes
schematically the changes described in seasonal relationships in the sub-decadal covariability between SNCEUR
and the AO before and after the mid 1980s.
Little is understood about snow–climate interactions that span from the warm to cold seasons; however, it
most probably differs from those attributed to autumn SNCEUR and the winter climate, since the troposphere

Tropospheric
Height Anomaly
ArcticOscillation

Before themid 1980s

Eurasiansnow cover

AutumnWinterSpring

Tropospheric
Height Anomaly
After themid 1980s

ArcticOscillation
Eurasiansnow cover: important season(s): season(s) representing the sub-decadal variability: seasonal SNC EUR persistence: influence of autumn SNC EUR on the winter AO: significant association between the AO and SNC EUR

Figure 6. Schematic diagram of changes of the seasonal relationships in the sub-decadal covariability between Eurasian snow coveranomaly (SNCEUR ) and the AO. Dark (light) shading denotes strong (modest) concurrent relation between the AO and NAO
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

42

K. SAITO, T. YASUNARI AND J. COHEN

0.2

-0.2

0.1

-0.1

0.0

0.0
0.1

-0.1
-0.2
1970

NAO

SNCNA

(a) 70--120-mo filtered (Marr): SNCNA and NAO

NAO
SNCNA
1975

0.2
1980

1985
Year

1990

1995

2000

1995

2000

1995

2000

(b) SNCNA
0.2
0.1
0.0
-0.1
-0.2
1970

JJA
DJF
Filtered
1975

1980

1990
1985
Year
(c) NAO

-0.2
-0.1
0.0
0.1
0.2
1970

MAM
JJA (SNCNA)
DJF
Filtered
1975

1980

1990
1985
Year

Figure 7. (a) As Figure 1(a), except for North American snow cover (solid; SNCNA ) and the NAO (dashed; scale at right). (b) As
Figure 1(b), but for SNCNA . (c) As (b), but for the NAO (thick) and its DJF (dash-dotted) and MAM (long-dashed) seasonal anomaly.
JJA SNCNA (dashed) is superimposed. This figure is available in colour online at http://www.interscience.wiley.com/ijoc

and stratosphere are uncoupled in the warm season. Instead, dynamical pathways are most probably confinedto the lower troposphere. Also, no clear evidence was found in this study that snow cover anomalies influence
the atmospheric circulation from winter to spring seasons.
Finally, we look at the sub-decadal variability of the NAO and its relation to the North American snow
cover anomaly (SNCNA ) in Figure 7. The bandpass-filtered NAO and SNCNA show two lead changes: thefirst in late 1970s and the second in the early 1990s (Figure 7(b)). In the first period, the NAO appears
to lead, whereas in the second it lags. In the third period, the NAO leads slightly, but they vary mostly
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

CHANGES IN SUB-DECADAL SNOW-ATMOSPHERE AVAILABILITY

43

in phase. Interannual variability of seasonal SNCNA resembles that of SNCEUR , in that winter SNCNA
dominates the interannual fluctuations in the earlier period but the warm season accounts for more of the
sub-decadal variations in the later period. Concerning the NAO and SNCNA covariability, the sub-decadal
NAO variability (Figure 7(c)), which shares similar characteristics with respect to seasonal dependence as
the AO (Figure 1(c)), is astonishingly well depicted by June to August (JJA) SNCNA (short dashed line in
Figure 7(c); cf. Figure 7(b)). This result suggests that JJA SNCNA variability is closely related to the subdecadal NAO variability, consistent with previous studies (Bamzai, 2002; Bojariu and Gimeno, 2003; Saito
and Cohen, 2003; Saunders et al., 2003).

4. CONCLUSION
We attempt in this study to gain a deeper insight into the covariability between continental snow cover extent
and the dominant mode of atmospheric variability, the AO or NAO, on sub-decadal time scales for the
remotely sensed period. We found that there were changes in the structure of the covariability. A significant
change occurred in the mid 1980s in the relation between SNCEUR and the AO. Before the mid 1980s theautumn SNCEUR –winter AO relation is considered to have been decisive in the sub-decadal covariability,with snow cover leading by several or more months. After the mid 1980s, loss of autumn–winter SNCEUR
persistence and the emergence of a strong AO–SNCEUR connection from winter to spring brought about anin-phase relationship on sub-decadal time scales. For the NAO and SNCNA , the bandpass-filtered series shows
two lead changes in the late 1970s and early 1990s, but the sub-decadal variability of the NAO is closely
related to that of summer SNCNA throughout the remotely sensed period.
The entire period examined in this study did experience a long-term increase of high-latitude Ts and adecrease of continental snow-covered area (cf. figures 3 and 10 in Serreze et al. (2000)). Also, the length of
the analysis period is limited. How these long-term trends may have also impacted the relationship between
snow and the dominant mode of winter atmospheric variability was not explored in our study. It remains forfurther research to determine whether the change detected in the sub-decadal covariability between continentalsnow cover anomalies and atmospheric variability of hemispheric scale is due in part to a long-term increasing
trend of high-latitude Ts and decreasing trend in snow-covered area in recent decades, or as part of internalinter-decadal fluctuations of the climate system. Well-organized general circulation model simulations and/or
observational analysis of more extended length are required to answer this question.
ACKNOWLEDGMENTS

This investigation was partly supported by NSF grant ATM-0124904.
REFERENCES
Bamzai AS. 2002. Relationship of snow cover variability and Arctic oscillation index on a hierarchy of time scales. International.
Journal of Climatology 23: 131–142.
Barnston AG, Livezey RE. 1987. Classification, seasonality and persistence of low-frequency atmospheric circulation patterns. Monthly
Weather Review 115: 1083–1126.
Bojariu R, Gimeno L. 2003. The role of snow cover fluctuations in multiannual NAO persistence. Geophysical Research Letters 30(4):DOI: 10.1029/2002GL015651.
Czaja A, Frankignoul C. 1999. Influence of the North Atlantic SST on the atmospheric circulation Geophysical Research Letters 26:
2969–2972.
Clark MP, Serreze MC, Robinson DA. 1999. Atmospheric controls on Eurasian snow extent. International Journal of Climatology 19:
27–40.
Cohen J, Entekhabi D. 1999. Eurasian snow cover variability and Northern Hemisphere climate predictability. Geophysical Research
Letters 26: 345–348.
Cohen J, Salstein D, Saito K. 2002. A dynamical framework to understand and predict the major Northern Hemisphere mode.
Geophysical Research Letters 29: DOI: 10.1029/2001GL014117.
Feldstein SB. 2002. The recent trend and variance increase of the annular mode. Journal of Climate 15: 88–94.
Gong G, Entekhabi D, Cohen J. 2003. Modeled Northern Hemisphere winter climate response to realistic Siberian snow anomalies.
Journal of Climate 16: 3917–3931.
Hurrell JW, van Loon H. 1997. Decadal variations associated with the North Atlantic oscillation. Climatic Change 36: 301–326.
Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)

44

K. SAITO, T. YASUNARI AND J. COHEN

Jones PD, Moberg A. 2003. Hemispheric and large-scale surface air temperature variations and an update to 2001. Journal of Climate
16: 206–223.
Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L, Iredell M, Saha S, White G, Woollen J, Zhu Y, Chelliah M,
Ebisuzaki W, Higgins W, Janowiak J, Mo KC, Ropelewski C, Wang J, Leetmaa A, Reynolds R, Jenne R, Joseph D. 1996. The
NCEP/NCAR 40-year reanalysis project. Bulletin of the American Meteorological Society 77: 437–471.
Marshall J, Kushnir Y, Battisti D, Chang P, Czaja A, Dickson R, McCartney M, Saravanan R, Visbeck M. 2001. North Atlantic climate
variability: phenomena, impacts and mechanisms. International Journal of Climatology 21: 1863–1898.
Plumb A. 1985. On the three-dimensional propagation of stationary waves. Journal of the Atmospheric Sciences 42: 217–229.
Rayner NA, Parker DE, Horton EB, Folland CK, Alexander LV, Rowell DP, Kent EC, Kaplan A. 2003. Global analysis of sea surface
temperature, sea ice, and night marine air temperature since the late nineteenth century. Journal of Geophysical Research: Atmosphere108(D14): DOI: 10.1029/2002JD002670.
Robertson AW. 2001. Influence of ocean–atmosphere interaction on the Arctic Oscillation in two general circulation models. Journalof Climate 14: 3240–3254.
Robinson DA, Dewey F, Heim Jr R. 1993. Northern Hemispheric snow cover: an update. Bulletin of the American MeteorologicalSociety 74: 1689–1696.
Saito K, Cohen J, Entekhabi D. 2001. Evolution of atmospheric response to early-season Eurasian snow cover anomalies. Monthly
Weather Review 129: 2746–2760.
Saito K, Cohen J. 2003. The potential role of snow cover in forcing interannual variability of the major Northern Hemisphere mode.
Geophysical Research Letters 30(6): DOI: 10.1029/2002GL016341.
Saunders MA, Qian B, Lloyd-Hughes B. 2003. Summer snow extent heralding of the winter North Atlantic oscillation. GeophysicalResearch Letters 30(7): DOI: 10.1029/2002GL016832.
Serreze MC, Walsh JE, Chapin III FS, Osterkamp T, Dyurgerov M, Romanovsky V, Oechel WC, Morison J, Zhang T, Barry RG. 2000.
Observational evidence of recent changes in the northern high-latitude environment. Climatic Change 46: 159–207.
Thompson DW, Wallace JM. 1998. The Arctic oscillation signature in the wintertime geopotential height and temperature fields.
Geophysical Research Letters 25: 1297–1300.
Torrence C, Compo GP. 1998. A practical guide to wavelet analysis. Bulletin of the American Meteorological Society 79: 61–78.
Trigo RM, Osborn TJ, Corte-Real JM. 2002. The North Atlantic oscillation influence on Europe: climate impacts and associated physicalmechanism. Climate Research 20: 9–17.
Wang J, Ikeda M. 2000. Arctic oscillation and Arctic sea-ice oscillation. Geophysical Research Letters 27: 1287–1290.
Zhou S, Miller AJ, Wang J, Angell JK. 2002. Downward-propagating temperature anomalies in the preconditioned polar stratosphere.
Journal of Climate 15: 781–792.

Copyright  2004 Royal Meteorological SocietyInt. J. Climatol. 24: 33–44 (2004)