Stable isotope record of ODP Sites 108-658 and 108-659

Ocean Drilling Program Site 658, cored below a major upwelling cell offshore Cap Blanc, contains a largely undisturbed hemipelagic sediment section spanning the Brunhes Chron and the early Quaternary and late Pliocene. The companion Site 659 recovered a complete and undisturbed Neogene profile further offshore that serves as a nonupwelling pelagic reference section. Oxygen and carbon isotope ratios in benthic (C. wuellerstorfi and in part Uvigerina sp.) and planktonic foraminifers (G. inflata) provide a climatic record of high resolution for the Brunhes Chron. At Site 658 the record extends back to the early Pleistocene and late Pliocene. The standard oxygen isotope record of the last 730,000 yr is markedly refined by a well-documented high-frequency variation (e.g., by a new "aborted" ice age at stage 13.2 and by Younger-Dryas style climatic setbacks during most terminations). In the late Pliocene, the numerical oxygen isotope stage taxonomy was extended back to stage 137 about 3.3 Ma ago. In comparison with published records, stage 114 at 2.7 Ma represents the first major glaciation event, when 18O was short-term enriched up to a middle Pleistocene glacial d18O level. About 3.17 Ma ago (stage 133), the interglacial oxygen isotope values of C. wuellerstorfi started to increase by 0.5 per mil until 2.7 Ma and then remained largely constant until the Holocene. Based on the d13C difference between C. wuellerstorfi and G. inflata, the dissolved CO2 in the ambient bottom water of Site 658 was dominated by the flux of particulate carbon from the overlying upwelling cell during the last 630,000 yr. In contrast, the advection of (upper) North Atlantic Bottom Water dominated in the control of the local CO2 content during the early Pleistocene and late Pliocene.

A new modern analogue reference dataset and its application to the 430-kyr pollen record from Lake Biwa

This study presents a newly compiled dataset of modern pollen and climate data from 798 sites across Japan and the Russian Far East. This comprehensive reference dataset combined with the modern analogue technique (MAT) provides a powerful tool for pollen-based reconstruction of the Quaternary Northwest Pacific climate. Pollen-derived reconstruction of the modern climate at the reference pollen-sampling sites matches well with the estimated modern climate values (R2 values vary between 0.79 and 0.95, and RMSEP values vary between 5.8 and 9.7% of the modern climatic range for all nine tested variables). The successful testing of the method encourages its application to the fossil pollen records. We used a coarse-resolution pollen record from Lake Biwa to reconstruct glacial-interglacial climate dynamics in central Japan since ~438 kyr and compared it to the earlier reconstruction based on a less representative reference dataset. The current and earlier results consistently demonstrate that the coldest glacial intervals experienced pronounced cooling in winter and moderate cooling in summer, supporting the growth of cool mixed forest (COMX) where warm mixed forest (WAMX) predominates today. During the last glacial, maximum (~24 kyr BP) mean temperatures of the coldest (MTCO) and warmest (MTWA) month were about -13 °C (RMSEP = 2.34 °C) and 21 °C (RMSEP = 1.66 °C) respectively, and annual precipitation (PANN) was about 800 mm (RMSEP = 158.06 mm). During the thermal optimums of the interglacial intervals, the temperatures of the coldest and warmest month were above 0 °C and 25 °C respectively, leading to the reconstruction of WAMX and temperate conifer forest (TECO). Although both these vegetation types grow in the southern part of Japan today, WAMX requires warmer space. The presence of WAMX during marine isotope stages (MIS) 11 and 1, and its absence during MIS 9 and MIS 5 contradict the marine isotope and Antarctic ice records, suggesting that the latter two interglacials were the warmest of the last 800 kyr. The apparent contradiction allows at least three different explanations including low temporal resolution of the pollen record; different trends in CO2 concentrations during 'short' and 'long' interglacials; and regional climate variability and non-linear response of different regions to the global forcing. More definitive conclusions will be possible on the basis of forthcoming high-resolution pollen records from central Japan.

Isotope characterisation of Villars Cave stalagmites

