Tropical Pacific Ocean HC in winter 2022
The Climate Prediction Center (CPC) issued a final La Niña advisory in March 2023, signaling the end of the unusually long cold event which lasted from winter 2020 to spring 2023[
8]. Because the large-scale tropical atmospheric and oceanic circulations are closely coupled to each other during ENSO[
9], this rare La Niña was featured by strong and long-lasting easterly anomalies over the equatorial western-central Pacific which, by inducing equatorward Sverdrup transport[
10], have maintained positive thermocline depth anomalies in the western equatorial Pacific (
Figure 1a-c).
Observations and classic ENSO theories suggest that when the thermocline depth in the western basin is deeper than normal in winter, an El Niño is likely to occur at the end of the following year[
4,
5,
11]. We therefore examine the thermocline depth denoted by the 20°C isotherm (D20) averaged from December 2022 to February 2023 (
Figure 1d). For comparison, the averaged thermocline depth from December of the preceding year (referred to as year (-1)) to February of the current year (year (0)) of 1982, 1997, and 2015 are also given (
Figure 1e-g). Deeper-than-normal thermocline was observed in the entire warm pool this winter, with an intensity much stronger than those preceding the onset of the three extreme El Niño events. Note that the spatial pattern of this winter’s thermocline depth anomaly bears a resemblance to that of 1982, a year when HC rather than HFPs was the deciding factor for the subsequent occurrence of super El Niño[
12].
The strong positive thermocline depth anomaly indicates that the western tropical Pacific has accumulated an abnormally large amount of warm water in the upper layer (Extended Data
Figure 1).
Figure 1h displays the integrated warm water volume (WWV) over the equatorial Pacific (120°E-80°W, 5°S-5°N) in D(-1)JF(0), an index often used as a precursor for ENSO event[
3], for every El Niño onset year since 1982, as well as for 2023 (Methods). While not all El Niño events were preceded by WWV/HC buildup, the strong ones certainly did. It is clear that the Pacific Ocean restored the largest HC in winter 2022, and such a large HC buildup appeared in all of the datasets we have analyzed (Extended Data
Figure 2 and
Figure 3).
Impact of HC on the upcoming El Niño
We next conduct a forecast experiment from March 1
st, 2023 using a state-of-the-art ENSO ensemble prediction system which consists of a global coupled model and a data assimilation system (Methods), termed CESM-G here. The NINO3.4 index, which is the average SST anomaly in the central-eastern equatorial Pacific (170°W-120°W, 5°S-5°N), is used to measure ENSO intensity. An event with NINO3.4 greater than 0.5°C is defined as an El Niño, which can be further divided into weak (0.5 to 0.9 °C), moderate (1.0 to 1.4°C), strong (1.5 to 1.9°C) and extreme (>2.0°C) categories[
13]. As shown in
Figure 2a, most ensemble members from CESM-G predict the coming of an El Niño in late 2023. The ensemble-mean of predicted NINO3.4 is ~1.17°C in D(0)JF(1) 2023, meeting the criterion of a moderate event. One notable feature is the level-off of the predictions after September 2023, lacking a distinct peak in winter that is often observed[
11]. Another feature is the rather large spread among 11 ensemble members, with predictions falling into all El Niño categories, and even a case of negative NINO3.4.
The uncertainty in the predicted ENSO intensity resides in the fact that ENSO is a nonlinear system and is affected by many processes at various spatiotemporal scales[
14,
15]. First, a slight initial perturbation may grow fast due to nonlinearity and take the prediction away from its unperturbed track. It is suggested that an initial perturbation close to the optimal growth pattern is more likely to trigger an ENSO event[
16]. Second, the multiplicative HFPs, especially those occurring in sensitive areas and seasons[
17], also have considerable impact on ENSO evolution and intensity[
7]. Finally, some processes out of the tropical Pacific, such as the Pacific meridional mode[
18] and the SST warming in the northern Atlantic[
19], are thought to influence ENSO development.
To focus on the impact of current tropical Pacific Ocean HC on the 2023/24 El Niño, we first try to reduce the uncertainty associated with the processes out of the tropical Pacific by replacing the surface wind stress outside of (100°E-70°W, 30°S-30°N) by the climatology. The predictions using this model setup, termed CESM-TP, are given in
Figure 2b. While considerable forecast spread still exists, the uncertainty is apparently reduced as compared to that of CESM-G. The NINO3.4 index in D(0)JF(1) 2023 is now ~1.69°C, meeting the criterion of a strong El Niño. The increase of El Niño intensity in CESM-TP is clearly a result of reduced active influences from outside of the tropical Pacific, which is probably a mixture of signal and noise[
20]. Considering the fact that the predicted 2023/24 El Niño from CESM-TP also shows a more realistic evolution pattern with a peak around December 2023, we tend to believe that the prediction of CESM-TP is likely to be more robust than that of CESM-G.
To further reduce the uncertainty and isolate the impact of HC buildup, we apply an online low-pass filtering scheme to the CESM-TP (Methods). This scheme is designed to filter out the internal HFPs produced by the model itself which are usually unrealistic[
21], thus letting the predicted ENSO depend on the influence from the low-frequency tropical ocean-atmosphere coupling alone[
22]. The predictions of NINO3.4 using this model setup, termed CESM-TP+LP, are given in
Figure 2c. The forecast spread is sharply reduced in this case, with a more prominent peak phase. The evolution of the predicted El Niño is in strong resemblance to the extreme El Niño events in history[
6]. The ensemble-mean NINO3.4 averaged in D(0)JF(1) 2023 is ~1.78°C, indicating that without the influences from outside of the tropical Pacific and HFPs, the tropical Pacific Ocean HC buildup in winter 2022 unambiguously sets the stage for the coming of a strong El Niño toward the end of 2023.
Forecast uncertainty due to HFPs
The occurrence and intensity of El Niño do not solely depend on the HC buildup. Many studies have shown that ENSO is strongly modulated by atmospheric HFPs, especially the westerly wind burst (WWB) in the western-central equatorial Pacific[
23]. As a matter of fact, every El Niño event during the past 50 years was accompanied by WWB activities[
24,
25]. In particular, the strong WWBs in spring 1997 and 2015 are found to be essential to the genesis of the 1997/98 and 2015/16 extreme El Niño events[
26,
27,
28]. Aside from WWB, the easterly wind surge (EWS) is also found to influence ENSO evolution substantially[
29]. For example, it is argued that the strong EWS in July 2014 played a key role in punctuating the 2014/15 El Niño[
30], though this EWS might be overestimated[
28]. For reliable predictions of El Niño, we have to account for the effects of HFPs[
7,
22,
25,
29], especially WWB and EWS, in addition to the HC buildup.
Figure 3 presents the scatterplots of the HC averaged in D(-1)JF(0) and the integrated HFPs in MAMJ(0) versus the NINO3.4 index averaged in D(0)JF(1) for every El Niño and La Niña event since 1982. Here we take HFPs in June into account because the strong EWS mentioned above occurred in late June 2014. From
Figure 3a, it is clear that although ENSO intensity is positively and significantly correlated to the HC, the degree of HC buildup is not the only determining factor. For example, while the 1982/83 and 1997/98 extreme El Niño events were preceded by a very large HC buildup in the previous winter, the even stronger 2015/16 extreme El Niño was preceded by a moderate HC buildup. On the other hand,
Figure 3b shows that ENSO intensity is also positively correlated with the spring and early summer HFPs. For instance, the WWBs in spring 2015 were strong enough to compensate for the weak HC buildup in the preceding winter, and to trigger and maintain the extreme 2015/16 El Niño event[
28]. The results shown in
Figure 3 can also be found using other datasets (Extended Data
Figure 4).
Given the importance of HFPs to ENSO, it is desirable to include them in forecast models. Unfortunately, present climate models can hardly simulate HFPs in a realistic manner due to their multiplicative nature [
7], as they are closely related to various synoptic and sub-seasonal processes. For example, WWBs are found to be associated with the western Pacific tropical cyclones[
31,
32], East Asian cold surge[
33,
34], and the Madden-Julian Oscillation[
35,
36]. Since WWBs are modulated by ENSO with stronger and more frequent bursts to occur during large El Niño events, it is possible to formulate semi-stochastic, environment-dependent parameterization for them[
37]. The response of ENSO to WWBs is sensitive to the timing of their occurrence[
17], and we have shown previously that the consecutive WWBs in spring and early summer are most effective drivers[
17,
22,
25].
Here we take a simple approach to investigate how HFPs influence the 2023/24 El Niño by randomly adding a set of HFPs to our forecast model in MAMJ(0) (Methods). The number of HFPs are within the range of observations related to ENSO events (Extended Data Table 1), and the location, duration and intensity of them are statistically consistent with observed values. In order to reduce the uncertainty associated with the internal HFPs in the model itself, the forced experiments are conducted using the aforementioned low-pass filtering scheme, and are termed as CESM-TP+LP+HFPs. Thus the impact of the stochastic forcing of HFPs that may occur in spring and early summer can be evaluated on top of the deterministic low-frequency atmosphere-ocean coupling started from the current tropical Pacific Ocean HC buildup. The experiments are carried out from three initial conditions which predict the strongest, weakest, and medium event in D(0)JF(1) 2023 in the CESM-TP+LP (
Figure 2d).
The predictions from the medium case are shown in
Figure 4a for different number and type of HFPs. In general, El Niño intensity increases with more WWBs and decreases with more EWSs. In the case where four WWBs occur in MAMJ(0), the ensemble-mean of NINO3.4 in D(0)JF(1) is ~2.36°C. For reference, the WWB forcing in MAMJ(0) 1997 was equivalent to 5.04 WWBs used this experiment (Extended Data Table 1), and the intensity of 1997/98 El Niño was ~2.24°C. Another point worth noting is that adding two EWSs in MAMJ(0) prolongs the negative NINO3.4 to July 2023, yet the model still predicts a moderate El Niño in D(0)JF(1). This is different from Hu and Fedorov[
30] who showed that the cease of the 2014/15 El Niño event was due to one shot of EWS. This is probably because the HC buildup in winter 2022 was much larger than that in winter 2013 so that even two strong EWSs cannot stop the warming tendency.
The changes of El Niño intensity from unforced to forced runs for all three groups of experiments are given in
Figure 4b. Just like what we have seen in
Figure 4a, adding WWBs (EWSs) increases (decreases) the 2023/24 El Niño intensity in all groups of model runs. However, the impact of WWBs decreases as the original intensity of El Niño increases. In other words, WWBs matter less when the HC buildup is large enough for a strong El Niño to come (as in 1982), but they are essential for the development of a strong El Niño when the HC buildup is relatively small (as in 2015)[
12]. Considering the unprecedentedly large HC buildup this winter, we can reasonably expect the coming of a strong El Niño in spite of what would happen with the HFPs, as long as their number and magnitude are within the range of historically observed values.
In summary, the analyses presented here reveal that the buildup of the equatorial Pacific Ocean HC in winter 2022 ranks the largest over the past 40 years (
Figure 1h). This HC buildup alone is sufficient to boost a strong El Niño event in late 2023, and such a daring prediction is not sensitive to the uncertainty in initial conditions (
Figure 2c). Other uncertainties mostly come from two sources. One is the effects of randomly appearing HFPs, which would enhance the intensity of El Niño and even push it into an extreme event in the case of WWB, and would lower the magnitude of El Niño in the case of EWS (
Figure 4). The other source of uncertainty is the influences from the processes outside of the tropical Pacific, which might degrade the 2023/24 El Niño to a moderate event if the prediction of our global model is to be trusted (
Figure 2a). At any rate, the stage is set and the curtain is drawn for an El Niño to make its grand appearance.