(Table 1) Stable isotope records of clathrite and water

Aragonitic clathrites are methane-derived precipitates that are found at sites of massive near-seafloor gas hydrate (clathrate) accumulations at the summit of southern Hydrate Ridge, Cascadia margin. These platy carbonate precipitates form inside or in proximity to gas hydrate, which in our study site currently coexists with a fluid that is highly enriched in dissolved ions as salts are excluded during gas hydrate formation. The clathrites record the preferential incorporation of 18O into the hydrate structure and hence the enrichment of 16O in the surrounding brine. We measured d18O values as high as 2.27 per mil relative to Peedee belemnite that correspond to a fluid composition of -1.18 per mil relative to standard mean ocean water. The same trend can be observed in Ca isotopes. Ongoing clathrite precipitation causes enrichment of the 44Ca in the fluid and hence in the carbonates. Carbon isotopes confirm a methane source for the carbonates. Our triple stable isotope approach that uses the three main components of carbonates (Ca, C, O) provides insight into multiple parameters influencing the isotopic composition of the pore water and hence the isotopic composition of the clathrites. This approach provides a tool to monitor the geochemical processes during clathrate and clathrite formation, thus recording the evolution of the geochemical environment of gas hydrate systems.

Stable isotope record of sediment core MD98-2177

We present a reconstruction of El Niño Southern Oscillation (ENSO) variability spanning the Medieval Climate Anomaly (MCA, A.D. 800-1300) and the Little Ice Age (LIA, A.D. 1500-1850). Changes in ENSO are estimated by comparing the spread and symmetry of d18O values of individual specimens of the thermocline-dwelling planktonic foraminifer Pulleniatina obliquiloculata extracted from discrete time horizons of a sediment core collected in the Sulawesi Sea, at the edge of the western tropical Pacific warm pool. The spread of individual d18O values is interpreted to be a measure of the strength of both phases of ENSO while the symmetry of the d18O distributions is used to evaluate the relative strength/frequency of El Niño and La Niña events. In contrast to previous studies, we use robust and resistant statistics to quantify the spread and symmetry of the d18O distributions; an approach motivated by the relatively small sample size and the presence of outliers. Furthermore, we use a pseudo-proxy approach to investigate the effects of the different paleo-environmental factors on the statistics of the d18O distributions, which could bias the paleo-ENSO reconstruction. We find no systematic difference in the magnitude/strength of ENSO during the Northern Hemisphere MCA or LIA. However, our results suggest that ENSO during the MCA was skewed toward stronger/more frequent La Niña than El Niño, an observation consistent with the medieval megadroughts documented from sites in western North America.

