Results of Heinrich event 1 simulation from the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model (CCSM3)

We investigated changes in tropical climate and vegetation cover associated with abrupt climate change during Heinrich Event 1 (HE1, ca. 17.5 ka BP) using two different global climate models: the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model version 3 (CCSM3). Tropical South American and African pollen records suggest that the cooling of the North Atlantic Ocean during HE1 influenced the tropics through a southward shift of the rain belt. In this study, we simulated the HE1 by applying a freshwater perturbation to the North Atlantic Ocean. The resulting slowdown of the Atlantic Meridional Overturning Circulation was followed by a temperature seesaw between the Northern and Southern Hemispheres, as well as a southward shift of the tropical rain belt. The shift and the response pattern of the tropical vegetation around the Atlantic Ocean were more pronounced in the CCSM3 than in the UVic ESCM simulation. For tropical South America, opposite changes in tree and grass cover were modeled around 10° S in the CCSM3 but not in the UVic ESCM. In tropical Africa, the grass cover increased and the tree cover decreased around 15° N in the UVic ESCM and around 10° N in the CCSM3. In the CCSM3 model, the tree and grass cover in tropical Southeast Asia responded to the abrupt climate change during the HE1, which could not be found in the UVic ESCM. The biome distributions derived from both models corroborate findings from pollen records in southwestern and equatorial western Africa as well as northeastern Brazil.

Short-term variations in radiolarian faunal dsitribution in the eastern equatorial Pacific

Two short time intervals centered at 2.3 and 4.7 Ma were studied to investigate short-term variations in surface-ocean processes as indicated by changes in the radiolarian microfossil population. These time intervals represent two different settings of late Neogene climate. The older interval represents a time when tropical circulation between the Pacific and Atlantic oceans was not blocked by the Isthmus of Panama, whereas the younger interval represents a time when Northern Hemisphere glaciation was present but did not display the dominance of the 100,000-yr cycle that characterizes the late Pleistocene. The younger time slice at 2.3 Ma was sampled at all Leg 138 sites except Site 844, where significant reworking was evident. All sites except 844, 853, and 854 were sampled for the older time slice. Samples were taken at 10- to 20-cm intervals at each site and spanned a GRAPE density maximum and minimum. Thus, it was possible to investigate whether the changes in carbonate content (as indicated by GRAPE density) were associated with changes in surface-ocean conditions (indicated by radiolarian assemblage variations). For both time slices, the radiolarian data indicate that intervals of decreased carbonate content are periods of cooler water conditions and possibly enhanced biogenic production. Times of increased carbonate content are associated with inferred warmer oceanographic conditions, as indicated by the dominance of tropical assemblages at 2.3 Ma and tropical and western Pacific assemblages during the time slice centered at 4.8 Ma. However, the spatial patterns of change during each time slice show a distinct difference in the mapped patterns of radiolarian assemblage dominance. The older time slice, representing a period before the closing of the Isthmus of Panama, shows more zonal patterns presumably associated with a more zonal character of equatorial circulation. After the closing of the isthmus, the shifts in faunal patterns between times of high and low carbonates are characterized by shifts in the dominance of the tropical and transitional assemblages, respectively, throughout the region.

(Table 1) Sea surface temperature reconstruction for eastern equatorial Pacific surface sediment samples

Significant uncertainties persist in the reconstruction of past sea surface temperatures in the eastern equatorial Pacific, especially regarding the amplitude of the glacial cooling and the details of the post-glacial warming. Here we present the first regional calibration of alkenone unsaturation in surface sediments versus mean annual sea surface temperatures (maSST). Based on 81 new and 48 previously published data points, it is shown that open ocean samples conform to established global regressions of Uk'37 versus maSST and that there is no systematic bias from seasonality in the production or export of alkenones, or from surface ocean nutrient concentrations or salinity. The flattening of the regression at the highest maSSTs is found to be statistically insignificant. For the near-coastal Peru upwelling zone between 11-15°S and 76-79°W, however, we corroborate earlier observations that Uk'37 SST estimates significantly over-estimate maSSTs at many sites. We posit that this is caused either by uncertainties in the determination of maSSTs in this highly dynamic environment, or by biasing of the alkenone paleothermometer toward El Niño events as postulated by Rein et al. (2005).

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.

