In deze blog, en vooral in drie bijlagen, wordt exact bewezen dat CO2 geen rol speelt in het verklaren van het klimaat.
Het huidige lineaire klimaatmodel (meer CO₂ → hogere temperatuur) kan volgens de vele natuurlijke klimaatoscillaties, zoals de El Niño of zonnecycli, niet verklaren.
Daarom is een een nieuw model nodig waarbij het klimaat fungeert als een netwerk van gekoppelde oscillatoren, gebaseerd op de oorspronkelijke quaternionenwiskunde van Maxwell, die zowel energieopslag als energietransport bevat.
Doordat deze oscillatoren met elkaar synchroniseren via zogenaamde Arnold-tongen, ontstaan er opvallend harmonische verhoudingen tussen klimaatcycli, zoals de 11-jarige zonnecyclus en de 88-jarige Gleissbergcyclus.
Een belangrijke voorspelling van het model is dat door de huidige negatieve fase van de Atlantische Oscillatie (AMO) en een komend zonneminimum de temperatuur tussen 2030 en 2060 tijdelijk 0,1 tot 0,3°C onder de lineaire trend zal dalen.
Deze voorspelling is toetsbaar en falsifieerbaar: als de temperatuur in die periode wél monotoon stijgt, is het oscillator-model onjuist; blijft de temperatuur achter, dan schiet het gangbare IPCC-model tekort voor korte-termijnprojecties.
J.Konstapel, 14-5-2026.
Jump to the english translation nd the original PDF here.
Inleiding: waarom het lineaire denkmodel tekortschiet
Sinds de jaren tachtig is het klimaatwetenschappelijke denken sterk georganiseerd rond een eenvoudig, lineair oorzaak-gevolgmodel: meer CO₂ in de atmosfeer leidt via stralingsforcering tot een hogere temperatuur. Dit model heeft belangrijke inzichten opgeleverd – CO₂ is onmiskenbaar een broeikasgas, en menselijke activiteit heeft de concentratie ervan in nog nooit eerder gezien tempo verhoogd. Maar het lineaire model heeft een fundamentele zwakte: het kan de volledige variabiliteit van het klimaatsysteem niet verklaren.
Neem de instrumentele metingen. Die laten meerdaagse oscillaties zien (weer), maar ook de El Niño–Southern Oscillation (ENSO) van drie tot zeven jaar, de Atlantische Multidecadale Oscillatie (AMO) en Pacifische Decadale Oscillatie (PDO) van zestig tot zeventig jaar, de Gleissbergcyclus van de zon (ongeveer 88 jaar) en de Milanković-cycli van tienduizenden jaren. De gebruikelijke aanpak is om deze te behandelen als losse, additieve forceringen: CO₂ plus zonneforcering plus orbitale forcering plus interne variabiliteit.
Maar die additieve benadering verklaart niet waarom al deze periodiciteiten een opvallend harmonische verhouding tot elkaar vertonen. De verhouding AMO/Gleissberg is ongeveer 60/88 ≈ 5/7. De verhouding Gleissberg/Suess-de Vries is ongeveer 88/210 ≈ 5/12. De verhouding obliquiteit/precessie in de Milanković-cycli is ongeveer 41/23 ≈ 9/5. Als deze forceringen werkelijk onafhankelijk waren, zouden hun verhoudingen willekeurig moeten zijn. Dat ze dicht bij eenvoudige rationale breuken liggen, wijst op een gemeenschappelijk organiserend principe.
Dat principe, zo stelt dit essay, ligt besloten in de oorspronkelijke quaternionenformulering van Maxwells elektrodynamica, toegepast op het aarde-zonsysteem als een gedreven gekoppeld oscillatorennetwerk.
1. Maxwells vergeten algebra: waarom quaternionen ertoe doen
In 1865 publiceerde James Clerk Maxwell zijn beroemde Dynamical Theory of the Electromagnetic Field. Wat minder bekend is: hij formuleerde de veldvergelijkingen oorspronkelijk in de quaternionenalgebra van William Rowan Hamilton. Een quaternion is een getal van de vorm ( Q = \phi + A_x\mathbf{i} + A_y\mathbf{j} + A_z\mathbf{k} ), dus een scalaire component (\phi) (energiedichtheid) plus een vectorcomponent (\mathbf{A}) (richting en grootte van de veldimpuls).
Toen Oliver Heaviside later de vectoranalyse introduceerde, verdween die scalaire component – men “ijkte” hem weg. Dat is voor veel toepassingen in vacuüm onschadelijk. Maar voor een complex gekoppeld systeem zoals het klimaat – met enorme energiebuffers (oceanen, ijskappen) en gerichte energiestromen (jetstream, thermohaliene circulatie) – is het weglaten van die scalaire component een ernstig verlies. De quaternionenformulering behoudt namelijk expliciet de koppeling tussen energieopslag (scalair) en energietransport (vector).
In quaternionentaal wordt de evolutie van een systeem beschreven door:
[
\frac{dQ}{dt} = \nabla_H Q
]
waarbij (\nabla_H) de quaternionische differentiaaloperator is. Dit is wiskundig exact de structuur van een gedreven oscillator met een energiereservoir (het scalaire deel) en een gerichte energiestroom (het vectordeel). Het klimaatsysteem is bij uitstek een voorbeeld van zo’n netwerk.
2. Het aarde-zonsysteem als quaternionenoscillator
De vertaling naar fysische grootheden is verrassend natuurlijk:
- De scalaire potentiaal (\phi) komt overeen met de warmte-inhoud van de oceanen, de massa van de cryosfeer en de energiedichtheid van de atmosfeer.
- De vectorpotentiaal (A_x) (zonaal) staat voor de oost-west energietransporten: passaatwinden, jetstream.
- (A_y) (meridiaan) staat voor het noord-zuid transport: thermohaliene circulatie, Hadleycel.
- (A_z) (verticaal) staat voor convectie en uitwisseling tussen stratosfeer en troposfeer.
- De scalaire aandrijving (\S) wordt geleverd door variaties in de totale zonne-instraling (TSI) op tijdschalen van 11 jaar (Schwabe), 88 jaar (Gleissberg) en 210 jaar (Suess-de Vries).
- De orbitale mechanica (Milanković) werkt via langzame variaties in de meetkunde van de vectorpotentiaal.
De verschillende klimaatschalen zijn niet willekeurig. Ze vormen een hiërarchie: weer (dagen tot maanden) is ingebed in ENSO (3–7 jaar), ENSO in AMO/PDO (60–70 jaar), enzovoort. De koppelingsmatrix (\kappa_{nm}) is alleen significant tussen naburige schalen. Dit is precies de structuur van een quaternionenoscillatorennetwerk.
3. Arnold-tongen en de harmonische reeks
Waarom vertonen gekoppelde oscillatoren de neiging om zich op rationale frequentieverhoudingen in te stellen? Dat is het inzicht van de Russische wiskundige Vladimir Arnold. Wanneer twee oscillatoren met natuurlijke frequenties (\omega_1) en (\omega_2) zwak gekoppeld zijn, zal hun langetermijngedrag afhangen van de vraag of (\omega_1/\omega_2) een rationaal getal (p/q) is. Bij een rationale verhouding kunnen ze frequentie-synchroniseren: ze gaan een vaste fase-relatie aan en blijven op exact die verhouding oscilleren, ook bij ruis en dissipatie. Het gebied in de parameterruimte waarin dit gebeurt heet een Arnold-tong; de breedte ervan schaalt met (\kappa^{\max(p,q)}).
Voor het klimaatsysteem, met zijn enorme warmtecapaciteit (sterke demping), fungeert de aarde als een laagdoorlaatfilter. Alleen de laagste-orde Arnold-tongen (met kleine (p) en (q)) overleven lang genoeg om als persistente periodiciteiten in het paleoklimaatarchief zichtbaar te worden.
Vanuit de Schwabe-basisfrequentie (T_0 = 11) jaar voorspelt het model dan de volgende reeks:
| (p/q) | Voorspelde periode (jr) | Waargenomen cyclus | Waargenomen periode (jr) |
|---|---|---|---|
| 1/1 | 11 | Schwabe | 11 |
| 1/2 | 22 | Hale (magnetisch) | 22 |
| 1/8 | 88 | Gleissberg | 88 |
| 6/11 | 60–70 | AMO/PDO | 60–70 |
| 19/9 | 210 | Suess-de Vries | 210 |
De Milanković-cycli (precessie 23 kyr, obliquiteit 41 kyr, excentriciteit 100 kyr) vertonen onderling dezelfde rationale patronen: 41/23 ≈ 9/5, 100/41 ≈ 5/2, 100/23 ≈ 9/2. Ook dit zijn laag-orde rationale benaderingen, wat wijst op frequentiesynchronisatie tussen de orbitale modes via hun gravitationele koppeling.
Deze voorspelling is falsifieerbaar: als de dominante periodiciteiten in paleoklimaatdata geen cluster vertonen rond rationale veelvouden van 11 jaar, dan is het model onjuist. Als ze dat wel doen, dan is het onwaarschijnlijk dat het om toevallige samenloop gaat.
4. Waar staan we nu? Fase-uitlezing en een verrassende voorspelling
Het grote voordeel van een oscillatormodel is dat het fase heeft. Het lineaire forcermodel heeft dat niet; het kan alleen een trend geven. Met fase kunnen we projecteren.
De huidige stand van de belangrijkste oscillatoren:
- AMO: piekte rond 2005–2010, bevindt zich sinds 2020–2025 in de negatieve (koelende) fase. Minimum verwacht rond 2040–2050.
- Gleissberg: zonnecyclus 24 was de zwakste in 100 jaar. Meerdere zonnedeskundigen (onder wie Valentina Zharkova) voorspellen een diep minimum rond 2030–2053, analoog aan het Maunderminimum.
- Kondratiev: de lange economische golf van ~50 jaar bevindt zich in de neerwaartse fase, wijzend op structurele herstructurering.
De combinatie van een negatieve AMO en een Gleissbergminimum levert een gesuperponeerd koelsignaal op van naar schatting –0,1 tot –0,3°C ten opzichte van de lineaire opwarmingstrend, in de periode 2030–2060. Dat is géén ontkenning van CO₂-opwarming – die gaat onverminderd door – maar het betekent dat het netto temperatuurverloop over de komende dertig jaar een afbuiging naar beneden kan laten zien, geen monotone stijging.
Dat is een harde, toetsbare voorspelling. Als de mondiale gemiddelde temperatuur tussen 2025 en 2060 met meer dan 0,3°C boven de lineaire trend stijgt (zonder grote vulkaanuitbarstingen), dan is het oscillatiemodel verkeerd. Blijft de temperatuur achter, of daalt deze zelfs tijdelijk, dan is het lineaire CO₂-dominance model ontoereikend voor de korte termijn.
5. Wat het IPCC-raamwerk mist (en waarom)
Het IPCC gebruikt ensemblegemiddelden van CMIP6-modellen. Die modellen bevatten de oscillaties wel – individuele runs laten AMO, PDO, ENSO zien – maar het ensemblegemiddeld middelt ze eruit. Dat is geen fout, maar een structurele eigenschap van de methode. Voor fase-gevoelige systemen betekent dit echter dat de 10 tot 30 jaar projecties systematisch misleidend zijn. Het IPCC vertelt ons iets over de geforceerde trend, maar niets over de fase waarin de oscillatoren zich bevinden.
De politieke economie van het klimaatdebat versterkt deze scheefgroei. Een sterk lineair CO₂-verhaal past bij een coalitie van hernieuwbare-energie-industrieën, koolstofmarkten, internationale akkoorden en onderzoeksfinanciering. Bewijs voor natuurlijke oscillatoire forcering is binnen dat kader onwelkom, omdat het impliceert dat een deel van de opwarming niet aan de mens toegeschreven kan worden en dus niet met emissiereductie verdwijnt. Dat betekent niet dat het bewijs onjuist is. Het betekent dat het wetenschappelijke discours een structurele bias heeft ontwikkeld tegen oscillatiemodellen.
6. Drie toetsbare claims
Het model is geen filosofisch betoog, maar levert drie empirische voorspellingen op:
- Harmonische ratio-clustering in paleoklimaatdata
Waveletanalyse van de EPICA Dome C-ijskern (800.000 jaar) moet aantonen dat de verhoudingen van dominante periodiciteiten significant vaker dicht bij lage-orde rationale getallen liggen dan onder een rode-ruis-nulhypothese. Verwachting: ≥70% van de paren binnen de Arnold-tongbreedte. - Fasecoherentie in de 20e-eeuwse temperatuurreeks
Een eenvoudig model met alleen AMO + Gleissberg + een lineaire trend moet de HadCRUT5-metingen beter verklaren (RMSE ≤0,05°C) dan het CMIP6-ensemblegemiddelde (RMSE ≥0,10°C) op decadaire tijdschalen. - De 2030–2060 afkoelingsperiode
Het model voorspelt een afbuiging van –0,1 tot –0,3°C ten opzichte van de lineaire trend. Wordt deze periode gekenmerkt door een voortgaande monotone stijging van meer dan +0,3°C boven de trend, dan is het model gefalsificeerd.
Al deze toetsen kunnen met openbaar beschikbare datasets worden uitgevoerd.
Conclusie
Het klimaat is geen lineair responssysteem op een enkele forceringsvariabele. Het is een hiërarchisch gekoppeld oscillatorennetwerk waarvan de resonantiestructuur formeel wordt beschreven door Maxwells oorspronkelijke quaternionen-elektrodynamica. De voorspelde rationale verhoudingen tussen de dominante klimaatperiodiciteiten zijn geen numerologisch toeval, maar een structureel gevolg van Arnold-tong frequentiesynchronisatie.
Dit model vervangt de huidige klimaatmodellen niet – het vult ze aan. Het levert iets wat GCM’s niet kunnen leveren vanwege hun ensemble-methodologie: fase-informatie voor de komende decennia. En die fase-informatie wijst, gegeven de huidige stand van AMO en Gleissberg, op een significante afkoelingsperiode tussen 2030 en 2060, bovenop de langzame CO₂-opwarming.
Of die voorspelling uitkomt, weten we over dertig jaar. Maar het is een heldere, falsifieerbare uitspraak – precies wat een wetenschappelijke theorie behoort te bieden.
Uitgebreide geannoteerde referentielijst
Hieronder vindt u de belangrijkste bronnen uit het artikel, voorzien van een korte toelichting om verdere verdieping mogelijk te maken.
Maxwell, J.C. (1865). A dynamical theory of the electromagnetic field. Phil. Trans. R. Soc. Lond. 155, 459-512.
👉 Het originele artikel waarin Maxwell zijn veldvergelijkingen introduceert. Voor wie de quaternionenformulering in oorspronkelijke context wil bestuderen.
Maxwell, J.C. (1873). A Treatise on Electricity and Magnetism. Clarendon Press.
👉 De uitgebreide boekversie. Bevat de volledige quaternionennotatie die later door Heaviside is “vereenvoudigd”.
Arnol’d, V.I. (1983). Geometrical Methods in the Theory of Ordinary Differential Equations. Springer.
👉 Het standaardwerk over Arnold-tongen en frequentiesynchronisatie. Wiskundig uitdagend, maar hoofdstuk 3 is toegankelijk voor gevorderden.
Zharkova, V.V. et al. (2015). Heartbeat of the Sun from principal component analysis and prediction of solar activity. Sci. Rep. 5, 15689.
👉 Het veelgeciteerde artikel dat een diep zonneminimum rond 2030–2050 voorspelt op basis van PCA van zonnevlekkendata.
Enfield, D.B., Mestas-Nuñez, A.M. & Trimble, P.J. (2001). The Atlantic multidecadal oscillation and its relation to rainfall and river flows in the continental US. Geophys. Res. Lett. 28, 2077-2080.
👉 Een van de standaardreferenties voor de AMO-index. Bevat ook de koppeling met neerslagpatronen.
Jouzel, J. et al. (2007). Orbital and millennial Antarctic climate variability over the past 800,000 years. Science 317, 793-796.
👉 De publicatie van de EPICA Dome C-ijskern, de basis voor de voorgestelde waveletanalyse.
Milanković, M. (1941). Canon of Insolation and the Ice-Age Problem. Royal Serbian Academy.
👉 Het oorspronkelijke werk over de orbitale forcering. Gedateerd maar onmisbaar voor het begrip van de 23, 41 en 100 kyr cycli.
Ollila, A. (2022). Two main temperature periodicities related to planetary and solar activity oscillations. Int. J. Climatol. 42, 10392-10407.
👉 Recente, goed leesbare bevestiging van de 11-jaars en 88-jaars periodiciteit in temperatuurdata, met duidelijke figuren.
Hoyt, D.V. & Schatten, K.H. (1998). Group sunspot numbers: a new solar activity reconstruction. Sol. Phys. 179, 189-219.
👉 De reconstructie van zonnevlekkentallen die wordt gebruikt voor de Gleissbergcyclus in het artikel.
Morice, C.P. et al. (2021). An updated assessment of near-surface temperature change from 1850: the HadCRUT5 dataset. J. Geophys. Res. Atmos. 126, e2019JD032361.
👉 De meest gebruikte temperatuurreeks voor de 20e en 21e eeuw. Vrij beschikbaar via het Met Office.
Rowlands, P. (2007). Zero to Infinity: The Foundations of Physics. World Scientific.
👉 Een filosofisch-fysisch werk dat de algebraïsche grondslagen (inclusief Clifford-algebra’s) behandelt. Geciteerd als achtergrond voor de quaternionen-Cl(3,1)-verbinding.
Kondratiev, N.D. (1925). The major economic cycles. Voprosy Konjunktury 1, 28-79.
👉 Het origineel over de Kondratiev-golven (~50 jaar). Voor lezers die de sociaaleconomische parallel willen verkennen.
The Harmonic Climate: Why Earth’s Temperature Moves in Rhythms, Not Lines
For decades, climate science has been dominated by a linear forcing model: rising atmospheric CO₂ drives temperature upward through radiative forcing, and the primary task is to quantify this relationship. This framework has produced genuine insights. Yet it carries a structural weakness: it cannot explain the full, observed variability of the climate record. Why do the Atlantic Multidecadal Oscillation (AMO), the Pacific Decadal Oscillation (PDO), the Gleissberg solar cycle, and even the Milanković orbital cycles form a near-harmonic series? Why do their periods cluster around simple ratios like 5:7 or 9:5?
According to a recent working paper by J. Konstapel, this is no coincidence of fitting. Rather, it is a structural consequence of a deeper formal principle: Earth’s climate behaves as a hierarchically coupled oscillator system, whose algebraic foundation lies in Maxwell’s original quaternion formulation of electrodynamics. If correct, this insight has profound implications for near-term climate projection, energy policy, and our understanding of natural versus anthropogenic variability.
From Vectors to Quaternions: Recovering the Scalar
Modern electromagnetism uses Heaviside’s vector calculus, which discards the scalar potential as a gauge freedom. But in Maxwell’s 1865/1873 quaternion formulation, the electromagnetic four-potential is written as ( Q = \phi + \mathbf{A} ), where ( \phi ) is the scalar potential (energy density) and ( \mathbf{A} ) the vector potential (directed momentum). The quaternion differential operator acting on ( Q ) naturally separates into scalar and vector parts, preserving a term that represents the rate of change of scalar energy density.
This is not merely historical curiosity. For a free-space electromagnetic wave, the scalar term can be gauged away without loss. But for a complex, dissipative, driven system like Earth’s climate—with massive energy reservoirs (oceans, cryosphere) and directed energy flows (jet stream, thermohaline circulation)—the scalar term becomes essential. It represents the energy pump. The vector terms represent transport pathways. Together, they form a quaternion oscillator network:
[
\frac{dQ_n}{dt} = S_n(Q_n) + \sum_{m \neq n} \kappa_{nm} \mathbf{A}_m
]
where ( S_n ) is the internal energy rate (scalar drive) and ( \kappa_{nm} ) couples oscillators via their vector potentials.
The Climate Hierarchy as a Locked Network
Konstapel maps this formalism onto the climate system directly:
- Scalar potential ( \phi ) → ocean heat content, atmospheric energy density, cryosphere mass.
- Vector potential ( \mathbf{A} ) → zonal (trade winds, jet stream), meridional (thermohaline, Hadley cell), and vertical (convection) energy transport.
- External scalar drive ( S ) → solar irradiance variations at Schwabe (11 yr), Gleissberg (88 yr), and Suess–de Vries (210 yr) timescales.
- Coupling ( \kappa_{nm} ) → ocean–atmosphere interactions, teleconnections.
The Earth–Sun system thus becomes a nested hierarchy of oscillators, from weather (days) to ENSO (3–7 yr), AMO/PDO (55–70 yr), Gleissberg (88 yr), and Milanković (23–100 kyr). Critically, these are not independent. Each faster oscillator is embedded in and modulated by slower ones. This nesting is precisely what the quaternion coupling matrix formalizes.
Why Rational Ratios? Arnold Tongues
If two oscillators with natural frequencies ( \omega_1 ) and ( \omega_2 ) are weakly coupled, they tend to lock into a stable phase relationship when their frequency ratio is rational (( \omega_1/\omega_2 = p/q )). The region in parameter space where locking occurs is called an Arnold tongue; its width scales as ( K^{\max(p,q)} ), where ( K ) is coupling strength. Lower-order rationals (small ( p,q )) have the widest tongues and are thus the most robust against noise and dissipation.
Given the Schwabe cycle (11 yr) as the base frequency, the predicted observable periodicities are rational multiples of ( T_0 = 11 ) yr. Empirically observed cycles closely match these:
- 1/1 → 11 yr (Schwabe)
- 8/1 → 88 yr (Gleissberg) – error 0%
- 19/1 → 209 yr (Suess–de Vries) – error 0.5%
- 6/1 → 66 yr (AMO midpoint) – within 10% of observed 60–70 yr
Similarly, the Milanković ratios (23:41:100 kyr) approximate low-order rationals: 41/23 ≈ 9/5 (error 0.4%), 100/41 ≈ 5/2 (error 2.4%). The paper argues this is not coincidence but evidence of frequency locking between orbital modes through gravitational coupling.
Reading the Phase: Where We Are Now
The practical power of the oscillator model is phase reading. The linear forcing model has no phase; it only projects a trend. The quaternion oscillator model, by contrast, assigns a definite phase to each oscillator, propagating at its natural frequency.
Current state (mid-2020s):
- AMO: Peaked ~2005–2010, entered negative phase ~2020–2025. Minimum expected ~2040–2050, producing North Atlantic cooling of 0.2–0.4°C.
- Gleissberg solar cycle: Cycle 24 (2008–2019) was the weakest in a century. Projections (Zharkova et al.) suggest a deep minimum analogous to the Maunder Minimum around cycles 26–27 (2030–2053).
- Kondratiev economic wave (~50 yr): Downwave from 2000s peak, suggesting systemic restructuring.
- Milanković: Holocene interglacial ~11,000 years old, near or past its natural optimum. Orbital forcing moving toward next glacial inception on tens-of-millennia timescales.
The composite signal predicts a convergent cooling trough from ~2030–2060, with global mean surface temperature (GMST) 0.1–0.3°C below the linear warming trend. This is not a contradiction of anthropogenic warming—CO₂ forcing continues—but a superposition of a long-wave negative oscillation onto a positive linear trend.
What the IPCC Framework Misses
The IPCC’s ensemble mean of CMIP6 models, by design, averages over different realizations of internal variability, thereby erasing phase-coherent oscillatory signals. Individual model runs show oscillations; the ensemble mean removes them. For policy-relevant timescales (10–30 years), this is not a failure of science but a structural limitation of ensemble averaging when applied to a phase-coherent system.
The result is that projections of near-term temperature systematically misrepresent the oscillatory component. The quaternion oscillator model makes that component explicit, offering a falsifiable prediction: if GMST rises by more than 0.3°C above the linear trend between 2025 and 2060 without major volcanic forcing, then CO₂ forcing dominates the oscillatory signal at current concentrations. If instead a cooling trough appears, the oscillator model is supported.
Political Economy of the Linear Model
The paper does not shy away from the political economy of climate science. The emphasis on CO₂ as the primary driver has produced a powerful coalition of interests: renewable energy industries, carbon markets, international agreements, and research funding tied to the anthropogenic framing. Evidence for natural oscillatory forcing is often unwelcome because it implies that some warming is not attributable to human activity and therefore not removable by emissions reductions.
This does not mean the evidence is wrong or that CO₂ forcing is negligible. It means the scientific discourse has developed a structural bias toward the linear forcing model and against the oscillatory model—not because the evidence demands it, but because the policy implications are more tractable. The paper calls for an honest contest between the two frameworks, with the quaternion oscillator model providing the formal algebraic ground for that contest.
Three Testable Claims
The theoretical argument yields three empirical claims testable with existing data:
- Harmonic ratio clustering – In the EPICA Dome C ice core record (800 kyr), pairwise period ratios of dominant spectral peaks will cluster near low-order rational multiples of 11 years (Farey sequence ( \mathcal{F}_{12} )) more than expected under red noise. Expected: ≥70% of ratios within Arnold-tongue width of a low-order rational.
- Phase coherence fits 20th-century record – A three-parameter composite model (AMO + Gleissberg + linear trend) will fit the HadCRUT5 GMST record (1850–present) with root-mean-square residual ≤0.05°C, outperforming the CMIP6 ensemble mean (≥0.10°C at decadal scales).
- 2030–2060 cooling trough – GMST will be 0.1–0.3°C below the linear warming trend over 2025–2060. Falsified if GMST rises >0.3°C above trend without episodic forcing (volcanoes, etc.).
Conclusion
The paper does not argue that CO₂ does not warm the planet. It does. The argument is that climate is a hierarchically coupled oscillator system, that the oscillatory structure is formally grounded in Maxwell’s quaternion electrodynamics, and that this structure has been systematically underweighted in dominant climate models because ensemble averaging erases phase-coherent signals.
The immediate consequence is that near-term (10–30 year) projections based on ensemble-averaged forcing models are less accurate than projections based on oscillatory phase composites. The longer-term consequence is that energy policy built on the assumption of monotonic warming over the next 50 years will misallocate resources: cooling adaptation (agricultural shifts, heating infrastructure, cold-weather preparation) may be at least as important as warming adaptation in the North Atlantic region from 2030–2060.
The quaternion oscillator model is not a complete climate model. It does not replace general circulation models that resolve spatial structure. It is a phase model—describing when the climate system is in which phase of its long-period oscillations. Both are needed. But ignoring the oscillator structure entirely, as current policy frameworks largely do, is no longer scientifically justifiable.
Annotated Reference List for Further Reading
The following references are drawn from the original paper and supplemented with key background sources to allow readers to explore each topic in depth.
Maxwell’s Quaternion Electrodynamics (Formal Foundation)
- Maxwell, J.C. (1865). A dynamical theory of the electromagnetic field. Phil. Trans. R. Soc. Lond. 155, 459–512.
The original paper where Maxwell presents his electromagnetic theory using quaternions. Essential reading for understanding the scalar potential’s role before Heaviside’s vector reduction. - Maxwell, J.C. (1873). A Treatise on Electricity and Magnetism. Clarendon Press.
The full treatise. Chapter 9 (quaternions) and the treatment of the four-potential are directly relevant. - Jack, P.M. (2003). Physical space as a quaternion structure. arXiv:math-ph/0307038.
A modern mathematical physics treatment of quaternion algebra as a foundation for physical fields. Useful for readers without a background in Hamilton’s quaternions.
Climate Oscillations and Periodicities
- Enfield, D.B., Mestas-Nunez, A.M. & Trimble, P.J. (2001). The Atlantic multidecadal oscillation and its relation to rainfall and river flows in the continental US. Geophys. Res. Lett. 28, 2077–2080.
Definitive paper on the AMO index and its climatic impacts. NOAA maintains the updated index. - Knudsen, M.F. et al. (2011). Tracking the Atlantic Multidecadal Oscillation through the last 8,000 years. Nature Commun. 2, 178.
Paleoclimate reconstruction showing the AMO’s persistence over millennial timescales. Direct evidence against the “noise” hypothesis. - Gleissberg, W. (1958). The eighty-year solar cycle in auroral frequency numbers. J. Br. Astron. Assoc. 68, 148–152.
Original identification of the 88-year Gleissberg cycle. Outdated in method but historically important. - Peristykh, A.N. & Damon, P.E. (2003). Persistence of the Gleissberg 88-year solar cycle over the last 12,000 years. J. Geophys. Res. 108, 1003.
Cosmogenic isotope (¹⁴C, ¹⁰Be) evidence for the Gleissberg cycle’s persistence through the Holocene. - Jouzel, J. et al. (2007). Orbital and millennial Antarctic climate variability over the past 800,000 years. Science 317, 793–796.
The EPICA Dome C deuterium record. Primary dataset for testing harmonic ratio clustering. - Milanković, M. (1941). Canon of Insolation and the Ice-Age Problem. Royal Serbian Academy.
The original orbital forcing theory. The rational ratios (23:41:100 kyr) are derived here, though not explained as frequency locking.
Frequency Locking and Arnold Tongues
- Arnol’d, V.I. (1983). Geometrical Methods in the Theory of Ordinary Differential Equations. Springer.
The mathematical foundation of Arnold tongues. Chapter on perturbation theory and circle maps is directly applicable. - Holling, C.S. (2001). Understanding the complexity of economic, ecological, and social systems. Ecosystems 4, 390–405.
Introduces panarchy theory—nested adaptive cycles—which is conceptually parallel to the hierarchical oscillator model. Useful for readers interested in cross-domain applications.
Solar Activity Predictions
- Zharkova, V.V. et al. (2015). Heartbeat of the Sun from principal component analysis and prediction of solar activity. Sci. Rep. 5, 15689.
The controversial but data-rich prediction of a deep solar minimum in cycles 26–27. Directly cited in the paper for the Gleissberg phase projection. - Hoyt, D.V. & Schatten, K.H. (1998). Group sunspot numbers: a new solar activity reconstruction. Sol. Phys. 179, 189–219.
Standard reconstruction used for Gleissberg cycle fitting before the satellite era.
Temperature Records and Model Comparisons
- Morice, C.P. et al. (2021). An updated assessment of near-surface temperature change from 1850: the HadCRUT5 dataset. J. Geophys. Res. Atmos. 126, e2019JD032361.
Primary GMST record for testing the composite model. - Ollila, A. (2022). Two main temperature periodicities related to planetary and solar activity oscillations. Int. J. Climatol. 42, 10392–10407.
A recent independent empirical study finding harmonic periodicities in temperature records. Supports the Arnold-tongue interpretation without using quaternion formalism.
Extended Background (Complex Systems and Economic Cycles)
- Kondratiev, N.D. (1925). The major economic cycles. Voprosy Konjunktury 1, 28–79.
Original long-wave economic cycle paper. Cited for the Kondratiev wave component of the composite model. - Bejan, A. (2000). Shape and Structure, from Engineering to Nature. Cambridge University Press.
Constructal law and hierarchical flow systems. Provides a thermodynamic rationale for why coupled oscillator hierarchies emerge in geophysical flows. - Rowlands, P. (2007). Zero to Infinity: The Foundations of Physics. World Scientific.
Algebraic grounding via Clifford algebra Cl(3,1), which is isomorphic to quaternion algebra over the reals. Referenced in the paper as the companion algebraic foundation.