B-isotope and Mg/Ca ratios of planktonic foraminifera from ODP Hole 108-668B and esimates of salinity, alkalinity and pCO2 (Table 1)

Knowledge of past atmospheric pCO2 is important for evaluating the role of greenhouse gases in climate forcing. Ice core records show the tight correlation between climate change and pCO2, but records are limited to the past ~900 kyr. We present surface ocean pH and pCO2 data, reconstructed from boron isotopes in planktonic foraminifera over two full glacial cycles (0-140 and 300-420 kyr). The data co-vary strongly with the Vostok pCO2-record and demonstrate that the coupling between surface ocean chemistry and the atmosphere is recorded in marine archives, allowing for quantitative estimation of atmospheric pCO2 beyond the reach of ice cores.

(Table A1) Isotopic ratios of Sr, Nd, and Pb for sediments in ODP Leg 138 sites

The provenance of eolian dust supplied to deep-sea sediments has the potential to offer insights into changes in past atmospheric circulation. Specifically, measuring temporal changes in dust provenance can shed light on changes in the mean position of the Intertropical Convergence Zone (ITCZ), a region acting as a barrier separating wind-blown material derived from northern versus southern hemisphere sources. Here we have analyzed Nd, Sr, and Pb isotope ratios in the operationally-defined detrital component extracted from deep-sea sediments in the eastern equatorial Pacific (EEP) along a meridional transect at 110°W from 3°S to 7°N (ODP Leg 138, sites 848-853). Sr isotope results show that barite Sr has a significant influence on 87Sr/86Sr isotope ratios of samples in the upwelling zone of the EEP. However, sites located >3° or more away from the equator (sites 852 and 853) are believed to not be affected by barite Sr and provide useful detrital Sr signals. 208Pb/206Pb and 207Pb/206Pb ratios in all cores fall into the Pb-isotope space of five potential dust sources (Asia, North and Central/South America, Sahara, and Australia), with no distinct isotopic fingerprinting of the dominant source(s). epsilon-Nd values were most valuable for discerning detrital source provenance, and their values at all sites, ranging from ~5.46 to ~3.25, were more unradiogenic for sediments deposited during the last glacial than for those deposited during the Holocene. There are distinct latitudinal trends in the epsilon-Nd values, with more radiogenic values further south and less radiogenic values further north, excluding site 848. This distinction holds true for both Holocene and last glacial periods. For the most southerly site, 848, we invoke, for the first time, a distinct southern hemisphere Australian source as being responsible for the unradiogenic Nd isotope ratios. Both average last glacial and Holocene epsilon-Nd values show similar sharp gradients along the transect between 5.29°N and 2.77°N, suggesting little movement of the glacial ITCZ in the EEP. However, during the deglacial, this gradient is stronger and shifted further north between 5.29°N and 7.21°N, suggesting a more northerly, possibly stronger, deglacial ITCZ.

Alkenone and boron based Oligocene pCO2 records

The Pliocene period is the most recent time when the Earth was globally significantly (~3°C) warmer than today. However, the existing pCO2 data for the Pliocene are sparse and there is little agreement between the various techniques used to reconstruct palaeo-pCO2. Moreover, the temporal resolution of the published records does not allow a robust assessment of the role of declining pCO2 in the intensification of the Northern Hemisphere Glaciation (INHG) and a direct comparison to other proxy records are lacking. For the first time, we use a combination of foraminiferal (delta11B) and organic biomarker (alkenone-derived carbon isotopes) proxies to determine the concentration of atmospheric CO2 over the past 5 Ma. Both proxy records show that during the warm Pliocene pCO2 was between 330 and 400 ppm, i.e. similar to today. The decrease to values similar to pre-industrial times (275-285 ppm) occurred between 3.2 Ma and 2.8 Ma - coincident with the INHG and affirming the link between global climate, the cryosphere and pCO2.

Stable isotope record, opal accumulation rate and geochemistry of middle Miocene sediments of IODP Site 321-U1338

During the middle Miocene, Earth's climate transitioned from a relatively warm phase (Miocene climatic optimum) into a colder mode with re-establishment of permanent ice sheets on Antarctica, thus marking a fundamental step in Cenozoic cooling. Carbon sequestration and atmospheric CO2 drawdown through increased terrestrial and/or marine productivity have been proposed as the main drivers of this fundamental transition. We integrate high-resolution (1-3 k.y.) benthic stable isotope data with XRF-scanner derived biogenic silica and carbonate accumulation estimates in an exceptionally well-preserved sedimentary archive, recovered at Integrated Ocean Drilling Program Site U1338, to reconstruct eastern equatorial Pacific productivity variations and to investigate temporal linkages between high- and low-latitude climate change over the interval 16-13 Ma. Our records show that the climatic optimum (16.8-14.7 Ma) was characterized by high amplitude climate variations, marked by intense perturbations of the carbon cycle. Episodes of peak warmth at (southern hemisphere) insolation maxima coincided with transient shoaling of the carbonate compensation depth and enhanced carbonate dissolution in the deep ocean. A switch to obliquity-paced climate variability after 14.7 Ma concurred with a general improvement in carbonate preservation and the onset of stepwise global cooling, culminating with extensive ice growth over Antarctica at ~13.8 Ma. We find that two massive increases in opal accumulation at ~14.0 and ~13.8 Ma occurred just before and during the final and most prominent cooling step, supporting the hypothesis that enhanced siliceous productivity in the eastern equatorial Pacific contributed to CO2 drawdown.

