Individual foraminifera variability from the western and eastern equatorial Pacific during the Late Holocene and Last Glacial Maximum

El Niño-Southern Oscillation (ENSO) is a major source of global interannual variability, but its response to climate change is uncertain. Paleoclimate records from the Last Glacial Maximum (LGM) provide insight into ENSO behavior when global boundary conditions (ice sheet extent, atmospheric partial pressure of CO2) were different from those today. In this work, we reconstruct LGM temperature variability at equatorial Pacific sites using measurements of individual planktonic foraminifera shells. A deep equatorial thermocline altered the dynamics in the eastern equatorial cold tongue, resulting in reduced ENSO variability during the LGM compared to the Late Holocene. These results suggest that ENSO was not tied directly to the east-west temperature gradient, as previously suggested. Rather, the thermocline of the eastern equatorial Pacific played a decisive role in the ENSO response to LGM climate.

Elemental, biomarker, palynology and maceral data from Seam 1, Schöningen and associated interbeds

Sphagnum moss is the dominant plant type in modern boreal and (sub)arctic ombrotrophic bogs and is of particular interest due to its sensitivity to climate and its important role in wetland biogeochemistry. Here we reconstruct the occurrence of Sphagnum moss - and associated biogeochemical change - within a thermally immature, early Paleogene (~55 Ma) lignite from Schöningen, NW Germany using a high-resolution, multi-proxy approach. Changes in the abundance of Sphagnum-type spores and the C23/C31n-alkane ratio indicate the expansion of Sphagnum moss within the top of the lignite seam. This Sphagnum moss expansion is associated with the development of waterlogged conditions, analogous to what has been observed within modern ombrotrophic bogs. The similarity between biomarkers and palynology also indicates that the C23/C31n-alkane ratio may be a reliable chemotaxonomic indicator for Sphagnum during the early Paleogene. The d13C value of bacterial hopanes and mid-chain n-alkanes indicates that a rise in water table is not associated with a substantial increase in aerobic methanotrophy. The absence of very low d13C values within the top of the seam could reflect either less methanogenesis or less efficient methane oxidation under waterlogged sulphate-rich conditions.

Radiogenic isotope record of Arctic Ocean circulation and weathering inputs of the past 15 million years

Lead (Pb), neodymium (Nd), and strontium (Sr) isotopic analyses were carried out on sediment leachates (reflecting the isotope composition of past seawater) and digests of the bulk residues (reflecting detrital continental inputs) of Integrated Ocean Drilling Program (IODP) Leg 302 and core PS2185 from the Lomonosov Ridge (Arctic Ocean). Our records are interpreted to reflect changes in continental erosion and oceanic circulation, driven predominantly by tectonic forcing on million-year timescales in the older (pre-2 Ma) part of the record and by climatic forcing of weathering and erosion of the Eurasian continental margin on thousand-year timescales in the younger (post-2 Ma) part. These data, covering the past ~15 Ma, show that continental inputs to the central Arctic Ocean have been more closely linked to glacial and hydrological processes occurring on the Eurasian margin than on continental North America and Greenland. The constancy of the detrital input signatures supports the early existence of an Arctic sea ice cover, whereas the major initiation of Northern Hemisphere glaciation at 2.7 Ma appears to have had little impact on the weathering regime of the Eurasian continental margin.

Stable oxygen isotope and Mg/Ca ratios on Globorotalia inflata

