Planktonic and benthic foraminiferal stable isotopes, SSTs and thermocline temperature, coccoliths and benthic foraminiferal assemblages, and productivity reconstruction for sediment core MD06-3067

We present sea surface, upper thermocline, and benthic d18O data, as well as temperature and paleoproductivity proxy data, from the International Marine Global Change Study Program (IMAGES) Core MD06-3067 (6°31'N, 126°30'E, 1575 m water depth), located in the western equatorial Pacific Ocean within the flow path of the Mindanao Current. Our records reveal considerable glacial-interglacial and suborbital variability in the Mindanao Dome upwelling over the last 160 kyr. Dome activity generally intensified during glacial intervals resulting in cooler thermocline waters, whereas it substantially declined during interglacials, in particular in the early Holocene and early marine oxygen isotope stage (MIS) 5e, when upwelling waters did not reach the thermocline. During MIS 3 and MIS 2, enhanced surface productivity together with remarkably low SST and low upper ocean thermal contrast provide evidence for episodic glacial upwelling to the surface, whereas transient surface warming marks periodic collapses of the Mindanao Dome upwelling during Heinrich events. We attribute the high variability during MIS 3 and MIS 2 to changes in the El Niño Southern Oscillation state that affected boreal winter monsoonal winds and upper ocean circulation. Glacial upwelling intensified when a strong cyclonic gyre became established, whereas El Niño-like conditions during Heinrich events tended to suppress the cyclonic circulation, reducing Ekman transport. Thus, our findings demonstrate that variations in the Mindanao Dome upwelling are closely linked to the position and intensity of the tropical convection and also reflect far-field influences from the high latitudes.

Stable carbon and oxygen isotope records of Cibicidoides wuellerstorfi in DSDP and ODP sites in the Atlantic and Pacific Ocean

Carbon isotopic records from benthic foraminifera are used to map patterns of deep ocean circulation between 3 and 2 million years ago, the interval when significant northern hemisphere glaciation began. The delta18O and delta13C data from four Atlantic sites (552, 607, 610, and 704) and one Pacific site (677) show that global cooling over this interval was associated with increased suppression of North Atlantic Deep Water (NADW) formation. However, the relative strength of NADW production was always greater than is observed during late Pleistocene glaciations when extreme decreases in NADW are observed in the deep North Atlantic. Our data indicate that an increase in the equator-to-pole temperature gradient associated with the onset of northern hemisphere glaciation did not intensify deepwater production in the North Atlantic but rather the opposite occurred. This is not unexpected as it is the "warm high-salinity" characteristic, rather than the "low temperature", of thermocline waters that is critical to the deepwater formation process in this region today.

Magnesium/Calcium ratios and sea surface temperature estimation for sediments of the Timor Sea

We use a multiproxy approach to monitor changes in the vertical profile of the Indonesian Throughflow as well as monsoonal wind and precipitation patterns in the Timor Sea on glacial-interglacial, precessional, and suborbital timescales. We focus on an interval of extreme climate change and sea level variation: marine isotope (MIS) 6 to MIS 5e. Paleoproductivity fluctuations in the Timor Sea follow a precessional beat related to the intensity of the Australian (NW) monsoon. Paired Mg/Ca and d18O measurements of surface- and thermocline-dwelling planktonic foraminifers (G. ruber and P. obliquiloculata) indicate an increase of >4°C in both surface and thermocline water temperatures during Termination II. Tropical sea surface temperature changed synchronously with ice volume (benthic d18O) during deglaciation, implying a direct coupling of high- and low-latitude climate via atmospheric and/or upper ocean circulation. Substantial cooling and freshening of thermocline waters occurred toward the end of Termination II and during MIS 5e, indicating a change in the vertical profile of the Indonesian Throughflow from surface- to thermocline-dominated flow.

Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Holocene and the Last Glacial Maximum