Alkenones and TEX86 temperature proxies for core-top sediments from the southern Adriatic Sea

The Mediterranean Sea is at the transition between temperate and tropical air masses and as such of importance for studying climate change. The Gulf of Taranto and adjacent SW Adriatic Sea are at the heart of this region. Their sediments are excellently suited for generating high quality environmental records for the last millennia with a sub-decadal resolution. The quality of these records is dependent on a careful calibration of the transfer functions used to translate the sedimentary lipid signals to the local environment. Here, we examine and calibrate the UK'37 and TEX86 lipid-based temperature proxies in 48 surface sediments and relate these to ambient sea surface temperatures and other environmental data. The UK'37-based temperatures in surface sediments reflect winter/spring sea surface temperatures in agreement with other studies demonstrating maximum haptophyte production during the colder season. The TEX86-based temperatures for the nearshore sites also reflect winter sea surface temperatures. However, at the most offshore sites, they correspond to summer sea surface temperatures. Additional lipid and environmental data including the distribution of the BIT index and remote-sensed chlorophyll-a suggest a shoreward increase of the impact of seasonal and spatial variability in nutrients and control of planktonic archaeal abundance by primary productivity, particle loading in surface waters and/or overprint by a cold-biased terrestrial TEX86 signal. As such the offshore TEX86 values seem to reflect a true summer signal to the effect that offshore UK'37 and TEX86 reconstruct winter and summer temperature, respectively, and hence provide information on the annual temperature amplitude.

Mass balance and acceleration of the Greenland ice sheet

Within the last decade, the Greenland ice sheet (GrIS) and its surroundings have experienced record high surface temperatures (Mote, 2007, doi:10.1029/2007GL031976; Box et al., 2010), ice sheet melt extent (Fettweis et al., 2011, doi:10.5194/tc-5-359-2011) and record-low summer sea-ice extent (Nghiem et al., 2007, doi:10.1029/2007GL031138). Using three independent data sets, we derive, for the first time, consistent ice-mass trends and temporal variations within seven major drainage basins from gravity fields from the Gravity Recovery and Climate Experiment (GRACE; Tapley et al., 2004, doi:10.1029/2004GL019920), surface-ice velocities from Inteferometric Synthetic Aperture Radar (InSAR; Rignot and Kanagaratnam, 2006, doi:10.1126/science.1121381) together with output of the regional atmospheric climate modelling (RACMO2/ GR; Ettema et al., 2009, doi:10.1029/2009GL038110), and surface-elevation changes from the Ice, cloud and land elevation satellite (ICESat; Sorensen et al., 2011, doi:10.5194/tc-5-173-2011). We show that changing ice discharge (D), surface melting and subsequent run-off (M/R) and precipitation (P) all contribute, in a complex and regionally variable interplay, to the increasingly negative mass balance of the GrIS observed within the last decade. Interannual variability in P along the northwest and west coasts of the GrIS largely explains the apparent regional mass loss increase during 2002-2010, and obscures increasing M/R and D since the 1990s. In winter 2002/2003 and 2008/2009, accumulation anomalies in the east and southeast temporarily outweighed the losses by M/R and D that prevailed during 2003-2008, and after summer 2010. Overall, for all basins of the GrIS, the decadal variability of anomalies in P, M/R and D between 1958 and 2010 (w.r.t. 1961-1990) was significantly exceeded by the regional trends observed during the GRACE period (2002-2011).

Stable carbon and oxygen isotope compositions of bulk rock samples analyzed from the Zumaia section

Astronomical tuning of sedimentary records to precise orbital solutions has led to unprecedented resolution in the geological time scale. However, the construction of a consistent astronomical time scale for the Paleocene is controversial due to uncertainties in the recognition of the exact number of 405-kyr eccentricity cycles and accurate correlation between key records. Here, we present a new Danian integrated stratigraphic framework using the land-based Zumaia and Sopelana hemipelagic sections from the Basque Basin and deep-sea records drilled during Ocean Drilling Program (ODP) Legs 198 (Shatsky Rise, North Pacific) and 208 (Walvis Ridge, South Atlantic) that solves previous discrepancies. The new coherent stratigraphy utilises composite images from ODP cores, a new whole-rock d13C isotope record at Zumaia and new magnetostratigraphic data from Sopelana. We consistently observe 11 405-kyr eccentricity cycles in all studied Danian successions. We achieve a robust correlation of bioevents and stable isotope events between all studied sections at the ~100-kyr short-eccentricity level, a prerequisite for paleoclimatic interpretations. Comparison with and subsequent tuning of the records to the latest orbital solution La2011 provides astronomically calibrated ages of 66.022 ± 0.040 Ma and 61.607 ± 0.040 Ma for the Cretaceous-Paleogene (K-Pg) and Danian-Selandian 105 (D-S) boundaries respectively. Low sedimentation rates appear common in all records in the mid-Danian interval, including conspicuous condensed intervals in the oceanic records that in the past have hampered the proper identification of cycles. The comprehensive interbasinal approach applied here reveals pitfalls in time scale construction, filtering techniques in particular, and indicates that some caution and scrutiny has to be applied when building orbital chronologies. Finally, the Zumaia section, already hosting the Selandian Global Boundary Stratotype Section and Point (GSSP), could serve as the global Danian unit stratotype in the future.