Abstract

The characteristics and possible impact factors of the South China Sea summer monsoon (SCSSM) evolution from onset to withdrawal before and after 1993/94 are investigated using ERA-Interim, CPC rainfall, and OLR data. During the late-onset period of 1979–1993, the SCSSM was characterized by stronger onset intensity and a gradual withdrawal, resulting in a continuous, strong preflood season in Southern China and a slower rain-belt retreat from north to south China in September. In addition, the rain-belt in the Yangtze River basin persisted much longer during summer. However, during the early-onset period in 1994–2016, the SCSSM is associated with a weaker onset intensity and comparatively faster retreat. The advanced preflood season lasted intermittently throughout May and the whole eastern China precipitation lasted until October when it retreated rapidly, making the rain-belt in Southern China persist for an extended duration. Further analysis indicates that a strong modulation of SCS intraseasonal oscillation (ISO) on the SCSSM evolution is observed. There are two active low-frequency oscillations over the SCS in summer during the late-onset period but three during the early-onset period. The wet ISO in the Northwest Pacific propagating northwestward into the SCS and enhanced SCSSM ISO activity may contribute to the early onset and faster withdrawal after 1993/94. The effect of warm western Pacific sea surface temperatures (SST) on the SCSSM evolution is also discussed.

1. Introduction

Asian monsoons are one of the main energy and moisture sources of the global atmosphere [1]. The SCSSM has a significant position in the Asian monsoon, forming a vital link between the East Asian and South Asian summer monsoons [24]. Climatologically, the SCSSM onset is generally accompanied by a continuous eastward retreat of the western North Pacific (WNP) subtropical high, convection enhancement, and the reversal of low-level zonal wind in the SCS [5, 6]. The synoptic-scale circulation systems, the disturbance of mid-latitude, and the ISO activities can trigger SCSSM onset [715]. In addition, the sea surface temperature anomaly (SSTA) is closely related to the interannual variability of the SCSSM onset. A warm (cold) El Niño-Southern Oscillation (ENSO) event in the previous winter and the warm (cold) tropical Indian Ocean (TIO) in the spring can delay (advance) the SCSSM onset [1619].

There are fewer studies that focus on SCSSM withdrawal when compared to the SCSSM onset [20]. Hu et al. [21] found that SCSSM withdrawal mainly occurred due to the westward intrusion of the WNP subtropical high, which is accompanied by the retreat of the weakening low-level intertropical convergence zone (ITCZ) and rain-belt in the SCS-WNP. The factors contributing to the interannual variability of the SCSSM withdrawal may include tropical cyclones (TCs) and ISO [22]. Luo and Lin [23] suggested that El Niño can delay the monsoon onset and advance the SCSSM withdrawal, thus shortening the length of the summer monsoon season.

Recent research has shown that the SCSSM onset date has advanced by about two weeks since 1994 [2429]. Kajikawa and Wang [25] suggested that the enhanced northwestward intraseasonal variability (ISV) and TC activity from WNP could be responsible for the interdecadal change in the SCSSM onset. Xiang and Wang [30] indicated that the interdecadal advance of the Asian summer monsoon onset may be attributed to a grand La Niña-like pattern change in the Pacific, while that in the South China Sea (SCS) is primarily determined by the abrupt SST warming near the Philippine Sea. Liu et al. [31] highlighted the impact of southern Indian Ocean (SIO) on the SCSSM onset advancement after 1993. Many efforts have been made to understand the cause driving interdecadal change of the SCSSM onset, and interdecadal change in the SCSSM withdrawal after the mid-2000 has been found recently [32]. Considering the integrity of SCSSM evolution, particularly from onset to withdrawal, the following questions can be raised: (1) What changes have happened in the SCSSM evolution from onset to withdrawal since 1993? (2) Previous studies have highlighted the role of SSTA and ISO during this advanced SCSSM onset. What are the effects of ISO and SSTAs on the SCSSM evolution from onset to withdrawal in different periods? These questions will be answered in this paper.

