(Table S3) Benthic foraminiferal Mg/Ca, Li/Ca and stable isotope ratios of ODP Site 122-761

The sensitivities of benthic foraminiferal Mg/Ca and Li/Ca to bottom water temperature and carbonate saturation state have recently been assessed. Here we present a new approach that uses paired Mg/Ca and Li/Ca records to calculate simultaneous changes in temperature and saturation state. Using previously published records, we first use this approach to document a cooling of deep ocean waters associated with the establishment of the Antarctic ice sheet at the Eocene-Oligocene climate transition. We then apply this approach to new records of the Middle Miocene Climate Transition from ODP Site 761 to estimate variations in bottom water temperature and the oxygen isotopic composition of seawater. We estimate that the oxygen isotopic composition of seawater varied by ~1 per mil between the deglacial extreme of the Miocene Climatic Optimum and the glacial maximum following the Middle Miocene Climate Transition, indicating large amplitude variations in ice volume. However, the longer-term change between 15.3 and 12.5 Ma is marked by a ~1°C cooling of deep waters, and an increase in the oxygen isotopic composition of seawater of ~0.6 per mil. We find that bottom water saturation state increased in the lead up to the Middle Miocene Climate Transition and decreased shortly after. This supports decreasing pCO2 as a driver for global cooling and ice sheet expansion, in agreement with existing boron isotope and leaf stomatal index CO2 records but in contrast to the published alkenone CO2 records.

Planktic and benthic foraminifreral isotope data for ODP Site 207-1258

The mid-Cretaceous is widely considered the archetypal ice-free greenhouse interval in Earth history, with a thermal maximum around Cenomanian-Turonian boundary time (ca. 90 Ma). However, contemporaneous glaciations have been hypothesized based on sequence stratigraphic evidence for rapid sea-level oscillation and oxygen isotope excursions in records generated from carbonates of questionable preservation and/or of low resolution. We present new oxygen isotope records for the mid-Cenomanian Demerara Rise that are of much higher resolution than previously available, taken from both planktic and benthic foraminifers, and utilizing only extremely well preserved glassy foraminifers. Our records show no evidence of glaciation, calling into question the hypothesized ice sheets and rendering the origin of inferred rapid sea-level oscillations enigmatic. Simple mass-balance calculations demonstrate that this Cretaceous sea-level paradox is unlikely to be explained by hidden ice sheets existing below the limit of d18O detection.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Quaternary isotope stratigraphy of DSDP Hole 90-593 from the Challenegr Plateau, South Tasman Sea

A moderate-resolution isotope stratigraphy (with an average of one sample per 17,500 yr.) derived from the benthic foraminifer Uvigerina (or Cibicides), the planktonic foraminifer Globigerina bulloides, and calcareous nannofossil concentrates is presented for the entire Quaternary (and latest Pliocene) section of mid-upper bathyal calcareous oozes from DSDP Site 593, western Challenger Plateau, south Tasman Sea. Superimposed on a trend of gradually increasing average delta18O values through the Pleistocene, reflecting the progressive buildup of polar ice sheets, is a record of highfrequency but generally low amplitude (0.5-1?) isotope fluctuations in the early Quaternary (1.9-1.0 m.y.), followed by a greatly increased intensity (1.5-2.0 ?) of glacial-interglacial fluctuations during the late Quaternary (< 1.0 m.y.). The standard late Quaternary isotope stages 1 to 24 are mainly resolvable. Significant excursions in both delta18O and delta13C values at various times during the Quaternary are suggested to be due to periodic, fundamental changes in ocean circulation properties over the plateau. For example, intensified upwelling of Antarctic Intermediate Waters during several glacial periods is indicated by the convergence of benthic and planktonic foraminiferal delta18O data, and productivity variations may account for certain delta13C spikes in the record. With increasingly higher resolution analysis this core will provide a useful Quaternary isotope reference section for southern temperate waters in the southwest Pacific, centered on New Zealand.

Paleocene radiolarian distribution and biostratigraphy from offshore eastern New Zealand

