Climate — Related Work

For the conclusion this page positions, see Climate Summary. For the empirical anchors, see Climate Formula and Insolation Null Test.

Closest precedent — Dutkiewicz et al. 2024

“Deep-sea hiatus record reveals orbital pacing by 2.4 Myr eccentricity grand cycles” — Nature Communications (2024). DOI: 10.1038/s41467-024-46171-5 

The closest published parallel to the model’s gravitational-coupling claim:

• Explicitly identifies Mars–Earth gravitational resonance (the g₄ − g₃ secular beat) as a direct climate driver — framed as gravitational forcing, not insolation.
• Uses 65 Myr of satellite-mapped deep-sea sediment hiatus data; finds gaps spaced at 2.4-Myr intervals matching the Mars–Earth resonance.
• Frames mechanism as “Mars’ gravitational pull is shifting Earth’s path around the sun”.

Where the model aligns: same physical mechanism (planet–planet secular coupling → climate variance), same direction of causation.

Where the model extends further: Dutkiewicz documents one cycle (2.4 Myr = g₄−g₃). The Holistic model documents 32 cycles unified under one fundamental period — the 2.4-Myr cycle is just one of them. It also documents the 1.2-Myr obliquity-band Grand cycle (s₄−s₃) and the 9-Myr 13H cycle, neither of which Dutkiewicz covers.

Direct revision of classical tenets — Zeebe & Lantink 2024

“A secular solar system resonance that disrupts the dominant cycle in Earth’s orbital eccentricity (g₂−g₅): Implications for astrochronology” — The Astronomical Journal 167(5), 204 (2024). DOI: 10.3847/1538-3881/ad32cf  · arXiv:2403.09332 

The paper states explicitly:

“The paradigm that the longest Milanković cycle dominates Earth’s astronomical forcing, is stable, and has a period of ~405 kyr requires revision.”

The mechanism: a secular resonance σ₁₂ = (g₁−g₂) + (s₁−s₂) can destabilise the g₂−g₅ beat over long timescales without major planetary-orbit changes. During σ₁₂-resonance episodes, the 405-kyr line is weak or absent.

Where the model aligns: the Holistic framework already treats 405-kyr as off-lattice (Climate Formula) — it is not an integer divisor of 8H. It enters the formula as the L2 carbon-thermostat fundamental, not as an L1 lattice integer. Structurally consistent with Zeebe–Lantink: if 405-kyr is dynamically unstable, it shouldn’t carry the metronome role classical astrochronology assigns it.

Where the model goes differently: the reason 405-kyr is off-lattice in this framework is integer-divisor structure, not dynamical instability. Both explanations may be correct simultaneously — the 405-kyr beat doesn’t fit the L1 fundamental-period lattice AND is dynamically unstable on long timescales.

The skeptical bookend — Wunsch 2003 / Roe 2006

The 2000s saw a substantive debate over whether orbital forcing explains climate at all.

Wunsch 2003 — “The spectral description of climate change including the 100 ky energy” — Climate Dynamics. DOI: 10.1007/s00382-002-0279-z 

Argued that climate records show “red-noise process or random walk” dominance, that Milankovitch frequencies contribute only a “small fraction of total climate variance”, and that the 100-kyr peak can be rationalised without invoking orbital forcing (threshold dynamics suffice).

Roe 2006 — “In defense of Milankovitch” — Geophysical Research Letters. DOI: 10.1029/2006GL027817 

Counter-argument: reframed Milankovitch as a rate-of-change relationship; ice-volume rate-of-change is in zero-lag phase with NH summer insolation; CO₂ lags ice-volume rate-of-change → ice melting precedes CO₂ change.

Position relative to this debate: the model’s R² = 0.87 on post-MPT LR04 is incompatible with Wunsch’s red-noise null and consistent with Roe’s defence of orbital forcing. It uses a different parameterisation basis than either (the L1 lattice rather than ε / e / ϖ insolation features), which makes the channel question (insolation vs broader gravitational coupling) explicit rather than implicit.

Grand-cycle precedents — Boulila 2019, Saillenfest 2020s

Boulila 2019 — “Coupling between Grand cycles and Events in Earth’s climate during the past 115 million years” — Scientific Reports 9, 327. DOI: 10.1038/s41598-018-36509-7 

Documented ~9 Myr and ~36 Myr Grand cycles in 115-Myr benthic foraminifera δ¹⁸O. Major Cenozoic climatic events (PETM, EOT, Mi-1, etc.) occur preferentially at extremes of these cycles.

Where the model aligns: the Holistic model has the 9-Myr cycle as 13H = 9.04 Myr. A stability test on 13H found it does not behave as a single coherent oscillator (amp CV 42–50%, circular phase std 97.9°), so it is not promoted to canonical L1/L2 — consistent with the Grand cycle being a multi-mechanism phenomenon.

Saillenfest et al. (a series of 2019–2021 papers on Saturn’s obliquity drift and Neptune resonance) work in the same secular-coupling-as-climate-driver framing.

Methodological neighbour — cyclostratigraphy (Hinnov, Acycle)

