This post proposes that regular variations in interplanetary dust between the Earth and the Sun is the primary cause of recent ice ages whose oscillations have been parameterised here. The Milankowitz theory of insolation is able to well explain the continuous 41,000y climate oscillation but it cannot explain why a 100,000y cycle became prevalent during the last million years, even though it is in phase with orbital eccentricity of the Earth’s orbit around the sun. A different origin for the emergence of the 100,000y cycle of major ice ages is proposed. This is based on the hypothesis that a large comet broke up close to Earth’s orbit about 3.3 million years ago, leading to an increase in inter-planetary dust whose density eventually phase locked to variations in Earth’s eccentricity. This process is able to explain the cooling trend starting 3 million years ago eventually leading to major eccentricity driven glaciations with persistent modulation in insolation from the obliqueness cycle. Evidence in support of regular cycles in inter-planetary dust density have been found in ocean sediments, and there is some evidence of a comet fragment impacting with Earth at 3.3 million years ago.
Milankowitch cycles inducing insolation changes in the polar regions is the generally accepted primary cause of Ice Ages. However the details for insolation changes due to the ellipticity of the Earth’s orbit don’t really work out. Although the insolation varies by 6%(today) up to 20% between perihelion(currently southern hemisphere summer) and aphelion(currently northern hemisphere summer), the net annual change for the Earth is insignificant. The poles themselves are unique places because they essentially experience just 1 day and 1 night per year. For the North Pole the sun rises at the spring equinox and sets again at the autumn equinox. At mid summer’s day the 24 hour insolation reaches average values higher than at the equator. The winter night then follows resulting in zero insolation for the next 6 months. This also explains the rapid melting and freezing cycles observed in polar regions.
Data used for insolation calculations at the poles are the three Milankowitz cycles of Orbital Eccentricity, Obliqueness and Precession originally calculated by  are taken from the Nasa site. The insolation is then calculated by tracking the orbit through the year on an hourly basis. Sample results are shown below.
The effects on maximum insolation, which effects summer melting rates of ice at the poles, is evident for both obliqueness and precession. However the effect of changes in the orbital ellipticity is more subtle. An increase in insolation for one pole close to perihelion is offset by a diminished insolation for the other pole. Later in the precession period when the poles are positioned at the equivalent of todays equinoxes, the polar insolation is similar to that of a circular orbit. Overall the net annual effect of eccentricity of the Earth’s insolation is negligible. For the last 900,000 years however, Ice Ages show a clear correlation with ellipticity at 100,000 years and this remains unexplained. The basic problem can be seen better by calculating long term insolation trends. Figure 2 shows the maximum and total insolation for the poles over the last 600,000 years.
The main effect of ellipticity is to modulate the effect of the precession terms at 23,000 years. Indeed if the orbit was circular then the precession of the equinoxes would have no effect whatsoever. However the direct effect on insolation of the increase in ellipticity of the Earth’s orbit averaged over a year is negligible. This is the paradox as to why regular glaciations and warm interglacials occur every 100,000 years in phase with the eccentricity.
The clearest insolation effect on climate is the obliqueness signal which when averaged out over the combined ellipticity and precession yields a regular variation of around 25 watts/m2 in maximum flux. The delta 18O data show a clear signal throughout the 5 million year period in phase with obliquity. The precession signal is directly dependent on eccentricity and is less evident in the data. However there is no explanation for the primary 100,000 year variation in phase with ellipticity. An excellent review of this problem has been given by Maya Elkibbin & JoseA. Rial  which highlights the questions that any complete theory of the evolution of Ice Ages must answer namely:
1. The observed gradual decrease in temperatures starting just over 3 million years ago.
2. The 43,000 year continuous signal dominated until about 1 million years ago when the onset of 70-100,00000 glacial cycles begin. Why did this happen ?
3. What is the cause of the recent 100,000y glacial cycles, since Insolation changes caused by ellipticity are too small to be an explanation ?
4. Why did the frequency change from about 70,000 years between 900-700,000 ybp to the current cycle of 100,000 years ?
5. Why is the larger 400,000 year eccentricity signal absent from recent data ?
Variations in insolation can explain the 43,000y obliquity and 23,000 y precession signals observed in the data. However another explanation is needed to explain the points made above. A possible solution is proposed as described below
Muller et al. have proposed an alternative explanation for the 100,000y cycle . They suggested that the driver is the change in inclination of the Earth’s orbit. This follows an 80,000 year cycle relative to the current ecliptic plane, but when projected to the invariate plane defined as the centre of mass plane of the solar system this is extended to 100,000 years. The physical process envisaged is that the interplanetary dust cloud is centerd on the invariate plane and the Earth passes through this during its orbit. Evidence from ocean sediments show a 100,000 year cycle of He3 deposits associated with cosmic dust . In effect, changes to orbital parameters of eccentricity, inclination and longitude of ascending node are all caused by the mutual gravitation of solar system planets as they orbit the Sun. Similarly all bodies in the solar system will follow such cycles including dust clouds. Let’s investigate correlations of inclination and glaciations. Figure 4 shows the data as given by  compared to the benthic foram stack.
The hypothesis that the Earth passes through a static interplanetary dust cloud aligned with the invariate plane attenuating sunlight is attractive. Furthermore there is direct evidence in ocean sediments of a regular 100,000 cycle of He3 deposits . Unfortunately the data do not seem to support inclination as being responsible, despite the similarity in timespans. One would expect colder temperatures at low alignment where presumably the Earth crosses the dust plane. The correlation with eccentricity is better with high eccentricity seemingly algning with warmer interglacial periods. So what is really happening and can we explain how the climate changed to regular glaciations ?
First we try and estimate if there is today an observable effect of interplanetary dust attenuating the solar constant as measured accurately by satellite.
The surface temperature of the sun is 5778 K and the mean distance of Earth to the Sun is 1.496×108 km. Applying Stefan Boltzman’s law we can calculate the radiation flux (solar constant) out side the atmosphere. Radius of the Sun = 6.96*10^4km
S0 = (sigmaT^4)*4piR^2/4pi*D^2 putting the numbers in gives:
S0 = 1368 watts/m2
The measured average solar constant over the orbit is 1361 watts/m2 . This implies therefore that currently during an interglacial maximim on average 7 watts/m2 is absorbed (or scattered) by interplanetary dust currently lying between the earth and the Sun. In addition accretion of dust by the atmosphere which is estimated at 40000 tons per year will also slightly reduce sunlight incident on earth. Therefore it is proposed that variations in the dust cloud aligned with the Earth – Sun ecliptic plane is the origin of recent Ice Ages. Until 1 million years ago, colder glacial periods were driven by variations in obliquness induced insolation. However this changed presumably when the dust cloud settled into harmony with the Earth’s orbit. For the last million years “Milankowitch” variations in the eccentricity of the orbit have also applied to the dust cloud inducing regular variations in density distributions thus causing variations in the solar constant. A hypothesis that can explain the slow cooling of the Earth, the onset of ice ages and their transition from 43,000 y to 100,000y dominant cycles is as follows.
Some 3.3 million years ago a large comet breaks up near the Earth creating a large irregular dust cloud lying between the earth and Sun. There is some direct evidence of a collision of a significant cometry fragment with the Earth at this time from Argentina . Initially this irregular dust cloud orbited the sun and acted to gradually reduce temperatures on Earth as it condensed into a Earth-Sun ecliptic orbit. The cloud gradually spread out into a disk in synchronous orbit with the Earth around the Sun. Variations in the orbital eccentricity of the Earth are caused by the regular combined gravitational forces of the outer solar system planets. Exactly the same forces will apply to the dust cloud now in synchronous orbit. However, the low mass of particles in the cloud accentuate the orbital effects especially in eccentricity. Thus we can envisage a current dust cloud whose density and shape is moulded by those same eccentricity forces as that on the Earth. Initially the diffuse cloud caused global cooling for 2 million years after breaking up. Slowly as the dust cloud settled into the ecliptic plane it became resonant with the orbital effects of the Earth. When the Earth is at a maximum of eccentricity in the 100,000y cycle the dust cloud is spread in an elliptical cloud with much of it lying outside the Earth Sun line of sight. As the Earth resorts back to low eccentricity and a more circular orbit so too does the dust cloud leading to an increase in dust density lying between the Earth and the Sun. This effect is shown schematically below.
If it is assumed that the net effect of the current dust density between Earth and the Sun is such as to reduce the solar constant by 7 watts/m2, as described above, then a doubling of this density as eccentricity drops to a minimum would cause a global temperature drop of about 2K leading to a new an Ice Age. This is based on a Stefan Boltzman respons of : DS = 4sigmaT^3 DT
This model of an interplanetary dust cloud resonant with the Earth’s eccentricity cycle predicts that the following phenomona should be observable.
- A slow decrease in the solar constant as the density of dust increases with falling eccentricity.
- Evidence of more dust external to the earth’s orbit along the aphelion axis compared to the aehlion axis.
- A higher density of dust between the Earth and the Sun at perihelion rather than at aphelion. This could be reflected in a slightly lower value of solar constant at perihelion than otherwise predicted. In addition the Zodiacal light should be more intense at perihelion than aphelion.
A doubling of CO2 levels in the atmosphere by the end of this century will cause a direct radiative forcing of 5.3Ln(2) watts/m2 or 3.67 watts/m2. IPCC models usually assume a further feedback mainly from water vapor of 2 watts/m2/degC. A reasonable assumption would be that AGW would induce a warming term of about 5-6 watts/m2 after which CO2 levels are expected to stabilize. The net effect of all human induced global warming could be to increase radiative forcing by something like between 4 and 8 watts/m2.
The eccentricity of the Earth’s orbit is currently decreasing with cooling beginning in 2000 years time leading eventually to another intense Ice age after 20,000 years and lasting a further 60,000 years. A new Ice Age would be disastrous for human and animal life at large latitudes. If this theory of a variable interplanetery dust cloud is correct, then AGW will likely delay the early onset of the next ice age. To fully offset the cooling effects of a thicker dust cloud, it would be necessary to keep CO2 levels above about 700 ppm for another 70,000 years.
A new possible cause of Ice Ages is proposed based on an increase in interplanetary dust for lower eccentricities of the Earth’s orbit for the last million years leading to large cooling of the Earth. This is proposed as the cause of the switch from glaciations driven by the 43,000 year variation of insolation due to obliquity to a 100,000 eccentricity driven cycle. The next cooling cycle should begin in 2000 years time reaching a glacial minium in about 10,000 years time lasting 70,000 years. At the same time anthropogenic global warming may delay the cooling trend if induced temperature rises of 2-3 degrees occur and defer the next glaciation if CO2 levels were to stabilize at ~700 ppm for 80,000 years.
 T.R. Quinn, S. Tremaine, & M.A. Duncan A three million year integration of the Earth’s orbit, Astronomical Journal 101, 2287Ð2305, 1991
 http://aom.giss.nasa.gov/srorbpar.html Andre L. Berger, 1978. Long Term Variations of Daily Insolation and Quaternary Climatic Changes. Journal of the Atmospheric Sciences, volume 35(12), 2362-2367
 Maya Elkibbi), Jose ? A. Rial, An outsider’s review of the astronomical theory of the climate: is the eccentricity-driven insolation the main driver of the ice ages? Earth-Science Reviews 56 Ž2001. 161–177
 Higgins et al. Sediment focusing creates 100-ka cycles in interplanetary dust accumulation on the Ontong Java Plateau, Earth and Planetary Science Letters 203 (2002) 383-397
 P. H. Schultz, M. Zarate, W. Hames, C. Camilión and J. King, A 3.3-Ma Impact in Argentina and Possible Consequences, Science 11 December 1998: Vol. 282 no. 5396 pp. 2061-2063
 Surface Temperature of the Sun http://en.wikipedia.org/wiki/Sun
 Kopp, G. and J. L. Lean, A new, lower value of total solar irradiance: Evidence and climate significance, Geophys. Res. Lett., 38
 “Glacial Cycles and Astrnomical Forcing” Richard Muller and Gordon Macdonald, SCIENCE Vol 277 Pages 215 – 218 11th July 1997