(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.

Stable isotope record of Globorotalia inflata

We sampled the upper water column for living planktic foraminifera along the SW-African continental margin. The species Globorotalia inflata strongly dominates the foraminiferal assemblages with an overall relative abundance of 70-90%. The shell delta18O and delta13C values of G. inflata were measured and compared to the predicted oxygen isotope equilibrium values (delta18O(eq)) and to the carbon isotope composition of the total dissolved inorganic carbon (delta13C(DIC)) of seawater. The delta18O of G. inflata reflects the general gradient observed in the predicted delta18O(eq) profile, while the delta13C of G. inflata shows almost no variation with depth and the reflection of the delta13C(DIC) in the foraminiferal shell seems to be covered by other effects. We found that offsets between delta18O(shell) and predicted delta18O(eq) in the surface mixed layer do not correlate to changes in seawater [CO3[2-]]. To calculate an isotopic mass balance of depth integrated growth, we used the oxygen isotope composition of G. inflata to estimate the fraction of the total shell mass that is grown within each plankton tow depth interval of the upper 500 m of the water column. This approach allows us to calculate the DELTA delta13C(interval added-DIC); i.e. the isotopic composition of calcite that was grown within a given depth interval. Our results consistently show that the DELTA delta13C(IA-DIC) correlates negatively with in situ measured [CO3[2-]] of the ambient water. Using this approach, we found DELTA delta13C(IA-DIC)/[CO3[2-]] slopes for G. inflata in the large size fraction (250-355 µm) of -0.013 per mil to 0.015 per mil (µmol/kg)**-1 and of -0.013 per mil to 0.017 per mil (µmol/kg)**-1 for the smaller specimens (150-250 µm). These slopes are in the range of those found for other non-symbiotic species, such as Globigerina bulloides, from laboratory culture experiments. Since the DELTA delta13C(IA-DIC)/[CO3[2-]] slopes from our field data are nearly identical to the slopes established from laboratory culture experiments we assume that the influence of other effects, such as temperature, are negligibly small. If we correct the delta13C values of G. inflata for a carbonate ion effect, the delta13C(shell) and delta13C(DIC) are correlated with an average offset of 2.11.

Mineralogical and stable isotope analyses of carbonates from the eastern Black Sea

Authigenic carbonates associated with cold seeps provide valuable archives of changes in the long-term seepage activity. To investigate the role of shallow-buried hydrates on the seepage strength and fluid composition we analysed methane-derived carbonate precipitates from a high-flux hydrocarbon seepage area ("Batumi seep area") located on the south-eastern Black Sea slope in ca. 850 m. In a novel approach, we combined computerized X-ray tomography (CT) with mineralogical and isotope geochemical methods to get additional insights into the three-dimensional internal structure of the carbonate build-ups. X-ray diffractometry revealed the presence of two different authigenic carbonate phases, i.e. pure aragonitic rims associated with vital microbial mats and high-Mg calcite cementing the hemipelagic sediment. As indicated by the CT images, the initial sediment has been strongly deformed, first plastic then brittle, leading to brecciation of the progressively cemented sediment. The aragonitic rims on the other hand, represent a presumably recent carbonate growth phase since they cover the already deformed sediment. The stable oxygen isotope signature indicates that the high-Mg calcite cement incorporated pore water mixed with substantial hydrate water amounts. This points at a dominant role of high gas/fluid flux from decomposing gas hydrates leading to the deformation and cementation of the overlying sediment. In contrast, the aragonitic rims do not show an influence of 18O-enriched hydrate water. The differences in d18O between the presumably recent aragonite precipitates and the older high-Mg cements suggest that periods of hydrate dissociation and vigorous fluid discharge alternated with times of hydrate stability and moderate fluid flow. These results indicate that shallow-buried gas hydrates are prone to episodic decomposition with associated vigorous fluid flow. This might have a profound impact on the seafloor morphology resulting e.g. in the formation of carbonate pavements and pockmark-like structures but might also affect the local carbon cycle.

(Supplementary Table 1) Proportions of methane (C1) to ethane (C2) and the stable isotope composition of methane (d13CCH4) of gas and oil collected by gas bubble sampler during cruise M 114 in the southern Gulf of Mexico

Hydrocarbon seepage is a widespread process at the continental margins of the Gulf of Mexico. We used a multidisciplinary approach, including multibeam mapping and visual seafloor observations with different underwater vehicles to study the extent and character of complex hydrocarbon seepage in the Bay of Campeche, southern Gulf of Mexico. Our observations showed that seafloor asphalt deposits previously only known from the Chapopote Knoll also occur at numerous other knolls and ridges in water depths from 1230 to 3150 m. In particular the deeper sites (Chapopopte and Mictlan knolls) were characterized by asphalt deposits accompanied by extrusion of liquid oil in form of whips or sheets, and in some places (Tsanyao Yang, Mictlan, and Chapopote knolls) by gas emission and the presence of gas hydrates in addition. Molecular and stable carbon isotopic compositions of gaseous hydrocarbons suggest their primarily thermogenic origin. Relatively fresh asphalt structures were settled by chemosynthetic communities including bacterial mats and vestimentiferan tube worms, whereas older flows appeared largely inert and devoid of corals and anemones at the deep sites. The gas hydrates at Tsanyao Yang and Mictlan Knolls were covered by a 5-to-10 cm-thick reaction zone composed of authigenic carbonates, detritus, and microbial mats, and were densely colonized by 1-2 m-long tube worms, bivalves, snails, and shrimps. This study increased knowledge on the occurrences and dimensions of asphalt fields and associated gas hydrates at the Campeche Knolls. The extent of all discovered seepage structure areas indicates that emission of complex hydrocarbons is a widespread, thus important feature of the southern Gulf of Mexico.

CO2 concentration and stable isotope ratios of three Antarctic ice cores covering the period from 149.4 - 1.5 kyr before 1950

We present new d13C measurements of atmospheric CO2 covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in d13C(atm) of 0.5 permil occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in d13C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO2]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of d13C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