We present a species-specific Mg/Ca-calcification temperature calibration for Globorotalia inflata from a suite of 38 core top samples from the South Atlantic (from 8° to 49°S). G. inflata is a deep-dwelling planktonic foraminifer commonly occurring in subtropical to subpolar conditions, which qualifies it for reconstructions of the permanent thermocline. Apparent calcification depths and calcification temperatures were determined by comparing measured d18O with equilibrium d18O of calcite based on water column properties. Based on our core top samples, G. inflata apparent calcification depth is constant throughout the South Atlantic mid-latitudes with a depth of 350-400 m within the permanent thermocline. The resulting Mg/Ca-calcification temperature calibration is Mg/Ca = 0.72 +/-0.045/0.042 exp (0.076 +0.006 calcification 2 temperature) (r2 = 0.81) and covers the temperature range 3.1-16.5°C. We applied our Mg/Ca calibration to gravity core PS2495-3 from the Mid-Atlantic Ridge at ca. 41°S to test its validity by reconstructing a low-resolution record covering the last two glacial-interglacial cycles. Our paleotemperature record reveals large changes in temperature for Terminations I and II, when permanent thermocline temperature increased by as much as 8°C. The G. inflata paleotemperature record suggests that oceanic fronts repeatedly migrated over the location of site PS2495-3 during the last 160 kyr. This study shows the potential of G. inflata Mg/Ca to reconstruct paleotemperatures in the permanent thermocline.

