ENSO phase space dynamics in CMIP models

ENSO phase space dynamics in CMIP models | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article ENSO phase space dynamics in CMIP models Priyamvada Priya, Dietmar Dommenget This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4727039/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract This study analyses the El Niño-Southern Oscillation (ENSO) phase space as simulated by the Coupled Model Intercomparison Projects 5 and 6 (CMIP5 and CMIP6) models. The ENSO phase space describes the ENSO cycle between the sea surface temperature (SST) anomaly in the eastern equatorial Pacific ( T ) and the equatorial mean thermocline depth anomaly ( h ). We find that the characteristics out-of-phase cross-correlation between T and h is shifted to negative values in CMIP models, suggesting that the coupling between T and h is regionally sifted to the east compared to the observed central Pacific. If we consider the CMIP models with an eastward shifted h then the models have better agreements with the observed characteristics. While the models can capture some of the non-linear aspects with high correlations, they do largely underestimate the strength of non-linear ENSO aspects. They underestimate the likelihood of extreme El Niño and discharge states, they cannot capture the enhanced growth rates during the recharge state, the enhanced decay after the discharge state nor the reduced phase transitions after the La Niña phases. Weaker than observed wind-SST feedback and weaker h variability are likely some of the reasons why models cannot fully capture the non-linear ENSO phase space dynamics. Further, we found no indication of significant improvements from the CMIP 5 to 6 ensemble, suggesting that the two ensembles are essentially the same in terms of their ENSO dynamics. There is, however, a large spread within the model ensembles, leading to models with quite different ENSO dynamics. El Niño-Southern Oscillation (ENSO) dynamics recharge oscillator model ENSO phase space ENSO events CMIP5 and CMIP6 models non-linear dynamics Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Figure 10 1. Introduction The El Niño-Southern Oscillation (ENSO) is a climate phenomenon that influences the global climate and weather patterns, ecosystems, and society (Tsonis et al. 2003 ; McPhaden et al. 2006 ; Latif and Keenlyside 2009 ; Timmermann et al. 2018 ; Alizadeh 2022 ). It is characterized by its two distinct phases which are the warm phase (El Niño) and cold phase (La Niña) in the equatorial east Pacific Ocean (Chen and Li 2021 ). The variations of ENSO are controlled by complex dynamics that include feedback from the atmosphere and ocean (Bellenger et al. 2014 ; Kim et al. 2014 ; Wang et al. 2017 ; Vijayeta and Dommenget 2018 ; Planton et al. 2021 ). The aim of this study is to contribute to a better understanding of these complex dynamics. Key ENSO dynamics can be determined by the interaction of ocean-atmosphere processes such as sea surface temperature (SST), wind stress, and thermocline depth, known as Bjerknes feedback (Bjerknes 1969 ). A simplified concept for combining these processes and for understanding ENSO dynamics is described by the recharge oscillator (ReOsc) model (Jin 1997 ; Burgers et al. 2005 ). The ReOsc model describes the variability of ENSO as an interaction between the tendencies of SST anomaly in the eastern equatorial Pacific ( T ) and equatorial mean thermocline depth anomaly ( h ), with the latter representing a proxy for the subsurface equatorial heat content. It has been used extensively in previous studies to analyse ENSO dynamics (Jin 1997 ; Burgers et al. 2005 ; Frauen and Dommenget 2010 ; Yu et al. 2016 ; Vijayeta and Dommenget 2018 ; Wengel et al. 2018 ; Dommenget and Vijayeta 2019 ; Crespo et al. 2022 ; Dommenget and Al-Ansari 2022 ; Dommenget et al. 2023 ; Priya et al. 2024 ). The idealised ReOsc model is an oscillation between T and h , presenting a 2-dimensional phase space, as known from physics for a harmonic oscillator or complex system analysis (Meinen and McPhaden 2000 ; Kessler 2002 ; Dommenget and Al-Ansari 2022 ; Dommenget et al. 2023 ; Priya et al. 2024 ). Detailed analysis of the characteristics of the ENSO phase can help to understand how the dynamics of ENSO control the amplitude, non-linear behaviour, and time evolution as functions of the ENSO phases. The study by Dommenget et al. ( 2023 ) showed that the ENSO phase space highlights strong non-linear aspects, with enhanced extremes in the El Niño to discharge states, larger growth rates in the transition from recharge to El Niño state, strong decay after the discharge states, faster phase speeds after El Niño states and slow phase speeds after La Niña states. These non-linear aspects of the ENSO cycle have also been reported in other studies, using different measures (An and Jin 2004 ; An et al. 2005 ; Dommenget et al 2013 ; Sun et al. 2016 ; Okumura 2019 ; An et al. 2020 ). The Coupled Model Intercomparison Projects (CMIP) models are widely used to simulate and forecast ENSO activity. They are highly capable of simulating many aspects of ENSO and do show some improvements in the consecutive generations of CMIP from phases 3 to 6 in metrics used to evaluate ENSO, such as seasonal cycle, ENSO patterns, and diversity of ENSO events (Bellenger et al. 2014 ; Kim et al. 2014 ; Planton et al. 2021 ). Despite their capabilities, the models do show some biases if compared against observations and they do show large variations within the model ensemble (Kim and Yu 2012 ; Guilyardi et al. 2012 ; Bellenger et al. 2014 ; Vijayeta and Dommenget 2018 ; Bayr et al. 2019 ; McKenna et al. 2020 ). The CMIP models exhibit biases in main atmospheric and oceanic feedback, which do control ENSO amplitude and stability (Guilyardi et al. 2012 ; Bellenger et al. 2014 ; Kim et al. 2014 ; Vijayeta and Dommenget 2018 ; Bayr et al. 2019 ; Beobide-Arsuaga et al. 2021 ). Bayr and Latif ( 2023 ) discussed how atmospheric feedback, specifically zonal wind feedback and heat flux feedback, are underestimated in CMIP model (CMIP5 and CMIP6) simulations of ENSO dynamics. Further, the models exhibit biases in their representation of the SST patterns, with a westward extension compared to observed data in CMIP5 and CMIP6 models (Weller and Cai 2013 ; Yang and Giese 2013 ; Bellenger et al. 2014 ; Capotondi et al. 2015 ; Wang et al. 2015 ; Guilyardi et al. 2020 ; Jiang et al. 2021 ; Planton et al. 2021 ). The models also exhibit biases in the equatorial thermocline depth, intensity, and zonal shape due to biases in the equatorial zonal wind stress, which affects upwelling and the zonal tilt of the equatorial thermocline (Guilyardi et al. 2020 ). About 80% of the variations in simulating ENSO amplitude within the CMIP model ensemble are due to variations in the linear growth rates of T and h , and in the strength of the stochastic noise forcing, with the latter also potentially being linked to non-linear processes (Wengel et al. 2018 ) (non-linear processes such as non-linear air-sea coupling, non-linear dynamical heating, non-linear ocean advection, non-linear relationship of precipitate-SST feedback, heat flux because of the tropical ocean instability waves, atmospheric non-linear response to the SST anomaly (wind-SST relation) (An and Jin 2004 ; An et al. 2005 ; Sun et al. 2016 ; Okumura 2019 ; An et al. 2020 )). Several studies on ENSO non-linearities point out that CMIP models tend to underestimate the non-linear aspects of ENSO (An and Jin 2004 ; An et al. 2005 ; Dommenget et al. 2013 ; Sun et al. 2016 ; Okumura 2019 ; Planton et al. 2021 ). In this study, we analyze the ENSO phase space statistics as simulated in the CMIP 5 and 6 simulations and compare them against the observed values. We examine how well the CMIP models can reproduce the observed ENSO phase space asymmetries to gain some understanding of how well models can simulate the non-linear behavior of ENSO as a function of the different phases of ENSO. We explore potential limitations in model simulations that may affect the model performance. Additionally, we are also interested in assessing the improvements from CMIP5 to CMIP6 models in representing the ENSO phase space characteristics. This paper is structured as follows: Section 2 describes the observational datasets used, CMIP simulations, phase space statistics methods, and the ReOsc model. This is followed by the main results section, in which we present the ENSO phase space statistics for CMIP simulations and observations. Here we also analyze the linear ReOsc model fit to CMIP models and discuss some aspects of the wind, thermocline, and SST cross-relations in the context of ENSO dynamics. The study is concluded with a summary and discussion section. 2. Data and Methods 2.1 Observational datasets Observational datasets for SST and the thermocline depth estimates have been taken from the Hadley Centre global sea Ice coverage and SST (HadISST) (Rayner et al. 2003 ), UK Met Office Hadley Centre 'EN' series observations datasets, version 4 (EN4) (Good et al. 2013 ) respectively from 1980 to 2021. The reanalysis datasets for SST and thermocline depth estimates are taken from Simple Ocean Data Assimilation version 3 (SODA3) (Carton and Giese 2008 ) (1980–2017), Ocean Reanalysis System 5 (ORAS5) (1980–2018) (Zuo et al. 2019 ), and CMCC Historical Ocean Reanalysis system: CHOR AS and CHOR RL (1980–2010) (Yang et al. 2017). To observe the variability more precisely, these 5 different datasets are concatenated to create a considerable time series of SST and thermocline depth repeated for the same year. Anomalies are calculated by removing the mean seasonal cycle (calculated over the entire period) from the data. T is the SST anomaly calculated for the NINO3 region (5°S to 5°N, 150°W to 90°W) and considered for the same in this study. Previous studies often use the depth of the 20 o C ( Z 20 ) as a proxy for thermocline depth ( h ), however, according to Dommenget et al. ( 2023 ), the maximum gradient in the temperature profile ( Z mxg ) is a better proxy of h , which is calculated by averaging Z mxg over the equatorial Pacific region (5°S to 5°N, 130°E to 80°W), for the analysis of the ENSO dynamics within the framework of the ReOsc model than the widely used Z 20 . The ENSO phase space for T index by using Z mxg as a proxy of h is closer to the idealized recharge oscillator model than Z 20 . This motivates us to use Z mxg to estimate the equatorial Pacific Ocean thermocline depth h , following the method of Dommenget et al. ( 2023 ). The zonal wind stress data are calculated from 10-m surface zonal wind datasets taken from the European Centre for Medium-Range Weather Forecasts (ECMWF-ERA5) reanalysis (1980–2021) (Hersbach et al. 2020 ) to analyse the wind dynamics of ENSO. Here, we only consider wind stress data from atmospheric reanalysis datasets, ERA5, following a similar approach to previous papers (Guilyardi et al. 2012 ; Bellenger et al. 2014 ; Vijayeta and Dommenget 2018 ; Guilyardi et al. 2020 ). Table 1 CMIP5 and CMIP6 models with the corresponding model numbers CMIP5 CMIP6 Model No. Model Name Model No. Model Name Model No. Model Name Model No. Model Name 1. ACCESS1-0 21. GISS-E2-H-CC 1. ACCESS-CM2 21. EC-Earth3-Veg 2. ACCESS1-3 22. GISS-E2-R 2. ACCESS-ESM1-5 22. EC-Earth3-Veg-LR 3. bcc-csm1-1 23. GISS-E2-R-CC 3. BCC-ESM1 23. FGOALS-f3-L 4. bcc-csm1-1-m 24. HadCM3 4. BCC-CSM2-MR 24. FGOALS-g3 5. CanESM2 25. HadGEM2-CC 5. CanESM5 25. FIO-ESM-2-0 6. CCSM4 26. HadGEM2-ES 6. CESM2 26. GFDL-CM4 7. CESM1-BGC 27. IPSL-CM5A-LR 7. CESM2-FV2 27. GISS-E2-1-G 8. CESM1-CAM5 28. IPSL-CM5A-MR 8. CESM2-WACCM 28. GISS-E2-1-G-CC 9. CESM1-FASTCHEM 29. IPSL-CM5B-LR 9. CESM2-WACCM-FV2 29. INM-CM4-8 10. CESM1-WACCM 30. MIROC5 10. CIESM 30. INM-CM5-0 11. CMCC-CESM 31. MPI-ESM-LR 11. CMCC-CM2-SR5 31. MIROC-ES2L 12. CMCC-CM 32. MPI-ESM-MR 12. CMCC-ESM2 32. MPI-ESM1-2-LR 13. CMCC-CMS 33. MPI-ESM-P 13. CNRM-CM6-1 33. MPI-ESM-1-2-HAM 14. CNRM-CM5 34. MRI-CGCM3 14. CNRM-ESM2-1 34. MRI-ESM2-0 15. CNRM-CM5-2 35. MRI-ESM1 15. E3SM-1-0 35. NESM3 16. CSIRO-Mk3L-1-2 36. NorESM1-M 16. E3SM-1-1 36. NorESM2-MM 17. GFDL-CM3 37. NorESM1-ME 17. E3SM-1-1-ECA 37. SAM0-UNICON 18. GFDL-ESM2G 18. EC-Earth3 38. TaiESM1 19. GFDL-ESM2M 19. EC-Earth3-AerChem 39. UKESM1-0-LL 20. GISS-E2-H 20. EC-Earth3-CC 2.1 CMIP Simulations For the CMIP simulations, we used the historical simulations of the CMIP5 from the years 1900 to 1999 (Taylor et al. 2012 ) and CMIP6 for the period 1850–2014 (Eyring et al. 2016 ) models, utilizing every model with the variables required for this analysis. We use a single realization of each model and use all models with the required data available. See Table 1 for the 37 CMIP5 and 39 CMIP6 models taken into consideration. 2.2 ENSO Phase space For the ENSO phase space analysis, we are following the approach by Dommenget and Al-Ansari ( 2022 ), Dommenget et al. ( 2023 ), and Priya et al. ( 2024 ). We normalise the monthly means of T and h by dividing them by their standard deviations, resulting in a non-dimensional presentation of the variables, T n and h n , respectively. The ENSO phase space is represented with T n on the x-axis and h n on the y-axis. The cartesian coordinate system has been transformed into a spherical coordinate system with phase angle, 𝜑 = 0° in h (y-direction) and 90° in T (x-direction), where the phase angle, 𝜑, rotates in a clockwise direction. The ENSO system anomaly, S , is represented by a vector with two cartesian coordinates T n and h n . When it is converted into the polar coordinate system, the magnitude of S remains constant for constant radius and is independent of phase angle (𝜑). The ENSO system is further described by the magnitude S and phase 𝜑. The ENSO system's tendencies in the polar coordinate system are separated into radial and tangential components. The phase dependent growth rate is calculated by dividing the radial component of the tendencies by the magnitude of S , which is defined as the tendency to move away from the origin (positive values) and towards the origin (negative values). The phase speed is calculated by dividing the tangential component of the tendencies by the magnitude of S and is defined as the system's tendency to circle around the origin. Positive values indicate clockwise motion, while negative values indicate anti-clockwise motion; additionally, smaller positive values indicate a slower transition in the ENSO cycle. The average time required to complete one complete cycle is determined through the integration of the phase speed from 0° to 360°. Note, that a proper analysis of the ENSO phase space statistics assumes that T and h are orthogonal to each other (Dommenget et al. 2023 ). Thus, they have a zero cross-correlation at time lag zero. If T and h are not orthogonal to each other (e.g., lag-zero cross-correlation is not zero), then the phase space statistics are skewed towards the diagonals between the T -axis and h -axis of the diagrams, which will result into skewed statistics of probabilities, growth rates, and phase speeds. 2.3 ReOsc Model The linear ReOsc model is based on two tendency equations with six parameters as given below (Burgers et al. 2005 ): $$\:\frac{dT\left(t\right)}{dt}={a}_{11}T\left(t\right)+\:{a}_{12}h\left(t\right)+{\xi\:}_{1}$$ 1 $$\:\frac{dh\left(t\right)}{dt}=\:{a}_{21}T\left(t\right)+\:{a}_{22}h\left(t\right)+\:{\xi\:}_{2}$$ 2 where a 11 and a 22 are the damping factor or growth rate of T and h respectively, a 12 is the coupling factor for T to h , a 21 is the coupling factor for h to T , ξ 1 is the noise forcing of T , and ξ 2 is the noise forcing of h . The model parameters of the ReOsc model are determined by carrying out a multivariate linear regression on the monthly mean tendencies of T and h in relation to their monthly mean values of T and h (Burgers et al. 2005 ; Jansen et al. 2009; Dommenget and Al-Ansari 2022 ; Dommenget et al. 2023 ; Priya et al. 2024 ). In reference to the normalised presentation of the ENSO phase space ( T n and h n ), we normalise all values of the ReOsc model. The ENSO phase space resulting from the linear ReOsc model fit can be estimated by integrating the Eqs. ( 1 – 2 ) with random white noise forcings ξ 1 and ξ 2 over a period of 10 4 . 3. CMIP simulations of ENSO dynamics Our main results section is divided into four subsections. We start the analysis in subsection 3.1 with investigating the cross-correlation between T and h , which is a fundamental statistic on which the ENSO phase space analysis is based on. This analysis will also focus on the regional shift in ENSO patterns, as it affects cross-correlation between T and h . In the following subsection 3.2, we analyse the ENSO phase space in CMIP simulations, which is the main focus of this study. This is followed by the analysis of the linear ReOsc model fitted to the CMIP simulations. We conclude our discussion by examining the effect of wind stress on ENSO dynamics and how well CMIP models capture this atmospheric feedback in the subsection 3.4. 3.1 Regional variations of T and h Fig. 1a shows the observed cross-correlation between T and h . It shows the characteristic out-of-phase relation between T and h with positive cross-correlation when h leads T , zero cross-correlation at lag zero, and strong negative correlations when T leads h . The CMIP 5 and 6 ensembles are very similar to each other, and both show significant deviations from the observed cross-correlation (Fig. 1 a). All CMIP 5 and 6 simulations have negative cross-correlation at lag zero, shifting the whole ensemble of simulations into clear negative correlations at lag zero. This shift also affects the out-of-phase cross-correlations with lower positive correlation when h leads T that tend to peak at longer (10mon.) lead times than observed (7mon.). In turn, the negative cross-correlations when T leads h are stronger in the CMIP models than in the observations and tend to have short lead times (5mon.) than observed (7mon.). The shift towards negative lag-zero cross-correlations between T and h in CMIP simulations could indicate an eastward shift in the thermocline depth variability associated with ENSO in the CMIP simulations. This would be consistent with observed regional variations in the SST cross-correlation with thermocline depth variability along the equatorial Pacific, which is known to be more positive in the eastern Pacific and more negative in the western Pacific (Burgers and Stephenson 1999 ; Dommenget et al. 2023 ). Given that the central assumption in the ENSO phase space analysis is that the cross-correlation between T and h is zero at lag zero, we now focus on the hypothesis that regional shifts in the simulated ENSO patterns could affect ENSO dynamics (e.g., the cross-correlation between T and h ) in CMIP simulations. We now analyse the equatorial patterns of the ENSO modes in the models, focussing on the SST pattern first and then on the more important regional shifts in the thermocline depth variability. The T index in the ReOsc model is most closely related to the empirical orthogonal function (EOF) mode-1 of the tropical Pacific SST, see Fig. 2 a. It is marked by a horseshoe pattern with warm SST anomalies in the eastern equatorial Pacific that transition into negative SST anomalies in the western equatorial and off-equatorial Pacific. This transition crosses zero equatorial SST anomalies at the longitude of about 159 o E (see dashed line in Fig. 2 a). In Figs. 2 b-d, we show three scenarios for EOF mode-1 of SST from CMIP simulations to highlight cases closest to the observed (Fig. 2 b; ACCESS-CM2), average (Fig. 2 c; MRI-ESM2-0) and worst (Fig. 2 d; NorESM2-MM) in relation to the observed transition line for zero equatorial SST anomalies. While the best CMIP simulation is close to the observed with some negative SST anomalies in the western equatorial Pacific, the worst model has essentially only positive SST anomalies over the whole equatorial Pacific. In general, the EOF mode-1 of tropical Pacific SST of CMIP simulations is similar to the observed, but the transition line for zero equatorial SST anomalies is, in ensemble mean, about 20 o further west for both CMIP 5 and 6 simulations (Fig. 3 a). This is consistent with the westward shift reported in previous studies (Weller and Cai 2013 ; Yang and Giese 2013 ; Bellenger et al. 2014 ; Capotondi et al. 2015 ; Wang et al. 2015 ; Guilyardi et al. 2020 ; Jiang et al. 2021 ; Planton et al. 2021 ). The observed EOF mode-1 of thermocline depth variability is a dipole pattern (tilting mode; not shown) that has an in-phase relation with the T index in the ReOsc model and thus does not represent the h index (Meinen and McPhaden 2000 ). The EOF mode-2 of the equatorial Pacific is similar to a basin-wide mode, but the pattern is somewhat shifted to the eastern Pacific, see Fig. 2 e. This mode is more closely correlated to h (correlation = 0.91). While the EOF mode-2 is mostly a monopole, it does have a transition from negative anomalies to positive anomalies in the far western equatorial Pacific at a longitude of about 146°E (see dashed line in Fig. 2 e). The EOF mode-2 of the CMCC-CM2-SR5 (best model) exhibits a similar spatial pattern, with near-zero values present around the same region with observation (Fig. 2 f). A medium CMIP6 model, ACCESS-ESM1-5, demonstrates that the near-zero values shift towards the eastern equatorial Pacific region, and we can see a dipole-like structure with deeper thermocline depth in the western and shallower thermocline depth in the central to eastern pacific region (Fig. 2 g). In addition, the CNRM-CM6-1 CMIP6 model (worst model) exhibits a clear dipole pattern in EOF mode-2, more similar to the tilting mode described in Meinen and McPhaden ( 2000 ; Fig. 2 h). In the ensemble of CMIP5 and 6 models, we can see that the EOF mode-2 patterns of thermocline depth ( h ) suggest an eastward shift in the equatorial Pacific region of about 20°E (Fig. 3 b). In summary, we find significant regional shifts in the SST and thermocline depth patterns associated with the ReOsc model. Interestingly, we find that the SST pattern is shifted to the west and the thermocline depth pattern is shifted to the east. At this point, we don't have a clear understanding of the reasons behind these biases, but it appears very likely that such regional shifts can affect the cross-correlation between T and h . Following the analysis of Dommenget et al. ( 2023 ) we know that an eastward shift of the h index region leads to a more positive lag-zero cross-correlation between T and h . This could compensate for the negative lag-zero cross-correlation between T and h found in CMIP models (Fig. 1 a). To fit into the ReOsc model and to follow the phase space concept, CMIP models should have an out-of-phase relationship between T and h (e.g., zero lag-zero cross-correlation). To define a revised h index for the CMIP models which has zero lag-zero cross-correlation between T and h , we shifted the index region for h eastward along the equator ( h shift ) until we found a region that most closely follows the observed zero lag-zero cross-correlation between T and h , for all CMIP simulations. We concluded that the out-of-phase characteristics for CMIP models in the ensemble mean are very similar to observation when h shift is considered for the east equatorial pacific region (5°S to 5°N, 190°E to 80°W) (Fig. 1 b). However, individual models do have some significant variations from this overall good fit. In addition to this, we also investigated various regions of T to understand the impact of the SST pattern shift on the ENSO phase space for CMIP models. However, we found that changing the T regions had no significant impact on the out-of-phase relationship between T and h. As a result, we have focused solely on the regional shift of h for CMIP models for the remaining analysis. 3.2 ENSO phase space analysis for CMIP models The observed ENSO phase space is depicted in Fig. 4 a. The main features show a clockwise rotation of the tendencies, with clear phase propagation through all four phases of the ENSO cycle. Unlike the phase space of an idealised linear ReOsc model, which is symmetric in all phases (Dommenget and Al-Ansari 2022 ), the observed variability is clearly skewed towards El Niño (positive T n values) to discharge states (negative h n values) (quarter Q2) and the tendencies are stronger in quarter Q2 and weaker in quarter Q4 (La Niña to recharge phases). A more detailed discussion is presented in Dommenget et al. ( 2023 ). Examples of CMIP models are shown for the phase space estimated based on h (left column in Fig. 4 ) and based on h shift (right column), see Fig. 4 . We selected a best (NorESM2-MM), medium (CESM2-WACCM), and worst (BCC-ESM1) model based on the correlation values between the mean phase space estimated based on h shift and the observed value. If we compare the left with the right column, we can notice that the phase spaces of the left column are pronounced along the diagonal of Q2 to Q4, which is a signature of a negative cross-correlation between T and h (Dommenget et al. 2023 ). The right column, in turn, is more equally distributed on all four quarters of the phase space and generally, closer to the observed behaviour. This highlights that the CMIP models compare better with the observed phase space when h shift is considered. The best model (NorESM2-MM) closely reproduces the observed ENSO phase space with high correlation in the mean phase space (0.9), clear clockwise propagation in all four phases, and variability skewed towards stronger positive T n (El Niño), and stronger negative h n (discharge phase). In turn, the worst model (BCC-ESM1) shows opposite asymmetries in the phase space to those observed, with a negative correlation in the mean phase space (-0.85), indicating that it has extreme recharge states instead of the observed extreme discharge states. However, it also shows clear clockwise propagation in all four phases. The mean statistical features in the ENSO phase space are shown for all CMIP models and the observations in Fig. 5 . Again, the left column shows the CMIP models based on h and the right column based on h shift . The observed values are shown in both columns for comparison (both based on h ). The CMIP models phase space based on h (left column) shows clear biases in all statistics. The mean phase space and the probability distribution of S > 2 are both more pronounced along the diagonal of Q2 to Q4, the growth rate is more pronounced along the T -axis and the phase speed along the diagonal of Q1 to Q3. All these biases directly result from the negative cross-correlation between T and h (Dommenget et al. 2023 ). As a result, the CMIP models phase space statistics based on h do not correlate well with the observed values. The phase space statistics of the CMIP models based on h shift show much less obvious biases and correlate much better with the observed values. The CMIP models generally do capture the asymmetry in the mean phase space, with the largest mean value in Q2 (El Niño and discharge state), smallest value in Q4 (La Niña and recharge state), and a correlation of the ensemble mean value in CMIP 5 and 6 of 0.88 and 0.81, respectively. However, the CMIP ensemble does not quite capture the intensity of the shift from Q4 to Q2, indicating the CMIP models are underestimating the non-linearity in the ENSO phase space. This also holds for the probability distribution of S > 2 (extreme values), which in observations show a clear shift to Q2 and away from Q4. This is captured by the CMIP models, but again not with the right intensity. In particular, the absence of observed extremes in Q4 is not fully captured by the CMIP models. Furthermore, we see a very large spread in CMIP ensemble members, indicating widely different behaviour between the models. The observed mean growth rate as a function of the ENSO phase shows positive values mostly around the recharge towards the El Niño state and negative values in the transition from discharge to a La Niña state (Fig. 5 g). This signature is somewhat captured in the CMIP ensemble, but only with a moderate correlation. The CMIP models capture the negative growth rate in the transition from discharge to a La Niña state better than the positive values during the recharge to the El Niño state. Again, the CMIP models show a wide spread between the ensemble members. Notably, we find several models with very strong negative growth rates around the discharge and recharge phases, suggesting that these ENSO states can collapse much faster in these models than observed. The observed phase speed is fastest, or most clear, in Q2 (after El Niño states) and slowest, or least clear in Q4 (after La Niña states; Fig. 5 h). The CMIP model ensemble can partly capture this signature, but only with a moderate positive correlation. In particular, the models have problems in capturing the smaller phase speed after La Niña events. This suggests that models tend to have a faster or clearer transition from La Niña states to the recharge state. The model ensemble also has a wide range of different phase speeds. Some models transition through some ENSO phases more than twice as fast than observed, and some models have phase speeds near zero in some ENSO phases, suggesting the models have no clear ENSO cycle in these phases. The phase speed analysis suggests that CMIP models are oscillating a bit faster than observed. Figure 6 shows the observed and CMIP simulated power spectrum of the T index. We can first notice that the ensemble mean power spectrum of CMIP 5 and 6 are in very good agreement with the observed, with no clear shift in the power to higher frequencies. However, both CMIP ensembles have a small tendency to peak at slightly shorter periods (~ 3yrs period) than observed (~ 4yrs period). Individual members of the ensemble can, however, behave quite differently, with some models having clear peaks at different periods and other models having less obvious peaks with a more continuous power spectrum. 3.3 Fitted linear ReOsc model for CMIP The ENSO phase space is based on the ReOsc model. It therefore does help to analyse the ReOsc model fitted to the data of the CMIP models to understand what is causing the characteristics of the ENSO phase space. Here we focus on the linear model, thus not considering non-linearities, and without considering any seasonality in the fitted parameters. Neglecting non-linearities is a strong limitation, as will be shown below, and is only reflecting that we yet do not know what non-linear model can describe the observed asymmetries in the ENSO phase space (Dommenget and Al-Ansari 2022 ; Dommenget et al. 2023 ). Figure 7 shows the fitted ReOsc model parameters for all CMIP models and observations for models based on h (left column) and based on h shift (right column). The observed estimate is the same in both columns. Here all parameters are normalized (as in T n and h n ) to allow a better discussion of the dynamical implication for the ENSO phase space. We can first-of-all note that there is quite some spread between models, with many models having parameters being far away from the observed values in both estimates. Further, we find that the two estimates of the ReOsc models are quite different from each other. The models based on h are strongly biased in the growth rates and the noise strength in h relative to the observed values but are much less biased in the coupling parameters. In turn, the models based on h shift are more strongly biased in the coupling parameters. The results for CMIP5 models are also somewhat different from Vijayeta and Dommenget ( 2018 ; hereafter VD18), which used Z 20 to estimate h instead of Z mxg . This does alter the ReOsc model parameter fits. However, the CMIP5 models tend to overestimate the damping (negative a 22 ) of h in both estimates present here and in VD18. We now focus on the ReOsc model fits based on h shift as they better represent the observed ENSO phase space. Here we are primarily interested in what dynamical elements of the ReOsc model are causing the simulated ENSO phase space characteristics. Any deviation of the ENSO phase space from a perfect cycle (a perfectly symmetric phase space), can be linked to asymmetries (deviations for the dashed lines in Fig. 7 ) in the ReOsc models parameter pairs: growth rates (a 11 , a 22 ), coupling (a 12 , a 21 ) or noise forcings (ξ 1 , ξ 2 ). In a non-linear model, they can also result from non-linearities in any of these aspects. Individual models do show strong deviations from the dashed lines in the growth rates and coupling parameters, but less so in the noise forcings. In the ensemble mean the largest asymmetries are found in the coupling parameters, where nearly all models are below the dashed line, indicating a stronger coupling of T to h (a 12 ) than the coupling of h to T (a 21 ), which is quite different from the observed symmetric coupling. The coupling parameters are most important for the period of ENSO (Wyrtki 1985 ). The Wyrtki-index, which measures the peak period of ENSO as a function of a 12 and a 21 , is for most models slightly shifted to shorter periods (Fig. 7 e), consistent with somewhat larger phase speed values seen in the phase space diagrams (Fig. 5 h). Integrating the linear ReOsc model fitted to each CMIP model allows us to compute the phase space statistics based on the resulting T n and h n values, see left column in Fig. 8 . The phase spaces resulting from the ReOsc model are in general similar to those of the CMIP models, but, due to the linear approach, the phase spaces are by construction symmetric (e.g., all structures are mirroring at the origin) and therefore lack all asymmetric variations. This is particularly important when compared against the observed phase space, because most of the interesting observed phase space characteristics are asymmetric. None of the observed asymmetries can be captured by the analysis of the linear ReOsc model, indicating that such structures must result from the non-linear processes. Thus, a first and important outcome of the linear ReOsc model fit is that the CMIP model mismatch to the observed phase space characteristics are likely to result from non-linear processes. Despite the limitation of the linear ReOsc model fit, we can still gain some understanding about the CMIP model ensemble spread and how different linear aspects of the ReOsc model affect the ENSO phase space statistics. The mean S , the probability of extreme S > 2 and the phase speed are all only varying along the diagonals of the phase space but have no variation along the T -axis and h -axis of the diagrams (left column in Fig. 8 ). The growth rates are only varying along T -axis and h -axis, but not along the diagonals. However, the actual data of the CMIP ensembles does show more complex variations in all statistical aspects (right column Fig. 5 ). This in turn suggests that variations not captured by the linear model parameters are a result of non-linear processes. We can further test the sensitivity of the ENSO phase space to the individual ReOsc model parameter variations in the CMIP ensembles, by varying only a subset of the parameters, while holding the other parameters at the ensemble mean values following the approach of VD18. This approach can allow us to determine which ReOsc model dynamics are causing asymmetries or variations in the ENSO phase space. The following three scenarios are presented: asymmetries in the growth rates (second column of Fig. 8 ), coupling (third column of Fig. 8 ), and forcing strength (right column of Fig. 8 ). The spread in the mean S values in the phase space of the linear ReOsc model fits (first row Fig. 8 ) is much smaller in each of the three sensitivity experiments, suggesting that the spread in mean S is resulting from a more complex combination of several parameters, but does not result from any of the three asymmetries tested. The ensemble mean S value is mostly affected by asymmetries in the coupling and forcings, which both have opposing effects. Here the asymmetries in linear coupling do have some similarities with the observed, but they lack the non-linear elements. Asymmetries in the growth rates have little impact on the mean S . The phase space variations of the extreme values ( S > 2) are affected by all their parameter asymmetries, but most strongly by the asymmetries in the coupling and forcings (second row Fig. 8 ). The phase space variations of the growth rates are clearly linked to the asymmetries in the growth rates, but asymmetries in the coupling and forcings do also have a significant effect (third row Fig. 8 ). The sensitivities of the phase speed to the linear ReOsc parameters show complex behaviour (last row Fig. 8 ). When we only consider the asymmetries in the growth rates or forcing strength, we get much larger phase speeds than observed in the CMIP ensemble. This suggests that the asymmetry in the coupling parameters has a strong impact on the phase speed. We can further notice that none of the three asymmetries reflect the overall variations in the phase speed, suggesting that more complex combinations of different ReOsc model parameters affect the phase speed. 3.4 Wind dynamics for CMIP models The surface wind response to SST anomalies is a key element of the Bjerknes feedback and is one of the main processes that control ENSO growth rate and period. While this relation does not directly relate to the ENSO phase space or the ReOsc model, it is implicitly related to both as it affects all parameters of the ReOsc model, in particular the growth rate of T . Figure 9 a shows the cross-correlation between the NINO4 region (5°S to 5°N, 160°E to150°W) zonal wind stress ( τ x ) and T for the observations and the CMIP models. The observed cross-correlation shows a strong and mostly in-phase relation between τ x and T , consistent with a strong response of τ x to variations in T . The relation is, however, strongest when τ x leads T by about one month, indicating also that variations in τ x are forcing variations in T . Overall, we see a close agreement between CMIP 5 and 6, and between the models and the observations. Despite the good agreement, there are some deviations in the CMIP models from observations. The CMIP models have a slightly low cross-correlation at lag zero, indicating a weaker response of τ x to T . The maximum cross-correlation in the CMIP ensembles is when τ x leads T by about two months, indicating a stronger forcing of τ x on variations in T than observed. The spread in these relations is quite large in the CMIP ensembles, indicating quite diverse behaviour in the individual models. Figure 9 b shows the regression value of τ x per T , which is often referred to as the atmospheric Bjerknes feedback (Bjerknes 1969 ). The observed atmospheric Bjerknes feedback is approximately 12.5*10 − 3 Nm − 2 / o C, whereas the CMIP ensemble averages are 7.4 (CMIP5) and 7.8*10 − 3 (CMIP6) Nm − 2 / o C. Thus, CMIP models underestimate this important feedback by approximately 40%, with none of the CMIP models reaching the observed value and no significant improvement from CMIP 5 to 6. There is a substantial spread within the CMIP model ensemble that is also of similar magnitude in both CMIP 5 and 6. These results are largely consistent with previous studies (Guilyardi et al. 2012 ; Bellenger et al. 2014 ; Vijayeta and Dommenget 2018 ; Guilyardi et al. 2020 ). To further examine the variability in T , τ x and h we looked at the standard deviation of these variables (Fig. 10 ). The standard deviation of T in the CMIP model averages closely resembles observations, with the CMIP5 model average being smaller than the CMIP6 value (Fig. 10 a). However, there are significant differences between individual CMIP models, with some models having significantly smaller and others larger standard deviations than the ensemble model or observed. The observed τ x has a significantly larger standard deviation than both CMIP5 and CMIP6 model averages (Fig. 10 b). This would be consistent with the weaker T forcing on τ x . Some CMIP models show standard deviations that are closer to the observation, while others have noticeably smaller standard deviations. The results suggest that the weaker τ x variations may be a cause or a reflection of the weaker τ x relation to T in the models. The CMIP6 ensemble mean standard deviation for h is very close to the observed, whereas the CMIP5 ensemble mean is a bit smaller, and some show large deviations for the ensemble mean (Fig. 10 c). However, more relevant for the ENSO phase space discussion is the variability of h shift (Fig. 10 d). This is significantly weaker than the observed h for both CMIP ensembles. Given a similar cross-correlation between h shift and T in the CMIP model as in the observations (Fig. 1 b) and the smaller standard deviation of h shift in the CMIP as shown above, we would expect larger values of a 12 in CMIP models than observed (see Fig. 7 e), since the regression values are divided the smaller standard deviations of h shift . Thus, per h shift anomaly in the CMIP models we get a larger T anomaly than in observations. It should be noted here that this does not need necessarily imply a larger sensitivity of T to h shift , because correlation does not allow to determine causality here. 4. Summary and discussions In this study, we examined how well the CMIP5 and CMIP6 models can simulate the observed ENSO phase space characteristics. In the context of the ReOsc model, the ENSO phase space presents ENSO as a cycle for which we can analyze important statistics, such as extreme value probabilities, growth rates, and phase speeds, as a function of the different phases of ENSO. This allows us to consider linear and non-linear aspects of the amplitude, growth rates, and phase transition speeds at the same time. A key starting point for the analysis of the ENSO phase space is the cross-correlation between T and h , which is assumed and observed to be out-of-phase (orthogonal) to allow for a proper evaluation of the phase space statistics, and it is an assumed relation for an idealized linear ReOsc model. We found that this relationship is clearly not out-of-phase for the CMIP models but is shifted to negative lag-zero cross-correlations. This shift in the cross-relation between T and h is linked to a regional shift in the ENSO patterns in CMIP models, with thermocline depth variability being shifted to the eastern equatorial Pacific. If the eastward shift in the thermocline depth variability ( h shift ) is considered in estimating the ENSO phase space, then we find a much better agreement between models and observations. This suggests that CMIP models do simulate the ENSO cycle with some regional displacements in the ocean-atmosphere interactions. The main characteristics of the observed phase space are mostly well simulated by the CMIP ensembles. The key feature of the observed higher likelihoods of extreme El Niño to discharge states and reduced likelihoods of extreme La Niña to recharge states are captured by the CMIP ensembles mean with a high correlation. However, the models strongly underestimate the magnitude of non-linearity in this aspect. This suggests that although models, in their ensemble mean, capture the non-linearity of phases of ENSO, the non-linear processes responsible for these characteristics are not strong enough (weaker non-linear processes in CMIP models than observation as suggested by An et al. 2005 ; Sun et al. 2016 ) or are too linear. The observed non-linearity in growth rates of the ENSO phase space are less well captured by the CMIP ensembles. While they do capture the enhanced growth rates during the El Niño states, they lack increased growth rates during the recharge state, and they lack strong decay (negative growth rates) of ENSO amplitudes in the transition from discharge to La Niña states. The observed variations in the ENSO phase speed are also captured by the CMIP ensembles, but only with low correlations. The models can capture the phase speed in the transition from El Niño to the discharge state, but the transition is too fast (or too clear) from the La Niña to the recharge state, where the observed transition is much slower (or less clear). Combined with the lack of variations in the growth rates this suggests that the ENSO cycle in CMIP models is too regular, with too little differences between the El Niño and La Niña phases. The ENSO phase space is strongly linked to the idea of the ReOsc model. It therefore does make sense to explain the observed and simulated characteristics of the phase space by the ReOsc model dynamics. We therefore fitted the linear ReOsc model parameters to the observed and simulated T and h data but encountered a number of problems. First, the regional shift in the ENSO pattern in CMIP models does strongly affect the ReOsc model estimates. Second, the most interesting aspects of the observed ENSO phase space must result from non-linear dynamics, where it is unclear what kind of dynamics these are (Dommenget et al. 2023 ). We further analyzed the relation between SST ( T ), wind stress ( τ x ) and thermocline depth ( h ) variability. Like previous studies, we find that the linear relation of τ x and T is much weaker than observed. This weaker relation is likely to contribute to the weaker non-linear ENSO phase space characteristics, since previous studies have indicated that the wind-SST relation is a key element of the non-linear behavior of ENSO. The variability of h in the eastward shifted region of the Pacific, which is most important for CMIP models, is also weaker than observed. This, combined with the weaker wind-SST relation indicates that the ENSO mode in the CMIP model is more wind-driven and less influenced by the thermocline variability. A previous study by Zhao and Sun 2022 found that the CMIP models have a shallower thermocline depth in the eastern equatorial Pacific than observed, which could explain the weaker response of the above-mentioned processes. However, further analysis into the SST, wind stress and thermocline depth interactions are needed to understand how they impact the ENSO phase space. It is worth noting that Planton et al. 2021 made no mention of thermocline depth biases in their ENSO metrics package, which is quite surprising, given that thermocline depth and their biases have a significant impact on the ENSO phase space characteristics and its dynamics as we found in our study. While the ensemble mean of the CMIP models does show good agreement with the observations in many aspects, there is a widespread within the model ensembles. Some models show opposite non-linear ENSO phase spaces from those observed, marking severe biases in the model behavior. Other models clearly outperform the ensemble mean, indicating that much more realistic ENSO simulations are possible within the CMIP model framework. In the comparison of the CMIP5 versus the CMIP6 ensembles, we find no substantial differences in the ENSO phase space dynamics, indicating no substantial improvements from CMIP5 to CMIP6. In particular, we find no differences in the skill of the model in simulating the significant non-linear ENSO behavior. This is a bit disappointing and worrying, as it indicates that no progress has been made from CMIP5 to CMIP6. This finding is largely consistent with Planton et al. 2021 . More future work is required to better understand the limitations of the CMIP models in simulating the non-linear ENSO phase space characteristics. A better understanding of how the ReOsc model can incorporate the non-linear aspects to simulate the observed ENSO phase space would be important to better evaluate CMIP models. The regional shift in the simulation of ENSO dynamics are also an important bias. We still don't know why these biases exist and what are the reason for these discrepancies, necessitating further analysis and a deeper understanding of this question, as well as opportunities for future studies to address these concerns. Declarations Competing Interests None Funding This research is funded by the Australian Research Council (ARC), discovery project “Improving projections of regional sea level rise and their credibility” (DP200102329) and the Centre of Excellence for Climate Extremes (CLEX) Grant Number: CE170100023. Authors Contributions Priyamvada Priya formulated the concept and designed this study. All authors contributed to the material preparation, collecting data, and analysis of this study. Priyamvada Priya wrote the first draft of the manuscript, and all authors commented and provided feedback on that. All authors read and approved the final version of the manuscript. Acknowledgements The authors would like to thank Monash University Melbourne, Australia, for providing all the resources and instruments required to conduct this research. This research is funded by the Australian Research Council (ARC) Centre of Excellence for Climate Extremes (CLEX), CLEX Grant Number: CE170100023, and the ARC discovery project “Improving projections of regional sea level rise and their credibility” (DP200102329). The authors would also like to thank Australia's computing facility National Computational Infrastructure (NCI) and CLEX for providing us with a wide range of facilities and computational tools for this study. Data Availability The data sources for this study are cited in the text. References Alizadeh O (2022) A review of the El Niño-Southern Oscillation in future. Earth Sci Rev 104246. https://doi.org/10.1016/j.earscirev.2022.104246 An SI, Jin FF (2004) Nonlinearity and asymmetry of ENSO. J Clim 17(12):2399–2412. https://doi.org/10.1175/1520-0442(2004)0172.0.CO;2 An SI, Ham YG, Kug JS, Jin FF, Kang IS (2005) El Niño–La Niña asymmetry in the coupled model intercomparison project simulations. J Clim 18(14):2617–2627. https://doi.org/10.1175/JCLI3433.1 An SI, Tziperman E, Okumura YM, Li T (2020) : ENSO irregularity and asymmetry. El Niño Southern Oscillation in a changing climate , 153–172, https://doi.org/10.1002/9781119548164.ch7 Bayr T, Domeisen DI, Wengel C (2019) The effect of the equatorial Pacific cold SST bias on simulated ENSO teleconnections to the North Pacific and California. Clim Dyn 53:3771–3789. https://doi.org/10.1007/s00382-019-04746-9 Bayr T, Latif M (2023) ENSO atmospheric feedbacks under global warming and their relation to mean-state changes. Clim Dyn 60(9–10):2613–2631. https://doi.org/10.1007/s00382-022-06454-3 Bellenger H, Guilyardi É, Leloup J, Lengaigne M, Vialard J (2014) ENSO representation in climate models: From CMIP3 to CMIP5. Clim Dyn 42(7):1999–2018. https://doi.org/10.1007/s00382-013-1783-z Beobide-Arsuaga G, Bayr T, Reintges A, Latif M (2021) Uncertainty of ENSO-amplitude projections in CMIP5 and CMIP6 models. Clim Dyn 56:3875–3888. https://doi.org/10.1007/s00382-021-05673-4 Bjerknes J (1969) : Atmospheric teleconnections from the equatorial Pacific. Monthly weather review , 97 (3), 163–172, https://doi.org/10.1175/1520-0493 (1969)097%3C0163:ATFTEP%3E2.3.CO;2 Burgers G, Stephenson DB (1999) The normality of el niño. Geophys Res Lett 26(8):1027–1030. https://doi.org/10.1029/1999GL900161 Burgers G, Jin FF, Van Oldenborgh GJ (2005) The simplest ENSO recharge oscillator. Geophys Res Lett 32(13). https://doi.org/10.1029/2005GL022951 Capotondi A, Wittenberg AT, Newman M, Di Lorenzo E, Yu JY, Braconnot P, … and, Yeh SW (2015) Understanding ENSO diversity. Bull Am Meteorol Soc 96(6):921–938. https://doi.org/10.1175/BAMS-D-13-00117.1 Carton JA, Giese BS (2008) A reanalysis of ocean climate using Simple Ocean Data Assimilation (SODA). Mon Weather Rev 136(8):2999–3017. https://doi.org/10.1175/2007MWR1978.1 Chen M, Li T (2021) ENSO evolution asymmetry: EP versus CP El Niño. Clim Dyn 56(11):3569–3579. https://doi.org/10.1007/s00382-021-05654-7 Crespo LR, Rodríguez-Fonseca MB, Polo I, Keenlyside N, Dommenget D (2022) Multidecadal variability of ENSO in a recharge oscillator framework. Environ Res Lett 17(7). p.074008, 10.1088/1748–9326/ac72a3 Dee DP, Uppala S (2009) Variational bias correction of satellite radiance data in the ERA-Interim reanalysis. Q J Royal Meteorological Society: J atmospheric Sci Appl Meteorol Phys Oceanogr 135(644):1830–1841. https://doi.org/10.1002/qj.493 Dommenget D, Bayr T, Frauen C (2013) Analysis of the non-linearity in the pattern and time evolution of El Niño southern oscillation. Clim Dyn 40:2825–2847. https://doi.org/10.1007/s00382-012-1475-0 Dommenget D, Vijayeta A (2019) Simulated future changes in ENSO dynamics in the framework of the linear recharge oscillator model. Clim Dyn 53(7):4233–4248. https://doi.org/10.1007/s00382-019-04780-7 Dommenget D, Al-Ansari M (2022) Asymmetries in the ENSO phase space. Clim Dyn 1–20. https://doi.org/10.1007/s00382-022-06392-0 Dommenget D, Priya P, Vijayeta A (2023) ENSO phase space dynamics with an improved estimate of the thermocline depth. Clim Dyn. https://doi.org/10.1007/s00382-023-06883-8 Eyring V, Bony S, Meehl GA, Senior CA, Stevens B, Stouffer RJ, Taylor KE (2016) Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci Model Dev 9(5):1937–1958. https://doi.org/10.5194/gmd-9-1937-2016 Frauen C, Dommenget D (2010) El Niño and La Niña amplitude asymmetry caused by atmospheric feedbacks. Geophys Res Lett 37(18). https://doi.org/10.1029/2010GL044444 Good SA, Martin MJ, Rayner NA (2013) EN4: Quality controlled ocean temperature and salinity profiles and monthly objective analyses with uncertainty estimates. J Geophys Research: Oceans 118(12):6704–6716. 10.1002/2013JC009067 Guilyardi E, Bellenger H, Collins M, Ferrett S, Cai W, Wittenberg A (2012) A first look at ENSO in CMIP5. Clivar Exchanges 17(58):29–32 Guilyardi E, Capotondi A, Lengaigne M, Thual S, Wittenberg AT (2020) : ENSO modeling: History, progress, and challenges. El Niño Southern Oscillation in a changing climate , 199–226, https://doi.org/10.1002/9781119548164.ch9 Hersbach H, Bell B, Berrisford P, Hirahara S, Horányi A, Muñoz-Sabater J, … and, Thépaut JN (2020) The ERA5 global reanalysis. Q J R Meteorol Soc 146(730):1999–2049. https://doi.org/10.1002/qj.3803 Jiang W, Huang P, Huang G, Ying J (2021) Origins of the excessive westward extension of ENSO SST simulated in CMIP5 and CMIP6 models. J Clim 34(8):2839–2851. https://doi.org/10.1175/JCLI-D-20-0551.1 Jin FF (1997) : An equatorial ocean recharge paradigm for ENSO. Part I: Conceptual model. Journal of the atmospheric sciences , 54 (7), 811–829, https://doi.org/10.1175/1520-0469 (1997)054%3C0811:AEORPF%3E2.0.CO;2 Kessler WS (2002) Is ENSO a cycle or a series of events? Geophys Res Lett 29(23) 40 – 1. https://doi.org/10.1029/2002GL015924 Kim ST, Yu JY (2012) The two types of ENSO in CMIP5 models. Geophys Res Lett 39(11). https://doi.org/10.1029/2012GL052006 Kim ST, Cai W, Jin FF, Yu JY (2014) ENSO stability in coupled climate models and its association with mean state. Clim Dyn 42:3313–3321. https://doi.org/10.1007/s00382-013-1833-6 Latif M, Sperber K, Arblaster J, Braconnot P, Chen D, Colman A, … and, Zebiak S (2001) ENSIP: the El Nino simulation intercomparison project. Clim Dyn 18:255–276. https://doi.org/10.1007/s003820100174 Latif M, Keenlyside NS (2009) : El Niño/Southern Oscillation response to global warming. Proceedings of the National Academy of Sciences , 106 (49), 20578–20583, https://doi.org/10.1073/pnas.0710860105 McKenna S, Santoso A, Gupta AS, Taschetto AS, Cai W (2020) Indian Ocean Dipole in CMIP5 and CMIP6: characteristics, biases, and links to ENSO. Sci Rep 10(1):11500. https://doi.org/10.1038/s41598-020-68268-9 McPhaden MJ, Zebiak SE, Glantz MH (2006) ENSO as an integrating concept in earth science. Science 314(5806):1740–1745. https://doi.org/10.1126/science.1132588 Meinen CS, McPhaden MJ (2000) Observations of warm water volume changes in the equatorial Pacific and their relationship to El Niño and La Niña. J Clim 13(20):3551–3559. https://doi.org/10.1175/1520-0442(2000)0132.0.CO;2 Okumura YM (2019) ENSO diversity from an atmospheric perspective. Curr Clim Change Rep 5:245–257. https://doi.org/10.1007/s40641-019-00138-7 Planton YY, Guilyardi E, Wittenberg AT, Lee J, Gleckler PJ, Bayr T, … and, Voldoire A (2021) Evaluating climate models with the CLIVAR 2020 ENSO metrics package. Bull Am Meteorol Soc 102(2):E193–E217. https://doi.org/10.1175/BAMS-D-19-0337.1 Priya P, Dommenget D, McGregor S (2024) The dynamics of the El Niño–Southern Oscillation diversity in the recharge oscillator framework. Clim Dyn 1–21. https://doi.org/10.1007/s00382-024-07158-6 Rayner NAA, Parker DE, Horton EB, Folland CK, Alexander LV, Rowell DP, Kent EC, Kaplan A (2003) Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. J Geophys Res Atmos 108(D14). https://doi.org/10.1029/2002JD002670 Sun Y, Wang F, Sun DZ (2016) Weak ENSO asymmetry due to weak nonlinear air–sea interaction in CMIP5 climate models. Adv Atmos Sci 33:352–364. https://doi.org/10.1007/s00376-015-5018-6 Taylor KE, Stouffer RJ, Meehl GA (2012) An overview of CMIP5 and the experiment design. Bull Am Meteorol Soc 93(4):485–498. https://doi.org/10.1175/BAMS-D-11-00094.1 Timmermann A, An SI, Kug JS, Jin FF, Cai W, Capotondi A, … and, Zhang X (2018) El Niño–southern oscillation complexity. Nature 559(7715):535–545. https://doi.org/10.1038/s41586-018-0252-6 Tsonis AA, Hunt AG, Elsner JB (2003) On the relation between ENSO and global climate change. Meteorol Atmos Phys 84(3–4) 229–242, 3579, 3579. https://doi.org/10.1007/s00703-003-0001-7 Vijayeta A, Dommenget D (2018) An evaluation of ENSO dynamics in CMIP simulations in the framework of the recharge oscillator model. Clim Dyn 51(5):1753–1771. https://doi.org/10.1007/s00382-017-3981-6 Wang G, Dommenget D, Frauen C, Wang (2015) An evaluation of the CMIP3 and CMIP5 simulations in their skill of simulating the spatial structure of SST variability. Clim Dyn 44:95–114. https://doi.org/10.1007/s00382-014-2154-0 Wang C, Deser C, Yu JY, DiNezio P, Clement A (2017) El Niño and southern oscillation (ENSO): a review. Coral reefs East Trop Pacific: Persistence loss dynamic Environ 85–106. https://doi.org/10.1007/978-94-017-7499-4_4 Weller E, Cai W (2013) Asymmetry in the IOD and ENSO teleconnection in a CMIP5 model ensemble and its relevance to regional rainfall. J Clim 26(14):5139–5149. https://doi.org/10.1175/JCLI-D-12-00789.1 Wengel C, Dommenget D, Latif M, Bayr T, Vijayeta A (2018) What controls ENSO-amplitude diversity in climate models? Geophys Res Lett 45(4):1989–1996. https://doi.org/10.1002/2017GL076849 Wyrtki K (1985) Water displacements in the Pacific and the genesis of El Niño cycles. J Geophys Research: Oceans 90(C4):7129–7132. https://doi.org/10.1029/JC090iC04p07129 Yang C, Giese BS (2013) El Niño Southern Oscillation in an ensemble ocean reanalysis and coupled climate models. J Geophys Research: Oceans 118(9):4052–4071. https://doi.org/10.1002/jgrc.20284 Yu Y, Dommenget D, Frauen C, Wang G, Wales S (2016) ENSO dynamics and diversity resulting from the recharge oscillator interacting with the slab ocean. Clim Dyn 46(5):1665–1682. https://doi.org/10.1007/s00382-015-2667-1 Zuo H, Balmaseda MA, Tietsche S, Mogensen K, Mayer M (2019) The ECMWF operational ensemble reanalysis–analysis system for ocean and sea ice: a description of the system and assessment. Ocean Sci 15(3):779–808. https://doi.org/10.5194/os-2018-154 Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4727039","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":326358611,"identity":"fb71934f-1596-475e-a7ad-215f985e7151","order_by":0,"name":"Priyamvada Priya","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABKklEQVRIiWNgGAWjYFACxgYIfQBEGDAw8CO4CURqkWwgqAUGDkBpgwMEtJhLH26T+JlzWJ7v+AG2xxUFh+WNrx0+9uDHnzsM/Ow5Bti0WPYltkn2bjtsOPNMArvhGYPDhttup6Ub9rY9Y5DseYNVi8EZxjYJ3m1pjBsOJLBJNhgcZtx2O8dMgrfhMIPBDey2gLRI/t2WZr/h/AOwFvvNs3PMJP/8Ocxgj0eLNO82m8QNNyC2JG6QzjGT5mED2iKBwy89jM3WsttskmfeeNhu2GCQnjzjdlqatGzbYR6JM88KsIYYD/vDm2+3Sdj2nU8+9rDhj7Vt/+zkY5Jv/hyW429P3oDVYQwMLBIQJmMbkGiGy/BgUw7VwvwBymYD4jpcCkfBKBgFo2AEAwDCJ2rfv/TS8wAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-1285-8093","institution":"Monash University","correspondingAuthor":true,"prefix":"","firstName":"Priyamvada","middleName":"","lastName":"Priya","suffix":""},{"id":326358612,"identity":"8f99dccd-69ee-49ab-bb12-275b37e38194","order_by":1,"name":"Dietmar Dommenget","email":"","orcid":"","institution":"Monash University","correspondingAuthor":false,"prefix":"","firstName":"Dietmar","middleName":"","lastName":"Dommenget","suffix":""}],"badges":[],"createdAt":"2024-07-11 23:47:56","currentVersionCode":1,"declarations":{"humanSubjects":false,"vertebrateSubjects":false,"conflictsOfInterestStatement":false,"humanSubjectEthicalGuidelines":false,"humanSubjectConsent":false,"humanSubjectClinicalTrial":false,"humanSubjectCaseReport":false,"vertebrateSubjectEthicalGuidelines":false},"doi":"10.21203/rs.3.rs-4727039/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4727039/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":61671557,"identity":"3ee9ad41-33bf-49c4-a876-f0ad2db053c8","added_by":"auto","created_at":"2024-08-02 18:41:45","extension":"jpeg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":259028,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003e(a) Cross-Correlation of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e vs. \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e and (b) Cross-Correlation of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e vs. \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003eshift\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e, with \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e defined for equatorial pacific region (5°S to 5°N, 130°E to 80°W) and \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003eshift\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e for east equatorial pacific region (5°S to 5°N, 190°E to 80°W). The observed (black) values in both panels are based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e. CMIP5 (blue), CMIP6 (red), CMIP5 individual models (light blue), and CMIP6 individual models (light red)\u003c/em\u003e\u003c/p\u003e","description":"","filename":"floatimage1.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/7b6cb860c0c12ab69f4853dd.jpeg"},{"id":61671313,"identity":"a686f786-73d9-42e3-8c20-a7508499797e","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"jpeg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":793100,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eEmpirical Orthogonal Function (EOF) mode-1 patterns of SST for (a) Observation, (b) ACCESS-CM2 CMIP6 model, (c) MRI-ESM2-0 CMIP6 model, (d) NorESM2-MM CMIP6 model, and EOF mode-2 patterns of thermocline depth for (e) Observation, (f) CMCC-CM2-SR5 CMIP6 model, (g) ACCESS-ESM1-5 CMIP6 model, (h) CNRM-CM6-1 CMIP6 model. The black dashed line represents the zonal transition lines for zero equatorial anomalies for observed SST and thermocline depth, and the blue dashed lines for models, respectively\u003c/em\u003e\u003c/p\u003e","description":"","filename":"floatimage2.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/16054764ad7afe90943ca774.jpeg"},{"id":61671312,"identity":"fc771e44-6511-4f5c-a866-62a9bff83345","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"jpeg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":120410,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eThe longitudes of zonal transition lines for zero equatorial anomalies of the (a) EOF mode-1 pattern for SST, and (b) EOF mode-2 pattern for thermocline depth (\u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e) for observation and CMIP models. Observed (black), CMIP5 (blue), CMIP6 (red), CMIP5 individual models (light blue), and CMIP6 individual models (light red)\u003c/em\u003e\u003c/p\u003e","description":"","filename":"floatimage3.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/2c997ac251912cc7bdcce327.jpeg"},{"id":61671316,"identity":"9ba2bd03-9716-4457-947a-bd7f535f4e31","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"jpeg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":890646,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eENSO phase space: (a) observed, left column three CMIP6 models based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e and right column based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003eshift\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e: (b, e) NorESM2-MM, (c, f) CESM2-WACCM, (d, g) BCC-ESM1. All data (grey), mean as a function of phase CMIP (red) and observed (black). The black arrows represent the tendency vectors of the mean tendencies of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003en\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e and \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003en\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e\u003cstrong\u003e \u003c/strong\u003e\u003c/em\u003e\u003cem\u003eover the range of ±0.4, with the reference point (origin of the vector) in phase space. A unit-length vector represents a tendency of 1mon\u003c/em\u003e\u003csup\u003e\u003cem\u003e-1\u003c/em\u003e\u003c/sup\u003e\u003cem\u003e, with the vector’s scale proportional to the magnitude of the tendencies. Vectors for \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eS \u003c/strong\u003e\u003c/em\u003e\u003cem\u003e\u0026lt;3.5 are only displayed here when the data is available. The correlation coefficients between CMIP6 models and observation are given for each panel (red numbers)\u003c/em\u003e\u003c/p\u003e","description":"","filename":"floatimage5.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/14fb0925ee9bb782fec86607.jpeg"},{"id":61671321,"identity":"50a966cd-7e03-4dfd-ac3b-f4c714313c1c","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"jpeg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":794065,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eStatistics of the ENSO phase space: (a, e) mean \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eS\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e, (b, f) probability distribution (\u003c/em\u003e\u003cem\u003e\u003cstrong\u003eS\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e \u0026gt;2), (c, g) growth rate, (d, h) phase speed. The left column is based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e (5°S to 5°N, 130°E to 80°W), and the right column is based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003eshift\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e (5°S to 5°N, 190°E to 80°W). The observed (black) values in both columns are based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e. Ensemble mean CMIP5 (blue), CMIP6 (red), and individual models CMIP5 (light blue), and CMIP6 (light red). In (c) and (g), the grey circle represents zero growth rates, values inside the circle represent negative, and values outside denote positive growth rates. The correlation between observation and CMIP5 (blue numbers), and CMIP6 (red numbers) are given for all diagrams\u003c/em\u003e\u003c/p\u003e","description":"","filename":"floatimage7.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/3f5fc98c5f97095bd79056be.jpeg"},{"id":61671318,"identity":"93d10b88-dc70-4b44-aaba-7f5b7d38759e","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"jpeg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":148153,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003ePower spectrum of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e for observed and CMIP models. Observed (black), CMIP5 (blue), CMIP6 (red), CMIP5 individual models (light blue), and CMIP6 individual models (light red). Dotted black lines are representing particular periods\u003c/em\u003e\u003c/p\u003e","description":"","filename":"floatimage8.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/1a2dd41eb7065d89a54354e5.jpeg"},{"id":61671558,"identity":"2d7ed038-9059-4607-bc60-97a9c2b36b3d","added_by":"auto","created_at":"2024-08-02 18:41:45","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":260875,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eParameters of the ReOsc model: (a) growth rates of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e and \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e (b) \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e coupling to \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e vs. \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e coupling \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e, (c) noise forcing of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e and \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e for observation and CMIP models. The left column is based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e and the right panel is based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003eshift\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e. The observed (black) values in both columns are based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e. The Wyrtki-index is shown as colored solid lines in panels b and e for specific periods assuming a\u003c/em\u003e\u003csub\u003e\u003cem\u003e11\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e = a\u003c/em\u003e\u003csub\u003e\u003cem\u003e22\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e = -0.06. Uncertainty in the observed parameters for 95% confidence interval is shown as a grey box. CMIP5 (blue), CMIP6 (red), CMIP5 individual models (light blue), CMIP6 individual models (light red), and one-to-one line (grey dashed)\u003c/em\u003e\u003c/p\u003e","description":"","filename":"fig7.png","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/43434a5fec3f3d3a3c913af0.png"},{"id":61671317,"identity":"01b4c53c-7c20-415c-aa2c-64f6dc3e53ba","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":170755,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eENSO phase space statistics of different linear ReOsc model fitted to CMIP models compared against the observed: (first row) mean \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eS\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e, (second row) probability distribution (\u003c/em\u003e\u003cem\u003e\u003cstrong\u003eS\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e\u0026gt;2), (third row) growth rate, (last row) phase speed. First column: complete ReOsc model fitted, and other columns for ReOsc models only considering asymmetries in the fitted parameters of: (second column) growth rates (a\u003c/em\u003e\u003csub\u003e\u003cem\u003e11\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e, a\u003c/em\u003e\u003csub\u003e\u003cem\u003e22\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e), (third column) coupling (a\u003c/em\u003e\u003csub\u003e\u003cem\u003e12\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e, a\u003c/em\u003e\u003csub\u003e\u003cem\u003e21\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e), and (last column) noise forcing strength (x\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e, x\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e). The observed (black) values in all columns are based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e. CMIP5 (blue), CMIP6\u003c/em\u003e\u003cem\u003e\u003cstrong\u003e \u003c/strong\u003e\u003c/em\u003e\u003cem\u003e(red), CMIP5 individual models (light blue), and CMIP6 individual models (light red). The grey circles in the third row panels indicate zero growth rates, values inside the circle indicate negative growth rates, and values outside the circle indicate positive growth rates. In the last row, the \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eS\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e axes scale varies to better highlight the model variations\u003c/em\u003e\u003c/p\u003e","description":"","filename":"Onlinefloatimage12.png","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/72de2a02b946e2026a6ce3ea.png"},{"id":61671320,"identity":"17082fda-f1a7-437e-a6f6-02e76efee4e7","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":120624,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003e(a) Cross-Correlation of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e vs. tau\u003c/em\u003e\u003csub\u003e\u003cem\u003ex\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e for NINO4 region (5°S to 5°N, 160°E to150°W), and (b) Atmospheric Bjerknes feedback (regression values of \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e vs. tau\u003c/em\u003e\u003csub\u003e\u003cem\u003ex\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e) for observed and CMIP models. Observed (black), CMIP5 (blue), CMIP6 (red), CMIP5 individual models (light blue), and CMIP6 individual models (light red)\u003c/em\u003e\u003c/p\u003e","description":"","filename":"Onlinefloatimage13.png","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/4bd74f7caf55e5df78b1658b.png"},{"id":61671319,"identity":"2ec336ed-0b7a-4d69-9e35-8bb0d4c9c120","added_by":"auto","created_at":"2024-08-02 18:33:45","extension":"png","order_by":10,"title":"Figure 10","display":"","copyAsset":false,"role":"figure","size":43830,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eStandard deviation of (a) \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eT\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e, (b) tau\u003c/em\u003e\u003csub\u003e\u003cem\u003ex \u003c/em\u003e\u003c/sub\u003e\u003cem\u003efor NINO4 region (5°S to 5°N, 160°E to150°W), (c) \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e, (d) \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003csub\u003e\u003cem\u003e\u003cstrong\u003eshift\u003c/strong\u003e\u003c/em\u003e\u003c/sub\u003e\u003cem\u003e. The observed values in (c) and (d) are based on \u003c/em\u003e\u003cem\u003e\u003cstrong\u003eh\u003c/strong\u003e\u003c/em\u003e\u003cem\u003e for reference. Observed (black), CMIP5 (blue), CMIP6 (red), CMIP5 individual models (light blue), CMIP6 individual models (light red), and observation reference line (black dashed)\u003c/em\u003e\u003c/p\u003e","description":"","filename":"Onlinefloatimage14.png","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/53692ca64c5cc2ceb0019158.png"},{"id":65049444,"identity":"a565d616-4b2d-4e49-8123-b44b22bf423b","added_by":"auto","created_at":"2024-09-23 05:37:37","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4637965,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4727039/v1/edfc246b-9bb4-42c0-853e-692adbec5a7a.pdf"}],"financialInterests":"","formattedTitle":"ENSO phase space dynamics in CMIP models","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eThe El Ni\u0026ntilde;o-Southern Oscillation (ENSO) is a climate phenomenon that influences the global climate and weather patterns, ecosystems, and society (Tsonis et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; McPhaden et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2006\u003c/span\u003e; Latif and Keenlyside \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Timmermann et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Alizadeh \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). It is characterized by its two distinct phases which are the warm phase (El Ni\u0026ntilde;o) and cold phase (La Ni\u0026ntilde;a) in the equatorial east Pacific Ocean (Chen and Li \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The variations of ENSO are controlled by complex dynamics that include feedback from the atmosphere and ocean (Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Kim et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Wang et al. \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Vijayeta and Dommenget \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The aim of this study is to contribute to a better understanding of these complex dynamics.\u003c/p\u003e \u003cp\u003eKey ENSO dynamics can be determined by the interaction of ocean-atmosphere processes such as sea surface temperature (SST), wind stress, and thermocline depth, known as Bjerknes feedback (Bjerknes \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e1969\u003c/span\u003e). A simplified concept for combining these processes and for understanding ENSO dynamics is described by the recharge oscillator (ReOsc) model (Jin \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e1997\u003c/span\u003e; Burgers et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2005\u003c/span\u003e). The ReOsc model describes the variability of ENSO as an interaction between the tendencies of SST anomaly in the eastern equatorial Pacific (\u003cb\u003eT\u003c/b\u003e) and equatorial mean thermocline depth anomaly (\u003cb\u003eh\u003c/b\u003e), with the latter representing a proxy for the subsurface equatorial heat content. It has been used extensively in previous studies to analyse ENSO dynamics (Jin \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e1997\u003c/span\u003e; Burgers et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Frauen and Dommenget \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Yu et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Vijayeta and Dommenget \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Wengel et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Dommenget and Vijayeta \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Crespo et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Dommenget and Al-Ansari \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Priya et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe idealised ReOsc model is an oscillation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e, presenting a 2-dimensional phase space, as known from physics for a harmonic oscillator or complex system analysis (Meinen and McPhaden \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2000\u003c/span\u003e; Kessler \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2002\u003c/span\u003e; Dommenget and Al-Ansari \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Priya et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). Detailed analysis of the characteristics of the ENSO phase can help to understand how the dynamics of ENSO control the amplitude, non-linear behaviour, and time evolution as functions of the ENSO phases.\u003c/p\u003e \u003cp\u003eThe study by Dommenget et al. (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) showed that the ENSO phase space highlights strong non-linear aspects, with enhanced extremes in the El Ni\u0026ntilde;o to discharge states, larger growth rates in the transition from recharge to El Ni\u0026ntilde;o state, strong decay after the discharge states, faster phase speeds after El Ni\u0026ntilde;o states and slow phase speeds after La Ni\u0026ntilde;a states. These non-linear aspects of the ENSO cycle have also been reported in other studies, using different measures (An and Jin \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; An et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Dommenget et al \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Sun et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Okumura \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; An et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe Coupled Model Intercomparison Projects (CMIP) models are widely used to simulate and forecast ENSO activity. They are highly capable of simulating many aspects of ENSO and do show some improvements in the consecutive generations of CMIP from phases 3 to 6 in metrics used to evaluate ENSO, such as seasonal cycle, ENSO patterns, and diversity of ENSO events (Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Kim et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Despite their capabilities, the models do show some biases if compared against observations and they do show large variations within the model ensemble (Kim and Yu \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Guilyardi et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Vijayeta and Dommenget \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Bayr et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; McKenna et al. \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). The CMIP models exhibit biases in main atmospheric and oceanic feedback, which do control ENSO amplitude and stability (Guilyardi et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Kim et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Vijayeta and Dommenget \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Bayr et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Beobide-Arsuaga et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Bayr and Latif (\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) discussed how atmospheric feedback, specifically zonal wind feedback and heat flux feedback, are underestimated in CMIP model (CMIP5 and CMIP6) simulations of ENSO dynamics.\u003c/p\u003e \u003cp\u003eFurther, the models exhibit biases in their representation of the SST patterns, with a westward extension compared to observed data in CMIP5 and CMIP6 models (Weller and Cai \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Yang and Giese \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Capotondi et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Wang et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Guilyardi et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Jiang et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The models also exhibit biases in the equatorial thermocline depth, intensity, and zonal shape due to biases in the equatorial zonal wind stress, which affects upwelling and the zonal tilt of the equatorial thermocline (Guilyardi et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAbout 80% of the variations in simulating ENSO amplitude within the CMIP model ensemble are due to variations in the linear growth rates of \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e, and in the strength of the stochastic noise forcing, with the latter also potentially being linked to non-linear processes (Wengel et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) (non-linear processes such as non-linear air-sea coupling, non-linear dynamical heating, non-linear ocean advection, non-linear relationship of precipitate-SST feedback, heat flux because of the tropical ocean instability waves, atmospheric non-linear response to the SST anomaly (wind-SST relation) (An and Jin \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; An et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Sun et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Okumura \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; An et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2020\u003c/span\u003e)). Several studies on ENSO non-linearities point out that CMIP models tend to underestimate the non-linear aspects of ENSO (An and Jin \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; An et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Dommenget et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Sun et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Okumura \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn this study, we analyze the ENSO phase space statistics as simulated in the CMIP 5 and 6 simulations and compare them against the observed values. We examine how well the CMIP models can reproduce the observed ENSO phase space asymmetries to gain some understanding of how well models can simulate the non-linear behavior of ENSO as a function of the different phases of ENSO. We explore potential limitations in model simulations that may affect the model performance. Additionally, we are also interested in assessing the improvements from CMIP5 to CMIP6 models in representing the ENSO phase space characteristics.\u003c/p\u003e \u003cp\u003eThis paper is structured as follows: Section 2 describes the observational datasets used, CMIP simulations, phase space statistics methods, and the ReOsc model. This is followed by the main results section, in which we present the ENSO phase space statistics for CMIP simulations and observations. Here we also analyze the linear ReOsc model fit to CMIP models and discuss some aspects of the wind, thermocline, and SST cross-relations in the context of ENSO dynamics. The study is concluded with a summary and discussion section.\u003c/p\u003e"},{"header":"2. Data and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Observational datasets\u003c/h2\u003e \u003cp\u003eObservational datasets for SST and the thermocline depth estimates have been taken from the Hadley Centre global sea Ice coverage and SST (HadISST) (Rayner et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2003\u003c/span\u003e), UK Met Office Hadley Centre 'EN' series observations datasets, version 4 (EN4) (Good et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) respectively from 1980 to 2021. The reanalysis datasets for SST and thermocline depth estimates are taken from Simple Ocean Data Assimilation version 3 (SODA3) (Carton and Giese \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2008\u003c/span\u003e) (1980\u0026ndash;2017), Ocean Reanalysis System 5 (ORAS5) (1980\u0026ndash;2018) (Zuo et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2019\u003c/span\u003e), and CMCC Historical Ocean Reanalysis system: CHOR AS and CHOR RL (1980\u0026ndash;2010) (Yang et al. 2017). To observe the variability more precisely, these 5 different datasets are concatenated to create a considerable time series of SST and thermocline depth repeated for the same year. Anomalies are calculated by removing the mean seasonal cycle (calculated over the entire period) from the data. \u003cb\u003eT\u003c/b\u003e is the SST anomaly calculated for the NINO3 region (5\u0026deg;S to 5\u0026deg;N, 150\u0026deg;W to 90\u0026deg;W) and considered for the same in this study.\u003c/p\u003e \u003cp\u003ePrevious studies often use the depth of the 20\u003csup\u003eo\u003c/sup\u003eC (\u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003e20\u003c/b\u003e\u003c/sub\u003e) as a proxy for thermocline depth (\u003cb\u003eh\u003c/b\u003e), however, according to Dommenget et al. (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), the maximum gradient in the temperature profile (\u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003emxg\u003c/b\u003e\u003c/sub\u003e) is a better proxy of \u003cb\u003eh\u003c/b\u003e, which is calculated by averaging \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003emxg\u003c/b\u003e\u003c/sub\u003e over the equatorial Pacific region (5\u0026deg;S to 5\u0026deg;N, 130\u0026deg;E to 80\u0026deg;W), for the analysis of the ENSO dynamics within the framework of the ReOsc model than the widely used \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003e20\u003c/b\u003e\u003c/sub\u003e. The ENSO phase space for \u003cb\u003eT\u003c/b\u003e index by using \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003emxg\u003c/b\u003e\u003c/sub\u003e as a proxy of \u003cb\u003eh\u003c/b\u003e is closer to the idealized recharge oscillator model than \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003e20\u003c/b\u003e\u003c/sub\u003e. This motivates us to use \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003emxg\u003c/b\u003e\u003c/sub\u003e to estimate the equatorial Pacific Ocean thermocline depth \u003cb\u003eh\u003c/b\u003e, following the method of Dommenget et al. (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe zonal wind stress data are calculated from 10-m surface zonal wind datasets taken from the European Centre for Medium-Range Weather Forecasts (ECMWF-ERA5) reanalysis (1980\u0026ndash;2021) (Hersbach et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) to analyse the wind dynamics of ENSO. Here, we only consider wind stress data from atmospheric reanalysis datasets, ERA5, following a similar approach to previous papers (Guilyardi et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Vijayeta and Dommenget \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Guilyardi et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eCMIP5 and CMIP6 models with the corresponding model numbers\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"8\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCMIP5\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCMIP6\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/th\u003e \u003c/tr\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eModel No.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eModel Name\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eModel No.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eModel Name\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eModel No.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eModel Name\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eModel No.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003eModel Name\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e1.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eACCESS1-0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e21.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eGISS-E2-H-CC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e1.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eACCESS-CM2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e21.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eEC-Earth3-Veg\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e2.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eACCESS1-3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e22.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eGISS-E2-R\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e2.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eACCESS-ESM1-5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e22.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eEC-Earth3-Veg-LR\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e3.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ebcc-csm1-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e23.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eGISS-E2-R-CC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e3.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eBCC-ESM1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e23.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eFGOALS-f3-L\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e4.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ebcc-csm1-1-m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e24.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHadCM3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e4.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eBCC-CSM2-MR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e24.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eFGOALS-g3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e5.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCanESM2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e25.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHadGEM2-CC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e5.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCanESM5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e25.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eFIO-ESM-2-0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e6.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCCSM4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e26.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHadGEM2-ES\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e6.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCESM2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e26.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eGFDL-CM4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e7.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCESM1-BGC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e27.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eIPSL-CM5A-LR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e7.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCESM2-FV2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e27.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eGISS-E2-1-G\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e8.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCESM1-CAM5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e28.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eIPSL-CM5A-MR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e8.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCESM2-WACCM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e28.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eGISS-E2-1-G-CC\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e9.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCESM1-FASTCHEM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e29.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eIPSL-CM5B-LR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e9.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCESM2-WACCM-FV2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e29.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eINM-CM4-8\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e10.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCESM1-WACCM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e30.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMIROC5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e10.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCIESM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e30.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eINM-CM5-0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e11.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCMCC-CESM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e31.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMPI-ESM-LR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e11.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCMCC-CM2-SR5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e31.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eMIROC-ES2L\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e12.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCMCC-CM\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e32.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMPI-ESM-MR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e12.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCMCC-ESM2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e32.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eMPI-ESM1-2-LR\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e13.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCMCC-CMS\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e33.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMPI-ESM-P\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e13.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCNRM-CM6-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e33.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eMPI-ESM-1-2-HAM\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e14.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCNRM-CM5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e34.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMRI-CGCM3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e14.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eCNRM-ESM2-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e34.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eMRI-ESM2-0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e15.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCNRM-CM5-2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e35.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMRI-ESM1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e15.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eE3SM-1-0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e35.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eNESM3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e16.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCSIRO-Mk3L-1-2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e36.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eNorESM1-M\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e16.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eE3SM-1-1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e36.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eNorESM2-MM\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e17.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGFDL-CM3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cb\u003e37.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eNorESM1-ME\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e17.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eE3SM-1-1-ECA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e37.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eSAM0-UNICON\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e18.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGFDL-ESM2G\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e18.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eEC-Earth3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e38.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eTaiESM1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e19.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGFDL-ESM2M\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e19.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eEC-Earth3-AerChem\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cb\u003e39.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eUKESM1-0-LL\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003e20.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eGISS-E2-H\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cb\u003e20.\u003c/b\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eEC-Earth3-CC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/td\u003e \u003ctd align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.1 CMIP Simulations\u003c/h2\u003e \u003cp\u003eFor the CMIP simulations, we used the historical simulations of the CMIP5 from the years 1900 to 1999 (Taylor et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2012\u003c/span\u003e) and CMIP6 for the period 1850\u0026ndash;2014 (Eyring et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) models, utilizing every model with the variables required for this analysis. We use a single realization of each model and use all models with the required data available. See Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e for the 37 CMIP5 and 39 CMIP6 models taken into consideration.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003e2.2 ENSO Phase space\u003c/h2\u003e \u003cp\u003eFor the ENSO phase space analysis, we are following the approach by Dommenget and Al-Ansari (\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), Dommenget et al. (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), and Priya et al. (\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). We normalise the monthly means of \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e by dividing them by their standard deviations, resulting in a non-dimensional presentation of the variables, \u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e, respectively. The ENSO phase space is represented with \u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e on the x-axis and \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e on the y-axis. The cartesian coordinate system has been transformed into a spherical coordinate system with phase angle, \u0026#120593; = 0\u0026deg; in \u003cb\u003eh\u003c/b\u003e (y-direction) and 90\u0026deg; in \u003cb\u003eT\u003c/b\u003e (x-direction), where the phase angle, \u0026#120593;, rotates in a clockwise direction. The ENSO system anomaly, \u003cb\u003eS\u003c/b\u003e, is represented by a vector with two cartesian coordinates \u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e. When it is converted into the polar coordinate system, the magnitude of \u003cb\u003eS\u003c/b\u003e remains constant for constant radius and is independent of phase angle (\u0026#120593;). The ENSO system is further described by the magnitude \u003cb\u003eS\u003c/b\u003e and phase \u0026#120593;.\u003c/p\u003e \u003cp\u003eThe ENSO system's tendencies in the polar coordinate system are separated into radial and tangential components. The phase dependent growth rate is calculated by dividing the radial component of the tendencies by the magnitude of \u003cb\u003eS\u003c/b\u003e, which is defined as the tendency to move away from the origin (positive values) and towards the origin (negative values). The phase speed is calculated by dividing the tangential component of the tendencies by the magnitude of \u003cb\u003eS\u003c/b\u003e and is defined as the system's tendency to circle around the origin. Positive values indicate clockwise motion, while negative values indicate anti-clockwise motion; additionally, smaller positive values indicate a slower transition in the ENSO cycle. The average time required to complete one complete cycle is determined through the integration of the phase speed from 0\u0026deg; to 360\u0026deg;.\u003c/p\u003e \u003cp\u003eNote, that a proper analysis of the ENSO phase space statistics assumes that \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e are orthogonal to each other (Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Thus, they have a zero cross-correlation at time lag zero. If \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e are not orthogonal to each other (e.g., lag-zero cross-correlation is not zero), then the phase space statistics are skewed towards the diagonals between the \u003cb\u003eT\u003c/b\u003e-axis and \u003cb\u003eh\u003c/b\u003e-axis of the diagrams, which will result into skewed statistics of probabilities, growth rates, and phase speeds.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003e2.3 ReOsc Model\u003c/h2\u003e \u003cp\u003eThe linear ReOsc model is based on two tendency equations with six parameters as given below (Burgers et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2005\u003c/span\u003e):\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$$\\:\\frac{dT\\left(t\\right)}{dt}={a}_{11}T\\left(t\\right)+\\:{a}_{12}h\\left(t\\right)+{\\xi\\:}_{1}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equ2\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e\n$$\\:\\frac{dh\\left(t\\right)}{dt}=\\:{a}_{21}T\\left(t\\right)+\\:{a}_{22}h\\left(t\\right)+\\:{\\xi\\:}_{2}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere a\u003csub\u003e11\u003c/sub\u003e and a\u003csub\u003e22\u003c/sub\u003e are the damping factor or growth rate of \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e respectively, a\u003csub\u003e12\u003c/sub\u003e is the coupling factor for \u003cb\u003eT\u003c/b\u003e to \u003cb\u003eh\u003c/b\u003e, a\u003csub\u003e21\u003c/sub\u003e is the coupling factor for \u003cb\u003eh\u003c/b\u003e to \u003cb\u003eT\u003c/b\u003e, ξ\u003csub\u003e1\u003c/sub\u003e is the noise forcing of \u003cb\u003eT\u003c/b\u003e, and ξ\u003csub\u003e2\u003c/sub\u003e is the noise forcing of \u003cb\u003eh\u003c/b\u003e. The model parameters of the ReOsc model are determined by carrying out a multivariate linear regression on the monthly mean tendencies of \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e in relation to their monthly mean values of \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e (Burgers et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Jansen et al. 2009; Dommenget and Al-Ansari \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Priya et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn reference to the normalised presentation of the ENSO phase space (\u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e), we normalise all values of the ReOsc model. The ENSO phase space resulting from the linear ReOsc model fit can be estimated by integrating the Eqs.\u0026nbsp;(\u003cspan refid=\"Equ1\" class=\"InternalRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan refid=\"Equ2\" class=\"InternalRef\"\u003e2\u003c/span\u003e) with random white noise forcings ξ\u003csub\u003e1\u003c/sub\u003e and ξ\u003csub\u003e2\u003c/sub\u003e over a period of 10\u003csup\u003e4\u003c/sup\u003e.\u003c/p\u003e \u003c/div\u003e"},{"header":"3. CMIP simulations of ENSO dynamics","content":"\u003cp\u003eOur main results section is divided into four subsections. We start the analysis in subsection 3.1 with investigating the cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e, which is a fundamental statistic on which the ENSO phase space analysis is based on. This analysis will also focus on the regional shift in ENSO patterns, as it affects cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e. In the following subsection 3.2, we analyse the ENSO phase space in CMIP simulations, which is the main focus of this study. This is followed by the analysis of the linear ReOsc model fitted to the CMIP simulations. We conclude our discussion by examining the effect of wind stress on ENSO dynamics and how well CMIP models capture this atmospheric feedback in the subsection 3.4.\u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Regional variations of T and h\u003c/h2\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eFig. 1a shows the observed cross-correlation between \u003cstrong\u003e\u003cem\u003eT\u003c/em\u003e\u003c/strong\u003e and \u003cstrong\u003e\u003cem\u003eh\u003c/em\u003e\u003c/strong\u003e. It shows the characteristic out-of-phase relation between \u003cstrong\u003e\u003cem\u003eT\u003c/em\u003e\u003c/strong\u003e and \u003cstrong\u003e\u003cem\u003eh\u003c/em\u003e\u003c/strong\u003e with positive cross-correlation when \u003cstrong\u003e\u003cem\u003eh\u003c/em\u003e\u003c/strong\u003e leads \u003cstrong\u003e\u003cem\u003eT\u003c/em\u003e\u003c/strong\u003e, zero cross-correlation at lag zero, and strong negative correlations when \u003cstrong\u003e\u003cem\u003eT\u003c/em\u003e\u003c/strong\u003e leads \u003cstrong\u003e\u003cem\u003eh\u003c/em\u003e\u003c/strong\u003e.\u003c/p\u003e\u003cp\u003eThe CMIP 5 and 6 ensembles are very similar to each other, and both show significant deviations from the observed cross-correlation (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea). All CMIP 5 and 6 simulations have negative cross-correlation at lag zero, shifting the whole ensemble of simulations into clear negative correlations at lag zero. This shift also affects the out-of-phase cross-correlations with lower positive correlation when \u003cb\u003eh\u003c/b\u003e leads \u003cb\u003eT\u003c/b\u003e that tend to peak at longer (10mon.) lead times than observed (7mon.). In turn, the negative cross-correlations when \u003cb\u003eT\u003c/b\u003e leads \u003cb\u003eh\u003c/b\u003e are stronger in the CMIP models than in the observations and tend to have short lead times (5mon.) than observed (7mon.).\u003c/p\u003e \u003cp\u003eThe shift towards negative lag-zero cross-correlations between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e in CMIP simulations could indicate an eastward shift in the thermocline depth variability associated with ENSO in the CMIP simulations. This would be consistent with observed regional variations in the SST cross-correlation with thermocline depth variability along the equatorial Pacific, which is known to be more positive in the eastern Pacific and more negative in the western Pacific (Burgers and Stephenson \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e1999\u003c/span\u003e; Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eGiven that the central assumption in the ENSO phase space analysis is that the cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e is zero at lag zero, we now focus on the hypothesis that regional shifts in the simulated ENSO patterns could affect ENSO dynamics (e.g., the cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e) in CMIP simulations. We now analyse the equatorial patterns of the ENSO modes in the models, focussing on the SST pattern first and then on the more important regional shifts in the thermocline depth variability.\u003c/p\u003e \u003cp\u003eThe \u003cb\u003eT\u003c/b\u003e index in the ReOsc model is most closely related to the empirical orthogonal function (EOF) mode-1 of the tropical Pacific SST, see Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea. It is marked by a horseshoe pattern with warm SST anomalies in the eastern equatorial Pacific that transition into negative SST anomalies in the western equatorial and off-equatorial Pacific. This transition crosses zero equatorial SST anomalies at the longitude of about 159\u003csup\u003eo\u003c/sup\u003eE (see dashed line in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea).\u003c/p\u003e \u003cp\u003eIn Figs.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb-d, we show three scenarios for EOF mode-1 of SST from CMIP simulations to highlight cases closest to the observed (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb; ACCESS-CM2), average (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec; MRI-ESM2-0) and worst (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed; NorESM2-MM) in relation to the observed transition line for zero equatorial SST anomalies. While the best CMIP simulation is close to the observed with some negative SST anomalies in the western equatorial Pacific, the worst model has essentially only positive SST anomalies over the whole equatorial Pacific.\u003c/p\u003e \u003cp\u003eIn general, the EOF mode-1 of tropical Pacific SST of CMIP simulations is similar to the observed, but the transition line for zero equatorial SST anomalies is, in ensemble mean, about 20\u003csup\u003eo\u003c/sup\u003e further west for both CMIP 5 and 6 simulations (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea). This is consistent with the westward shift reported in previous studies (Weller and Cai \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Yang and Giese \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Capotondi et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Wang et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Guilyardi et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Jiang et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe observed EOF mode-1 of thermocline depth variability is a dipole pattern (tilting mode; not shown) that has an in-phase relation with the \u003cb\u003eT\u003c/b\u003e index in the ReOsc model and thus does not represent the \u003cb\u003eh\u003c/b\u003e index (Meinen and McPhaden \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2000\u003c/span\u003e). The EOF mode-2 of the equatorial Pacific is similar to a basin-wide mode, but the pattern is somewhat shifted to the eastern Pacific, see Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee. This mode is more closely correlated to \u003cb\u003eh\u003c/b\u003e (correlation\u0026thinsp;=\u0026thinsp;0.91). While the EOF mode-2 is mostly a monopole, it does have a transition from negative anomalies to positive anomalies in the far western equatorial Pacific at a longitude of about 146\u0026deg;E (see dashed line in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe EOF mode-2 of the CMCC-CM2-SR5 (best model) exhibits a similar spatial pattern, with near-zero values present around the same region with observation (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ef). A medium CMIP6 model, ACCESS-ESM1-5, demonstrates that the near-zero values shift towards the eastern equatorial Pacific region, and we can see a dipole-like structure with deeper thermocline depth in the western and shallower thermocline depth in the central to eastern pacific region (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eg). In addition, the CNRM-CM6-1 CMIP6 model (worst model) exhibits a clear dipole pattern in EOF mode-2, more similar to the tilting mode described in Meinen and McPhaden (\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2000\u003c/span\u003e; Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eh).\u003c/p\u003e \u003cp\u003eIn the ensemble of CMIP5 and 6 models, we can see that the EOF mode-2 patterns of thermocline depth (\u003cb\u003eh\u003c/b\u003e) suggest an eastward shift in the equatorial Pacific region of about 20\u0026deg;E (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn summary, we find significant regional shifts in the SST and thermocline depth patterns associated with the ReOsc model. Interestingly, we find that the SST pattern is shifted to the west and the thermocline depth pattern is shifted to the east. At this point, we don't have a clear understanding of the reasons behind these biases, but it appears very likely that such regional shifts can affect the cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e.\u003c/p\u003e \u003cp\u003eFollowing the analysis of Dommenget et al. (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) we know that an eastward shift of the \u003cb\u003eh\u003c/b\u003e index region leads to a more positive lag-zero cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e. This could compensate for the negative lag-zero cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e found in CMIP models (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea). To fit into the ReOsc model and to follow the phase space concept, CMIP models should have an out-of-phase relationship between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e (e.g., zero lag-zero cross-correlation).\u003c/p\u003e \u003cp\u003eTo define a revised \u003cb\u003eh\u003c/b\u003e index for the CMIP models which has zero lag-zero cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e, we shifted the index region for \u003cb\u003eh\u003c/b\u003e eastward along the equator (\u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e) until we found a region that most closely follows the observed zero lag-zero cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e, for all CMIP simulations. We concluded that the out-of-phase characteristics for CMIP models in the ensemble mean are very similar to observation when \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e is considered for the east equatorial pacific region (5\u0026deg;S to 5\u0026deg;N, 190\u0026deg;E to 80\u0026deg;W) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb). However, individual models do have some significant variations from this overall good fit.\u003c/p\u003e \u003cp\u003eIn addition to this, we also investigated various regions of \u003cb\u003eT\u003c/b\u003e to understand the impact of the SST pattern shift on the ENSO phase space for CMIP models. However, we found that changing the \u003cb\u003eT\u003c/b\u003e regions had no significant impact on the out-of-phase relationship between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh.\u003c/b\u003e As a result, we have focused solely on the regional shift of \u003cb\u003eh\u003c/b\u003e for CMIP models for the remaining analysis.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003e3.2 ENSO phase space analysis for CMIP models\u003c/h2\u003e \u003cp\u003eThe observed ENSO phase space is depicted in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea. The main features show a clockwise rotation of the tendencies, with clear phase propagation through all four phases of the ENSO cycle. Unlike the phase space of an idealised linear ReOsc model, which is symmetric in all phases (Dommenget and Al-Ansari \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), the observed variability is clearly skewed towards El Ni\u0026ntilde;o (positive \u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e values) to discharge states (negative \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e values) (quarter Q2) and the tendencies are stronger in quarter Q2 and weaker in quarter Q4 (La Ni\u0026ntilde;a to recharge phases). A more detailed discussion is presented in Dommenget et al. (\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eExamples of CMIP models are shown for the phase space estimated based on \u003cb\u003eh\u003c/b\u003e (left column in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e) and based on \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e (right column), see Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e. We selected a best (NorESM2-MM), medium (CESM2-WACCM), and worst (BCC-ESM1) model based on the correlation values between the mean phase space estimated based on \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e and the observed value. If we compare the left with the right column, we can notice that the phase spaces of the left column are pronounced along the diagonal of Q2 to Q4, which is a signature of a negative cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e (Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The right column, in turn, is more equally distributed on all four quarters of the phase space and generally, closer to the observed behaviour. This highlights that the CMIP models compare better with the observed phase space when \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e is considered.\u003c/p\u003e \u003cp\u003eThe best model (NorESM2-MM) closely reproduces the observed ENSO phase space with high correlation in the mean phase space (0.9), clear clockwise propagation in all four phases, and variability skewed towards stronger positive \u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e (El Ni\u0026ntilde;o), and stronger negative \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e (discharge phase). In turn, the worst model (BCC-ESM1) shows opposite asymmetries in the phase space to those observed, with a negative correlation in the mean phase space (-0.85), indicating that it has extreme recharge states instead of the observed extreme discharge states. However, it also shows clear clockwise propagation in all four phases.\u003c/p\u003e \u003cp\u003eThe mean statistical features in the ENSO phase space are shown for all CMIP models and the observations in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e. Again, the left column shows the CMIP models based on \u003cb\u003eh\u003c/b\u003e and the right column based on \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e. The observed values are shown in both columns for comparison (both based on \u003cb\u003eh\u003c/b\u003e). The CMIP models phase space based on \u003cb\u003eh\u003c/b\u003e (left column) shows clear biases in all statistics. The mean phase space and the probability distribution of \u003cb\u003eS\u003c/b\u003e\u0026thinsp;\u0026gt;\u0026thinsp;2 are both more pronounced along the diagonal of Q2 to Q4, the growth rate is more pronounced along the \u003cb\u003eT\u003c/b\u003e-axis and the phase speed along the diagonal of Q1 to Q3. All these biases directly result from the negative cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e (Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). As a result, the CMIP models phase space statistics based on \u003cb\u003eh\u003c/b\u003e do not correlate well with the observed values.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe phase space statistics of the CMIP models based on \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e show much less obvious biases and correlate much better with the observed values. The CMIP models generally do capture the asymmetry in the mean phase space, with the largest mean value in Q2 (El Ni\u0026ntilde;o and discharge state), smallest value in Q4 (La Ni\u0026ntilde;a and recharge state), and a correlation of the ensemble mean value in CMIP 5 and 6 of 0.88 and 0.81, respectively. However, the CMIP ensemble does not quite capture the intensity of the shift from Q4 to Q2, indicating the CMIP models are underestimating the non-linearity in the ENSO phase space.\u003c/p\u003e \u003cp\u003eThis also holds for the probability distribution of \u003cb\u003eS\u003c/b\u003e\u0026thinsp;\u0026gt;\u0026thinsp;2 (extreme values), which in observations show a clear shift to Q2 and away from Q4. This is captured by the CMIP models, but again not with the right intensity. In particular, the absence of observed extremes in Q4 is not fully captured by the CMIP models. Furthermore, we see a very large spread in CMIP ensemble members, indicating widely different behaviour between the models.\u003c/p\u003e \u003cp\u003eThe observed mean growth rate as a function of the ENSO phase shows positive values mostly around the recharge towards the El Ni\u0026ntilde;o state and negative values in the transition from discharge to a La Ni\u0026ntilde;a state (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eg). This signature is somewhat captured in the CMIP ensemble, but only with a moderate correlation. The CMIP models capture the negative growth rate in the transition from discharge to a La Ni\u0026ntilde;a state better than the positive values during the recharge to the El Ni\u0026ntilde;o state. Again, the CMIP models show a wide spread between the ensemble members. Notably, we find several models with very strong negative growth rates around the discharge and recharge phases, suggesting that these ENSO states can collapse much faster in these models than observed.\u003c/p\u003e \u003cp\u003eThe observed phase speed is fastest, or most clear, in Q2 (after El Ni\u0026ntilde;o states) and slowest, or least clear in Q4 (after La Ni\u0026ntilde;a states; Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eh). The CMIP model ensemble can partly capture this signature, but only with a moderate positive correlation. In particular, the models have problems in capturing the smaller phase speed after La Ni\u0026ntilde;a events. This suggests that models tend to have a faster or clearer transition from La Ni\u0026ntilde;a states to the recharge state. The model ensemble also has a wide range of different phase speeds. Some models transition through some ENSO phases more than twice as fast than observed, and some models have phase speeds near zero in some ENSO phases, suggesting the models have no clear ENSO cycle in these phases.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe phase speed analysis suggests that CMIP models are oscillating a bit faster than observed. Figure\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e shows the observed and CMIP simulated power spectrum of the \u003cb\u003eT\u003c/b\u003e index. We can first notice that the ensemble mean power spectrum of CMIP 5 and 6 are in very good agreement with the observed, with no clear shift in the power to higher frequencies. However, both CMIP ensembles have a small tendency to peak at slightly shorter periods (~\u0026thinsp;3yrs period) than observed (~\u0026thinsp;4yrs period). Individual members of the ensemble can, however, behave quite differently, with some models having clear peaks at different periods and other models having less obvious peaks with a more continuous power spectrum.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Fitted linear ReOsc model for CMIP\u003c/h2\u003e \u003cp\u003eThe ENSO phase space is based on the ReOsc model. It therefore does help to analyse the ReOsc model fitted to the data of the CMIP models to understand what is causing the characteristics of the ENSO phase space. Here we focus on the linear model, thus not considering non-linearities, and without considering any seasonality in the fitted parameters. Neglecting non-linearities is a strong limitation, as will be shown below, and is only reflecting that we yet do not know what non-linear model can describe the observed asymmetries in the ENSO phase space (Dommenget and Al-Ansari \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eFigure 7 shows the fitted ReOsc model parameters for all CMIP models and observations for models based on \u003cstrong\u003e\u003cem\u003eh\u003c/em\u003e\u003c/strong\u003e (left column) and based on\u0026nbsp;\u003cstrong\u003e\u003cem\u003eh\u003csub\u003eshift\u003c/sub\u003e\u0026nbsp;\u003c/em\u003e\u003c/strong\u003e(right column). The observed estimate is the same in both columns. Here all parameters are normalized (as in \u003cstrong\u003e\u003cem\u003eT\u003csub\u003en\u003c/sub\u003e\u003c/em\u003e\u003c/strong\u003e and \u003cstrong\u003e\u003cem\u003eh\u003csub\u003en\u003c/sub\u003e\u003c/em\u003e\u003c/strong\u003e) to allow a better discussion of the dynamical implication for the ENSO phase space.\u003c/p\u003e \u003cp\u003eWe can first-of-all note that there is quite some spread between models, with many models having parameters being far away from the observed values in both estimates. Further, we find that the two estimates of the ReOsc models are quite different from each other. The models based on \u003cb\u003eh\u003c/b\u003e are strongly biased in the growth rates and the noise strength in \u003cb\u003eh\u003c/b\u003e relative to the observed values but are much less biased in the coupling parameters. In turn, the models based on \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e are more strongly biased in the coupling parameters.\u003c/p\u003e \u003cp\u003eThe results for CMIP5 models are also somewhat different from Vijayeta and Dommenget (\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; hereafter VD18), which used \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003e20\u003c/b\u003e\u003c/sub\u003e to estimate \u003cb\u003eh\u003c/b\u003e instead of \u003cb\u003eZ\u003c/b\u003e\u003csub\u003e\u003cb\u003emxg\u003c/b\u003e\u003c/sub\u003e. This does alter the ReOsc model parameter fits. However, the CMIP5 models tend to overestimate the damping (negative a\u003csub\u003e22\u003c/sub\u003e) of \u003cb\u003eh\u003c/b\u003e in both estimates present here and in VD18.\u003c/p\u003e \u003cp\u003eWe now focus on the ReOsc model fits based on \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e as they better represent the observed ENSO phase space. Here we are primarily interested in what dynamical elements of the ReOsc model are causing the simulated ENSO phase space characteristics. Any deviation of the ENSO phase space from a perfect cycle (a perfectly symmetric phase space), can be linked to asymmetries (deviations for the dashed lines in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e) in the ReOsc models parameter pairs: growth rates (a\u003csub\u003e11\u003c/sub\u003e, a\u003csub\u003e22\u003c/sub\u003e), coupling (a\u003csub\u003e12\u003c/sub\u003e, a\u003csub\u003e21\u003c/sub\u003e) or noise forcings (ξ\u003csub\u003e1\u003c/sub\u003e, ξ\u003csub\u003e2\u003c/sub\u003e). In a non-linear model, they can also result from non-linearities in any of these aspects.\u003c/p\u003e \u003cp\u003eIndividual models do show strong deviations from the dashed lines in the growth rates and coupling parameters, but less so in the noise forcings. In the ensemble mean the largest asymmetries are found in the coupling parameters, where nearly all models are below the dashed line, indicating a stronger coupling of \u003cb\u003eT\u003c/b\u003e to \u003cb\u003eh\u003c/b\u003e (a\u003csub\u003e12\u003c/sub\u003e) than the coupling of \u003cb\u003eh\u003c/b\u003e to \u003cb\u003eT\u003c/b\u003e (a\u003csub\u003e21\u003c/sub\u003e), which is quite different from the observed symmetric coupling. The coupling parameters are most important for the period of ENSO (Wyrtki \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e1985\u003c/span\u003e). The Wyrtki-index, which measures the peak period of ENSO as a function of a\u003csub\u003e12\u003c/sub\u003e and a\u003csub\u003e21\u003c/sub\u003e, is for most models slightly shifted to shorter periods (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003ee), consistent with somewhat larger phase speed values seen in the phase space diagrams (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eh).\u003c/p\u003e \u003cp\u003eIntegrating the linear ReOsc model fitted to each CMIP model allows us to compute the phase space statistics based on the resulting \u003cb\u003eT\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003en\u003c/b\u003e\u003c/sub\u003e values, see left column in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e. The phase spaces resulting from the ReOsc model are in general similar to those of the CMIP models, but, due to the linear approach, the phase spaces are by construction symmetric (e.g., all structures are mirroring at the origin) and therefore lack all asymmetric variations. This is particularly important when compared against the observed phase space, because most of the interesting observed phase space characteristics are asymmetric. None of the observed asymmetries can be captured by the analysis of the linear ReOsc model, indicating that such structures must result from the non-linear processes. Thus, a first and important outcome of the linear ReOsc model fit is that the CMIP model mismatch to the observed phase space characteristics are likely to result from non-linear processes.\u003c/p\u003e \u003cp\u003eDespite the limitation of the linear ReOsc model fit, we can still gain some understanding about the CMIP model ensemble spread and how different linear aspects of the ReOsc model affect the ENSO phase space statistics. The mean \u003cb\u003eS\u003c/b\u003e, the probability of extreme \u003cb\u003eS\u003c/b\u003e\u0026thinsp;\u0026gt;\u0026thinsp;2 and the phase speed are all only varying along the diagonals of the phase space but have no variation along the \u003cb\u003eT\u003c/b\u003e-axis and \u003cb\u003eh\u003c/b\u003e-axis of the diagrams (left column in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). The growth rates are only varying along \u003cb\u003eT\u003c/b\u003e-axis and \u003cb\u003eh\u003c/b\u003e-axis, but not along the diagonals. However, the actual data of the CMIP ensembles does show more complex variations in all statistical aspects (right column Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e). This in turn suggests that variations not captured by the linear model parameters are a result of non-linear processes.\u003c/p\u003e \u003cp\u003eWe can further test the sensitivity of the ENSO phase space to the individual ReOsc model parameter variations in the CMIP ensembles, by varying only a subset of the parameters, while holding the other parameters at the ensemble mean values following the approach of VD18. This approach can allow us to determine which ReOsc model dynamics are causing asymmetries or variations in the ENSO phase space. The following three scenarios are presented: asymmetries in the growth rates (second column of Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e), coupling (third column of Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e), and forcing strength (right column of Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe spread in the mean \u003cb\u003eS\u003c/b\u003e values in the phase space of the linear ReOsc model fits (first row Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e) is much smaller in each of the three sensitivity experiments, suggesting that the spread in mean \u003cb\u003eS\u003c/b\u003e is resulting from a more complex combination of several parameters, but does not result from any of the three asymmetries tested. The ensemble mean \u003cb\u003eS\u003c/b\u003e value is mostly affected by asymmetries in the coupling and forcings, which both have opposing effects. Here the asymmetries in linear coupling do have some similarities with the observed, but they lack the non-linear elements. Asymmetries in the growth rates have little impact on the mean \u003cb\u003eS\u003c/b\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe phase space variations of the extreme values (\u003cb\u003eS\u003c/b\u003e\u0026thinsp;\u0026gt;\u0026thinsp;2) are affected by all their parameter asymmetries, but most strongly by the asymmetries in the coupling and forcings (second row Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). The phase space variations of the growth rates are clearly linked to the asymmetries in the growth rates, but asymmetries in the coupling and forcings do also have a significant effect (third row Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe sensitivities of the phase speed to the linear ReOsc parameters show complex behaviour (last row Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). When we only consider the asymmetries in the growth rates or forcing strength, we get much larger phase speeds than observed in the CMIP ensemble. This suggests that the asymmetry in the coupling parameters has a strong impact on the phase speed. We can further notice that none of the three asymmetries reflect the overall variations in the phase speed, suggesting that more complex combinations of different ReOsc model parameters affect the phase speed.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003e3.4 Wind dynamics for CMIP models\u003c/h2\u003e \u003cp\u003eThe surface wind response to SST anomalies is a key element of the Bjerknes feedback and is one of the main processes that control ENSO growth rate and period. While this relation does not directly relate to the ENSO phase space or the ReOsc model, it is implicitly related to both as it affects all parameters of the ReOsc model, in particular the growth rate of \u003cb\u003eT\u003c/b\u003e.\u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003ea shows the cross-correlation between the NINO4 region (5\u0026deg;S to 5\u0026deg;N, 160\u0026deg;E to150\u0026deg;W) zonal wind stress (\u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e) and \u003cb\u003eT\u003c/b\u003e for the observations and the CMIP models. The observed cross-correlation shows a strong and mostly in-phase relation between \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eT\u003c/b\u003e, consistent with a strong response of \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e to variations in \u003cb\u003eT\u003c/b\u003e. The relation is, however, strongest when \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e leads \u003cb\u003eT\u003c/b\u003e by about one month, indicating also that variations in \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e are forcing variations in \u003cb\u003eT\u003c/b\u003e.\u003c/p\u003e \u003cp\u003eOverall, we see a close agreement between CMIP 5 and 6, and between the models and the observations. Despite the good agreement, there are some deviations in the CMIP models from observations. The CMIP models have a slightly low cross-correlation at lag zero, indicating a weaker response of \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e to \u003cb\u003eT\u003c/b\u003e. The maximum cross-correlation in the CMIP ensembles is when \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e leads \u003cb\u003eT\u003c/b\u003e by about two months, indicating a stronger forcing of \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e on variations in \u003cb\u003eT\u003c/b\u003e than observed. The spread in these relations is quite large in the CMIP ensembles, indicating quite diverse behaviour in the individual models.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003eb shows the regression value of \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e per \u003cb\u003eT\u003c/b\u003e, which is often referred to as the atmospheric Bjerknes feedback (Bjerknes \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e1969\u003c/span\u003e). The observed atmospheric Bjerknes feedback is approximately 12.5*10\u003csup\u003e\u0026minus;\u0026thinsp;3\u003c/sup\u003e Nm\u003csup\u003e\u0026minus;\u0026thinsp;2\u003c/sup\u003e/\u003csup\u003eo\u003c/sup\u003eC, whereas the CMIP ensemble averages are 7.4 (CMIP5) and 7.8*10\u003csup\u003e\u0026minus;\u0026thinsp;3\u003c/sup\u003e (CMIP6) Nm\u003csup\u003e\u0026minus;\u0026thinsp;2\u003c/sup\u003e/\u003csup\u003eo\u003c/sup\u003eC. Thus, CMIP models underestimate this important feedback by approximately 40%, with none of the CMIP models reaching the observed value and no significant improvement from CMIP 5 to 6. There is a substantial spread within the CMIP model ensemble that is also of similar magnitude in both CMIP 5 and 6. These results are largely consistent with previous studies (Guilyardi et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Bellenger et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Vijayeta and Dommenget \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Guilyardi et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eTo further examine the variability in \u003cb\u003eT\u003c/b\u003e, \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eh\u003c/b\u003e we looked at the standard deviation of these variables (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e10\u003c/span\u003e). The standard deviation of \u003cb\u003eT\u003c/b\u003e in the CMIP model averages closely resembles observations, with the CMIP5 model average being smaller than the CMIP6 value (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e10\u003c/span\u003ea). However, there are significant differences between individual CMIP models, with some models having significantly smaller and others larger standard deviations than the ensemble model or observed.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe observed \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e has a significantly larger standard deviation than both CMIP5 and CMIP6 model averages (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e10\u003c/span\u003eb). This would be consistent with the weaker \u003cb\u003eT\u003c/b\u003e forcing on \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e. Some CMIP models show standard deviations that are closer to the observation, while others have noticeably smaller standard deviations. The results suggest that the weaker \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e variations may be a cause or a reflection of the weaker \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e relation to \u003cb\u003eT\u003c/b\u003e in the models.\u003c/p\u003e \u003cp\u003eThe CMIP6 ensemble mean standard deviation for \u003cb\u003eh\u003c/b\u003e is very close to the observed, whereas the CMIP5 ensemble mean is a bit smaller, and some show large deviations for the ensemble mean (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e10\u003c/span\u003ec). However, more relevant for the ENSO phase space discussion is the variability of \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e10\u003c/span\u003ed). This is significantly weaker than the observed \u003cb\u003eh\u003c/b\u003e for both CMIP ensembles. Given a similar cross-correlation between \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eT\u003c/b\u003e in the CMIP model as in the observations (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb) and the smaller standard deviation of \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e in the CMIP as shown above, we would expect larger values of a\u003csub\u003e12\u003c/sub\u003e in CMIP models than observed (see Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003ee), since the regression values are divided the smaller standard deviations of \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e. Thus, per \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e anomaly in the CMIP models we get a larger \u003cb\u003eT\u003c/b\u003e anomaly than in observations. It should be noted here that this does not need necessarily imply a larger sensitivity of \u003cb\u003eT\u003c/b\u003e to \u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e, because correlation does not allow to determine causality here.\u003c/p\u003e \u003c/div\u003e"},{"header":"4. Summary and discussions","content":"\u003cp\u003eIn this study, we examined how well the CMIP5 and CMIP6 models can simulate the observed ENSO phase space characteristics. In the context of the ReOsc model, the ENSO phase space presents ENSO as a cycle for which we can analyze important statistics, such as extreme value probabilities, growth rates, and phase speeds, as a function of the different phases of ENSO. This allows us to consider linear and non-linear aspects of the amplitude, growth rates, and phase transition speeds at the same time.\u003c/p\u003e \u003cp\u003eA key starting point for the analysis of the ENSO phase space is the cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e, which is assumed and observed to be out-of-phase (orthogonal) to allow for a proper evaluation of the phase space statistics, and it is an assumed relation for an idealized linear ReOsc model. We found that this relationship is clearly not out-of-phase for the CMIP models but is shifted to negative lag-zero cross-correlations. This shift in the cross-relation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e is linked to a regional shift in the ENSO patterns in CMIP models, with thermocline depth variability being shifted to the eastern equatorial Pacific.\u003c/p\u003e \u003cp\u003eIf the eastward shift in the thermocline depth variability (\u003cb\u003eh\u003c/b\u003e\u003csub\u003e\u003cb\u003eshift\u003c/b\u003e\u003c/sub\u003e) is considered in estimating the ENSO phase space, then we find a much better agreement between models and observations. This suggests that CMIP models do simulate the ENSO cycle with some regional displacements in the ocean-atmosphere interactions.\u003c/p\u003e \u003cp\u003eThe main characteristics of the observed phase space are mostly well simulated by the CMIP ensembles. The key feature of the observed higher likelihoods of extreme El Ni\u0026ntilde;o to discharge states and reduced likelihoods of extreme La Ni\u0026ntilde;a to recharge states are captured by the CMIP ensembles mean with a high correlation. However, the models strongly underestimate the magnitude of non-linearity in this aspect. This suggests that although models, in their ensemble mean, capture the non-linearity of phases of ENSO, the non-linear processes responsible for these characteristics are not strong enough (weaker non-linear processes in CMIP models than observation as suggested by An et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Sun et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) or are too linear.\u003c/p\u003e \u003cp\u003eThe observed non-linearity in growth rates of the ENSO phase space are less well captured by the CMIP ensembles. While they do capture the enhanced growth rates during the El Ni\u0026ntilde;o states, they lack increased growth rates during the recharge state, and they lack strong decay (negative growth rates) of ENSO amplitudes in the transition from discharge to La Ni\u0026ntilde;a states.\u003c/p\u003e \u003cp\u003eThe observed variations in the ENSO phase speed are also captured by the CMIP ensembles, but only with low correlations. The models can capture the phase speed in the transition from El Ni\u0026ntilde;o to the discharge state, but the transition is too fast (or too clear) from the La Ni\u0026ntilde;a to the recharge state, where the observed transition is much slower (or less clear). Combined with the lack of variations in the growth rates this suggests that the ENSO cycle in CMIP models is too regular, with too little differences between the El Ni\u0026ntilde;o and La Ni\u0026ntilde;a phases.\u003c/p\u003e \u003cp\u003eThe ENSO phase space is strongly linked to the idea of the ReOsc model. It therefore does make sense to explain the observed and simulated characteristics of the phase space by the ReOsc model dynamics. We therefore fitted the linear ReOsc model parameters to the observed and simulated \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e data but encountered a number of problems. First, the regional shift in the ENSO pattern in CMIP models does strongly affect the ReOsc model estimates. Second, the most interesting aspects of the observed ENSO phase space must result from non-linear dynamics, where it is unclear what kind of dynamics these are (Dommenget et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eWe further analyzed the relation between SST (\u003cb\u003eT\u003c/b\u003e), wind stress (\u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e) and thermocline depth (\u003cb\u003eh\u003c/b\u003e) variability. Like previous studies, we find that the linear relation of \u003cb\u003eτ\u003c/b\u003e\u003csub\u003e\u003cb\u003ex\u003c/b\u003e\u003c/sub\u003e and \u003cb\u003eT\u003c/b\u003e is much weaker than observed. This weaker relation is likely to contribute to the weaker non-linear ENSO phase space characteristics, since previous studies have indicated that the wind-SST relation is a key element of the non-linear behavior of ENSO. The variability of \u003cb\u003eh\u003c/b\u003e in the eastward shifted region of the Pacific, which is most important for CMIP models, is also weaker than observed. This, combined with the weaker wind-SST relation indicates that the ENSO mode in the CMIP model is more wind-driven and less influenced by the thermocline variability. A previous study by Zhao and Sun 2022 found that the CMIP models have a shallower thermocline depth in the eastern equatorial Pacific than observed, which could explain the weaker response of the above-mentioned processes. However, further analysis into the SST, wind stress and thermocline depth interactions are needed to understand how they impact the ENSO phase space. It is worth noting that Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e made no mention of thermocline depth biases in their ENSO metrics package, which is quite surprising, given that thermocline depth and their biases have a significant impact on the ENSO phase space characteristics and its dynamics as we found in our study.\u003c/p\u003e \u003cp\u003eWhile the ensemble mean of the CMIP models does show good agreement with the observations in many aspects, there is a widespread within the model ensembles. Some models show opposite non-linear ENSO phase spaces from those observed, marking severe biases in the model behavior. Other models clearly outperform the ensemble mean, indicating that much more realistic ENSO simulations are possible within the CMIP model framework.\u003c/p\u003e \u003cp\u003eIn the comparison of the CMIP5 versus the CMIP6 ensembles, we find no substantial differences in the ENSO phase space dynamics, indicating no substantial improvements from CMIP5 to CMIP6. In particular, we find no differences in the skill of the model in simulating the significant non-linear ENSO behavior. This is a bit disappointing and worrying, as it indicates that no progress has been made from CMIP5 to CMIP6. This finding is largely consistent with Planton et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2021\u003c/span\u003e.\u003c/p\u003e \u003cp\u003eMore future work is required to better understand the limitations of the CMIP models in simulating the non-linear ENSO phase space characteristics. A better understanding of how the ReOsc model can incorporate the non-linear aspects to simulate the observed ENSO phase space would be important to better evaluate CMIP models. The regional shift in the simulation of ENSO dynamics are also an important bias. We still don't know why these biases exist and what are the reason for these discrepancies, necessitating further analysis and a deeper understanding of this question, as well as opportunities for future studies to address these concerns.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e \u003ch2\u003eCompeting Interests\u003c/h2\u003e \u003cp\u003eNone\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eThis research is funded by the Australian Research Council (ARC), discovery project \u0026ldquo;Improving projections of regional sea level rise and their credibility\u0026rdquo; (DP200102329) and the Centre of Excellence for Climate Extremes (CLEX) Grant Number: CE170100023.\u003c/p\u003e\u003ch2\u003eAuthors Contributions\u003c/h2\u003e \u003cp\u003ePriyamvada Priya formulated the concept and designed this study. All authors contributed to the material preparation, collecting data, and analysis of this study. Priyamvada Priya wrote the first draft of the manuscript, and all authors commented and provided feedback on that. All authors read and approved the final version of the manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgements\u003c/h2\u003e \u003cp\u003eThe authors would like to thank Monash University Melbourne, Australia, for providing all the resources and instruments required to conduct this research. This research is funded by the Australian Research Council (ARC) Centre of Excellence for Climate Extremes (CLEX), CLEX Grant Number: CE170100023, and the ARC discovery project \u0026ldquo;Improving projections of regional sea level rise and their credibility\u0026rdquo; (DP200102329). The authors would also like to thank Australia's computing facility National Computational Infrastructure (NCI) and CLEX for providing us with a wide range of facilities and computational tools for this study.\u003c/p\u003e\u003ch2\u003eData Availability\u003c/h2\u003e \u003cp\u003eThe data sources for this study are cited in the text.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eAlizadeh O (2022) A review of the El Ni\u0026ntilde;o-Southern Oscillation in future. Earth Sci Rev 104246. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.earscirev.2022.104246\u003c/span\u003e\u003cspan address=\"10.1016/j.earscirev.2022.104246\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAn SI, Jin FF (2004) Nonlinearity and asymmetry of ENSO. J Clim 17(12):2399\u0026ndash;2412. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/1520-0442(2004)017\u0026lt;2399:NAAOE\u0026gt;2.0.CO;2\u003c/span\u003e\u003cspan address=\"10.1175/1520-0442(2004)017%3C2399:NAAOE%3E2.0.CO;2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAn SI, Ham YG, Kug JS, Jin FF, Kang IS (2005) El Ni\u0026ntilde;o\u0026ndash;La Ni\u0026ntilde;a asymmetry in the coupled model intercomparison project simulations. J Clim 18(14):2617\u0026ndash;2627. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/JCLI3433.1\u003c/span\u003e\u003cspan address=\"10.1175/JCLI3433.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAn SI, Tziperman E, Okumura YM, Li T (2020) : ENSO irregularity and asymmetry. \u003cem\u003eEl Ni\u0026ntilde;o Southern Oscillation in a changing climate\u003c/em\u003e, 153\u0026ndash;172, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/9781119548164.ch7\u003c/span\u003e\u003cspan address=\"10.1002/9781119548164.ch7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBayr T, Domeisen DI, Wengel C (2019) The effect of the equatorial Pacific cold SST bias on simulated ENSO teleconnections to the North Pacific and California. Clim Dyn 53:3771\u0026ndash;3789. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-019-04746-9\u003c/span\u003e\u003cspan address=\"10.1007/s00382-019-04746-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBayr T, Latif M (2023) ENSO atmospheric feedbacks under global warming and their relation to mean-state changes. Clim Dyn 60(9\u0026ndash;10):2613\u0026ndash;2631. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-022-06454-3\u003c/span\u003e\u003cspan address=\"10.1007/s00382-022-06454-3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBellenger H, Guilyardi \u0026Eacute;, Leloup J, Lengaigne M, Vialard J (2014) ENSO representation in climate models: From CMIP3 to CMIP5. Clim Dyn 42(7):1999\u0026ndash;2018. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-013-1783-z\u003c/span\u003e\u003cspan address=\"10.1007/s00382-013-1783-z\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBeobide-Arsuaga G, Bayr T, Reintges A, Latif M (2021) Uncertainty of ENSO-amplitude projections in CMIP5 and CMIP6 models. Clim Dyn 56:3875\u0026ndash;3888. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-021-05673-4\u003c/span\u003e\u003cspan address=\"10.1007/s00382-021-05673-4\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBjerknes J (1969) : Atmospheric teleconnections from the equatorial Pacific. \u003cem\u003eMonthly weather review\u003c/em\u003e, \u003cem\u003e97\u003c/em\u003e(3), 163\u0026ndash;172, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/1520-0493\u003c/span\u003e\u003cspan address=\"10.1175/1520-0493\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e(1969)097%3C0163:ATFTEP%3E2.3.CO;2\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBurgers G, Stephenson DB (1999) The normality of el ni\u0026ntilde;o. Geophys Res Lett 26(8):1027\u0026ndash;1030. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/1999GL900161\u003c/span\u003e\u003cspan address=\"10.1029/1999GL900161\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBurgers G, Jin FF, Van Oldenborgh GJ (2005) The simplest ENSO recharge oscillator. Geophys Res Lett 32(13). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/2005GL022951\u003c/span\u003e\u003cspan address=\"10.1029/2005GL022951\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCapotondi A, Wittenberg AT, Newman M, Di Lorenzo E, Yu JY, Braconnot P, \u0026hellip; and, Yeh SW (2015) Understanding ENSO diversity. Bull Am Meteorol Soc 96(6):921\u0026ndash;938. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/BAMS-D-13-00117.1\u003c/span\u003e\u003cspan address=\"10.1175/BAMS-D-13-00117.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCarton JA, Giese BS (2008) A reanalysis of ocean climate using Simple Ocean Data Assimilation (SODA). Mon Weather Rev 136(8):2999\u0026ndash;3017. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/2007MWR1978.1\u003c/span\u003e\u003cspan address=\"10.1175/2007MWR1978.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen M, Li T (2021) ENSO evolution asymmetry: EP versus CP El Ni\u0026ntilde;o. Clim Dyn 56(11):3569\u0026ndash;3579. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-021-05654-7\u003c/span\u003e\u003cspan address=\"10.1007/s00382-021-05654-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCrespo LR, Rodr\u0026iacute;guez-Fonseca MB, Polo I, Keenlyside N, Dommenget D (2022) Multidecadal variability of ENSO in a recharge oscillator framework. Environ Res Lett 17(7). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ep.074008, 10.1088/1748\u0026ndash;9326/ac72a3\u003c/span\u003e\u003cspan address=\"p.074008, 10.1088/1748\u0026ndash;9326/ac72a3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDee DP, Uppala S (2009) Variational bias correction of satellite radiance data in the ERA-Interim reanalysis. Q J Royal Meteorological Society: J atmospheric Sci Appl Meteorol Phys Oceanogr 135(644):1830\u0026ndash;1841. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/qj.493\u003c/span\u003e\u003cspan address=\"10.1002/qj.493\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDommenget D, Bayr T, Frauen C (2013) Analysis of the non-linearity in the pattern and time evolution of El Ni\u0026ntilde;o southern oscillation. Clim Dyn 40:2825\u0026ndash;2847. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-012-1475-0\u003c/span\u003e\u003cspan address=\"10.1007/s00382-012-1475-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDommenget D, Vijayeta A (2019) Simulated future changes in ENSO dynamics in the framework of the linear recharge oscillator model. Clim Dyn 53(7):4233\u0026ndash;4248. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-019-04780-7\u003c/span\u003e\u003cspan address=\"10.1007/s00382-019-04780-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDommenget D, Al-Ansari M (2022) Asymmetries in the ENSO phase space. Clim Dyn 1\u0026ndash;20. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-022-06392-0\u003c/span\u003e\u003cspan address=\"10.1007/s00382-022-06392-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDommenget D, Priya P, Vijayeta A (2023) ENSO phase space dynamics with an improved estimate of the thermocline depth. Clim Dyn. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-023-06883-8\u003c/span\u003e\u003cspan address=\"10.1007/s00382-023-06883-8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEyring V, Bony S, Meehl GA, Senior CA, Stevens B, Stouffer RJ, Taylor KE (2016) Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci Model Dev 9(5):1937\u0026ndash;1958. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.5194/gmd-9-1937-2016\u003c/span\u003e\u003cspan address=\"10.5194/gmd-9-1937-2016\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFrauen C, Dommenget D (2010) El Ni\u0026ntilde;o and La Ni\u0026ntilde;a amplitude asymmetry caused by atmospheric feedbacks. Geophys Res Lett 37(18). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/2010GL044444\u003c/span\u003e\u003cspan address=\"10.1029/2010GL044444\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGood SA, Martin MJ, Rayner NA (2013) EN4: Quality controlled ocean temperature and salinity profiles and monthly objective analyses with uncertainty estimates. J Geophys Research: Oceans 118(12):6704\u0026ndash;6716. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1002/2013JC009067\u003c/span\u003e\u003cspan address=\"10.1002/2013JC009067\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuilyardi E, Bellenger H, Collins M, Ferrett S, Cai W, Wittenberg A (2012) A first look at ENSO in CMIP5. Clivar Exchanges 17(58):29\u0026ndash;32\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuilyardi E, Capotondi A, Lengaigne M, Thual S, Wittenberg AT (2020) : ENSO modeling: History, progress, and challenges. \u003cem\u003eEl Ni\u0026ntilde;o Southern Oscillation in a changing climate\u003c/em\u003e, 199\u0026ndash;226, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/9781119548164.ch9\u003c/span\u003e\u003cspan address=\"10.1002/9781119548164.ch9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHersbach H, Bell B, Berrisford P, Hirahara S, Hor\u0026aacute;nyi A, Mu\u0026ntilde;oz-Sabater J, \u0026hellip; and, Th\u0026eacute;paut JN (2020) The ERA5 global reanalysis. Q J R Meteorol Soc 146(730):1999\u0026ndash;2049. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/qj.3803\u003c/span\u003e\u003cspan address=\"10.1002/qj.3803\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJiang W, Huang P, Huang G, Ying J (2021) Origins of the excessive westward extension of ENSO SST simulated in CMIP5 and CMIP6 models. J Clim 34(8):2839\u0026ndash;2851. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/JCLI-D-20-0551.1\u003c/span\u003e\u003cspan address=\"10.1175/JCLI-D-20-0551.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJin FF (1997) : An equatorial ocean recharge paradigm for ENSO. Part I: Conceptual model. \u003cem\u003eJournal of the atmospheric sciences\u003c/em\u003e, \u003cem\u003e54\u003c/em\u003e(7), 811\u0026ndash;829, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/1520-0469\u003c/span\u003e\u003cspan address=\"10.1175/1520-0469\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e(1997)054%3C0811:AEORPF%3E2.0.CO;2\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKessler WS (2002) Is ENSO a cycle or a series of events? Geophys Res Lett 29(23) 40\u0026thinsp;\u0026ndash;\u0026thinsp;1. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/2002GL015924\u003c/span\u003e\u003cspan address=\"10.1029/2002GL015924\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKim ST, Yu JY (2012) The two types of ENSO in CMIP5 models. Geophys Res Lett 39(11). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/2012GL052006\u003c/span\u003e\u003cspan address=\"10.1029/2012GL052006\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKim ST, Cai W, Jin FF, Yu JY (2014) ENSO stability in coupled climate models and its association with mean state. Clim Dyn 42:3313\u0026ndash;3321. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-013-1833-6\u003c/span\u003e\u003cspan address=\"10.1007/s00382-013-1833-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLatif M, Sperber K, Arblaster J, Braconnot P, Chen D, Colman A, \u0026hellip; and, Zebiak S (2001) ENSIP: the El Nino simulation intercomparison project. Clim Dyn 18:255\u0026ndash;276. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s003820100174\u003c/span\u003e\u003cspan address=\"10.1007/s003820100174\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLatif M, Keenlyside NS (2009) : El Ni\u0026ntilde;o/Southern Oscillation response to global warming. \u003cem\u003eProceedings of the National Academy of Sciences\u003c/em\u003e, \u003cem\u003e106\u003c/em\u003e(49), 20578\u0026ndash;20583, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1073/pnas.0710860105\u003c/span\u003e\u003cspan address=\"10.1073/pnas.0710860105\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcKenna S, Santoso A, Gupta AS, Taschetto AS, Cai W (2020) Indian Ocean Dipole in CMIP5 and CMIP6: characteristics, biases, and links to ENSO. Sci Rep 10(1):11500. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41598-020-68268-9\u003c/span\u003e\u003cspan address=\"10.1038/s41598-020-68268-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcPhaden MJ, Zebiak SE, Glantz MH (2006) ENSO as an integrating concept in earth science. Science 314(5806):1740\u0026ndash;1745. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1126/science.1132588\u003c/span\u003e\u003cspan address=\"10.1126/science.1132588\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMeinen CS, McPhaden MJ (2000) Observations of warm water volume changes in the equatorial Pacific and their relationship to El Ni\u0026ntilde;o and La Ni\u0026ntilde;a. J Clim 13(20):3551\u0026ndash;3559. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/1520-0442(2000)013\u0026lt;3551:OOWWVC\u0026gt;2.0.CO;2\u003c/span\u003e\u003cspan address=\"10.1175/1520-0442(2000)013%3C3551:OOWWVC%3E2.0.CO;2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOkumura YM (2019) ENSO diversity from an atmospheric perspective. Curr Clim Change Rep 5:245\u0026ndash;257. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s40641-019-00138-7\u003c/span\u003e\u003cspan address=\"10.1007/s40641-019-00138-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePlanton YY, Guilyardi E, Wittenberg AT, Lee J, Gleckler PJ, Bayr T, \u0026hellip; and, Voldoire A (2021) Evaluating climate models with the CLIVAR 2020 ENSO metrics package. Bull Am Meteorol Soc 102(2):E193\u0026ndash;E217. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/BAMS-D-19-0337.1\u003c/span\u003e\u003cspan address=\"10.1175/BAMS-D-19-0337.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePriya P, Dommenget D, McGregor S (2024) The dynamics of the El Ni\u0026ntilde;o\u0026ndash;Southern Oscillation diversity in the recharge oscillator framework. Clim Dyn 1\u0026ndash;21. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-024-07158-6\u003c/span\u003e\u003cspan address=\"10.1007/s00382-024-07158-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRayner NAA, Parker DE, Horton EB, Folland CK, Alexander LV, Rowell DP, Kent EC, Kaplan A (2003) Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. J Geophys Res Atmos 108(D14). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/2002JD002670\u003c/span\u003e\u003cspan address=\"10.1029/2002JD002670\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSun Y, Wang F, Sun DZ (2016) Weak ENSO asymmetry due to weak nonlinear air\u0026ndash;sea interaction in CMIP5 climate models. Adv Atmos Sci 33:352\u0026ndash;364. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00376-015-5018-6\u003c/span\u003e\u003cspan address=\"10.1007/s00376-015-5018-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTaylor KE, Stouffer RJ, Meehl GA (2012) An overview of CMIP5 and the experiment design. Bull Am Meteorol Soc 93(4):485\u0026ndash;498. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/BAMS-D-11-00094.1\u003c/span\u003e\u003cspan address=\"10.1175/BAMS-D-11-00094.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTimmermann A, An SI, Kug JS, Jin FF, Cai W, Capotondi A, \u0026hellip; and, Zhang X (2018) El Ni\u0026ntilde;o\u0026ndash;southern oscillation complexity. Nature 559(7715):535\u0026ndash;545. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41586-018-0252-6\u003c/span\u003e\u003cspan address=\"10.1038/s41586-018-0252-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTsonis AA, Hunt AG, Elsner JB (2003) On the relation between ENSO and global climate change. Meteorol Atmos Phys 84(3\u0026ndash;4) 229\u0026ndash;242, 3579, 3579. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00703-003-0001-7\u003c/span\u003e\u003cspan address=\"10.1007/s00703-003-0001-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVijayeta A, Dommenget D (2018) An evaluation of ENSO dynamics in CMIP simulations in the framework of the recharge oscillator model. Clim Dyn 51(5):1753\u0026ndash;1771. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-017-3981-6\u003c/span\u003e\u003cspan address=\"10.1007/s00382-017-3981-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang G, Dommenget D, Frauen C, Wang (2015) An evaluation of the CMIP3 and CMIP5 simulations in their skill of simulating the spatial structure of SST variability. Clim Dyn 44:95\u0026ndash;114. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-014-2154-0\u003c/span\u003e\u003cspan address=\"10.1007/s00382-014-2154-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang C, Deser C, Yu JY, DiNezio P, Clement A (2017) El Ni\u0026ntilde;o and southern oscillation (ENSO): a review. Coral reefs East Trop Pacific: Persistence loss dynamic Environ 85\u0026ndash;106. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/978-94-017-7499-4_4\u003c/span\u003e\u003cspan address=\"10.1007/978-94-017-7499-4_4\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWeller E, Cai W (2013) Asymmetry in the IOD and ENSO teleconnection in a CMIP5 model ensemble and its relevance to regional rainfall. J Clim 26(14):5139\u0026ndash;5149. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1175/JCLI-D-12-00789.1\u003c/span\u003e\u003cspan address=\"10.1175/JCLI-D-12-00789.1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWengel C, Dommenget D, Latif M, Bayr T, Vijayeta A (2018) What controls ENSO-amplitude diversity in climate models? Geophys Res Lett 45(4):1989\u0026ndash;1996. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/2017GL076849\u003c/span\u003e\u003cspan address=\"10.1002/2017GL076849\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWyrtki K (1985) Water displacements in the Pacific and the genesis of El Ni\u0026ntilde;o cycles. J Geophys Research: Oceans 90(C4):7129\u0026ndash;7132. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1029/JC090iC04p07129\u003c/span\u003e\u003cspan address=\"10.1029/JC090iC04p07129\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang C, Giese BS (2013) El Ni\u0026ntilde;o Southern Oscillation in an ensemble ocean reanalysis and coupled climate models. J Geophys Research: Oceans 118(9):4052\u0026ndash;4071. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/jgrc.20284\u003c/span\u003e\u003cspan address=\"10.1002/jgrc.20284\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYu Y, Dommenget D, Frauen C, Wang G, Wales S (2016) ENSO dynamics and diversity resulting from the recharge oscillator interacting with the slab ocean. Clim Dyn 46(5):1665\u0026ndash;1682. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00382-015-2667-1\u003c/span\u003e\u003cspan address=\"10.1007/s00382-015-2667-1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZuo H, Balmaseda MA, Tietsche S, Mogensen K, Mayer M (2019) The ECMWF operational ensemble reanalysis\u0026ndash;analysis system for ocean and sea ice: a description of the system and assessment. Ocean Sci 15(3):779\u0026ndash;808. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.5194/os-2018-154\u003c/span\u003e\u003cspan address=\"10.5194/os-2018-154\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"El Niño-Southern Oscillation (ENSO) dynamics, recharge oscillator model, ENSO phase space, ENSO events, CMIP5 and CMIP6 models, non-linear dynamics","lastPublishedDoi":"10.21203/rs.3.rs-4727039/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4727039/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThis study analyses the El Ni\u0026ntilde;o-Southern Oscillation (ENSO) phase space as simulated by the Coupled Model Intercomparison Projects 5 and 6 (CMIP5 and CMIP6) models. The ENSO phase space describes the ENSO cycle between the sea surface temperature (SST) anomaly in the eastern equatorial Pacific (\u003cb\u003eT\u003c/b\u003e) and the equatorial mean thermocline depth anomaly (\u003cb\u003eh\u003c/b\u003e). We find that the characteristics out-of-phase cross-correlation between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e is shifted to negative values in CMIP models, suggesting that the coupling between \u003cb\u003eT\u003c/b\u003e and \u003cb\u003eh\u003c/b\u003e is regionally sifted to the east compared to the observed central Pacific. If we consider the CMIP models with an eastward shifted \u003cb\u003eh\u003c/b\u003e then the models have better agreements with the observed characteristics. While the models can capture some of the non-linear aspects with high correlations, they do largely underestimate the strength of non-linear ENSO aspects. They underestimate the likelihood of extreme El Ni\u0026ntilde;o and discharge states, they cannot capture the enhanced growth rates during the recharge state, the enhanced decay after the discharge state nor the reduced phase transitions after the La Ni\u0026ntilde;a phases. Weaker than observed wind-SST feedback and weaker \u003cb\u003eh\u003c/b\u003e variability are likely some of the reasons why models cannot fully capture the non-linear ENSO phase space dynamics. Further, we found no indication of significant improvements from the CMIP 5 to 6 ensemble, suggesting that the two ensembles are essentially the same in terms of their ENSO dynamics. There is, however, a large spread within the model ensembles, leading to models with quite different ENSO dynamics.\u003c/p\u003e","manuscriptTitle":"ENSO phase space dynamics in CMIP models","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-08-02 18:33:40","doi":"10.21203/rs.3.rs-4727039/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"5e41d313-621c-416e-9058-5c8c0b4ffaf2","owner":[],"postedDate":"August 2nd, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2024-11-15T06:36:35+00:00","versionOfRecord":[],"versionCreatedAt":"2024-08-02 18:33:40","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4727039","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4727039","identity":"rs-4727039","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}