Rather than fishing statistical predictors in a variety of oceanic and atmospheric fields, we examine only SSTAs since they are the dominated lower boundary forcing conditions. Following Li and Wang (2018), we focus on two types of precursors, i.e., seasonal mean SSTAs in spring and SST tendency defined as the difference in SSTAs between spring and winter (the former minus the latter). To determine the most meaningful predictors for JA NECR, we calculate correlation coefficients between the SSTs and JA NECRI. Particular attention is paid to the regions where the correlations are significant at the 95% confidence level.
4.1 SST predictors
To explore potential precursors of JA NECR, we investigate correlations of the NECRI with the SSTAs in the preceding winter (December–January–February, DJF) and spring (March–April–May, MAM), and the SST tendency between MAM and DJF. To establish an empirical prediction model, the entire period of 1963–2019 (N=57 years) was separated into a calibration period of 1963–2002 (n1=40) and an independent validation period of 2003–2019 (n2=17 years). Fig. 3 shows the correlation maps between the JA NECRI and the SSTs in DJF and MAM during 1963–2002. Three empirical predictors with definite physical bases are discovered.
In DJF, there are significant positive correlation coefficients (CCs) in the TNIO, South China Sea, western part of the Maritime Continent, and negative CCs in the tropical central–Eastern Pacific Ocean (Fig. 3a). In the following MAM, enhanced positive CCs remain over the TNIO, while the prominent positive CCs in the tropical central–Eastern Pacific Ocean disappear (Fig. 3b). Thus, we extracted the first predictor index; that is, the MAM TNIO index (TNIOI, SST averaged over 0–15°N, 50–120°E), which has a significant CC of 0.63 with the JA NECRI. Moreover, the positive relationship between the MAM TNIOI and the JA NECRI is best represented when the values of the NECRI and TNIOI are both greater than 0.5 standard deviation (Fig. 4). For the 40-year calibration period of interest, a high (low) TNIOI in MAM corresponds to 7 (6) extreme positive (negative) NECRI, while the opposite example is only one in JA (Fig. 4a). This finding suggests that the MAM TNIO SSTAs in MAM contribute importantly to the interannual variations of the JA NECR.
Previous studies have revealed that SST tendency can be considered as an important predictor of short-term climate (Wang et al. 2013; Yim et al. 2014, 2015; Wen et al. 2017). Two additional tendency predictors are selected for the NECRI. The first tendency is a long-term dipolar-like SST tendency with warming in the West Pacific and cooling in the tropical central–Eastern Pacific Ocean from DJF to MAM (Fig. 3c). Therefore, based on Fig. 3c, an index that represents this dipolar-like SST tendency is defined as WEPDT index (WEPDTI, Table 1), which has a significant CC of 0.70 with the JA NECRI (Table 1). For the 40-year calibration period of interest, a high (low) WEPDTI corresponds to 8 (7) extreme wet (dry) JA NECR, while there are no years in which the opposite is true (Fig. 4b). This finding shows that the dipolar-like SST tendency between the Western and Eastern Pacific also contribute to the interannual variations of JA NECR.
The second tendency predictor is the tripolar SSTA over the North Atlantic Ocean from DJF to MAM (Fig. 3c). This tripolar-like pattern is accompanied with warming tendencies in the North tropical Atlantic (0°–20°N) and Northeast Atlantic (46°–60°N) and cooling tendencies in the subtropical Atlantic (26°–34°N). Therefore, based on Fig. 3c, an index reflecting this tripolar SST tendency is defined as NAT index (NATI, Table 1), which has a significant CC of 0.53 with the JA NECRI. For the 40-year calibration period of interest, a high (low) NATI corresponds to 7 (7) extreme wet (dry) JA NECR, while the opposite example is only one year (Fig. 4c).
4.2 Discussion of physical linkage between the predictors and JA NECR
For the empirical prediction, it is vital to comprehend the sources of predictability, i.e., the physical processes and mechanisms linking the predictors and predictand. In this section, we give a brief discussion of their physical linkage to the interannual variations of JA NECR by both observation analysis and numerical modeling experiments.
4.2.1 Spring TNIO SST and NECR
To illustrate the influence of TNIO on JA NECR, regressions of the SST anomalies, stream function and horizontal wind anomalies at 850 hPa, and precipitation anomalies in JA against the MAM TNIO index are computed for the period of 1963–2002. It is clear that the TNIO SSTAs strongly persist from MAM to JJA, and the warm SSTAs can spread eastward to the eastern Philippine Sea (Fig. 5a). These warm SSTAs give rise to easterly anomalies over the eastern Indian Ocean–western Pacific (Fig. 5b) via the Gill-Matsuno mechanism (Gill 1980), which further result in an anomalous anticyclone over the Philippine Sea (i.e., WNPAC). The associated southwesterly anomalies along the northwestern flank of the WNPAC transport more warm moisture into NEC, leading to strengthened rainfall over NEC (Fig. 5c). This is consistent with Wu et al. (2010), which revealed that the TNIO warming plays a dominant role in modulating the WNPAC variability in late summer.
To confirm the possible mechanisms and assess the contribution of TNIO SSTAs to the JA NECR variability, we perform a set of numerical experiments using the CAM5.3 model forced with SSTAs plus the climatological mean SST in the tropical Indian Ocean (30°S–30°N, 30°–110°E; referred to at the TNIO run). Here, the SSTA forcing is derived by regressing the monthly SST anomalies against the MAM TNIOI during 1963–2002. The atmospheric circulation response to the TNIO SSTA forcing is represented as the ensemble-mean difference between the TNIO and control runs (former minus latter).
Figure 6 shows the responses of horizontal wind and stream function at 850 hPa, precipitation to the TNIO SSTA forcing in JA. There are positive precipitation responses in most parts of the TNIO (Fig. 6b), indicating that the warm SSTAs in the TNIO tend to enhance the local convection activity. Significant easterly anomalies prevail over the northern tropical Indian Ocean–Western Pacific as a Gill-Matsuno-type response, leading to anti-cyclonic flows to the north of the easterly (Fig. 6a). In particular, a strong anti-cyclonic response occurs over the Philippine Sea (Fig. 6a), consistent with its observed counterpart (Fig. 5b). The southwesterly anomalies at the western flank of the anomalous anticyclone favor more water vapor from the ocean to NEC, resulting in enhanced precipitation over there (Fig. 6b), though the amplitude of the simulated precipitation anomalies is weaker than that of the observations. Nevertheless, similarity between the model simulations and observations confirms that the TNIO warming can exert a significant impact on the NEC precipitation variability through enhancing the WNPAC in late summer. Also, our result is supported by the previous studies that TNIO warming can intensify the WNPAC through the eastward propagation of Kelvin waves and/or modulating the Hadley circulation (Sui et al. 2007; Yang et al. 2007; Wu and Zhou 2008; Xie et al. 2009; Wang et al. 2013). It should be noted that the TNIO warming tends to be lagged the El Niño events by one season, and therefore the spring TNIO warming plays a role in linking the previous winter ENSO and following summer weather and climate over East Asia (Xie et al. 2009). In addition to ENSO, ocean dynamics and local ocean–atmosphere interaction also contribute to the TNIO SST variability (Huang and Kinter 2002; Du et al. 2009).
4.2.2 SST tendency index 1—WEPDTI and NECR
Same as the TNIO index, regressions of large-scale lower boundary anomalies in JA against the WEPDTI are computed to examine influence of the dipolar-like SST tendency in the Pacific Ocean from DJF to MAM. It is indicated that the warming SSTA tendency in the Western Pacific and the cooling SSTA tendency in the tropical Eastern Pacific associated with the WEPDTI can endure until summer (Fig. 7a), which corresponds to the developing phase of La Niña. This east-west dipole SST tendency can strengthen the equatorial easterly anomalies through zonal thermal gradients, thus leading to central–eastern Pacific cooling in the following JA season (Fig. 7a). Such central–eastern Pacific cooling generates westward propagating descending Rossby waves through suppressed convection over the western North Pacific (Li and Wang 2015). This Rossby wave response enhances the WNPAC, resulting in stronger southwesterlies that transport water vapor further north (Fig. 7b), thereby enhancing NECR (Fig. 7c).
To confirm the influence of tropical Pacific SSTAs to the JA NECR variability, we perform another set of numerical experiments using the CAM5.3 model forced with SSTAs plus the climatological mean SST in the tropical Pacific Ocean (20S°–20°N, 110°–270°E; referred to at the WEP run). The SSTA forcing is derived by regressing the monthly SST anomalies against the WEPDTI during 1963–2002. Figure 8 shows the differences in the precipitation, horizontal wind and stream function at 850 hPa in response to the tropical Pacific SSTAs forcing in JA. It is indicated that the atmospheric circulation responses are strong in the tropics in the WEP run. Especially, there is a pair of anomalous low-level cyclones over the eastern Pacific–Atlantic region and a pair of anomalous low-level anti-cyclones over the western Pacific (Fig. 8a), and therefore leads to more precipitation anomalies over NEC (Fig. 8b), which are consistent with the results obtained from observations (Fig. 7b). Therefore, the model simulations support the significant influence of the WEPDTI-related SSTAs on the WNPAC and thus the NEC precipitation variability in late summer
4.2.3 SST tendency index 2—NATI and NECR
Figure 9 displays regressions of SST and atmospheric circulation anomalies in JA against the NATI. It is shown that the North Atlantic tripolar-like SST tendency, particularly the warming SST tendency in the tropical North Atlantic (0–20°N), can persist from spring into late summer (Fig. 9a). The delayed impact of Atlantic SSTAs on JA NECR is possibly related to the long persistence of SSTAs from spring to summer (Hong et al. 2014; Zuo et al. 2019) and atmosphere–ocean interactions along the Pacific Intertropical Convergence Zone (Ham et al. 2013). Associated with the high-index phase of the NATI, there is a prominent anomalous anticyclone over the subtropical western North Pacific and a cyclone over the eastern North Pacific–Atlantic (Fig. 9b). Previous studies have suggested that the North Atlantic–WNPAC connection is related to an anomalous zonally overturning circulation extending westward from the tropical Atlantic to the tropical Pacific (Hong et al. 2014, 2015; Zuo et al. 2019). When the tropical Atlantic SSTAs are warmer than normal, enhanced local convection activity (Fig. 9d), low-level convergence and upper-level divergence (Fig. 9c) are clearly observed over the heating forcing. Consequently, enhanced upper-level convergence and low-level divergence are observed over the central tropical Pacific, which suppress the local convection activity and in turn trigger an anomalous low-level anticyclone over the western North Pacific as a Rossby-wave response. Moreover, an anomalous cyclone is observed over Northeast China (Fig. 9b). This anomalous cyclone over Northeast China and the enhanced WNPAC favor intensified NECR in JA (Fig. 9d).
To confirm the linkage between the NAT and JA NECR, numerical experiments are conducted with the CAM5.3 model forced with SSTAs plus the climatological mean SST in the North Atlantic Ocean (0°–70°N, 10°–70°W; referred to at the NAT run). Here, the SSTA forcing is derived by regressing the monthly SST anomalies against the NATI during 1963–2002. Figure 10 shows the ensemble-mean differences of the 850-hPa stream function and horizontal wind, 200-hPa velocity potential, and precipitation in JA between the NAT run and control run (former minus latter). It is indicated that the atmospheric circulation and precipitation in JA responses to the NAT SST forcing are strong, which are highly consistent with the observations. Particularly, there are apparent positive precipitation responses over the tropical North Atlantic (Fig. 10c), indicating that the warm SSTAs over there tend to enhance local convection activity. Correspondingly, a pair of anomalous low-level cyclonic flow responses occurs over the subtropical eastern Pacific–western Atlantic , and a pair of anomalous low-level anticyclone occurs over the western North Pacific in the experiments (Fig. 10a), the pattern of which is similar to its observed counterpart (Fig. 9b).
Moreover, there is significant anomalous low-level convergence and upper-level divergence over the tropical North Atlantic and opposite anomalies over the central tropical Pacific in the high-index phase of the NATI (Fig. 10b). Therefore, the tropical North Atlantic warming appears to be accompanied by an anomalous zonally overturning circulation extending westward from the tropical Atlantic to the central tropical Pacific in JA, which is associated with anomalous ascending motion over the tropical North Atlantic and opposite anomalies over the central tropical Pacific. The enhanced subsidence over the central tropical Pacific suppresses the local convection activity, which in turn triggers an anomalous low-level anticyclone over the western North Pacific as a Gill-type Rossby-wave response (Gill 1980). This result is consistent with previous numerical experiment result (Zuo et al. 2019), which revealed that warm SSTAs over the tropical North Atlantic favor the intensification of the WNPAC via modulating the Walker circulation over the Atlantic-Pacific Oceans in summer. Therefore, consistency between the simulations and observations confirms the impact of the North Atlantic SSTAs on the WNPAC and thus the precipitation anomalies over NEC in JA (Fig. 10c).
Previous studies have revealed that the tripole-like SSTAs in the North Atlantic can also exert an impact on the East Asian summer climate by triggering a set of wave trains propagating from the North Atlantic downward to Eurasian mid–high latitudes (Wu et al. 2009; Wu et al. 2012; Zuo et al. 2013; Li et al. 2018). As seen in Fig. 11a, which displays the regressed 500-hPa geopotential height anomalies in JA against the NATI, there exists a set of wave trains propagating from the North Atlantic middle latitudes downward to Northeast Asia. A similar wave-train pattern is observed in the ensemble-mean difference in 500-hPa geopotential height between the NAT and control runs (Fig. 11b). But it should be noted that amplitude of both the observed and simulated wave-train patterns is rather weak in the Eurasian mid–high latitudes.