2. Data and Methods

The datasets used in this paper include (1) the daily mean interpolated outgoing longwave radiation (OLR) data from the National Oceanic and Atmospheric Administration (NOAA) satellite with resolution of 2.5° latitude by 2.5° longitude from 1979 to 2016; (2) the daily mean ECMWF Interim Re-Analysis dataset (ERA-Interim) with resolution of 2.5° latitude by 2.5° longitude from 1979 to 2016 (the variables used in this paper include temperature and three-dimensional wind fields); (3) the monthly mean Hadley Centre sea surface temperature (SST) dataset for the period of 1979–2016 [33]; and (4) daily precipitation set from the CPC Unified Precipitation Project which is underway at NOAA Climate Prediction Centre (CPC).

In the present study, the SCSSM onset is identified by the seasonal transition of the mid- to upper-tropospheric meridional temperature gradient (MTG, ) from winter to summer, which is defined by Mao et al. [34] and Liu et al. [31]. Specifically, the timing of the SCSSM onset is defined by the moment when the pentadly area-averaged MTG in the mid-to-upper troposphere (500–200 hPa) changes from negative to positive and remains positive for at least three pentads over the SCS (10–20°N, 110–120°E). Then the first pentad from negative to positive phase is defined as the onset of the SCSSM. They suggest that MTG between 500 hPa and 200 hPa is more abrupt in assessing monsoon onset and consistent with the main climatological characteristics of the SCSSM onset, such as area-averaged zonal wind shear reversal, zonal wind reversal, and convection enhancement. There are 72 pentads per year. However, the 850 hPa zonal wind and convection are the two common and significant indices for identifying the SCSSM onset [21, 35]. In order to take low-level features for a more comprehensive consideration into account the, we corrected the onset date by combining circulation and the onset date defined by the National Climate Centre. The National Climate Centre defines the SCSSM onset by the first pentad when 850 hPa zonal wind in the SCS changes to westerly wind for two constant pentads stably and the average potential pseudoequivalent temperature is greater than 340 K. As a result, we corrected the onset date of 1982, 1997, 2000, 2008, and 2009. The SCSSM withdrawal pentad for 1979–2016 is also defined by MTG. The MTG-defined SCSSM onset and withdrawal originate from thermal differences between north and south, which is the essence of monsoon formation. In order to extract the more obvious signal caused by the SCSSM interdecadal transition, the years with later than average onset are regarded as the typical late-onset years during 1979–1993, which are 1981, 1982, 1985, 1987, 1991, and 1993. The years with an earlier than average onset are regarded as the typical early-onset years during 1994–2016, which are 1994, 1996, 2000, 2001, 2002, 2005, 2006, 2008, 2012, 2013, and 2014.

3. The Decadal Shift of the SCSSM before and after 1993/94

3.1. Difference of the SCSSM Onset Characteristics during the Two Epochs

To show the reliability of the MTG-defined SCSSM, the time series of the SCSSM onset date is plotted in Figure 1, which shows a significant advancement of the SCSSM onset date after 1993/1994 and the SCSSM onset time we defined is basically consistent with that defined by the National Climate Centre and its correlation is 0.58. The mean SCSSM onset occurs in about pentad 29 during 1979–1993, while during 1994–2016 it occurs in about pentad 27. The change is about two pentads and the standard deviation (STD) of the SCSSM onset time is 1.67. This interdecadal change of the onset date is consistent with previous studies [2426] and passes the sliding t-test (not shown). Based on the interdecadal change of SCSSM onset, the circulation and convection differences during the SCSSM onset should be researched further.

