Basic characteristics of JSN. As shown in Fig. 1a, JS Niña often develops as early as April-June (AMJ) along the south coast of Java, peaks in July-September (JAS) as the cold SSTAs slowly shift to the southwest coast of Sumatra, and dissipates in October-December (OND). Although weak, there is a secondary peak of JS Niña in the preceding January-March (JFM; hereafter, pre-season). This pre-season cooling pattern is closely tied to a basin-wide cooling of the tropical Indian Ocean that is largely modulated by the Indian Ocean Basin (IOB) mode, which tends to occur during the dissipation phase of La Niña in boreal winter and spring35, as discussed later in this study. Interestingly, relatively weak warm SSTAs develop in the central tropical South Indian Ocean in JAS, and then propagate westward to produce a positive IOD-like pattern in OND.
During the most active months (JAS), JS Niña is typically characterized by cold SSTAs, a shallower thermocline, southeasterly surface wind anomalies along the Java-Sumatra coast, a deeper thermocline farther westward in the tropical South Indian Ocean, and easterly surface wind anomalies along the eastern equatorial Indian Ocean (Fig. 1b). These anomalous atmosphere-ocean patterns indicate that local Bjerknes feedback plays an important role in the development of JS Niña23. More specifically, the cold SSTAs off the Java-Sumatra coast suppress atmospheric convection aloft. The associated heat sink induces an anomalous anticyclone to its west (i.e., a descending Rossby wave response34), producing southeasterly wind anomalies and anomalous upwelling off the Java-Sumatra coast, which in turn reinforce the initial cooling23. The heat sink also drives equatorial easterly wind anomalies (i.e., an equatorial divergent wind response34), thus producing anomalous upwelling in the eastern equatorial Indian Ocean and promoting an expansion of the cold SSTAs toward the equator along the Sumatra coast (Fig. 2).
JS Niño is characterized by nearly opposite patterns of anomalous SSTs, thermocline depth and surface winds compared to JS Niña (Fig. 1d,e). However, there are several noticeable differences that break the antisymmetry. JS Niño tends to develop and decay slightly earlier than JS Niña in boreal summer and fall, and the associated atmosphere-ocean anomalies are generally much weaker36. In contrast, the pre-season peak in JFM, which is tied to a positive phase of the IOB mode and the dissipation phase of El Niño35, has a much larger amplitude, being on par with the peak amplitude in June-August. Additionally, unlike JS Niña, JS Niño tends to remain a monopole throughout boreal summer and fall.
JS Niña events are linked to decreased rainfall over the Maritime Continent (i.e., maritime Southeast Asia) and Malay Peninsula, and vice versa for JS Niño events (Fig. 1c,f). Some JS Niña/Niño events are also linked to increased rainfall in South Asia and mainland Southeast Asia (hereafter, South-Southeast Asia). However, these and other basic characteristics of JSN derived from the composite mean fields (Fig. 1a-f) are not applicable to all events. Indeed, the spatiotemporal evolution patterns are quite different between individual events in terms of the timing, zonal pattern, and amplitude of their onset, peak, and dissipation (Supplementary Figs. 1 and 2).
Spatiotemporal diversity of JSN. As described in Methods, we carry out a spatiotemporal EOF analysis to identify the first two preferred modes of the observed JS Niño/Niña events. Each EOF mode for warm/cold events represents two contrasting JS Niño/Niña varieties or flavors (PC changing from − 1 to 1) that correspond to adding and subtracting the EOF pattern to the composite mean, leading to four sets of the most frequently recurring JS Niño/Niña varieties. See Methods and Supplementary Fig. 3 for more details.
One of the leading JS Niño/Niña varieties describes warm/cold events that form in JAS and thus is referred to as the on-season JS Niño/Niña variety (Fig. 3f,h). The 1989 warm event and the 1994 cold event are good examples of the on-season JS Niño/Nina variety (Fig. 3i,j; see also Supplementary Figs. 1 and 2). Other JS Niño/Niña varieties describe warm/cold events that appear as early as JFM and then, after a short pause, reappear in JAS (hereafter, the double-peak JS Niño/Niña variety; Fig. 3e,g) or almost completely dissipate before boreal summer (hereafter, the pre-season JS Niño/Niña variety; Fig. 3a,c). The 2016 warm event is an example of the former, while the 1974 cold event and the 1983 warm event are examples of the latter (see also Supplementary Figs. 1 and 2).
Occasionally, strong cold events may develop in boreal summer, persist throughout SON and spawn robust warm SSTAs in the central tropical South Indian Ocean, which in turn propagate westward and peak in OND. The spatiotemporal evolution of such strong and persistent cold events (e.g., the 2019 JS Niña) is well captured by the persistent JS Niña variety (Fig. 3b). Finally, although rare and very weak in amplitude, a pre-season warm event may switch to a cold event in the boreal summer (hereafter, the transition JS Niño variety; Fig. 3d). It is worth noting that the persistent JS Niña variety displays a positive IOD-like pattern in OND, whereas the double-peak JS Niño variety prevails largely as a monopole throughout boreal summer and fall. Both the double-peak JS Niña and on-season JS Niño varieties display an IOD-like pattern in boreal summer, while the on-season JS Niña variety has a more complex zonal tripole structure.
Interactions with ENSO. The pre-season and double-peak JS Niño/Niña varieties are strongly connected to the dissipating phase of El Niño/La Niña in boreal spring (Fig. 4a,c,e,g). Specifically, the pre-season JS Niño/Niña variety is favored if El Niño/La Niña dissipates into an ENSO-neutral phase in boreal summer, whereas the double-peak JS Niño/Niña variety is favored if El Niño/La Niña transitions to the opposite phase and then intensifies during boreal summer and fall. However, some double-peak JS Niño/Niña events occurred under a near-neutral ENSO condition in JAS (e.g., the 2016 JS Niño; Supplementary Figs. 4 and 5). In contrast, there is no apparent ENSO signal linked to the on-season JS Niño/Niña variety in boreal summer and fall. This suggests that the on-season JS Niño/Niña variety is predominantly driven by atmosphere-ocean processes internal to the Indian Ocean. Similarly, the transition JS Niño variety seems to have little connection to ENSO.
Perhaps the most interesting case is the persistent JS Niña variety, which is characterized by the greatest intensity of cold SSTAs off the Java-Sumatra coast. During these events, strong warm SSTAs develop in the southwestern equatorial Indian Ocean around OND, so this variety can also be described as strong positive IOD. The persistent JS Niña variety is association with the developing phase of El Niño, but the regressed ENSO signal is relatively weak, especially during the onset phase of the persistent JS Niña variety around May-July (Fig. 4b). Therefore, it is unlikely that ENSO is required for the onset of the persistent JS Niña variety. Instead it is likely driven by atmosphere-ocean processes internal to the Indian Ocean.
To further study the potential role of ENSO in driving the JSN varieties, longitude-time composite maps of SSTAs along the equatorial Pacific and tropical South Indian Oceans are constructed for observed El Niño and La Niña events during 1949–2019 (Fig. 4i,j). There is an overall weak connection between the developing phase of ENSO and SSTAs along the Java-Sumatra coast. An IOD-like pattern of SSTAs emerges around SON from the southeast equatorial Indian Ocean to the western equatorial Indian Ocean. However, during the onset phase of JSN (i.e., prior to July), equatorial Pacific SSTAs associated with the developing phase of ENSO are relatively weak; thus, it is unlikely that ENSO drives the onset of JSN. The remote influence of ENSO during the onset, peak and dissipation phases of JSN, and the associated regional rainfall patterns, are further investigated next.
Regional rainfall patterns linked to JSN varieties. Figure 5 shows the rainfall and 850 hPa wind anomalies associated with El Niño and the leading JS Niña varieties for the onset (AMJ), peak (JAS) and dissipation (OND) phases of JS Niña. During the peak phase of JS Niña, the rainfall anomalies associated with El Niño show a predominantly zonal see-saw pattern with decreased rainfall over the Maritime Continent and increased rainfall in the western tropical Pacific (Fig. 5b). The anomalous low-level westerly winds along the western tropical Pacific clearly indicate a weakened Pacific Walker cell. The rainfall and low-level wind anomalies during the peak phase of the double-peak JS Niña variety are largely consistent with those of El Niño (Fig. 5e), suggesting that the remote influence of El Niño via changes in the Pacific Walker cell largely dictates the regional rainfall anomalies associated with the double-peak JS Niña variety.
In contrast, the rainfall anomalies during the peak phase of the persistent JS Niña variety show a meridional see-saw pattern between the Maritime Continent and South-Southeast Asia (Fig. 5h). The associated anomalous low-level winds suggest a Gill-type cross-hemispheric atmospheric response to diabatic cooling anomalies aloft in the southeastern equatorial Indian Ocean34. The Gill-type atmospheric circulation carries extra moisture from the southeastern equatorial Indian Ocean region to South-Southeast Asia and can be interpreted as a strengthening of the South and Southeast Asian summer monsoons. The on-season JS Niña variety during its peak phase also shows a similar meridional see-saw pattern of rainfall anomalies and the Gill-type atmospheric response, though the rainfall and low-level wind anomalies are generally weaker than those of the persistent JS Niña variety (Fig. 5k). Additionally, the teleconnection to southeastern Australia, which is known to occur in boreal summer during some IOD events14, is observed only during the peak phase of the persistent JS Niña variety (Fig. 5h).
During the dissipation phase (OND) of JS Niña, the rainfall anomalies associated with El Niño and the double-peak, persistent and on-season JS Niña varieties all show a zonal tripole pattern with increased rainfall from eastern Africa to the western tropical Indian Ocean, decreased rainfall over and around the Maritime Continent, and increased rainfall over the western tropical Pacific Ocean (Fig. 5c,f,i,l). The zonal tripole pattern of rainfall anomalies that appears during the dissipation phase of the persistent JS Niña variety is particularly strong, with large rainfall anomalies over eastern Africa, northern Australia and the Maritime Continent, whereas the rainfall response to the on-season JS Niña variety is generally very weak in OND.
Figure 6 shows the rainfall and 850 hPa wind anomalies associated with La Niña and the leading JS Niño varieties for the onset, peak and dissipation phases of JS Niño. The rainfall and low-level wind anomalies associated with La Niña are almost mirror images of those associated with El Niño, indicating a zonal see-saw pattern of rainfall anomalies and a strengthened Pacific Walker cell during the peak phase of JS Niño (Fig. 6b). The rainfall and low-level wind anomalies associated with the peak phase of the double-peak JS Niño variety are largely consistent with those of La Niña (Fig. 6e). Interestingly, Gill-type atmospheric circulation anomalies still appear during the peak phase in response to diabatic heating anomalies over the southeastern equatorial Indian Ocean. However, these anomalies are not robust enough to produce a meridional see-saw pattern of rainfall anomalies between the Maritime Continent and South-Southeast Asia. Similarly, the anomalous low-level winds associated with the peak phase of the on-season JS Niño variety show a Gill-type atmospheric response to diabatic heating anomalies over the southeastern equatorial Indian Ocean (Fig. 6e). The associated regional rainfall anomalies are somewhat indicative of a meridional see-saw pattern of regional rainfall anomalies. However, the rainfall response to the on-season JS Niño variety is generally weak.
During the dissipation phase (OND) of JS Niño, the rainfall anomalies associated with La Niña and the double-peak and on-season JS Niño varieties are characterized by a zonal tripole pattern with decreased rainfall from eastern Africa to the western tropical Indian Ocean, increased rainfall over the Maritime Continent and over northern and eastern Australia, and decreased rainfall over the western tropical Pacific Ocean (Fig. 6c,f,i). However, the rainfall response to the on-season JS Niño variety is very weak in OND.
Onset mechanisms of JSN varieties. As discussed earlier, the pre-season and double-peak JS Niño/Niña, and persistent JS Niña varieties are connected to ENSO, while the on-season JS Niño/Niña variety is largely independent of ENSO. Here, we explore in more detail the onset mechanisms of leading JSN varieties. First, we examine the anomalous SSTs, thermocline depth and surface winds associated with El Niño/La Niña and then compare them with those of the leading JS Niña/Niño varieties during the onset, peak and dissipation phases of JSN (Figs. 7 and 8).
SSTAs in the tropical Indian Ocean associated with El Niño are generally weak (Fig. 7a-c). Although weak cold SSTAs appear in a broad region of the southeastern tropical Indian Ocean during the peak phase (JAS) of JS Niña, the surface wind anomalies are too weak to trigger or sustain JS Niña. A positive IOD-like pattern of weak SSTAs develops during the dissipation phase (OND) of JS Niña, with corresponding surface wind anomalies (i.e., southeasterly) that are favorable for coastal upwelling along the Java-Sumatra coast and downwelling in the central tropical South Indian Ocean. However, since the southeasterly monsoonal wind weakens and can reverse its direction as the Australian summer monsoon develops37, the seasonal mean thermocline along Java and Sumatra is too deep to sustain local Bjerknes feedback in OND (Supplementary Fig. 6).
For the double-peak JS Niña variety, cold SSTAs appear in the southeastern tropical Indian Ocean during its onset phase (AMJ), with almost no surface wind anomalies in the region (Fig. 7d). This suggests that the cold SSTAs are the residual of basin-wide cold SSTAs forced remotely during the dissipating phase of La Niña (Fig. 4e). Although still weak, the anomalous SSTs, thermocline depth and surface winds are generally larger than those of El Niño during the peak and dissipation phases of JS Niña (Fig. 7e,f), suggesting an involvement of local Bjerknes feedback in the double-peak JS Niña variety, possibly enhanced by the residual of cold SSTAs formed during the pre-season.
Unlike El Niño or the double-peak JS Niña variety, a robust pattern of upwelling favorable surface wind anomalies appears along the equatorial and southeastern equatorial Indian Ocean during the onset phase of the persistent JS Niña variety (Fig. 7g). Although weak, cold SSTAs appear along the south coast of Java, and the thermocline becomes shallower from the south coast of Java to the northeast equatorial region off Sumatra. During the peak phase, the anomalous surface winds strengthen greatly, and the cold SSTAs intensify and expand to the southwest coast of Sumatra (Fig. 7h). Additionally, the thermocline becomes shallower along the Java-Sumatra coast and in the eastern equatorial region, and the anomaly expands to the west. In contrast, warm SSTAs and a deeper thermocline develop in the central tropical South Indian Ocean. These patterns of amplifying surface wind, SST and thermocline depth anomalies from the onset to peak phase strongly suggest an active role of local Bjerknes feedback.
During the dissipation phase of the persistent JS Niña variety, the cold SSTAs off the Java-Sumatra coast persist but weaken, while the anomalous patterns of surface winds and thermocline depth further intensify and propagate westward (Fig. 7i). These surface wind and thermocline depth anomalies are accompanied by a pronounced basin-wide see-saw pattern of SSTAs between the southeast and southwestern tropical Indian Ocean, thus producing a strong positive IOD. Although the surface wind anomalies are still strong and upwelling favorable, the cold SSTs along the Java-Sumatra coast start to weaken in OND as the southeasterly monsoonal wind weakens and the seasonal mean thermocline along the Java-Sumatra coast becomes too deep to influence the coastal SSTs. These also explain how the weakening and directional shifts in monsoonal winds in boreal winter prevent JSN from continuing beyond OND (Supplementary Fig. 6).
The spatial patterns and amplitudes of anomalous SSTs, thermocline depth and surface winds associated with the on-season JS Niña variety are very similar to those of the persistent JS Niña variety during the onset and peak phases (Fig. 7j,k). However, unlike the persistent JS Niña variety, the anomalous SSTs and surface winds diminish rapidly in OND (Fig. 7l). It is not entirely clear why some JS Niña events are persistent (i.e., the persistent JS Niña variety) and others dissipate more quickly (i.e., the on-season JS Niña variety). One possibility is that some JS Niña events are enhanced and prolonged by remote forcing from the Pacific. Another possibility is that an exceptionally strong JS Niña event may greatly weaken the Indo-Pacific Walker cells, producing warm SSTAs in the neighboring regions of the tropical Indian and Pacific Oceans. The warm SSTAs in the tropical Indian and Pacific Oceans may in turn further weaken the Indo-Pacific Walker cells to prolong the JS Niña event in OND. Consistent with this hypothesis, the persistent JS Niña variety is linked to warm SSTAs in the southwestern tropical Indian Ocean and to weak El Niño conditions (Fig. 4b), but there are almost no such warm SSTAs for the on-season JS Niña variety (Fig. 4f). This hypothesis is also in line with the idea that a strong IOD event may induce a weak (or pseudo) ENSO event38–41.
Anomalous SSTs, surface winds and thermocline depth associated with La Niña are more or less mirror images of those associated with El Niño (Fig. 8a-c). In other words, the anomalies are too weak to trigger or sustain JS Niño during the onset (AMJ) and peak (JAS) phases of JS Niño. During the dissipation phase (OND) of JS Niño, a negative IOD-like pattern of weak SSTAs develops with corresponding surface wind anomalies (i.e., northwesterly) that are favorable for coastal downwelling along the Java-Sumatra coast. However, as discussed earlier, due to the weakening of the southeasterly monsoonal wind and the associated deepening of the coastal thermocline after JAS (supplementary Fig. 6), the anomalous northwesterly winds along the Java-Sumatra coast in OND do not produce a pronounced JS Niño event.
Unlike La Niña, a robust pattern of downwelling favorable surface wind anomalies emerges along the equatorial and southeastern equatorial Indian Ocean during the onset phase of the double-peak JS Niño variety (Fig. 8d). As a result, the thermocline is deeper from the south coast of Java to the northeast equatorial region off Sumatra and shallower immediately southwest of the Java-Sumatra coast. Also during this early period (AMJ), warm SSTAs already prevail widely across the tropical Indian Ocean including in the southeastern equatorial region. The early season warm SSTAs are the residual of pre-season warming, forced remotely during the dissipation phase of El Niño in boreal winter and spring (Fig. 4g). These surface wind and thermocline depth anomalies and the warm SSTAs along the Java-Sumatra coast intensify during the peak phase (Fig. 8e). This suggests that local Bjerknes feedback still plays an active role in the double-peak JS Niño variety, given that the remote influence of La Niña is generally weak during the onset and peak phases of JS Niño (Fig. 8a,b). Consistently, some double-peak JS Niño events occurred under a near-neutral ENSO condition in JAS. For instance, the 2016 JS Niño event, which is the strongest warm event during 1949–2019, occurred under a very weak and short-lived La Niña event and brought extremely high rainfall over Indonesia and Australia in June-September41 (Supplementary Figs. 2 and 5). During the dissipation phase of the double-peak JS Niño variety, warm SSTAs off the Java-Sumatra coast quickly diminish, while the anomalous patterns of surface winds and thermocline depth further intensify and propagate westward (Fig. 8f). Additionally, warm SSTAs develop in the central tropical South Indian Ocean south of 10°S. It is likely that these atmosphere-ocean anomalies in OND are, to a certain extent, enhanced by La Niña through a strengthened Walker cell in the tropical Indian Ocean and the associated anomalous wind stress curl and Ekman pumping.
Similar to the double-peak JS Niño variety, during the onset phase of the on-season JS Niño variety, downwelling favorable surface wind anomalies appear along the equatorial and southeastern equatorial Indian Ocean (Fig. 8g). Although weak, warm SSTAs also appear along the south coast of Java. During the peak phase of the on-season JS Niño variety, the anomalous surface winds intensify; thus, the thermocline is deeper along the Java-Sumatra coast and shallower in the central tropical South Indian Ocean (Fig. 8h). The anomalous SSTs and surface winds quickly diminish during the dissipation phase since the seasonal mean thermocline along the Java-Sumatra coast is too deep to sustain local Bjerknes feedback in OND (Supplementary Fig. 6).
Potential triggering mechanisms of local Bjerknes feedback. It is widely accepted that ENSO is an important forcing mechanism for the IOD17–22, which tends to peak in October1. However, as shown and discussed in this study, during the onset phase of JSN (AMJ), the remote influence of ENSO is generally not robust enough to directly force the development of JSN. Thus, for most of the leading JSN varieties (i.e., the on-season JS Niña/Niño, double-peak JS Niño and persistent JS Niña varieties), local Bjerknes feedback23 is a key ingredient for their emergence. As suggested in Li et al.23, it is likely that JSN is a weakly damped oscillation, which by definition requires a finite amplitude perturbation to trigger and sustain the local Bjerknes feedback.
Previous studies suggested that westward propagating oceanic Rossby waves across the equatorial South Indian Ocean could directly affect atmospheric convection in the western equatorial region to produce the IOD29,42. The oceanic Rossby waves may also reflect from the western boundary to produce eastward propagating equatorial Kelvin waves. These oceanic waves may further disturb the Wyrtki jets, the eastward flowing ocean currents along the equatorial Indian Ocean appearing in boreal spring and fall, and thus potentially trigger the IOD41,43−52. In line with these hypotheses, during the onset phase of the double-peak JS Niño variety, the thermocline across the Seychelles Dome in the southwestern tropical Indian Ocean is deeper (Fig. 8d). The deepened thermocline is remotely forced from the Pacific during the pre-season53–56 (Fig. 4g). It is possible that the anomaly later propagates to the west and reflects from the western boundary to produce eastward propagating equatorial Kelvin waves, initiating local Bjerknes feedback that kick-starts the second peak. It is also possible that the residual of warm SSTAs in the southeastern tropical Indian Ocean formed during the pre-season may persist and enhance the local Bjerknes feedback (Fig. 8d). However, in some cases the pre-season anomalies completely dissipate before boreal summer (i.e., the pre-season JS Niño/Niña variety). In addition, there are little-to-no pre-season anomalies for the on-season JS Niña/Niña variety or the persistent JS Niña variety. Therefore, further study is needed to understand the extent to which the pre-season anomalies of thermocline depth and SSTs can actually trigger the local Bjerknes feedback.
Other proposed onset mechanisms of the IOD involve subseasonal atmospheric processes, such as the Madden-Julian Oscillation and boreal summer intraseasonal oscillation25,57−59, a delayed response to the Pacific ENSO through the Indonesian Throughflow12,60, recharge oscillation61 and inter-hemispheric sea level pressure gradient62. Yet, some other studies point to atmosphere-ocean anomalies originating from the subtropical South Indian Ocean as a potential driver of the IOD24,27. In agreement with this hypothesis, prior to the onset phase of the persistent and on-season JS Niña varieties, an anticyclonic sea level pressure anomaly appears over the subtropical South Indian Ocean in March (Supplementary Fig. 7a,d). This sea level pressure anomaly slowly shifts (or expands) toward the equator, producing anomalous southeasterly winds off the Java-Sumatra coast as early as April (Supplementary Fig. 7b,e). This may in turn serve as a trigger for local Bjerknes feedback. Similarly, prior to the onset phase of the double-peak and on-season JS Niño varieties, a cyclonic sea level pressure anomaly forms over the subtropical South Indian Ocean and then slowly migrates toward the equator (Supplementary Fig. 7g-l). Further study is needed to understand how these subtropical sea level pressure anomalies form and migrate equatorward.