Constraining the Last Glacial Maximum climate by data-model ( i LOVECLIM ) comparison using oxygen stable isotopes

Introduction Conclusions References


Introduction
Oxygen stable isotopes (δ 18 O) are among the most widely used/common tools in palaeoclimatology-palaeoceanography.δ 18 O constitutes an important tracer of the hydrological cycle for the different components of the climatic system (ocean, atmosphere, ice sheets) but processes that control the recorded δ 18 O signal are various and complex.Simulation of climate and its associated isotopic signal allow an investigation of these various and complex processes (e.g.Roche et al., 2004a;Lewis et al., 2010).
Water isotopes have been implemented in numerous atmospheric and oceanic general circulation models (Jouzel et al., 1987;Joussame and Jouzel, 1993;Hoffmann et al., 1998;Schmidt, 1998;Paul et al., 1999;Delaygue et al., 2000;Werner et al., 2000Werner et al., , 2011;;Noone and Simmonds, 2002;Mathieu et al., 2002;Lee et al., 2007;Yoshimura et al., 2008;Zhou et al., 2008;Tindall et al., 2009;Risi et al., 2010;Xu et al., 2012).However, they have very seldom been used in coupled climate simulations with water isotopes in both the atmospheric and oceanic components, due to computational costs (Schmidt et al., 2007;LeGrande and Schmidt, 2011).We chose here to use an intermediate complexity model to circumvent that limitation, while retaining a full oceanic general circulation model to allow investigating the details of the oceanic response where numerous palaeodata are available.i LOVECLIM (Roche, 2013;Roche and Caley, 2013;Caley and Roche, 2013).In the present study, we focus on the Last Glacial Maximum (LGM) climate, a cold climatic extreme.We will refer to the LGM as the period around 21 thousand years before present because (1) it is coherent with the commonly used time period in previous climate simulations and ( 2) it has been defined as the period between 23 and 19 thousand years before present based on a thorough examination of palaeoclimatic data (Mix et al., 2001).To our knowledge, the present study is the first attempt at a global atmosphere -ocean coupled simulation of the Last Glacial Maximum with an oceanic general circulation model.One earlier coupled model study including oxygen isotopes used a spatially greatly simplified model (Roche et al., 2004b).
The LGM has become a standard target period to constrain climate sensitivity and evaluate model's capability to simulate a climate that is drastically different from that of the present-day.Efforts have been made to improve multi-model intercomparisons under glacial boundary conditions to diagnose the range of possible responses, to a given change in forcing.These efforts have been realized by the pioneering work of the Palaeoclimate Modeling Intercomparison Project (PMIP).PMIP used Atmospheric General Circulation models (Joussaume and Taylor, 2000) in its first phase, and Coupled Atmospheric-Ocean models in its second phase (PMIP2) (Crucifix et al., 2005;Braconnot et al., 2007a, b).PMIP is now in its third phase (CMIP5/PMIP3 palaeo-simulations) and uses a coherent framework between past, present and future simulations to provide a systematic and quantitative assessment of the realism of the models used to predict the future (Braconnot, 2012;Schmidt et al., 2013).
The evaluations of model results are supported by environmental proxy data that have been compiled in several projects aimed at the reconstruction of the LGM surface conditions, like Climate Long-Range Investigation Mapping and Prediction (CLIMAP, 1981), Glacial Atlantic Mapping and Prediction (GLAMAP2000, Sarnthein et al., 2003) and more recently The Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO, Kucera et al., 2005;MARGO Project Members, 2009).However, all environmental proxies are influenced by non-climatic factors and/or Introduction

Conclusions References
Tables Figures

Back Close
Full further climatic factors.A promising way of exploiting palaeo-environmental data for climate-model evaluation is to use models that explicitly simulate the sensor.A decisive advantage of isotope-enabled models is their ability to directly simulate measured quantities (in the present case, δ 18 O), so that their results can be directly compared with data from the different climatic archives.
In this study, we present a comparison between LGM simulated and measured oxygen isotopes.The results are presented in terms of anomaly between the LGM and the present.The use of anomaly renders absolute values irrelevant; it concerns a purely relative change.We can therefore ignore complications such as speciesspecific climate variable relationships, vital effect offsets, and calibrations of values measured relative to the VPDB standard to values on the VSMOW scale (Rohling and Cooke, 1999).We consider both the atmospheric and oceanic component and discuss the agreement between data and model results, as well as the processes driving the isotopic signal recorded.

