Zhaolu Hou

El Niño and the Southern Oscillation (ENSO) have a worldwide impact on seasonal to yearly climate. However, there are decadal variations in the seasonal prediction skill of ENSO in dynamical and statistical models; in particular, ENSO prediction skill has declined since 2000. The shortcomings of models mean that it is very important to study ENSO seasonal predictability and its decadal variation using observational/reanalysis data. Here we quantitatively estimate the seasonal predictability limit (PL) of ENSO from 1900 to 2015 using Nonlinear local Lyapunov exponent (NLLE) theory with an observational/reanalysis dataset and explore its decadal variations. The mean PL of sea surface temperature (SST) is high in the central/eastern tropical Pacific and low in the western tropical Pacific, reaching 12–15 and 7–8 months, respectively. The PL in the tropical Pacific varies on a decadal timescale, with an interdecadal standard deviation of up to 2 months in the central tropical Pacific that has similar spatial structure to the mean PL. Taking the PL of SST in the Niño 3.4 region as representative of the PL in the central/eastern tropical Pacific, there are clearly higher values in the 1900s, mid-1930s, mid-1960s, and mid-1990s, and lower values in the 1920s, mid-1940s, and mid-2010s. Meanwhile, the PL of SST in the Niño 6 region—whose average value is 7 months—is in good agreement with the PL of most regions in the western tropical Pacific, with higher values in the 1910s, 1940s, and 1980s and lower values in the 1930s, 1950s, and mid-1990s.In the framework of NLLE theory, the PL is determined by the error growth rate (representing the dissipation rate of the predictable signal) and the saturation value of relative error (representing predictable signal intensity). We reveal that the spatial structure of the mean PL in the tropical Pacific is determined mainly by the error growth rate. The decadal variability of PL is affected more by the variation of the saturation value of relative error in the equatorial Pacific, whereas the error growth rate cannot be ignored in the PL of some regions. As an important source of predictability in ENSO dynamics, the relationship between warm water volume and SST in the Niño 3.4 region has a critical role in the decadal variability of PL in the tropical Pacific through the error growth rate and saturation value of relative error. This strong relationship reduces the error growth rate in the initial period and increases the saturated relative error, contributing to the high PL.

This study investigates the inter-seasonal variability of predictability limits (PLs) for South China seasonal precipitation (SCSP) using reanalysis data, emphasizing the crucial role of the western North Pacific anomalous anticyclone (WNPAC) in modulating SCSP predictability through ENSO-driven processes. Incorporating the relative global attractor radius as the saturation threshold into the nonlinear local Lyapunov exponent method to address SCSP’s oscillatory error growth after saturation reveals that its PL peaks in winter (3.72 months) and drops to its lowest during the autumn (3.12 months). ENSO is identified as the primary source of SCSP predictability, with its PLs generally extending longer and exhibiting distinct inter-seasonal variability compared to those of SCSP. As a crucial system linking ENSO and SCSP, WNPAC mirrors with SCSP in both PLs’ values and inter-seasonal variability. Physical analysis shows that seasonal shift in the WNPAC-ENSO relationship exerts a strong influence on SCSP predictability.

The Indian Ocean dipole (IOD) is a remarkable interannual variability in the tropical Indian Ocean. The improved prediction of IOD is of a great value because of its large socioeconomic impacts. Previous studies reported that both El Ni˜ no-Southern Oscillation (ENSO) and South China Sea summer monsoon (SM) play a dominant role in the western and eastern pole of the IOD, respectively. They can be used as predictors of the IOD at 3 month lead beyond self-persistence. Here, we develop an empirical model of multi-factors in which the western pole is predicted by ENSO and persistence and the eastern pole is predicted by SM and persistence. This new empirical model outperforms largely the average level of the dynamical models from the North American multi-model ensemble (NMME) project in predicting the peak IOD in boreal autumn, with a correlation coefficient of ∼0.86 and a root mean square error of ∼0.24°C. Furthermore, the hit rate of positive culminated IOD in this new empirical model is equivalent to that in current NMME models (above 65%), much higher than that for negative culminated IOD. This improvement of skill using the empirical model suggests a perspective for better understanding and predicting the IOD.

The North Pacific sea surface temperature (SST) exerts profound climatic influence. El Niño-Southern Oscillation (ENSO) significantly impacts North Pacific SST, yet the influence from ENSO’s distinct phases on SST predictability remains unclear. Overcoming model limitations, this study assesses SST predictability under diverse ENSO phases using reanalysis. Quantifying predictability limits (PL), results unveil asymmetry: El Niño PL at 5.5 months, La Niña at 8.4 months, and Neutral at 5.9 months. This asymmetry mirrors contemporary multimodal prediction skills. Error growth dynamics reveal La Niña's robust signal strength with slow error growth rate, contrasting El Niño's weaker signal and faster error growth. Neutral exhibits intermediate signal strength and elevated error growth. Physically, predictability signal strength aligns with SST variability, whereas error growth rate correlates with atmospheric-ocean heating anomalies. La Niña, inducing positive heating anomalies, minimizes atmospheric noise impact, resulting in lower error growth. The results are beneficial for improving North Pacific SST predictions.