Figure 2 shows the composite OLR and 850 hPa winds of the SCSSM onset from before and after the onset pentad in each epoch. Lower OLR corresponds to heavy rainfall and strong convection. In the late-onset years of 1979–1993, before the SCSSM onset pentad (Figure 2(a)), a weak anticyclone is centered over the SCS, which makes the SCS be controlled by easterly wind. The active convection area is mainly located in the Bay of Bengal and Indochina Peninsula, extending to the south of China and the south of Japan. When the SCSSM onsets (Figure 2(b)), the developing convection from the Bay of Bengal expands to the northern part of the SCS along with weakening convection over the southern Indian Ocean, indicating that the ITCZ stretched northwards [25]. Southwesterly wind from the Bay of Bengal reaches the SCS resulting in the anticyclone of the Northwest Pacific retreating eastward. After the onset (Figure 2(c)), the active convection over the Arabian Sea and Indochina Peninsula is expanded into the SCS. The convection in the Bay of Bengal is significantly stronger than that in 1994–2016 and the southwesterly wind in the SCS is more intense.

In early-onset years of 1994–2016, the convection in the eastern equator of Indian Ocean is weaker than that in 1979–1993, before the onset pentad (Figure 2(d)). Another distinct convection area is over the Indonesian Maritime Continent. This convection is associated with tropical disturbances, such as easterly wave and MJO, which are favorable for the onset of the SCSSM [22]. When the SCSSM builds up (Figure 2(e)), both active convection over the Indonesian Maritime Continent and the eastern equator of Indian Ocean are expanded northward. The wind direction over the SCS abruptly changes to westerly direction. After the onset (Figure 2(f)), the convection over the Indonesian Maritime Continent moves northward to the Philippines and enters south of the SCS, but the southwesterly wind in the SCS weakens.

According to the above analysis, the convection source favoring monsoon onset has changed since 1993/94. In addition, it is clearly shown that the reversal and increment of low-level wind in SCS are more noticeable during 1979–1993 when the SCSSM onset is delayed. Figure 3(a) shows the composite 850 hPa zonal wind evolution over the SCS for typical years and climatology. In 1979–1993, the monsoon builds up late, corresponding to the zonal wind turning positive late. On the other hand, the zonal wind growth is peaking faster in late June. However, in 1994–2016, the monsoon builds earlier; the zonal wind turns positive in early May and grows slowly. The first peak is in mid-June. Thus, what is the relationship between monsoon onset date and intensity?

Due to the change of the zonal wind in the SCS, we expressed the zonal wind growth intensity by zonal wind over the SCS after the SCSSM onset pentad minus the SCSSM onset pentad, which objectively depicts the abrupt onset intensity of SCSSM. Time series of standardized OLR in onset pentad and zonal wind growth intensity in the SCS shows that positive zonal wind growth intensity generally corresponds to negative OLR (Figure 3(b)). This means that a strong zonal wind intensity increase in the SCS generally corresponds to an enhanced convection. Therefore, the convection and wind field in the SCS have obvious coupling characteristics. They show that the SCSSM onset process is a combination of dynamic and thermal. To reveal the relationship between the SCSSM onset date and intensity, we calculated the correlation coefficients between the regional averaged OLR and the zonal wind growth intensity with the SCSSM onset date in the period of 1979–1993 as well as in the period of 1979–2016 (Table 1). The SCSSM onset date is significantly related to the SCSSM onset intensity. There is a positive correlation between the zonal wind growth intensity and the onset date of SCSSM, and there is a negative correlation between the OLR and the onset date of SCSSM. However, the correlation coefficient of two variations during 1994–2016 is much higher and exceeds the 0.01 confidence level, while that during the 1979–1993 epoch fails in passing the 0.05 confidence level and this phenomenon is more obvious during 1994–2016. However, different from Hu et al. [20] and Torrence and Webster [36], sliding correlation between the SCSSM onset date with OLR and zonal wind growth intensity shows no significant interdecadal shift in the relationship between zonal wind growth intensity and OLR in the mid-1990s. This result illustrates that the late onset of the SCSSM is likely to correspond to a strengthened zonal wind and enhanced convection activity, which means that the SCSSM onset intensity is stronger. Therefore, the later (earlier) the SCSSM onset is, the stronger (weaker) the intensity of the SCSSM onset is. Between 1979 and 1993 (1994 and 2016), the SCSSM onset is generally late (early), so the zonal wind growth intensity is stronger (more weakened) and convection is more active (suppressive) in the SCS.

