Pliocene-Pleistocene stack of globally distributed benthic stable oxygen isotope records

We present a 5.3-Myr stack (the ''LR04'' stack) of benthic d18O records from 57 globally distributed sites aligned by an automated graphic correlation algorithm. This is the first benthic delta18O stack composed of more than three records to extend beyond 850 ka, and we use its improved signal quality to identify 24 new marine isotope stages in the early Pliocene. We also present a new LR04 age model for the Pliocene-Pleistocene derived from tuning the delta18O stack to a simple ice model based on 21 June insolation at 65 N. Stacked sedimentation rates provide additional age model constraints to prevent overtuning. Despite a conservative tuning strategy, the LR04 benthic stack exhibits significant coherency with insolation in the obliquity band throughout the entire 5.3 Myr and in the precession band for more than half of the record. The LR04 stack contains significantly more variance in benthic delta18O than previously published stacks of the late Pleistocene as the result of higher resolution records, a better alignment technique, and a greater percentage of records from the Atlantic. Finally, the relative phases of the stack's 41- and 23-kyr components suggest that the precession component of delta18O from 2.7-1.6 Ma is primarily a deep-water temperature signal and that the phase of d18O precession response changed suddenly at 1.6 Ma.

Sea surface temperature reconstruction based on alkenones of sediment core M35003-4

Evidence for abrupt climate changes on millennial and shorter timescales is widespread in marine and terrestrial climate records (Dansgard et al., 1993, doi:10.1038/364218a0; Bond et al., 1993, doi:10.1038/365143a0; Charles et al., 1996, doi:10.1016/0012-821X(96)00083-0, Bard et al., 1997, doi:10.1038/385707a0). Rapid reorganization of ocean circulation is considered to exert some control over these changes (Broecker et al., 1985, doi:10.1038/315021a0), as are shifts in the concentrations of atmospheric greenhouse gases (Broecker, 1994, doi:10.1038/372421a0). The response of the climate system to these two influences is fundamentally different: slowing of thermohaline overturn in the North Atlantic Ocean is expected to decrease northward heat transport by the ocean and to induce warming of the tropical Atlantic (Crowley, 1992, doi:10.1029/92PA01058; Manabe and Stouffer, 1997, doi:10.1029/96PA03932), whereas atmospheric greenhouse forcing should cause roughly synchronous global temperature changes (Manabe et al., 1991, doi:10.1175/1520-0442(1991)004<0785:TROACO>2.0.CO;2). So these two mechanisms of climate change should be distinguishable by the timing of surface-water temperature variations relative to changes in deep-water circulation. Here we present a high-temporal-resolution record of sea surface temperatures from the western tropical North Atlantic Ocean which spans the past 29,000 years, derived from measurements of temperature-sensitive alkenone unsaturation in sedimentary organic matter. We find significant warming is documented for Heinrich event H1 (16,900-15,400 calendar years bp) and the Younger Dryas event (12,900-11,600 cal. yr bp), which were periods of intense cooling in the northern North Atlantic. Temperature changes in the tropical and high-latitude North Atlantic are out of phase, suggesting that the thermohaline circulation was the important trigger for these rapid climate changes.

Figure 2. Downcore variations of total diatom concentration in core MD02-2529

The modern eastern equatorial Pacific (EEP) is a major natural source for atmospheric carbon dioxide and is thought to be connected to high-latitude ocean dynamics by oceanic teleconnections on glacial-interglacial timescales. A wealth of sedimentary records aiming at reconstructing last Quaternary changes in primary productivity and nutrient utilization have been devoted to understanding those linkages between the EEP and other distant oceanic areas. Most of these records are, however, clustered in the pelagic EEP cold tongue, with comparatively little attention devoted to coastal areas. Here we present downcore measurements of the composition and concentration of the diatom assemblage together with opal (biogenic silica) concentration at site MD02-2529 recovered in the coastal Panama Basin. Piston core MD02-2529, collected in an area affected by a multitude of processes, provides evidence for strong variations in diatom production at the millennial timescale during the last glacial cycle. The maxima in total diatom concentration occurred during the early marine isotopic stage (MIS) 4 as well as during the MIS 4/3 transition and MIS 3. Rapid changes in diatom concentrations during the MIS 3 mimics Bond cycles as independently recorded by the SSS estimation derived from planktonic foraminifera from the same core. Such patterns indicate a clear linkage between diatom production in the coastal EEP and rapid climate changes in the high-latitude North Atlantic. In parallel, the long-term succession of the diatom community from coastal diatoms, predominantly thriving during MIS 5 and 4, towards pelagic diatoms, dominant during MIS 3 and 2, points to a long-term change in the surface hydrology. During Heinrich Events, diatoms strongly reduced their production, probably due to enhanced stratification in the upper water column. After the last glacial maximum, diatom production and valve preservation strongly decreased in response to the advection of nutrient (H2SiO4)-depleted, warmer water masses. Our high-resolution record highlights how regional climatic processes can modulate rapid changes in siliceous primary production as triggered by wind-induced local upwelling, indicating that millennial climatic variability can overtake other prominent hydrological processes such as those related to silicic acid leakage.

Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Caribbean and North Atlantic