Hinnov 2018 — “Cyclostratigraphy and Astrochronology in 2018” — Stratigraphy & Timescales 3:1–80. DOI: 10.1016/bs.sats.2018.08.004 

Cyclostratigraphy is the established discipline of fitting multiple orbital sinusoids to paleoclimate records — the methodological neighbour of the model’s L1 lattice fit. The Acycle software, the Cyclostratigraphy Intercomparison Project, and 40+ years of multi-tone-fit papers (Imbrie & Imbrie 1980, Hays-Imbrie-Shackleton 1976, Huybers 2007 ) routinely reach R² 0.6–0.9 on post-MPT LR04 and similar records.

Where the model aligns: same method class — sinusoidal regression at known orbital frequencies, sequential-fit decomposition.

Where the model differs: cyclostratigraphy fits ~10–15 individual Laskar eigenmode beats as independent oscillators. The model fits 32 integer divisors of a single fundamental period (8H) — a stronger structural claim. The 405-kyr metronome is treated as the universal anchor in cyclostratigraphy; the model treats it as off-lattice L2.

The channel question — Munk et al. 2002

Munk, Dzieciuch & Jayne 2002 — “Millennial Climate Variability: Is There a Tidal Connection?” — Journal of Climate 15:370. Link 

Munk made an observation the model implicitly answers:

“The tide community is concerned with the relatively rapid gravitational forces (periods up to 18.6 yr) and the climate community with the long-period Milankovitch insolation terms (exceeding 20,000 yr).”

The gravitational vs insolation channel split is a known disciplinary gap. The two communities don’t talk to each other much because their forcings are at different periods. The model’s framing implicitly closes the gap: gravitational coupling is the common source for both channels.

Climate sensitivity — paleo references

The model’s frequency-resolved decomposition of LR04 ice-volume variance into ice-albedo versus CO₂/GHG forcing per orbital band (see Climate Formula §Ice-Albedo Decomposition) yields a Charney equilibrium climate sensitivity estimate. Marginalized over the four calibration constants (Monte Carlo, N = 5,000 draws), the result is Charney ECS = 3.63 K, 90 % CI [3.01, 4.31] K — consistent with all four major published reference values within their respective uncertainties.

Where the model aligns: the Holistic value sits within all four published reference ranges. Hansen 2013’s 3.0 ± 0.5 K central value, Sherwood et al. 2020’s 2.6–3.9 K Bayesian-synthesis range, PALAEOSENS / Köhler 2017’s 3.0–4.5 K multi-proxy paleo range, and IPCC AR6’s likely 2.5–4.0 K range all bracket the Holistic 3.63 K estimate.

Where the model differs methodologically: the Holistic estimate comes from an orthogonal direction — a fixed-lattice frequency-domain decomposition rather than the standard LGM-to-Holocene time-domain comparison (Hansen 2013) or Bayesian synthesis over independent lines of evidence (Sherwood 2020). The lattice-based estimate also produces a frequency-resolved ice-albedo decomposition that none of the surveyed references produce: 57 % ice-share at the 100-kyr band, 64 % at obliquity, 69 % at precession, 46 % at long periods (>130 kyr). The same per-band decomposition computed pre-MPT shows a +31 percentage-point obliquity-band ice-share elevation, quantitatively confirming Willeit et al. 2019 ’s mid-Pleistocene-transition framing.

Novelty assessment

Novel claims (not found in surveyed literature)

Aligned with active 2024 research

Differences from the conventional view

An observational angle: the 405-kyr phase position of the current era

A fitted output of the framework’s L2 layer (Climate Summary §5) is that the 405-kyr eccentricity-driven carbon thermostat is near a warm peak around the present era. Two sides of where this sits in the literature:

• Mainstream Holocene attribution emphasises the 23.7-kyr precession cycle: NH summer insolation peaked at the Holocene Climate Optimum (~10 kyr BP) and has been declining since, implying a cooling natural trend. The IPCC AR4 paleoclimate chapter and follow-ups discuss this extensively. Lorenz et al. 2006  found orbital forcing has been “up to four times larger than the 1.6 W m⁻² net anthropogenic forcing since 1750” — in the same precession-emphasising tradition.
• The 405-kyr cycle’s current phase position is rarely cited in modern-era attribution discussions. The cycle is studied extensively in deep-time cyclostratigraphy (Hinnov, Laskar, Zeebe–Lantink) — but typically for events tens or hundreds of millions of years old, not for what the present-era phase implies about pre-industrial baseline. The “Holocene Temperature Conundrum ” debate touches this question but does not resolve the role of the 405-kyr phase position.

The framework’s fitted L2 phase provides one quantitative element to the Holocene-temperature debate that the literature treats less prominently than the precession signal. This is a positioning observation, not an attribution claim — the framework’s rate-of-change implications cannot account for industrial-era warming (see Climate Summary §5 for the rate vs level distinction).

See also

• Climate Summary — the conclusion this page situates
• Climate Formula — canonical L1+L2+L3 architecture
• L1 Attribution — per-integer Berger vs Holistic mapping
• Insolation Null Test — empirical anchor for the ΔR² = 0 claim