3.2. Differences of the SCSSM Withdrawal Characteristics during the Two Epochs

The circulation evolution and abrupt characteristics of the SCSSM onset of the two epochs are discussed above. After the SCSSM onset, the low-level zonal winds appear with a two-peak pattern in summer (Figure 3(a)). In the late-onset years of 1979–1993, the low-level zonal wind over the SCS grew rapidly from early June until it reached the first peak, and then it began to decrease. Subsequently, the zonal wind grew again and reached a second peak at the end of July. In the early-onset years of 1994–2016, the zonal wind peaked firstly in mid-June, with a weak intensity. Then, in early July, it started to grow rapidly until the end of July. After summer, what changes have happened to the SCSSM withdrawal since 1993/94?

To find the difference in SCSSM withdrawal characteristics between the two epochs, we must first figure out the climatological circulation evolution of 850 hPa wind and OLR during the SCSSM withdrawal. Figure 4 shows the latitudinal-time and longitudinal-time cross section of low-level wind and OLR averaged from August 15 to the end of October. Different from the abrupt features of the SCSSM onset, it takes longer for the SCSSM to withdraw. The southwesterly wind and convection retreat from the SCS do not occur simultaneously during the SCSSM withdrawal. Meridionally, the southwesterly wind retreats completely from the northern SCS to the southern SCS in mid-September, while the convection retreats gradually from the SCS until mid-October (Figure 4(a)). Different times for the wind and convection to retreat to the SCS can also be seen in the zonal variation (Figure 4(b)). This nonsimultaneous retreat of westerlies and convection in the SCS is due to the activities of TCs and the synoptic systems embedded in the basic flow of easterly winds, which can contribute to the active convection around the SCS [21]. However, the TC and other synoptic scale systems are also greatly modulated by the WNP anticyclonic/cyclonic anomalies [35], which are closely related to the SCSSM withdrawal. Therefore, it is necessary to consider the progressive retreat convection in the process of SCSSM withdrawal, especially in the determination of the SCSSM retreat date.

According to the above analysis, we use pentadly area-averaged MTG in the mid-to-upper troposphere (500–200 hPa) when it changes from positive to negative to calculate the date of the SCSSM withdrawal. Considering the inconsistency of convection and wind field retreat during the SCSSM withdrawal, the persistence of MTG turning negative cannot be maintained for three continuous pentads; therefore, the large-scale circulations such as wind field and convection are also taken into account and every pentad in which the MTG changes from positive to negative from September is checked and corrected. Figure 5 shows the time series of the SCSSM withdrawal pentad. The mean and standard deviation (SD) of the SCSSM retreat are pentads 55.3 (about October 8th) and 3.4, while the interannual fluctuation of the SCSSM retreat is stronger than that of the onset, which is consistent with Luo and Lin [23] and Hu et al. [37] who provided two different sets of SCSSM withdrawal dates characteristics. The correlation between the SCSSM withdrawal dates defined by MTG and by Hu et al. [21] is 0.4, which exceeds the 99% confidence level. We define pentads 53–57 as the withdrawal stages to describe the interannual variability of the SCSSM retreat process. However, the average withdrawal time from 1979 to 1993 is pentad 54.5, while the average withdrawal time from 1994 to 2016 is pentad 55.9. The retreat time was later in the later epoch but it failed to pass the Mann–Kendall test. In Hu et al.’s work [32], there was a significant delay of SCSSM withdrawal between the periods of 1995–2005 and 2006–2016, but the periods in this research were 1979–1993 and 1994–2016. This may be because the study period is different. But the SCSSM withdrawal also exhibits a delay after the mid-2000s if it was defined by our approach. For example, the average SCSSM withdrawal date from 1995 to 2005 was 54.4 pentads, and the average SCSSM withdrawal date from 2006 to 2016 was 57.7 pentads, which is almost consistent with Hu et al.’s result. From the composite results, the earlier the SCSSM onset, the later the SCSSM withdrawal after 1993.

