REVIEW 3 major objections 4 minor 16 references
The ENSO–monsoon link is strong and stationary when measured with a global subsurface index, because the monsoon lags a synchronized chaotic partner by 18 months.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
A D20-based Global-ENSO index shows a stable ~0.8 correlation with Indian monsoon rainfall at 18-month lead via lag synchronization of chaotic oscillators, independent of global warming.
T0 review reviewed 2026-07-13 challenge →
load-bearing objection Solid multi-dataset evidence that a D20-based Global-ENSO index yields a high, relatively stationary 18-month correlation with ISMR; the chaotic lag-synchronization framing is an untested re-labeling of that correlation. the 3 major comments →
Resolving the Paradox of Changing El Ni\~no-Monsoon Relation through Synchronization of Chaotic Oscillators
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
When the ENSO–monsoon relationship is measured with the D20-based Global-ENSO index Dp rather than with Pacific SST, the true relationship is strong (long-term correlation ~0.8–0.85), stationary across the historical record and selected CMIP6 projections, and arises as lag synchronization between ISMR and Dp at an 18-month lead.
What carries the argument
Global-ENSO index Dp: the projection of tropical (30°S–30°N) 20°C isotherm depth anomalies onto the long-term correlation pattern with Indian summer monsoon rainfall. At 18-month lead this index both maximizes the correlation and exhibits elevated phase-locking with the monsoon, which the authors interpret as lag synchronization of chaotic oscillators.
Load-bearing premise
That high lagged correlation and rising phase-locking value prove lag synchronization of two deterministic chaotic oscillators, rather than shared multi-decadal forcing or residual analysis artifacts in the subsurface data.
What would settle it
Show that the 18-month lead Dp–ISMR correlation and the associated phase-locking value collapse once off-equatorial multi-decadal signals (AMV/PDO imprint) are removed from D20, or that independent ocean reanalyses and fully coupled models fail to reproduce a stable high correlation at that lead.
If this is right
- Apparent multi-decadal weakening of the ENSO–monsoon relationship is an artifact of Pacific SST indices and does not imply declining monsoon predictability.
- Long-lead (18-month) skill for Indian summer monsoon rainfall is physically grounded and remains available under moderate greenhouse-gas forcing.
- Analogous Global-ENSO predictors constructed for other tropical monsoon systems should likewise reveal long-lead predictability if the same lag-synchronization mechanism operates.
- Predictability estimates should be taken as the highest correlation at any lead, not only the simultaneous Pacific-SST correlation.
Where Pith is reading between the lines
- If lag synchronization is the operative mechanism, the same construction of a subsurface Global-ENSO index should unlock multi-year predictability for other tropical monsoons and for ENSO itself.
- The result suggests that slow ocean memory can systematically overcome the spring predictability barrier once the correct spatial filter is applied.
- A practical next test is whether deep-learning or dynamical models that ingest tropical D20 at 18-month lead can convert the reported potential skill into real forecast skill over independent decades.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript argues that the well-known multi-decadal fluctuations in the ENSO–Monsoon relationship (EMR) are largely an artifact of using noisy Pacific SST indices (Niño 3.4, Sp). It constructs a Global-ENSO index Dp by projecting tropical D20 anomalies onto the long-term (1875–2010) correlation map with ISMR, shows that Dp–ISMR correlation reaches ~0.8–0.85 at 18-month lead and is nearly stationary across historical data and selected CMIP6 projections, and attributes this stationarity and long-lead skill to lag synchronization of two chaotic oscillators (ISMR and Dp). Supporting evidence includes multi-dataset moving correlations, regression maps of SST and monsoon Hadley circulation, scatter plots, and phase-locking values (PLV rising from 0.44 at lag 0 to 0.75 at lag 18).
Significance. If the empirical stationarity of the Dp–ISMR teleconnection holds, the paper would resolve a long-standing controversy about the reliability of monsoon predictability under climate change and would strengthen the case for subsurface ocean predictors in seasonal-to-interannual forecasting. The multi-dataset comparison of Pacific SST indices and the demonstration that a carefully constructed D20-based index remains skillful at long leads are useful contributions. The claim of a new physical basis for predictability beyond the conventional chaos limit is potentially high-impact, but currently rests on correlative diagnostics rather than a completed dynamical demonstration of chaotic lag synchronization.
major comments (3)
- Section 3.1 and Abstract: The central mechanistic claim that stationarity of EMR is a natural consequence of lag synchronization of deterministic chaotic oscillators is not yet established. The paper itself states that a dynamical demonstration is currently under investigation and that the long-lead filtering of sampling noise is only an intuitive understanding (Sec. 3.2). Linear scatter (Fig. 2j) and elevated PLV are consistent with shared low-frequency forcing (AMV/PDO imprint on off-equatorial D20, also discussed in Sec. 3.2) and do not uniquely diagnose chaotic lag synchronization (Rosenblum et al.). Either supply Lyapunov spectra, mutual-information or controlled-coupling diagnostics that isolate the oscillators from AMV/PDO, or reframe the claim as high lagged correlation consistent with lag synchronization pending further dynamical tests.
- Text S1 and construction of Dp: Dp is obtained by projecting D20 anomalies onto the long-term (1875–2010) correlation map with ISMR itself. Using the same map both to define the predictor and to evaluate its skill introduces a circularity burden that inflates the reported long-term correlation (~0.8). A fully independent construction (e.g., correlation map from a non-overlapping early period, or a leave-one-decade-out scheme) and out-of-sample skill scores are needed before the stationarity and predictability claims can be considered robust.
- Section 3.3 and Figs. 3g–h / S3: The claim that EMR remains strong under future GHG forcing is supported by only a subset of CMIP6 models (at most three under RCP4.5). Inter-model spread is large and several models show weak or insignificant correlations. The text should quantify how many models pass a pre-defined skill threshold, discuss systematic biases in ENSO–monsoon coupling and multidecadal variability, and avoid generalizing from the best-performing models to the broader ensemble.
minor comments (4)
- Figure 1 panels (c–f) and (g–h) are dense; adding panel labels and a clearer legend for the ensemble-mean line would improve readability.
- Notation for lead times is inconsistent (lead 0, −5, −18, 18-month lead). Standardize throughout text and figures.
- Table S1 reports no mean for Kaplan Niño 3.4; either supply the value or note why it is omitted.
- Several recent references on ENSO–monsoon non-stationarity and subsurface predictors could be added for completeness; the existing citation list is otherwise adequate.
Circularity Check
Dp is constructed by projecting D20 onto the long-term ISMR correlation map; the reported high skill and stationarity of EMR are therefore partly by construction of the predictor itself.
specific steps
-
self definitional
[Text S1 (G-ENSO Indices) and Introduction (construction of Dp)]
"To construct Dp, first we obtain correlation maps between D20 and ISMR anomalies over a long time period (1875-2010) and retain correlations with >95% significance (Text S1). ... Next, using Hadamard (location-by-location) multiplication of the spatial pattern of D20 anomalies with the correlation map and then summing up across all locations, we obtain the Dp time series. ... For any given lead, the index is generated by projecting the anomaly field ... onto a long-term correlation pattern obtained between the corresponding anomaly field and ISMR anomalies during a reference period."
Dp is defined as the projection of D20 onto the very correlation pattern that maximises its association with ISMR over the same historical window used for skill evaluation. The subsequent claim that ISMR-Dp correlation reaches ~0.8-0.85 and is stationary is therefore partly tautological: the index is constructed to be the linear combination of D20 that best matches ISMR, so high correlation (and reduced epochal scatter once the map is fixed) is expected by construction rather than independently discovered.
-
fitted input called prediction
[Section 3.2 / Figures 3d-f and 1g-h; Abstract claim of 'strong and stable correlation (~0.8)']
"Contrary to previous findings, ISMR shows a strong and stable correlation (~0.8) with Dp at 18-month lead during the historical period. ... At 18-month lead the EMR stabilizes, showing minimal epochal variability and consistently strong positive correlations across datasets (Figure 3f). The long-term correlation increases from 0.6 at zero lead to 0.85 at 18 months."
The 'prediction' skill and the stationarity of EMR are evaluated with the same Dp that was fitted (via the 1875-2010 correlation map) to maximise association with ISMR. The 18-month peak and the disappearance of multidecadal fluctuations are therefore not free predictions; they inherit the spatial filter that was chosen precisely because it correlates with the target. Moving-window correlations still use the fixed map derived from the full period, so the reported stability is not an out-of-sample discovery.
-
self citation load bearing
[Introduction and Section 3.1 (prior establishment of Dp skill)]
"Having established Dp as a superior ISMR predictor (Sharma, Das, Chakraborty, et al., 2025; Sharma et al., 2022) (Figure 1g), here, we provide a physical basis for the 18-month lead predictability of ISMR. ... deep learning model demonstrates the feasibility of achieving a correlation skill of 0.65 at 18-month lead over 44 years of independent ISMR forecasts (Sharma, Das, Chakraborty, et al., 2025)."
The claim that Dp is already a 'superior' predictor whose 18-month skill needs only a physical interpretation is load-bearing on the authors' own prior papers that introduced the identical construction. Those earlier results are not independently re-derived or externally validated here; they are imported to justify treating the high lagged correlation as given, after which the present paper re-labels it as lag synchronization.
-
renaming known result
[Section 3.1-3.2 (lag synchronization claim)]
"the high correlation and PLV at longer leads signifies that these two chaotic systems are lag synchronized rather than being independent. ... We discover that the stationarity in the El Nino-Monsoon relationship emerges as a natural consequence of chaotic synchronization between Indian summer monsoon rainfall and 20 degree Celsius isotherm at an 18-month lead. ... Although the implications of chaotic synchronization between ISMR and Dp at 18-month lead ... is currently under investigation"
The observables offered (linear scatter Y(t)≈X(t-τ) and PLV rising from 0.44 to 0.75) are simply re-descriptions of the high lagged correlation already produced by the self-defined Dp. The paper renames this correlation 'lag synchronization of chaotic oscillators' while conceding that a dynamical demonstration remains under investigation and that the multidecadal memory itself arises from shared AMV/PDO forcing on off-equatorial D20. No new diagnostic (Lyapunov spectra, mutual information, controlled coupling) is supplied; the known empirical pattern is re-labelled.
full rationale
The paper's central empirical claim (stable EMR ~0.8 at 18-month lead, independent of warming) rests on the G-ENSO index Dp. Text S1 and the Introduction define Dp by projecting D20 anomalies onto the long-term (1875-2010) correlation map between D20 and ISMR itself, retaining only >95% significant locations. The same map that defines the predictor is then used to evaluate its correlation with ISMR (Figs. 1g-h, 3d-f). This is a classic self-definitional / fitted-input construction: the index is the spatial pattern that maximises contemporaneous association with the target, so elevated lagged skill and reduced epochal variability are partly forced by that choice rather than discovered independently. The lag-synchronization interpretation (Y(t)≈X(t-τ), PLV rising to 0.75) is not algebraically forced by the same map, but the paper itself states that a dynamical demonstration is 'currently under investigation' and that the long-lead filtering of sampling noise is only an 'intuitive understanding'; the synchronization claim therefore inherits the circular construction of Dp without independent verification (Lyapunov spectra, controlled coupling, isolation from AMV/PDO). Self-citations to Sharma et al. (2022, 2025) supply the prior skill numbers but are not the sole load-bearing step; the definitional circularity of Dp is. Score 6 reflects partial circularity of the main empirical result without rendering the entire paper vacuous.
Axiom & Free-Parameter Ledger
free parameters (3)
- 95 % significance threshold for correlation-map mask
- 21-year / 31-year moving-correlation window lengths
- 18-month lead selected as the synchronization lag
axioms (4)
- ad hoc to paper High lagged linear correlation plus elevated phase-locking value between two climate time series implies lag synchronization of deterministic chaotic oscillators.
- domain assumption Subsurface D20 anomalies are essentially free of the high-frequency climate noise that contaminates SST, and therefore provide a cleaner memory of ENSO precursors.
- domain assumption Off-equatorial tropical D20 (10–30° latitude) carries a quasi-stationary multi-decadal imprint of AMV/PDO that phase-locks with ISMR multi-decadal variability at long leads.
- domain assumption Linear projection of anomaly fields onto a fixed long-term correlation map yields a physically meaningful Global-ENSO index.
invented entities (2)
-
Global-ENSO (G-ENSO) indices Sp, Dp, Hp
no independent evidence
-
Lag-synchronized chaotic oscillator pair (ISMR, Dp at 18-month lead)
no independent evidence
Cite this review
Pith. "Pith review of Resolving the Paradox of Changing El Ni\~no-Monsoon Relation through Synchronization of Chaotic Oscillators." pith.science (2026). https://pith.science/paper/AHEMI5ZU
@misc{pith2026260316346,
author = {Pith},
title = {Pith review of: Resolving the Paradox of Changing El Ni\~no-Monsoon Relation through Synchronization of Chaotic Oscillators},
year = {2026},
howpublished = {\url{https://pith.science/paper/AHEMI5ZU}},
note = {Machine review of arXiv:2603.16346}
}
read the original abstract
For over a century, the relationship between Indian summer monsoon rainfall and El Nino-Southern Oscillation has been the foundation of 'long-range' prediction of Indian monsoon. This relation is estimated from correlations between Indian summer monsoon rainfall and a Pacific sea surface temperature-based index of El Nino-Southern Oscillation. However, a prominent multi-decadal variability in the correlation raises doubts on the realism of El Nino-Monsoon relation and stability of Indian monsoon predictability. Previous studies discussed that Pacific-based El Nino-Southern Oscillation indices do not represent El Nino's global influence completely, making their correlation with Indian monsoon unreliable. To address this limitation, a Global El Nino-Southern Oscillation framework based on the depth of the 20 degree Celsius isotherm is developed, integrating subsurface signal from all three tropical ocean basins and maximizing Indian monsoon teleconnections. Contrary to previous findings, the 20 degree Celsius isotherm-based Global El Nino-Southern Oscillation exhibits a strong and stable correlation (greater than 0.8) with Indian monsoon at an 18-month lead. Through a re-examination of the El Nino-Monsoon relationship with the superior Global El Nino-Southern Oscillation predictor, we show that the true relationship is independent of global warming and stationary in time. We discover that the stationarity in the El Nino-Monsoon relationship emerges as a natural consequence of chaotic synchronization between Indian summer monsoon rainfall and 20 degree Celsius isotherm at an 18-month lead. Such synchronization between chaotic climate variables provides a new physical basis for climate predictability beyond the conventional deterministic limit set by chaos. Our findings provide the foundation and renewed confidence in long range climate prediction.
Reference graph
Works this paper leans on
-
[1]
1 Resolving the Paradox of Changing ENSO-Monsoon Relation through Global-ENSO Devabrat Sharma1,2, Shruti Tandon1,2, Gaurav Chopra3, R. I. Sujith1,2, and B. N. Goswami4,* 1Department of Aerospace Engineering, Indian Institute of Technology Madras, Chennai- 600036, India 2Centre of Excellence for Studying Critical Transitions in Complex Systems, Indian Inst...
2000
-
[2]
While ENSO is considered the dominant predictor of ISMR (Gowariker et al., 1989; Rajeevan, 2001; Rajeevan et al., 2004; Rajeevan & Pai, 2007; Thapliyal & Kulshrestha S
(Figure 1c–f). While ENSO is considered the dominant predictor of ISMR (Gowariker et al., 1989; Rajeevan, 2001; Rajeevan et al., 2004; Rajeevan & Pai, 2007; Thapliyal & Kulshrestha S. M., 1992), the reported post-1980s (K. K. Kumar et al., 1999; Xavier et al., 2007), weakening of EMR has raised concerns about ENSO’s predictive role. Several studies have a...
1989
-
[3]
Kumar et al., 2007; X
appears weaker (Figure 1c-f), remaining within the multidecadal variability range, and even suggests modest strengthening post-2000s (P. Kumar et al., 2007; X. Yang & Huang, 2021). Furthermore, when EMR is evaluated using non-traditional predictor indices of ISMR, such as the thermodynamic Indian Summer Monsoon index (TISM), the results show weak multidec...
2007
-
[4]
First, ISMR predictability is always underestimated due to ‘climate noise’ in Niño 3.4, making it unsuitable for long-lead predictions
identified two major issues in using Niño 3.4 as a predictor to assess EMR. First, ISMR predictability is always underestimated due to ‘climate noise’ in Niño 3.4, making it unsuitable for long-lead predictions. As a result, the ENSO-ISMR teleconnection strength at longer forecast lead times cannot be estimated using Niño 3.4. Second, the influence of ENS...
2022
-
[5]
Further, we show that not only is ISMR more strongly coupled to G-ENSO than to Pacific ENSO, but also that EMR is not changing as perceived earlier
(Figure 1g), here, we provide a physical basis for the 18-month lead predictability of ISMR. Further, we show that not only is ISMR more strongly coupled to G-ENSO than to Pacific ENSO, but also that EMR is not changing as perceived earlier. Remarkably, we find that EMR was strong in the past and likely to remain strong in the foreseeable future, even und...
2016
-
[6]
7 The Niño 3.4 index, a measure of estimating the strength of the Pacific ENSO is measured using SST anomalies averaged over 170°W-120°W, 5°S-5°N
which is based on a fixed network of 1384 rain gauge stations, and (iv) the IMD Climate Monitoring Portal (referred to as IMD2), also covering 1901–2023 (IMD, 2023). 7 The Niño 3.4 index, a measure of estimating the strength of the Pacific ENSO is measured using SST anomalies averaged over 170°W-120°W, 5°S-5°N. We show the time series of four Niño 3.4 ind...
1901
-
[7]
In this backdrop, high potential predictability at 18-month lead may seem counterintuitive
is an evidence of influence of chaos even on the seasonal mean. In this backdrop, high potential predictability at 18-month lead may seem counterintuitive. However, studies in nonlinear dynamics have shown that when two chaotic oscillators are coupled, their dynamics can synchronize such that one oscillator follows a time-lagged version of the other, a re...
1997
-
[8]
(g-j) Scatter plot of ISMR with Sp(0), Dp(0), Dp(-5), and Dp(-18), respectively
(f) Same as (e) but for ISMR from IITM data and Dp at 5- and 18-month lead. (g-j) Scatter plot of ISMR with Sp(0), Dp(0), Dp(-5), and Dp(-18), respectively. 3.2 ENSO-Monsoon Relationship and ISMR predictability in a Warming Climate The high ISMR-Dp correlation at 18-month lead is a result of lag synchronization (Pikovsky et al., 2010; Rosenblum et al.,
2010
-
[9]
Lag synchronization occurs when the states of two dynamical systems become nearly identical but are separated by a finite time delay, such that Y(t) ≈ X(t-�)
of interannual variability. Lag synchronization occurs when the states of two dynamical systems become nearly identical but are separated by a finite time delay, such that Y(t) ≈ X(t-�). Thus, two key signatures are expected: (i) a delayed correspondence between the states of the systems, and (ii) a consistent phase relation between their oscillatory comp...
2020
-
[10]
Figure 3: G-ENSO-ISMR relationship in a warming climate
can generate quasi- 14 stationary off-equatorial oceanic Rossby waves, introducing multidecadal variability in the off-equatorial thermocline, which is captured by Dp at 18-month lead. Figure 3: G-ENSO-ISMR relationship in a warming climate. (a) 21-year moving simultaneous correlation between JJAS Sp and ISMR from (i) IITM, (ii) ERA5, (iii) IMD, and (iv) ...
2023
-
[11]
on the ENSO-ISMR relationship, it has been estimated using equatorial Pacific SST indices, treating them as ideal predictors. Regional SST-based predictors such as the IOD, PDO, and Atlantic Niño, often assumed independent of ENSO, have also been used as ISMR predictors at various leads. Our study show that these assumptions are premature, as regional SST...
-
[12]
https://doi.org/10.1175/JCLI-D-19-0395.1 Kinter, J
Journal of Climate, 33(4), 1351–1379. https://doi.org/10.1175/JCLI-D-19-0395.1 Kinter, J. L., Miyakoda, K., & Yang, S. (2002). Recent Change in the Connection from the Asian Monsoon to ENSO. Journal of Climate, 15, 1203–1215. Retrieved from http://www.cpc.ncep.noaa.gov/ IMD. (2023). Southwest Monsoon Rainfall Data (1901 – 2023). Climate Monitoring Portal,...
-
[13]
Current Sciences, 86, 422–431. Retrieved from https://www.tropmet.res.in/~kolli/MOL/Monsoon/year2005/imdlrfmodel.pdf Rajeevan, M., Bhate, J., & Jaswal, A. K. (2008). Analysis of variability and trends of extreme rainfall events over India using 104 years of gridded daily rainfall data. Geophysical Research Letters, 35(18), L18707. https://doi.org/https://...
-
[14]
https://doi.org/10.1029/2002jd002670 Rosenblum, M. G., Pikovsky, A. S., & Kurths, J. (1996). Phase Synchronization of Chaotic Oscillators. Physical Review Letters, 76(11), 1804–1807. Sahastrabuddhe, R., Ghausi, S. A., Joseph, J., & Ghosh, S. (2023). Indian Summer Monsoon Rainfall in a changing climate: a review. Journal of Water and Climate Change, 14(4),...
-
[15]
The conclusions qualitatively remain unchanged if we use different SST and rainfall data sets
and SST from COBE (Hirahara et al., 2014). The conclusions qualitatively remain unchanged if we use different SST and rainfall data sets. Text S2: Phase Locking Value The Phase Locking Value (Lachaux et al.,
2014
-
[16]
The dark contours represent correlation above 95% confidence level. � 28 Figure S2: Power spectrum of 15-month running D20 anomaly averaged over the tropical region between (a) 5°S-5°N, (b) 20°S-20°N. ���������������������������������� ������������������ � 29 Figure S3: Correlation between Dp and ISMR from CMIP6 model simulations with (a) RCP2.6 and (b) R...
This paper was first reviewed by grok-4.5 on July 13, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.