A new interpretation of the two-step δ 18 O signal at the Eocene-Oligocene boundary

The most marked step in the global climate transition from “Greenhouse” to “Icehouse” Earth occurred at the Eocene-Oligocene (E − O) boundary, 33.7 Ma. Evidence for climatic changes comes from many sources, including the marine benthic δ 18 O record, showing an increase by 1.2–1.5‰ at this time. This positive excursion is characterised 5 by two steps, separated by a plateau. The increase in δ 18 O values has been attributed to rapid glaciation of the Antarctic continent, previously ice-free. Simultaneous changes in the δ 13 C record are indicative of a greenhouse gas control on climate. Previous studies show that a decline in p CO 2 beyond a certain threshold value may have initiated the growth of a Southern Hemispheric ice sheet. These studies were not able to con-10 clusively explain the remarkable two-step proﬁle in δ 18 O. Furthermore, they did not address the potential role of changes in ocean circulation in the E − O transition. Here a new interpretation of the δ 18 O signal is presented, based on model simulations using a simple coupled 8-box-ocean, 4-box-atmosphere model with an added land ice component. The model was forced with a slowly decreasing atmospheric carbon dioxide 15 concentration. It is argued that the ﬁrst step in the δ 18 O represents a shift in meridional overturning circulation from a Southern Ocean to a bipolar source of deep-water formation, which is associated with a cooling of the deep sea. This shift can be initiated by a small density perturbation in the model, although there is also a parameter regime for which the shift occurs spontaneously. The second step in the δ 18 O proﬁle occurs