Now that the SCSSM withdrawal dates have been defined, we will explore the circulation evolution during the retreat period. Figure 6 shows the climatological compositions of low-level wind and OLR during the SCSSM withdrawal period from 1979 to 2016. It can be seen that active convection is located in the Bay of Bengal to the Northwest Pacific before the SCSSM withdrawal pentad, and it retreats gradually from northeast to southwest (Figures 6(a) and 6(b)) and then completely withdraws from the SCS, since the SCSSM withdrawal pentad and convection centers shift to the Indonesia area (Figures 6(c) and 6(d)). Meanwhile, the low-level wind is already easterly before the SCSSM withdrawal pentad, but it is obvious that the northeasterly wind has become stronger on the SCSSM withdrawal pentad. To further objectively investigate the variability of circulation factors during the SCSSM withdrawal, the pentad-mean of 850 hPa zonal wind and OLR averaged over the SCS evolution are composited in Figure 7. There appears to be a significant decrease in zonal wind and a weakened convection whose value rises remarkably from about 220 W/m2 to around 240 W/m2 over the SCS in the SCSSM withdrawal pentad. Zonal wind and convection almost develop at the same time when the SCSSM onsets, but, during the SCSSM retreat, the westerly wind and convection retreat from the SCS do not occur simultaneously [21]. It was zonal wind shifting to easterly wind firstly and then convection retreat from the SCS (Figure 4). The time difference of them is about 20 days. The SCSSM withdrawal dates defined in this paper take into account not only the wind field but also OLR; therefore, the 850 hPa zonal wind in pentad (−3) is already easterly wind. It reveals the “abrupt change” of the SCSSM onset and the “gradual change” of the SCSSM retreat from the dynamics and thermal structure of atmosphere. Thus, the characteristics of the wind field and convection during the SCSSM withdrawal show that the retreat time defined by MTG is reliable. However, based on the interdecadal change of the SCSSM, what characteristics have changed on the SCSSM withdrawal?

To further study the differences of the SCSSM withdrawal before and after 1993, we composite the anomalies of 850 hPa wind and OLR during the SCSSM withdrawal period for typical years in Figure 8. During the period of 1979–1993, before the SCSSM withdrawal pentad, the low-level and convection anomalies are not significant in the SCS (Figure 8(a)). On the SCSSM retreat pentad, an enhanced easterly wind anomaly and increased OLR anomaly, meaning weakened convection, appear over the SCS and subtropical western North Pacific. The easterly wind anomaly and weakened convection extend eastward to the Bay of Bengal (Figures 8(b) and 8(c)). Comparison of the difference before and after the SCSSM withdrawal pentad (Figure 8(d)) shows that the SCSSM withdrawal leads to remarkable increases in easterly wind from the SCS to the Bay of Bengal, as well as an anomalous anticyclone over the western North Pacific, which accompanies a weakened convection over the east of Philippines and east of the SCS.

In the period of 1994–2016, before the SCSSM withdrawal pentad, strong westerly anomalies are found over the Bay of Bengal extending through the South China Sea to the east of the Philippines, which transports abundant warm moisture to favor the active convection over the SCS and Northwest Pacific (Figure 8(e)). On the withdrawal pentad, the low-level wind abruptly transits to a northeasterly anomaly in the SCS and persists into the next pentad, which is stronger than that in the period of 1979–1993. On the other hand, a suppressed convection appears in the SCS and eventually extends eastward to the Bay of Bengal on the next pentad of the SCSSM withdrawal (Figures 8(f) and 8(g)). Comparing the differences before and after the SCSSM withdrawal pentad, the increased easterly wind in the SCS is sharper than that in the period of 1979–1993, and the convection over the Northwest Pacific, SCS, and the Bay of Bengal is weakened along with the SCSSM withdrawal, which is associated with a southward migration of the seasonal march of the ITCZ (Figure 8(h)). It is suggested that the SCSSM withdrawal rate during the period of 1994–2016 is faster than that of the period of 1979–1993.