We measured oxygen isotopes and Mg/Ca ratios in the surface-dwelling planktonic foraminifer Globigerinoides ruber (white s.s.) and the thermocline dweller Pulleniatina obliquiloculata to investigate upper ocean spatial variability in the Indo-Pacific Warm Pool (IPWP). We focused on three critical time intervals: the Last Glacial Maximum (LGM; 18-21.5 ka), the early Holocene (8-9 ka), and the late Holocene (0-2 ka). Our records from 24 stations in the South China Sea, Timor Sea, Indonesian seas, and western Pacific indicate overall dry and cool conditions in the IPWP during the LGM with a low thermal gradient between surface and thermocline waters. During the early Holocene, sea surface temperatures increased by ~3°C over the entire region, indicating intensification of the IPWP. However, in the eastern Indian Ocean (Timor Sea), the thermocline gradually shoaled from the LGM to early Holocene, reflecting intensification of the subsurface Indonesian Throughflow (ITF). Increased surface salinity in the South China Sea during the Holocene appears related to northward displacement of the monsoonal rain belt over the Asian continent together with enhanced influx of saltier Pacific surface water through the Luzon Strait and freshwater export through the Java Sea. Opening of the freshwater portal through the Java Sea in the early Holocene led to a change in the vertical structure of the ITF from surface- to thermocline-dominated flow and to substantial freshening of Timor Sea thermocline waters.

Late Quaternary planktonic foraminifera of ODP Hole 175-1087A

Planktonic foraminiferal assemblages from the upper Pleistocene part of Hole 1087A (0 to 12.1 meters below seafloor) are investigated to assess the role of global and local climate changes on surface circulation in the southern Benguela region. The benthic stable isotope record indicates that the studied interval is representative of the last four climatic cycles, that is, down to marine isotope Stage (MIS) 12. The species assemblages bear a clear transitional to subpolar character, with Neogloboquadrina pachyderma (d), Globorotalia inflata, and Globigerina bulloides, in order of decreasing abundance, as the dominant taxa. This species association presently characterizes the mixing domain of old upwelled and open ocean waters, seaward of the Benguela upwelling cells. Abundance variation of the dominant foraminiferal species roughly follows a glacial-interglacial pattern down to MIS 8, suggesting an alternation of upwelling strength and associated seaward extension of the belt of upwelled water as a response to global climate changes. This pattern is interrupted from ~250 ka down to MIS 12, where the phase relationship with global climate is ill defined and might be interpreted as a local response of the southern Benguela region to the mid-Brunhes event. Of particular interest is a single pulse of newly upwelled waters at the location of Site 1087 during early MIS 9 as indicated by a peak abundance of sinistral N. pachyderma (s). Variable input of warm, salty Indian Ocean thermocline waters into the southeast Atlantic, a key component of the Atlantic heat conveyor, is indicated by abundance changes of the tropical taxon Globorotalia menardii. From this tracer, we suggest that interocean exchange was hardly interrupted throughout the last 460 k.y., but was most effective at glacial terminations, particularly during Terminations I and II, as well as during the upper part of MIS 12. This maximum input of Indian Ocean waters around the southern tip of Africa is associated with the reseeding of G. menardii in the tropical Atlantic.

Benthic foraminiferal d18O and Mg/Ca studies of sediment cores off tropical NW Africa

Benthic foraminiferal d18O and Mg/Ca of sediment cores off tropical NW Africa are used to study the properties of Atlantic central waters during the Last Glacial Maximum (LGM) and Heinrich Stadial 1 (HS1). We combined our core top data with published results to develop a new Mg/Ca-temperature calibration for Planulina ariminensis, which shows a Mg/Ca-temperature sensitivity of 0.19 mmol/mol per °C. Estimates of the LGM and HS1 thermocline temperatures are comparable to the present-day values between 200 and 400 m water depth, but were 1.2-1.5°C warmer at 550-570 m depth. The HS1 thermocline waters (200-570 m depth) did not show any warming relative to the LGM. This is in contrast to previous climate model studies, which concluded that tropical Atlantic thermocline waters warmed significantly when Atlantic meridional overturning circulation was reduced. However, our results suggest that thermocline temperatures of the northeastern tropical Atlantic show no pronounced sensitivity to changes in the thermohaline circulation during glacial periods. In contrast, we find a significant increase in thermocline-water salinity during the LGM (200-550 m depth) and HS1 (200-400 m depth) with respect to the present-day, which we relate to changes in the wind-driven circulation. We infer that the LGM thermocline (200-550 m depth) and the HS1 upper thermocline (200-400 m depth) in the northeastern tropical Atlantic was ventilated by surface waters from the North Atlantic rather than the southern-sourced waters. This suggests that the frontal zone between the modern South Atlantic and North Atlantic Central Waters was probably shifted southward during the LGM and HS1.

Age determination and stable isotope record of sediment cores at 9°S in the Indian Ocean

