a b s t r a c tThe Last Interglacial (LIG, 129e116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain.Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG.We have also compiled four data-based time slices with temperature anomalies (compared to presentday conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions. Crown

Abstract. Understanding natural and anthropogenic climate change processes involves using computational models that represent the main components of the Earth system: the atmosphere, ocean, sea ice, and land surface. These models have become increasingly computationally expensive as resolution is increased and more complex process representations are included. However, to gain robust insight into how climate may respond to a given forcing, and to meaningfully quantify the associated uncertainty, it is often required to use either or both ensemble approaches and very long integrations. For this reason, more computationally efficient models can be very valuable tools. Here we provide a comprehensive overview of the suite of climate models based around the HadCM3 coupled general circulation model. This model was developed at the UK Met Office and has been heavily used during the last 15 years for a range of future (and past) climate change studies, but has now been largely superseded for many scientific studies by more recently developed models. However, it continues to be extensively used by various institutions, including the BRIDGE (Bristol Research Initiative for the Dynamic Global Environment) research group at the University of Bristol, who have made modest adaptations to the base HadCM3 model over time. These adaptations mean that the original documentation is not entirely representative, and several other relatively undocumented configurations are in use. We therefore describe the key features of a number of configurations of the HadCM3 climate model family, which together make up HadCM3@Bristol version 1.0. In order to differentiate variants that have undergone development at BRIDGE, we have introduced the letter B into the model nomenclature. We include descriptions of the atmosphere-only model (HadAM3B), the coupled model with a low-resolution ocean (HadCM3BL), the high-resolution atmosphere-only model (HadAM3BH), and the regional model (HadRM3B). These also include three versions of the land surface scheme. By comparing withPublished by Copernicus Publications on behalf of the European Geosciences Union. observational datasets, we show that these models produce a good representation of many aspects of the climate system, including the land and sea surface temperatures, precipitation, ocean circulation, and vegetation. This evaluation, combined with the relatively fast computational speed (up to 1000 times faster than some CMIP6 models), motivates continued development and scientific use of the HadCM3B family of coupled climate models, predominantly for quantifying uncertainty and for long multi-millennial-scale simulations.

It is thought that the Northern Hemisphere experienced only ephemeral glaciations from the Late Eocene to the Early Pliocene epochs (about 38 to 4 million years ago), and that the onset of extensive glaciations did not occur until about 3 million years ago. Several hypotheses have been proposed to explain this increase in Northern Hemisphere glaciation during the Late Pliocene. Here we use a fully coupled atmosphere-ocean general circulation model and an ice-sheet model to assess the impact of the proposed driving mechanisms for glaciation and the influence of orbital variations on the development of the Greenland ice sheet in particular. We find that Greenland glaciation is mainly controlled by a decrease in atmospheric carbon dioxide during the Late Pliocene. By contrast, our model results suggest that climatic shifts associated with the tectonically driven closure of the Panama seaway, with the termination of a permanent El Niño state or with tectonic uplift are not large enough to contribute significantly to the growth of the Greenland ice sheet; moreover, we find that none of these processes acted as a priming mechanism for glacial inception triggered by variations in the Earth's orbit.

A Community Climate System Model, Version 3 (CCSM3) simulation for 125 ka during the Last Interglacial (LIG) is compared to two recent proxy reconstructions to evaluate surface temperature changes from modern times. The dominant forcing change from modern, the orbital forcing, modified the incoming solar insolation at the top of the atmosphere, resulting in large positive anomalies in boreal summer. Greenhouse gas concentrations are similar to those of the pre-industrial (PI) Holocene. CCSM3 simulates an enhanced seasonal cycle over the Northern Hemisphere continents with warming most developed during boreal summer. In addition, year-round warming over the North Atlantic is associated with a seasonal memory of sea ice retreat in CCSM3, which extends the effects of positive summer insolation anomalies on the high-latitude oceans to winter months. The simulated Arctic terrestrial annual warming, though, is much less than the observational evidence, suggesting either missing feedbacks in the simulation and/or interpretation of the proxies. Over Antarctica, CCSM3 cannot reproduce the large LIG warming recorded by the Antarctic ice cores, even with simulations designed to consider observed evidence of early LIG warmth in Southern Ocean and Antarctica records and the possible disintegration of the West Antarctic Ice Sheet. Comparisons with a HadCM3 simulation indicate that sea ice is important for understanding model polar responses. Overall, the models simulate little global annual surface temperature change, while the proxy reconstructions suggest a global annual warming at LIG (as compared to the PI Holocene) of approximately 1 ° C, though with possible spatial sampling biases. The CCSM3 SRES B1 (low scenario) future projections suggest high-latitude warmth similar to that reconstructed for the LIG may be exceeded before the end of this century.

Abstract. Ice thickness and bedrock topography are essential boundary conditions for numerical modelling of the evolution of the Greenland ice-sheet (GrIS). The datasets currently in use by the majority of GrIS modelling studies are over two decades old and based on data collected from the 1970s and 80s. We use a newer, high-resolution Digital Elevation Model of the GrIS and new temperature and precipitation forcings to drive the Glimmer ice-sheet model offline under steady state, present day climatic conditions. Comparisons are made of ice-sheet geometry between these new datasets and older ones used in the EISMINT-3 exercise. We find that changing to the newer bedrock and ice thickness makes the greatest difference to Greenland ice volume and ice surface extent. When all boundary conditions and forcings are simultaneously changed to the newer datasets the ice-sheet is 33% larger in volume compared with observation and 17% larger than that modelled by EISMINT-3.We performed a tuning exercise to improve the modelled present day ice-sheet. Several solutions were chosen in order to represent improvement in different aspects of the GrIS geometry: ice thickness, ice volume and ice surface extent. We applied these new parameter sets for Glimmer to several future climate scenarios where atmospheric CO 2 concentration was elevated to 400, 560 and 1120 ppmv (compared with 280 ppmv in the control) using a fully coupled General Circulation Model. Collapse of the ice-sheet was found to occur between 400 and 560 ppmv, a threshold substantially lower than previously modelled using the standard EISMINT-3 setup. This work highlights the need to assess carefully boundary conditions and forcings required by Correspondence to: E. J. Stone (emma.j.stone@bristol.ac.uk) ice-sheet models, particularly in terms of the abstractions required for large-scale ice-sheet models, and the implications that these can have on predictions of ice-sheet geometry under past and future climate scenarios.