In short, the distinct characteristics of the SCSSM are that, during the period of 1979–1993 (1994–2016), the SCSSM onset is generally later (earlier) and the intensity of the SCSSM onset is stronger (weaker). On the other hand, there is no obvious difference in the SCSSM retreat date, but the speed of the SCSSM retreat is slower (faster) and the duration is longer (shorter). This raises the following question: what is the effect of this change on rainfall in China?

3.3. Distinct Rainfall Distribution from the Onset to Withdrawal of the SCSSM during the Two Periods

Many previous studies have confirmed that the evolution of the SCSSM has a broad impact on the spatial and temporal distribution of rainfall in East Asia during the period from May to September. Particularly, the SCCSSM onset has been regarded as the beginning of the East Asian rainy season [12, 3841]. It is necessary to figure out the link between the different characteristics of the SCSSM in two epochs and the rainfall distribution in Eastern China. Seasonal evolution of latitudinal (110°–120°E) averaged rainfall between two epochs is shown in Figure 9. In the late-onset years of 1979–1993, the maximum rainfall belt (>10 mm·day−1), corresponding to the first rainy season in Southern China, begins at the end of May (Figure 9(a)). Moreover, for the late-onset period, the SCSSM intensity is stronger, leading to more warm moisture being constantly transported into Southern China. The rainfall persists consistently until the beginning of June until the major rain-belt moves northward to the middle and lower reaches of the Yangtze River basin (around 30°N), which is associated with the Meiyu front. By the end of July, the maximum rainfall belt moves over Southern China again, which can be regarded as the beginning of the last flood season in Southern China. The rainfall in Eastern China reduces from 35°N north to 25°N south from early September to early October.

In the period of 1994–2016, the SCSSM onset is generally early. The first rainy season in South China begins at the beginning of May and continues intermittently to the end of June until the major rain-belt moves northward to the middle and lower reaches of the Yangtze River basin (Figure 9(b)). The Meiyu period during 1994–2016 is shorter than that during 1979–1993. By mid-July, the maximum rainfall belt moves to the Southern China again and it decreases with Eastern China rainfall simultaneously from late September to early October, consistent with the faster and shorter retreat of the SCSSM. On the other hand, the two migrations of the rain-belt are closely related to the two-peak pattern of the zonal wind over the SCS in summer, which is consistent with previous studies [42].

4. Effects of ISO on the SCSSM Evolution from Onset to Withdrawal during the Two Periods

The above analyses have shown different characteristics of the SCSSM from onset to withdrawal before and after 1993/94 and the effect of the SCSSM change on rainfall distribution. Thus, a question to be further addressed is, what are the mechanisms responsible for this characteristic change? Several studies have discussed the importance of two dominant modes of ISOs on the triggers of the SCSSM onset and withdrawal [4346]. It is found that 10–25-day and 30–80-day oscillations also control the SCSSM activity from May when the SCSSM onsets to October when the SCSSM retreats through power spectrum analysis (fig. not shown). To detect the effects of ISO on the SCSSM evolution from onset to withdrawal during the two epochs, we calculated standard deviations of the 10–25- and 30–80-day filtered OLR from 15 May to 15 October, which cover the 5 months from climatological SCSSM onset to withdrawal. Figure 10 shows the differences (1994–2016 to 1979–1993) of mean standard deviation between the two epochs. The difference of the 30–80-day oscillation mode over the SCS is significant and the strength of the 30–80-day oscillations is large in the SCS during the period of 1994–2016. Thus, the effect of 30–80-day ISO on the SCSSM evolution during the two periods will be studied further.