Organic-geochemical proxies (UK'37 and TEXH86) records in surface sediments of the tropical Eastern Indian Ocean

In this study we reconstruct sea surface temperatures (SSTs) using two lipid-based biomarker proxies (alkenone unsaturation index UK'37 and TEX86 index based on glycerol dibiphytanyl glycerol tetraethers) in 36 surface sediment samples from the Indonesian continental margin off west Sumatra and south of Java and the Lesser Sunda Islands. Comparison of measured temperatures (World Ocean Atlas 09) to reconstructed temperatures suggests that SST-UK'37 reflects the SE monsoon SST in the upwelling area south of Java and the Lesser Sunda Islands, whereas Temp-TEXH86 estimates are up to 2°C lower than SST-UK'37. This offset is possibly related to either one or a combination of two factors: i) the depth habitats of the source organisms; ii) different seasonal production and/ or seasonality of export associated with phytoplankton blooming triggered by primary productivity. In the non-upwelling area off west Sumatra, the alkenone-based SSTs are cooler than measured temperatures during the entire year, likely due to the reduced sensitivity of the UK'37 proxy beyond 28°C. However, reconstructed temperatures based on TEXH86 are consistent with mean annual SST, implying that the Temp-TEXH86 reflects the mean annual SST in the non-upwelling area of the tropical Eastern Indian Ocean.

(Table 1) Lipid biomarkers in surface sediments from the Gulf of Genoa

A series of molecular organic markers were determined in surface sediments from the Gulf of Genoa (Ligurian Sea) in order to evaluate their potential for palaeo-environmental reconstructions. Allochthonous input can be characterized by the distributions of n-C29 and n-C31 alkanes, n-C26 and n-C28 alkanols and branched glycerol dialkyl glycerol tetraethers (GDGTs), whose concentrations are generally highest near the river mouths. In the open basin however, terrestrial n-alkanes and n-alkanols may have an additional, aeolian source. Autochthonous input is represented by crenarchaeol and isoprenoid GDGTs. Their concentrations are highest in the open basin showing the preference of Thaumarchaeota for oligotrophic waters. Indications of a significant degradation of sterols and C37 alkenones exclude these lipids as reliable productivity proxies. Using terrestrial and aquatic lipids as end-members allows estimating the percentage of terrestrial organic matter between 20% and 58% in the coastal area decreasing to 1 to 30% in the deep basin. The spatial distribution of sea surface temperature (SST) estimates using the alkenone-based UK'37 index is very similar to the autumnal (November) mean satellite-based SST distribution. Conversely, TEXH86-derived SST estimates are close to winter SSTs in the coastal area and summer SSTs in the open basin. This pattern reflects presumably a shift in the main production of Thaumarchaeota from the coastal area in winter to the open basin in summer. This study represents a major prerequisite for the future application of lipid biomarkers on sediment cores from the Gulf of Genoa.

Oligocene to early Miocene benthic foraminifera in the eastern Equatorial Pacific

We investigated Oligocene and early Miocene benthic foraminiferal faunas (> 105 µm in size) from Ocean Drilling Program (Leg 199) Site 1218 (4826 m water depth and ~3300 to ~4000 m paleo-water depth) and Site 1219 (5063 m water depth and ~4200 to ~4400 m paleo-water depth) to understand the response of abyssal benthic foraminifera to mid-Oligocene glacial events in the eastern Equatorial Pacific Ocean. Two principal factor assemblages were recognized. The Factor 1 assemblage (common Nuttallides umbonifer) is related to either an influx of the Southern Component Water (SCW), possibly carbonate undersaturated, or a decrease in seasonality of the food supply from the surface ocean. The Factor 2 assemblage is characterized by typical deep-sea taxa living under variable trophic conditions, possibly with a seasonal component in food supply. The occurrence of abyssal benthic foraminifera faunas during the mid-Oligocene depends on either the effect of SCW or the seasonality of food resources. The Factor 1 assemblage was most common near 76Ol-C11r, 73Ol-C10rn and 67Ol-C9n (ca. 30.2, 29.1 and 26.8 Ma respectively by Pälike et al. (2006, doi:10.1126/science.1133822)). This indicates that the effect of SCW increased or the seasonal input of food from the surface ocean to benthic environments was weakened close to these glacial events. In contrast, the huge export flux of small biogenic carbonate particles close to these glacial events might be responsible for carbonate-rich sediments buffering carbonate undersaturation. Changes in deep-water masses or the periodicity of food supply from the surface ocean and variation in surface carbonate production affected by orbital forcing had an impact on the mid-Oligocene faunas of abyssal benthic foraminifera around the intervals of glacial events in the eastern Equatorial Pacific Ocean. The Factor 1 assemblage decreased sharply at ? 30 Ma (29.8 Ma by Pälike et al. (2006), 30.0 Ma by CK95) and returned to dominance after ? 29 Ma (28.6 Ma by Pälike et al. (2006), 28.8 Ma by CK95). It is likely that the effect of SCW (possibly carbonate undersaturated) has intensified since the late Oligocene. The faunal transition of benthic foraminifera in the eastern Equatorial Pacific Ocean at ~29 Ma might be attributable to the influence of Northern Component Water (NCW) input to the Southern Ocean and the subsequent formation of SCW at about that time.

Sea surface temperatures reconstructed from sediment core U1338

This study presents a new alkenone-derived Sea Surface Temperature (SST) record and d18ONoelaerhabdaceae data of the 2-5 µm carbonate fractions from the IODP site U1338 located in the Eastern Equatorial Pacific (EEP), over the Miocene-Pliocene. Our data and those available from other sites of the same area show the establishment of a cold tongue during the early Pliocene (4.4 - 3.6 Ma). SST and d18ONoelaerhabdaceae time-series indicate periods of significant salinity variations. Comparison with the d18Obenthic curve from sediment cores of the Equatorial Pacific Ocean allow us to distinguish between global changes and local surface salinity variations in the EEP. Ice sheet growth and evaporation-precipitation are then discussed as possible drivers of such changes, as well as the role of Central American and the Indonesian seaway restriction. Our data suggest a shallowing of the thermocline in the EEP, between 6.8 and 6 Ma, and its shoaling between 4.8 and 4.0 Ma, suggesting the appearance of the cold tongue (Steph et al., 2010). The Pliocene climate transition would therefore not be primarily driven by Northern hemisphere glaciation.

Chemical concentrations and fluxes from EPICA ice cores EDML and EDC

Continuous sea salt and mineral dust aerosol records have been studied on the two EPICA (European Project for Ice Coring in Antarctica) deep ice cores. The joint use of these records from opposite sides of the East Antarctic plateau allows for an estimate of changes in dust transport and emission intensity as well as for the identification of regional differences in the sea salt aerosol source. The mineral dust flux records at both sites show a strong coherency over the last 150 kyr related to dust emission changes in the glacial Patagonian dust source with three times higher dust fluxes in the Atlantic compared to the Indian Ocean sector of the Southern Ocean (SO). Using a simple conceptual transport model this indicates that transport can explain only 40% of the atmospheric dust concentration changes in Antarctica, while factor 5-10 changes occurred. Accordingly, the main cause for the strong glacial dust flux changes in Antarctica must lie in environmental changes in Patagonia. Dust emissions, hence environmental conditions in Patagonia, were very similar during the last two glacials and interglacials, respectively, despite 2-4 °C warmer temperatures recorded in Antarctica during the penultimate interglacial than today. 2-3 times higher sea salt fluxes found in both ice cores in the glacial compared to the Holocene are difficult to reconcile with a largely unchanged transport intensity and the distant open ocean source. The substantial glacial enhancements in sea salt aerosol fluxes can be readily explained assuming sea ice formation as the main sea salt aerosol source with a significantly larger expansion of (summer) sea ice in the Weddell Sea than in the Indian Ocean sector. During the penultimate interglacial, our sea salt records point to a 50% reduction of winter sea ice coverage compared to the Holocene both in the Indian and Atlantic Ocean sector of the SO. However, from 20 to 80 ka before present sea salt fluxes show only very subdued millennial changes despite pronounced temperature fluctuations, likely due to the large distance of the sea ice salt source to our drill sites.

Mid-Pleistocene (800-400 ka) records of benthic and planktic foraminfer d18O and d13C, planktic foraminifer based SST, and IRD of IODP Site 306-U1314

Stable carbon and oxygen isotopes from benthic and planktic foraminifers, planktic foraminifer assemblages and ice rafted debris from the North Atlantic Site U1314 (Integrated Ocean Drilling Program Expedition 306) were examined to investigate orbital and millennial-scale climate variability in the North Atlantic and its impact on global circulation focusing on the development of glacial periods during the mid-Pleistocene (ca 800-400 ka). Glacial initiations were characterized by a rapid cooling (6-10 °C in less than 7 kyr) in the mean annual sea surface temperature (SST), increasing benthic d18O values and high benthic d13C values. The continuous increase in benthic d18O suggests a continuous ice sheet growth whereas the positive benthic d13C values indicate that the flow of the Iceland Scotland Overflow water (ISOW) was vigorous. Strong deep water formation in the Norwegian Greenland Sea promoted a high transfer of freshwater from the ocean to the continents. However, low SSTs at Site U1314 suggest a subpolar gyre cooling and freshening that may have reduced deep water formation in the Labrador Sea during glacial initiations. Once the 3.5 per mil threshold in the benthic d18O record was exceeded, ice rafting started and ice sheet growth was punctuated by millennial-scale waning events which returned to the ocean part of the freshwater accumulated on the continents. Ice-rafting events were associated with a rapid reduction in the ISOW (benthic d13C values dropped 0.5-1 per mil) and followed by millennial-scale warmings. The first two millennial-scale warm intervals of each glacial period reached interglacial temperatures and were particularly abrupt (6-10 °C in ~3 kyr). Subsequent millennial-scale warm events were cooler probably because the AMOC was rather reduced as suggested by the low benthic d13C values. These two abrupt warming events that occurred at early glacial periods were also observed in the Antarctic temperature and CO2 records, suggesting a close correlation between both Hemispheres. The comparison of the sea surface proxies with the benthic d18O record (as the Southern sign) indicates the presence of a millennial-scale seesaw pattern similar to that seen during the Last Glacial period.

The Antarctic ice core chronology (AICC2012)

Abstract of Bazin et al. (2013): An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120-800 ka. In this framework, new measurements of d18Oatm over Marine Isotope Stage (MIS) 11-12 on EDC and a complete d18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from d18Oatm, dO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one. Abstract of Veres et al. (2013): The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on d15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores.

Model - data comparison of western African rainfall evolutions during the Last Interglacial (130-115 ka)

This dataset characterizes the evolution of western African precipitation indicated by marine sediment geochemical records in comparison to transient simulations using CCSM3 global climate model throughout the Last Interglacial (130-115 ka). It contains (1) defined tie-points (age models), newly published stable isotopes of benthic foraminifera and Al/Si log-ratios of eight marine sediment cores from the western African margin and (2) annual and seasonal rainfall anomalies (relative to pre-industrial values) for six characteristic latitudinal bands in western Africa simulated by CCSM3 (two transient simulations: one non-accelerated and one accelerated experiment).