3.1 Climatology of ISMR and SST
The climatological features of rainfall in the Indian subcontinent and SST in the Indo-Pacific domain during the summer monsoon season (June to September) for the period of 1901 to 2020 are given in Fig. 1. Figure 1a depicts high rainfall over the windward side of the Western Ghats and the northeastern regions of the country, where the monsoon rainfall exceeds 16 mm day− 1. Over central India, the rainfall is about 10 mm day− 1 and it is less than 5 mm day− 1 over southeast, north and northwest Indian regions. The Western Ghats in the west coast region receives high rainfall, with an average annual rainfall of about 3800 mm year− 1 in the windward side (Rao., 1976). North eastern region also gets high amounts of rainfall compared to the north and central part of India. Rain shadow region of south India receives less rainfall during the monsoon season.
In Fig. 1b, the equatorial Pacific Ocean experiences high SST greater than 27o C, which is concentrated over the western side between 120o E and 160o E. Eastern side experiences relatively low SST between 21o C and 25o C. SST decreases towards higher latitudes beyond 20o N. The equatorial Indian Ocean shows SST more than 27o C. The PDO region is comparatively cooler than the Niño region. The spatial distribution of rainfall and SST varies considerably during the occurrence of ENSO and PDO, moreover the co-occurrence of both the events can amplify the modification in their spatial pattern. Hence a detailed understanding of the variability in different space time domains and relationship of PDO and ENSO with ISMR is needed.
3.2 Interannual variation and relationship of PDO and ENSO
Interannual variations between PDO and ENSO along with linear trend for the period before and after 1976 (1930–1975 and 1976–2020) for a period from 1901 to 2020 during southwest monsoon are given in Fig. 2a. PDO and ENSO interannual variations have witnessed major alterations before and after the 1970s. The standard deviation of Niño 3.4 index is 0.60 and PDO index is 0.92, both are significant at more than 99% confidence level. In the decadal variations, before the 1970s (pre-shift), the phase transitions of PDO took 15–20 years but after the 1970s (post-shift), only 12–15 years were taken for it. In recent years after the 1970s, the warm (cold) phase of PDO tends to occur with the warm (cold) phase of ENSO. Earlier years showed a less intense ENSO and PDO where magnitude of indices were in between − 1.5 and + 1.5. But recent years have phases with high intensity which are varying in between − 2.5 to + 2.5. During 1930 to 1975, interannual variations of PDO index and Niño 3.4 index are inhomogeneous. Occurrence of El Niño (La Niña) is not coinciding with positive PDO (negative PDO). Intensity of ENSO is smaller ( in between − 1 to + 1) than that of PDO (in between − 2 to + 2). During 1976 to 2020, the occurrence of El Niño (La Niña) coincided with the occurrence of positive PDO (negative PDO). Both indices have similar intensities (ranging from − 3 to + 3) and interannual fluctuations. These observable changes indicate that a climate shift occurred in 1976, which is consistent with previous studies (Graham 1994; Chaluvadi et al., 2021, Miller et al., 1994). It should be noticed that the ENSO trend has changed from negative to positive (trends are not significant), however the PDO trend has shown significant changes in their trend values (from − 0.29 to -0.01). Global warming and abrupt changes in SST conditions can be a reason for the changes occurring in the prevailing climatic conditions (Philander and Philander, 2008; Schuldt et al., 2011).
Figure 2b is showing the regression analysis of the Niño 3.4 index and PDO index during the southwest monsoon from 1901 to 2020. In the Fig. 2b, the first quadrant represents (warm) positive PDO (pPDO) and El Niño phases (positive values for PDO index and Niño 3.4 index), second quadrant represents (cold) negative PDO (nPDO) and El Niño phases (negative values for PDO index and positive values for Niño 3.4 index), third quadrant represents nPDO and La Niña (negative values for PDO index and Niño 3.4 index) and fourth quadrant represents pPDO and La Niña (positive values for PDO index and negative values for Niño 3.4 index). There exist years having warm (cold) phases of PDO and ENSO. Mixed cases are also seen but the number of occurrences is less. Regression line is showing a positive relationship. Points are staggered and the linear correlation coefficient is 0.5 and regression coefficient of linear relation is 0.3 which denotes an in-phase relationship. Even though they are having one to one relation, lots of staggered points are there which gives an importance to study its properties on quadrant basis.
3.3 Combined effect of PDO and ENSO on SST and ISMR anomalies
To substantiate the relation between ISMR and combination of PDO and ENSO, composites of summer monsoon seasonal anomalies of rainfall and SST are shown in Fig. 3 and Fig. 4, respectively. The number of pPDO and El Niño, pPDO and La Niña, nPDO and El Niño and nPDO and La Niña years included in the composite analysis are 32, 27, 21 and 39, respectively. The nPDO and La Niña composite of SST is shown in Fig. 3a indicating negative SST anomalies in the eastern part along the west coast of North America and stronger positive SST anomalies in the central and western parts of the north Pacific. This gives strong nPDO signatures. A pattern of severe negative anomalies exists in the eastern and central equatorial Pacific, with a neutral situation in the western Pacific indicating La Niña signatures. The pPDO and El Niño composite (Fig. 3c) shows a pattern similar to the nPDO and La Niña composite but with SST anomalies of opposite signs. Figure 3c shows positive SST anomalies in the eastern Pacific and negative SST anomalies in the North western and central Pacific. The equatorial Pacific is having positive anomalies of SST in the eastern side and weak negative anomalies in the western side. Both of the above mentioned in-phase examples, the Indian Ocean appears to be a mirror image of the equatorial Pacific Ocean.
Figure 3b shows the SST anomalies of nPDO and El Niño. It has positive anomalies in the equatorial eastern Pacific and weak negative anomalies in the western Pacific, consistent with El Niño conditions. Positive anomalies in the north Pacific with slightly weak negative anomalies of SST consistent with the nPDO phase. However, the anomalies are weak compared to the in-phase composites. Thus, the nPDO phase and El Niño counteract each other to minimise their proper signals. Similar to the case of pPDO and La Niña composite in Fig. 3d, which shows negative anomalies in the equatorial eastern Pacific as well as northwest Pacific and weak positive anomalies in northeast Pacific. In the case of out-of-phase examples mentioned above, the Indian Ocean is a mirror image of the equatorial Pacific, just as in in-phase ones, but with different spatial distribution and homogeneity.
The previously described SST features influence southwest monsoon rainfall patterns. Spatial distribution of rainfall anomalies of southwest monsoon during different combinations of PDO and ENSO in Fig. 4 is used to analyse this concept. The negative PDO and La Niña composite of rainfall shows positive anomalies of rainfall all over India except the North eastern region (Fig. 4a). Higher anomalies (> 3 mm day− 1) are observed in the northern portion of the Konkan coast. This enhanced flood condition over almost the entire India indicates that La Niña and nPDO, both of which produce flood conditions over India, are complementing each other. pPDO and El Niño composite (Fig. 4c) shows below normal rainfall patterns over central India and western coast of the country (Krishnan et al., 2003; Kanamitsu and Krishnamurti, 1978; Krishnamurthy and Goswami, 2000). However the Northeastern region is showing positive anomalies of rainfall. This enhanced signature in the rainfall anomalies seems to be the result of pPDO and El Niño, both associated with negative rainfall anomalies over India, complementing each other. Figure 4b shows nPDO and El Niño composite of rainfall in which India is experiencing a combination of drought and neutral conditions, with a domination of drought conditions. In a pPDO and La Niña composite (Fig. 4d), floods occur in the north and central regions of India, whereas negative rainfall anomalies occur in the southern India and western regions of Gujarat. Hence out-of-phase cases having a mixed pattern of rainfall with ENSO domination.
3.4 Comparison of SST and rainfall using pure conditions
To further substantiate the relation between monsoon and ENSO and to understand such relation relative to PDO, composites of summer monsoon seasonal anomaly of rainfall and SST of pure ENSO and PDO cases are shown in Fig. 5 and Fig. 6, respectively. Figure 5a illustrates La Niña conditions with a negative SST anomaly in the eastern and central equatorial Pacific and weak positive anomalies in the western equatorial Pacific. The Indian Ocean is likewise witnessing a negative SST anomaly, particularly in the Arabian Sea. Rainfall composite during La Niña period (Fig. 6a) exhibits an above-normal precipitation pattern across India, with the exception of the central and eastern parts, which show below-normal precipitation. El Niño is causing a horseshoe pattern of positive SST anomalies in the central and eastern equatorial Pacific and negative anomalies in the western equatorial Pacific (Fig. 5b). The positive SST anomaly can also be detected in the western Indian Ocean. Rainfall composite during El Niño (Fig. 6b) depicts the drought situation in India. A positive SST anomaly can be evident in the northeastern Pacific during a warm or pPDO (Fig. 5c), but a significant negative SST anomaly can be detected in the central and western Pacific. A positive SST anomaly can also be found across the Indian Ocean. Rainfall composite of pPDO (Fig. 6c) indicating a minor precipitation deficit. A negative SST anomaly exists in the northeastern Pacific during nPDO (Fig. 5d), while positive SST anomalies exist in the central and western Pacific. The Indian Ocean is experiencing neutral SST conditions which leads to the above normal rainfall pattern in the rainfall composite of nPDO (Fig. 6d).
Drought conditions are experienced when the pPDO and El Niño is present. Drought severity is greater in El Niño situations than in pPDO. Flooding is the result of both a nPDO and La Niña. As noted previously, the La Niña situation causes more positive rainfall anomalies than negative PDO situations. It can be concluded that the ISMR has a negative correlation with PDO and ENSO. Spatial correlation of ISMR with PDO and Niño 3.4 indices can further substantiate this result. Figure 7 is the spatial correlation coefficient of Indian summer monsoon rainfall with PDO and Niño 3.4 indices during southwest monsoon. Contours represent regions with 95% significance level.
The two indices show regionally varying relationships with ISMR. The correlation is significant over most parts of India, except in the northeast, where the correlation is insignificantly out-of-phase as noted in earlier studies (e.g., Guhathakurta and Rajeevan, 2008; Varikoden and Babu 2015; Nair et al., 2018; Shukla and Mooley 1987). The correlation values are more significant in the Niño 3.4 index than PDO index. In some parts of India correlation of ISMR with PDO index is becoming positive. The ENSO (Fig. 8b) is having a strong negative correlation with ISMR compared to PDO. Most of the regions are having a significance of more than 95% confidence level. This can be identified using contours drawn in the figure. When comparing years with a combination of phases and pure ones, it is noted that if they are in-phase, strong and well defined SST features of both events occurring in the site. When they are out-of-phase, their intensity drops and various places experience neutral conditions. The horseshoe pattern of SST, which is an important feature of the ENSO SST pattern, is well characterised in in-phase conditions but not in out-of-phase conditions. Rainfall intensifies (weakens) when both are in warm (cold) phase, ie, during in-phase conditions. In the out-of-phase scenario, ENSO patterns are more prominent than PDO patterns.
3.5 Climate shift and change in the behaviour of teleconnections
Climate shifts in the Pacific Ocean have a significant impact on the phenomena that occur there (meehl et al., 2009; Mantua and Hare, 2002). During the shift basic characteristics and properties of all the phenomena in the Pacific Ocean changed (Graham, 1994). Considering the changes in the Pacific ocean, It is assumed that phenomena like ENSO and PDO, which are produced by the SST and SLP variations in the Pacific Ocean can also be changed along with this.
To recognize changes in the SST anomalies of the Pacific Ocean under distinct ENSO and PDO phase combinations, spatial distribution of SST (oC) anomalies during PDO and ENSO were drawn for pre-shift and post-shift (Fig. 8). Negative SST anomalies can be detected in the northwestern and central Pacific during pPDO and El Niño situations (Fig. 8a), with positive anomalies near the western coast of north America. Positive anomalies can also be seen in the central and eastern equatorial Pacific, while negative anomalies can be seen in the western equatorial Pacific. Positive anomalies detected around the western coast of North America, as well as in the eastern and central Pacific. However, the same relationship exists in the post-shift period with higher magnitudes (Fig. 8e). Negative anomalies in the western Pacific region are diminishing. The Indian Ocean has negative anomalies and neutral temperatures before the shift, but positive anomalies after the shift.
Negative SST anomalies can be detected in the northwestern and central Pacific during positive PDO and El Niño situations (Fig. 8b), with positive anomalies near the western coast of North America. Positive anomalies can also be seen in the central and eastern equatorial Pacific, while negative anomalies can be seen in the western equatorial Pacific. Positive anomalies are found around the western coast of North America, as well as the eastern and central Pacific, become stronger after the shift (Fig. 8f). Negative anomalies in the western Pacific region are diminishing. The Indian Ocean has negative anomalies and neutral temperatures in the pre-shift condition, however, it turned to positive anomalies during the post-shift period.
In nPDO and El Niño conditions (Fig. 8c), positive SST anomalies in the north Pacific with a weak negative anomaly around it show nPDO characteristics. The horseshoe pattern of positive SST anomalies across central and eastern Pacific with weak negative anomalies at the western Pacific shows El Niño signature. However, the anomalies are weak compared to the in-phase composites. Thus, the nPDO phase and El Niño counteract each other to minimise their proper signals. El Niño characters are visible compared to PDO. The Indian Ocean is having a negative SST anomaly. In the post-shift phase (Fig. 8g), a stronger positive anomaly of SST is observed in the central and western Pacific, with a weaker negative anomaly along north America's western coast. Positive anomaly is present in the equatorial Pacific, reducing El Niño signatures which depicts clear dominance of PDO characteristics. Figure 8d, which shows nPDO and La Niña, has negative SST anomalies over the western coast of North America and positive SST anomalies along the western and central Pacific, both are having strong PDO signatures. Negative SST over the equatorial Pacific is a strong indicator of La Niña. SST traces of both events are observed. Positive anomalies are evident everywhere over the Oceanic region after the shift (Fig. 8h), with the exception of the eastern side of the Pacific, which has weak negative SST anomalies. Both phenomena's signature traits are missing, despite the increased positive anomaly of SST. The Indian Ocean is having negative SST anomalies in pre-shift however positive anomalies in post-shift.
To investigate variations of ISMR and to identify differences in ISMR pattern before and after climate shift, spatial distribution of rainfall (mm day− 1) anomalies during PDO and ENSO periods before climate shift and after climate shift is shown in Fig. 9. In the positive PDO and El Niño case, the pre-shift condition (Fig. 9a) is showing a negative anomaly of rainfall with some positive rainfall anomaly western coast and northern parts of India. whereas the aftershift (Fig. 9e) is showing a more intensive negative anomaly of rainfall. pPDO and La Niña case is having positive rainfall anomalies favourable for La Niña condition in the pre-shift period (Fig. 9b), whereas negative rainfall anomalies which is a favourable condition for positive PDO are identified in post-shift (Fig. 9f). During nPDO and El Niño, negative anomalies of rainfall (El Niño favourable situation) can be seen in the pre-shift case (Fig. 9c), whereas positive anomalies (nPDO favourable situation) can be seen in the after shift phase (Fig. 9g). nPDO and La Niña show an increased amount of rainfall anomaly in pre-shift (Fig. 9d) but during post-shift (Fig. 9h), positive rainfall anomaly weakens and can see some negative rainfall anomalies in north eastern and central portions of India. Thus the analysis shows that during in-phase conditions, SST anomalies for ENSO and PDO combinations are showing their individual SST signatures. Warm (cold) phase has characteristics of El Niño (La Niña) and pPDO (nPDO) in pre-shift. In post-shift, negative (positive) anomalies become weaker (stronger). Rainfall decreases (increases) when warm (cold) phases coexist during pre-shift. After the climate shift, negative (positive) anomalies of rainfall become stronger (weaker). Out-of-phase conditions counteract each other and minimise their characteristic SST signals. But ENSO features are prominent during the pre-shift and PDO features are prominent during after shift. Rainfall patterns are coinciding with this result.
3.6 Enhancement of changes using dynamical parameters
Understanding the low level circulation anomalies helps us to identify the regions where the circulation pattern is favourable for rainfall. Figure 10 shows circulation anomalies at 850 hPa during the pre-shift and the post-shift periods. In pPDO and El Niño conditions, the pre-shift case (Fig. 10a) shows the westerly component of easterly trade winds strengthens over the central Pacific region. This strengthening of westerly components modifies the ENSO - monsoon relationship. Anomalous westerlies are also observed at 20°N in the Pacific Ocean, which indicates the weakening trade winds due to the high SST/ low SLP produced by pPDO and El Nino conditions. This combined effect decreases rainfall in India. The divergence pattern of the wind vector which leads to descending motion can also be attributed to decreased rainfall. Figure 10e shows the post-shift condition of pPDO and El Niño, decreased intensity of anomalous westerlies is seen in both Niño and PDO regions. The divergence pattern of wind vectors is prominent over India which further decreases rainfall than that of the pre-shift condition. The pPDO and La Niña during the pre-shift situation (Fig. 10b), show a convergence of the wind, thereby enhanced rainfall over India. Easterlies are prominent in the Niño region and the Indian Ocean region. Weakened easterlies originated from the western side of the Niño region due to the lower SST and higher SLP due to the climatic conditions caused by the La Niña. In the post-shift case (Fig. 10f), there is divergence all over central India, and strong easterlies are evident, decreasing rainfall. The rainfall is further weakened by westerlies across the Indian Ocean.
The pre-shift (Fig. 10c) has weak anomalous westerlies over the Indian region, which strengthen the prevalent trade winds over the Indian subcontinent that leads to deficit rainfall during nPDO and El Niño conditions. During the post-shift period (Fig. 10g), strong westerlies led to convergence over the Indian regions, enhancing convection and thereby enhancing rainfall over the Indian subcontinent. In nPDO and La Niña conditions (Fig. 10d), a convergence is observed over India in association with very weak easterlies in the western Pacific. This convergence weakens during the post-shift situation (Fig. 10h), which weakens the rainfall. Strong easterlies are observed in the western central Pacific region along with strong westerlies in the central eastern and north Pacific regions. Analysis of wind vectors are consistent with the results obtained from the analysis of SST and rainfall.
Further analysis of velocity potential and stream function is done to understand large scale features of horizontal circulation at 850 hPa levels (Fig. 11). During the pre-shift case (Fig. 11a), a tendency of divergence across India is noticed that reduces rainfall over there. The PDO region is also showing higher convergence and thereby ascending motion of air. In the post-shift period (Fig. 11e), a positive anomaly of velocity potential intensified, increasing the descending motion through an increase of divergence over India and decreasing rainfall. In the pre-shift phase of pPDO and La Niña (Fig. 11b), intensified convergence over India increases the amount of rainfall. In the post-shift (Fig. 11f) phase, there is a divergence, which reduces rainfall. Intensified divergence patterns are observed in between 20o N and 40o N in the post-shift period due to the increase in the temperature at that region after the shift.
During the nPDO and El Niño, the divergence pattern caused by the relocation of the convergence centre towards the Indian Ocean resulted in deficit rainfall over India in the pre-shift case (Fig. 11c). The post-shift case (Fig. 11g) shows a small convergence over north India, causing an above normal rainfall. Increased convergence near the Indian subcontinent enhances positive rainfall anomalies in the pre-shift situation of nPDO and La Niña (Fig. 11d). But the post-shift situation (Fig. 11h) shows a slightly divergent pattern over India, which weakens the positive rainfall anomaly. Thus in-phase conditions show an increased positive (negative) rainfall anomalies when warm (cold) phases of PDO and ENSO coexist. In the case of out-of-phase conditions, the rainfall pattern shows similar to the rainfall pattern of ENSO during the pre-shift case; however, the rainfall pattern shows as that of the PDO pattern during the post-shift period.
In order to compare large scale circulation features at lower level (850 hPa) between pre- and post-shift periods, anomalies of stream function are given in Fig. 12. During pPDO and El Niño during the pre-shift (Fig. 12a), Above the Indian subcontinent, a positive anomaly in stream function was seen, resulting in anticyclonic circulation and descending motion. As a result, the amount of rainfall falls. In the post-shift (Fig. 12e) case, positive anomalies shift towards the Arabian sea which inhibits convection hence further decreases the rainfall. In the pre-shift case (Fig. 12b), pPDO and La Niña shows a feeble negative anomaly of stream function over India, which increases rainfall as per property of La Niña but post-shift (Fig. 12f) shows just the opposite, which decreases the rainfall pattern as in pPDO. In the nPDO and El Niño pre-shift case (Fig. 12c), a positive anomaly of stream function causes the rainfall to decrease as in the case of pure El Niño. These positive anomalies result in the descending motion and anticyclonic circulation over India, which further suppresses the rainfall. But after shift (Fig. 12g) shows a negative anomaly of stream function over India with increased rainfall by the counterclockwise circulation pattern and ascending motion of air as in the pure PDO rainfall pattern. The nPDO and La Niña (Fig. 12d) case is showing negative anomalies with enhanced rainfall in the pre-shift but post-shift shows a weak positive anomaly (Fig. 12h), and it can be attributed to deficit rainfall. In general, the results of the previous analysis are also consistent with SST and rainfall observations, as well as dynamical characteristics like horizontal wind velocity and velocity potential.