Variations in the strength of the North Atlantic Ocean thermohaline circulation have been linked to rapid climate changes during the last glacial cycle through oscillations in North Atlantic Deep Water formation and northward oceanic heat flux. The strength of the thermohaline circulation depends on the supply of warm, salty water to the North Atlantic, which, after losing heat to the atmosphere, produces the dense water masses that sink to great depths and circulate back south. Here we analyse two Caribbean Sea sediment cores, combining Mg/Ca palaeothermometry with measurements of oxygen isotopes in foraminiferal calcite in order to reconstruct tropical Atlantic surface salinity during the last glacial cycle. We find that Caribbean salinity oscillated between saltier conditions during the cold oxygen isotope stages 2, 4 and 6, and lower salinities during the warm stages 3 and 5, covarying with the strength of North Atlantic Deep Water formation. At the initiation of the Bølling/Allerød warm interval, Caribbean surface salinity decreased abruptly, suggesting that the advection of salty tropical waters into the North Atlantic amplified thermohaline circulation and contributed to high-latitude warming.

(Table 1) Stable oxygen isotope ratios of Globigerinoides sacculifer of ODP Hole 130-806B

A detailed Pliocene oxygen isotope record from the Ontong Java Plateau, based on measurements of the surface-dwelling planktonic foraminifer Globigerinoides sacculifer, was produced for the period from 5 to 2 Ma. The record documents major long- and short-term climate changes. The results show periods of enhanced ice volume at 4.6 to 4.3 Ma and after 2.85 Ma, a long-term warming trend from 4.1 to 3.7 Ma, and a distinct cooling trend that was initiated at 3.5 Ma and progressed through the initiation of large-scale Northern Hemisphere glaciation after 2.85 Ma (according to the time scale of Shackleton and others proposed in 1990). Periods of high average ice volumes also show the highest d18O amplitudes. The pattern of climate cyclicity changed markedly at about 2.85 Ma. Earlier times were marked by high-frequency variability at the precessional frequencies or even higher frequencies, pointing to low-latitude processes as a main controlling factor driving planktonic d18O variability in this period. The high-frequency variability is not coherent with insolation and points to strong nonlinearity in the way the climate system responded to orbital forcing before the onset of large scale Northern Hemisphere glaciation. After 3 Ma, stronger 41-k.y. cyclicity appears in the record. The shift in pattern is clearest around 2.85 Ma (according to the time scale proposed by Shackleton and others in 1990), 100-200 k.y. before the most dramatic spread of Northern Hemisphere ice sheets. This indicates that high-latitude processes from this point on began to take over and influence most strongly the d18O record, which now reflects ice-volume fluctuations related to the climatic effects of obliquity forcing on the seasonality of high-latitude areas, most probably in the Northern Hemisphere. The general Pliocene trend is that high-latitude climate sensitivity and instability was increasing, and the causal factors producing the intensified glacial cyclicity during the Pliocene must be factors that enhance cooling and climate sensitivity in the subarctic areas.

Age determination and stable isotope ratios of planktonic foraminifera from DSDP Hole 94-609

As shown by the work of Dansgaard and his colleagues, climate oscillations of one or so millennia duration punctuate much of glacial section of the Greenland ice cores. These oscillations are characterized by 5°C air temperature changes, severalfold dust content changes and 50 ppm CO2 changes. Both the temperature and CO2 change are best explained by changes in the mode of operation of the ocean. In this paper we provide evidence which suggests that oscillations in surface water conditions of similar duration are present in the record from a deep sea core at 50°N. Based on this finding, we suggest that the Greenland climate changes are driven by oscillations in the salinity of the Atlantic Ocean which modulate the strength of the Atlantic's conveyor circulation.

(Appendix 1) Coral analysis and isostatic rebound effects from different Holes of IODP Expedition 310

Past sea-level records provide invaluable information about the response of ice sheets to climate forcing. Some such records suggest that the last deglaciation was punctuated by a dramatic period of sea-level rise, of about 20 metres, in less than 500 years. Controversy about the amplitude and timing of this meltwater pulse (MWP-1A) has, however, led to uncertainty about the source of the melt water and its temporal and causal relationships with the abrupt climate changes of the deglaciation. Here we show that MWP-1A started no earlier than 14,650 years ago and ended before 14,310 years ago, making it coeval with the Bølling warming. Our results, based on corals drilled offshore from Tahiti during Integrated Ocean Drilling Project Expedition 310, reveal that the increase in sea level at Tahiti was between 12 and 22 metres, with a most probable value between 14 and 18 metres, establishing a significant meltwater contribution from the Southern Hemisphere. This implies that the rate of eustatic sea-level rise exceeded 40 millimetres per year during MWP-1A.