LGM boundary conditions
We use the boundary conditions defined in/by the PMIP2 protocol to simulate the LGM climate.Lowered levels of atmospheric greenhouse gas concentrations (CO 2 = 185 ppm, CH 4 = 350 ppb and NO 2 = 200 ppb) are used in agreement with icecore measurements (Fluckiger et al., 1999;Dällenbach et al., 2000;Monnin et al., 2001).Ice-sheet topography changes are taken from Peltier (2004) and the surface albedo is set accordingly.Orbital parameters correspond to 21 000 yr before present (Berger and Loutre, 1992).To account for the ∼ 130 m decrease in sea level relative to present day, the land-sea mask and the oceanic bathymetry are modified (Lambeck and Chappell, 2001).Some variations exist among the PMIP simulations, mainly for the Northern Hemisphere, in the handling of changes in the river basins (Weber et al., Introduction Conclusions References Tables Figures

Back Close
Full 2007), i.e. changes in river routing due to the presence of ice-sheets.In our LGM simulation we included changes in the water routing from the Laurentide ice-sheet over North America and from the Fennoscandian ice sheet over Eurasia.We applied the above forcings to the model's equilibrium state and integrated it until a new equilibrium was reached, after 5000 yr of integration.
Our choice of using the PMIP2 boundary conditions instead of the more recent PMIP3 protocol arise from several considerations: (1) having a state readily comparable to the already published LGM state of an earlier version of the model (Roche et al., 2007), (2) possibility in a future study to intercompare our atmospheric results to already published PMIP2 LGM atmospheric general circulation model, (3) the fact that the main difference between the PMIP2 and PMIP3 protocols is in the ice-sheet elevation (cf.PMIP, 2013) which, interpolated on our atmospheric coarse resolution grid, is not very large.

Global datasets compilation
Global oxygen isotopic datasets for the atmospheric and oceanic components have been already compiled for the Late Holocene (LH) and have been compared and discussed with i LOVECLIM results (Caley and Roche, 2013) (Rozanski et al., 1993).S1 and S2).Working with anomalies instead of LGM and LH values would reduce uncertainties related to potential changes in planktic foraminifera seasonality or depth habitat (Waelbroeck et al., 2005).We selected marine cores for which the chronology is essentially based on radiogenic dating ( 14 C).Chronostratigraphic quality has been defined following the MARGO project definition (Kucera et al., 2005) (Tables S1 and S2).
Concerning planktic foraminifera, we worked with δ 18 O data of the most commonly measured planktic foraminifer species: Globigerinoides ruber white and pink, Globigerinoides sacculifer, Globigerina bulloides, and Neogloboquadrina pachyderma sinistral for the majority of the sites.Two exceptions occur at sites MD96-2077 and OCE326-GGC5 for which Globorotalia inflata is used (Table S1).We extended our data set with 10 calcite δ 18 O measurements of planktic foraminifera Neogloboquadrina pachyderma sinistral compiled by Meland et al. (2005).We only selected records for which the LGM level was determined by AMS 14 C and for which LH data (last 3000 yr) was available.These authors used a definition for the LGM chronozone slightly different (18-21.5 ka) but the data allow comparison with model results in the Nordic seas, an area poorly documented in our compilation (Table S1).
Concerning benthic foraminifera δ 18 O, we extended our data set with 22 calcite

