  The Earth's revolution is modified by changes in inclination of its rotation
axis. Despite the fact that the gravity field is central, the Earth's
trajectory is not closed and the equinoxes drift. Milankovic (1920) argued that
the shortest precession period of solstices is 20,7kyr: the Summer solstice in
one hemisphere takes place alternately every 11kyr at perihelion and at
aphelion. We have submitted the time series for the Earth's pole of rotation,
global mean surface temperature and ephemeris to iterative Singular Spectrum
Analysis. iSSA extracts from each a trend, a 1yr and a 60yr component. Both the
apparent drift of solstices of Earth around the Sun and the global mean
temperature exhibit a strong 60yr oscillation. The "fixed dates" of solstices
actually drift. Comparing the time evolution of the Winter and Summer solstices
positions of the rotation pole and the first iSSA component (trend) of the
temperature allows one to recognize some common features. A basic equation from
Milankovic links the derivative of heat received at a given location on Earth
to solar insolation, known functions of the location coordinates, solar
declination and hour angle, with an inverse square dependence on the Sun-Earth
distance. We have translated the drift of solstices as a function of distance
to the Sun into the geometrical insolation theory of Milankovic. Shifting the
inverse square of the 60yr iSSA drift of solstices by 15 years with respect to
the first derivative of the 60yr iSSA trend of temperature, that is exactly a
quadrature in time, puts the two curves in quasi-exact superimposition. The
probability of a chance coincidence appears very low. Correlation does not
imply causality when there is no accompanying model. Here Milankovic's equation
can be considered as a model that is widely accepted. This paper identifies a
case of agreement between observations and a mathematical formulation.