Temperature measured in permafrost borehole at ICDP site 5011-3 in El'gygytgyn Crater, 2008-2011

This study focuses on the temperature field observed in boreholes drilled as part of interdisciplinary scientific campaign targeting the El'gygytgyn Crater Lake in NE Russia. Temperature data are available from two sites: the lake borehole 5011-1 located near the center of the lake reaching 400 m depth, and the land borehole 5011-3 at the rim of the lake, with a depth of 140 m. Constraints on permafrost depth and past climate changes are derived from numerical simulation of the thermal regime associated with the lake-related talik structure. The thermal properties of the subsurface needed for these simulations are based on laboratory measurements of representative cores from the quaternary sediments and the underlying impact-affected rock, complemented by further information from geophysical logs and data from published literature. The temperature observations in the lake borehole 5011-1 are dominated by thermal perturbations related to the drilling process, and thus only give reliable values for the lowermost value in the borehole. Undisturbed temperature data recorded over more than two years are available in the 140 m deep land-based borehole 5011-3. The analysis of these observations allows determination of not only the recent mean annual ground surface temperature, but also the ground surface temperature history, though with large uncertainties. Although the depth of this borehole is by far too insufficient for a complete reconstruction of past temperatures back to the Last Glacial Maximum, it still affects the thermal regime, and thus permafrost depth. This effect is constrained by numerical modeling: assuming that the lake borehole observations are hardly influenced by the past changes in surface air temperature, an estimate of steady-state conditions is possible, leading to a meaningful value of 14 ± 5 K for the post-glacial warming. The strong curvature of the temperature data in shallower depths around 60 m can be explained by a comparatively large amplitude of the Little Ice Age (up to 4 K), with low temperatures prevailing far into the 20th century. Other mechanisms, like varying porosity, may also have an influence on the temperature profile, however, our modeling studies imply a major contribution from recent climate changes.

Sea surface temperature estimation for the Iberian Margin

Centennial climate variability over the last ice age exhibits clear bipolar behavior. High-resolution analyses of marine sediment cores from the Iberian margin trace a number of associated changes simultaneously. Proxies of sea surface temperature and water mass distribution, as well as relative biomarker content, demonstrate that this typical north-south coupling was pervasive for the cold phases of climate during the past 420,000 years. Cold episodes after relatively warm and largely ice-free periods occurred when the predominance of deep water formation changed from northern to southern sources. These results reinforce the connection between rapid climate changes at Mediterranean latitudes and century-to-millennial variability in northern and southern polar regions.

(Table 1) Organic carbon and radiogenic isotopes analysis of DSDP Hole 22-218 sediments

Modern erosion of the Himalaya, the world's largest mountain range, transfers huge dissolved and particulate loads to the ocean. It plays an important role in the long-term global carbon cycle, mostly through enhanced organic carbon burial in the Bengal Fan. To understand the role of past Himalayan erosion, the influence of changing climate and tectonic on erosion must be determined. Here we use a 12 Myr sedimentary record from the distal Bengal Fan (Deep Sea Drilling Project Site 218) to reconstruct the Mio-Pliocene history of Himalayan erosion. We use carbon stable isotopes (d13C) of bulk organic matter as paleo-environmental proxy and stratigraphic tool. Multi-isotopic - Sr, Nd and Os - data are used as proxies for the source of the sediments deposited in the Bengal Fan over time. d13C values of bulk organic matter shift dramatically towards less depleted values, revealing the widespread Late Miocene (ca. 7.4 Ma) expansion of C4 plants in the basin. Sr, Nd and Os isotopic compositions indicate a rather stable erosion pattern in the Himalaya range during the past 12 Myr. This supports the existence of a strong connection between the southern Tibetan plateau and the Bengal Fan. The tectonic evolution of the Himalaya range and Southern Tibet seems to have been unable to produce large re-organisation of the drainage system. Moreover, our data do not suggest a rapid change of the altitude of the southern Tibetan plateau during the past 12 Myr. Variations in Sr and Nd isotopic compositions around the late Miocene expansion of C4 plants are suggestive of a relative increase in the erosion of High Himalaya Crystalline rock (i.e. a simultaneous reduction of both Transhimalayan batholiths and Lesser Himalaya relative contributions). This could be related to an increase in aridity as suggested by the ecological and sedimentological changes at that time. A reversed trend in Sr and Nd isotopic compositions is observed at the Plio-Pleistocene transition that is likely related to higher precipitation and the development of glaciers in the Himalaya. These almost synchronous moderate changes in erosion pattern and climate changes during the late Miocene and at the Plio-Pleistocene transition support the notion of a dominant control of climate on Himalayan erosion during this time period. However, stable erosion regime during the Pleistocene is suggestive of a limited influence of the glacier development on Himalayan erosion.