Stable carbon and oxygen isotope ratios of planktonic foraminifera from the equatorial Atlantic

Carbon isotopic records of nutrient-depleted surface water place constraints on the past fertility of the oceans and on past atmospheric pCO2 levels. The best records of nutrient-depleted delta13C are obtained from planktonic foraminifera living in the thick mixed layers of the western equatorial and tropical Atlantic Ocean. We have produced a composite, stacked Globigerinoides sacculifer delta13C record from the equatorial Atlantic, which exhibits significant spectral power at the 100,000- and 41,000-year Milankovitch periods, but no power at the 23,000-year period. Similar to the record presented by Shackleton and Pisias [1985], surface-deep ocean Delta delta13C produced with the G. sacculifer record leads the delta18O ice volume record. However, the glacial-interglacial amplitudes of Delta delta13C differ between our record and Shackleton and Pisias [1985] record. Although large changes in Delta delta13C occur in the equatorial Atlantic during early stages of the last three glacial cycles, surface-deep Delta delta13C at glacial maxima (18O stage 2, late stage 6, and late stage 8) was only about 0.2? greater than during the subsequent interglacial. Our results imply that nutrient-driven pCO2 changes account for about one third of the pCO2 decrease observed in ice cores, and consequently, Delta delta13C should not be used as a proxy pCO2 index. Enough variance in the ice core pCO2 records remains to be explained that conclusions about pCO2 and ice volume phase relationships should also be reexamined. As much as 40 ppm pCO2 change still has not been accounted for by models of past physics and chemistry of the ocean.

Barite accumulation rates throughout Pliocene-Pleistocene in Eastern Equatorial Pacific at ODP Site 138-849

Export production is an important component of the carbon cycle, modulating the climate system by transferring CO2 from the atmosphere to the deep ocean via the biological pump. Here we use barite accumulation rates to reconstruct export production in the eastern equatorial Pacific over the past 4.3 Ma. We find that export production fluctuated considerably on multiple time scales. Export production was on average higher (51 g C/m**2/yr) during the Pliocene than the Pleistocene (40 g C/m**2/yr), decreasing between 3 and 1 Ma (from more than 60 to 20 g C/m**2/yr) followed by an increase over the last million years. These trends likely reflect basin-scale changes in nutrient inventory and ocean circulation. Our record reveals decoupling between export production and temperatures on these long (million years) time scale. On orbital time scales, export production was generally higher during cold periods (glacial maxima) between 4.3 and 1.1 Ma. This could be due to stronger wind stress and higher upwelling rates during glacial periods. A shift in the timing of maximum export production to deglaciations is seen in the last ~1.1 million years. Results from this study suggest that, in the eastern equatorial Pacific, mechanisms that affect nutrient supply and/or ecosystem structure and in turn carbon export on orbital time scales differ from those operating on longer time scales and that processes linking export production and climate-modulated oceanic conditions changed about 1.1 million years ago. These observations should be accounted for in climate models to ensure better predictions of future climate change.

Surface wind stress modeled for the Panama Seaway with the Community Climate System Model 3

The early Pliocene warm phase was characterized by high sea surface temperatures and a deep thermocline in the eastern equatorial Pacific. A new hypothesis suggests that the progressive closure of the Panamanian seaway contributed substantially to the termination of this zonally symmetric state in the equatorial Pacific. According to this hypothesis, intensification of the Atlantic meridional overturning circulation (AMOC) - induced by the closure of the gateway - was the principal cause of equatorial Pacific thermocline shoaling during the Pliocene. In this study, twelve Panama seaway sensitivity experiments from eight ocean/climate models of different complexity are analyzed to examine the effect of an open gateway on AMOC strength and thermocline depth. All models show an eastward Panamanian net throughflow, leading to a reduction in AMOC strength compared to the corresponding closed-Panama case. In those models that do not include a dynamic atmosphere, deepening of the equatorial Pacific thermocline appears to scale almost linearly with the throughflow-induced reduction in AMOC strength. Models with dynamic atmosphere do not follow this simple relation. There are indications that in four out of five models equatorial wind-stress anomalies amplify the tropical Pacific thermocline deepening. In summary, the models provide strong support for the hypothesized relationship between Panama closure and equatorial Pacific thermocline shoaling.

