AMERICAN METEOROLOGICAL SOCIETY Journal of Climate
EARLY ONLINE RELEASE This is a preliminary PDF of the author-produced manuscript that has been peer-reviewed and accepted for publication. Since it is being posted so soon after acceptance, it has not yet been copyedited, formatted, or processed by AMS Publications. This preliminary version of the manuscript may be downloaded, distributed, and cited, but please be aware that there will be visual differences and possibly some content differences between this version and the final published version. The DOI for this manuscript is doi: 10.1175/JCLI-D-16-0214.1 The final published version of this manuscript will replace the preliminary version at the above DOI once it is available. If you would like to cite this EOR in a separate work, please use the following full citation: Zhang, R., R. Zhang, and Z. Zuo, 2017: Impact of Eurasian spring snow decrement on East Asian summer precipitation. J. Climate. doi:10.1175/JCLI-D16-0214.1, in press. © 2017 American Meteorological Society
Manuscript (non-LaTeX)
Click here to download Manuscript (non-LaTeX) manuscript(second revision).docx
Impact of Eurasian spring snow decrement on East Asian summer precipitation
Ruonan Zhang Chinese Academy of Meteorological Sciences, Beijing, China National Climate Center, China Meteorological Administration, Beijing, China
Renhe Zhang* Chinese Academy of Meteorological Sciences, Beijing, China CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing, China
Zhiyan Zuo Chinese Academy of Meteorological Sciences, Beijing, China
Submitted to the Journal of Climate on March 12, 2016 Revised on January 6, 2017
*
Corresponding author:
Dr. Renhe Zhang No. 46, Zhong-Guan-Cun South Avenue, Beijing 100081, China E-mail:
[email protected]
1
1
ABSTRACT: In this study, the relationship between Eurasian spring snow decrement (SSD) and
2
East Asian summer precipitation and related mechanisms were investigated using observational
3
data and the Community Atmospheric Model version 3.1 (CAM3.1). The results show that a
4
west-east dipole pattern in Eurasian SSD anomalies, with a negative center located in the region
5
between eastern Europe and the West Siberia Plain (EEWSP) and a positive center located around
6
Baikal Lake (BL), is significantly associated with East Asian summer precipitation via triggering
7
an anomalous mid-latitude Eurasian wave train. Reduced SSD over EEWSP corresponds to
8
anomalously dry local soil conditions from spring to the following summer, thereby increasing
9
surface heat flux and near-surface temperatures. Similarly, the increase in SSD over BL is
10
accompanied by anomalously low near-surface temperatures. The near-surface thermal anomalies
11
cause an anomalous meridional temperature gradient, which intensifies the lower-level
12
baroclinicity and causes an acceleration of the subtropical westerly jet stream, leading to an
13
enhanced and maintained Eurasian wave train. Additionally, the atmospheric response to changed
14
surface thermal conditions tends to simultaneously increase the local 1000–500 hPa thickness,
15
which further enhances the Eurasian wave train. Consequently, significant wave activity flux
16
anomalies spread from eastern Europe eastward to East Asia and significantly influence the
17
summer precipitation over China, with more rainfall over northeastern China and the Yellow
18
River valley and less rainfall over Inner Mongolia and southern China.
19 20
KEYWORDS: spring
21
activity flux
snow decrement, East Asian summer precipitation, Eurasian wave train, wave
22
2
23
1. Introduction
24
Variations in Eurasian snow cover in winter and spring play a major role in the East Asian climate
25
and precipitation during the following seasons (Gong et al., 2003, 2004; Wu and Kirtman, 2007;
26
Zhang et al., 2008; Wu et al., 2009a; Zuo et al., 2011, 2012; Mu and Zhou, 2012, 2015; Zhang et
27
al., 2013; and many others). The associated mechanisms can be attributed to albedo and
28
hydrological effects, with the latter considered the more important (Barnett et al., 1989; Cohen
29
and Rind, 1991; Gong et al., 2004). Previous studies have used various snow parameters and
30
representations to indicate the influence of Eurasian snow cover on circulation anomalies in the
31
Northern Hemisphere and associated dynamic mechanisms. For example, Yasunari et al. (1991)
32
found that excessive Eurasian snow mass during March may trigger the spreading of Rossby
33
waves from the Eurasian continent to North America, leading to the positive phase in the Pacific
34
North America (PNA) teleconnection. Chen and Song (2000a, b) found that the early spring
35
retreat of snow cover extent in Eurasia absorbs large amounts of heat from the atmosphere,
36
leading to thermal anomalies in the troposphere, which ultimately induce anomalous circulation
37
patterns over the mid-high latitudes of the Eurasian continent.
38
Several land factors, such as soil moisture, soil temperature, and frozen soil have been
39
considered to be the links with the long memories of land surface hydrological and thermal
40
conditions and related feedback to climate (Shinoda, 2001). Some studies have investigated the
41
role of spring soil moisture (Matsuyama et al., 1998; Zuo and Zhang, 2007, 2016; Zhao et al.,
42
2007; Zhang and Zuo, 2011) or frozen soil (Li et al., 2009) in the subsequent summer
43
precipitation over East Asia and provide a theoretical basis for constructing the linkages. Chen
44
and Song (2000a, b) suggested that the persistent soil moisture anomaly produced by early snow
45
cover retreat from spring to summer appears to decrease the geopotential height field, which is 3
46
not favorable for the development of blocking activity. Meanwhile, the soil moisture anomaly
47
could also excite anomalous wave trains propagating from western Europe to East Asia through
48
the geopotential height field, with opposite signs between northern and southern China. Using
49
snow water equivalent (SWE) data, Wu et al. (2009a) examined the wave train response to
50
reduced Eurasian spring SWE. Kripalani et al. (2002) and Wu et al. (2014) found a west-east
51
contrast pattern of the snow cover over Eurasian continent is closely connected to summer
52
precipitation anomalies over East Asia. Furthermore, Mu and Zhou (2012, 2015) demonstrated
53
that the soil temperature and frozen soil anomalies associated with the extent of fresh snow in
54
Eurasia might increase the northerly wind over the Baikal Lake (BL) region, leading to relatively
55
colder conditions in northern East Asia, a strengthened subtropical upper-level jet stream, and a
56
weaker Asian summer monsoon (Bamzai and Marx, 2000; Dash et al., 2005).
57
The aforementioned studies are mainly based on snow cover variability itself. As a matter of
58
fact, spring is a transition season from winter to summer. Spring snow decrement (SSD) is closely
59
related with snowmelt and latent heat sink and simultaneously reflects changes in soil moisture
60
and the hydrological cycle, which influences the local thermal conditions and ultimately leads to
61
climate anomalies both locally and remotely (Walsh et al., 1985; Groisman et al., 1994). Although
62
these possible physical processes are identified, the detailed and quantitative climate impacts of
63
Eurasian SSD on the East Asian climate have not been convincingly and deeply understood. Little
64
attention has been paid to variations in the SSD and its relationship with the East Asian summer
65
climate.
66
The purpose of the present study is to investigate the effect of Eurasian SSD on East Asian
67
summer precipitation and to explore the associated dynamic and thermodynamic processes. The
68
remainder of this paper is organized as follows. In section 2, we describe the data and methods 4
69
used in the present study. Section 3 presents the relationship between Eurasian SSD and East
70
Asian summer rainfall. Section 4 explores the associated dynamic and thermodynamic
71
mechanisms. Section 5 discusses the role of El Niño and Arctic sea ice in the relationship between
72
Eurasian SSD and East Asia summer rainfall, and conclusions are reported in section 6.
73 74
2. Data, methods and model
75
The datasets used in this study include (1) the monthly mean atmospheric data from the National
76
Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP/NCAR)
77
Reanalysis I with a horizontal resolution of 2.5º×2.5º available after 1948 (Kalnay et al., 1996), (2)
78
the monthly mean sea ice concentration (SIC) and sea surface temperature (SST) from the Met
79
Office Hadley Centre with a horizontal resolution of 1º×1º available since 1870 (Rayner et al.,
80
2003), (3) the monthly mean soil moisture from the ERA-Interim reanalysis with a horizontal
81
resolution of 0.75º×0.75º and four layers of 0-7, 7-28, 28-100, 100-289 cm in depth for the period
82
1979–2013 (Balsamo et al., 2009), which can well represent the features of observed soil moisture
83
over East Asia (Zuo and Zhang, 2009; Liu et al., 2014), and (4) the monthly mean snow water
84
equivalent (SWE) data of GlobSnow V2 from the Finnish Meteorological Institute for the period
85
1979–2013 (http://www.globsnow.info/swe/), which are obtained using a combination of
86
ground-based data and satellite microwave radiometer-based measurements and each monthly
87
matrix having 721×721 grid (Solberg et al., 2009). We have converted the raw GlobSnow V2
88
SWE data to regular 1º×1º grid for our analysis. For the majority of summer (June–September),
89
measurable SWE is absent; thus, the data in summer are not credible and not utilized. In addition,
90
SWE data are lacking south of 35ºN; thus, the climate impact of SWE over the Tibetan Plateau is
91
not discussed in our analysis. 5
92
In this study the time period is taken from 1979 to 2013. The methods used in the present study
93
include the empirical orthogonal function (EOF) as well as correlation and regression analyses.
94
The Student’s t-test is used to check the statistical significance for the correlation and regression
95
analyses.
96
The monthly snow decrement can be considered as the difference in SWE between two
97
successive months. Positive and negative snow decrements represent increase and decrease SWE
98
in succeeding month with respect to the previous one. Figure 1 shows the climatological annual
99
cycle of SWE and snow decrement over the Eurasian continent (35º–76ºN, 0º–150ºE) during the
100
period 1980–2013, with the absence of June–September SWE data. The monthly SWE rises
101
considerably from October to the peak in February−March (approximately 61 mm) and then
102
declines sharply in the following months to approximately 18 mm in May. Corresponding to the
103
annual cycle of SWE, snow enhancement can be observed in cold seasons starting in October to
104
the following February, followed by a hiatus in March; finally, the most significant snow
105
reduction can be observed in April (-21 mm) and May (-24 mm). Therefore, the most severe snow
106
reduction occurs mainly in the warming seasons in spring due to the gradually increasing
107
temperature and decreasing snowfall. Therefore, in order to quantitatively depict the variations in
108
snow decrement during the spring, in the present study, the Eurasian spring snow decrement (SSD)
109
is defined as the SWE difference between early spring (February−March) and late spring (May).
110
In this study, the NCAR Community Atmospheric Model version 3.1 (CAM3.1), with a
111
horizontal resolution of T42 (approximately 2.8º×2.8º) and 26 vertical levels (up to 3.7 hPa)
112
(Collins et al., 2006), was employed to better understand the dynamic and thermodynamic
113
influences of Eurasian SSD on circulation anomalies over Eurasia, particularly East Asia and
114
China. As shown in Collins et al. (2006) and Hack et al. (2006), this model is able to reproduce 6
115
remarkably realistic atmospheric circulation and broad climatological characteristics of East Asia
116
summer climate.
117 118
3. Results
119
3.1 Eurasian SSD and its relationship with East Asian summer precipitation
120
An empirical orthogonal function (EOF) analysis was applied to the Eurasian SSD. The first three
121
EOF modes of the SSD (see Supplementary Fig. S1) account for 22.2%, 14.3% and 7.4% of the
122
total variance, respectively. According to North et al. (1982), all these three modes are well
123
separated with each other. Figure 2 shows the spatiotemporal features of the second EOF mode
124
(EOF2) of Eurasian SSD. The spatial pattern of EOF2 is characterized by a west-east dipole
125
pattern, with a negative center located in the region between eastern Europe and the West Siberia
126
Plain (EEWSP) and a positive center in the vicinity of Baikal Lake (BL). Additionally, two other
127
positive centers over eastern Europe and the North Siberia Plain are also identified (Figure 2a).
128
This dipole pattern somewhat resembles those of previous studies that are based on spring snow
129
cover fraction (Yim et al., 2010; Wu et al., 2014). Wu and Kirtman (2007) also considered the West
130
Siberia snow cover anomaly during spring to be a notable precursor for China seasonal
131
precipitation. As shown by the principal component of EOF2 (PC2) (Figure 2b), this pattern
132
exhibits significant interannual and interdecadal variations. To quantify the Eurasian SSD
133
variability and its association with the East Asian summer climate, we defined a SSD index (SSDI)
134
as the difference in area-averaged SSD between the domains with positive values (46º–58ºN, 96º–
135
138ºE) and negative values (54º–68ºN, 48º–84ºE). The interannual and interdecadal variations of
136
the SSDI are highly consistent with those of the PC2, and their correlation coefficient is as high as
137
0.91, indicating that the SSDI index can provide a feasible representation of the spatiotemporal 7
138
variations in Eurasian SSD during 1980–2013.
139
Corresponding to the variations in Eurasian SSD, the related East Asian summer precipitation
140
anomaly is investigated by correlation analysis with the SSDI index. As shown in Figure 3a, the
141
correlation coefficient pattern reveals a meridional quadrupole structure prevailing in the regions
142
between northern and southern East Asia, with excessive precipitation over regions west of BL,
143
northeastern China and the Yellow River valley and deficient precipitation over Inner Mongolia
144
and southern China. We checked the first three EOF modes of the East Asian summer
145
precipitation (Supplementary Fig. S2). The first three modes account for 12.9%, 11.1% and 8.0%
146
of the total variance, respectively, and they are well separated with each other according to North
147
et al. (1982). The meridional quadrupole pattern strongly resembles that of the spatial distribution
148
of the third EOF mode (EOF3) of East Asian summer precipitation (Figure 3b), and their spatial
149
correlation coefficient is as high as 0.67. Furthermore, the variations in PC3 of the East Asian
150
summer precipitation is similar to that of the SSDI index during 1980–2013 (Figure 3c), with
151
correlation coefficient being 0.45, which exceeds the 95% confidence level. The aforementioned
152
indicates the significant relationship between the variations in Eurasian SSD and the East Asian
153
summer precipitation.
154
However, there are also some discrepancies between the two precipitation patterns. The
155
negative regions over Inner Mongolia in the EOF3 pattern appear to extend more northward
156
toward BL and are more significant than those of the correlation pattern (Figure 3b). These
157
discrepancies illustrate that other factors may also have certain contributions, or their linkage may
158
derive from the combined effects of sea surface temperature (Zhang et al., 1999) and Arctic sea ice
159
(Wu et al., 2009a). For instance, in 1998, the East Asian summer precipitation has negative
160
anomalies larger than two standard deviations, whereas the SSD anomaly is insignificant (Figure 8
161
3c), which may be related to Tibetan Plateau thermal forcing and East Pacific sea surface
162
temperature anomalies (Chinese National Climate Center, 1998; Chen, 2001; Lau and Weng,
163
2001).
164
The above analysis reveals a statistical linkage between the Eurasian SSD and East Asian
165
summer precipitation. In the next section, we will illustrate the physical mechanisms responsible
166
for this linkage and associated causal relationship through exploring the possible dynamic and
167
thermodynamic mechanisms associated with Eurasian SSD and East Asian summer precipitation,
168
with a focus on surface thermal forcings.
169
3.2 Possible physical mechanisms
170
Previous studies suggest that Eurasian snowmelt may directly increase the soil moisture and affect
171
soil temperature over some regions of Siberia during warming seasons, thereby leading to the
172
reduction of surface air temperature and tropospheric diabatic heating (Cohen and Rind, 1991;
173
Matsuyama et al., 1998; Saito and Cohen, 2003; Zuo et al., 2011). Figure 4 displays the regressed
174
100-289 cm soil moisture anomalies during late spring (AMJ) and summer (JJA) onto the SSDI
175
index. Figure 4a illustrates that the soil moisture anomalies feature an out-of-phase structure over
176
the middle and high latitudes of the Eurasian continent during the late spring. In comparison, the
177
western SSD anomalies are accompanied by stronger negative in situ moisture anomalies over
178
EEWSP while the eastern SSD anomalies by positive soil moisture anomalies over regions around
179
BL. The east-west contrast in soil moisture anomalies during late spring is consistent with the
180
east-west contrast in the SSD anomaly pattern; thus, the soil moisture anomalies may be attributed
181
to the effect of local spring SSD. The colocation of SSD and soil moisture anomalies supports the
182
clear snow hydrological effect during the period 1980–2013, with excessive (less) SSD leading to
9
183
more (less) snowmelt and thus wetter (drier) soil. Notably, the regions north of 70ºN feature
184
notable negative soil moisture anomalies, which is not in accordance with the excessive SSD
185
during spring. In fact, these regions have the latest melt season due to their higher latitudes and
186
the longest snowy seasons. Correspondingly, the soil moisture anomalies may not correspond well
187
to the snow decrement during spring. In addition to the mid-high latitudes, the majority of regions
188
south of 35ºN show coherent soil moisture anomalies. Because of the long memory of surface
189
hydrologic conditions in Eurasia (Shinoda, 2001), the Eurasian soil moisture structure persists to
190
the subsequent summer. As shown in Figure 4b, the pattern of the soil moisture anomalies in
191
summer is quite similar to that in spring. In this study, we focused on the 100-289 cm layer soil
192
moisture; other shallower soil levels display the coherent features but with relatively weaker
193
signals (figures not shown), which may be attributed to the longer memory of soil moisture
194
anomalies in deeper layer.
195
The soil moisture in summer (Figure 4b) can affect the surface thermal balance. Zhang and Zuo
196
(2011) examined the impact of soil moisture over East Asia on the surface heat balance and
197
pointed out that a wetter (drier) soil corresponds to a lower (higher) surface air temperature. To
198
analyze the effect of soil moisture on surface thermal condition in summer, we investigate the
199
surface solar radiation flux, longwave radiation flux, sensible heat flux, latent heat flux, and net
200
heat flux in regression onto the SSDI index, and the results are shown in Figure 5. The solar
201
radiation at the surface is enhanced over EEWSP and weakened in the vicinity of BL (Figure 5a).
202
Meanwhile, the longwave radiation varies oppositely (Figure 5b). These may result from locally
203
changed cloud cover in the two regions (Zampieri et al. 2009). The changed surface energy
204
further results in more sensible heat but less latent heat in EEWSP, and less sensible heat and
205
latent heat in the surrounding of BL (Figures 5c and 5d). A clear picture comes out in the net heat 10
206
flux (Figure 5e), which illustrates a similar anomalous pattern to that of the solar and sensible heat
207
flux. The similarity of the anomalous net heat flux pattern to the distribution of the Eurasian SSD
208
in spring as well as the soil moisture anomalies in spring and summer suggests a change of the
209
local energy balance by the Eurasian SSD anomalies.
210
In Figure 6 we show the atmospheric thermal features over Eurasian continent associated with
211
the SSD. From Figure 6a it can be seen that in summer over the dryer area in EEWSP significant
212
higher surface temperatures is observed, whereas over the wetter area around BL significant lower
213
temperatures appears. These anomalous surface thermal forcings can warm or cool the lower
214
troposphere. As shown in Figure 6b, the near-surface thermal conditions associated with SSD
215
substantially change the local atmospheric thickness of the lower-middle troposphere between
216
500 hPa and 1000 hPa (Z500-Z1000), with a notably thickened atmosphere over EEWSP and
217
thinner atmosphere over regions west of BL. The west-east dipole pattern in atmospheric
218
thickness is quite similar to those in surface temperature and soil moisture. The meridional heat
219
flux anomalies over East Asia shown in Figure 5e form the north-south dipole structure in
220
temperature anomalies over East Asia. This dipole further strengthens the meridional temperature
221
gradient and atmospheric baroclinicity over the mid-latitude of East Asia at 700 hPa and
222
simultaneously decreases them to the north and south (Figure 6c). These results reflect the
223
thermal processes associated with anomalous Eurasian SSD, and illustrate the possible
224
thermodynamic influences on atmospheric circulation over East Asia in summer.
225
According to the thermal wind adjustment theorem, the intensified atmospheric baroclinicity
226
can accelerate the 200 hPa westerly winds over the mid-latitudes of the Eurasian continent,
227
leading to a strengthened subtropical jet stream over northern China and easterly wind anomalies
228
north and south of the subtropical jet stream (Figure 7a). The strengthened subtropical jet tends to 11
229
act as a connection between near-surface thermal forcings and the development of meridional
230
circulation anomalies, as suggested by Mu and Zhou (2015). Figure 7b illustrates that at the 500
231
hPa geopotential heights, an anomalous cyclonic pattern and anomalous anticyclonic pattern
232
emerge to the north and south, respectively, of the exit region of the East Asian subtropical jet
233
stream. In addition to this dipole pattern, an evident mid-latitude Eurasian wave train also prevails
234
over the regions from the West Siberia Plain to the Northwest Pacific, with a negative center
235
located over Mongolia and positive centers over EEWSP and the Northwest Pacific, respectively.
236
This Eurasian wave train can be interpreted as a manifestation of a single stationary Rossby wave
237
train propagation. Meanwhile, coexisted with the mid-latitude wave train, in the higher latitudes
238
there is a relatively weaker polar wave train with weakened heights emerged over regions around
239
western Siberian coast and strengthened heights over eastern Siberian coast. These results imply
240
that the anomalous SSD-induced upper-level westerly anomalies can provide favorable dynamic
241
conditions for the development of a north-south dipole pattern and Eurasian wave train pattern
242
over the lower-middle and higher latitudes of East Asia, which tend to act as the atmospheric
243
bridge linking the surface thermal forcings and East Asian summer climate.
244
In addition to the influences of the wave trains discussed above, the SSD and/or soil moisture
245
anomalies may also directly trigger the wave trains during late spring (AMJ) and summer. The
246
origins and development of the Eurasian wave trains are demonstrated in Figure 8, which presents
247
the 500 hPa stream function anomalies and associated wave activity flux (WAF, Takaya and
248
Nakamura, 2001) during late spring and summer. Over EEWSP, significant positive stream
249
function anomalies associated with reduced SSD and strengthened thermal heating occur in situ
250
during late spring. Meanwhile, there are relatively weak WAFs spreading from extratropical North
251
Atlantic northeastward to eastern Europe along the sub-polar belt. This wave train is then 12
252
markedly reinforced over eastern Europe, propagating northeastward to the West Siberia Plain,
253
and then splitting into two branches, with one turning to the higher latitudes and another
254
southeastward to East Asian (Figure 8a). This reinforced propagation pattern can be attributed to
255
the thermal forcing associated with SSD anomalies. In summer, similar and significant
256
mid-latitude WAFs emanated from eastern Europe and propagated directly eastward to BL and
257
further to East Asia, exhibiting predominantly zonally oriented. The horizontal propagation of
258
quasi-stationary wave train in summer (Figure 8b) resemble the Silk-Road and polar wave trains
259
as suggested by Orsolini et al. (2015) that have substantial impacts on precipitation over northeast
260
China. Compared with those during late spring, in summer the upper-stream wave train across the
261
extratropical North Atlantic is much weaker and the wave propagation is more zonally oriented.
262
This implies that the Eurasian wave trains are essentially triggered by SSD anomaly over EEWSP,
263
and the persistence of the wave propagation is mainly related to the consecutive surface thermal
264
forcings.
265
Some recent studies have emphasized the interference between the forced and climatological
266
wave trains (Smith et al. 2011, 2012; Peings and Magnusdottir 2014; Orsolini et al. 2015). These
267
studies suggest the importance of both forced and climatological wave trains in affecting the
268
atmospheric circulation. For instance, the atmospheric circulation is stronger when the forced
269
wave train is in phase with the climatological one, and vice versa. Figure 9 plots the forced
270
anomalous wave train regressed on the SSDI index along with the climatological wave train,
271
which are represented by 500 hPa geopotential heights. The Fourier decomposition is used to
272
decompose the stationary waves from wave 1 to wave 4 according to zonal wavenumber. In the
273
original field (Figure 9a), the climatological wave trains in summer exhibit the co-existing of the
274
polar and mid-latitude wave trains. The forced and climatological wave trains are in phase over 13
275
EEWSP and regions south of BL, and out of phase over regions north of BL and around Okhotsk
276
Sea, with a spatial correlation coefficient of 0.17. Over EEWSP the spatial correlation coefficient
277
is as high as 0.90, which denotes a constructive linear interference between the forced and
278
climatological wave train. Concerning wave 1 and wave 3 (Figures 9b and 9d), the forced wave
279
train is out of phase with the climatological one with the spatial correlation coefficients being
280
-0.53 and -0.44 for waves 1 and 3, respectively. These destructive interferences indicate that the
281
atmospheric circulation modified by wave 1 and wave 3 is opposite to the climatological wave
282
train. As for the wave 2 and wave 4 (Figures 9c and 9e), their spatial correlation coefficients with
283
the climatological wave train are 0.36 and 0.41, respectively, suggesting that the anomalous
284
forced waves should be linearly interfered with and enhance the climatological wave train. Thus,
285
the constructive linear interferences exhibited in wave 2 and 4 would reinforce the climatological
286
wave train and then prompt the influence on East Asian summer climate.
287
These anomalous circulation patterns provide favorable conditions for summer precipitation
288
over East Asia, particularly over China, with excessive precipitation over regions west of BL,
289
northeastern China and the Yellow River valley, and deficient precipitation over Inner Mongolia
290
and southern China (Figure 3a). These results indicate that the Eurasian SSD-induced diabatic
291
heating forcings may make prominent contributions to East Asian summer precipitation through
292
the formation and development of the Eurasian wave trains. However, these results depend largely
293
on observation analyses and cannot guarantee a causal relationship. Considering that the
294
numerical model simulations are essential to confirming the diagnostic results, the simulations are
295
addressed in the next section.
14
296
3.3 Model simulations
297
To elucidate the influence of anomalous Eurasian SSD on East Asian summer climate, we
298
performed a set of control experiments with CAM3.1 to derive different initial fields. The model
299
is spun up for 20 years and then integrates forward for 50 years. The data for December 1 of the
300
last 50 years are used as initial fields, and the 50 initial conditions are used for ensemble
301
simulation. Then, the control and sensitivity experiments, each with 50 members, are integrated
302
from the beginning of December to the end of August in the following year. In the sensitivity
303
experiments, we refer to Barnett et al. (1989) for the experimental design. To quantitatively mimic
304
the observed Eurasian SSD anomalies (as shown in Figure 2), we reduce the snowfall rate by 50%
305
over EEWSP and simultaneously increase the snowfall rate by 50% over regions around eastern
306
Europe, North Siberia Plain, and BL. Such modification is applied to the entire period of the
307
sensitivity simulations. Finally, the simulation results are obtained as the ensemble difference
308
between the sensitivity and the control experiments. To increase the robustness and reliability of
309
the numerical experiments, a 50-member ensemble mean is conducted to reduce the uncertainties
310
arising from differing initial conditions. Corresponding to the change in snowfall rate, other
311
hydrological and energy variables in CAM3.1 simulations change accordingly (see
312
Supplementary Fig. S3). In both late spring and summer, the deficient snowfall is accompanied by
313
negative anomalies of 100-289 cm layer soil moisture (Fig. S3a) and higher temperature over
314
EEWSP (Fig. S3c), while more snowfall by positive soil moisture anomalies (Fig. S3a) and lower
315
temperature (Fig. S3c) over regions around BL. These anomalous patterns are similar to the
316
Eurasian SWE and SSD anomalies over Eurasia (Figure 10).
317
To substantiate the conjecture of our designed experimental schemes, we primarily examined
318
the model-simulated response of spring SWE and SSD to the anomalous snowfall rate. As shown 15
319
in Figure 10a, two regions of large and significant negative anomalies are induced in the surface
320
SWE field during spring, one located over EEWSP and the other over the Tibetan Plateau, with
321
the maximum value occurring in the latter region and exceeding -8 mm. This reflects the in-phase
322
variation in SWE over the two regions and offsets the absence of observed lower-latitude SWE.
323
Similarly, three positive SWE anomalous regions are also triggered, with centers emerging over
324
eastern Europe, North Siberia Plain and BL. The maximum positive SWE center occurs in regions
325
west of BL, with the center value exceeding 8 mm. Corresponding to the simulated anomalous
326
SWE, the SSD reveals a similar pattern, namely, more SSD over eastern Europe and BL and less
327
SSD over EEWSP and the northern Tibetan Plateau (Figure 10b). However, there are some
328
differences between SWE and SSD. One is that the largest negative SSD center occurs over
329
EEWSP at approximately 60ºN, which differs from that of the SWE at approximately 70ºN.
330
Another is that the negative Tibetan Plateau SSD center is much weaker than that of the SWE and
331
even positive center appears, implying weak snowmelt and snow enhancement over the northern
332
and southern Tibetan Plateau, respectively, possibly due to the lower temperature over the Plateau.
333
Moreover, an insignificant negative SSD anomaly is simulated over the North Siberia Plain
334
maybe due to the persistent snowfall, even during early summer. These simulations largely
335
reproduce the observed anomalous SSD patterns (Figure 2a), indicating that our numerical
336
schemes of the snowfall rate changes are representative of the SSD anomalies in the CAM3.1.
337
To verify whether the anomalous SSD can induce the Eurasian wave train, we simulated the
338
physical processes that determine the circulation in response to snowfall rate changes. Figure 11
339
presents the simulated composite differences between sensitivity and control experiments in those
340
fields parallel to reanalysis (Figure 6-8). The surface air temperature field (Figure 11a) shows a
341
positive difference over eastern Europe and a negative difference over BL, which reproduce the 16
342
west-east dipole pattern in reanalysis. In the atmospheric thickness field (Figure 11b), two
343
enhanced anomalous thickness centers appear over eastern Europe and central China, whereas a
344
decreased thickness center emerges over regions north of BL. Meanwhile, in the 700 hPa
345
meridional temperature gradient field (Figure 11c), an anomalous baroclinicity center appear to
346
the south of BL. These patterns respond directly to the thermal forcing of Eurasian SSD. Then, a
347
strengthened subtropical jet stream (Figure 11d), an Eurasian mid-latitude wave train and a
348
north-south dipole structure over East Asia with a decreased height over BL and an enhanced
349
height over southern China (Figure 11e) are subsequently simulated. Figure 11f shows weak
350
WAFs propagating from extratropical North Atlantic eastward to East Asia, which is markedly
351
reinforced in eastern Europe, indicating the active role of SSD anomaly over EEWSP in triggering
352
a strong Rossby wave train.
353
The linear interference between the model forced and climatological wave trains are further
354
examined in Figure 12. In the total wave (Figure 12a), the forced and climatological waves are in
355
phase over EEWSP with the spatial correlation coefficient being 0.83 and out of phase over
356
regions north of BL. Concerning the decomposed waves, the forced wave is out of phase with the
357
climatological wave for the wave 1 and wave 3 with spatial correlation coefficients being -0.23
358
and -0.36, respectively, and in phase for the wave 2 and wave 4 with spatial correlation
359
coefficients being 0.23 and 0.39, respectively. These responses further confirm that the linear
360
interference between the forced and climatological waves in wave 2 and wave 4 intensifies the
361
influence on East Asian summer climate.
362
In general, these simulations reproduce the thermodynamic and dynamic responses of
363
atmospheric circulation to anomalous Eurasian SSD as identified in the analysis using observation
364
data, and demonstrate the previous argument that the anomalous Eurasian SSD can excite a 17
365
Eurasian wave train propagating from eastern Europe eastward to East Asia. However, there are
366
still some discrepancies in the simulated geopotential height response for the polar wave train
367
compared to that in reanalysis. For example, the total wave response (Figure 12a) is further north
368
than in the reanalysis (Fig 9a). Examining Figure 12d, it seems that the model has a stronger
369
forced response at high latitudes than the reanalysis (Fig 9d). In Fig 12e, the response is slightly
370
shifted north compared to reanalysis (Fig 9e). These discrepancies may be due to the weak SSD
371
forcing over high-latitudes in CAM3.1 or contributions from other factors such as North Atlantic
372
SST (Wu et al., 2011) and Arctic sea ice (Wu et al., 2013).
373
The associated summer precipitation anomalies are also shown in Figure 13. Positive
374
precipitation anomalies are identified over regions from BL eastward to northeastern China and
375
the lower Yellow River valley, and negative anomalies over western Inner Mongolia and most
376
regions of southern China. When compared with the observed anomalous precipitation patterns in
377
Figures 3a and 3b, opposite signs appear in regions from northern China to the middle Yellow
378
River valley. These discrepancies may arise from the deficiency of CAM3.1 in simulating
379
variations in East Asia summer precipitation, as described in Hack et al. (2006) and Wei et al.
380
(2011). Nevertheless, the numerical simulations demonstrate the significant role of Eurasian SSD
381
in influencing the East Asian summer precipitation through modulating atmospheric circulation
382
pattern.
383 384
4. Discussion on other factors’ effect
385
Although the close connection between Eurasian SSD and East Asian summer climate has been
386
examined with both observations and numerical model simulations in the above sections, other
387
atmospheric factors and lower boundary forcings, such as the ENSO (Zhang et al., 1996, 1999; 18
388
Ferranti and Molteni, 1999; Shaman and Tziperman, 2005) and Arctic sea ice (Wu et al., 2009b,
389
2009c; Zhang and Wu, 2011; Bader et al., 2011; Liu et al., 2012), can act as important drivers of
390
variations in both Eurasian snow cover and East Asian summer climate. In this section, we will
391
employ partial correlation analyses (Spiegel, 1988) to examine if the relationship between
392
Eurasian SSD and East Asian summer climate revealed in our present study is affected by the
393
ENSO and Arctic sea ice.
394
Figure 14 shows the partial correlation coefficient patterns between the SSDI index and East
395
Asian summer precipitation after the linear parts related to Nino3.4 index and SIC index are
396
removed, respectively. The Nino3.4 index is defined as the average of eastern equatorial Pacific
397
sea surface temperature anomalies in the area -5ºS–5ºN, 120º–170ºW. The SIC index refers to the
398
regional averaged SIC north of 60ºN. The partial correlation patterns feature similar structures as
399
that of the EOF3 pattern of East Asian summer precipitation (Figure 3a) and the correlation
400
coefficient pattern (Figure 3b), with anomalous meridional quadrupole precipitation centers
401
extending from northern to southern East Asia. To further examine their resemblance, we
402
perform a spatial correlation analysis to confirm the resemblances between the partial correlation
403
coefficient patterns and EOF3 pattern of East Asian summer precipitation. The spatial correlation
404
coefficients between the EOF3 pattern of East Asian summer precipitation and the partial
405
correlation coefficient patterns with ENSO and SIC removed are 0.65 and 0.68, respectively,
406
which are close to 0.67 between the correlation pattern (Figure 3a) and the EOF3 pattern of East
407
Asian summer precipitation (Figure 3b). Therefore, the relationship between Eurasian SSD and
408
East Asian summer precipitation reveal in our present study should be independent of ENSO and
409
Arctic sea ice concentration.
410
5. Conclusions
411
In the present study, we investigated the relationship between Eurasian spring snow decrement
412
(SSD) and East Asian summer precipitation and the related thermodynamic and dynamic 19
413
mechanisms using both observational data and the CAM3.1 model. The results show that the
414
second EOF mode of Eurasian SSD exhibits a west-east dipole pattern, with a negative center
415
located over eastern Europe and the West Siberia Plain (EEWSP) and a positive center over the
416
area around Baikal Lake (BL). This anomalous SSD pattern is significantly associated with the
417
third EOF mode of East Asian summer rainfall through triggering an anomalous mid-latitude
418
Eurasian wave train. The reduced SSD over EEWSP tends to decrease the local soil moisture from
419
spring to the following summer, thereby increase the surface heat flux and near-surface
420
temperatures. Similarly, the increase in SSD over BL is accompanied by anomalously low
421
near-surface temperatures. Changes in near-surface temperatures further intensify the meridional
422
temperature gradient and lower-level baroclinicity, leading to the acceleration of the upper-level
423
subtropical westerly jet stream. At the 500 hPa geopotential heights, an anomalous cyclone and
424
anomalous anticyclone emerge to the north and south, respectively, of the exit region of East
425
Asian subtropical jet stream. Meanwhile, the changed surface thermal conditions enhance the
426
local 1000–500 hPa thickness over EEWSP while decreasing it over BL. These factors both create
427
favorable physical conditions for the maintenance and enhancement of the anomalous
428
mid-latitude Eurasian wave train prevailing over the regions from eastern Europe eastward to the
429
Northwest Pacific. We further explored the origin of the Eurasian wave train and found that there
430
are zonally oriented WAFs spreading from eastern Europe eastward to East Asia. This finding
431
demonstrates the role of anomalous SSD in triggering the Eurasian wave train. These circulation
432
patterns ultimately significantly influence the precipitation over East Asia, especially over China,
433
with excessive precipitation over regions west of BL, northeastern China and the Yellow River
434
valley and deficient precipitation over Inner Mongolia and southern China. Therefore, Our study
435
confirms the significant role of Eurasian SSD in influencing East Asian summer precipitation. 20
436
Our model results demonstrate that the CAM3.1 can reproduce the positive summer
437
precipitation anomalies over most of northern East Asia and negative anomalies over southern
438
China. Corresponding to the imposed anomalous Eurasian SSD forcings in CAM3.1, positive
439
surface air temperature and atmospheric thickness responses over EEWSP and negative responses
440
over BL are primarily simulated. In the simulation the subtropical jet stream over East Asia is
441
strengthened, the Eurasian mid-latitude wave train is formed, and reinforced WAFs propagate
442
from eastern Europe to East Asia. These simulated dynamic and thermodynamic processes result
443
in according precipitation anomalies to occur over East Asia.
444 445
Acknowledgments
446
The authors are very grateful for the constructive comments from anonymous reviewers. This
447
research was funded by the National Key Research and Development Program of China (Grant
448
2016YFA0601501, 2016YFA0600602), the National Basic Research Program of China (Grant
449
2013CB430203), the Basic Scientific Research and Operation Foundation of the CAMS (Grant
450
2015Z001), the National Natural Science Foundation of China (NSFC, Grants 41661144017,
451
41205059, and 41505053), and the Public Sector (Meteorology) Research and Special Project of
452
China (Grant GYHY20120617).
21
453
References
454
Barnett, T. P., and Coauthors, 1989: The effect of Eurasian snow cover on regional and global
455
climate variations. J. Atmos. Sci., 46, 661–685.
456
Bamzai, A.S., and L. Marx, 2000: COLA AGCM simulation of the effect of anomalous spring
457
snow over Eurasia on the Indian summer monsoon. Quart. J. Roy. Meteor. Soc., 126, 2575–
458
2584.
459
Balsamo, G., and Coauthors, 2009: A revised hydrology for the ECMWF model: Verification
460
from field site to terrestrial water storage and impact in the Integrated Forecast System. J.
461
Hydrometeorol, 10, 623–643.
462
Bader, J., M. D. S. Mesquita, K. I. Hodges, N. Keenlyside, S. Osterhus, and M. Miles, 2011: A
463
review on Northern Hemisphere sea-ice, storminess and the North Atlantic Oscillation:
464
Observations and projected changes. Atmos. Res., 101, 809–834.
465
Chen, L. T., 2001: The role of the anomalous snow cover over the Qinghai-Xizang Plateau and
466
ENSO in the great floods of 1998 in the Changjiang River valley. Chin. J. Atmos. Sci., 2, 004.
467
(in Chinese)
468
Chen, X. F., and W. L. Song, 2000a: Circulation analysis of different influence of snow cover over
469
the Tibetan Plateau and Eurasia in winter on summertime droughts and floods of China. Atmos.
470
Sci., 24, 585–592.
471
Chen, X. F., and W. L. Song, 2000b: Analysis of relationship between snow cover on Eurasia and
472
Qinghai-Xizang Plateau in winter and summer rainfall in China and application to prediction.
473
Plateau Meteor., 19, 214–224. (in Chinese)
474 475
Chinese
National
Climate
Center, 1998: The
1998
Abnormality. Meteorology Press, 137pp. (in Chinese) 22
Flood
in
China
and
Climatic
476
Cohen, J., and D. Rind, 1991: The effect of snow cover on the climate. J. Climate, 4, 89–706.
477
Collins, W. D., C. M. Bitz, and M. L. Blackmon, 2006: The Community Climate System Model
478
version 3 (CCSM3). J. Climate, 19, 2122–2143.
479
Dash, S. K., G. P. Singh, M. S. Shekhar, and A. D. Vernekar, 2005: Response of the Indian
480
summer monsoon circulation and rainfall to seasonal snow depth anomaly over Eurasia. Clim.
481
Dyn., 24, 1–10.
482
Ferranti, L., and F. Molteni, 1999: Ensemble simulations of Eurasian snow–depth anomalies and
483
their influence on the summer Asian monsoon. Quart. J. Roy. Meteor. Soc., 125, 2597–2610.
484
Gong, G., D. Entekhabi, and J. Cohen, 2003: Modeled Northern Hemisphere winter climate
485 486 487 488 489
response to realistic Siberian snow anomalies. J. Climate, 16, 3917–3931. Gong, G., D. Entekhabi, J. Cohen, and D. Robinson, 2004: Sensitivity of atmospheric response to modeled snow anomaly characteristics. J. Geophys. Res., 109, D06107. Groisman, P. Y., and Coauthors, 1994: Changes of snow cover, temperature, and radiative heat balance over the Northern Hemisphere. J. Climate, 7, 1633–1656.
490
Hack, J. J., and Coauthors, 2006: Simulation of the Global Hydrological Cycle in the CCSM
491
Community Atmosphere Model Version 3 (CAM3): Mean Features. J. Climate, 19, 2199–2221.
492
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-year reanalysis project. Bull. Amer.
493 494 495
Meteor. Soc., 77, 437–471. Kripalani, R. H., B. J. Kim, J. H. Oh, and S. E. Moon, 2002: Relationship between Soviet snow and Korean rainfall. Int. J. Climatol., 22, 1313–1325.
496
Lau, K. M., and H. Y. Weng, 2001: Coherent modes of global SST and summer rainfall over
497
China: An assessment of the regional impacts of the 1997-98 El Niño. J. Climate, 14,
498
1294-1308. 23
499 500 501 502 503 504
Li, Z. K., B. Y. Wu, W. J. Zhu, 2009: Numerical simulation on effect of spring Eurasian snow cover on summer rainfall in China. Adv. Clim. Change, 5, 196–201. (in Chinese) Liu, L., R. H. Zhang, and Z. Y. Zuo, 2014: Intercomparison of spring soil moisture among multiple reanalysis data sets over eastern China. J Geophys. Res.: Atmos., 119, 54-64. Liu, J. P., J. A. Curry, H. J. Wang, M. R. Song, and R. M. Horton, 2012: Impact of declining Arctic sea ice on winter snowfall. Proc. Natl. Acad. Sci., 109, 4074–4079.
505
Matsuyama, H., and K. Masuda, 1998: Seasonal/interannual variations of soil moisture in the
506
former USSR and Its relationship to Indian summer monsoon rainfall. J. Climate, 11, 652–658.
507
Mu, S. N., and G. Q. Zhou, 2012: Mechanism for the correlation of winter fresh snow extent over
508
Northern Eurasia and summer climate anomalies in China: Anomalous seasonal transition of
509
land as a bond. Chin. J. Atmos. Sci., 36, 297–315. (in Chinese)
510
Mu, S. N., and G. Q. Zhou, 2015: Mechanism for the correlation of winter fresh snow extent over
511
Northern Eurasia and summer climate anomalies in China: The important role of Eurasian
512
winter fresh snow extent. Chin. J. Atmos. Sci., 39, 611–633. (in Chinese)
513 514
North, G. R., T. L. Bell, R. F. Cahalan, and F. J. Moeng, 1982: Sampling errors in the estimation of EOFs. Mon. Wea. Rev., 110, 699–706.
515
Orsolini, Y. J., L. Zhang, D. H. W. Peters, et al., 2015: Extreme precipitation events over north
516
China in August 2010 and their link to eastward-propagating wave-trains across Eurasia:
517
observations and monthly forecasting. Quart. J. Roy. Meteor. Soc., 141, 3097–3105.
518
Peings, Y., and G. Magnusdottir, 2013: Response of the wintertime Northern Hemisphere
519
atmospheric circulation to current and projected Arctic sea ice decline: A numerical study with
520
CAM5. J. Climate, 27, 244-264.
521
Rayner, N. A., and Coauthors, 2003: Global analyses of sea surface temperature, sea ice, and 24
522 523 524
night marine air temperature since the late nineteenth century. J. Geophys. Res., 108, 4407. Saito, K., and J. Cohen, 2003: The potential role of snow cover in forcing interannual variability of the major Northern Hemisphere mode. Geophys. Res. Lett., 30, 1302.
525
Shaman, J., and E. Tziperman, 2005: The Effect of ENSO on Tibetan Plateau Snow Depth: A
526
Stationary Wave Teleconnection Mechanism and Implications for the South Asian Monsoons. J.
527
Climate, 18, 2067–2079.
528 529 530
Shinoda, M., 2001: Climate memory of snow mass as soil moisture over central Eurasia. J. Geophys. Res., 106, 393–403. Smith, K. L., and P. J. Kushner, 2012: Linear interference and the initiation of extratropical
531
stratosphere-troposphere interactions.
532
2012JD017587.
533 534
J. Geophys.
Res.,
117, D13107, doi:10.1029/
Smith, K. L., P. J. Kushner, and J. Cohen, 2011: The role of linear interference in northern annular mode variability associated with Eurasian snow cover extent. J. Climate, 24, 6185–6202.
535
Solberg, R., J. Amlien, H. Koren, B. Wangensteen, K. Luojus, J. Pulliainen, M. Takala, J.
536
Lemmetyinen, T. Nagler, H. Rott, F. Muller, C. Derksen, S. Metsamaki, and K. Bottcher, 2009:
537
Global Snow Monitoring for Climate Research - Design Justification File. European Space
538
Agency Study Contract Report, 21703/08/I-EC Deliverable 1.7. 246 pages.
539 540
Spiegel, M. R., 1988: Theory and Problems of Statistics, Schaum’s Outline Series, McGraw-Hill, 5th Printing, 359 pp.
541
Takaya, K., and H. Nakamura, 2001: A formulation of a phase-independent wave-activity flux for
542
stationary and migratory quasigeostrophic eddies on a zonally varying basic flow. J. Atmos. Sci.,
543
58, 608–627.
544
Walsh, J. E., D. R. Tucek, and M. R. Peterson, 1982: Seasonal snow cover and short–term 25
545 546 547 548 549 550 551 552 553 554 555 556 557 558 559
climatic fluctuations over the United States. Mon. Wea. Rev., 110, 1474–1485. Wei, T., L. Wang, Dong W, and J. Zhang, 2011: A comparison of East Asian summer monsoon simulations from CAM3.1 with three dynamic cores. Theor. Appl. Climatol., 106, 295–306. Wu, B. Y., K. Yang, R. H. Zhang, 2009a: Eurasian snow cover variability and its association with summer rainfall in China. Adv. Atmos. Sci., 26, 31–44. Wu, B. Y., R. H. Zhang, B. Wang, and R. D'Arrigo, 2009b: On the association between spring Arctic sea ice concentration and Chinese summer rainfall. Geophys. Res. Lett., 36, L09501. Wu, B. Y., R. H. Zhang, and B. Wang, 2009c: On the association between spring arctic sea ice concentration and Chinese summer rainfall: A further study. Adv. Atmos. Sci., 26, 666–678. Wu, B., R. Zhang, R. D'Arrigo, and J. Su, 2013: On the relationship between winter sea ice and summer atmospheric circulation over Eurasia. J. Climate, 26, 5523-5536. Wu, R. G., and B. P. Kirtman, 2007: Observed relationship of spring and summer East Asian rainfall with winter and spring Eurasian snow. J. Climate, 20, 1285–1304. Wu, R., S. Yang, S. Liu, L. Sun, Y. Lian, and Z. Gao, 2011: Northeast China summer temperature and North Atlantic SST. J. Geophys. Res., 116, D16116, doi: 10.1029/2011JD015779.
560
Wu, R. G., G. Liu, and P. Zhao, 2014: Contrasting Eurasian spring and summer climate anomalies
561
associated with western and eastern Eurasian spring snow cover changes. J. Geophy. Res., 119,
562
7410–7424.
563
Yasunari, T., A. Kitoh, and T. Tokioka, 1991: Local and remote responses to excessive snow cover
564
mass over Eurasia appearing in the northern spring and summer climate–A study with the
565
MRI.GCM. J. Meteor. Soc. Japan, 69, 473–485.
566
Yim, S. Y., J. G. Jhun, R. Lu, and B. Wang, 2010: Two distinct patterns of spring Eurasian snow
567
cover anomaly and their impacts on the East Asian summer monsoon. J. Geophys. Res., 115, 26
568
D22113.
569
Zampieri, M., F. D'Andrea, R. Vautard, P. Ciais, N. D. Nobletducoudré, and P. Yiou, 2009: Hot
570
European summers and the role of soil moisture in the propagation of mediterranean drought.
571
J. Climate, 22, 4747–4758.
572 573 574 575
Zhang, R. H., A. Sumi and M. Kimoto, 1996: Impact of El Niño on the East Asia Monsoon: A diagnostic study of the '86/87 and '91/92 events. J. Meteor. Soc. Japan, 74, 49–62. Zhang, R. H., A. Sumi, M. Kimoto, 1999: A diagnostic study of the impact of El Nino on the precipitation in China. Adv. Atmos. Sci., 16, 229–241.
576
Zhang, R. H., B. Y. Wu, P. Zhao, and J. P. Han, 2008: The decadal shift of the summer climate in
577
the late 1980s over eastern China and its possible causes. Acta. Meteor. Sinica, 22, 435–445.
578
Zhang, R. H., and Z. Y. Zuo, 2011: Impact of spring soil moisture on surface energy balance and
579
summer monsoon circulation over East Asia and precipitation in East China. J. Climate, 24,
580
3309–3322.
581
Zhang, R. H., B. Y. Wu, J. P. Han, and Z. Y. Zuo, 2013: Effects on summer monsoon and rainfall
582
change over China due to Eurasian snow cover and ocean thermal conditions. In: Clim. Change
583
- Realities, Impacts over Ice Cap, Sea Level and Risks (Ed. B. R. Singh), InTech, Rijeka, 227–
584
250.
585
Zhang, R. N, and B. Y. Wu, 2011: The Northern Hemisphere Atmospheric Response to Spring
586
Arctic Sea Ice Anomalies in CAM3.0 Model. Chin. J. Atmos. Sci., 35, 847–862. (in Chinese)
587
Zhao, P., R. H. Zhang, J. P. Liu, X. Zhou, and J. He, 2007: Onset of southwesterly wind over
588 589 590
eastern China and associated atmospheric circulation and rainfall. Clim. Dyn., 28, 797–811. Zuo, Z. Y., and R. H. Zhang, 2007: The spring soil moisture and the summer rainfall in eastern China. Chin. Sci. Bull., 52, 3310–3312. 27
591 592
Zuo, Z. Y., and R. H. Zhang, 2009: Temporal and spatial features of the soil moisture in boreal spring in eastern China. Sci. China Series D, 52, 269–278.
593
Zuo, Z. Y., S. Yang, W. Q. Wang, A. Kumar, Y. Xue, and R. H. Zhang, 2011: Relationship between
594
anomalies of Eurasian snow and southern China rainfall in winter. Environ. Res. Lett., 6,
595
045402.
596
Zuo, Z. Y., R. H. Zhang, B. Y. Wu, and X. Y. Rong, 2012: Decadal variability in springtime snow
597
over Eurasia: Relation with circulation and possible influence on springtime rainfall over China.
598
Int. J. Climatol., 32, 1336–1345.
599 600
Zuo, Z. Y., and R. H. Zhang, 2016: Influence of soil moisture in eastern China on East Asian summer monsoon. Adv. Atmos. Sci., 33, 151–163.
28
601
Figure captions:
602
Figure 1. Climatological annual cycle of monthly snow water equivalent (SWE, light blue bars)
603
and snow decrement (SD, orange bars) over the Eurasian continent (35º–76ºN, 0º–
604
150ºE) during 1980–2013 (Units: mm). The data in June-September are absent.
605
Positive and negative SDs represent increase and decrease SWE in succeeding month
606
with respect to the previous one.
607
Figure 2. (a) Spatial pattern (shadings, units: mm) of the second EOF (EOF2) and (b) detrended
608
time series of the second principal component (PC2, black line) of SSD over the
609
Eurasian continent during 1980–2013. The blue line in (b) shows the SSDI index
610
defined as the difference of SSD between the rectangles with positive values (46º–
611
58ºN, 96º–138ºE) and negative values (54º–68ºN, 48º–84ºE) shown in (a). Purple and
612
black dotted areas in (a) denote values with statistical significance exceeding the 90%
613
and 95% confidence levels, respectively.
614
Figure 3. (a) Correlation coefficients between the SSDI index and the subsequent summer
615
precipitation, (b) spatial pattern of the third EOF mode of East Asian summer
616
precipitation, and (c) detrended time series of SSDI and summer precipitation (PC3).
617
R represents the correlation coefficient between the two indices. Purple and black
618
dotted areas denote values with statistical significance exceeding the 90% and 95%
619
confidence levels, respectively.
620
Figure 4. (a) Late spring (AMJ) and (b) summer (JJA) 100-289 cm soil moisture (SM) anomalies
621
(shadings; units: m3/m3) regressed onto the SSDI index. Purple and black dotted areas
622
denote values with statistical significance exceeding the 90% and 95% confidence
623
levels, respectively. The rectangles are same as those in Figure 2. 29
624
Figure 5. Summer (JJA) (a) surface solar radiation flux, (b) longwave radiation flux, (c) sensible
625
heat flux, (d) latent heat flux, and (e) net heat flux anomalies (shadings; units: W/m2)
626
regressed onto the SSD index. Purple and black dotted areas denote values with
627
statistical significance exceeding the 90% and 95% confidence levels, respectively.
628
The rectangles are same as those in Figure 2.
629
Figure 6. (a) Surface temperature (shadings; units: °C), (b) atmospheric thickness between 500
630
hPa and 1000 hPa (Z500-Z1000, units: gpm), and (c) 700 hPa meridional temperature
631
gradient (shadings; units: °C/grid) anomalies during summer regressed onto the SSDI
632
index. Purple and black dotted areas denote values with statistical significance
633
exceeding the 90% and 95% confidence levels, respectively. The rectangles are same
634
as those in Figure 2.
635
Figure 7. (a) 200 hPa zonal wind (shadings; units: m/s) and (b) 500 hPa geopotential height
636
(shadings; units: gpm) anomalies regressed onto the SSDI index. The climatological
637
zonal winds with speed greater than 20 m/s are enclosed by the black lines in (a).
638
Purple and black dotted areas denote values with statistical significance exceeding the
639
90% and 95% confidence levels, respectively. The rectangles are same as those in
640
Figure 2.
641
Figure 8. 500 hPa stream function (shading; units: m2/s) and the associated wave activity flux
642
(vector; units: m2/s2) anomalies regressed onto the SSDI index during (a) late spring
643
(AMJ) and (b) summer (JJA). The rectangles are same as those in Figure 2.
644
Figure 9. SSD forced (contours) and climatological (shading) stationary waves represented by
645
summer 500 hPa geopotential heights for (a) total wave and (b–e) zonal wavenumbers
646
1–4. The stationary waves are obtained by removing the zonal mean of the 30
647
geopotential heights. R denotes the spatial coefficients between the SSD forced and
648
climatological waves. Contour intervals are 10 and 4 gpm for the climatological and
649
SSD forced total waves, respectively, and 5 and 2 gpm for climatological and forced
650
waves of zonal wavenumbers 1–4, respectively. The rectangles are same as those in
651
Figure 2.
652
Figure 10. CAM3.1 simulations of (a) spring SWE and (b) SSD (shadings; units: mm) anomalies
653
in response to the anomalous snowfall rate over the Eurasian continent. The thin and
654
thick black lines denote values with statistical significance exceeding the 90% and 95%
655
confidence levels, respectively.
656
Figure 11. CAM3.1 simulations of summer (a) surface air temperature (shadings, units: °C), (b)
657
atmospheric thickness (Z500-Z1000), (c) 700 hPa meridional temperature gradient
658
(shadings; units: °C/grid), (d) 200 hPa zonal wind (shadings, units: m/s), (e) 500 hPa
659
geopotential height (shadings, units: gpm) anomalies and (f) 500 hPa stream function
660
(shadings, units: m2/s) and associated wave activity flux (vectors, units: m/s)
661
anomalies. The thin and thick black lines denote values with statistical significance
662
exceeding the 90% and 95% confidence levels, respectively. The rectangles are same
663
as those in Figure 2.
664
Figure 12. Same as Figure 9 but for CAM3.1 simulations.
665
Figure 13. CAM3.1 simulations of summer precipitation anomalies (shadings, units: mm). The
666
thin and thick black lines denote values with statistical significance exceeding the 90%
667
and 95% confidence levels, respectively.
668
Figure 14. Partial correlation coefficients of SSDI index with East Asian summer precipitation
669
when the linear parts related to the (a) Nino3.4 and (b) Arctic sea ice concentration 31
670
(SIC) are removed. The Nino3.4 and Arctic SIC are detrended before calculating the
671
correlation coefficients. Purple and black dotted areas denote values with statistical
672
significance exceeding the 90% and 95% confidence levels, respectively.
32
673 674
Figure 1. Climatological annual cycle of monthly snow water equivalent (SWE, light blue bars)
675
and snow decrement (SD, orange bars) over the Eurasian continent (35º–76ºN, 0º–150ºE) during
676
1980–2013 (Units: mm). The data in June-September are absent. Positive and negative SDs
677
represent increase and decrease SWE in succeeding month with respect to the previous one.
678
33
679 680
Figure 2. (a) Spatial pattern (shadings, units: mm) of the second EOF (EOF2) and (b) detrended
681
time series of the second principal component (PC2, black line) of SSD over the Eurasian
682
continent during 1980–2013. The blue line in (b) shows the SSDI index defined as the difference
683
of SSD between the rectangles with positive values (46º–58ºN, 96º–138ºE) and negative values
684
(54º–68ºN, 48º–84ºE) shown in (a). Purple and black dotted areas in (a) denote values with
685
statistical significance exceeding the 90% and 95% confidence levels, respectively.
686
34
687 688
Figure 3. (a) Correlation coefficients between the SSDI index and the subsequent summer
689
precipitation, (b) spatial pattern of the third EOF mode of East Asian summer precipitation, and (c)
690
detrended time series of SSDI and summer precipitation (PC3). R represents the correlation
691
coefficient between the two indices. Purple and black dotted areas denote values with statistical
692
significance exceeding the 90% and 95% confidence levels, respectively.
693
35
694 695
Figure 4. (a) Late spring (AMJ) and (b) summer (JJA) 100-289 cm soil moisture (SM) anomalies
696
(shadings; units: m3/m3) regressed onto the SSDI index. Purple and black dotted areas denote
697
values with statistical significance exceeding the 90% and 95% confidence levels, respectively.
698
The rectangles are same as those in Figure 2.
699
36
700 701
Figure 5. Summer (JJA) (a) surface solar radiation flux, (b) longwave radiation flux, (c) sensible
702
heat flux, (d) latent heat flux, and (e) net heat flux anomalies (shadings; units: W/m2) regressed
703
onto the SSD index. Purple and black dotted areas denote values with statistical significance
704
exceeding the 90% and 95% confidence levels, respectively. The rectangles are same as those in
705
Figure 2.
706
37
707 708
Figure 6. (a) Surface temperature (shadings; units: °C), (b) atmospheric thickness between 500
709
hPa and 1000 hPa (Z500-Z1000, units: gpm), and (c) 700 hPa meridional temperature gradient
710
(shadings; units: °C/grid) anomalies during summer regressed onto the SSDI index. Purple and
711
black dotted areas denote values with statistical significance exceeding the 90% and 95%
712
confidence levels, respectively. The rectangles are same as those in Figure 2.
38
713 714
Figure 7. (a) 200 hPa zonal wind (shadings; units: m/s) and (b) 500 hPa geopotential height
715
(shadings; units: gpm) anomalies regressed onto the SSDI index. The climatological zonal winds
716
with speed greater than 20 m/s are enclosed by the black lines in (a). Purple and black dotted
717
areas denote values with statistical significance exceeding the 90% and 95% confidence levels,
718
respectively. The rectangles are same as those in Figure 2.
719
39
720 721
Figure 8. 500 hPa stream function (shading; units: m2/s) and the associated wave activity flux
722
(vector; units: m2/s2) anomalies regressed onto the SSDI index during (a) late spring (AMJ) and (b)
723
summer (JJA). The rectangles are same as those in Figure 2.
724
40
725 726
Figure 9. SSD forced (contours) and climatological (shading) stationary waves represented by
727
summer 500 hPa geopotential heights for (a) total wave and (b–e) zonal wavenumbers 1–4. The
728
stationary waves are obtained by removing the zonal mean of the geopotential heights. R denotes
729
the spatial coefficients between the SSD forced and climatological waves. Contour intervals are
730
10 and 4 gpm for the climatological and SSD forced total waves, respectively, and 5 and 2 gpm
731
for climatological and forced waves of zonal wavenumbers 1–4, respectively. The rectangles are
732
same as those in Figure 2.
733
41
734
735 736
Figure 10. CAM3.1 simulations of (a) spring SWE and (b) SSD (shadings; units: mm) anomalies
737
in response to the anomalous snowfall rate over the Eurasian continent. The thin and thick black
738
lines denote values with statistical significance exceeding the 90% and 95% confidence levels,
739
respectively.
740
42
741 742
Figure 11. CAM3.1 simulations of summer (a) surface air temperature (shadings, units: °C), (b)
743
atmospheric thickness (Z500-Z1000) (shadings, units: gpm), (c) 700 hPa meridional temperature
744
gradient (shadings; units: °C/grid), (d) 200 hPa zonal wind (shadings, units: m/s), (e) 500 hPa
745
geopotential height (shadings, units: gpm), and (f) 500 hPa stream function (shadings, units: m2/s)
746
and associated wave activity flux (vectors, units: m/s) anomalies. The thin and thick black lines
747
denote values with statistical significance exceeding the 90% and 95% confidence levels,
748
respectively. The rectangles are same as those in Figure 2.
749
43
750 751
Figure 12. Same as Figure 9 but for CAM3.1 simulations.
752
44
753 754
Figure 13. CAM3.1 simulations of summer precipitation anomalies (shadings, units: mm). The
755
thin and thick black lines denote values with statistical significance exceeding the 90% and 95%
756
confidence levels, respectively.
757
45
758 759
Figure 14. Partial correlation coefficients of SSDI index with East Asian summer precipitation
760
when the linear parts related to the (a) Nino3.4 and (b) Arctic sea ice concentration (SIC) are
761
removed. The Nino3.4 and Arctic SIC are detrended before calculating the correlation coefficients.
762
Purple and black dotted areas denote values with statistical significance exceeding the 90% and
763
95% confidence levels, respectively.
46