Oxygen stable isotopes in precipitation
In the following we only briefly discuss the annual mean distribution of δ 18 O in precipitation, since a complete description and intercomparison with previously published LGM δ 18 O in precipitation simulations will be undertaken in another study/article in preparation.
Isotopic values are expressed relative to a standard (Sharp, 2007): where R is the abundance ratio of the heavy and light isotopes (e.g., N( 18 O) p /N( 16 O) p ) for substance P and δ is commonly reported in units of parts per thousand (‰).The standard for carbonate is Pee Dee Belemnite (PDB) (Craig, 1957) and that for water is Standard Mean Ocean Water (SMOW) (Baertschi, 1976).
Model results for the annual mean precipitation δ 18 O anomaly show large negative values in the northern and southern high latitudes at the LGM compared to the present (Fig. 1).Particularly visible are the area of very depleted δ 18 O precipitation over the imposed LGM ice-sheet in North America and northern Eurasia.In contrast, the midlatitudes are only slightly depleted and the tropical regions slightly enriched.Over the oceans and with the notable exception of the North Atlantic, the anomaly is relatively symmetric with respect to the equator.Some areas over Siberia exhibit a surprisingly large positive anomaly, despite their increased continentality and the negative δ 18 O anomaly in Western Europe.We do not consider this result as valid but rather, as a result of a model deficiency at very low moisture content, a deficiency that has already been noted for present-day (Roche, 2013).Introduction

Conclusions References
Tables Figures

