Milankovitch orbital forcing: F(t) = sum A_k cos(omega_k t); spectrum on 3 frequencies

Milankovitch spectral-forcing framework. Setup: Earth's insolation at latitude phi is modulated on ~10 kyr timescales by three orbital parameters: eccentricity (~100 kyr), obliquity (~41 kyr), and precession (~23 and ~19 kyr) (Milankovitch…