1. Introduction
The spring persistence barrier (SPB) and the associated spring predictability barrier have received considerable attention in El Niño–Southern Oscillation (ENSO) prediction studies. It has long been noticed in the observations of the tropical Pacific Ocean sea surface temperature (SST) anomaly (e.g., Niño-3.4; 5°S–5°N, 170°–120°W; Ren et al. 2016) variability and the variability in sea level pressure (Troup 1965; Webster and Yang 1992) and rainfall (Walker and Bliss 1932; Wright 1979). The major feature of this SPB is that a band of maximum decline of monthly autocorrelation occurs at a fixed phase (or month), as seen from the monthly autocorrelation of SST anomaly variability and its lag gradient (e.g., Liu et al. 2019). This maximum decline of persistence of fixed phase shows that, regardless of the initial months, a damped persistence forecast loses its predictability most rapidly in the following May or June, forming the SPB of ENSO. Recent studies also suggested that there are distinct features of persistence barrier in the two types of ENSO, the eastern Pacific (EP) and central Pacific (CP) types (Ren et al. 2016). By using a statistical model, Ren et al. (2019) found that the summer barrier in the ENSO persistence for the CP type is closely related to the summer predictability barrier.
Despite intensive studies of the SPB of ENSO, the nature of the SPB has not been fully understood. Some studies suggested that the spring predictability barrier is related to systematic errors of the forecasting scheme and therefore is not an intrinsic feature of ENSO (Chen et al. 2004; Johnson et al. 2000; Yu et al. 2012), whereas some other studies demonstrated that the SPB or spring predictability barrier is an intrinsic feature in the real climate system related to the seasonally varying initial error growth or ENSO growth rate of SST anomaly, which is in turn related to the background of the tropical Pacific (Blumenthal 1991; Moore and Kleeman 1996; Torrence and Webster 1998; Mu et al. 2007; Levine and McPhaden 2015). Zhu et al. (2015) indicated that, over the eastern equatorial Pacific, the connection between the thermocline and SST is the weakest during boreal spring, which may be one of causes for ENSO SPB. Recently, Liu et al. (2019, hereinafter LJ19) derived a simple analytical solution for the SPB as an asymptotic solution of the damped persistence model with seasonal growth rate and noise forcing—that is, a Langevin equation of seasonally varying coefficients. This solution demonstrated explicitly that once the amplitude of the seasonal growth rate exceeds a threshold, the seasonal persistence barrier occurs. It showed theoretically that the persistence barrier is an intrinsic feature in a seasonal forcing system, at least in a system described by the Langevin equation. By linking the seasonal growth rate in this equation with the seasonally varying Bjerknes instability index (Jin et al. 2006), Jin et al. (2019) indicated that seasonal thermodynamic damping and thermocline positive feedback play an important role in determining the ENSO SPB. The role of ENSO growth rate in SPB strength (i.e., a weaker ENSO growth rate strengthens SPB) is further identified in two different ENSO regimes (Jin et al. 2020). Overall, the role of ENSO growth rate in SPB has been discussed by recent studies.
One major deficiency of the theory of LJ19 and Jin et al. (2019, 2020), however, is the lack of the effect of the natural ENSO period on SPB, because the Langevin equation does not include an ENSO period. Because the ENSO period has been suggested to be important for its predictability (Chen et al. 1997; Chen and van den Dool 1997), it is interesting to examine the role of an inherent ENSO period in SPB. More specifically, we are interested in the following question: How does the ENSO period influence the strength and phase of the SPB?
This paper is an attempt to understand the modulation of ENSO period on SPB mainly in the context of the neutral recharge oscillator (NRO) model (Stein et al. 2014). Both analytical and numerical solutions of the autocorrelation function (ACF) and SPB are derived in the NRO model. The features of SPB caused by ENSO period are also observed and interpreted in the damped regime, self-exciting regime, and ENSO observations.
The paper is arranged as follows. The recharge oscillator model, the definition of ENSO SPB strength and timing, and the reanalysis data will be presented in section 2. In sections 3 and 4, we will explore the effect of ENSO period on SPB in the NRO model in analytical and numerical solutions, respectively. Section 5 will further examine the SPB for ENSO in the damped and self-exciting regime as well as in the observations. A summary and discussion will be given in section 6.
2. Model, method, and data
a. The parametric recharge oscillator model
An approximate analytical solution of SPB is obtained from the neutrally stable, unforced case of the recharge model with R0 = 0, σ = 0, c = 0, and b = 0. Previous studies have shown that the approximate analytical solution of variance in this neutral model is able to explain the ENSO phase locking of variance (Stein et al. 2014). Here we will derive the corresponding ACF and in turn study the SPB. The robustness of the results of the NRO model is then confirmed numerically in the recharge model in the damped regime with noise forcing (R0 > 0, σ > 0, c = 0, and b = 0) or self-exciting unstable regime (R0 < 0, σ = 0, c > 0, and b = 0). All of the numerical results are from the last 500 years of 1000 model years. The numerical model Eqs. (2.1)–(2.3) are solved at a time step of 4 h.
b. Definition of SPB strength and timing
c. The reanalysis dataset
The reanalysis dataset of SST used here in this study is the monthly Hadley Centre Global Sea Ice and Sea Surface Temperature (HadISST; 1° × 1°). Here we use HadISST v1.1 (https://climatedataguide.ucar.edu/climate-data/sst-data-hadisst-v11) to identify some general features of SPB. All monthly data used here with their climatological seasonal cycle and linear trends being removed.
3. Effect of ENSO period on SPB in the NRO model: Analytical solution
In this section, the role of ENSO period in SPB is studied by deriving an approximate analytical solution in the NRO model.
The (a) numerical and (b) analytical solutions of the persistence map (ACF) for λε = 1/π and ω0 = (2π/48) month−1 in the NRO model. The black circles on the persistence map mark the lag month of maximum autocorrelation decline for different initial months.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
The solution of Eq. (3.11) shows explicitly that, when τB increases, the timing m* + τB is no longer a fixed phase; instead, it shifts earlier as |E/τB| becomes smaller. Therefore, in contrast to the Langevin equation in which the timing of SPB is strictly phase locked to a fixed month (Liu et al. 2019), the timing of SPB in NRO model varies in a range. That is, γ is no longer constant, which can be seen clearly in section 4. This may be caused by the interaction between the annual frequency and lower frequency (2–7 yr; ENSO period).
We examine the SPB for an entire annual cycle of initial months to explore the phase range of SPB. This feature indicates that although SPB occurs in May or June (Liu et al. 2019) for initial month January, when we start to forecast in October, the forecast skill will drop dramatically in the early spring, shorter than we expect. We will explain this feature through the numerical solutions.
4. Effect of ENSO period on SPB in the NRO model: Numerical solution
In this section, the features of SPB in light of the analytical solution derived above, and furthermore in the numerical solution of NRO model, are discussed.
a. Impact of ENSO period on the intensity of SPB
We first show that a shorter ENSO period strengthens the SPB. According to Eq. (3.7) (gradient of the ACF), ω* sinω*τ* is proportional to ω*. The last two terms, which contain the seasonal growth rate, are the main reasons for SPB. As
The numerical solution of the persistence map (ACF) for λε = 1.2/π, in the NRO model. (a) The ω0 = (2π/36) month−1 case. The black lines on the persistence map mark the lag phase of maximum autocorrelation decline for different initial phases [i.e., τB(m)]. The gray line indicates that the phase of SPB (γ) is constant. Here the daily data are used to calculate the persistence map. (b) The relationship between γ and the maximum persistence decline τB. The blue plus signs indicate γ for the ω0 = (2π/36) month−1 case, and the black line suggests γ for the Langevin equation. (c),(d) As in (a) and (b), but for the ω0 = (2π/60) month−1 case.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
To confirm the impact of ENSO period on SPB strength systematically, we examine the modulation of the SPB strength by varying the ENSO period from 2.5 to 6 years (Fig. 3). In the NRO model, there is a parametric resonance when ENSO period is two years (Stein et al. 2014) such that we study the ENSO period from 2.5 years. For a small amplitude of seasonal cycle of damping rate λε = 0.6/π, both numerical solutions (blue triangles in Fig. 3) and approximate analytical solutions (green triangles in Fig. 3) show that a shorter ENSO period strengthens ENSO SPB. The difference between the analytical and numerical solutions increases modestly when the ENSO period becomes shorter, because the approximate analytical solution is derived in the limit of ENSO period much larger than the annual cycle. This weakening effect of a longer ENSO period on SPB is also robust for an increased amplitude of seasonal cycle: λε = 1.8/π (squares in Fig. 3). The difference between the numerical solution and analytical solution is larger than that of λε = 0.6/π case, because the analytical solution is derived in the limit of small amplitude of seasonal cycle λε ≪ 1. The comparison of the cases of λε = 0.6/π and 1.8/π shows that the SPB is stronger for larger λε (blue triangles vs blue squares or green triangles vs green squares in Fig. 3). This is consistent with the Langevin equation where the SPB intensifies with the amplitude of the seasonal period (LJ19). Note here that previous studies have shown the relationship between the SPB intensity and ENSO phase locking (Torrence and Webster 1998; Tian et al. 2019). According to Eq. (3.3), ENSO phase locking strength is directly related to the seasonal forcing amplitude λε. Moreover, the strength of ENSO SPB is also controlled by the seasonal forcing amplitude (triangles vs squares in Fig. 3; Liu et al. 2019). Both ENSO SPB and phase locking strength are related to this amplitude, so both will show similar changes.
The relationship between ENSO period (yr) and SPB strength for λε = 0.6/π (for clarity, we use λ0 = 0.1 month−1 in the figure legend; triangles) and λε = 1.8/π (here we use λ0 = 0.3 month−1 in the figure legend; squares); the blue symbols are the numerical solution, and the green symbols indicate the results from the analytical solution of ACF.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
b. Impact of ENSO period on the phase of SPB
We now show that the timing of SPB in the NRO model is no longer phase locked to a fixed month. Instead, the SPB varies with a range of about 2 months, for some typical parameters. This is in contrast to the Langevin equation, where the timing of SPB is phase locked to a fixed month once SPB occurs. This change of the phase of SPB is caused by the presence of the ENSO period. As seen in Eq. (3.11), γ is no longer a constant in general, which can be shown in the examples in Fig. 2. For an ENSO period of 3 years (Fig. 2a), the first SPB lag τB0 is 75° at the initial phase of 145°. According to Eq. (2.6), the timing of SPB is 75° + 145° = 220°. With the increasing lag τB to, say, 190° at the initial month of 0°, the timing of the SPB is therefore 0° + 190° =190° and is shifted earlier by 30° (~1 month). With the further increased τB, the phase of the SPB is shifted even earlier. In particular, the relationship between the initial month m and maximum persistence decline τB is not linear (gray line in Fig. 2a), bending downward (black line in Fig. 2a). This change of SPB month can be seen more clearly in Fig. 2b (blue plus sign). For the 3-yr ENSO case, at τB = τB0, the SPB phase is γ(τB0) = 210°. As τB increases (m* is modulated as well), γ decreases in the NRO model, in contrast to the Langevin equation in which the SPB month maintains at 210° (black line in Fig. 2b), and the timing of SPB decreases by ~60° (from 210° to 150°; ~2 months). This feature of changing SPB timing is also robust for the case of a 5-yr ENSO period (Fig. 2c, and blue plus sign in Fig. 2d), although the phase range is reduced to ~40° (from 213° to 173°; ~1.3 months). Indeed, the Langevin equation can be treated as the case of ENSO period of infinitely long period (m* = 0), and the phase range is then reduced to zero (i.e., an exact phase locking), as shown in LJ19.
The phase range of SPB with ENSO period can be seen more systematically in the NRO model in both analytical and numerical solutions (Fig. 4), in which the ENSO period changes from 2.5 to 6 years. For the case of λε = 0.6/π, when the ENSO period is lengthened, the phase range of the SPB decreases monotonically in the numerical solution (blue triangles in Fig. 4). This decrease of phase range is well approximated by the analytical solution (except for very short ENSO periods; green triangle in Fig. 4). This is, in fact, straightforward to understand. When the ENSO period is longer (ω* → 0), the NRO model behaves more like the Langevin equation, in which the phase range vanishes (LJ19). The phase changes of SPB with the modulation of ENSO period may be caused by the interaction between annual frequency and lower frequency. When ENSO period is shorter, this interaction is stronger and shifts the timing of SPB. Physically, the shortened ENSO period may be caused by the thermocline and zonal advective feedbacks (Lu et al. 2018). As such, these feedbacks may also affect the timing of SPB. In the case of a larger amplitude of seasonal cycle λε = 1.8/π, the phase range is reduced (green square vs green triangle or blue square vs blue triangle in Fig. 4). This may be understood as follows. A larger λε represents a stronger seasonal forcing. As indicted by LJ19, the SPB is caused by this seasonal forcing and intensified with the seasonal cycle. When the seasonal forcing is intensified in the NRO model, the phase locking of the SPB to a fixed phase becomes more distinct, with the range of SPB phase decreased.
As in Fig. 3, but for the phase range of SPB.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
In sum, our study of the NRO model here shows that a shorter ENSO period tends to intensify the SPB and shift the timing of the SPB. These features will also be seen in the next section for ENSOs in the damped persistence regime, self-exciting regime, and observation.
5. Impact of ENSO period on the SPB features in the different ENSO regimes and observation
We now examine SPB response to ENSO period in the recharge oscillator model in the damped regime (R0 > 0, σ > 0, c = 0, and b = 0), self-exciting regime (R0 < 0, σ = 0, c > 0, and b = 0), and observation.
a. SPB for ENSO in the damped-persistence regime
In the damped regime, a shorter ENSO period will also strengthen the SPB. Parameters set for the observational ENSO are (Chen and Jin 2020) R0 = −0.1 month−1, σ = 1/9 month−1, and d = 0.66 month−1. Here we set λε = 1.8/π to show distinct ENSO SPB features. When ENSO period is 3 years, the SPB strength is 4.76 (Fig. 5a), which is significantly stronger than the case when ENSO period is 5 years (2.61; Fig. 5b). It shows that a shorter ENSO period will also strengthen SPB in the damped regime.
The SPB features in the damped recharge oscillator model. The numerical solution of persistence map for ENSO period is (a) 3 and (b) 5 years. The red line indicates an SPB timing, and the black line suggests another SPB timing.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
Furthermore, ENSO SPB is no longer phase locked to a fixed month in the damped regime as in the neutral case before. Here we use monthly data in the damped recharge oscillator model to directly compare it with the observation. When initial month is April (Fig. 5a), the maximum persistence decline τB occurs in the lag-2 month, which indicates that SPB occurs in June (red line in Fig. 5a). However, as the τB increases, for initial month January, SPB happens in May (black line in Fig. 5a). The timing of ENSO SPB moves forward, which is consistent with SPB features in NRO model. The sudden jump of the SPB month is due to this monthly SST data used here. This feature of SPB also occurs when ENSO period is 5 years (Fig. 5b). It is noted that, as compared with the neutral oscillation in NRO model, damping produces one realistic feature of the SPB: the emergence of the SPB from the initial lag (i.e., τB0 = 0). This is consistent with the study in the Langevin equation, because the decline of persistence in the initial lags depends critically on the damping of the system (LJ19).
To confirm the impact of ENSO period on SPB features systematically in the damped regime, we vary the ENSO period from 2.5 to 6 years, as shown in Fig. 6. We find that when the ENSO period is shorter, the SPB strength is stronger (squares or triangles in Fig. 6a). Meanwhile, ENSO SPB is stronger for a larger seasonal amplitude λε (squares vs triangles in Fig. 6a). It is consistent with the results in NRO model. We also find that a phase range of SPB in the damped regime exists (1–2 months). When λε is larger, the phase range of SPB is small (square vs triangle in Fig. 6b), which is the same as the results in the NRO model.
(a),(b) As in Figs. 3 and 4, respectively, but for the numerical solution in the damped recharge oscillator model. The unit in (b) is degrees, which is consistent with Fig. 4 (30° represents one month).
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
b. SPB for ENSO in the self-exciting regime
In the self-exciting regime, we set the parameter for realistic ENSO as R0 = 0.1 month−1, c = 1/62.3 month−1, and b = 0 (Chen and Jin 2020). Here we set λε = 1.8/π to be consistent with the damped regime. When the ENSO period is 3 years (Fig. 7a), the SPB strength is 4.59, which is significantly stronger for the case when ENSO period is 5 years (1.69; Fig. 7b). It shows that the SPB weakens with the increased period of ENSO period. Meanwhile, we find both in two cases, the timing of SPB is not phase locked (Fig. 7). For initial month May (Fig. 7a), the maximum persistence decline τB occurs in lag 2-month such that the month of SPB is July (red line in Fig. 7a). When τB is increased, for initial month January, τB is 5 months such that the month of SPB is June (black line in Fig. 7a), which is one month earlier for the case of initial month May. We also find similar SPB features when ENSO period is 5 years (Fig. 7b). In all, the phase of SPB moves forward as τB is increased in the unstable regime.
As in Fig. 5, but for the SPB features in the unstable ENSO regime.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
We also examine the SPB features systematically in the unstable regime by varying the ENSO period from 2.5 to 6 years (Fig. 8). It is found that the SPB is strengthened with decreased ENSO period (squares or triangles in Fig. 8a). Meanwhile, the SPB exists for a phase range in the unstable regime (square or triangle in Fig. 8b). When seasonal forcing amplitude is decreased, ENSO SPB strength is weaker (squares vs triangles in Fig. 8a) and the phase range of SPB is larger (squares vs triangles in Fig. 8b). In all, the SPB features caused by ENSO period are consistent with the results in NRO model.
As in Fig. 6, but for the SPB features in the unstable ENSO regime.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
c. SPB phase shift for observed ENSO
The earlier shifting of SPB can be identified in the observational ENSO. During the year 1958–80, for the initial lags of 2–3 months, which corresponds to the initial months of March to April, the timing of the SPB occurs in June (red line in Fig. 9a). However, for increased lag τB and the initial months of May–December, the SPB is shifted one month earlier to May (black line in Fig. 9a). This change of SPB month can also be seen for the years 1980–2016, when the timing of SPB also shifts earlier for increased lags, from June for the initial months of January to April (red line in Fig. 9b) to May for initial months of May to November (black line in Fig. 9b). This shift of the SPB month is not caused by the sampling error, since for the late half year (e.g., Fig. 9b) the timing of SPB moves forward systematically. Instead, our theory suggests that this shift is an intrinsic feature of the observed ENSO and it identifies that the timing of SPB will be shifted earlier from the early summer (May–June) to late spring (April–May), corresponding to the initial months of the early half year and later half year, respectively. This feature indicates that for later half year, ENSO predictability decreases as the cause of ENSO period.
As in Fig. 5, but for observational data during (a) 1958–80 and (b) 1980–2016.
Citation: Journal of Climate 34, 6; 10.1175/JCLI-D-20-0540.1
In this section, we identify that the features of the SPB in the NRO model are robust for the damped persistence regime, self-exciting regime, and observational ENSO. We find that a shorter ENSO period will strengthen ENSO SPB. Moreover, the timing of SPB is no longer phase locked to a specific month. As τB is increased, the timing of SPB shifts earlier, indicating that ENSO predictability decreases for later initial months.
6. Summary and discussion
As a follow up of LJ19, this paper attempts to understand the role of ENSO period in SPB. In particular, on the basis of the approximate analytical solution of SST anomaly variance in the NRO model (Stein et al. 2014), we derive the ACF and develop a theory that sheds light on the influence of ENSO period on SPB. We find two effects of ENSO period on the SPB. First, a shorter ENSO period strengthens the SPB, implying that the prediction skill will drop more significantly if the ENSO period is shortened. Second, in contrast to the SPB in the Langevin equation (LJ19), the timing of SPB is no longer phase locked exactly to a fixed month. Instead, the timing of the SPB shifts earlier by 1–2 months. Therefore, the maximum drop of ENSO prediction skill will depend on the forecast lead and initial forecast month. The forecast skill will drop dramatically in late spring for small forecast leads when the predictions start from the early half of the year, whereas the prediction skill will drop significantly in early spring for long-term leads of forecast when the prediction starts from the late half of the year. These two impacts of ENSO period on SPB are further found to be robust in the recharge model for ENSOs in the damped regime, self-exciting regime, and observation.
The physical mechanism of influencing ENSO SPB strength is complicated according to the analytical solution of NRO model and Bjerknes instability index (Jin et al. 2006). On the basis of the recharge oscillator framework, Lu et al. (2018) indicated that the thermocline and zonal advective feedbacks multiplied by the efficiency factor of discharging/recharging of the equatorial content driven by ENSO wind stress anomalies [ω0 in Eqs. (2.1) and (2.2), respectively] control ENSO linear periodicity. As demonstrated by Eqs. (2.4) and (2.5), a larger ω0 leads to a shorter ENSO period and a strengthened SPB. Accordingly, the thermocline and zonal advective feedbacks may affect the strength of SPB by controlling the ENSO period. On the other hand, the Bjerknes instability index is identified to play an important role in the strength of SPB (Jin et al. 2019, 2020). The thermocline and zonal advective feedbacks can also affect SPB strength by controlling this Bjerknes instability index. Therefore, the physical mechanism of affecting SPB strength still needs further study.
This study may have implications for ENSO prediction under a global warming scenario. As suggested by Collins et al. (2010), year-to-year ENSO variability is controlled by a delicate balance of amplifying and damping feedbacks, and it is not yet possible to determine whether ENSO period will change. As such, ENSO SPB, and in turn ENSO predictability, may not be strengthened or weakened as ENSO period will not change dramatically.
Overall, our theory provides an explanation of the role of ENSO period on SPB. We suggest that it can serve as a null hypothesis for the ENSO SPB.
Acknowledgments
This work is supported by Chinese MOST 2017YFA0603801 and NSFC41630527 and U.S. National Science Foundation AGS-1656907.
APPENDIX
REFERENCES
An, S.-I., 2008: Interannual variations of the tropical ocean instability wave and ENSO. J. Climate, 21, 3680–3686, https://doi.org/10.1175/2008JCLI1701.1.
An, S.-I., and F.-F. Jin, 2011: Linear solutions for the frequency and amplitude modulation of ENSO by the annual cycle. Tellus, 63A, 238–243, https://doi.org/10.1111/j.1600-0870.2010.00482.x.
Blumenthal, M. B., 1991: Predictability of a coupled ocean–atmosphere model. J. Climate, 4, 766–784, https://doi.org/10.1175/1520-0442(1991)004<0766:POACOM>2.0.CO;2.
Chen, D., S. E. Zebiak, M. A. Cane, and A. J. Busalacchi, 1997: Initialization and predictability of a coupled ENSO forecast model. Mon. Wea. Rev., 125, 773–788, https://doi.org/10.1175/1520-0493(1997)125<0773:IAPOAC>2.0.CO;2.
Chen, D., M. A. Cane, A. Kaplan, S. E. Zebiak, and D. Huang, 2004: Predictability of El Niño over the past 148 years. Nature, 428, 733–736, https://doi.org/10.1038/nature02439.
Chen, H.-C., and F.-F. Jin, 2020: Fundamental behavior of ENSO phase locking. J. Climate, 33, 1953–1968, https://doi.org/10.1175/JCLI-D-19-0264.1.
Chen, W. Y., and H. M. van den Dool, 1997: Atmospheric predictability of seasonal, annual, and decadal climate means and the role of the ENSO cycle: A model study. J. Climate, 10, 1236–1254, https://doi.org/10.1175/1520-0442(1997)010<1236:APOSAA>2.0.CO;2.
Collins, M., and Coauthors, 2010: The impact of global warming on the tropical Pacific Ocean and El Niño. Nat. Geosci., 3, 391–397, https://doi.org/10.1038/ngeo868.
Jin, F.-F., 1997a: An equatorial ocean recharge paradigm for ENSO. Part I: Conceptual model. J. Atmos. Sci., 54, 811–829, https://doi.org/10.1175/1520-0469(1997)054<0811:AEORPF>2.0.CO;2.
Jin, F.-F., 1997b: An equatorial ocean recharge paradigm for ENSO. Part II: A stripped-down coupled model. J. Atmos. Sci., 54, 830–847, https://doi.org/10.1175/1520-0469(1997)054<0830:AEORPF>2.0.CO;2.
Jin, F.-F., S. T. Kim, and L. Bejarano, 2006: A coupled-stability index for ENSO. Geophys. Res. Lett., 33, L23708, https://doi.org/10.1029/2006GL027221.
Jin, Y., Z. Liu, Z. Lu, and C. He, 2019: Seasonal period of background in the tropical Pacific as a cause of ENSO spring persistence barrier. Geophys. Res. Lett., 46, 13 371–13 378, https://doi.org/10.1029/2019GL085205.
Jin, Y., Z. Lu, and Z. Liu, 2020: Controls of spring persistence barrier strength in different ENSO regimes and implications for 21st century changes. Geophys. Res. Lett., 47, e2020GL088010, https://doi.org/10.1029/2020GL088010.
Johnson, S. D., D. S. Battisti, and E. Sarachik, 2000: Empirically derived Markov models and prediction of tropical Pacific sea surface temperature anomalies. J. Climate, 13, 3–17, https://doi.org/10.1175/1520-0442(2000)013<0003:EDMMAP>2.0.CO;2.
Levine, A. F., and M. J. McPhaden, 2015: The annual cycle in ENSO growth rate as a cause of the spring predictability barrier. Geophys. Res. Lett., 42, 5034–5041, https://doi.org/10.1002/2015GL064309.
Liu, Z., Y. Jin, and X. Rong, 2019: A theory for the seasonal predictability barrier: Threshold, timing, and intensity. J. Climate, 32, 423–443, https://doi.org/10.1175/JCLI-D-18-0383.1.
Lu, B., F.-F. Jin, and H.-L. Ren, 2018: A coupled dynamic index for ENSO periodicity. J. Climate, 31, 2361–2376, https://doi.org/10.1175/JCLI-D-17-0466.1.
Moore, A. M., and R. Kleeman, 1996: The dynamics of error growth and predictability in a coupled model of ENSO. Quart. J. Roy. Meteor. Soc., 122, 1405–1446, https://doi.org/10.1002/qj.49712253409.
Mu, M., H. Xu, and W. Duan, 2007: A kind of initial errors related to “spring predictability barrier” for El Niño events in Zebiak–Cane model. Geophys. Res. Lett., 34, L03709, https://doi.org/10.1029/2006GL027412.
Ren, H.-L., F.-F. Jin, B. Tian, and A. A. Scaife, 2016: Distinct persistence barriers in two types of ENSO. Geophys. Res. Lett., 43, 10 973–10 979, https://doi.org/10.1002/2016GL071015.
Ren, H.-L., J. Zuo, and Y. Deng, 2019: Statistical predictability of Niño indices for two types of ENSO. Climate Dyn., 52, 5361–5382, https://doi.org/10.1007/s00382-018-4453-3.
Stein, K., A. Timmermann, N. Schneider, F.-F. Jin, and M. F. Stuecker, 2014: ENSO seasonal synchronization theory. J. Climate, 27, 5285–5310, https://doi.org/10.1175/JCLI-D-13-00525.1.
Tian, B., H.-L. Ren, F.-F. Jin, and M. F. Stuecker, 2019: Diagnosing the representation and causes of the ENSO persistence barrier in CMIP5 simulations. Climate Dyn., 53, 2147–2160, https://doi.org/10.1007/s00382-019-04810-4.
Torrence, C., and P. J. Webster, 1998: The annual period of persistence in the El Niño/Southern Oscillation. Quart. J. Roy. Meteor. Soc., 124, 1985–2004, https://doi.org/10.1002/qj.49712455010.
Troup, A. J., 1965: The “southern oscillation.” Quart. J. Roy. Meteor. Soc., 91, 490–506, https://doi.org/10.1002/qj.49709139009.
Walker, G. T., and E. W. Bliss, 1932: World weather. V. Mem. Roy. Meteor. Soc., 4, 53–84.
Webster, P., and S. Yang, 1992: Monsoon and ENSO: Selectively interactive systems. Quart. J. Roy. Meteor. Soc., 118, 877–926, https://doi.org/10.1002/qj.49711850705.
Wright, P. B., 1979: Persistence of rainfall anomalies in the central Pacific. Nature, 277, 371–374, https://doi.org/10.1038/277371a0.
Yu, Y., M. Mu, and W. Duan, 2012: Does model parameter error cause a significant “spring predictability barrier” for El Niño events in the Zebiak–Cane model? J. Climate, 25, 1263–1277, https://doi.org/10.1175/2011JCLI4022.1.
Zhu, J., A. Kumar, and B. Huang, 2015: The relationship between thermocline depth and SST anomalies in the eastern equatorial Pacific: Seasonality and decadal variations. Geophys. Res. Lett., 42, 4507–4515, https://doi.org/10.1002/2015GL064220.