Natural remanent magnetization (NRM) from ODP Site 1263 covering magnetochron C20r and C21n and high-resolution bulk carbon isotope (d13C) records from Sites 702 and 1263

To explore cause and consequences of past climate change, very accurate age models such as those provided by the astronomical timescale (ATS) are needed. Beyond 40 million years the accuracy of the ATS critically depends on the correctness of orbital models and radioisotopic dating techniques. Discrepancies in the age dating of sedimentary successions and the lack of suitable records spanning the middle Eocene have prevented development of a continuous astronomically calibrated geological timescale for the entire Cenozoic Era. We now solve this problem by constructing an independent astrochronological stratigraphy based on Earth's stable 405 kyr eccentricity cycle between 41 and 48 million years ago (Ma) with new data from deep-sea sedimentary sequences in the South Atlantic Ocean. This new link completes the Paleogene astronomical timescale and confirms the intercalibration of radioisotopic and astronomical dating methods back through the Paleocene-Eocene Thermal Maximum (PETM, 55.930 Ma) and the Cretaceous-Paleogene boundary (66.022 Ma). Coupling of the Paleogene 405 kyr cyclostratigraphic frameworks across the middle Eocene further paves the way for extending the ATS into the Mesozoic.

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.

Planktic foraminifer-based sea-surface temperature reconstruction for sediment cores JM97-948/2A and MD95-2011 from the Vøring Plateau, eastern Norwegian Sea

The early to mid-Holocene thermal optimum is a well-known feature in a wide variety of paleoclimate archives from the Northern Hemisphere. Reconstructed summer temperature anomalies from across northern Europe show a clear maximum around 6000 years before present (6 ka). For the marine realm, Holocene trends in sea-surface temperature reconstructions for the North Atlantic and Norwegian Sea do not exhibit a consistent pattern of early to mid- Holocene warmth. Sea-surface temperature records based on alkenones and diatoms generally show the existence of a warm early to mid-Holocene optimum. In contrast, several foraminifer and radiolarian based temperature records from the North Atlantic and Norwegian Sea show a cool mid- Holocene anomaly and a trend towards warmer temperatures in the late Holocene. In this paper, we revisit the foraminifer record from the Vøring Plateau in the Norwegian Sea. We also compare this record with published foraminifer based temperature reconstructions from the North Atlantic and with modelled (CCSM3) upper ocean temperatures. Model results indicate that while the seasonal summer warming of the seasurface was stronger during the mid-Holocene, sub-surface depths experienced a cooling. This hydrographic setting can explain the discrepancies between the Holocene trends exhibited by phytoplankton and zooplankton based temperature proxy records.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Seawater carbonate chemistry and calcification during experiments with a coral Madracis auretenra, 2010

Physiological data and models of coral calcification indicate that corals utilize a combination of seawater bicarbonate and (mainly) respiratory CO2 for calcification, not seawater carbonate. However, a number of investigators are attributing observed negative effects of experimental seawater acidification by CO2 or hydrochloric acid additions to a reduction in seawater carbonate ion concentration and thus aragonite saturation state. Thus, there is a discrepancy between the physiological and geochemical views of coral biomineralization. Furthermore, not all calcifying organisms respond negatively to decreased pH or saturation state. Together, these discrepancies suggest that other physiological mechanisms, such as a direct effect of reduced pH on calcium or bicarbonate ion transport and/or variable ability to regulate internal pH, are responsible for the variability in reported experimental effects of acidification on calcification. To distinguish the effects of pH, carbonate concentration and bicarbonate concentration on coral calcification, incubations were performed with the coral Madracis auretenra (= Madracis mirabilis sensu Wells, 1973) in modified seawater chemistries. Carbonate parameters were manipulated to isolate the effects of each parameter more effectively than in previous studies, with a total of six different chemistries. Among treatment differences were highly significant. The corals responded strongly to variation in bicarbonate concentration, but not consistently to carbonate concentration, aragonite saturation state or pH. Corals calcified at normal or elevated rates under low pH (7.6-7.8) when the seawater bicarbonate concentrations were above 1800 µm. Conversely, corals incubated at normal pH had low calcification rates if the bicarbonate concentration was lowered. These results demonstrate that coral responses to ocean acidification are more diverse than currently thought, and question the reliability of using carbonate concentration or aragonite saturation state as the sole predictor of the effects of ocean acidification on coral calcification.