Winter precipitation anomalies in South China (SC) frequently result in severe disasters. However, the evaluation of prediction performance and distinctions between positive precipitation anomaly events (PPA, wet condition) and negative precipitation anomaly events (NPA, dry condition) in current operational models remains incomplete. This study employed the Climate Forecast System version 2 (CFSv2) to assess winter precipitation prediction accuracy in SC from 1983 to 2021. Differences in predicting PPA and NPA events and the underlying physical mechanisms were explored. The results indicate that CFSv2 can effectively predict interannual variations in winter precipitation in SC, as there is a significant time correlation coefficient of 0.68 (0.62) between observations and predictions, with a lead time of 0 (3) months. The model revealed an intriguing asymmetry in prediction skills: PPA outperformed NPA in both deterministic and probabilistic prediction. The higher predictability of PPA, as indicated by the perfect model correlation and signal-to-noise ratio, contributed to its superior prediction performance when compared to NPA. Physically, tropical signals from the ENSO and extratropical signals from the Arctic sea ice anomaly, were found to play pivotal roles in this asymmetric feature. ENSO significantly impacts PPA events, whereas NPA events are influenced by a complex interplay of factors involving ENSO and Arctic sea ice, leading to low NPA predictability. The capability of the model to replicate Arctic sea ice signals is limited, but it successfully predicts ENSO signals and reproduces their related circulation responses. This study highlights the asymmetrical features of precipitation prediction, aiding in prediction models improvement.

The North Pacific sea surface temperature (SST) has a profound climatic influence. The El Niño‐Southern Oscillation (ENSO) significantly impacts the North Pacific SST; however, the influence of the distinct phases of ENSO on SST predictability remains unclear. To overcome the model limitations, this study assessed SST predictability under diverse ENSO phases using reanalysis. The predictability limit of the North Pacific SST under La Niña (8.4 months) is longer than that under Neutral (5.9 months) and El Niño (5.5 months) conditions, which unveils asymmetry. This asymmetry mirrors contemporary multimodal prediction skills. Error growth dynamics reveal La Niña's robust signal strength with a slow error growth rate, in contrast to El Niño's weaker signal and faster error growth. There exhibits intermediate signal strength and elevated error growth in Neutral condition. Physically, predictability signal strength aligns with SST variability, whereas the error growth rate correlates with atmospheric‐ocean heating anomalies. La Niña, which induces positive heating anomalies, minimizes the impact of atmospheric noise, resulting in lower error growth. The result is beneficial for improving North Pacific SST predictions.

Drought monitoring is crucial for various sectors, and soil moisture data play a pivotal role, especially in agricultural contexts. This study focuses on the recent CYGNSS Level 3 soil moisture data derived from the NASA Cyclone Global Navigation Satellite System (CYGNSS), notable for its wide coverage and rapid revisit times, yet underexplored in drought research. Spanning from 1 January 2018 to 31 December 2022, this research analyzed daily CYGNSS soil moisture data, comparing them with the ERA5, SMAP, and GLDAS-NOAH datasets. It was found that the average and standard deviation (std) of CYGNSS soil moisture exhibited spatial patterns largely similar to other datasets, although some regions showed discrepancies (std differences reached up to 0.05 in some regions). The correlation coefficients and RMSE values between CYGNSS and other datasets depended on climate and land cover types. Four drought indicators from different soil moisture datasets were compared with the improved monthly Standardized Precipitation Evapotranspiration Index (SPEI). The drought indicators based on CYGNSS data demonstrate the capacity to describe drought extent and intensity. The correlation coefficients between certain drought indicators obtained from CYGNSS and SPEI reached 0.27 for drought percentage and 0.16 for drought intensity. Further investigations with selected extreme drought cases revealed that the indicator from CYGNSS data is relatively weak, influenced by the selected regions, times, and drought indicators. The results of this study provide insights into the potential application of CYGNSS soil moisture data in drought monitoring, offering a foundation for future research and practical implementation with current and future improved products.

Numerical seasonal forecasts in Earth science always contain forecast errors that cannot be eliminated by improving the ability of the numerical model. Therefore, correction of model forecast results is required. Analog correction is an effective way to reduce model forecast errors, but the key question is how to locate analogs. In this paper, we updated the local dynamical analog (LDA) algorithm to find analogs and depicted the process of model error correction as the LDA correction scheme. The LDA correction scheme was first applied to correct the operational seasonal forecasts of sea surface temperature (SST) over the period 1982–2018 from the state-of-the-art coupled climate model named NCEP Climate Forecast System, version 2. The results demonstrated that the LDA correction scheme improves forecast skill in many regions as measured by the correlation coefficient and root-mean-square error, especially over the extratropical eastern Pacific and tropical Pacific, where the model has high simulation ability. El Niño–Southern Oscillation (ENSO) as the focused physics process is also improved. The seasonal predictability barrier of ENSO is in remission, and the forecast skill of central Pacific ENSO also increases due to the LDA correction method. The intensity of the ENSO mature phases is improved. Meanwhile, the ensemble forecast results are corrected, which proves the positive influence from this LDA correction scheme on the probability forecast of cold and warm events. Overall, the LDA correction scheme, combining statistical and model dynamical information, is demonstrated to be readily integrable with other advanced operational models and has the capability to improve forecast results.