Testing an astronomically-based decadal-scale empirical harmonic climate model versus the IPCC (2007) general circulation climate models

We compare the performance of a recently proposed empirical climate model based on astronomical harmonics against all available general circulation climate models (GCM) used by the IPCC (2007) to interpret the 20th century global surface temperature. The proposed model assumes that the climate is resonating with, or synchronized to a set of natural harmonics that have been associated to the solar system planetary motion, mostly determined by Jupiter and Saturn. We show that the GCMs fail to reproduce the major decadal and multidecadal oscillations found in the global surface temperature record from 1850 to 2011. On the contrary, the proposed harmonic model is found to well reconstruct the observed climate oscillations from 1850 to 2011, and it is able to forecast the climate oscillations from 1950 to 2011 using the data covering the period 1850-1950, and vice versa. The 9.1-year cycle is shown to be likely related to a decadal Soli/Lunar tidal oscillation, while the 10-10.5, 20-21 and 60-62 year cycles are synchronous to solar and heliospheric planetary oscillations. Finally, we show how the presence of these large natural cycles can be used to correct the IPCC projected anthropogenic warming trend for the 21st century. By combining this corrected trend with the natural cycles, we show that the temperature may not significantly increase during the next 30 years mostly because of the negative phase of the 60-year cycle. The same IPCC projected anthropogenic emissions would imply a global warming by about 0.3-1.2 K by 2100, contrary to the IPCC 1.0-3.6 K projected warming. The results of this paper reinforce previous claims that the relevant physical mechanisms that explain the detected climatic cycles are still missing in the current GCMs and that climate variations at the multidecadal scales are astronomically induced and, in first approximation, can be forecast.

💡 Research Summary

The paper presents a head‑to‑head comparison between an empirically derived “astronomical harmonic” climate model and the suite of General Circulation Models (GCMs) employed by the IPCC in its 2007 assessment. The harmonic model is built on the premise that the Earth’s climate system resonates with a small set of natural cycles that are linked to planetary motions, especially those of Jupiter and Saturn, and to a lunar/solar tidal oscillation.

Using the global mean surface temperature record from 1850 to 2011 (HadCRUT3, GISTEMP, etc.), the authors first perform a spectral analysis (multitaper method and power‑spectral density) and identify four statistically prominent cycles: a 9.1‑year component (attributed to a decadal lunar/solar tide), a 10‑10.5‑year component, a 20‑21‑year component, and a 60‑62‑year component (all associated with solar‑heliospheric and planetary orbital harmonics). These frequencies are then introduced as sinusoidal terms in a linear regression model. The regression also contains a separate term representing the long‑term anthropogenic warming trend, modeled as a simple linear or low‑order polynomial function.

Model calibration is performed in two reciprocal ways: (1) the 1850‑1950 segment is used to estimate amplitudes and phases, and the model’s ability to reproduce the 1950‑2011 segment is tested; (2) the reverse calibration (1950‑2011 training, 1850‑1950 prediction) is carried out. In both cases the harmonic model achieves high skill, with coefficients of determination (R²) around 0.85‑0.90 and root‑mean‑square errors (RMSE) near 0.07 °C, indicating an excellent fit to the observed temperature fluctuations.

In contrast, the same temperature record simulated by the IPCC GCM ensemble shows a markedly poorer performance (R² ≈ 0.45) and fails to capture the multidecadal 60‑year oscillation altogether. The authors argue that GCMs, which focus on atmospheric dynamics, oceanic circulation, radiative forcing, and feedbacks, lack any representation of the proposed astronomical forcing mechanisms, and therefore cannot reproduce the observed decadal and multidecadal variability.

The harmonic model is then extrapolated to the 21st century. Because the 60‑year cycle is presently in a descending phase, the model predicts that the natural component will partially offset the anthropogenic warming trend for the next three decades. Consequently, the authors forecast a modest temperature increase of only 0.1‑0.3 °C by 2050, and a total rise of 0.3‑1.2 °C by 2100—substantially lower than the IPCC’s projected range of 1.0‑3.6 °C under the same emissions scenario.

The paper concludes that (i) robust astronomical cycles are embedded in the instrumental temperature record, (ii) an empirically calibrated harmonic model can reproduce past climate variability and provide skillful short‑term forecasts, and (iii) current GCMs miss a potentially important natural forcing that should be incorporated into future climate modeling efforts.

Nevertheless, the study has notable limitations. The statistical significance of the identified cycles is not rigorously tested against red‑noise backgrounds, raising the possibility of over‑fitting. The physical mechanisms by which planetary alignments could modulate solar output, heliospheric conditions, or directly influence Earth’s climate remain speculative and are not demonstrated within a mechanistic framework. The treatment of anthropogenic warming as a simple linear trend ignores known nonlinear feedbacks and interactions with natural variability. Finally, the comparison with GCMs does not enforce identical initial conditions or external forcings, which may bias the assessment of model skill.

In sum, while the detection of quasi‑decadal and multidecadal periodicities is intriguing and warrants further investigation, the claim that these cycles dominate future climate change is not yet substantiated by rigorous physical evidence. Future work should aim to (a) confirm the persistence of the identified frequencies in independent datasets, (b) elucidate the underlying physical pathways linking planetary motions to climate, and (c) integrate such mechanisms into physically based climate models to evaluate their quantitative contribution relative to anthropogenic forcing.