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.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

(Table S1) Age determination of ODP Site 202-1240

Interannual-decadal variability in the equatorial Pacific El Niño-Southern Oscillation (ENSO) induces climate changes at global scale, but its potential influence during past global climate change is not yet well constrained. New high-resolution eastern equatorial Pacific proxy records of thermocline conditions present new evidence of strong orbital control in ENSO-like variability over the last 275,000 years. Recurrent intervals of saltier thermocline waters are associated with the dominance of La Niña-like conditions during glacial terminations, coinciding with periods of low precession and high obliquity. The parallel dominance of d13C-depleted waters supports the advection of Antarctic origin waters toward the tropical thermocline. This "oceanic tunneling" is proposed to have reinforced orbitally induced changes in ENSO-like variability, composing a complex high- and low-latitude feedback during glacial terminations.

(Table 2) Age determination of sediment cores from the Panama Basin

Application of the 230Th normalization method to estimate sediment burial fluxes in six cores from the eastern equatorial Pacific (EEP) reveals that bulk sediment and organic carbon fluxes display a coherent regional pattern during the Holocene that is consistent with modern oceanographic conditions, in contrast with estimates of bulk mass accumulation rates (MARs) derived from core chronologies. Two nearby sites (less than 10 km apart), which have different MARs, show nearly identical 230Th-normalized bulk fluxes. Focusing factors derived from the 230Th data at the foot of the Carnegie Ridge in the Panama Basin are >2 in the Holocene, implying that lateral sediment addition is significant in this part of the basin. New geochemical data and existing literature provide evidence for a hydrothermal source of sediment in the southern part of the Panama Basin and for downslope transport from the top of the Carnegie Ridge. The compilation of core records suggests that sediment focusing is spatially and temporally variable in the EEP. During oxygen isotope stage 2 (OIS 2, from 13-27 ka BP), focusing appears even higher compared to the Holocene at most sites, similar to earlier findings in the eastern and central equatorial Pacific. The magnitude of the glacial increase in focusing factors, however, is strongly dependent on the accuracy of age models. We offer two possible explanations for the increase in glacial focusing compared to the Holocene. The first one is that the apparent increase in lateral sediment redistribution is partly or even largely an artifact of insufficient age control in the EEP, while the second explanation, which assumes that the observed increase is real, involves enhanced deep sea tidal current flow during periods of low sea level stand.

Stable isotope record of Cibicidoides wuellerstorfi from Pleistocene sediments of DSDP Site 86-502 in the Caribbean Sea

Variations in the stable isotopic composition of benthic foraminifera from Deep-Sea Drilling Project (DSDP) site 502B in the Caribbean Sea are used to reconstruct Atlantic intermediate water circulation variability over the last 1.2 m.y. Comparison of this record with other North Atlantic benthic isotope records indicates that Atlantic intermediate water circulation was relatively enhanced during glacial maxima when North Atlantic deep water (NADW) production was reduced. However, a simple, compensatory relationship between intermediate and deepwater circulation is not apparent. Geochemical models have shown that such changes in ocean circulation can affect atmospheric CO2 levels by changing vertical nutrient and alkalinity profiles. The Delta delta13C difference between Caribbean site 502B and deep equatorial Pacific site 677 is highly coherent and in phase with ice volume. Like the delta18O record, there is an increase in amplitude (40%) and a large increase in 100 kyr power after 0.7 Ma. The 1.2? Delta delta13C amplitude scales to 70 ppm V in atmospheric CO2 using Boyle's (1986) box model result. The implied increase in CO2 amplitude after 0.7 Ma may suggest a positive feedback role in effecting the higher-amplitude climatic fluctuations which characterize the last 0.7 m.y.