A continuous 10-m-long section consisting of roughly two thirds Ethmodiscus rex (a diatom) and one third mixed planktonic foraminifera was identified in a core from 3800 m depth at 9°S on the Indian Ocean's 90°E Ridge. Radiocarbon dates place the onset of deposition of this layer at >30,000 years B.P. and its termination at close to 11,000 years B.P. However, precise dating of the foraminifera from the Ethmodiscus layer itself proved to be impossible owing to the presence of secondary calcite presumably precipitated from the pore waters. During the Holocene, high calcium carbonate content ooze free of diatoms was deposited at this locale. As the site currently lies beneath the pathway taken by upper ocean waters entering the Indian Ocean from the Pacific (via the Indonesian Straits), it appears that during glacial time, thermocline waters moving along this same path provided the silica and other nutrients required by these diatoms.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Southern Ocean

Recent evidence that dissolved organic carbon (DOC) is a significant component of the organic carbon flux below the photic layer of the ocean (1), together with verification of high respiration rates in the dark ocean (2), suggests that the downward flux of DOC may play a major role in supporting respiration there. Here we show, on the basis of examination of the relation between DOC and apparent oxygen utilization (AOU), that the DOC flux supports ~10% of the respiration in the dark ocean. The contribution of DOC to pelagic respiration below the surface mixed layer can be inferred from the relation between DOC and apparent oxygen utilization (AOU, µM O2), a variable quantifying the cumulative oxygen consumption since a water parcel was last in contact with the atmosphere. However, assessments of DOC/AOU relations have been limited to specific regions of the ocean (3, 4) and have not considered the global ocean. We assembled a large data set (N = 9824) of concurrent DOC and AOU observations collected in cruises conducted throughout the world's oceans (fig. S1, table S1) to examine the relative contribution of DOC to AOU and, therefore, respiration in the dark ocean. AOU increased from an average (±SE) 96.3 ± 2.0 µM at the base of the surface mixed layer (100 m) to 165.5 ± 4.3 µM at the bottom of the main thermocline (1000 m), with a parallel decline in the average DOC from 53.5 ± 0.2 to 43.4 ± 0.3 µM C (Fig. 1). In contrast, there is no significant decline in DOC with increasing depth beyond 1000 m depth (Fig. 1), indicating that DOC exported with overturning circulation plays a minor role in supporting respiration in the ocean interior (5). Assuming a molar respiratory quotient of 0.69, the decline in DOC accounts for 19.6 ± 0.4% of the AOU within the top 1000 m (Fig. 1). This estimate represents, however, an upper limit, because the correlation between DOC and AOU is partly due to mixing of DOC-rich warm surface waters with DOC-poor cold thermocline waters (6). Removal of this effect by regressing DOC against AOU and water temperature indicates that DOC supports only 8.4 ± 0.3% of the respiration in the mesopelagic waters.

Age model for ODP Leg 175 sites