We present a new mid-latitude speleothem record of millennial-scale climatic variability during OIS3 from the Villars Cave that, combined with former published contemporaneous samples from the same cave, gives a coherent image of the climate variability in SW-France between ~55 ka and ~30 ka. The 0.82 m long stalagmite Vil-stm27 was dated with 26 TIMS U-Th analyses and its growth curve displays variations that are linked with the stable isotopes, both controlled by the climatic conditions. It consists in a higher resolved replicate of the previously published Vil-stm9 and Vil-stm14 stalagmites where Dansgaard-Oeschger (DO) events have been observed. The good consistency between these three stalagmites and the comparison with other palaoeclimatic reconstructions, especially high resolution pollen records (ODP 976 from the Alboran Sea, Monticchio Lake record from southern Italy) and the nearby MD04-2845 Atlantic Ocean record, permits to draw a specific climatic pattern in SW-France during the OIS3 and to see regional differences between these sites. Main features of this period are: 1) warm events corresponding to Greenland Interstadials (GIS) that are characterized by low speleothem d13C, high temperate pollen percentages, warm temperatures and high humidity; among these events, GIS#12 is the most pronounced one at Villars characterized by an abrupt onset at ~46.6 ka and a duration of about 2.5 ka. The other well individualized warm event coincides with GIS#8 which is however much less pronounced and occurred during a cooler period as shown by a lower growth rate and a higher d13C; 2) cold events corresponding to Greenland Stadials (GS) that are clearly characterized by high speleothem d13C, low temperate pollen abundance, low temperature and enhanced dryness, particularly well expressed during GS coinciding with Heinrich events H5 and H4. The main feature of the Villars record is a general cooling trend between the DO#12 event ~45.5 ka and the synchronous stop of the three stalagmites at ~30 ka ±1, with a first well marked climatic threshold at ~41 ka after which the growth rate and the diameter of all stalagmites slows down significantly. This climatic evolution differs from that shown at southern Mediterranean sites where this trend is not observed. The ~30 ka age marks the second climatic threshold after which low temperatures and low rainfalls prevent speleothem growth in the Villars area until the Lateglacial warming that occurred at ~16.5 ± 0.5 ka. This 15 ka long hiatus, as the older Villars growth hiatus that occurred between 67.4 and 61 ka, are linked to low sea levels, reduced ocean circulation and a southward shift of the Polar Front that likely provoked local permafrost formation. These cold periods coincide with both low summer 65°N insolation, low atmospheric CO2 concentration and large ice sheets development (especially the Fennoscandian).

Stable isotope record and planktonic foraminiferal fragmentation of the late Pleistocene section of ODP Hole 134-828A

A high-resolution record of foraminiferal fragmentation (a dissolution indicator) for the last 250 k.y. (isotopic Stages 1 to 7) is identified in the upper 61.9 m of Ocean Drilling Program (ODP) Hole 828A, west Vanuatu. This record is comparable in detail to the atmospheric CO2 record and the d18O stack. Phase shifts between preservation spikes and maximum ice volumes (d18O of Globigerinoides sacculifer) are analogous to those on Ontong Java Plateau. Mass spectrometer (AMS14C) dating of a sample taken at the base of dissolution cycle B1 and the position of the last glacial maximum indicates a lag in time of ~8 k.y. in the Vanuatu region for the last glacial termination. When dissolution spikes are compared with minimum ice volumes there is no phase shift for the last two glacial terminations. The difference between Vanuatu and Ontong Java Plateau may be explained by local CO2 sinks and the interplay between intermediate and deep water masses. Terrigenous input increasingly affected sediment of Hole 828A on the North d'Entrecasteaux Ridge (NDR) as it approached Espiritu Santo Island. Mud and silt suspended in mid-water flows become important after 125 ka, while turbidites bypass the New Hebrides Trench only towards the last glacial maximum (LGM). Terrigenous supply seems to affect the lysocline profile that changed from an "open ocean" to a "near continent" type, thus favoring dissolution. Fragmentation of planktonic foraminifers is a more sensitive indicator of lysocline variations than is foraminiferal susceptibility to dissolution, the foraminiferal dissolution index, the abundance of benthic foraminifers, or CaCO3 content. A modern foraminiferal lysocline for the neighboring area (between 10°S and 30°S, and 160°E and 180°E) is found at 3.1 km below sea level, compared to west Vanuatu where it is shallower. The past lysocline level was deeper than 3086 m during intervals of dissolution minima, and ranged from ~2550 to 3000 m during intervals of dissolution maxima. The high sedimentation rates (in the order of 10 to 50 cm/k.y.) found in Hole 828A offer a great potential for future high-resolution studies either in this hole or other western localities along the NDR. Areas of high sedimentation near continental regions have been discarded for paleoceanographic and/or paleoclimatic studies. Nonetheless, conditions analogous to those found in Hole 828A are expected to occur in many trench areas around the world where mid-water flows have preserved as yet undiscovered fine high-resolution sedimentary records.

Geochemistry and Nd-Pb isotope ratios of sediments from the Lesser Antilles arc