Stable isotope record of planktonic foraminifera of sediment core TR163-19

Stable isotope data from eastern equatorial Pacific (EEP) core TR163-19 (2°15'N, 90°57'W, 2348 m) are presented for the surface-dwelling foraminifers Globigerinoides ruber and G. sacculifer and thermocline-dwelling Globorotalia menardii and Neogloboquadrina dutertrei. Using species-specific normalization factors derived from experimental and plankton tow data, we reconstruct a 360 kyr record of water column hydrography across the past three glacial cycles. We demonstrate that G. ruber maintains a mixed layer habitat throughout the entire record, while G. sacculifer records a mixture of thermocline and mixed layer conditions and G. menardii and N. dutertrei record thermocline properties. We conclude that G. sacculifer is not appropriate for paleoceanographic applications in regions with steep vertical hydrographic gradients. Results suggest that this region of the EEP had a thicker mixed layer and deeper d13CDIC boundary between the surface and equatorial undercurrent during the last two glacial periods. A shift in N. dutertrei and G. sacculifer geochemistry prior to ~185 kyr suggests water column structure and chemocline gradients changed, possibly due to a shift in the position of the undercurrent relative to this site. The timing and magnitude of glacial-interglacial d13C variations between species indicates that near-surface carbon chemistry is controlled by changes in productivity, atmospheric circulation, and advected intermediate water sources north of the Antarctic polar front. These results demonstrate that when properly calibrated for species differences, multispecies geochemical data sets can be invaluable for reconstructing water column structure and properties in the past.

Middle Miocene evolution and stable isotope ratios of Globorotalia (Fohsella) in the western equatorial Pacific

Evolution of the planktic foraminiferal lineage Globorotalia (Fohsella) occurred during the Miocene between 23.7 and 11.8 Ma and forms the basis for stratigraphic subdivision of the early middle Miocene (Zones N 10 through N 12). Important morphologic changes within the G. (Fohsella) lineage included a marked increase in test size, a transition from a rounded to an acute periphery, and the development of a keel in later forms. We found that the most rapid changes in morphology of G. (Fohsella) occurred between 13 and 12.7 Ma and coincided with an abrupt increase in the delta18O ratios of shell calcite. Comparison of isotopic results of G. (Fohsella) with other planktic foraminifers indicate that delta18O values of the lineage diverge from surface-dwelling species and approach deep-dwelling species after 13.0 Ma, indicating a change in depth habitat from the surface mixed layer to intermediate depth near the thermocline. Isotopic and faunal evidence suggests that this change in depth stratification was associated with an expansion of the thermocline in the western equatorial Pacific. After adapting to a deeper water habitat at 13.0 Ma, the G. (Fohsella) lineage became extinct abruptly at 11.8 Ma during a period when isotopic and faunal evidence suggest a shoaling of the thermocline. Following the extinction of G. (Fohsella), the ecologic niche of the lineage was filled by the Globorotalia (Menardella) group, which began as a deep-water form and later evolved to an intermediate-water habitat. We suggest that the evolution of G. (Fohsella) and G. (Menardella) were tightly linked to changes in the structure of the thermocline in the western equatorial Pacific.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.

Stable isotope records of corals from the northern Red Sea