(Appendix A) Abundance of selected coccolithophores in sediment core DS97-2P, North Atlantic

The calcareous nannofossil assemblages from sediment core DS97-2P from the Reykjanes Ridge have been investigated to document oceanographic changes in surface water during the Holocene. The recorded variations in coccolithophore species assemblages and accumulation rates indicate that the region was subjected to rapid changes of surface water masses throughout the entire Holocene. Coccolithophore assemblages generally are of low species diversity and consist mainly of Emiliania huxleyi and Coccolithus pelagicus ssp. pelagicus. Two major events occurred at 8.5-7 ka and at 4.5-3.5 ka, showing higher coccolith accumulation rates, suggesting that the influence of relatively warm Atlantic waters via the Irminger Current was strong in the investigated area. The coccolithophore assemblages have been compared with diatom, foraminifer and sedimentological records within the same core. These data, supported by a comparison with previously published proxy records, add credit to the hypothesis that Holocene changes did not occur uniformly across the North Atlantic. The results have highlighted the Holocene pattern in the North Atlantic, as a period influenced by strong regionalism with discrepancies in the hydrographical trends and in the distribution of the planktonic proxies.

Oxygen consumption in Arctic sediments

Total sediment oxygen consumption rates (TSOC or Jtot), measured during sediment-water incubations, and sediment oxygen microdistributions were studied at 16 stations in the Arctic Ocean (Svalbard area). The oxygen consumption rates ranged between 1.85 and 11.2 mmol m**-2 d**-1, and oxygen penetrated from 5.0 to >59 mm into the investigated sediments. Measured TSOC exceeded the calculated diffusive oxygen fluxes (Jdiff) by 1.1-4.8 times. Diffusive fluxes across the sediment-water interface were calculated using the whole measured microprofiles, rather than the linear oxygen gradient in the top sediment layer. The lack of a significant correlation between found abundances of bioirrigating meiofauna and high Jtot/Jdiff ratios as well as minor discrepancies in measured TSOC between replicate sediment cores, suggest molecular diffusion, not bioirrigation, to be the most important transport mechanism for oxygen across the sediment-water interface and within these sediments. The high ratios of Jtot/Jdiff obtained for some stations were therefore suggested to be caused by topographic factors, i.e. underestimation of the actual sediment surface area when one-dimensional diffusive fluxes were calculated, or sampling artifacts during core recovery from great water depths. Measured TSOC correlated to water depth raised to the -0.4 to -0.5 power (TSOC = water depth**-0.4 to -0.5) for all investigated stations, but they could be divided into two groups representing different geographical areas with different sediment oxygen consumption characteristics. The differences in TSOC between the two areas were suggested to reflect hydrographic factors (such as ice coverage and import/production of reactive particulate organic material) related to the dominating water mass (Atlantic or polar) in each of the two areas. The good correlation between TSOC and water depth**-0.4 to -0.5 rules out any of the stations investigated to be topographic depressions with pronounced enhanced sediment oxygen consumption.

Eemian temperature reconstruction for European sites