The Lesser Antilles arc is a particularly interesting island arc because it is presently very active, it is located perpendicular to the South American continent and its chemical and isotopic compositions display a strong north-south gradient. While the presence in the south of a thick pile of sedimentary material coming from the old South American continent has long been suspected to explain the geochemical gradient, previous studies failed to demonstrate unambiguously a direct link between the arc lava compositions and the subducted sediment compositions. Here, we present new Nd, Sm, Th, U and Pb concentrations and Nd-Pb isotopic data for over 60 sediments from three sites located in the fore arc region of the Lesser Antilles arc. New data for DSDP Site 543 drill core located east of Dominica Island complement the data published by White et al. (1985, doi:10.1016/0016-7037(85)90082-1) and confirm their relatively uniform isotopic compositions (i.e., 206Pb/204Pb between 19.13 and 19.53). In contrast, data obtained on DSDP Site 144 located further south, on the edge of the South American Rise and on sediments from Barbados Island are much more variable (206Pb/204Pb ranges from 18.81 to 27.69). The very radiogenic Pb isotopic compositions are found in a 60 m thick black shale unit, which has no age equivalent in the Site 543 drill core. We interpret the peculiar composition of the southern sediments as being due to two factors, (a) the proximity of the South American craton, which contributes coarse grain old detrital material that does not travel far from the continental shelf, and (b) the presence of older sediments including the thick black shale unit formed during Oceanic Anoxic events 2 and 3. The north-south isotopic change known along the Lesser Antilles arc can be explained by the observed geographical changes in the composition of the subducted sediments. About 1% contamination of the mantle wedge by Site 543 sediments explains the composition of the northern islands while up to 10% sediments like those of Site 144 is required in the source of the southern island lavas. The presence of black shales in the subducted pile provides a satisfactory explanation for the very low Delta8/4 values that characterize the Lesser Antilles arc.

Middle to late Miocene stable isotope record of benthic foraminifera from ODP Site 184-1146

We present high-resolution (2-3 kyr) benthic foraminiferal stable isotopes in a continuous, well-preserved sedimentary archive from the West Pacific Ocean (Ocean Drilling Program Site 1146), which track climate evolution in unprecedented resolution over the period 12.9 to 8.4 Ma. We developed an astronomically tuned chronology over this interval and integrated our new records with published isotope data from the same location to reconstruct long-term climate and ocean circulation development between 16.4 and 8.4 Ma. This extended perspective reveals that the long eccentricity (400 kyr) cycle is prominently encoded in the d13C signal over most of the record, reflecting long-term fluctuations in the carbon cycle. The d18O signal closely follows variations in short eccentricity (100 kyr) and obliquity (41 kyr). In particular, the obliquity cycle is prominent from ~14.6 to 14.1 Ma and from ~9.8 to 9.2 Ma, when high-amplitude variability in obliquity is congruent with low-amplitude variability in short eccentricity. The d18O curve is additionally characterized by a series of incremental steps at ~14.6, 13.9, 13.1, 10.6, 9.9, and 9.0 Ma, which we attribute to progressive deep water cooling and/or glaciation episodes following the end of the Miocene climatic optimum. On the basis of d18O amplitudes, we find that climate variability decreased substantially after ~13 Ma, except for a remarkable warming episode at ~10.8-10.7 Ma at peak insolation during eccentricity maxima (100 and 400 kyr). This transient warming, associated with a massive negative carbon isotope shift, is reminiscent of intense global warming events at eccentricity maxima during the Miocene climatic optimum.

Eocene sedimentary calcium carbonate contents and stable isotope composition of benthic foraminifera

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

Planktonic foraminifera oxygen isotope and trace metal data ODP Hole 172-1056A, 46-54 ka BP

Salinity increase in the subtropical gyre system may have pre-conditioned the North Atlantic Ocean for a rapid return to stronger overturning circulation and high-latitude warming following meltwater events during the Last Glacial period. Here we investigate the Gulf Stream - subtropical gyre system properties over Dansgaard-Oeschger (DO) cycles 14 to 12, including Heinrich ice-rafting event 5. During the Holocene and Last Glacial Maximum a positive gradient in surface dwelling planktonic foraminifera d18O (Globigerinoides ruber) can be observed between the Gulf Stream and subtropical gyre, due to decreasing temperature, increasing salinity, and a change from summer to year-round occurrence of G. ruber. We assess whether this gradient was a common feature during stadial-interstadial climate oscillations of Marine Isotope Stage 3, by comparing existing G. ruber d18O from ODP Site 1060 (subtropical gyre location) and new data from ODP Site 1056 (Gulf Stream location) between 54 and 46 ka. Our results suggest that this gradient was largely absent during the period studied. During the major warm DO interstadials 14 and 12 we infer a more zonal and wider Gulf Stream, influencing both ODP Sites 1056 and 1060. A Gulf Stream presence during these major interstadials is also suggested by the large vertical d18O gradient between shallow dwelling planktonic foraminifera species, especially G. ruber, and the deep dwelling species Globorotalia inflata at site 1056, which we associate with strong summer stratification and Gulf Stream presence. A major reduction in this vertical d18O gradient from 51 ka until the end of Heinrich event 5 at 48.5 ka suggests site 1056 was situated within the subtropical gyre in this mainly cold period, from which we infer a migration of the Gulf Stream to a position nearer to the continental shelf, indicative of a narrower Gulf Stream with possibly reduced transport.