A 100-m-thick Paleocene sequence of mainly pelagic sediments at ODP Site 1121, on the eastern flanks of the Campbell Plateau, contains few to common radiolarians of relatively low diversity in the lower 40 m (Early to early Late Paleocene) and abundant, diverse radiolarian assemblages in the upper 60 m (mid-Late Paleocene). The 150 taxa recorded from the entire Paleocene interval are thought to under-represent the actual species diversity by at least one half as many morphotypes have not been differentiated below the level of genus. Assemblages in the lower 40 m are similar to those described from onland New Zealand and DSDP Site 208 (northern Lord Howe Rise); they are correlated with South Pacific radiolarian zones RP4 and RP5. Assemblages in the upper 60 m differ from other known Late Paleocene assemblages in the great abundance of plagiacanthids and cycladophorids. Similarities are noted with later Cenozoic cool-water assemblages. This upper interval is correlated with South Pacific zone RP6, as revised herein, based on comparison with faunas from Site 208 and Marlborough, New Zealand. The interval is also correlated with the upper part of North Atlantic zone RP6 (RP6b-c) based on the presence of Aspis velutochlamydosaurus, Plectodiscus circularis and Pterocodon poculum. Other species, such as Buryella tetradica and Buryella pentadica, are valuable for local correlation but exhibit considerable diachroneity between the Pacific, Indian and Atlantic Oceans. An age model for the Paleocene interval at Site 1121, based on well-constrained nannofossil and radiolarian datums, indicates that the rate of compacted sediment accumulation doubles from 15 to 30 mm/ka at the RP5/RP6 zonal boundary. In large part this is due to a sudden and pronounced increase in accumulation rates for all siliceous fossils; radiolarians and larger diatoms increase from <100 to >10 000 specimens/cm2/ka. This apparent increase in biosiliceous productivity is age-equivalent to a mid-Paleocene cooling event (57-59 Ma) identified from global stable isotope records that is associated with the heaviest delta13C values for the entire Cenozoic.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Oxygen and carbon isotope stratigraphy of ODP Holes 107-653A and 107-654A

Stable isotope analysis of two species (or groups of species) of planktonic foraminifers: Globigerinoides ruber (or G. obliquus and G. obliquus extremus) and Globigerina bulloides (or G. falconensis and G. obesa) from ODP Hole 653A and Site 654 in the Tyrrhenian basin, records the Pliocene-Pleistocene glacial history of the Northern Hemisphere. The overall increase in mean d18O values through the interval 4.6-0.08 Ma is 1.7 per mil for G. bulloides and 1.5 per mil for G. ruber. The time interval 3.1-2.5 Ma corresponds to an important phase of 18O enrichment for planktonic foraminifers. In this interval, glacial d18O values of both species G. bulloides and G. ruber increase by about l per mil, this increase being more progressive for G. ruber than for G. bulloides. The increase of interglacial d18O values is higher for G. bulloides (1.5 per mil) than for the Gruber group (1 per mil). These data suggest a more pronounced seasonal stratification of the water masses during interglacial phases. Large positive d18O fluctuations of increasing magnitude are also recorded at 2.25 and 2.15 Ma by G bulloides and appear to be diachronous with those of Site 606 in the Atlantic Ocean. Other events of increasing d18O values are recorded between 1.55 and 1.3 Ma, at 0.9 Ma, 0.8 Ma, and near 0.34 Ma. In the early Pliocene the d18O variability recorded by the planktonic species G. bulloides was higher in the Mediterranean than in the Atlantic at the same latitude. This suggests that important cyclic variations in the water budget of the Mediterranean occurred since that time. Step increases in the d18O variability are synchronous with those of the open ocean at 0.9 and 0.34 Ma. The higher variability as well as the higher amplitude of the peaks of 18O enrichment may be partly accounted for by increase of dryness over the Mediterranean area. In particular the high amplitude d18O fluctuations recorded between 3.1 and 2.1 Ma are correlated with the onset of a marked seasonal contrast and a summer dryness, revealed by pollen analyses. Strong fluctuations towards d13C values higher than modern ones are recorded by the G. ruber group species before 1.7 Ma and suggest a high production of phytoplankton. When such episodes of high primary production are correlated with episodes of decreasing 13C content of G. bulloides, they are interpreted as the consequence of a higher stratification of the upper water masses resulting itself from a marked seasonality. Such episodes occur between 4.6 and 4.05 Ma, 3.9 and 3.6 Ma, and 3.25 and 2.66 Ma. The interval 2.66-1.65 Ma corresponds to a weakening of the stratification of the upper water layers. This may be related to episodes of cooling and increasing dryness induced by the Northern Hemisphere Glaciations. The Pleistocene may have been a less productive period. The transition from highly productive to less productive surface waters also coincides with a new step increase in dryness and cooling, between 1.5 and 1.3 Ma. The comparison of the 13C records of G ruber and G. bulloides in fact suggests that a high vertical convection became a dominant feature after 2.6 Ma. Increases in the nutrient input and the stratification of the upper water masses may be suspected, however, during short episodes near 0.86 Ma (isotopic stage 25), 0.57-0.59 Ma (isotopic stage 16), 0.49 Ma (isotopic stage 13), 0.4-0.43 Ma (isotopic stage 11), and 0.22 and 0.26 Ma (part of isotopic stage 7 and transition 7/8). In fact, changes in the C02 balance within the different water masses of the Tyrrhenian basin as well as in the local primary production did not follow the general patterns of the open ocean.

