(Table S1) Boron isotope and trace element data and reconstructed carbonate system parameters, ODP Holes 122-761B and 154-926A

The middle Miocene Climatic Optimum (17-15 Ma; MCO) is a period of global warmth and relatively high CO2 and is thought to be associated with a significant retreat of the Antarctic Ice Sheet (AIS). We present here a new planktic foraminiferal d11B record from 16.6 to 11.8 Ma from two deep ocean sites currently in equilibrium with the atmosphere with respect to CO2. These new data demonstrate that the evolution of global climate during the middle Miocene (as reflected by changes in the cyrosphere) was well correlated to variations in the concentration of atmospheric CO2. What is more, within our sampling resolution (~1 sample per 300 kyr) there is no evidence of hysteresis in the response of ice volume to CO2 forcing during the middle Miocene, contrary to what is understood about the Antarctic Ice Sheet from ice sheet modelling studies. In agreement with previous data, we show that absolute levels of CO2 during the MCO were relatively modest (350-400 ppm) and levels either side of the MCO are similar or lower than the pre-industrial (200-260 ppm). These new data imply the presence of either a very dynamic AIS at relatively low CO2 during the middle Miocene or the advance and retreat of significant northern hemisphere ice. Recent drilling on the Antarctic margin and shore based studies indicate significant retreat and advance beyond the modern limits of the AIS did occur during the middle Miocene, but the complete loss of the AIS was unlikely. Consequently, it seems that ice volume and climate variations during the middle Miocene probably involved a more dynamic AIS than the modern but also some component of land-based ice in the northern hemisphere.

(Table 1) Age determination of ODP Sites 124-769 and 124-768

The Sulu Sea is located in the 'warm pool' of the western Pacific Ocean, where mean annual temperatures are the highest of anywhere on Earth. Because this large heat source supplies the atmosphere with a significant portion of its water vapour and latent heat, understanding the climate history of the region is important for reconstructing global palaeoclimate and for predicting future climate change. Changes in the oxygen isotope composition of planktonic foraminifera from Sulu Sea sediments have previously been shown to reflect changes in the planetary ice volume at glacial-interglacial and millenial timeseales, and such records have been obtained for the late Pleistocene epoch and the last deglaciation (Linsley and Thunell, 1990, doi:10.1029/PA005i006p01025; Lindley and Dunbar, 1994, doi:10.1029/93PA03216; Kudrass et al., 1991, doi:10.1038/349406a0). Here I present results that extend the millenial time resolution record back to 150,000 years before present. On timescales of around 10,000 years, the Sulu Sea oxygen-isotope record matches changes in sea level deduced from coral terraces on the Huon peninsula (Chappell and Shackleton, doi:10.1038/324137a0). This is particularly the case during isotope stage 3 (an interglacial period 23,000 to 58,000 years ago) where the Sulu Sea oxygen-isotope record deviates from the SPECMAP deep-ocean oxygen-isotope record (Imbrie et al., 1984). Thus these results support the idea (Chappell and Shackleton, doi:10.1038/324137a0; Shackleton, 1987, doi:10.1016/0277-3791(87)90003-5) that there were higher sea levels and less continental ice during stage 3 than the SPECMAP record implies and that sea level during this interglacial was just 40-50 metres below present levels. The subsequent rate of increase in continental ice volume during the return to full glacial conditions was correspondingly faster than previously thought.

Sea-surface temperature reconstruction in the Southwest Pacific

Uniquely in the Southern Hemisphere the New Zealand micro-continent spans the interface between a subtropical gyre and the Subantarctic Circumpolar Current. Its 20° latitudinal extent includes a complex of submerged plateaux, ridges, saddles and basins which, in the present interglacial, are partial barriers to circulation and steer the Subtropical (STF) and Subantarctic (SAF) fronts. This configuration offers a singular opportunity to assess the influence of bottom topography on oceanic circulation through Pleistocene glacial - interglacial (G/I) cycles, its effect on the location and strength of the fronts, and its ability to generate significant differences in mixed layer thermal history over short distances. For this study we use new planktic foraminiferal based sea-surface temperature (SST) estimates spanning the past 1 million years from a latitudinal transect of four deep ocean drilling sites. We conclude that: 1. the effect of the New Zealand landmass was to deflect the water masses south around the bathymetric impediments; 2. the effect of a shallow submerged ridge on the down-current side (Chatham Rise), was to dynamically trap the STF along its crest, in stark contrast to the usual glacial-interglacial (G-I) meridional migration that occurs in the open ocean; 3. the effect of more deeply submerged, downstream plateaux (Campbell, Bounty) was to dynamically trap the SAF along its steep southeastern margin; 4. the effects of saddles across the submarine plateaux was to facilitate the development of jets of subtropical and subantarctic surface water through the fronts, forming localized downstream gyres or eddies during different phases in the G-I climate cycles; 5. the deep Pukaki Saddle across the Campbell-Bounty Plateaux guided a branch of the SAF to flow northwards during each glacial, to form a strong gyre of circumpolar surface water in the Bounty Trough, especially during the mid-Pleistocene Climate Transition (MIS 22-16) when exceptionally high SST gradients existed across the STF; 6. the shallower Mernoo Saddle, at the western end of the Chatham Rise, provided a conduit for subtropical water to jet southwards across the STF in the warmest interglacial peaks (MIS 11, 5.5) and for subantarctic water to flow northwards during glacials; 7. although subtropical or subantarctic drivers can prevail at a particular phase of a G-I cycles, it appears that the Antarctic Circumpolar Current is the main influence on the regional hydrography. Thus complex submarine topography can affect distinct differences in the climate records over short distances with implications for using such records in interpreting global or regional trends. Conversely, the local topography can amplify the paleoclimate record in different ways in different places, thus enhancing its value for the study of more minor paleoceanographic influences that elsewhere are more difficult to detect. Such sites include DSDP 594, which like some other Southern Ocean sites, has the typical late Pleistocene asymmetrical saw-tooth G-I climate pattern transformed to a gap-tooth pattern of quasi-symmetrical interglacial spikes that interrupt extended periods of minimum glacial temperatures.

Pliocene-Pleistocene stack of globally distributed benthic stable oxygen isotope records

We present a 5.3-Myr stack (the ''LR04'' stack) of benthic d18O records from 57 globally distributed sites aligned by an automated graphic correlation algorithm. This is the first benthic delta18O stack composed of more than three records to extend beyond 850 ka, and we use its improved signal quality to identify 24 new marine isotope stages in the early Pliocene. We also present a new LR04 age model for the Pliocene-Pleistocene derived from tuning the delta18O stack to a simple ice model based on 21 June insolation at 65 N. Stacked sedimentation rates provide additional age model constraints to prevent overtuning. Despite a conservative tuning strategy, the LR04 benthic stack exhibits significant coherency with insolation in the obliquity band throughout the entire 5.3 Myr and in the precession band for more than half of the record. The LR04 stack contains significantly more variance in benthic delta18O than previously published stacks of the late Pleistocene as the result of higher resolution records, a better alignment technique, and a greater percentage of records from the Atlantic. Finally, the relative phases of the stack's 41- and 23-kyr components suggest that the precession component of delta18O from 2.7-1.6 Ma is primarily a deep-water temperature signal and that the phase of d18O precession response changed suddenly at 1.6 Ma.