Era, Earth's climate experienced a transition from a "Greenhouse" world to the current "Icehouse" world, the most marked step of which occurred at the Eocene-Oligocene (E−O) boundary, 33.7 Ma.The Eocene climate was relatively warm, with polar surface and deep water temperatures up to 10 • C warmer than today.The Antarctic continent -though in its present position -was ice-free and lushly vegetated.At the E−O boundary, global temperatures dropped and a semi-permanent ice-sheet formed on Antarctica.An overwhelming amount of evidence from fossils and proxies indicates that climate change at the E-O boundary was global and involved all parts of the climate system (Lear et al., 2000;DeConto and Pollard, 2003a;Dockery and Lozouet, 2003;Hay et al., 2005;Zachos and Kump, 2005;Coxall and Pearson, 2007;Liu et al., 2009).
The most pronounced evidence for relatively rapid climatic changes comes from the marine δ 18 O record.Numerous cores taken from the ocean floor display an abrupt increase in δ 18 O of 1.2-1.5‰across the E−O boundary.The shift lasts approximately 500 kyr and is characterised by a remarkable two-step profile in which two 40-kyr steps are separated by a plateau of 200 kyr.The E−O transition terminates with a sustained maximum persisting for 400 kyr, followed by several stepped decreases.Eventually δ 18 O values stabilise ∼31 Ma to a value ∼1‰ higher than before the transition (Zachos et al., 1996;Coxall et al., 2005;Zachos and Kump, 2005;Coxall and Pearson, 2007).The rapid shift in δ 18 O towards a sustained maximum has been associated with a rapid and extensive glaciation of at least the Antarctic continent.Glaciomarine sediments in the vicinity of the Antarctic continent support this interpretation.Although Northern Hemispheric ice rafted debris dating back to the middle Eocene has also been found, data are insufficient to determine the extent of a Northern Hemisphere glaciation at this time (Coxall and Pearson, 2007, and references therein).It is highly likely that a rapid glaciation of the Antarctic continent was accompanied by a (local) decrease in temperature and there are many indications that the δ 18 O signal comprises both an ice volume and a temperature change (Oerlemans, 2004;Coxall and Pearson, 2007).Analysis of the Mg/Ca-ratio in benthic foraminifera indicates a cooling of up to Figures 2.5 • C in the tropics (Lear et al., 2000(Lear et al., , 2004(Lear et al., , 2008)).U K 37 and TEX 86 analyses have lead to the conclusion that high-latitude sea surface temperatures (SSTs) decreased by 4.8-5.4 • C from the late Eocene to the early and mid-Oligocene.Model simulations show that such a decrease in SSTs corresponds to a decrease in deep-water temperatures of 3-5 • C. Combining this result with the δ 18 O record produces maximum ice volume estimates of 40-120% of modern Antarctic ice volume (Liu et al., 2009).Changes in the global carbon cycle are reflected in the δ 13 C record, which displays a rapid and stepwise increase of 1‰ similar in shape and duration to the shift observed in δ 18 O values.The increase in δ 13 C is thought to be indicative of enhanced marine export production (Zachos and Kump, 2005;Coxall and Pearson, 2007).The sharp rise of carbon isotope ratios has also been associated with a two-step deepening of the carbon compensation depth (CCD) by more than 1 km (Rea and Lyle, 2005;Zachos and Kump, 2005).The CCD is linked to ocean acidity, which is in turn linked to atmospheric CO 2 concentration (Coxall et al., 2005;Merico et al., 2008).Proxy records show a general decline in pCO 2 throughout the Cenozoic from 2-5 times pre-industrial atmospheric levels (PAL) in the mid-Eocene, to near-modern values in the lower Miocene, although sufficient detail to accurately determine the CO 2 decline across the E−O boundary is lacking (Pagani et al., 2005;Coxall and Pearson, 2007;Pearson et al., 2009).
During the Eocene, the Southern Ocean and even mid-latitudes are thought to have played a greater role in deep-water formation.Sedimentary and seismic evidence in combination with analysis of neodymium (Nd) isotopes shows that there was a transition from a southern sinking state to a bipolar source of deep-water formation at the beginning of the Oligocene (Thomas et al., 2003;Thomas, 2004;Via and Thomas, 2006;Coxall and Pearson, 2007).Other implied climatic changes include an increased Figures

Back Close
Full Precise timing of the opening of these gateways has, however, proven to be problematic (Livermore et al., 2005;Scher and Martin, 2006;Coxall and Pearson, 2007).
Evidence to support a vigourous ACC before the mid-Miocene is scarce (Coxall and Pearson, 2007).Furthermore, DeConto and Pollard (2003a) argue that the influence of colder Southern Ocean SSTs on Antarctic meteorology is poorly resolved, and Huber and Nof (2006) demonstrate that ocean heat transport would have had to decrease enormously for any significant changes in Antarctic continental temperatures to occur.Huber and Nof (2006) argue additionally, that the major changes in productivity during the E−O transition are strong evidence in favour of a greenhouse gas control on climate.It has been proposed that, once some CO 2 threshold (∼750 ppm, Pearson et al., 2009) was reached, feedbacks related to snow/ice-albedo and ice-sheet height/massbalance could have initiated rapid ice-sheet growth (DeConto and Pollard, 2003a,b;Coxall and Pearson, 2007).DeConto and Pollard (2003a) ran a global climate model coupled to a dynamic ice-sheet model, forcing it with a CO 2 decline from 4 to 2 times PAL in 10 Myr.They found that relatively small ice caps form due to high levels of winter precipitation, which start to grow rapidly beyond some threshold CO 2 concentration and eventually coalesce to form one continental scale ice-sheet.In their simulations with a non-dynamical, 50-m slab ocean, the opening of Drake Passage only affects the timing of ice-sheet growth.It is thought that the initiation and rapid growth of the ice-sheet is aided by the specific orbital settings at the E−O boundary: low eccentricity and low-amplitude changes in obliquity favour dampened seasonality and cold summers (DeConto and Pollard, 2003a,b;Coxall et al., 2005;Coxall and Pearson, 2007) Thus far, model studies have largely focussed on the role of individual parts of the climate system.Since data show that climate change at the E−O boundary was global and involved all parts of the climate system, it is to be expected that interactions and feedbacks between the different climate system components are important.In addition, previous model studies have failed to address the potential role of changes in the global meridional ocean circulation (MOC) in the E−O climate transition.
Here a new interpretation of the δ 18 O record will be presented, in which several parts of the climate system are involved and switches in the MOC play a dominant role.The transition from "Greenhouse" to "Icehouse" is simulated using an adapted version of the model developed by Gildor andTziperman (2000, 2001) and Gildor et al. (2002)

The model
The model that is used is an adapted version of the coupled ocean-atmosphere model with added land ice and sea ice components that was created by Gildor andTziperman (2000, 2001) to explain the 100-kyr period observed in glacial-interglacial oscillations for the past 400 kyr.Added to this model are the biogeochemistry component developed by Gildor et al. (2002) and a module that we developed to compute changes in Figures In this section the most important characteristics of the model are discussed; a full description of the equations can be found in Gildor and Tziperman (2001) and Gildor et al. (2002).An overview of the parameters that were altered with respect to Gildor et al. (2002) is presented in Table 1.
The ocean is represented by four surface boxes and four deep boxes whose latitudinal extends are South Pole to 45 • S, 45 • S to equator, equator to 45 • N and 45 • N to North Pole.The meridional overturning circulation is buoyancy driven; advection and diffusion of temperature and salinity are balanced by surface fluxes from the atmosphere and land ice components.The atmospheric model consists of four vertically averaged boxes representing the same latitude bands as the ocean boxes.The surface below the atmospheric boxes is a combination of ocean, land, land ice and sea ice, each with its specified albedo.The averaged potential temperature in the box is calculated on the basis of the energy balance in the box, consisting of: incoming solar radiation, using box albedo and with the possibility to include Milankovitch variations; outgoing long-wave radiation; heat-flux into the ocean; and meridional atmospheric heat transport.The meridional atmospheric moisture transport is proportional to the meridional temperature gradient and the humidity of the box to which the flux is directed.
The land ice model assumes a simple, zonally symmetric ice sheet of perfect plasticity.The ice sheet grows due to precipitation in the box, which is assumed to be distributed evenly.It decreases as a result of ablation, evaporation and calving processes.This sink term is taken to be constant in time, as it is expected that the source term dependence on temperature dominates that of the sink term.A known problem of many models, including this one, is a too high sensitivity to the sea ice-albedo feed-Introduction

Conclusions References
Tables Figures

Back Close
Full in the Eocene-Oligocene transition (Huber and Sloan, 2001) -the sea ice component is left out of the study entirely to allow land fractions to take on realistic values.Hence, there will be no sea ice present in both polar boxes.
An additional module was added to the model that computes variations in δ 18 O.
The oxygen isotopic composition of sea-water δ w is calculated using the volumetric changes in land ice and ocean and conservation of total δ 18 O: where V w 0 is the volume of the Eocene ocean, and V w is the volume of the ocean after ice formation: V w =V w 0 −V ice .δ ice is the isotopic composition of land ice and δ w 0 and δ w are the isotopic compositions of the Eocene ocean and the ocean after ice formation, respectively.The values that were used for these variables can be found in Table 1.
The sea-water isotopic composition δ w together with the simulated temperature T (in • C) of the four lower ocean boxes is then used to calculate the calcite isotopic composition δ c according to the calcite-temperature relation of Erez and Luz (1983): The prognostic equations in the model are solved using a leapfrog time scheme with a uniform time step of 6 h.All variables are averaged over 200 years (Gildor and Tziperman, 2001).For the present study biogeochemical feedback processes are ignored.

The Eocene reference state
First the model was adapted to paleoclimatic boundary conditions representing those of the Eocene climate system.The largest differences between the climatic boundary conditions for which the model was developed (Gildor andTziperman, 2000, 2001; Introduction

Conclusions References
Tables Figures

Back Close
Full Gildor et al., 2002) and those of the Eocene climate, are the position of the continents and the background atmospheric greenhouse gas concentrations.Therefore in the model the land fraction per box and pCO 2 were adjusted.Representative values for the land fractions were obtained from Markwick et al. (2000) and can be found in Table 1.A CO 2 concentration was sought for which the land ice would disappear (almost) entirely, while not (excessively) exceeding the upper-limit of 5 times PAL as set by Pagani et al. (2005).
An overview of our modelled state of the Eocene climate system at a pCO 2 of 1500 ppm is presented in Fig. 2. At this concentration, the Southern Hemispheric ice sheet has disappeared entirely, and the ice sheet in the northern box is non-perennial and a factor 10 4 smaller than it is presently (Williams and Ferrigno, 2005), see Fig. 2f.A southern sinking pattern of the MOC (Fig. 2d) was found to be a stable solution of the model for these Eocene boundary conditions.Such a state is in agreement with proxy data and model reconstructions (Thomas et al., 2003;Thomas, 2004;Via and Thomas, 2006;Coxall and Pearson, 2007).
The yearly averaged atmospheric temperature (Fig. 2a) is slightly over 30 • C in the equatorial boxes, around 13 • C in the southern box and just below freezing in the northern box.Subsequently, the atmospheric meridional temperature gradient (Fig. 2b) in the Northern Hemisphere is quite high, while taking on a more moderate value in the South.A similar pattern can be observed in the yearly averaged surface ocean temperatures (Fig. 2c).The oceanic meridional temperature gradients are therefore comparable to the atmospheric ones.In the polar boxes the seasonality is much higher (10-15 • C) than in the equatorial boxes ( 5• C).The yearly averaged freshwater flux into the ocean (Fig. 2e) is positive in the polar boxes and negative in the equatorial boxes; i.e., there is evaporation at the equator and precipitation at the poles.The seasonality of the freshwater flux is much higher at the equator than at the poles.
To check the validity of this reference state, the ocean temperature results are compared with Eocene data from Bijl et al. (2009) and Liu et al. (2009), and GCM model results from Huber and Sloan (2001) Oberhuber (1988).These are of the same order of magnitude, so the model determines realistic values of the buoyancy flux.

Critical threshold for ice formation
In order to test the hypothesis of ice sheet inception after lowering of atmospheric greenhouse gas concentrations below a certain threshold value, the model was forced with linearly decreasing CO 2 concentrations.Following the measurements of Pagani et al. (2005), pCO 2 was chosen to decrease at a rate of 750 ppm per 200 kyr from 1500 ppm to 100 ppm.The carbon dioxide concentration and the corresponding Southern Hemispheric ice volume are plotted in Fig. 3a and b.The threshold CO 2 value for Antarctic ice growth is 270 ppm, at which the ice sheet starts to grow exponentially and reaches its full size within 100 kyr.In the northern polar box, no significant changes in ice sheet volume occur.The response of the climate system to decreasing carbon dioxide concentrations differs per box.At first, when no ice sheet is formed yet, the temperature decrease as plotted in Fig. 3c shows a similar pattern across the globe: slowly at first, but more rapidly as pCO 2 continues to fall.The inception of the Southern Hemispheric ice sheet lowers atmospheric temperatures in the most southern polar box and the southern equatorial box additionally through the ice-albedo feedback, while northern temperatures are barely affected.During the transition, the northern meridional temperature difference decreases, while the southern difference increases, especially after the ice sheet has started to grow.The amount of equatorial evaporation and polar precipitation decreases, see Fig.It was also tested whether the climate system returns to its initial state when pCO 2 is raised back to 1500 ppm.The Antarctic ice sheet starts to melt at 605 ppm, a value well above the CO 2 concentration of ice inception.The difference in critical pCO 2 levels is due to the ice-albedo feedback which lowers the atmospheric temperature as long as the ice sheet is present.It takes ∼300 kyr for the ice sheet to melt entirely, and this is a linear rather than an exponential process.
Changes in temperature and freshwater flux in the northern boxes under increasing CO 2 are a mirrored version of the response to decreasing CO 2 described above: temperatures increase at a decreasing rate.The response in the southern boxes follows the melting of the Antarctic ice sheet and is therefore linear.Once the CO 2 concentration has returned to 1500 ppm and the climate system is restored to equilibrium, the global temperature and freshwater distributions are exactly the same as they were before the transition.Throughout the transition the southern sinking state of the MOC remains unaffected.When the rate of pCO 2 decrease is doubled or halved the results remain qualitatively the same.

Existence of multiple climate equilibria
Differential heating of Earth's surface leads to cool, fresh polar surface waters and warm, salty equatorial surface waters.Temperature and salinity therefore have opposing roles in determining the location of deep water formation.This suggests that multiple MOC patterns are possible.This elementary problem was studied by Stommel (1961) using a density driven two-box model.It was shown that in this model indeed two different states of the MOC exist.In an extended version of this model, with three boxes (Welander, 1986;Thual and McWilliams, 1992), four solutions exist: sinking in the North and upwelling in the South (NPP), sinking in the South and upwelling in the North (SPP), sinking at both poles (TH), and upwelling at both poles (SA).Whether the solutions are stable depends on the boundary conditions (i.e.parameter settings) of the model.Since our ocean model is similar to the models described here, it is expected that multiple steady state solutions might exist in our Eocene climate.Indeed, 1401 Introduction

Conclusions References
Tables Figures

Back Close
Full for the Eocene boundary conditions a TH state (Fig. 4d) and NPP state (Fig. 5d) are also found to be stable solutions of the model, in addition to the SPP state found earlier (Fig. 2d).
The state of the MOC affects the global temperature distribution: compare Fig. 2 to Figs. 4 and 5.In the NPP and TH states there is more oceanic heat transport towards the North than in the SPP state; consequently, atmospheric and sea surface temperatures there are higher, leading to the disappearance of the Arctic ice sheet.Northerly SSTs are in these states more in accordance with the Bijl et al. (2009) and Liu et al. (2009) data and Huber and Sloan (2001) GCM results.Similarly, southerly atmospheric and sea surface temperatures in the NPP state are lower as a result of a decrease in heat transport.This lowering is insufficient to initiate the growth of an Antarctic ice sheet.Deep sea temperatures are coldest in the TH state (13.6 • C, Fig. 4c) and warmest in the NPP state (17.2 • C, Fig. 5c).
The simulations with changing pCO 2 (Sect.3.2) showed that even if global atmospheric temperatures change drastically, the SPP state of the Eocene reference state persists.To investigate the stability of all the MOC states under decreasing CO 2 concentrations, the model was run for each state with pCO 2 decreasing with 100 ppm increments.For every CO 2 concentration the model was run until the system was in equilibrium.The pattern and strength of the MOC can be deduced from the value of q N −q S , where q N is the volume transport from the surface layer to the deep ocean in the northern box, and q S in the southern box, see Fig. 2d.A q N −q S that is positive indicates a NPP state, a negative q N −q S indicates a SPP state and if q N −q S is approximately zero the ocean is in a TH state.The resulting bifurcation diagram is plotted in Fig. 6a and shows that the SPP and the TH state of the overturning circulation coexist for all values of pCO 2 .This is not the case for the NPP state, which ceases to exist for CO 2 concentrations lower than 1400 ppm.At these concentrations the NPP state can no longer be sustained due to a decrease in the freshwater flux in the northern polar box.Introduction

Conclusions References
Tables Figures

Back Close
Full Ice formation on the Southern Hemisphere in the TH state occurs at a CO 2 concentration of 300 ppm; for the SPP state ice formation starts at 200 ppm.For concentrations lower than 100 ppm oceanic temperatures become unphysical, causing the overturning circulation to collapse into the SA state.

Transitions between equilibria
The fact that the SPP and TH states coexist for the entire range of CO 2 concentrations between 100 and 1500 ppm indicates that it is possible for the climate system to switch between these states.Such a shift would be in accordance with data that show that the MOC at the E−O boundary changed from a Southern Ocean to a bipolar source of deep-water formation (Thomas et al., 2003;Thomas, 2004;Via and Thomas, 2006;Coxall and Pearson, 2007).Therefore it is tested whether it is possible to initiate such a transition in the model.
A transition from the SPP state to the TH state is induced at a CO 2 concentration of 800 ppm by adding a density perturbation to the surface layer in the northern polar box.This perturbation lasts for a short time (200 year, i.e. one averaging time-step) and is of small magnitude (1 kg m −3 ) compared to the background density.The immediate result is that water in this box starts to sink.Subsequently, heavier water dominating the artificial perturbation is advected and causes the MOC to change to the TH state.When the CO 2 concentration is lowered further, the system remains on the TH branch (Fig. 6a).The ability of the system to switch between MOC states exists for all CO 2 concentrations.Transitions between MOC states as a result of a finite amplitude density perturbation can always occur due to noise in the climate system.
The corresponding δ 18 O profile (Fig. 6b) consists of two rapid increases, separated temperature decrease of 2.3 • C.This is in accordance with the temperature decrease found by Liu et al. (2009).The second step represents Southern Hemispheric ice growth.The length of the plateau depends on the timing of the density perturbation with respect to the critical value of pCO 2 for Antarctic ice formation.As long as the pCO 2 of the MOC transition is higher than the pCO 2 of Southern Hemispheric ice inception, the δ 18 O profile associated with a decreasing CO 2 concentration will always display a two-step signature.The magnitude of the second δ 18 O shift is 3.5‰, of which 2.6‰ is due to ice growth.Since data indicate that the total increase in δ 18 O at the E−O boundary was of the order of 1.5‰, the magnitude of this step is too large.Inception of ice growth on the Southern Hemisphere in our model occurs at CO 2 concentrations below 400 ppm.This is very low compared to the measurements of Pearson et al. (2009) and model studies by DeConto and Pollard (2003a,b), who estimated Antarctic ice growth to start at 750 ppm.A possible explanation for the relatively low CO 2 concentration of ice inception in the model is that the ice sheet only grows if the temperature in the polar box rather than at the pole is below 0 • C. Since the box temperature is an average, it will already be far below freezing at the pole when this threshold is reached.Therefore, the pCO 2 of ice inception might be raised by increasing the inception temperature of ice in the model.The first test that was done involved raising the inception temperature of ice to 3 • C. Ice growth now occurs for CO 2 concentrations lower than 800 ppm.It is again possible to initiate a transition from the SPP to the TH state by perturbing the density.This can be done at any pCO 2 ≥800 ppm; here we chose to place the perturbation at 1000 ppm.The resulting bifurcation diagram is shown in Fig. 7a.Again, the δ 18 O shift (Fig. 7b) exhibits a two-step profile.Both steps are comparable in magnitude to the steps in Fig. 6b.These results show that it is indeed possible to improve the pCO 2 of ice inception in the model.When the inception temperature is raised even higher, to 5 a CO 2 concentration lower than 1500 ppm and higher than the pCO 2 of Antarctic ice inception.For this inception temperature, the transition from the SPP to the TH stateand hence the two-step δ 18 O profile -will be a spontaneous outcome of the model.

Summary and discussion
Global cooling and the initiation of the Antarctic ice sheet at the E−O boundary 33.Gildor andTziperman (2000, 2001) and Gildor et al. (2002) was used to simulate the response of the full climate system to a decrease in pCO 2 .First the Eocene reference state was found by changing the land fractions in the model and raising the CO 2 concentration to 1500 ppm.The MOC displayed a southern sinking (SPP) state, in accordance with proxy data.Under these boundary conditions, a northern sinking (NPP) and bipolar sinking (TH) state are also stable solutions.When the model is forced with decreasing pCO 2 it naturally becomes colder.At low carbon dioxide concentrations (270 ppm) an Antarctic ice sheet is formed.The SPP state of the MOC remains unaffected by the climatic changes.Therefore it was tested whether all states of the MOC exist for all CO 2 concentrations lower than 1500 ppm.
The NPP state ceases to exist for CO 2 concentrations lower than 1300 ppm.Both the SPP and the TH pattern exist for all levels of pCO 2 between 100 and 1500 ppm.This suggests that it is possible for the system to switch between these MOC states.
A transition from the SPP to the TH state was triggered by means of a small density perturbation.The switch in overturning state leads to a cooling of the deep sea and Introduction

Conclusions References
Tables Figures

Back Close
Full The CO 2 concentration of ice inception was increased by increasing the inception temperature of ice in the model.At an inception temperature of 3 • C ice growth is initiated at 800 ppm.This is in good agreement with data and modelling studies.For an ice inception temperature of 5 • C, the Eocene reference state with the MOC in a SPP pattern ceases to exist.This suggests that an inception temperature between 3 and 5 • C exists for which the transition from the SPP to the TH state under decreasing pCO 2 is a spontaneous outcome of the model.Despite the fact that a model of very low complexity was used to simulate the Eocene-Oligocene transition, our new interpretation of the δ 18 O record appears to be coherent and complies qualitatively well with proxy data.The exact timing of the modelled shifts and intermittent plateau is still undetermined and is largely dependent on the exact pattern of pCO 2 change.Quantitatively the first modelled step -representing a transition in overturning circulation -also agrees well with the δ 18 O proxy records.
Only the magnitude of the second step decreases the quantitative resemblance between model results and data: the modelled shift is too large.
In our model simulations, the Southern Hemispheric ice sheet was allowed to grow to its full size, 9×10 16 m 3 .This is three times larger than the present-day size of the Antarctic ice sheet (Williams and Ferrigno, 2005) and therefore much higher than the maximum E−O ice volume estimates (40-120% of modern Antarctic ice volume, Liu et al., 2009).It takes the ice sheet slightly over 100 kyr to grow to this size, while data indicate that the second step only lasted for 40 kyr.Indeed, after 40 kyr, the ice sheet in our model has grown to 9×10 15 m 3 , 30% of modern Antarctic ice volume.The δ 18 O excursion would then be 0.67‰, of which 0.24‰ is due to ice growth.The magnitude of the two shifts combined would now be 1.17‰, agreeing well with the observed increase of 1.2-1.5‰.Therefore there must have been a mechanism that stopped the growth of the Antarctic ice sheet after 40 kyr.This mechanism might be related to the CO 2 record, which is known to have displayed variable rates of changes around the E−O Introduction

Conclusions References
Tables Figures

Back Close
Full boundary.According to measurements by Pearson et al. (2009), there might even have been a temporary increase in CO 2 concentrations during the second step, which could have stabilised the growth of the ice sheet.Another plausible mechanism is phase locking with the obliquity.Northern polar temperatures in our modelled Eocene reference state were found to be too cold, and deep sea temperatures were found to be too warm.These discrepancies between modelled temperatures and data are a shortcoming of the model.In our box model, the surface flow from the North Pole to the equator is in one direction (southward) only, lowering northern polar temperatures significantly.In reality, a cold southward surface current would be (partially) matched by a warm northward current, as with the Gulf Stream today.This would decrease the total heat transport and thus moderate northern polar temperatures.Too high deep sea temperatures are a known problem of ocean box models.
Although changes in MOC and ice volume under decreasing CO 2 concentrations prove to be a plausible explanation for the observed δ 18 O signal, possible mechanisms behind the pCO 2 decline remain unaddressed.In our model study atmospheric pCO 2 was a prescribed forcing and not an output variable.Proposed causes of the CO 2 decline in the literature are reduced sea-floor spreading, increased silicate weathering, increased marine production and release of organic carbon from geological reservoirs (Zachos and Kump, 2005;Coxall and Pearson, 2007;Pearson et al., 2009).Positive feedback mechanisms exist between these processes and the changes in MOC and ice volume found in our simulations.The role and strength of these mechanisms could be tested with the biogeochemical module developed by Gildor et al. (2002) or the more comprehensive models of Zachos and Kump (2005) and Merico et al. (2008), but are outside the scope of this paper.
In this model study only a limited number of feedback processes has been incorporated.These processes have proven to be sufficient to newly interpret the δ 18 O record at the E−O boundary.Nonetheless, the results of this research will be even more robust when they can be reproduced using more complex models.The role of several Introduction

Conclusions References
Tables Figures

Back Close
Full potentially important processes should be investigated in the future.To fully understand the influence of the opening of oceanic gateways, a model should be used in which Southern Ocean dynamics are included and an ACC can be formed.It will also be interesting to see how the changes described here affect atmospheric circulation patterns, and vice versa.Although the land ice module used in this study is very simple, its results are in agreement with the research by DeConto and Pollard (2003a,b).
The precise details of Antarctic ice sheet formation can be better understood from their simulations.Changes in orbital configuration were not taken into account in our simulations, even though they are thought to have played a significant role.Full atmospheric heat transport coefficient 15.5×10 −20 13.5×10 Full long-term variations in our present-day climate are paced by the waxing and waning of the Northern Hemispheric ice sheet, Earth's surface was (almost) entirely ice-free at the beginning of the Cenozoic Era, 65.5 Ma.In the course of this Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | declining CO 2 values below a certain threshold value as the main driver.The opening of Drake Passage between Antarctica and South America and of Tasmanian Passage between Antarctica and Australia has facilitated the organisation of the Antarctic Circumpolar Current (ACC), reducing southward oceanic heat transport and cooling Southern Ocean SSTs.Since the opening of these passages occurred around the E−O boundary, it might have played a dominant role in the glaciation of Antarctica.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | δ 18 O.By changing the parameters of the model, it can be adapted to paleoclimatic boundary conditions (see below).
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | by a plateau.The first of these shifts represents the transition in the MOC from the SPP to the TH state.Deep sea temperatures in the TH state are colder than in the SPP state (compare Fig. 4c to Fig. 2c), because deep water is formed in both polar regions instead of one.At colder temperatures, the calcite δ 18 O value (δ c ) is higher.During the shift, the modelled increase in δ 18 O is 0.5‰, which corresponds to a deep sea Discussion Paper | Discussion Paper | Discussion Paper | • C, the Eocene reference state with the MOC in the SPP pattern ceases to exist.Therefore there must exist an inception temperature between 3 and 5 • C for which the SPP state ceases to exist at Discussion Paper | Discussion Paper | Discussion Paper |

7
Ma are thought to be the result of decreasing atmospheric carbon dioxide levels below a certain threshold value.Thus far, the remarkable two-step profile of the δ 18 O record at this boundary has remained unexplained.Proxy data indicate a change in the meridional overturning circulation (MOC) from a southern sinking to a bipolar sinking state.The role of these changes during the E−O climate transition has not yet been investigated.Therefore an adapted version of the simple coupled 8-box-ocean, 4-boxatmosphere model with added land ice component by Discussion Paper | Discussion Paper | Discussion Paper | hence a positive excursion of δ 18 O values.The second step in the δ 18 O record can be interpreted as representing the rapid glaciation of the Antarctic ice sheet.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |
. In the southern and equatorial boxes, modelled Figures SSTs are in reasonable agreement with measurements.However, proxy data as well as the GCM results indicate that both poles were relatively warm (up to 20 • C) during the Eocene.Simulated northern box temperatures are therefore too cold.Liu et al. (2009) measured deep-sea temperatures of ∼10 • C. Deep-sea temperatures in our modelled reference state are therefore too warm.Values for the freshwater flux into the ocean were compared against modern-day observations by Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | This model study has provided us with a new interpretation of the two-step δ 18 O profile at the E−O boundary.It was found that the first step is the result of a switch in meridional overturning circulation from a southern to a bipolar sinking state.This transition leads to a cooling of the deep sea and thus to an increase in marine benthic δ 18 O values.The second step represents rapid Antarctic ice growth, caused by lowering of the atmospheric carbon dioxide levels below a certain threshold level.The two-step signature will always occur when the MOC transition -which can occur spontaneously due to random density perturbations -occurs at higher pCO 2 levels than the critical pCO 2 level for Southern Hemispheric ice formation.The mechanisms that are proposed here are in good agreement with proxy data, and so is the order of magnitude of the two modelled δ 18 O steps.Introduction