Botanical data are widely used as terrestrial proxy data for climate reconstructions. Using a newly established method based on probability density functions (pdf-method), the temperature development throughout the last interglacial, the Eemian, is reconstructed for the two German sites Bispingen and Grobern and the French site La Grande Pile. The results are compared with previous reconstructions using other methods. After a steep increase in January as well as July temperatures in the early phase of the interglacial, the reconstructed most probable climate appears to be slightly warmer than today. While the temperature is reconstructed as relatively stable throughout the Eemian, a certain tendency towards cooler January temperatures is evident. January temperatures decreased from approx. 2-3° C in the early part to approx. -3° C in the later part at Bispingen, and from approx. 2° C to approx. -1° C at Grobern and La Grande Pile. A major drop to about -8° C marks the very end of the interglacial at all three sites. While these results agree well with other proxy data and former reconstructions based on the indicator species method, the results differ significantly from reconstructions based on the modern pollen analogue technique ("pollen transfer functions"). The lack of modern analogues is assumed to be the main reason for the discrepancies. It is concluded that any reconstruction method needs to be evaluated carefully in this respect if used for periods lacking modern analogous plant communities.

Palaeoclimate reconstructions of the Oligocene to Pliocene from eight sites in northwest Germany using the Coexistence Approach

The first detailed reconstruction of the continental palaeoclimate evolution of the Northwest German Tertiary (Late Oligocene to Pliocene) is presented. The paleoclimate data are derived from the paleobotanical record using the coexistence approach, a method recently introduced that employs climatic requirements of the Nearest Living Relatives of a fossil flora. Twenty six megafloras (fruits and seeds, leaves, woods) from the Tertiary succession of the Lower Rhine Basin and neighboring areas are analyzed with respect to ten meteorological parameters. Additionally, two sample sets from Late Miocene to Early Pliocene sediments comprising 396 palynofloras are analyzed by the same method providing a higher temporal resolution. The temperature curves show a comparatively cooler phase in the Late Oligocene, a warm interval the Middle Miocene, and a cooling starting at 14 Ma. The cooling trend persisted until Late Pliocene with a few higher frequency temperature variations observed. From the beginning of Late Miocene to the present, the seasonality increases and climate appears to have been less stable. As indicated by the precipitation data, a Cfa climate with wet summers persisted in NW Germany from Late Oligocene to Late Pliocene. The results obtained are well in accordance with regional and global isotope curves derived from the marine record, and allow for a refined correlation of the Tertiary succession in the Lower Rhine Basin with the international standard. It is shown that the reconstructed data are largely consistent with the continental climate record for the Northern Hemisphere, as reported by various authors. Discrepancies with previous reconstructions are discussed in detail.

Temperature affects the morphology and calcification of Emiliania huxleyi strains

The global warming debate has sparked an unprecedented interest in temperature effects on coccolithophores. The calcification response to temperature changes reported in the literature, however, is ambiguous. The two main sources of this ambiguity are putatively differences in experimental setup and strain-specificity. In this study we therefore compare three strains isolated in the North Pacific under identical experimental conditions. Three strains of Emiliania huxleyi type A were grown under non-limiting nutrient and light conditions, at 10, 15, 20 and 25 ºC. All three strains displayed similar growth rate versus temperature relationships, with an optimum at 20-25 ºC. Elemental production (particulate inorganic carbon (PIC), particulate organic carbon (POC), total particulate nitrogen (TPN)), coccolith mass, coccolith size, and width of the tube elements cycle were positively correlated with temperature over the sub-optimum to optimum temperature range. The correlation between PIC production and coccolith mass/size supports the notion that coccolith mass can be used as a proxy for PIC production in sediment samples. Increasing PIC production was significantly positively correlated with the percentage of incomplete coccoliths in one strain only. Generally, coccoliths were heavier when PIC production was higher. This shows that incompleteness of coccoliths is not due to time shortage at high PIC production. Sub-optimal growth temperatures lead to an increase in the percentage of malformed coccoliths in a strain-specific fashion. Since in total only six strains have been tested thus far, it is presently difficult to say whether sub-optimal temperature is an important factor causing malformations in the field. The most important parameter in biogeochemical terms, the PIC:POC, shows a minimum at optimum growth temperature in all investigated strains. This clarifies the ambiguous picture featuring in the literature, i.e. discrepancies between PIC:POC-temperature relationships reported in different studies using different strains and different experimental setups. In summary, global warming might cause a decline in coccolithophore's PIC contribution to the rain ratio, as well as improved fitness in some genotypes due to less coccolith malformations.