Marine osmium isotope record across the Cretaceous-Tertiary boundary

A composite late Maastrichtian (65.5 to 68.5 Ma) marine osmium (Os) isotope record, based on samples from the Southern Ocean (ODP Site 690), the Tropical Pacific Ocean (DSDP Site 577), the South Atlantic (DSDP Site 525) and the paleo-Tethys Ocean demonstrates that subaerially exposed pelagic carbonates can record seawater Os isotope variations with a fidelity comparable to sediments recovered from the seafloor. New results provide robust evidence of a 20% decline in seawater 187Os/188Os over a period of about 200 kyr early in magnetochron C29r well below the Cretaceous-Paleogene Boundary (KPB), confirming previously reported low-resolution data from the South Atlantic Ocean. New results also confirm a second more rapid decline in 187Os/188Os associated with the KPB that is accompanied by a significant increase in Os concentrations. Complementary platinum (Pt) and iridium (Ir) concentration data indicate that the length scale of diagenetic remobilization of platinum group elements from the KPB is less than 1 m and does not obscure the pre-KPB decline in 187Os/188Os. Increases in bulk sediment Ir concentrations and decreases in bulk carbonate content that coincide with the Os isotope shift suggest that carbonate burial flux may have been lower during the initial decline in 187Os/188Os. We speculate that diminished carbonate burial rate may have been the result of ocean acidification caused by Deccan volcanism.

Strontium and Oxygen isotope compositions from DSDP Holes 54-424A, 54-424B and 70-509B

Strontium and O isotope compositions of green clay minerals from sediment cores of three boreholes drilled into (sites 424A and 509B) and close to a hydrothermal mound (site 424B) near the Galapagos Spreading Center (DSDP Legs 54 and 70) were determined. The green clays consist mostly of a transition from Fesmectite (nontronite) to glauconite. 87Sr/86Sr ratios were measured on clay size-fractions after gentle acid leaching and on the recovered leachates from different samples. The 87Sr/86Sr ratios of the clay residues from both the 424A and B sites are well below the modern seawater value, which points consistently to precipitation from hydrothermal fluids that contained variable amounts of seawater, even away from mound. However, most of the clay residues from mound site 509B have 87Sr/86Sr ratios significantly above the seawater value, suggesting the occurrence of a detrital component together with the new authigenic particles. The clay minerals of the hydrothermal mound are mixed with detrital components, and that of the sample taken outside but near the mound as a reference for the surrounding oceanic environment, yields a hydrothermal signature. Crystallization temperatures of the clays range from 32 to 63 °C assuming a d18O value of +2.2 per mil for the mineralizing fluids. Hydrothermal fluids generated in the underlying oceanic crust, mixed in varied proportions with ambient seawater and migrated into beds of the mound in a sequence of recurrent processes that ultimately resulted in the formation of the observed clay minerals. No significant temperature differences were detected for crystallization of the K-rich glauconite and K-depleted nontronite. The 87Sr/86Sr ratios of the Sr leached off the clay particles are near the value of modern seawater, inferring a progressive replacement of the hydrothermal fluids by seawater in the pore space of the mound sediments.

Natural gamma ray and stable isotope record of ODP Site 181-1119

Ocean Drilling Program Site 1119 is ideally located to intercept discharges of sediment from the mid-latitude glaciers of the New Zealand Southern Alps. The natural gamma ray signal from the site's sediment core contains a history of the South Island mountain ice cap since 3.9 million years ago (Ma). The younger record, to 0.37 Ma, resembles the climatic history of Antarctica as manifested by the Vostok ice core. Beyond, and back to the late Pliocene, the record may serve as a proxy for both mid-latitude and Antarctic polar plateau air temperature. The gamma ray signal, which is atmospheric, also resembles the ocean climate history represented by oxygen isotope time series.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.