Upwelling along the western coast of Africa south of the equator may be partitioned into three major areas, each having its own dynamics and history: (1) the eastern equatorial region, comprising the Congo Fan and the area of Mid-Angola; (2) the Namibia upwelling system, extending from the Walvis Ridge to Lüderitz; and (3) the Cape Province region, where upwelling is subdued. The highest nutrient contents in thermocline waters are in the northern region, the lowest in the southern one. Wind effects are at a maximum near the southern end of the Namibia upwelling system, and maximum productivity occurs near Walvis Bay, where the product between upwelling rate and nutrient content of upwelled waters is at a maximum. In the Congo/Angola region, opal tends to follow organic carbon quite closely in the Quaternary record. However, organic carbon has a strong precessional component, while opal does not. Despite relatively low opal content, sediments off Angola show the same patterns as those off the Congo; thus, they are part of the same regime. The spectrum shows nonlinear interference patterns between high- and low-latitude forcing, presumably tied to thermocline fertility and wind. On Walvis Ridge, as in the Congo-Angola region, the organic matter record behaves normally; that is, supply is high during glacial periods. In contrast, interglacial periods are favorable for opal deposition. The pattern suggests reduction in silicate content of the thermocline during glacial periods. The reversed phase (opal abundant during interglacials) persists during the entire Pleistocene and can be demonstrated deep into the Pliocene, not just on Walvis Ridge but all the way to the Oranje River and off the Cape Province. From comparison with other regions, it appears that silicate is diminished in the global thermocline, on average, whenever winds become strong enough to substantially shorten the residence time of silicate in upper waters (Walvis Hypothesis, solving the Walvis Paradox of reversed phase in opal deposition). The central discovery during Leg 175 was the documentation of a late Pliocene opal maximum for the entire Namibia upwelling system (early Matuyama Diatom Maximum [MDM]). The maximum is centered on the period between the end of the Gauss Chron and the beginning of the Olduvai Chron. A rather sharp increase in both organic matter deposition and opal deposition occurs near 3 Ma in the middle of the Gauss Chron, in association with a series of major cooling steps. As concerns organic matter, high production persists at least to 1 Ma, when there are large changes in variability, heralding subsequent pulsed production periods. From 3 to 2 Ma, organic matter and opal deposition run more or less parallel, but after 2 Ma opal goes out of phase with organic matter. Apparently, this is the point when silicate becomes limiting to opal production. Thus, the MDM conundrum is solved by linking planetary cooling to increased mixing and upwelling (ramping up to the MDM) and a general removal of silicate from the upper ocean through excess precipitation over global supply (ramping down from the MDM). The hypothesis concerning the origin of the Namibia opal acme or MDM is fundamentally the same as the Walvis Hypothesis, stating that glacial conditions result in removal of silicate from the thermocline (and quite likely from the ocean as a whole, given enough time). The Namibia opal acme, and other opal maxima in the latest Neogene in other regions of the ocean, marks the interval when a cooling ocean selectively removes the abundant silicate inherited from a warm ocean. When the excess silicate is removed, the process ceases. According to the data gathered during Leg 175, major upwelling started in the late part of the late Miocene. Presumably, this process contributed to the drawing down of carbon dioxide from the atmosphere, helping to prepare the way for Northern Hemisphere glaciation.

Foraminiferal Mg/Ca and Mn/Ca ratios across the Eocene-Oligocene transition

Constraining the magnitude of high-latitude temperature change across the Eocene-Oligocene transition (EOT) is essential for quantifying the magnitude of Antarctic ice-sheet expansion and understanding regional climate response to this event. To this end, we constructed high-resolution stable oxygen isotope (d18O) and magnesium/calcium (Mg/Ca) records from planktic and benthic foraminifera at four Ocean Drilling Program (ODP) sites in the Southern Ocean. Planktic foraminiferal Mg/Ca records from the Kerguelen Plateau (ODP Sites 738, 744, and 748) show a consistent pattern of temperature change, indicating 2-3 °C cooling in direct conjunction with the first step of a two-step increase in benthic and planktic foraminiferal d18O values across the EOT. In contrast, benthic Mg/Ca records from Maud Rise (ODP Site 689) and the Kerguelen Plateau (ODP Site 748) do not exhibit significant temperature change. The contrasting temperature histories derived from the planktic and benthic Mg/Ca records are not reconcilable, since vertical d18O gradients remained nearly constant at all sites between 35.0 and 32.5 Ma. Based on the coherency of the planktic Mg/Ca records from the Kerguelen Plateau sites and complications with benthic Mg/Ca paleothermometry at low temperatures, the planktic Mg/Ca records are deemed the most reliable measure of Southern Ocean temperature change. We therefore interpret a uniform cooling of 2-3 °C in both deep surface (thermocline) waters and intermediate deep waters of the Southern Ocean across the EOT. Cooling of Southern Ocean surface waters across the EOT was likely propagated to the deep ocean, since deep waters were primarily sourced on the Antarctic margin throughout this time interval. Removal of the temperature component from the observed foraminiferal d18O shift indicates that seawater d18O values increased by 0.6 ± 0.15 per mil across the EOT interval, corresponding to an increase in global ice volume to a level equivalent with 60-130% modern East Antarctic ice sheet volume.

(Table 1) Spacing of diatom abundance index and opal index at ODP Site 175-1084