Late Quaternary benthic foraminiferal record of the Bounty Trough, east of New Zealand

The Bounty Trough, east of New Zealand, lies along the southeastern edge of the present-day Subtropical Front (STF), and is a major conduit via the Bounty Channel, for terrigenous sediment supply from the uplifted Southern Alps to the abyssal Bounty Fan. Census data on 65 benthic foraminiferal faunas (>63 µm) from upper bathyal (ODP 1119), lower bathyal (DSDP 594) and abyssal (ODP 1122) sequences, test and refine existing models for the paleoceanographic and sedimentary history of the trough through the last 150 ka (marine isotope stages, MIS 6-1). Cluster analysis allows recognition of six species groups, whose distribution patterns coincide with bathymetry, the climate cycles and displaced turbidite beds. Detrended canonical correspondence analysis and comparisons with modern faunal patterns suggest that the groups are most strongly influenced by food supply (organic carbon flux), and to a lesser extent by bottom water oxygen and factors relating to sediment type. Major faunal changes at upper bathyal depths (1119) probably resulted from cycles of counter-intuitive seaward-landward migrations of the Southland Front (SF) (north-south sector of the STF). Benthic foraminiferal changes suggest that lower nutrient, cool Subantarctic Surface Water (SAW) was overhead in warm intervals, and higher nutrient-bearing, warm neritic Subtropical Surface Water (STW) was overhead in cold intervals. At lower bathyal depths (594), foraminiferal changes indicate increased glacial productivity and lowered bottom oxygen, attributed to increased upwelling and inflow of cold, nutrient-rich, Antarctic Intermediate Water (AAIW) and shallowing of the oxygen-minimum zone (upper Circum Polar Deep Water, CPDW). The observed cyclical benthic foraminiferal changes are not a result of associations migrating up and down the slope, as glacial faunas (dominated by Globocassidulina canalisuturata and Eilohedra levicula at upper and lower bathyal depths, respectively) are markedly different from those currently living in the Bounty Trough. On the abyssal Bounty Fan (1122), faunal changes correlate most strongly with grain size, and are attributed to varying amounts of mixing of displaced and in-situ faunas. Most of the displaced foraminifera in turbiditic sand beds are sourced from mid-outer shelf depths at the head of the Bounty Channel. Turbidity currents were more prevalent during, but not restricted to, glacial intervals.

Pliocene stable oxygen and carbon isotope record of benthic foraminifera from ODP Site 138-846 in the eastern equatorial Pacific