To explore the 30–80-day ISO activities of the SCSSM, an East Asian summer monsoon index (IM), which is proposed by Mao and Chan [44], has been applied. That is, , where , R is OLR, and are the standard deviations of and OLR, respectively, and and are the multiyear averages of and OLR, respectively. The advantage of this index is that it takes into account the dynamic and thermodynamic characteristics of the monsoon. In this study, the 30–80-day band-pass filtering of IM over the SCS (110°–120°E, 10°–20°N) is applied as the daily ISO index of the SCSSM. Figure 11 shows the typical daily ISO evolution from April to November in the SCS between two epochs. A strong modulation of SCS ISOs on the SCSSM evolution is observed. In the late-onset period of 1979–1993, the SCS ISO has been in a dry phase, postponing the monsoon onset due to low-frequency easterly wind anomalies and weakening convection until late May. The ISO then turns to the positive phase and develops to the maximum amplitude of the whole monsoon period. The low-frequency convection and westerly wind are significantly strengthened and favor a strong SCSSM onset. The transformation time is very near pentad 29 (May 26–31). There are two active low-frequency oscillations in summer, namely, “South China Sea monsoon surge,” which has been discussed in depth by Li et al. [46], corresponding to the two rain-belt migrations. After September, the ISO gradually weakened and disappeared, corresponding to the weakening and slowing retreat of the SCSSM.

In the early-onset years of 1994–2016, the ISO of the SCSSM turns to the positive phase around the beginning of May. The first positive peak value of ISO is smaller than that in the period of 1979–1993, corresponding to a weaker SCSSM onset intensity. Then there are three active low-frequency oscillations from June to September, with stable amplitude. In early October, the ISO transforms rapidly to the negative phase, contributing to the shorter, faster retreat of the SCSSM. It is inferred that the 30–80-day ISO was crucial for the SCSSM from onset to withdrawal between two epochs.

We have made a preliminary analysis for the effect of local SCS 30–80-day ISO on the SCSSM evolution difference; how does each ISO propagate and where is the ISO propagation source? Figure 12 shows the zonal and meridional propagation of ISO along with 10°–20°N and 110°–120°E during the two periods. In the period of 1979–1993, the wet ISO propagates eastward from tropical Indian Ocean to the SCS along with 10°–20°N in early June and mid-August (Figure 12(c)). Moreover, when ISO activity in the Indian Ocean is relatively strong, the propagation to the SCS is also strong. In meridional direction, the positive ISO propagates northward into the SCS. In September, there is ISO propagating from Indian Ocean but not entering the SCS. Different from the late-onset period, in the period of 1994–2016, the first two wet ISOs propagate westward from the Northwest Pacific in May and June, and the last two propagate eastward from the Indian Ocean in July and September. In zonal direction, the ISOs also propagate northward into the SCS. Therefore, the wet ISO in the Northwest Pacific propagating northwestward into the SCS may contribute to the SCSSM evolution difference after 1993/94 and the SCSSM ISO intensity can be predicted by monitoring ISO anomalies over the northwest Pacific and Indian oceans.

5. Effects of SSTAs on the SCSSM Evolution from Onset to Withdrawal during the Two Periods

In addition to the effects of internal atmospheric dynamic mechanism, SST, an important external force, also plays an important role in the SCSSM evolution. Previous studies have discovered that the delayed onset of SCSSM may result from a La Niña event [47]. The early onset and delayed retreat of the SCSSM can be ascribed to warm SST anomalies in the Philippine Sea [23, 30]. A warm SST in the tropical Indian Ocean during the previous winter [19] can delay the SCSSM onset by suppressing the convection over the western Pacific. It is also interesting to find that the spatial pattern of the SST evolution exhibits specific interdecadal characteristics during two periods. As shown in Figure 13 on the left, significant positive SST anomalies appear around the equatorial Eastern Pacific in the late-onset years of the period of 1979–1993 in May and continue developing until November. Meanwhile, the spatial pattern of the SST changes is similar to the horseshoe-shaped SST warming during El Niño events. In theory, the warm tropical eastern Pacific SST anomaly may induce an anticyclone anomaly in the northwest Pacific due to the excitation of ascending Rossby waves. This anticyclone anomaly connects to the low-level easterly anomaly in the SCS, which possibly postpones the onset of SCSSM [48]. In the early-onset years of the period of 1994–2016, significant warm SST anomalies appear in the Philippine Sea from April to September, and the strongest warm SST appears in May and then gradually weakens. A few weak negative SST anomalies appear around the equatorial Eastern Pacific (Figure 13 on the right). Thus, the SST warming in the Philippine Sea from April to September may be a crucial external force for the SCSSM evolution in the period of 1994–2016. However, a deeper understanding of the mechanism of this SST change on the SCSSM evolution during two epochs calls for further study.