We present a study based on X-ray chronologies and the stable isotopic composition of fossil Porites spp. corals from the northern Gulf of Aqaba (Red Sea) covering the mid-Holocene period from 5750 to 4450 14C years BP (before present). The stable oxygen and carbon isotopic compositions of five specimens reveal regular annual periodicities. Compared with modern Porites spp. from the same environment, the average seasonal delta18O amplitude of the fossil corals is higher (by ca. 0.35-0.60?), whereas annual growth rates are lower (by ca. 3.5 to 2 mm/year). This suggests stronger seasonality of sea surface temperatures and increased variability of the oxygen isotopic composition of the sea water due to changes in the precipitation and evaporation regime during the mid-Holocene. Most likely, summer monsoon rains reached the northern end of the Red Sea at that time. Average annual coral growth rates are diminished probably due to an increased input and resuspension of terrestrial debris to the shallow marine environment during more humid conditions. Our results corroborate published reports of paleodata and model simulations suggesting a northward migration of the African monsoon giving rise to increased seasonalities during the mid-Holocene over northeastern Africa and Arabia.

Stable isotope analysis and temperature reconstruction data from DSDP Hole 94-609B and ODP Hole 162-984B

The Greenland ice sheet is accepted as a key factor controlling the Quaternary glacial scenario. However, the origin and mechanisms of major Arctic glaciation starting at 3.15 Ma and culminating at 2.74 Ma are still controversial. For this phase of intense cooling Ravelo et al. proposed a complex gradual forcing mechanism. In contrast, our new submillennial-scale paleoceanographic records from the Pliocene North Atlantic suggest a far more precise timing and forcing for the initiation of northern hemisphere glaciation (NHG), since it was linked to a 2-3 °C surface water warming during warm stages from 2.95 to 2.82 Ma. These records support previous models, claiming that the final closure of the Panama Isthmus (3.0- ~2.5 Ma induced an increased poleward salt and heat transport. Associated strengthening of North Atlantic Thermohaline Circulation and in turn, an intensified moisture supply to northern high latitudes resulted in the build-up of NHG, finally culminating in the great, irreversible climate crash at marine isotope stage G6 (2.74 Ma). In summary, there was a two-step threshold mechanism that marked the onset of NHG with glacial-to-interglacial cycles quasi-persistent until today.

Stable isotope record of planktonic foraminifera from sediment core MD96-2085

A high-resolution (~4-5cm/kyr) giant piston core record (MD962085) retrieved during an IMAGES II-NAUSICAA cruise from the continental slope of the southeast Atlantic Ocean reveals striking variations in planktonic foraminifer faunal abundances and sea-surface temperatures (SST) during the past 600 000 yr. The location and high-quality sedimentary record of the core provide a good opportunity to assess the variability of the Benguela Current system and associated important features of the ocean-climate system in the southeast Atlantic. The planktonic foraminifer faunal abundances of the core are dominated by three assemblages: (1) Neogloboquadrina pachyderma (right coiling) + Neogloboquadrina dutertrei, (2) Globigerina bulloides, and (3) Globorotalia inflata. The assemblage of N. pachyderma (right coiling) + N. dutertrei shows distinctive abundance changes which are nearly in-phase with glacial-interglacial variations. The high abundances of this assemblage are associated with major glacial conditions, possibly representing low SST/high nutrient level conditions in the southwestern Africa margin. In contrast, the G. bulloides and G. inflata assemblages show greater high-frequency abundance change patterns, which are not parallel to the glacial-interglacial changes. These patterns may indicate rapid oceanic frontal movements from the south, and a rapid change in the intensity of the Benguela upwelling system from the east. A single episode of maximum abundances of a polar water species N. pachyderma (left coiling) occurred in the beginning of stage 9 (~340-330 kyr). The event of the maximum occurrence of this species shown in this record may indicate instability in the Benguela coastal upwelling, or the Antarctic polar front zone position. A winter season SST estimate using transfer function techniques for this record shows primarily glacial-interglacial variations. The SST is maximal during the transitions from the major glacial to interglacial stages (Terminations I, II, IV, V), and is associated with the abundance maxima of a warm water species indicator Globigerinoides ruber. Cross-spectral analyses of the SST record and the SPECMAP stack reveal statistically significant concentrations of variance and coherencies in three major orbital frequency bands. The SST precedes changes in the global ice volume in all orbital frequency bands, indicating a dominant southern Hemispheric climate effect over the Benguela Current region in the southeast Atlantic.