Oxygen and carbon isotope ratios were measured in benthic foraminifers from the entire Pliocene and latest Miocene sections of Site 846, a 180-m section, at a sampling interval of 10 cm. This provides a temporal resolution of about 2500 yr. The documented continuity of the record is excellent. Using the time scale that was developed on the basis of orbital tuning of GRAPE density records, we observed a fairly constant phase relationship between delta18O and variations in the obliquity of Earth's rotational axis. A new numbering scheme for Pliocene isotope stages is proposed. This high-resolution delta18O record clarifies several interesting aspects of late Neogene climatic evolution, including a "glacial" event that may have caused the final Messinian desiccation of the Mediterranean Sea; one or more "interglacial" events that might have caused refilling of the Mediterranean; a well-resolved couplet of glacial events at about the age of the Sidujfall Subchron; interglacial extremes in the early part of the Gauss that could have resulted from either significant deglaciation on Antarctica or from warming of deep water; and a gradual ramp of increasingly extreme "glacial" events, starting at about the Kaena Subchron and culminating with delta18O stage 100 in the earliest Matuyama.

Temperatures, salinity, and stable oxygen isotope values of bottom water and of the mollusc shell Arctica islandica from Nantucket Shoals

The carbonate shell of the bivalve Arctica islandica has been recognized, for more than a decade, as a potentially important marine geochemical biorecorder owing to this species' great longevity (200+ years) and wide geographic distribution throughout the northern North Atlantic Ocean, a region vital to global climate and ocean circulation. However, until now this potential has not been realized owing to the difficulty of precisely sampling the shell of this slow growing species. Using newly available automated microsampling techniques combined with micromass stable isotope mass spectrometry, a stable oxygen isotope record (1956-1957 and 1961-1970) has been obtained from a live-captured, 38-year-old A. islandica specimen collected near the former position of the Nantucket Shoals Lightship (41°N. 69°W). The shell's delta18O signal is compared with an expected signal derived from ambient bottom temperature and salinity data recorded at the lightship for the same period. The results show that A islandica's delta18O record (1) is in phase with its growth banding, confirming the annual periodicity of this species' growth bands, (2) is in oxygen isotopic equilibrium with the ambient seawater, (3) shows a consistent shell growth shutdown temperature of ~6°C. which translates into an ~8-month (May-December) shell growth period at this location, and (4) records the ambient bottom temperature with a precision of ~ +/-1.2°C. These results add important information on the life history of this commercially important shellfish species and demonstrate that A. islandica shells can be used to reconstruct inter- and intra-annual records of the continental shelf bottom temperature.

Nd and Pb isotope data for DSDP and ODP sediments from the Bengal Fan

Correlation of new multichannel seismic profiles across the upper Indus Fan and Murray Ridge with a dated industrial well on the Pakistan shelf demonstrates that ~40% of the Indus Fan predates the middle Miocene, and ~35% predates uplift of the Murray Ridge (early Miocene, ~22 Ma). The Arabian Sea, in addition to the Makran accretionary complex, was therefore an important repository of sediment from the Indus River system during the Paleogene. Channel and levee complexes are most pronounced after the early Miocene, coincident with an increase in sedimentation rates. Middle Eocene sandstones from Deep Sea Drilling Project Site 224 on the Owen Ridge yield K-feldspars whose Pb isotopic composition, measured by in situ ion microprobe methods, indicates an origin in, or north of, the Indus suture zone. This observation requires that India-Asia collision had occurred by this time and that an Indus River system, feeding material from the suture zone into the basin, was active soon after collision. Pleistocene provenance was similar to that during the Eocene, albeit with greater contribution from the Karakoram. A mass balance of the erosional record on land with deposition in the fan and associated basins suggests that only ~40% of the Neogene sediment in the fan is derived from the Indian plate.

Response of nannoplankton to major changes in sea-surface temperature at the South Chatham Rise, southeastern New Zealand