Sedimentology, age models and stable isotope ratios of sediment cores from the equatorial Pacific

Based on detailed reconstructions of global distribution patterns, both paleoproductivity and the benthic d13C record of CO2, which is dissolved in the deep ocean, strongly differed between the Last Glacial Maximum and the Holocene. With the onset of Termination I about 15,000 years ago, the new (export) production of low- and mid-latitude upwelling cells started to decline by more than 2-4 Gt carbon/year. This reduction is regarded as a main factor leading to both the simultaneous rise in atmospheric CO2 as recorded in ice cores and, with a slight delay of more than 1000 years, to a large-scale gradual CO2 depletion of the deep ocean by about 650 Gt C. This estimate is based on an average increase in benthic d13C by 0.4-0.5 per mil. The decrease in new production also matches a clear 13C depletion of organic matter, possibly recording an end of extreme nutrient utilization in upwelling cells. As shown by Sarnthein et al., [1987], the productivity reversal appears to be triggered by a rapid reduction in the strength of meridional trades, which in turn was linked via a shrinking extent of sea ice to a massive increase in high-latitude insolation, i.e., to orbital forcing as primary cause.

(Table S2) The reconstructed sea surface temperature (SST) over the eastern equatorial Pacific

Abundant hydroclimatic evidence from western Amazonia and the adjacent Andes documents wet conditions during Heinrich Stadial 1 (HS1, 18-15 ka), a cold period in the high latitudes of the North Atlantic. This precipitation anomaly was attributed to a strengthening of the South American summer monsoon due to a change in the Atlantic interhemispheric sea surface temperature (SST) gradient. However, the physical viability of this mechanism has never been rigorously tested. We address this issue by combining a thorough compilation of tropical South American paleorecords and a set of atmosphere model sensitivity experiments. Our results show that the Atlantic SST variations alone, although leading to dry conditions in northern South America and wet conditions in northeastern Brazil, cannot produce increased precipitation over western Amazonia and the adjacent Andes during HS1. Instead, an eastern equatorial Pacific SST increase (i.e., 0.5-1.5 °C), in response to the slowdown of the Atlantic Meridional Overturning Circulation during HS1, is crucial to generate the wet conditions in these regions. The mechanism works via anomalous low sea level pressure over the eastern equatorial Pacific, which promotes a regional easterly low-level wind anomaly and moisture recycling from central Amazonia towards the Andes.

Dinocyst assemblage constraints on oceanographic and atmospheric processes in the eastern equatorial Atlantic over the last 44 kyr

A new 44 kyr long record of dinoflagellate (phytoplanktonic organisms) cysts (dinocysts) is presented from a marine sediment core collected on the Congolese margin with the aim of reconstructing past hydrological changes in the equatorial eastern Atlantic Ocean since Marine Isotopic Stage (MIS) 3. Our high-resolution dinocyst record indicates that significant temperature and moisture variations occurred across the glacial period, the last deglaciation and the Holocene. The use of specific dinocyst taxa, indicative of fluvial, upwelling and Benguela Current past environments for instance, provides insights into the main forcing mechanisms controlling palaeohydrological changes on orbital timescales. In particular, we are able, for the last 44 kyr, to correlate fluvial-sensitive taxa to monsoonal mechanisms related to precession minima–obliquity maxima combinations. While upwelling mechanisms appear as the main drivers for dinoflagellate productivity during MIS 2, dissolved nutrient-enriched Congo River inputs to the ocean also played a significant role in promoting dinoflagellate productivity between approximately 15.5 and 5 ka BP. Finally, this high-resolution dinocyst study permits us to precisely investigate the suborbital timing of the last glacial–interglacial termination, including an atypical warm and wet oceanic LGM signature, northern high-latitude abrupt climate change impacts in the equatorial eastern Atlantic, as well as a two-step decrease in moisture conditions during the Holocene at around 7–6 and 4–3.5 ka BP.