An important discovery during Ocean Drilling Program Leg 175, when investigating the record of upwelling off Namibia, was the finding of a distinct Late Pliocene diatom maximum spanning the lower half of the Matuyama reversed polarity chron (MDM, Matuyama Diatom Maximum) and centered around 2.6-2.0 Ma. This maximum was observed at all sites off southwestern Africa between 20°S and 30°S, and is most strongly represented in sediments of Site 1084, off Lüderitz, Namibia. The MDM is characterized by high biogenic opal content, high numbers of diatom valves, and a diatom flora rich in Southern Ocean representatives (with Thalassiothrix antarctica forming diatom mats) as well as coastal upwelling components. Before MDM time, diatoms are rare until ca. 3.6 Ma. After the MDM, in the Pleistocene, the composition of the diatom flora points to increased importance of coastal upwelling toward the present, but is accompanied by a general decrease in opal and diatom deposition. Here we present a simple conceptual model as a first step in formalizing a possible forcing mechanism responsible for the record of opal deposition in the upwelling system off Namibia. The model takes into account Southern Ocean oceanography, and a link with deepwater circulation and deepwater nutrient chemistry which, in turn, are coupled to the evolution of North Atlantic Deep Water (NADW). The model proposes that between the MDM and the Mid-Pleistocene climate revolution, opal deposition off Namibia is not directly tied to glacial-interglacial fluctuations (as seen in the global d18O record), but that, instead, a strong deepwater link exists with increased NADW production (as seen in the deepwater d13C record) accounting for higher supply of silicate to the thermocline waters that feed the upwelling process. The opal record of Site 1084 shows affinity to eccentricity on the 400-kyr scale but not for the 100-kyr scale. This points toward long-term geologic processes for delivery of silica to the ocean.

(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.

Stable isotope record, Mg/Ca ratios and sea surface reconstruction for sediment core MD06-3067

We present sea surface, upper thermocline, and benthic d18O data, as well as temperature and paleoproductivity proxy data, from the International Marine Global Change Study Program (IMAGES) Core MD06-3067 (6°31'N, 126°30'E, 1575 m water depth), located in the western equatorial Pacific Ocean within the flow path of the Mindanao Current. Our records reveal considerable glacial-interglacial and suborbital variability in the Mindanao Dome upwelling over the last 160 kyr. Dome activity generally intensified during glacial intervals resulting in cooler thermocline waters, whereas it substantially declined during interglacials, in particular in the early Holocene and early marine oxygen isotope stage (MIS) 5e, when upwelling waters did not reach the thermocline. During MIS 3 and MIS 2, enhanced surface productivity together with remarkably low SST and low upper ocean thermal contrast provide evidence for episodic glacial upwelling to the surface, whereas transient surface warming marks periodic collapses of the Mindanao Dome upwelling during Heinrich events. We attribute the high variability during MIS 3 and MIS 2 to changes in the El Niño Southern Oscillation state that affected boreal winter monsoonal winds and upper ocean circulation. Glacial upwelling intensified when a strong cyclonic gyre became established, whereas El Niño-like conditions during Heinrich events tended to suppress the cyclonic circulation, reducing Ekman transport. Thus, our findings demonstrate that variations in the Mindanao Dome upwelling are closely linked to the position and intensity of the tropical convection and also reflect far-field influences from the high latitudes.

Pliocene-Pleistocene carbon isotope record of foraminifera from DSDP Site 89-586

Oceanographic changes in the western equatorial Pacific during the past 6 m.y. are inferred from carbon isotopic analyses of planktonic and benthic foraminifers from Ontong Java Plateau (DSDP Site 586). Sample spacing is 1.5 m (ca. 35,000-75,000 yr). An overall trend of d13C toward lighter values is evident for the last 5 m.y. in all four foraminiferal taxa analyzed (G. sacculifer, Pulleniatina, P. wuellerstorfi, and O. umbonatus). This trend is interpreted as an enrichment of the global ocean with 12C, because of the addition of carbon from organic carbon reservoirs (or lack of removal of carbon to such reservoirs), as a consequence of an overall drop in sea level. Differences between shallow- and deep-water d13C decrease slightly during this time interval, suggesting a moderate drop in productivity. This drop is not sufficient to explain the drop in sedimentation rate, however, much of which apparently must be ascribed to winnowing effects. A marked convergence in the d13C values of planktonic taxa exists within the last 2 m.y. We propose that this convergence indicates nutrient depletion in thermocline waters, caused by the vigorous removal of phosphate in marginal upwelling regions, or by the stripping of intermediate waters in their source regions. No large shifts are seen in the carbon isotope record of the last 6 m.y., in contrast to the oxygen isotope record. Some indication of cyclicity is present, with a period between 0.5 and 1.0 m.y. (especially in the earlier portion of the record).