Upper ocean vertical supply: A neglected primary factor controlling the distribution of neodymium concentrations of open ocean surface waters?

Neodymium (Nd) isotopes are an important geochemical tool to trace the present and past water mass mixing as well as continental inputs. The distribution of Nd concentrations in open ocean surface waters (0�100 m) is generally assumed to be controlled by lateral mixing of Nd from coastal surface currents and by removal through reversible particle scavenging. However, using 228Ra activity as an indicator of coastal water mass influence, surface water Nd concentration data available on key oceanic transects as a whole do not support the above scenario. From a global compilation of available data, we find that more stratified regions are generally associated with low surface Nd concentrations. This implies that upper ocean vertical supply may be an as yet neglected primary factor in determining the basin-scale variations of surface water Nd concentrations. Similar to the mechanism of nutrients supply, it is likely that stratification inhibits vertical supply of Nd from the subsurface thermocline waters and thus the magnitude of Nd flux to the surface layer. Consistently, the estimated required input flux of Nd to the surface layer to maintain the observed concentrations could be nearly two orders of magnitudes larger than riverine/dust flux, and also larger than the model-based estimation on shelf-derived coastal flux. In addition, preliminary results from modeling experiments reveal that the input from shallow boundary sources, riverine input, and release from dust are actually not the primary factors controlling Nd concentrations most notably in the Pacific and Southern Ocean surface waters.

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.

Biogeochemical cycling of dissolved zinc along the GEOTRACES South Atlantic transect GA10 at 40°S

The biogeochemical cycle of zinc (Zn) in the South Atlantic, at 40°S, was investigated as part of the UK GEOTRACES program. To date there is little understanding of the supply of Zn, an essential requirement for phytoplankton growth, to this highly productive region. Vertical Zn profiles displayed nutrient-like distributions with distinct gradients associated with the watermasses present. Surface Zn concentrations are among the lowest reported for theworld’s oceans (<50 pM). A strong Zn-Si linear relationshipwas observed (Zn (nM)= 0.065 Si (μM), r2=0.97, n = 460). Our results suggest that the use of a global Zn-Si relationship would lead to an underestimation of dissolved Zn in deeper waters of the South Atlantic. By utilizing Si* and a new tracer Zn* our data indicate that the preferential removal of Zn in the Southern Ocean prevented a direct return path for dissolved Zn to the surface waters of the South Atlantic at 40°S and potentially the thermocline waters of the South Atlantic subtropical gyre. The importance of Zn for phytoplankton growth was evaluated using the Zn-soluble reactive phosphorus (SRP) relationship. We hypothesize that the low Zn concentrations in the South Atlantic may select for phytoplankton cells with a lower Zn requirement. In addition, a much deeper kink at ~ 500m in the Zn:SRP ratio was observed compared to other oceanic regions.

Changes in the intermediate water mass formation rates in the global ocean for the Last Glacial Maximum, mid-Holocene and pre-industrial climates

The paleoclimate version of the National Center for Atmospheric Research Community Climate System Model version 3 (NCAR-CCSM3) is used to analyze changes in the water formation rates in the Atlantic, Pacific, and Indian Oceans for the Last Glacial Maximum (LGM), mid-Holocene (MH) and pre-industrial (PI) control climate. During the MH, CCSM3 exhibits a north-south asymmetric response of intermediate water subduction changes in the Atlantic Ocean, with a reduction of 2 Sv in the North Atlantic and an increase of 2 Sv in the South Atlantic relative to PI. During the LGM, there is increased formation of intermediate water and a more stagnant deep ocean in the North Pacific. The production of North Atlantic Deep Water (NADW) is significantly weakened. The NADW is replaced in large extent by enhanced Antarctic Intermediate Water (AAIW), Glacial North Atlantic Intermediate Water (GNAIW), and also by an intensified of Antarctic Bottom Water (AABW), with the latter being a response to the enhanced salinity and ice formation around Antarctica. Most of the LGM intermediate/mode water is formed at 27.4 < sigma(theta) < 29.0 kg/m(3), while for the MH and PI most of the subduction transport occurs at 26.5 < sigma(theta) < 27.4 kg/m(3). The simulated LGM Southern Hemisphere winds are more intense by 0.2-0.4 dyne/cm(2). Consequently, increased Ekman transport drives the production of intermediate water (low salinity) at a larger rate and at higher densities when compared to the other climatic periods.

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.