1 SEPTEMBER 2016
SONG ET AL.
6363
Coupling Modes of Climatological Intraseasonal Oscillation in the East Asian Summer Monsoon ZEHAO SONG, CONGWEN ZHU, JINGZHI SU, AND BOQI LIU State Key Laboratory of Severe Weather, and Institute of Climate System, Chinese Academy of Meteorological Sciences, Beijing, China (Manuscript received 7 November 2015, in final form 1 June 2016) ABSTRACT The present study used harmonic and multivariate empirical orthogonal function (MV-EOF) analyses to identify the existence of climatological intraseasonal oscillation (CISO) in the diabatic heating, precipitation, and circulation of the East Asian summer monsoon (EASM). The strongest CISO signals are found in the north of the western North Pacific, possibly because of the horizontal gradient of diabatic heating induced by the seasonal land–sea thermal contrast. Further, the phase relationship between the diabatic heating components maintains the EASM CISO. The first two coupling modes of EASM CISO in the circulation are robust during May through August, with a period of 40–80 days, and exhibit phase locking to the stepwise establishment of the EASM, which reveals the coaction of the Mongolian cyclone (MC) around Lake Baikal at 850 hPa, the western North Pacific subtropical high (WNPSH) at 500 hPa, and the South Asian high (SAH) over the Tibetan Plateau (TP) at 200 hPa. The first mode shows that the jointly enhanced MC, WNPSH, and SAH correspond to a tripole rainfall anomaly with strong mei-yu and baiu fronts over East Asia. The second mode, however, indicates the eastward and northwestward propagation of MC and WNPSH, respectively, with suppressed SAH, as well as a dipole rainfall anomaly over East Asia. Both the observations and numerical simulation verify the importance of daily diabatic heating and SST in maintaining the CISO modes over the WNP, where the condensation heating related to atmospheric forcing determines the local intraseasonal air–sea interaction.
1. Introduction The East Asian summer monsoon (EASM) includes the tropical South China Sea (SCS) and the western North Pacific (WNP), along with the subtropical summer monsoon north of 208N over the East Asia–Pacific region (Lau and Li 1984; Zhu et al. 1986; Tao and Chen 1987; Wang and LinHo 2002; He et al. 2007; Zhu et al. 2011; Wang et al. 2012; Chi et al. 2015). The rain belts of the EASM with stronger low-level southwesterly winds span from the subtropical WNP to the higher latitudes of eastern Russia and become the most northern component in the global monsoon system (Wang et al. 2012). A number of studies have suggested that the WNP subtropical high (WNPSH) in the mid-to-lower troposphere (Figs. 1a,b) and the upper-level South Asian high (SAH)
Corresponding author address: Dr. Congwen Zhu, Institute of Climate System, Chinese Academy of Meteorological Sciences, Zhongguancun South St. 46, Beijing 100081, China. E-mail:
[email protected] DOI: 10.1175/JCLI-D-15-0794.1 Ó 2016 American Meteorological Society
over the southern Tibetan Plateau (Fig. 1c) are the major circulations of the EASM to regulate its rainfall. During the boreal summer, the enhanced WNPSH transports more warm moisture by strengthening the southwesterly winds, thus resulting in heavier summer rainfall in the mid-to-lower reaches of the Yangtze River in China (i.e., mei-yu fronts), in conjunction with the stronger ascent induced by the strengthened SAH (Tao and Zhu 1964; Ren et al. 2007; Lu et al. 2014). Stronger mei-yu fronts are also closely associated with enhanced rainstorm frequency and intensity due to the cold and warm airmass interaction between the tropical and mid-to-high latitudes (Tao 1980; Lau and Li 1984; Chen et al. 2001; Zhu et al. 2003; Lu et al. 2014), which is associated with the Mongolian cyclone (MC) and cold vortex over northern China (Li and Wang 2003; Miao et al. 2006a,b; He et al. 2007). In climatology, the EASM circulation and rainfall begin in early April, peak in August, and prevail over the East Asia–Pacific region until finally withdrawing in October (Zhu et al. 2005; Zhao et al. 2007; He et al. 2007;
6364
JOURNAL OF CLIMATE
VOLUME 29
FIG. 1. (a) Spatial pattern of the summer winds at 850 hPa (vector; m s–1) and precipitation (shading; mm day–1), averaged from May to August (MJJA); MJJA-averaged winds at (b) 500 and (c) 200 hPa (vector; m s 21 ); (d) seasonal variation of precipitation (contour) and its AC (shading) component along 1158E.
Zhu et al. 2011; Chi et al. 2015). The EASM rain belts exhibit a stepwise and subsequent northward jump with the seasonal migration of the EASM (Fig. 1d), thus resulting in the peak heavy rainfall in southern China, the Yangtze River, and northern China in early April, midJune, and late July, respectively (He et al. 2008; Tian and Yasunari 1998; Zhu et al. 2011). The timing of each enhanced rainfall fluctuation indicates the subtropical EASM, SCS, and tropical WNP monsoon onset, and the rapid northward shift of the rainband marking the start of the summer rainy season in southern China, mei-yu front, and northern China, respectively (e.g., Lau et al. 1988; Zhu et al. 2011). The seasonal progression of the EASM relative to its annual cycle exhibits subseasonal variation, which manifests as seasonal fluctuations of the EASM rainbands during May–August (Fig. 1d). This subseasonal variation of the EASM strongly modulates the waxing and waning of summer monsoon rainfall (e.g., Lau and Li 1984; Chen and Jin 1984; Lau et al. 1988; Wang and Rui 1990; Wang and Xu 1997; He et al. 2015; Chi et al. 2015). Wang and Xu (1997), by using climatological pentad mean outgoing longwave radiation (OLR) and the ECMWF winds have revealed that the subseasonal
variation in the Asian summer monsoon is directly associated with the activity of the climatological intraseasonal oscillation (CISO). They focused on the time series of OLR at 158N, 1408E and have confirmed the statistical significance of the CISO signal relative to the annual cycle in the WNP domain. Their results suggest that the summer monsoon CISO results from phase locking between the transient ISO and the annual cycle, which presents a dynamically coherent structure between enhanced convection and lowerlevel convergent (upper-level divergent) cyclonic (anticyclonic) circulation. The extreme phases of the CISO characterize monsoon singularities (i.e., monsoon events occurring on a fixed pentad with usual regularity; Wang and Xu 1997), thus implying the potential ability of CISO to predict the subseasonal variation of the EASM. Although the CISO of the EASM has been observed since the 1990s, its existence and detailed structure remain uncertain. Because the general circulation, particularly in the tropics, is driven by atmospheric diabatic heating (Yanai et al. 1976; Hantel and Baader 1978), which receives feedback from the circulation through atmospheric instabilities (Lau and Peng 1987; Emanuel
1 SEPTEMBER 2016
6365
SONG ET AL.
et al. 1994; Liu et al. 2004; Hazra and Krishnamurthy 2015), it has been argued that the annual and diurnal variations of radiative forcing with shortwave radiative heating and longwave radiative cooling produce the most significant climatological modes of the atmospheric circulation at the same time scale. However, it remains unclear whether the shortwave radiative heating and the longwave cooling intrinsically have an intraseasonal component. If so, can the CISO be treated as the atmospheric response to the intraseasonal variation of radiative heating and cooling? Moreover, because the EASM-related circulation exhibits both internal vertical and meridional coupling, with evident seasonal land–sea thermal contrast over the EASM, it is interesting to further investigate the coupling modes of the EASM CISO, their effects on the subseasonal variations of summer rainfall, and the possible impact from local intraseasonal air–sea interaction. The present study will first discuss the existence of CISO in the EASM by investigating the relationship among atmosphere diabatic heating, precipitation, and large-scale monsoon circulation over the East Asia– Pacific region. Subsequently, the coupling modes of EASM CISO and their possible impacts on the subseasonal rainfall activity will be further described. In addition, we will explore the possible regulation of air– sea interaction by the EASM CISO. The paper is organized as follows. The data and methods used are described in section 2. The existence of the CISO EASM modes is presented in section 3. Section 4 demonstrates the coupling modes of the EASM CISO in circulation, as revealed by the multivariate empirical orthogonal function (MV-EOF) analysis, and their impacts on the activity of subseasonal rainfall. The role of subseasonal air–sea interaction over the WNP in maintaining the EASM CISO is explored in section 5. Finally, a summary and discussion are provided in section 6.
2. Data and methods The data utilized in this study include the atmospheric components for the period 1981 to 2010 from the National Centers for Environmental Prediction–U.S. Department of Energy (NCEP–DOE) AMIP-II reanalysis (R-2; Kanamitsu et al. 2002), along with the pentadaveraged precipitation from the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) (Xie and Arkin 1997), with a horizontal resolution of 2.58 3 2.58. The daily NOAA high-resolution SST data (Reynolds et al. 2007) were provided by the NOAA/OAR/ESRL Physical Sciences Division (PSD), Boulder, Colorado, from their website (http://www.esrl. noaa.gov/psd/).
The vertically integrated atmospheric heat source Q1 (Hantel and Baader 1978; Yanai et al. 1992; Yanai and Tomita 1998; Liu et al. 2004) was calculated using the R-2 data as follows: Q1 5 SH 1 CO 1 RC,
(1)
where SH, CO, and RC represent the surface sensible heating, condensation heating released by convection, and net radiative heating of the air column, respectively. Here, SH was extracted from the NCEP–DOE reanalysis, and CO and RC were calculated using CO 5 Pr 3 Lw 3 r
and
RC 5 R‘ 2 R0 5 (SY‘
2 S[‘ ) 2 (SY0
(2) 2 S[0 )
2 (F0Y 2 F0[ ) 2 F‘ .
(3)
In Eq. (2), Pr, Lw 5 2.5 3 10 6 J Kg 21 , and r 5 1.0 3 10 3 kg m 23 are the precipitation amount, condensation heat coefficient, and water density, respectively. In Eq. (3), R‘ and R0 are the net radiation values at the top of atmosphere (TOA) and the land surface, respectively. The variables S and F denote the shortwave (SW) and longwave (LW) radiation fluxes, their subscripts ‘ and 0 denote the TOA and the ground surface, and the superscripts Y and [ represent downward and upward transport, respectively. The climatology was defined as the 30-yr arithmetic average value from 1981 to 2010, and the CISO was defined as the 20–90-day harmonic components relative to the annual cycle. Following previous studies (Wang and Xu 1997; Wheeler and Hendon 2004), we applied harmonic analysis to extract the CISO component for each grid as follows: A(t) 5 A 1 SVC(t) 1 ISVC(t) 1 R(t) , 3 2pm t , and SVC(t) 5 å Am exp i T0 m51 18 2pm ISVC(t) 5 å Am exp i t , T0 m54
(4) (5) (6)
where A(t) and A indicate the climatological daily mean and annual average for each element, respectively. The seasonal variation component SVC(t) was defined as the sum of the first three harmonics (m 5 1, . . ., 3) and corresponds to the annual and semiannual components. The term ISVC(t) represents the sum of harmonics 4 through 18 and has a period of 20–90 days. The residual in the harmonic analysis R(t) was for the transient components with periods less than 20 days. Here, we focused on the CISO in the Asia–Pacific region of 08–608N, 608–1608E during boreal summer (April–October). To identify the dominant modes of
6366
JOURNAL OF CLIMATE
VOLUME 29
FIG. 2. (a) The first harmonic-indicated distribution of the climatological annual cycle mode of diabatic heating, represented by its phase (vector) and amplitude (shading; W m22); (b) distribution of averaged Q1 during May– August (contour; W m22) and the standard deviation (shading; W m22) of CISO; (c) seasonal variation of the Q1 AC between 1158 (shading) and 1408E (contour); and (d) the difference of the AC (shading) and CISO (contour) in terms of Q1 between 1158 and 1408E (1158E minus 1408E). The phase of the cycle in (a) is shown as a 12-month clock, with a northward arrow indicating maximum heating in July. The arrows rotate clockwise with eastward, southward, and westward arrows indicating October, January, and April, respectively, and the arrow length defines the amplitude of the cycle. The contour interval of the Q1 AC is 50 W m22, and the positive (solid) and negative (dashed) green lines of the CISO interval in (d) is 30 W m22, where the amplitude below 30 W m22 is omitted.
vertical and meridional circulation interaction in the EASM, we used MV-EOF analysis (Wang 1992; Wang et al. 2008; He et al. 2015) on the wind fields at various levels (i.e., 850, 500, and 200 hPa). In addition, nondimensional eigenvectors are used to standardize the combined meteorological fields.
3. Existence of CISO in EASM a. CISO in diabatic heating The annual modes of Q1 indicated by the first harmonic component can be considered to be the atmospheric response to the annual solar radiation forcing. However, the phase peak of the annual cycle (AC) of Q1 is different from region to region, even at the same latitude, and it is very distinct between the land and ocean, owing to their different thermal capacities (Fig. 2a). For instance, to the
east of 1058E in the East Asia–Pacific region, Q1 AC over the continent reaches its peak in April and May at approximately 258 and 408N in southern and northern China, respectively. The peak maximum of the Q1 AC to the south of 308N over Asia generally appears before July, but over the Indian and Pacific Oceans, it occurs after July in August, with a 1–2-month lag. The positive amplitudes of the Q1 AC are centered to the south of 208N, with three branches over the Indian subcontinent, where it adjoins oceanic regions, Indo-China, and the eastern Philippine Sea (Fig. 2a). During May–August, the spatial pattern of the climatological Q1 is similar to that of the Q1 AC (Figs. 2a,b). Specifically, to the west of 1058E, the maximum standard deviation (STD) of CISO, located over the tropical Indian Ocean, is accompanied by the climatological Q1 center, whereas in the western Pacific, the CISO STD is centered over the SCS and the WNP regions, without the maximum of the climatological Q1.
1 SEPTEMBER 2016
SONG ET AL.
6367
FIG. 3. Seasonal variation of (a) AC and (b) AC-removed residuals [sum of CISO and residual in Eq. (3)] of Q1 (W m22) averaged over the WNP region (208–408N, 1158–1408E). The SW and LW radiation is multiplied by a factor of 5 in (a) and (b). (c) Lead–lag correlation coefficients between Q1 and its components without the annual cycle from May to October; (d) as in (c), but for the CISO results.
Moreover, the amplitude of CISO depends on the horizontal Q1 gradient during the seasonal match of the EASM. Figure 2c shows that the seasonal evolution of the Q1 AC between 1158 and 1408E is asymmetric in the meridional direction. The maximum of the Q1 AC at 1158E is located to the north of 208N during March– August, but it appears over the tropical ocean at 1408E to the south of 208N in May–December (Fig. 2c). Thus, the negative meridional gradient of the Q1 AC (›Q1/›y , 0) becomes evident in May–August over East Asia and the ocean to its east. Meanwhile, the negative zonal gradient of the Q1 AC (2›Q1/›x . 0) exists over the East Asia–Pacific region north of 208N from March to August, whereas the negative counterparts appear over the tropical western Pacific latitudes during all seasons (Fig. 2d). Therefore, in May–August, a negative vorticity source is formed to the north and east of the Q1 AC maximum, giving rise to the anticyclonic vorticity in the air column over East Asia and the surrounding ocean to the north of 208N (Liu et al. 2001). Subsequently, the air column contracts to enhance the local ascending and the Q1 AC, as required by the geostrophic relationship and the atmospheric continuous equation. Figures 3c,d also indicate that the CISO in Q1 is active in May–August over East Asia, with enhanced ascending and Q1 AC, where the CISO STD is centered (Fig. 2b). This result suggests that the existence of CISO depends partly on the distinct seasonal evolution of the Q1 AC between
East Asia and its surrounding oceans due to the land–sea thermal contrast and that the CISO is closely associated with the vertical motion over East Asia. In addition to the land–sea thermal contrast, the CISO in Q1 of the EASM is maintained by the phase relationship between the Q1 components. We applied a direct method to estimate each component separately (Newell et al. 1969, 1975). Our results suggest that the subseasonal variation of Q1 primarily derives from the AC-removed residuals of CO heating during May– October. The large STDs of Q1 and CO are equivalent over the WNP (208–408N, 1158–1408E) from May to August (figures not shown), but the contributions of SH, SW, and LW are less than 10% of the total. Similarly, we decomposed each component of Q1 into AC, CISO, and the residual terms by applying Eqs. (4)–(6). As shown in Fig. 3a, the AC of SW heating and LW cooling are evident with their peaks occurring in mid-June and early August, respectively. The maximum of SH occurs in the boreal winter, whereas the CO heating reaches its maximum in August. There is a lead–lag relationship between the AC of the Q1 components. The ACremoved variation of CO exhibits a larger fluctuation than do those of SW and LW, particularly during May– October (Fig. 3b). The subseasonal variation of SW and LW follows that of SH and CO, which is quite different from the case in AC mode (Fig. 3a). Specifically, the SH varies initially, followed by the CO variation;
6368
JOURNAL OF CLIMATE
VOLUME 29
FIG. 4. (a) Time series of pentad rainfall and its AC component at 308N, 1158E (mm day21); (b) the CISO and AC-removed residuals (mm day21); (c) the CISO and AC-removed residuals at 408N, 1158E (mm day21); (d) seasonal variation of EASMI (refer to the text for definition) and its AC component; (e) time series of CISO and the AC-removed residual of the EASM; and (f) wavelet analysis (Torrence and Compo 1998) of the residual, where the shading indicates regions of .95% confidence for a red-noise process, and the thick dashed line denotes the cone of influence, outside of which the edge effects become important.
subsequently, changes in SW and LW occur over the EASM region at a subseasonal time scale (Fig. 3c). Such a lead–lag relationship is better organized in the correlations of the CISO (Fig. 3d). In physics, the increase in SH is able to enhance the atmospheric convective instability, and the monsoon convection releases CO heating, thus resulting in increased cloud, which in turn modulates the SW and LW in situ. As a result, the subseasonal variation of SW and LW lags behind the CO heating, following the SH oscillation. Wavelet analysis of AC-removed Q1 and its components shows a significant 40–80-day oscillation during May–October (figures not shown). This periodicity is prominent in the SW, LW, and CO, but it is weak in the SH. The 40–80-day oscillation is consistent with the dominant periodicity of the CISO components, particularly the CO heating; this finding implies that the significance of CISO in Q1 can be mainly ascribed to the CO heating over the WNP and EASM region during the boreal summer. As a result, the evident CISO in Q1 is maintained by the lead–lag phase
relationship between the Q1 components, in which the feedback between SH and CO variations is dominant. In addition, the variations of SW heating and LW cooling are possible responses to the CO heating, instead of the origin of CISO in Q1.
b. CISO of EASM rainfall and circulation The CISO over the EASM region not only exists in Q1 but also is embedded in the monsoon rainfall and circulation. The summer rainfall related to the mei-yu front in the mid-to-low reaches of the Yangtze River indicates the remarkable activity of the EASM, where the rainfall at 308N, 1158E exhibits an outstanding subseasonal fluctuation relative to its AC mode (Fig. 4a). The time series of the AC-removed residual exhibits an intraseasonal oscillation that coincides with its CISO component during May–August; its positive and negative phases indicate the activity and cessation of EASM rainfall, respectively (Fig. 4b). A similar feature is also presented in the subseasonal variation of rainfall at
1 SEPTEMBER 2016
SONG ET AL.
408N, though with a smaller amplitude (Fig. 4c). To verify the existence of CISO in the EASM circulation, we adopted an EASM index (EASMI), which is defined on the basis of land–sea thermal contrasts by the composite of the normalized standard difference of sea level pressure (SLP; SLP at 1608E minus SLP at 1108E) and the vertical zonal wind shear (U850hPa 2 U200hPa) in the region of 08–108N, 1008–1308E (Zhu et al. 2005). The positive and negative phases of the AC in EASMI indicate the seasonal alternation between the summer and winter monsoon and the reversal of the land–sea thermal contrast over the East Asia–Pacific region. The seasonal variation of EASMI presents a stepwise enhancement with subseasonal fluctuation relative to its AC component (Fig. 4d). The EASMI without AC exhibits an intraseasonal oscillation with a significant 40– 80-day period during May–August, which is consistent with its CISO component (Figs. 4e,f). Consequently, the CISO does exist in the EASM and manifests as the residual between the daily climatology and its AC, with a significant 40–80-day period from May to August. As a result, the CISO of the EASM is closely associated with the subseasonal variation of CO heating ascribed to the zonal and meridional land–sea thermal contrast in the boreal summer monsoon season, and it is maintained via a lead–lag relationship between the Q1 components.
4. Coupling modes of EASM CISO and its impact on rainfall a. Coupling CISO modes in EASM The significant ISO in the diabatic heating, rainfall, and EASMI implies the possible existence of CISO modes in the vertical circulation over the East Asia– Pacific region. As shown in Fig. 1, the EASM rainfall is closely associated with the vertical coupling of circulation at different levels (i.e., the coupling modes). Because the CISO is significant in EASM rainfall, it is speculated that the coupling modes should also play an important role in the EASM CISO. To reveal the coupling modes of CISO in EASM circulation, we applied the MV-EOF analysis to the standard normalized components of CISO in the wind fields at 850, 500, and 200 hPa from April to October. Figure 5 shows the first two leading modes of the MV-EOFs, which account for 17.1% and 14.6% of the total variance, respectively. These two modes are statistically distinguished by the criterion of North et al. (1982). The first MV-EOF mode shows a stronger MC around Lake Baikal at 850 hPa, an enhanced WNPSH over the northwestern Pacific at 500 hPa, and a strengthened SAH over the Tibetan Plateau (TP) at 200 hPa (Figs. 5a, 5c, and 5e, respectively).
6369
This pattern represents a cooperatively enhanced MC, WNPSH, and SAH. The MC and WNPSH are highly pronounced, with well-organized barotropic structures. The second MV-EOF mode, however, is featured by an enhanced MC at 850 hPa to the east of Lake Baikal and a northwestward shift of WNPSH at 850 and 500 hPa, showing an eastward shift of MC over northeastern China in the vertical levels and an inland northwestward propagation of WNPSH over the WNP (Figs. 5b and 5d, respectively). In contrast, the weakened SAH manifests as a southward cyclone over the TP at 200 hPa (Fig. 5f). Thus, the first mode of the EASM’s CISO is characterized by the simultaneous enhancement of the MC, WNPSH, and SAH, whereas the second mode shows an eastward and northwestward shift of the MC and WNPSH accompanied by the weakening SAH. Therefore, the first two modes compose the coupling mode of CISO in the EASM circulation at the intraseasonal time scale in boreal summer, which features the interaction among MC, WNPSH, and SAH. Figure 6a presents the phase relationship between the first two principal components (PCs) of CISO in the EASM. PC1 and PC2 exhibit an ISO with an alternative appearance in terms of the minimum and maximum phases. In the boreal summer, the maximum PCs occur in June, when the mei-yu reaches its peak phase, whereas the minimum PCs occur in late July, corresponding to the rainfall season over northeastern China. Wavelet analysis suggests that PC1 and PC2 exhibit a significant 40–80-day period in boreal summer (figures not shown), where the lag correlation between PC1 and PC2 indicates that PC1 leads PC2 by approximately 20 days (Fig. 6b). The phase difference between PC1 and PC2 suggests that the CISO state can be depicted in a two-dimensional phase space by using the two PCs (Fig. 6c). This pattern shows that the sequential days move clockwise around the origin and that the CISO becomes active with a large amplitude from June to August, thus indicating CISO’s critical role in modulating the EASM’s intraseasonal variability. Furthermore, three circulation indices are defined according to the vorticity of the 850-hPa MC (408–508N, 1008–1208E), 500-hPa WNPSH (258–37.58N, 1208– 1408E), and 200-hPa SAH (258–37.58N, 758–1008E) to examine their association with the coupling mode of EASM CISO. Figure 7 shows the normalized time series of the seasonal cycle, annual cycle, and CISO of the MC index (MCI), WNPSH index (WNPSHI), and SAH index (SAHI), respectively. In daily climatological and AC modes, the MCI first becomes positive in April, followed by the SAHI in early May and the WNPSHI in June, during the seasonal transition from boreal winter
6370
JOURNAL OF CLIMATE
VOLUME 29
FIG. 5. Spatial patterns of the (left) first and (right) second MV-EOF modes of the intraseasonal component of daily climatological winds at (a),(b) 850, (c),(d) 500, and (e),(f) 200 hPa. Anticyclone and cyclone over the crucial circulation areas are denoted by blue and red, respectively. All of the quantities are nondimensional because they were derived from the correlation coefficient matrix. The values in parentheses above each panel indicate the covariance, explained by each MV-EOF mode. The red and blue arrows indicate the regions of the positive and negative vorticity corresponding to the active region of MC, WNPSH, and SAH at 850, 500, and 200 hPa during boreal summer.
to summer (Figs. 7a,b). These three indices subsequently reach their peak in July and weaken and become negative in September, November, and December, respectively. Therefore, the evolution of the critical circulations in the EASM exhibits a stepwise seasonal
development, which is consistent with the subsequent onset of the subtropical, SCS, and WNP summer monsoons of the EASM (Zhu et al. 2011). Moreover, the withdrawal of the EASM is related to the phase change in the annual cycle of these three circulation indices. The
1 SEPTEMBER 2016
SONG ET AL.
6371
FIG. 6. (a) PCs of the first and second MV-EOF modes; (b) lag correlation between PC1 and PC2 and with itself, with the two red dashed lines representing a 95% confidence level based on Student’s t test; and (c) phase space (PC1, PC2) changes for days from 1 April to 31 October, quadrants are labeled with the four defined regions of the phase space.
CISOs of the MCI, WNPSHI, and SAHI are evident and present large amplitudes from April to August that alternately reach their peak during the boreal summer (Fig. 7c). This feature persists until late September. Moreover, there are phase lags between the circulation indices during the EASM CSIO cycle (Figs. 7d,e). The maximum of the SAHI leads PC1 by 2 days, whereas the counterparts of the MCI and WNPSHI lag PC1 by 0 and 15 days, respectively (Fig. 7d). Meanwhile, the MCI leads PC2 by 10 days, whereas the SAHI exhibits a significant relationship with PC2 (Fig. 7e). Thus, the intraseasonal variation of the MCI, WNPSHI, and SAHI indicates that the MC, WNPSH, and SAH are the three critical circulations in the dynamical wheel of the EASM and that their various phases determine the coupling modes of the EASM’s CISO.
b. Impacts on rainfall The EASM rainfall is prominent between 1008 and 1208E over East Asia (Fig. 1d), where the CISO of rainfall is remarkable in the latitudinal belt of 208–358N (Figs. 4b,c). The influences of the coupling modes in the EASM CISO on rainfall are identified by the winds and rainfall regressed against the first two PCs (Fig. 8) and by
their correlation with the three circulation indices (Fig. 9). For the positive PC1, stronger southwesterly winds are located over southern China in the lower troposphere, as are the strengthened WNPSH and an anticyclone over the Okhotsk Sea in far-eastern Russia (Fig. 8a). The PC1-regressed circulation in the upper troposphere presents a 200-hPa anticyclone centered over the TP (Fig. 8c). Meanwhile, the PC1-regressed rainfall is enhanced along the Yangtze River and is suppressed over the northeastern China and WNP regions (Fig. 8a), consistently with the synthesis between WNPSHI- and SAHI-associated precipitation (Figs. 9c,d). The PC1-regressed rainfall mode is similar to that associated with the Pacific–Japan (P-J) or East Asia– Pacific (EAP) pattern (Nitta 1987; Huang and Li 1988), thus showing the robust coupling between the WNPSH and SAH in the EASM’s CISO. In contrast, the PC2-regressed winds and rainfall anomalies exhibit a dipole mode that is characterized by the weakening of the MC and WNPSH with enhanced and suppressed rainfall to the north and the south of the Yangtze River, respectively (Fig. 8b). However, the SAH in the upper troposphere is weakened in the PC2regressed wind anomalies at 200 hPa (Fig. 8d). The
6372
JOURNAL OF CLIMATE
VOLUME 29
FIG. 7. Variations in the normalized time series critical circulation indices. The red line represents the MCI, the blue line represents the WNPSHI, and the green line represents the SAHI for (a) seasonal variation, (b) annual cycle, and (c) CISO. The lag correlations of PCs (CSIO wind and SST–Q1 modes) and each CISO circulation index (MCI, WNPSHI, and SAHI) during April–October are shown. Lag correlation coefficients with (d) PC1 and (e) PC2 (CISO winds mode). The WNPSHI and SAHI have been multiplied by a factor of 21 for comparison.
rainfall and upper-level circulation regressed against PC2 is opposite to the counterpart regressed against the SAHI (Fig. 9d), indicating the decoupling of the SAH to the MC and WNPSH in the mid-to-lower troposphere. The circulation regressed against the PC2 shows the comprehensive effect of MC and WNPSH on the CISO in terms of EASM rainfall. Thus, the CISO of tripole precipitation over the EASM region with its center over the mid-tolower reaches of the Yangtze River can be explained by the joint impacts of the MC, WNPSH, and SAH at intraseasonal time scale. In particular, the CISO of SAH is closely related to a tripole mode of rainfall over the WNP and the area to its north, implicating the importance of SAH in the first coupling mode of EASM CISO. In addition, the second mode of EASM CISO featured by the
coupling between the MC and WNPSH results in a dipole mode of rainfall anomalies (Figs. 9b,c). We further conducted EOF analysis on the CISO rainfall during April– October to verify whether the PC-regressed rainfall would resemble the primary characteristics of monsoon precipitation. Notably, the PC1- and PC2-regressed rainfall anomalies resemble the first two leading EOF modes of rainfall (shading in Figs. 8c,d), and the correlation coefficients between the first two PCs of circulation and their rainfall counterparts are 10.67 and 10.66, respectively. The PC1 and PC2 of rainfall anomalies also exhibit a 40–80-day period, and PC2 lags behind PC1 by 15 days, with a maximum correlation of 10.6. Therefore, the first two leading EOF modes of CISO of EASM rainfall suggest snapshots during the northeastward shift
1 SEPTEMBER 2016
SONG ET AL.
6373
FIG. 8. Anomalies of intraseasonal mean state associated with the first two MV-EOF analyses of winds: (a),(c) PC1 and (b),(d) PC2. (top) Regressed intraseasonal precipitation (shading; mm day21) and winds at 850 hPa (vectors; m s21). (bottom) The spatial distribution of the first two leading intraseasonal precipitation EOF modes in summer (shading) and winds at 200 hPa (vectors; m s21).
of the EASM’s rain belts, which are directly regulated by the EASM’s CISO, as indicated by the coupling of the MC, WNPSH, and SAH.
5. Linkages to the air–sea interaction a. SST–Q1 modes Because of the tropical ocean’s high heat capacity, the sea surface temperature (SST) exhibits a slow change during its seasonal cycle; therefore, the oscillation in the atmosphere over a single month can be considered to be a response to the SST forcing in the tropical ocean (Hoskins 2013). The gradual variation in SST may couple with Q1 and modulate the EASM’s CISO. Here, the MV-EOF approach is again applied to the CISO components of the SST and Q1 over the region 08–408N, 408–1608E in boreal summer. The results represent the feedback between atmospheric circulation and SST over the EASM region. Figure 10 shows the first two leading modes of CISO in the SST and Q1, which account for 37.4% and 28.3% of the total variance, respectively. In the first MV-EOF
mode of the SST and Q1, anomalous diabatic heating occurs over continental Asia, with its center settling over southern India, the southern edge of the TP, and Japan (Fig. 10a). Conversely, the diabatic cooling anomalies mainly occur over the Arabian Sea and the WNP, with three cooling centers situated over western Sumatra, the SCS, and the open ocean near southern Japan. Further, tropical SST anomalies exhibit a tripole mode, manifesting as cooling to the south of 108N and around the Japan Sea and warming over the Bay of Bengal (BOB), SCS, and WNP (Fig. 10c). However, the tropical SST anomalies display a zonal dipole pattern in the northern Indian Ocean (08–208N, 458–908E). The configuration between the SST and Q1 is quite different over the WNP versus the Arabian Sea. Because tropical diabatic heating is primarily contributed by condensation heating, the diabatic heating (cooling) over the tropical ocean suggests local enhanced (suppressed) convection. Over the WNP, the warm (cold) SST is overlaid by diabatic cooling (warming), thus suggesting that the local convection is suppressed (encouraged) to warm (cool) the SST by
6374
JOURNAL OF CLIMATE
VOLUME 29
FIG. 9. (a) The time–latitude diagrams of CISO of precipitation (mm day21). Spatial distribution of the correlations during April–October between (b) MCI, (c) WNPSHI, and (d) SAHI and precipitation (shading; mm day21) and winds (vectors; m s21) on an intraseasonal scale on corresponding levels, with red (blue) shading denoting significant positive (negative) correlations.
increasing (decreasing) the shortwave radiation and decreasing (increasing) the surface latent heat flux. However, diabatic cooling is superimposed on the cold SST over the Arabian Sea. This pattern demonstrates that the suppressed convection related to diabatic cooling should result from the local cold SST in the Arabian Sea, which decreases the supply of energy to the atmosphere and the convective instability via the air–sea interaction. As a result, the anomalous tropical convection is the reason for rather than the result of the local SST anomalies over the WNP, whereas the SST anomalies determine the convection over the Arabian Sea. In the second SST–Q1 mode (Figs. 10b,d), the anomalous diabatic cooling and heating centered over land straddle the Yangtze River over East Asia (Fig. 10b), whereas over the ocean, the diabatic cooling is centered over the Arabian Sea and the SCS. In addition, warm SST anomalies are distributed from the Arabian Sea to the SCS crossing the BOB (Fig. 10d). A distinct zonal gradient of SST anomalies is then generated between the SCS and the Philippines. The SST anomalies exhibit a
meridional alternative pattern, but their phase differs from that of the first mode in the northwestern Pacific Ocean. Anomalous diabatic cooling is collocated with the warm SST over the Arabian Sea and the SCS, thus implying a local SST response to atmospheric forcing. The interaction between the EASM CISO and the surrounding SST anomalies can be indicated by their leading modes. Figures 10e,f show the first two PCs of CISO with respect to the EASM and SST–Q1 anomalies. The correlation coefficients of the first two PCs are 10.64 and 10.71, respectively, with statistical confidence at the 95% level. Notably, the EASM CISO leads to changes in SST for several days during June–September, implicating the EASM CISO modulation on the surrounding SST anomalies, which, in turn, affect the underlying condition of the EASM. During the positive PC1 phase, the coupled mode of the EASM generally exhibits an enhanced WNPSH with suppressed rainfall over the WNP, leading to local clear skies and thus greater shortwave radiation to warm the SST. To the south and north of the WNPSH, however, abundant rainfall and heavy cloud cover may reduce
1 SEPTEMBER 2016
SONG ET AL.
6375
FIG. 10. Spatial patterns of the two leading MV-EOF modes of the CISO of (a),(b) daily apparent heating source (W m22) and (c),(d) sea surface temperature (K). The (normalized) (e) PC1 and (f) PC2 time series, with solid lines depicting the Q1 and SST modes and dotted lines depicting intraseasonal wind modes. All of the quantities are nondimensional because they were derived from the correlation coefficient matrix. The percentages in parentheses above (a),(c) and (b),(d) indicate the covariance explained by MV-EOF1 and MV-EOF2, respectively. The correlation coefficient between the first (second) leading SST–Q1 modes and the PC1 (PC2) shown in Fig. 6a is presented in the parentheses above (e) and (f), respectively.
shortwave radiation reaching the sea surface. Subsequently, the warming and cooling SST in the WNP may not only enhance the meridional SST gradient (Fig. 10c) but also modulate the summer rainfall over the Yangtze River. We noted that the first two PCs of EASM CISO always lead the counterparts of CISO in the SST–Q1 mode, implying atmospheric forcing on the SST at the intraseasonal time scale. Although the PCs of the EASM CISO lead the SST–Q1 CISO by 10 days, the relationship between WNPSHI and SST–Q1 CISO is distinct in the two leading modes (Figs. 11a,b). The SAHI and MCI lead the PC1 of SST–Q1 CISO by 15 and 5 days, whereas the WNPSHI lags the PC1 by 5 days (Fig. 11a). The SAHI, MCI, and WNPSHI lead the PC2 of SST–Q1 CISO by 30, 20, and 5 days, respectively (Fig. 11b). This result indicates that the feedback between the WNPSH associated with the EASM CISO and the air–sea coupling system over East Asia may support the life cycle of the EASM CISO in boreal summer. The SAH and MC may cause the WNPSH to regulate the
local air–sea interaction over the WNP. Specifically, the enhanced SAH and MC may enhance the midlatitude westerly winds, which strengthen the WNPSH to shift northward. Subsequently, the sea surface would receive more shortwave radiation to warm the SST because of the lower cloud cover and weaker convection related to the enhanced WNPSH. In addition, the air–sea coupling system over East Asia may affect the strength of the EASM CISO by changing Q1 over the WNP. If Q1 is suppressed with warmer SSTs over the WNP, a positive vorticity source (›Q1/›y . 0) would emerge to the north of the WNP, where the cyclonic vorticity of the air column in situ is increased, thereby weakening the ascending and EASM convection. Finally, the EASM CISO would be dampened, leading to enhanced convection and colder SSTs over the WNP.
b. Model simulation An atmospheric general circulation model (AGCM) of ECHAM5.4 (Roeckner et al. 2003) is used to further
6376
JOURNAL OF CLIMATE
VOLUME 29
FIG. 11. The lag correlations of PCs (CSIO wind and SST–Q1 modes) and each CISO circulation index (MCI, WNPSHI, and SAHI) during April–October. Lag correlation coefficients with (a) PC1 and (b) PC2 (SST–Q1 modes).
examine the role of SST forcing in maintaining the EASM CISO. We used the version with a spectral T63 horizontal resolution and 31 vertical levels (T63L31); the vertical domain of the model covered the surface level to 10 hPa. The daily SSTs used to drive this model were the same as those in the data analysis. The 30-yr (1982–2011) averaged circulation, forced by the daily SSTs, was defined as the climatology of the model simulation. Figure 12 shows the model-simulated climatological mean winds and rainfall during May–August and the
seasonal variation of rainfall at 1158E. In contrast, the model is able to accurately simulate the observed wind patterns, the tropical low-level westerly jet, the WNPSH at 850 hPa, and the tropical rainbands. The WNPSH is also clearly observed at 500 hPa in the simulation results (Fig. 12b). The SAH, characterized by a large-scale anticyclone centered on the TP, is present in the upper troposphere (Fig. 12c). The simulated seasonal variation of subtropical rainfall along 1158E exhibits a northward shift during May–August, which is much weaker than
FIG. 12. As in Fig. 1, but for the model simulations.
1 SEPTEMBER 2016
SONG ET AL.
6377
FIG. 13. As in Fig. 4, but for the model simulations. (f) Wavelet analysis (Torrence and Compo 1998) of the residual (shading), where the enclosed contours indicate regions of 95% confidence for a red-noise process, and the thick blue dashed line denotes the cone of influence, outside of which the edge effects become important.
the observations and is limited to south of 308N (Figs. 1d and 12d). However, a subtropical rainband originating from 258N is still present, and it rapidly shifts northward after mid-May and reaches its maximum in August. Relative to the AC component, the seasonal variation of subtropical and tropical rainfall exhibits a clear fluctuation during this period. The EASM CISO was also simulated. The simulated seasonal variation of rainfall at 308N, 1158E presents two peaks occurring in late June and August, respectively (Fig. 13a), corresponding to the double AC peaks, which are distinct compared with the single peak in the observations. In addition, the observed maximum rainfall is approximately 12 mm day21, but it is only 4 mm day21 in the simulation (i.e., only one-third that of the observations). The simulated maximum rainfall occurs in August, which is later than in the observations (late June). Therefore, this model fails to simulate the observed seasonal variation of rainfall in the lower reaches of the Yangtze River. However, the enhanced rainfall CISO is captured by the model (Fig. 13b), indicating the enhancement of local rainfall at the intraseasonal time scale during the mei-yu season. Although the time series of the AC-removed rainfall at 408N, 1158E exhibits a CISO during May–August (Fig. 13c), its amplitude and evolution
are quite different from the observations (Fig. 4c). In contrast, the seasonal variation and the AC of EASMI are well reproduced by the model (Fig. 13d), indicating a strong ability in simulating the seasonal transition between winter and summer. Moreover, the time series of EAMI without AC shows a marked CISO (Fig. 13e), but its 40– 80-day period is insignificant in the wavelet analysis. The coupling modes of CISO in EASM can also be partly captured by the AMIP run. Figure 14 shows the first two model-simulated MV-EOF modes of the EASM’s CISO. In simulation, the first mode of circulation at different levels (Figs. 14a,c,e) is similar to that of the observations, which is characterized by the jointly enhanced MC, WNPSH, and SAH (Figs. 5a, 5c, and 5e, respectively). In the second mode, however, the simulated 500-hPa WNPSH is weaker and the 200-hPa SAH is stronger than the observations (Figs. 14d,f), and both are located farther to the north (Figs. 5d,f). Therefore, the model has a spatial bias and fails to simulate the second mode of the observed EASM CISO. Figure 15 shows the first two PCs and the time series of the MCI, WNPSHI, and SAHI in the observations and the model simulation, respectively. The high correlation coefficient between the simulated and observed PC1 suggests that the SST forcing can accurately reproduce
6378
JOURNAL OF CLIMATE
VOLUME 29
FIG. 14. As in Fig. 5, but for the model output winds.
the first leading mode of the EASM CISO. However, the AGCM cannot accurately represent the observed PC2. For each component of the EASM CISO, the MCI could be well simulated, whereas the AGCM cannot capture the observed evolution of the WNPSHI and SAHI. Because the SST in the observed second mode of SST– Q1 acts as a response of SST to the atmospheric circulation, the AGCM cannot accurately represent the
observed air–sea interaction, resulting in the failure in simulating PC2 of the EASM CISO.
6. Summary and discussion a. Summary The climatological intraseasonal oscillation (CISO) of the East Asian summer monsoon (EASM) is characterized
1 SEPTEMBER 2016
SONG ET AL.
6379
FIG. 15. Variation in circulation indices calculated from the observations (black line) and the model simulation results (red line) during the summer monsoon. The values in parentheses are for the correlation coefficient between observed and simulated circulation indices.
by the vertical and meridional interactions of monsoon circulation, with a stepwise northward shift of frontrelated rain belts during boreal summer, particularly from May to August. To reveal the vertical structure and the internal modes of the EASM CISO, as well as their interaction with the surrounding SSTs, we conducted harmonic and multivariate empirical orthogonal function (MV-EOF) analyses on the climatological daily winds, rainfall, diabatic heating, and SST for the period 1981–2010. The EASM CISO exists not only in diabatic heating, primarily contributed by condensation heating, but also in rainfall and circulation. The existence of the EASM’s CISO is determined by both the different timing of diabatic heating Q1 over East Asia due to the land–sea thermal contrast and the distinct phase relationship between the Q1 components. During the EASM season, the Q1 phase associated with condensation heating is centered over the WNP and East Asia. Thus, the horizontal gradient of Q1 produces a negative vorticity source to the north and east of the Q1 centers, favoring the local background ascending that triggers the CISO.
Moreover, the EASM CISO initiates with the change in surface sensible heating, which modulates the atmospheric convective instability, followed by the variation of monsoon convection, releasing condensation heating, as well as the regulation of radiative heating. The MV-EOF results reveal that the first two modes of the EASM CISO are mainly characterized by the coupling of the Mongolian cyclone (MC) around Lake Baikal at 850 hPa, the WNP subtropical high (WNPSH) at 500 hPa, and the South Asian high (SAH) over the Tibetan Plateau (TP) at 200 hPa. The first leading mode shows a jointly enhanced MC, WNPSH, and SAH accompanied by a tripole rainfall anomaly of the strong mei-yu and baiu fronts around the lower reaches of the Yangtze River and southern Japan, whereas the rainfall was suppressed over northeastern Asia and the WNP. The second leading mode, which indicates the eastward and northwestward propagation of the enhanced MC and WNPSH with the weakened SAH, is associated with a dipole of the rainfall anomaly, with abundant and deficient rainfall over the northeastern and southeastern Asia–Pacific regions, respectively. The first two CISO
6380
JOURNAL OF CLIMATE
modes of daily SST and atmospheric Q1 are closely associated with the counterparts of the EASM CISO, except for a phase difference. Warming SST caused by the suppressed convection would enhance the WNPSH and cause the triple rain belts over East Asia through the coupling modes of the EASM CISO. In this study, we also applied ECHAM5.4, with a version at T63L31, to perform the AGCM simulations forced by daily SST from 1981 to 2010. The results show that the simulated first CISO leading mode is consistent with the observed mode, suggesting the critical role of seasonal variations of daily SST in the EASM CISO over the Asia–Pacific region. However, the AGCM runs fail to realistically resemble the second CISO leading mode. The MC CISO is well reproduced by the AGCM, but the simulated CISO of the WNPSH and SAH is poor compared with the observations owing to the lack of feedback of atmosphere forcing.
b. Discussion Since the MJO is the most prominent ISO in the tropics, one might instinctively consider it as the main source of the EASM CISO. Note that the MJO events are distinct in the individual year, and the effects of MJO on the climate-mean state exist possibly due to their phase lock to the annual cycle. To examine the phase locking of MJO to the annual cycle, we have applied the statistical method of a box-and-whisker plot (Tukey 1977) to the daily real-time multivariate MJO indices (RMM1 and RMM2), MJO intensity, and phase during April–October for the period 1981–2010 (Wheeler and Hendon 2004). There are three MJO-like events observed during May– August, but its period is less than 40 days. In addition, less than 40% of OLR CISO over the southern WNP could be linked to the MJO events, implicating the limited contribution of MJO to the EASM CISO modes. The possible mechanism of EASM CISO could be attributed to the SST–Q1 CISO over the WNP, which exhibits a remarkable northward propagation of anomalous SST and Q1 as shown in Figs. 10a and 10b, respectively. Previous studies have ascribed the intraseasonal northward shifting of SST–rainfall pattern over the WNP to the local air–sea interaction (Lau and Nath 2009). A similar process is observed in the SST–Q1 CISO mode over the WNP (figures not shown). Since the basic flow over the WNP is featured by the strong southeasterly over the southern WNP and weak northeasterly over the northern WNP during boreal summer (May–October), the more-than-normal rainfall over the cold SST anomalies in the WNP could induce an abnormal near-surface cyclone in situ. It could enhance (weaken) the evaporation latent heating over the northern (southern) WNP via the wind–evaporation–
VOLUME 29
SST feedback (Xie and Philander 1994), which supports the northward propagation of SST CISO mode. In the meantime, the zonal land–sea thermal contrast over East Asia is modified to alter the EASM rainfall and circulation as required by the thermal wind relationship (He et al. 2008; Zhu et al. 2011; He and Liu 2016). The robust cycle of EASM CISO starts with the changes in WNPSH and MC in the lower troposphere, followed by the upper-level SAH variation when the rainfall anomalies over East Asia are enhanced. Finally, the EASM CISO and the SST–Q1 CISO over the WNP are coupled together to move northward. Wang and Xu (1997) have speculated that the regulation of seasonal stepwise planetary-scale circulation related to the atmospheric heating may be responsible for the phase-locking behavior of the CISO to an annual cycle. This appears to be the result of complex nonlinear interactions among atmosphere, ocean, and land. In our results, we found that the CISO is most evident during May–August, instead of the whole boreal summer season, and that it exhibits phase locking to the EASM establishment and the subsequent summer monsoon onset in the Asia–Pacific region. During this period, the annual cycle maximum of diabatic heating over the land and sea exhibits distinct phase lags (Fig. 2a), which may be responsible for the generation of the CISO signal in the summer monsoon region. Therefore, to better understand the mechanism of CSIO, further research should focus on the impacts of local air–sea interactions on the CISO over the Asia–Pacific region. In addition, the CISO can be observed over other summer monsoon regions in boreal summer (figure not shown), but its strength is weaker than that of the EASM CISO. The features of CISO over other monsoon regions and the mechanism for their formation also require further investigation. Acknowledgments. We thank three anonymous reviewers for their valuable comments that helped to improve our manuscript. This study is jointly supported by the National Basic Research Program of China (973 Program; 2015CB453202), the China NSF (41475057, 41221064, and 41505049), the special fund for public welfare industry (Meteorology; Grant GYHY201406019), and the CAMS key program (2015Z001). REFERENCES Chen, L.-X., and Z. Jin, 1984: The medium-range variations of the summer monsoon circulation system over East Asia. Adv. Atmos. Sci., 2, 224–233, doi:10.1007/BF02678135. ——, C. Zhu, W. Wang, and P. Zhang, 2001: Analysis of the characteristics of 30–60 day low-frequency oscillation over Asia during 1998 SCSMEX. Adv. Atmos. Sci., 18, 623–638, doi:10.1007/s00376-001-0050-0.
1 SEPTEMBER 2016
SONG ET AL.
Chi, Y., and Coauthors, 2015: Correlation between the onset of the East Asian subtropical summer monsoon and the eastward propagation of the Madden–Julian oscillation. J. Atmos. Sci., 72, 1200–1214, doi:10.1175/JAS-D-14-0038.1. Emanuel, K. A., J. D. Neelin, and C. S. Bretherton, 1994: On largescale circulations in convecting atmospheres. Quart. J. Roy. Meteor. Soc., 120, 1111–1143, doi:10.1002/qj.49712051902. Hantel, M., and H.-R. Baader, 1978: Diabatic heating climatology of the zonal atmosphere. J. Atmos. Sci., 35, 1180–1189, doi:10.1175/ 1520-0469(1978)035,1180:DHCOTZ.2.0.CO;2. Hazra, A., and V. Krishnamurthy, 2015: Space–time structure of diabatic heating in monsoon intraseasonal oscillation. J. Climate, 28, 2234–2255, doi:10.1175/JCLI-D-14-00280.1. He, J., and B. Q. Liu, 2016: The East Asian subtropical summer monsoon: Recent progress. J. Meteor. Res., 30, 135–155, doi:10.1007/s13351-016-5222-z. ——, and Coauthors, 2007: A review of recent advances in research on Asian monsoon in China. Adv. Atmos. Sci., 24, 972–992, doi:10.1007/s00376-007-0972-2. ——, and Coauthors, 2008: Discussion of some problems as to the East Asian subtropical monsoon. Acta Meteor. Sin., 22, 419–434. ——, L.-Y. Chang, and H. Chen, 2015: Meridional propagation of the 30- to 60-day variability of precipitation in the East Asian subtropical summer monsoon region: Monitoring and prediction. Atmos.–Ocean, 53, 251–263, doi:10.1080/ 07055900.2015.1017798. Hoskins, B., 2013: The potential for skill across the range of the seamless weather-climate prediction problem: A stimulus for our science. Quart. J. Roy. Meteor. Soc., 139, 573–584, doi:10.1002/qj.1991. Huang, R., and W. Li, 1988: Influence of heat source anomaly over the western tropical Pacific on the subtropical high over East Asia and its physical mechanism (in Chinese). Chin. J. Atmos. Sci., 12, 107–116. Kanamitsu, M., W. Ebisuzaki, J. Woollen, S.-K. Yang, J. Hnilo, M. Fiorino, and G. Potter, 2002: NCEP–DOE AMIP-II Reanalysis (R-2). Bull. Amer. Meteor. Soc., 83, 1631–1643, doi:10.1175/BAMS-83-11-1631. Lau, K.-M., and M.-T. Li, 1984: The monsoon of East Asia and its global associations—a survey. Bull. Amer. Meteor. Soc., 65, 114–125, doi:10.1175/1520-0477(1984)065,0114: TMOEAA.2.0.CO;2. ——, and L. Peng, 1987: Origin of low-frequency (intraseasonal) oscillations in the tropical atmosphere. Part I: Basic theory. J. Atmos. Sci., 44, 950–972, doi:10.1175/ 1520-0469(1987)044,0950:OOLFOI.2.0.CO;2. ——, G. Yang, and S. H. Shen, 1988: Seasonal and intraseasonal climatology of summer monsoon rainfall over East Asia. Mon. Wea. Rev., 116, 18–37, doi:10.1175/1520-0493(1988)116,0018: SAICOS.2.0.CO;2. Lau, N.-C., and M. J. Nath, 2009: A model investigation of the role of air–sea interaction in the climatological evolution and ENSO-related variability of the summer monsoon over the South China Sea and western North Pacific. J. Climate, 22, 4771–4792, doi:10.1175/2009JCLI2758.1. Li, J., and J. X. Wang, 2003: A new North Atlantic Oscillation index and its variability. Adv. Atmos. Sci., 20, 661–676, doi:10.1007/BF02915394. Liu, Y. M., G. X. Wu, H. Liu, and P. Liu, 2001: Condensation heating of the Asian summer monsoon and the subtropical anticyclone in the Eastern Hemisphere. Climate Dyn., 17, 327– 338, doi:10.1007/s003820000117. ——, ——, and R. C. Ren, 2004: Relationship between the subtropical anticyclone and diabatic heating. J. Climate,
6381
17, 682–698, doi:10.1175/1520-0442(2004)017,0682: RBTSAA.2.0.CO;2. Lu, R., and Coauthors, 2014: The 30–60-day intraseasonal oscillations over the subtropical western North Pacific during the summer of 1998. Adv. Atmos. Sci., 31, 1–7, doi:10.1007/ s00376-013-3019-x. Miao, C., and Coauthors, 2006a: The anomalous features of the northeast cold vortex during the first flood period in the last 50 years and its correlation with rainfall in south China (in Chinese). Chin. J. Atmos. Sci., 30, 1249–1256. ——, Z. Wu, and J. He, 2006b: Relationship among the Northern Hemisphere annual mode, the northeast China cold vortex and precipitation during the first yearly rainy period in south China (in Chinese). J. Trop. Meteor., 22, 593–599. Newell, R. E., and Coauthors, 1969: The energy balance of the global atmosphere. The Global Circulation of the Atmosphere, G. A. Corby, Ed., Royal Meteorological Society, 42–90. ——, J. W. Kidson, D. G. Vincent, and G. J. Boer, 1975: The General Circulation of the Tropical Atmosphere and Interactions with Extratropical Latitudes. Vol. 2. MIT Press, 371 pp. Nitta, T., 1987: Convective activities in the tropical western Pacific and their impact on the Northern Hemisphere summer circulation. J. Meteor. Soc. Japan, 65, 373–390. North, G. R., T. L. Bell, R. F. Cahalan, and F. J. Moeng, 1982: Sampling errors in the estimation of empirical orthogonal functions. Mon. Wea. Rev., 110, 699–706, doi:10.1175/ 1520-0493(1982)110,0699:SEITEO.2.0.CO;2. Ren, R. C., Y. M. Liu, and G. X. Wu, 2007: Impact of South Asia high on the short-term variation of the subtropical anticyclone over western Pacific in July 1998 (in Chinese). Acta Meteor. Sin., 65, 183–197. Reynolds, R. W., T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. G. Schlax, 2007: Daily high-resolution-blended analyses for sea surface temperature. J. Climate, 20, 5473–5496, doi:10.1175/2007JCLI1824.1. Roeckner, E., and Coauthors, 2003: The atmospheric general circulation model ECHAM5. Part I: Model description. Max Planck Institute Rep. 349, 140 pp. [Available online at https://www.mpimet.mpg.de/fileadmin/publikationen/ Reports/max_scirep_349.pdf.] Tao, S., 1980: Heavy Rainfalls in China. Science Press, 225 pp. ——, and F. Zhu, 1964: The variation of 100 mb circulation over South Asia in summer and its association with march and withdraw of west Pacific subtropical high. Acta Meteor. Sin., 34, 385–395. ——, and L. X. Chen, 1987: A review of recent research on the East Asia summer monsoon in China. Monsoon Meteorology, C. P. Chang and T. N. Krishnamurti, Eds., Oxford University Press, 60–92. Tian, S.-F., and T. Yasunari, 1998: Climatological aspects and mechanism of spring persistent rains over central China. J. Meteor. Soc. Japan, 76, 57–71. Torrence, C., and G. P. Compo, 1998: A practical guide to wavelet analysis. Bull. Amer. Meteor. Soc., 79, 61–78, doi:10.1175/ 1520-0477(1998)079,0061:APGTWA.2.0.CO;2. Tukey, J. W., 1977: Exploratory Data Analysis. Addison-Wesley, 688 pp. Wang, B., 1992: The vertical structure and development of the ENSO anomaly mode during 1979–1989. J. Atmos. Sci., 49, 698–712, doi:10.1175/1520-0469(1992)049,0698:TVSADO.2.0.CO;2. ——, and H. Rui, 1990: Synoptic climatology of transient tropical intraseasonal convection anomalies: 1975–1985. Meteor. Atmos. Phys., 44, 43–62, doi:10.1007/BF01026810.
6382
JOURNAL OF CLIMATE
——, and X. Xu, 1997: Northern Hemisphere summer monsoon singularities and climatological intraseasonal oscillation. J. Climate, 10, 1071–1085, doi:10.1175/1520-0442(1997)010,1071: NHSMSA.2.0.CO;2. ——, and LinHo, 2002: Rainy season of the Asian–Pacific summer monsoon. J. Climate, 15, 386–398, doi:10.1175/1520-0442 (2002)015,0386:RSOTAP.2.0.CO;2. ——, Z. Wu, J. Li, J. Liu, C.-P. Chang, Y. Ding, and G. Wu, 2008: How to measure the strength of the East Asian summer monsoon. J. Climate, 21, 4449–4463, doi:10.1175/ 2008JCLI2183.1. ——, J. Liu, H.-J. Kim, P. J. Webster, and S.-Y. Yim, 2012: Recent change of the global monsoon precipitation (1979–2008). Climate Dyn., 39, 1123–1135, doi:10.1007/ s00382-011-1266-z. Wheeler, M. C., and H. H. Hendon, 2004: An all-season real-time multivariate MJO index: Development of an index for monitoring and prediction. Mon. Wea. Rev., 132, 1917–1932, doi:10.1175/1520-0493(2004)132,1917:AARMMI.2.0.CO;2. Xie, P., and P. A. Arkin, 1997: Global precipitation: A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs. Bull. Amer. Meteor. Soc., 78, 2539–2558, doi:10.1175/1520-0477(1997)078,2539: GPAYMA.2.0.CO;2. Xie, S.-P., and S. G. H. Philander, 1994: A coupled ocean-atmosphere model of relevance to the ITCZ in the eastern Pacific. Tellus, 46A, 340–350, doi:10.1034/j.1600-0870.1994.t01-1-00001.x.
VOLUME 29
Yanai, M., and T. Tomita, 1998: Seasonal and interannual variability of atmospheric heat sources and moisture sinks as determined from NCEP–NCAR reanalysis. J. Climate, 11, 463–482, doi:10.1175/1520-0442(1998)011,0463:SAIVOA.2.0.CO;2. ——, C. Li, and Z. Song, 1992: Seasonal heating of the Tibetan Plateau and its effects on the evolution of the Asian summer monsoon. J. Meteor. Soc. Japan, 70, 319–351. ——, J. H. Chu, T. E. Stark, and T. Nitta, 1976: Response of deep and shallow tropical maritime cumuli to largescale processes. J. Atmos. Sci., 33, 976–991, doi:10.1175/ 1520-0469(1976)033,0976:RODAST.2.0.CO;2. Zhao, P., and Coauthors, 2007: Onset of southwesterly wind over eastern China and associated atmospheric circulation and rainfall. Climate Dyn., 28, 797–811, doi:10.1007/s00382-006-0212-y. Zhu, C., T. Nakazawa, J. Li, and L. Chen, 2003: The 30–60 day intraseasonal oscillation over the western North Pacific Ocean and its impacts on summer flooding in China during 1998. Geophys. Res. Lett., 30, 1952, doi:10.1029/2003GL017817. ——, W. Lee, H. Kang, and C. Park, 2005: A proper monsoon index for seasonal and interannual variations of the East Asian monsoon. Geophys. Res. Lett., 32, L02811, doi:10.1029/2004GL021295. ——, X. Zhou, P. Zhao, L. Chen, and J. He, 2011: Onset of East Asian subtropical summer monsoon and rainy season in China. Sci. China Earth Sci., 54, 1845–1853, doi:10.1007/s11430-011-4284-0. Zhu, Q., J. He, and P. Wang, 1986: A study of circulation differences between East-Asian and Indian summer monsoons with their interaction. Adv. Atmos. Sci., 3, 466–477, doi:10.1007/BF02657936.