
to deflect atmospheric planetary waves (Huybers and Tziper-
man, 2007). This changed the storm path across the North
Atlantic Ocean and prevented Gulf Stream penetration as far
north as today. This surface ocean change and increased melt-
water in the Nordic Seas and Atlantic Ocean ultimately reduced
the production of deepwater. This in turn reduced the amount
of warm water pulled northwards. All of these led to increased
cooling in the Northern Hemisphere and expansion of the ice
sheets.
There are other feedback mechanisms that also helped drive
the system towards maximum possible glacial conditions
(Raymo and Huybers, 2008). These include changes in the car-
bon cycle that reduced both atmospheric carbon dioxide and
methane. Ridgwell et al. (2003) discuss the main controls on
the carbon cycle and the possible causes of reduced glacial
atmospheric carbon dioxide. They speculate that glacial-inter-
glacial changes in atmospheric carbon dioxide may be primar-
ily driven by changes in oceanography in the Southern Ocean.
This could have altered the nutrient supply to the surface
waters, and hence surface water productivity. Dust-stimu lated
increased glacial surface water productivity would have drawn
down atmospheric carbon dioxide into the surface water to pro-
duce organic matter through photosynthesis. However, the con-
trols on the glacial-interglacial carbon cycle are still very
poorly understood. Glacial periods are also by their very nature
drier, which reduces atmospheric water vapor. For example,
Lea et al. (2000) provide clear evidence that the water vapor
production of the equatorial Pacific zone was greatly curtailed
during the last five glacial periods. CO
2
,CH
4
and water vapor
are all crucial greenhouse gases and any reduction in them
leads to general global cooling (Figure Q5), which in turn
furthers glaciation.
These feedbacks are prevented from becoming a runaway
effect by “moisture limitation.” As warm surface water is
forced further and further south, the moisture supply that is
required to build ice sheets decreases. Moreover, it is very clear
that ice sheets are naturally unstable and these feedbacks are
constantly altering direction during the whole glacial period,
hence it takes 80 ka to achieve the maximum ice extent during
the LGM and that period is characterized by rapid oscillations
such as Heinrich events and D-O cycles (see below).
The natural instability of ice sheets means that deglaciation
is much quicker than glaciation. For example, the last deglacia-
tion or Termination I, endured a maximum of 4 ka including
the brief return to glacial conditions called the Younger Dryas
period. The increase in summer 65
N insolation (e.g., Imbrie
et al., 1993) led to the initial melting of the northern ice sheets.
The resultant raised sea level then undercut the ice sheets adja-
cent to the oceans, which in turn increased sea level. This sea
level feedback mechanism was extremely rapid. Once the ice
sheets retreated then the other feedback mechanisms discussed
for glaciation reversed (see Figure Q5), leading to a massive
increase in the greenhouse gases CO
2
,CH
4
and water vapor.
These feedbacks were prevented from creating a runaway effect
by the limit of how much heat the North Atlantic could steal
from the South Atlantic to maintain the interglacial deep-water
overturning rate.
The latest data concerning the link between orbital forcing
and these feedbacks suggests that changes in the carbon cycle
lead to changes in global ice volume (Shackleton, 2000; Ridgwell
et al., 2003; Ruddiman and Raymo, 2003; Pagani et al., 2005;
Raymo and Huybers, 2008). This suggests that orbital forcing
has a direct effect on atmospheric carbon dioxide and methane,
which may drive global temperatures, and in turn change glo-
bal ice volume. Thus, the carbon cycle, instead of being a
response to glacial-interglacial cycles, could indeed be one of
the key driving forces (see Figure Q5).
The Mid-Pleistocene revolution
Debunking the eccentricity myth?
The Mid-Pleistocene Revolution (MPR) is the term used to
denote the marked prolongation and intensification of the glo-
bal glacial-interglacial climate cycles which occurred between
900 and 650 ka (Berger and Jansen, 1994). Prior to the MPR,
since at least the intensification of Northern Hemisphere glacia-
tion (2.75 Ma), global climate appears to have responded to
obliquity variations (Imbrie et al., 1992), i.e., the glacial-intergla-
cial cycles occurred with a frequency of 41 kyr. After about
800 ka, glacial-interglacial cycles became more pronounced
and occurred with a frequency of approximately 100 ka. The
MPR, therefore, marks a dramatic sharpening of the contrast
between warm and cold periods. Mudelsee and Stattegger
(1997) used advanced methods of time-series analysis to
review the deep-sea evidence spanning the MPR and summar-
ized the salient features (Figure Q6). The first transition
occurred between 942 and 892 ka when global ice volume sig-
nificantly increased. However, the 41 ka climate forcing contin-
ued. This situation persisted until about 650–725 ka when the
climate system entered a bi-modal phase and the strong 100 ka
climate cycles began (Mudelsee and Stattegger, 1997).
Imbrie et al. (1993) called the MPR the “100 ka problem,”
as the eccentricity signal is by far the weakest of the orbital
parameters. Hence, it has been assumed that there must have
been a change from a linearly to a non-linearly forced climate
system (Imbrie et al., 1992, 1993). Many different theories
have thus been postulated to produce this critical non-linear
transition, including: critical size of the Northern Hemisphere
ice sheets, global cooling trend, erosion of regolith beneath
the Laurentide Ice Sheet, o rbital inclination, Greenland-
Scotland submarine ridge and carbon cycle and atmospheric
CO
2
(see Maslin and Ridgwell, 2005 for detailed references).
There are a number of problems with associating the 100 ka
cycles with eccentricity. The first is that eccentricity has spec-
tral peaks at 95 ka, 125 ka and 400 ka. The spectral analysis
of benthic foraminifera oxygen isotopes (SPECMAP) – a
proxy for global ice volume – shows a consistent single peak
of 100 ka (Maslin and Ridgwell, 2005). However, the lengths
of the last 8 glacial-interglacial cycles vary from 87 to 119 ka,
averaging about 100 ka. Ridgwell et al. (1999) demonstra-
ted that this strong 100 ka peak could not be produced by
eccentricity, but rather by a simple saw-tooth pattern based on
a long glaciation period followed by a short deglaciation
(Figure Q7). This suggests that the 100 ka cycle seen in the
ice volume records is dominated by deglaciations that occur
every 100 ka after the MPR, on average. Note that summer
insolation at 65
N is also strongly influenced by precession;
hence, the saw-toothed pattern of rapid deglaciation occurs
every fourth or fifth precessional cycle over the last 600 ka.
A similar saw-tooth pattern can be produced from obliquity
using every second or third cycle (Ridgwell et al., 1999;
Huybers and Wunsch, 2005; Raymo et al., 2006; Raymo and
Huybers, 2008), but the spectral analysis does not match the
global ice volume signal as well (Ridgwell et al., 1999; Maslin
and Ridgwell, 2005). Models have also been used to re-create
the 100 ka cycle and the most successful is the LLN 2D model,
QUATERNARY CLIMATE TRANSITIONS AND CYCLES 847