Ice Ages revisited

In 2006 Gerard Roe published a paper “In defense of Milankovitz” [1] which showed a good qualitative agreement between the rate of change in ice volume with Arctic summer melting as measured by maximum insolation at 65N. The ice volume (temperature proxy) as measured by the Benthic DO16 stack is not correlated directly with insolation but instead interglacials are better correlated with the eccentricity cycle.

I have re-anlaysed data from a previous analysis to look for a correlation with the differential of Bentic DO16 data [2]. I use an orbital calculation of maximum insolation at the North pole to look for correlations. The result if shown in Figure 1.

Figure 1: Comparison of Ice Volume (Temperature proxy) from Benthic Fora data. Upper graph Rate of Change DC/DT overlayed with maximum 'm-d summer' insolation at North Pole (watts/m2). Shown below are the Ice Volume data

Figure 1: Comparison of Ice Volume (Temperature proxy) from Benthic Fora data. Upper graph Rate of Change DC/DT overlayed with maximum ‘m-d summer’ insolation at North Pole (watts/m2). Shown below are the Ice Volume data

There is a clear tendency that with Milankowitz increases in summer insolation at the North pole induces decreases in Ice volume. However it can’t explain why 100,000 year interglacial periods have prevailed for the last 600,000 years.

Fig2 - same as Fig8 in  ref [3]

Fig2 – same as Fig8 in ref [3]

Ruddimen [3] argues that CO2 changes induced by Ocean temperatures should be seen as a feedback on Ice Volume rather than a cause direct cause. The 41,000 year changes in tilt effect the latitude range insolation. Interglacials seem to occur at multiples 41,000 years. For 2 million years until 900,000 years ago interglacials were in phase with obliquity. Then as the climate continued to cool they switched first to 2 obliquity cycles and then to 3 cycles from 600,000 years ago. Ruddimen argues that the current saw tooth oscillations are caused by non-linear CO2 feedbacks. – see Figure 2.

CO2 is therefore not a cause of glaciations, but instead acts as a feedback on the primary cause which is Milankowitz orbital variations in insolation. Summer insolation reduces Ice volume and ocean temperatures rise. This increases CO2 levels in the atmosphere. Any enhanced greenhouse effect is a logarithmic response to solar forcing depending on temperature. It seems to need 3 tilt cycles before the induced CO2 response is large enough to trigger an interglacial.


1) Gerard Roe, In defense of Milankovitch, Geophysical Research Letters, Vol 33, L24703, 2006

2) Lisiecki, L. E., and M. E. Raymo (2005), A Pliocene-Pleistocene stack of 57 globally distributed benthic d18O records, Paleoceanography, 20, PA1003

3) W.F Ruddiman, Ice-Driven CO2 feedback on ice volume, Climates of the Past Discussions, 2 43-78, 2006

About Clive Best

PhD High Energy Physics Worked at CERN, Rutherford Lab, JET, JRC, OSVision
This entry was posted in Climate Change, climate science, Ice Ages, Science. Bookmark the permalink.

6 Responses to Ice Ages revisited

  1. AJ says:

    “The ice volume (temperature proxy)… “.

    Wouldn’t ice volume be a proxy for “energy in the climate system”, which in turn would be the integral of the surface air temperature? It seems to me that ice will melt fastest when temperatures are highest. CO2 would be also be an energy proxy. If this is true, then CO2 would lag temperatures by a much longer period.

    This paper, “Physical Basis for the Temperature-Based Melt-Index Method” might be interesting:

    Perhaps I’m off base? That is, due to coupling, temperature and energy might be nearly perfectly correlated. At this point, however, I doubt it.

    • Clive Best says:

      By CO2 being an energy proxy do you mean that it follows sea surface temperature ?

      Ice melting clearly depends on air temperature which rises as the 4th power of insolation. The highest insolation at the north pole occurs fo largest obliquity and when summer coincides with the orbital perihelion.

  2. AJ says:

    That should have been: “That is, due to coupling, temperature and energy might be nearly in-phase”

  3. AJ says:

    Then again, maybe the 65N and tropical temperatures are out of phase. Maybe ice volume is a good proxy for tropical temperatures and thus global temperatures. It’d be nice if someone somewhere could explain the ice volume = temperature relationship. I must be missing something. It seems I’m the only one interested in this question.

    • Clive Best says:

      I think the relationship between Ice volume and temperature must be non-linear. Something strange happens to trigger an interglacial. Currently inter-glacials are triggered every 3 cycles of obliquity. 3 million to 1 million years ago they were triggered every cycle of obliquity (41K years). Currently the earth has cooled to a stage where saw-tooth oscillations occur roughly every 100,000 years. This means that the ice melts fast once an instability occurs. This could be due to an enhancement above Milankovitz forcing due to albedo/ clouds or something else.

  4. John Reid says:

    Hi Clive

    I have just discovered your blog – nice to see something on climate by a fellow physicist.

    I think comparing peaks and troughs between time series is not a good way to go – it is too easy to fool yourself. Not only that but there are a lot of people out there with pet theories about “cycles” in climate and other fields. Short lengths of random time series with red spectra often give the impression of being deterministically cyclic even when they are not. See my unpublished paper on this at
    I am working this up into a full stochastic model of ice-age climate. I was researching radiation transfer vs lapse rate when I came across your blog.



Leave a Reply