A high-resolution history of paleoceanographic changes in the subpolar waters of the southern margin of the Subtropical Convergence Zone during the last 130 kyr, is present in foraminiferal assemblages of DSDP Site 594. The foraminifera indicate that sea-surface temperatures during the Last Interglacial Climax were warmer than today, and that between substage 5d through to the end of isotope stage 2, temperatures were mostly cooler than Holocene temperatures. The paleotemperatures suggest that (1) the Subtropical Convergence was located over the site during substage 5e, later moving further north, then moving southwards to near the site during the Holocene, and (2) the Polar Front was positioned over the Site during glacial stages 6, 4, 2 and possibly parts of stage 3. Several major events are indicated by the nannofloral assemblages during these large changes in sea-surface temperature and associated reorganization of ocean circulation. First, the time-progressive trends between E. huxleyi and medium to large Gephyrocupsa are unique to this site, with E. huxleyi dominating over medium Gephyrocupsa during stages 5c-a, middle part of stage 4 and after the middle point of stage 3. This unusual trend may (at least partly) be caused by the shift of the Polar Front across the site. Second, upwelling flora (E. huxleyi and small placoliths) increase in abundance during stages 1, 3 and 5, suggesting that upwelling or disturbance of water stratification took place during the interglacials. Thirdly, there are no significant differences between the distribution patterns of the various morphotypes of medium to large Gephyrocupsu, and the combined value of all medium Gephyrocupsu increases in abundance during glacials (stages 2 and 4 and the end of stage 6), similar to the abundance trends in benthic foraminifera. Finally, subordinate nannofossil taxa also show distinctive climatic trends during the last glacial cycle: (1) Syrucosphaera spp. are present in increased abundance during warmer extremes in climate (substages 5e, 5a, and stage 1); (2) Coccolithus pelagicus and Culcidiscus leptoporus dominate the subordinate nannofossil taxa, and their relative proportions seem to provide a useful paleoceanographic index, with C. pelagicus dominating when the Polar Front Zone is over the site (stages 6, 4 and 2), whilst C. leptoporus is relatively more abundant when the STC is positioned over the site (stages 1 and 5e). Increased abundance of C. pelagicus also can indicate intensified coastal upwelling.

Os and Re concentrations and Os isotope ratios for DSDP and ODP metalliferous sediments from the Pacific Ocean (Table 1)

We report new 187Os/186Os data and Re and Os concentrations in metalliferous sediments from the Pacific to construct a composite Os isotope seawater evolution curve over the past 80 m.y. Analyses of four samples of upper Cretaceous age yield 187Os/186Os values of between 3 and 6.5 and 187Re/186Os values below 55. Mass balance calculations indicate that the pronounced minimum of about 2 in the Os isotope ratio of seawater at the K-T boundary probably reflects the enormous input of cosmogenic material into the oceans by the K-T impactor(s). Following a rapid recovery to 187Os/186Os of 3.5 at 63 Ma, data for the early and middle part of the Cenozoic show an increase in 187Os/186Os to about 6 at 15 Ma. Variations in the isotopic composition of leachable Os from slowly accumulating metalliferous sediments show large fluctuations over short time spans. In contrast, analyses of rapidly accumulating metalliferous carbonates do not exhibit the large oscillations observed in the pelagic clay leach data. These results together with sediment leaching experiments indicate that dissolution of non-hydrogenous Os can occur during the hydrogen peroxide leach and demonstrate that Os data from pelagic clay leachates do not always reflect the Os isotopic composition of seawater. New data for the late Cenozoic further substantiate the rapid increase in the 187Os/186Os of seawater during the past 15 Ma. We interpret the correlation between the marine Sr and Os isotope records during this time period as evidence that weathering within the drainage basin of the Ganges-Brahmaputra river system is responsible for driving seawater Sr and Os toward more radiogenic isotopic compositions. The positive correlation between 87Sr/86Sr and U concentration, the covariation of U and Re concentrations, and the high dissolved Re, U and Sr concentrations found in the Ganges-Brahmaputra river waters supports this interpretation. Accelerating uplift of many orogens worldwide over the past 15 Ma, especially during the last 5 Ma, could have contributed to the rapid increase in 187Os/186Os from 6 to 8.5 over the past 15 Ma. Prior to 15 Ma the marine Sr and Os record are not tightly coupled. The heterogeneous distribution of different lithologies within eroding terrains may play an important role in decoupling the supplies of radiogenic Os and Sr to the oceans and account for the periods of decoupling of the marine Sr and Os isotope records.