Back Close
Full To compare our model results with isotopic data, we use the precipitation δ 18 O signal from ice cores.We focus on Greenland and Antarctic ice cores (Table 1).
For the Greenland ice-sheet, i LOVECLIM model result and data anomalies are in very good agreement (R 2 = 0.7) (Fig. 2).A gradient from weak negative δ 18 O anomalies from the North-West (NEEM, Camp Century) through the centre (GRIP, NGRIP) towards large negative δ 18 O anomalies in the South-East (Dye-3) is visible in both data and model results.  1 and Fig. 3).
Part of the mismatch between the model and the data at those sites may arise from the numerical humidity issue already outlined in Roche (2013).Some tropical ice cores exist in the Himalayan (Dunde, Guliya) and Andes regions (Huascaran, Sajama and Illimani).The tropics are marked by slightly positive δ 18 O anomalies or no change in i LOVECLIM (Fig. 1).The model fails to simulate the depletions ranging from −2 to −5 ‰ inferred from tropical ice cores (Thompson et al., 1989(Thompson et al., , 1995(Thompson et al., , 1998(Thompson et al., , 2000;;Ramirez et al., 2003;Risi et al., 2010).Part of the explanation could reside in the extreme altitude (around 5000-6000 m) of tropical ice cores, hampering a good data-model comparison.Indeed, the altitude in our coarse resolution model for these regions is lower (300 m for the Andes and 4000 m for the Himalayas).Another reason might be the complex precipitation setting of the Andes, Introduction

Conclusions References
Tables Figures

Back Close
Full with a combination of moisture from the Pacific and recycled moisture from inland regions (Risi et al., 2010).

Oxygen stable isotopes in continental carbonates (speleothems)
Under equilibrium conditions, the δ 18 O of continental carbonates (speleothems) depends on both the temperature through its control on equilibrium fractionation between water and calcite (Hendy, 1971;Kim and O'Neil, 1997) and the isotopic composition of drip water, from the cave site in which the speleothem grew, itself linked to the δ 18 O in precipitation.This relationship between the δ 18 O of carbonates, the δ 18 O of the water, and temperature is expressed by the equation of Kim and O'Neil (1997) for synthetic calcite: where T is the temperature in Kelvin.
Although recent works suggest that calcite speleothems do not precipitate under equilibrium conditions (Mickler et al., 2006;Daeron et al., 2011), the use of anomaly calculations between the LGM and LH would limit such potential bias.
As the atmospheric temperature is an important control on the calcite δ 18 O signal, we need to assess the modelled atmospheric temperature.To do so, we used the pollen-based continental climate reconstructions (Bartlein et al., 2011) (Fig. 4).A general cooling during the LGM is observed in both data and model results.Close to the Northern Hemisphere ice sheets the cooling is largest and well reproduced by the model.The LGM cooling indicated by pollen data is more pronounced in southern Europe in comparison to model results.This aspect was already noted (Vandenberghe et al., 2012) and was shown to reflect a too little southward extension of winter sea-ice along the western European coast.It could be attributed to (1) the low resolution of the atmospheric component and (2) the absence of Gibraltar in the oceanic part of the model (Roche and Caley, 2013), promoting warm waters from the Mediterranean along the European coast (Roche et al., 2007;Vandenberghe et al., 2012).In the tropics, the Introduction

Conclusions References
Tables Figures

Back Close
Full LGM and present is then compared to speleothem data (Fig. 5).Considering the error bars on the data (Table 1), we observe overall positive calcite δ 18 O anomalies except for Solufar cave (Fig. 5).As shown previously, mean annual atmospheric temperature are globally cooler during the LGM (Fig. 4).Precipitation reconstruction derived from subfossil pollens and plant macrofossils for the LGM suggests a significant decrease in precipitation compared to present over the Eurasia, Africa and North America (Bartlein et al., 2011).Both the significant cooling and drying of the LGM climate can explain the overall pattern toward positive calcite δ 18 O anomalies.Overall, for nine sites over the ten compiled, there is a very good agreement between data and model results (Fig. 5).In a previous data model comparison study for the late Holocene (Caley and Roche, 2013), we concluded that limitation of the model together with the processes operating in the atmosphere, soil zone, epikarst and cave system hampered a good quantitative data model comparison for the continental calcite δ 18 O signal.The better agreement between data and model results in term of annual mean anomaly suggests that this approach allows us to reduce complications with the atmospheric, soil and cave processes and that the model is capable to reproduce the right amplitude of changes.This offers large perspectives for the understanding of speleothems records covering glacial-interglacial time scale.Indeed, long-term transient simulation of water isotopes with i LOVECLIM could be realized and the relationship between the δ 18 O precipitation signal and climate variables such as temperature and precipitation rates could be investigated.Introduction

Conclusions References
Tables Figures

Back Close
Full

Oxygen stable isotopes in surface ocean carbonates (planktic foraminifera)
The carbonate isotopic concentration from various organisms such as foraminifera is mainly controlled by temperature and by the isotopic composition of seawater (δ 18 Osw) during shell formation (Urey, 1947;Shackleton, 1974).
The temperature dependence of the equilibrium fractionation of inorganic calcite precipitation around 16.9 • C is given in Shackleton (1974) as: As the calcite δ 18 O signal is controlled by temperature and δ 18 Osw, it is important to discuss and assess these variables in our model.The assessment can be carried out for sea surface temperature (SST) using reconstruction of LGM SST derived from different microfossil proxies (MARGO Project Members, 2009).On the contrary, there is currently no method to directly reconstruct surface water δ 18 O in the past.
We observe a very good agreement between simulated and measured SST anomalies between the LGM and the present (Fig. 6). Figure 6b illustrates the datamodel agreement or disagreement taking into accounts the uncertainties on LGM SST reconstructions (MARGO Project Members, 2009).Model results confirm that the strongest annual mean cooling occurred in the mid-latitude North Atlantic (MARGO Project Members, 2009).However, some discrepancies between SST reconstructions and model results occur south of Iceland and Greenland.We will discuss in detail this point later in the manuscript.In the tropical band (30 • S-30 • N) our model results are in excellent agreement with data and therefore confirm that the tropical cooling is more extensive than that proposed by CLIMAP (MARGO Project Members, 2009).Interbasin Introduction

Conclusions References
Tables Figures

Back Close
Full differences as well as west-east gradients within each basin, although much weaker in the model, mark the equatorial oceans in agreement with MARGO reconstructions.After subtraction of the LGM ice-sheets contribution (∼ 1 ‰) (Schrag et al., 1996;Duplessy et al., 2002), modelled surface water δ18 O anomalies exhibit negative values in the North Atlantic region (between 30 and 60 • N) (Fig. 7).This probably reflects the changes in ice-sheets distribution and their impact on surface water δ 18 O through depleted water discharges from rivers.Both sides of the Greenland ice-sheet are marked by positive anomalies (Fig. 7) that reflect the change from present day seasonal sea-ice to LGM permanent sea ice condition in agreement with proxy reconstruction (De vernal et al., 2006).Changes in ∆δ 18 Osw are null in the Southern Ocean and the rest of the oceans (tropical area) is mainly marked by slight positive anomalies, probably reflecting the more enriched δ 18 O precipitation signal during the LGM (Fig. 1).We computed annual mean calcite δ 18 O anomalies from simulated δ 18 Osw and SST and compared the results with deep-sea core data (Fig. 8).We chose a depth habitat of 0-50 m to calculate calcite δ 18 O anomalies as we previously demonstrated that it was suitable for a comparison with a global and varied dataset composed of different species of foraminifera (Caley and Roche, 2013).Although ecological effects can also play a role and are more expressed when individual species are considered (Caley and Roche, 2013) our strategy based on anomaly calculation limits such potential biases.We observe a good qualitative agreement between data and model results (Fig. 8a).
We compare these trends with a global latitudinal transect of ∆SST and the same latitudinal transect for the simulated ∆δ 18 Osw (Fig. 9b and c).We conclude that observed calcite δ 18 O gradients in the north Atlantic are mainly an effect of SST changes with latitude.Indeed, large positive calcite δ 18 O anomalies are associated with negative δ 18 Osw anomalies and colder temperatures indicating the dominant role of SST in driving the calcite δ 18 O signal (Fig. 9).
The data-model comparison also reveals large positive calcite δ 18 O anomalies in the Southern Ocean, between 45 and 50 • S (1.5-2 ‰).In a previous study, we argue that a data-model comparison for calcite δ 18 O in past climate could constitute an interesting way for mapping the potential shifts of the frontal systems and circulation changes through time, in particular in the Southern region (Caley and Roche, 2013).The large values observed are linked to a large negative SST anomaly in the region with weak changes in ∆δ 18 Osw observed in our model (Fig. 9).The cooling could be directly linked to the reorganisation of frontal systems during glacial periods as documented in many studies (Peeters et al., 2004;Bard and Rickaby, 2009;Caley et al., 2012).
However the driver of this potential fronts reorganisation is far from being understood.
It could be related to SH westerlies changes, although recent data-model comparison works have no clear conclusion on the behaviour of westerlies during the LGM (Kohfeld et al., 2013;Sime et al., 2013).It is interesting to note that the Southern Ocean, subtropical South Atlantic and Pacific are regions where large disagreements occur among the latest coupled-GCM LGM simulations (Braconnot et al., 2007a).The fact that i LOVECLIM reproduces well the observed cooling and the main pattern of the calcite δ 18 O signal in the Southern region illustrates how data-model comparison for oxygen isotopes can serve to evaluate model's capability to simulate a climate that is drastically different from that of the present-day.

Conclusions References
Tables Figures

Back Close
Full exception occurs in the North Indian Ocean for which we observed a higher calcite δ 18 O anomaly (2-3 ‰).This anomaly can not be explained by the cooling observed in the region but rather by a higher ∆δ 18 Osw (Figs. 9 and 10).Indeed, the cooling in the tropical Atlantic is more pronounced than in the north Indian Ocean but the calcite δ 18 O anomaly is not larger (Fig. 10).This large positive anomaly observed in data for the North Indian Ocean is overestimated in the model (Figs. 8,9 and 10).This is probably due to the anomalously low δ 18 Osw signal simulated by the model for the present day in the North Indian region (Roche and Caley, 2013).Explanations for the observed North Indian Ocean enrichment in δ 18 Osw at the LGM could be (1) a contraction of the Indian subtropical gyre and reduction of Agulhas leakage salty water (Caley et al., 2011a) and/or (2) an overall reduction of the hydrological cycle over the western and northern Asian region, in agreement with numerous Indo-Asian monsoonal reconstructions (Schultz et al., 1998;Iwamoto and Inouchi, 2007;Cheng et al., 2009;Guo et al., 2009;Caley et al., 2011b, c;Chabangborn et al., 2013).Also interesting is the low calcite δ 18 O anomaly observed in the China Sea (Figs. 8 and 10).This signal cannot be explained by a temperature effect as we observe a cooling more important in the China Sea in comparison to the North Indian Ocean (Figs. 6 and 10).Therefore, we hypothesize an important decrease of the ∆δ 18 Osw, a pattern exhibited in our model (Fig. 10c).The cause for such important decrease of the ∆δ 18 Osw is not completely clear because the monsoon in East Asia is rather reduced during the LGM (Iwamoto and Inouchi, 2007;Cheng et al., 2009;Guo et al., 2009).Nonetheless, some studies argue for substantial precipitation during the LGM in South China sea (Sun et al., 2000;Colin et al., 2010;Chabangborn et al., 2013).Indeed, part of the explanation could reside in the negative δ 18 O anomaly observed in precipitation over the China Sea (Fig. 1).
The use of calcite δ 18 O anomalies in the tropical regions (30 • N-30 • S) do not allow the confirmation of the presence of west-east SST gradients within each basin as the amount of data is rather limited.

Conclusions References
Tables Figures

Back Close
Full As discussed previously, the north Atlantic region exhibits steep calcite δ 18 O gradients that are also reflected in SST conditions (Figs. 8 and 9).However, large discrepancies remain with respect to reconstructed glacial temperatures based on different microfossil proxies (dinocysts and planktic foraminifera transfer function, alkenones and Mg/Ca ratio measurement on planktic foraminifera) in this region and there was no objective way to reconcile the divergent proxy results (MARGO Project Members, 2009).Our data-model comparison study of the calcite δ 18 O signal can shed light on these discrepancies.We therefore focussed on the North Atlantic and realized data-model comparison for calcite δ 18 O and SST signals (Fig. 11).The main difference between MARGO SST reconstructions and i LOVECLIM results corresponds to proxy reconstructions showing a positive SST anomaly in the region located south of the Iceland and Greenland ice-sheet.This positive anomaly disagrees with our model results but also with other proxy reconstructions (foraminifera transfer function) indicating an important cooling in the region during the LGM (Fig. 11a).
The calcite δ 18 O signal is strongly influenced by SST changes in the region and there is very good agreement between model and planktic foraminifera calcite data (Fig. 11).Gradients in calcite δ 18 O thus mirror SST changes.
In the south Iceland and Greenland region, data and model indicate large positive calcite δ 18 O anomalies (> 2 ‰) (Fig. 11b).These anomalies can solely be explained by a large negative SST anomaly.Therefore, the data indicating a positive anomaly or a weak cooling in the region are probably biased.After investigation, these data correspond to alkenones and, to a lesser extent, dinocyst reconstructions (Fig. 11a).Positive anomalies derived from dinocyst transfer functions have been interpreted as representing a "no-analogue" situation under partial sea ice coverage and glacial wind fields or a fine stratified layer that becomes warmer in summer (de Vernal et al., 2006).Concerning alkenones, several hypotheses can be proposed to explain the observed discrepancy.First, reconstructions are biased toward a specific season and can not be directly compared with our mean annual results (Rosell-Melé and Prahl, 2013).Second, culture studies of G. oceanica and E. huxleyi (Conte et al., 1998)  together with measurements of sinking particulate matter from Bermuda (Conte et al., 2001) suggest that the real shape of the alkenones-SST is probably sigmoidal, i.e. the calibration of the alkenone index converges asymptotically toward 0 and 1, for low and high temperatures, respectively (Conte et al., 2006).This could explain a bias with alkenone reconstructions at very low temperature.Finally, several early studies have noticed that alkenone-SST records are affected by laterally advected allochtonous input (Benthien and Müller, 2000;Mollenhauer et al., 2006;Rühlemann and Butzin, 2006).Anomalously warm SST during the last glacial period were observed in a marine sediment core recovered from the South East Indian Ridge (SEIR) at the location of the modern Subantarctic Front and attributed to a strong advection of detrital alkenones produced in warmer surface waters from the Agulhas region to SEIR (Sicre et al., 2005).These different hypotheses require further investigation in the region.Based on datamodel comparison for oxygen stable isotopes, we argue for a strong mean annual temperature anomaly between the LGM and present in the south of the Iceland and Greenland region (> 6 • C), supporting the foraminifera transfer function reconstruction.

Oxygen stable isotopes in deep ocean carbonates (benthic foraminifera)
Modelled results for δ 18 Osw changes in the ocean are presented on Fig. 12 for the deep Atlantic and Pacific Oceans after correction of the LGM ice-sheet contribution (1 ‰) (Schrag et al., 1996;Duplessy et al., 2002).Spatial differences in term of δ 18 Osw anomaly can be observed, suggesting that changes are not homogenous in the deep ocean.This is in agreement with reconstructions derived from pore fluids in deep-sea sediments (Adkins et al., 2002;Schrag et al., 2002;Malone et al., 2004) (Fig. 12a).
Considering the uncertainties on pore fluids reconstructions (0.1 ‰) (Adkins et al., 2002;Schrag et al., 2002), the data point that exhibits a negative δ 18 Osw anomaly in the deep Atlantic (∼ 4500 m) is the only measurement in significant disagreement with the model results (Fig. 12a).The value of this data point is surprising as it would imply a vigorous and/or deeper extension of the north Atlantic deep water at the LGM, in 121 Introduction

Conclusions References
Tables Figures

Back Close
Full anomaly is influenced by temperature.The modelled deep-water temperature for the present is around 2 • C lower than data, a pattern particularly marked in the Southern Ocean (Caley and Roche, 2013).This introduces bias in our data-model comparison.
The uses of anomalies only slightly limit such bias because deep LGM temperatures are close to the freezing point of seawater at the ocean's surface (Adkins et al., 2002).Reconstructions of past deep temperature are rather rare and localised.The north Atlantic, South Indian and equatorial Pacific are marked by negative temperature anomalies of ∼ 4, 3 and 2 • C respectively between the LGM and present (Adkins et al., 2002;Waelbroeck et al., 2002;Siddall et al., 2010) (Fig. 13e and f).Similarly, reconstruction of deep temperature anomaly in the South Pacific yields ∼ 2 • C (Malone et al., 2004;Elderfield et al., 2012) (Fig. 13f).The modelled deep temperature anomalies are in good agreement for the North Atlantic at 2000 but not at 4500 m as for the δ 18 Osw anomaly.For the Southern Ocean and deep equatorial Pacific, the negative anomalies between the LGM and present are about 1 • C too weak in the model (Fig. 13f).
Although there are some discrepancies between modelled and measured deepwater temperature, our data model comparison support a heterogeneous cooling of few degrees (2-4 • C) in the LGM Ocean.Introduction

Conclusions References
Tables Figures

Back Close
Full is in very good agreement with ice core records, whereas simulated values are in less good agreement with data over Antarctica.We observe a general pattern toward more enriched calcite δ 18 O in the model over the continents at the LGM, in agreement with speleothem data.This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation, as confirmed by reconstructions derived from pollens and plant macrofossils (Bartlein et al., 2011).The good agreement between data and model results in term of annual mean calcite δ 18 O anomaly offers large perspectives for the understanding of speleothems records covering glacial-interglacial time scale.Long-term transient simulations of water isotopes with i LOVECLIM are planned in the near future and the relationship between precipitation δ 18 O and climate variables such as temperature and precipitation will be investigated.
Data-model comparison for sea surface temperature indicates that i LOVECLIM is capable to satisfyingly simulate oceanic surface conditions at the LGM, whereas the majority of AO-GCM simulations experience some difficulties (Braconnot et al., 2007a).Large discrepancies with respect to glacial temperatures recorded by different microfossil proxies remain in the North Atlantic region and there was no objective way to reconcile the divergent proxy results (MARGO Project Members, 2009).
Water isotopes have been implemented and validated against data for the present day climate in the global three-dimensional model of intermediate complexity 107 Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

δ
18 O measurements compiled byZarriess and Mackensen (2011).These authors used a definition for the LGM chronozone (18.3-23.5 ka) that encompasses the one we used (19-23 ka).We also considered four calcite δ 18 O measurements from Adkins et al. (2002) and two calcite δ 18 O measurements from Malone et al. (2004) that have been combined in the same cores with pore fluids measurements in deep-sea sediments to reconstruct deep ocean temperature.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | cooling is reduced compared to the high northern latitudes.With the exception of South Africa where a lot of scatter is observed in the data from a large cooling of −10 • C to a moderate −2 • C cooling, i LOVECLIM reproduces very well the main features of atmospheric temperature between the LGM and present.Simulated atmospheric temperature and precipitation δ 18 O are used to calculate calcite δ 18 O.The simulated annual mean anomaly between the Discussion Paper | Discussion Paper | Discussion Paper | Figure8billustrates the data-model agreement or disagreement taking into accounts the uncertainties (2σ) on LGM and LH calcite δ 18 O reconstructions (TableS1).Overall, we observe quantitative good agreement between data and model except in the North Indian region.Although calcite δ 18 O anomalies are larger in the model than in the data, Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | disagreement with reconstructions based on δ 13 C and Cd/Ca proxies (Lynch-Stieglitz et al., 2007 and references therein) and with our model results.Qualitative and quantitative comparison of calcite δ 18 O anomaly, calculated with the model and benthic foraminifera data with associated error bars (2σ) (TableS2), are visible on Fig.13a-dfor the deep Atlantic and Pacific Oceans.Some discrepancies can be observed between data and model and are particularly marked in the equatorial Atlantic at ∼ 1000 m, the deep Southern Ocean and central Pacific.The calcite δ 18 O fully coupled atmosphere-ocean three-dimensional model of intermediate complexity i LOVECLIM to simulate the climate and oxygen stable isotopes (δ 18 O) in the atmospheric and oceanic component during the LGM.We also realized a careful compilation of global oxygen isotopic datasets to assess the model performance and constrain the LGM climate.Model results for the annual mean precipitation δ 18 O show more depleted values in the northern and southern high latitudes during the LGM than at present.The simulated spatial gradient in precipitation δ 18 O over Greenland Our datamodel comparison for planktic foraminiferal calcite δ 18 O indicates that a strong mean annual cooling characterized the LGM with respect to the present (> 6 • C), supporting Discussion Paper | Discussion Paper | Discussion Paper | the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions.The data-model comparison also reveals large positive calcite δ 18 O anomalies in the Southern Ocean linked to an important cooling that could be linked to a reorganization of frontal systems during the LGM.Nonetheless, the exact driver of this pattern remains unclear.From our data-model comparison of planktic foraminifer oxygen stable isotopes and SST we deduced a large positive/negative δ 18 Osw anomaly for the north Indian Ocean/China Sea between the LGM and present which may be explained by changes in the hydrological cycle over the region.Our simulation of the deep ocean suggests that changes in δ 18 Osw between the LGM and present are not spatially homogenous.This is supported by reconstructions derived from pore fluids in deep-sea sediments.The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite δ 18 O data.Some experiments are planned in a near future to try modifying this aspect in the model.Nonetheless, our data model comparison support a heterogeneous cooling of few degrees (2-4 • C) in the LGM Ocean.i LOVECLIM reproduces well the δ 18 O signals between the LGM and present and therefore illustrates how data-model comparison for oxygen isotopes can serve to evaluate model's capability to simulate a climate that is drastically different from that of the present-day.Supplementary material related to this article is available online at http://www.clim-past-discuss.net/10/105/2014/cpd-10DiscussionPaper | Discussion Paper | Discussion Paper |

Table 1
also includes calcite δ 18 O data from 10 published speleothem records.Calcite δ 18 O anomalies are reported as the difference between averaged δ 18 O values computed over the period of 20 000-22 000 yr BP and over the last 1000 yr of each record, using published U/Th chronologies.
The maximum depletion occurs in our model over Baffin Bay, with a strong effect of glacial/interglacial sea-ice changes.Concerning Antarctica, i LOVECLIM results are not in good agreement (R 2 = 0.19) (Fig.3, right panel) with the small data set used.Looking in more details at individual ice-core sites, we infer that the model does reproduce qualitatively the east -west gradient in Antarctica, with large negative precipitation δ 18O anomalies in western Antarctica (Epica Dronning Maud Land, Byrd, WDC, Siple Dome) and smaller negative anomalies in eastern Antarctica.The δ 18 O anomaly at Talos Dome ice core record is well reproduced in the model.Similarly, the values simulated at the Vostok and Dome B sites are not regionally inconsistent.The only ice-core sites that are really at odd with our results are Dome Fuji and Epica Dome C/Dome C data (Table