6. Summary and Discussion

In this study, we investigate the different characteristics and possible impact factors of the SCSSM evolution from onset to withdrawal before and after 1993/94 when an interdecadal change occurs. The results suggest that, during the period of 1979–1993, the SCSSM onsets later, corresponding to abrupt enhancement of zonal wind and convection activity, results in a stronger SCSSM onset intensity. Delayed onset also results in a continuous, strong rain-belt appearing in Southern China at the end of May. After the SCSSM onset, the low-level zonal winds appear with a two-peak pattern in summer, and two abrupt increases of the zonal winds correspond to the beginning of the Meiyu period in Yangtze River valley and the last flood season in Southern China, respectively. The monsoon retreats from September until October and exhibits an extended, gradual retreat, which causes the withdrawal of the rain-belt to begin in September from north to south China and completely vacate at the beginning of October.

However, after 1993/94, the SCSSM onsets earlier and is characterized by a weaker onset intensity and more rapid retreat. The preflood season advanced and lasted intermittently throughout May. The rain-belt in Southern China persisted for a longer time in summer. The monsoon started to retreat from October, resulting in the rapid retreat of the rain-belt from the whole eastern China.

Further analyses indicate that the enhanced 30–80-day SCS ISO activities are an important contributor for SCSSM evolution change. There are two active low-frequency oscillations over the SCS in summer during the late-onset period but three ISOs in the early-onset period. Prior to 1993/94, the SCS ISO had been in a dry phase, which postponed the monsoon onset until late May. Then a late and strong onset is triggered by enhanced wet ISO propagating northeastward from the equatorial Indian Ocean. After September, the ISO gradually weakened and disappeared, corresponding to the weakening and slowing retreat of the SCSSM. After 1993/94, the wet ISO from the Northwest Pacific propagated northwestward to the SCS in early May, which modulated the mid-May onset of SCSSM. After October, the SCS ISO transformed rapidly to the negative phase, contributing to the shorter duration and faster retreat of the SCSSM.

Moreover, the distinct SCSSM evolution is closely related to the SST change. Significantly positive SST anomalies appear around the equatorial Eastern Pacific in the late-onset years of the period of 1979–1993 in May and continue developing until November. The spatial patterns of the SST changes are similar to the horseshoe-shaped SST warming during El Niño events. However, in the early-onset years of the period of 1994–2016, warm SST anomalies appear in the Philippine Sea from April to September. Many researchers have noted the effect of ENSO-related western North Pacific SST and Philippine Sea temperature on the SCSSM onset in May and withdrawal in September [30, 32, 49]. The specific physical mechanism for the impact of SST on the monsoon evolution from May to September should thus be studied further. On the other hand, air-sea interaction can provide an additional northward propagation mechanism for 30–60-day Northwest Pacific ISO and ISO source [50]. The influence of SSTAs on the SCSSM ISO should also be studied as it may provide a basis for predicting the subsequent SCSSM evolution.

Data Availability

Daily precipitation data and OLR were derived from https://www.esrl.noaa.gov/psd/data/gridded/. The HadISST data were from http://www.metoffice.gov.uk/hadobs/hadisst/data/download.html. The ERA-Interim data were obtained from https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=pl/.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was jointly supported by the National Key R&D Program of China (2019YFC1510004) and the National Natural Science Foundation of China (41975085).