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.

Habitat characteristics and energy budget of the ocean quahog (Arctica islandica)

We compared lifetime and population energy budgets of the extraordinary long-lived ocean quahog Arctica islandica from 6 different sites - the Norwegian coast, Kattegat, Kiel Bay, White Sea, German Bight, and off northeast Iceland - covering a temperature and salinity gradient of 4-10°C (annual mean) and 25-34, respectively. Based on von Bertalanffy growth models and size-mass relationships, we computed organic matter production of body (PSB) and of shell (PSS), whereas gonad production (PG) was estimated from the seasonal cycle in mass. Respiration (R) was computed by a model driven by body mass, temperature, and site. A. islandica populations differed distinctly in maximum life span (40 y in Kiel Bay to 197 y in Iceland), but less in growth performance (phi' ranged from 2.41 in the White Sea to 2.65 in Kattegat). Individual lifetime energy throughput, as approximated by assimilation, was highest in Iceland (43,730 kJ) and lowest in the White Sea (313 kJ). Net growth efficiency ranged between 0.251 and 0.348, whereas lifetime energy investment distinctly shifted from somatic to gonad production with increasing life span; PS/PG decreased from 0.362 (Kiel Bay, 40 y) to 0.031 (Iceland, 197 y). Population annual energy budgets were derived from individual budgets and estimates of population mortality rate (0.035/y in Iceland to 0.173/y in Kiel Bay). Relationships between budget ratios were similar on the population level, albeit with more emphasis on somatic production; PS/ PG ranged from 0.196 (Iceland) to 2.728 (White Sea), and P/B ranged from 0.203-0.285/y. Life span is the principal determinant of the relationship between budget parameters, whereas temperature affects net growth efficiency only. In the White Sea population, both growth performance and net growth efficiency of A. islandica were lowest. We presume that low temperature combined with low salinity represent a particularly stressful environment for this species.

Age model and stable isotope record of DSDP holes 48-401 and 86-577

Early Paleogene warm climates may have been linked to different modes and sources of deepwater formation. Warm polar temperatures of the Paleocene and Eocene may have resulted from either increased atmospheric trace gases or increased heat transport through deep and intermediate waters. The possibility of increasing ocean heat transport through the production of warm saline deep waters (WSDW) in the Tethyan region has generated considerable interest. In addition, General Circulation Model results indicate that deepwater source regions may be highly sensitive to changing basin configurations. To decipher deepwater changes, we examined detailed benthic foraminiferal faunal and isotopic records of the late Paleocene through the early Eocene (~60 to 50 Ma) from two critical regions: the North Atlantic (Bay of Biscay Site 401) and the Pacific (Shatsky Rise Site 577). These records are compared with published data from the Southern Ocean (Maud Rise Site 690, Islas Orcadas Rise Site 702). During the late Paleocene, similar benthic foraminiferal delta18O values were recorded at all four sites. This indicates uniform deepwater temperatures, consistent with a single source of deep water. The highest delta13C values were recorded in the Southern Ocean and were 0.5 per mil more positive than those of the Pacific. We infer that the Southern Ocean was proximal to a source of nutrient-depleted deep water during the late Paleocene. Upper Paleocene Reflector Ab was cut on the western Bermuda Rise by cyclonically circulating bottom water, also suggesting a vigorous source of bottom water in the Southern Ocean. A dramatic negative excursion in both carbon and oxygen isotopes occurred in the latest Paleocene in the Southern Ocean. This is a short-term (<100 kyr), globally synchronous event which also is apparent in both the Atlantic and Pacific records as a carbon isotopic excursion of approximately 1 per mil. Faunal analyses from the North Atlantic and Pacific sites indicate that the largest benthic foraminiferal faunal turnover of the Cenozoic was synchronous with the isotopic excursion, lending support to the hypothesis that the extinctions were caused by a change in deepwater circulation. We speculate that the Southern Ocean deepwater source was reduced or eliminated at the time of the excursion. During the early Eocene, Southern Ocean delta13C values remained enriched relative to the North Atlantic and Pacific. However, the Southern Ocean was also enriched in delta18O relative to these basins. We interpret that these patterns indicate that although the Southern Ocean was proximal to a source of cool, nutrient-depleted water, the intermediate to upper deep water sites of the North Atlantic and Pacific were ventilated by a different source that probably originated in low latitudes, i.e., WSDW.

Geochemistry of sediments in the central equatorial Pacific Ocean

In order to quantify changes in export production and carbonate dissolution over the past 1 Myr in the central equatorial Pacific Ocean we analyzed Ba, P, Al Ti, and Ca in 1106 samples from five piston cores gathered from 5°S to 4°N at 140°W. We focused on Ba/Ti, Al/Ti, and P/Ti ratios as export proxies and employed areally integrated time slice as well as time series strategies. Carbonate maxima from 0-560 kyr are characterized by 15-30% greater export than carbonate minima. The increases in export fall on glacial delta18O transitions rather than glacial maxima. From 560-800 kyr, overlapping with the mid-Pleistocene transition, there is a very large increase in total export yet no glacial-interglacial variability. The highest latitudes (5°S and 4°N) record minimal absolute export change from glacials to interglacials and yet record the most extreme minima in percent CaCO3, indicating that carbonate records there are dominated by dissolution, whereas near the equator they are more influenced by changes in export.

Oligocene strontium isotope stratigraphy of Indian Ocean sediments

Strontium isotope (87Sr/86Sr) ages have been established for Oligocene samples of Leg 119 Site 744, Leg 120 Sites 747 and 748, and Leg 121 Sites 756 and 757. Ages were determined using the strontium isotope age equation of Miller et al. (1988) and preliminary correlations have been made with available nannofossil biostratigraphy. The strontium isotope ages calculated here augment biostratigraphy, which for the Oligocene is characterized by long biozones, and provide additional detail where the paleomagnetic record is not clear (Sites 756 and 757). Results from the lower latitude Ninetyeast Ridge sites where standard calcareous nannofossil datums are present are compared to those of the higher latitude Kerguelen Plateau sites in order to examine biostratigraphic events across latitude in the Indian Ocean. The 87Sr/86Sr determined ages are used here as a tool for correlation.

Global Ocean Particulate Organic Carbon flux merged with satellite parameters

The efficiency of the biological pump of carbon to the deep ocean depends largely on the biologically mediated export of carbon from the surface ocean and its remineralization with depth. Global satellite studies have primarily focused on chlorophyll concentration and net primary production (NPP) to understand the role of phytoplankton in these processes. Recent satellite retrievals of phytoplankton composition now allow for the size of phytoplankton cells to be considered. Here, we improve understanding of phytoplankton size structure impacts on particle export, remineralization and transfer. Particulate organic carbon (POC) flux observations from sediment traps and 234Th are compiled across the global ocean. Annual climatologies of NPP, percent microplankton, and POC flux at four time series locations and within biogeochemical provinces are constructed, and sinking velocities are calculated to align surface variables with POC flux at depth. Parameters that characterize POC flux vs. depth (export flux ratio, labile fraction, remineralization length scale) are then fit to the aligned dataset. Times of the year dominated by different size compositions are identified and fit separately in regions of the ocean where phytoplankton cell size showed enough dynamic range over the annual cycle. Considering all data together, our findings support the paradigm of high export flux but low transfer efficiency in more productive regions and vice versa for oligotrophic regions. However, when parsing by dominant size class, we find periods dominated by small cells to have both greater export flux and lower transfer efficiency than periods when large cells comprise a greater proportion of the phytoplankton community.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

Macrobenthic abundance, biomass, productivity and production in the deep Arctic ocean

The few existing studies on macrobenthic communities of the deep Arctic Ocean report low standing stocks, and confirm a gradient with declining biomass from the slopes down to the basins as commonly reported for deep-sea benthos. In this study we have further investigated the relationship of faunal abundance (N), biomass (B) as well as community production (P) with water depth, geographical latitude and sea ice concentration. The underlying dataset combines legacy data from the past 20 years, as well as recent field studies selected according to standardized quality control procedures. Community P/B and production were estimated using the multi-parameter ANN model developed by Brey (2012). We could confirm the previously described negative relationship of water depth and macrofauna standing stock in the Arctic deep-sea. Furthermore, the sea-ice cover increasing with high latitudes, correlated with decreasing abundances of down to < 200 individuals/m**2, biomasses of < 65 mg C/m**2 and P of < 75 mg C/m**2/y. Stations under influence of the seasonal ice zone (SIZ) showed much higher standing stock and P means between 400 - 1400 mg C/m**2/y; even at depths up to 3700 m. We conclude that particle flux is the key factor structuring benthic communities in the deep Arctic ocean, explaining both the low values in the ice-covered Arctic basins and the high values along the SIZ.

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.

Model results of sensitivity experiments for marine nitrogen cycle in four NetCDF files

The marine nitrogen (N) inventory is thought to be stabilized by negative feedback mechanisms that reduce N inventory excursions relative to the more slowly overturning phosphorus inventory. Using a global biogeochemical ocean circulation model we show that negative feedbacks stabilizing the N inventory cannot persist if a close spatial association of N2 fixation and denitrification occurs. In our idealized model experiments, nitrogen deficient waters, generated by denitrification, stimulate local N2 fixation activity. But, because of stoichiometric constraints, the denitrification of newly fixed nitrogen leads to a net loss of N. This can enhance the N deficit, thereby triggering additional fixation in a vicious cycle, ultimately leading to a runaway N loss. To break this vicious cycle, and allow for stabilizing negative feedbacks to occur, inputs of new N need to be spatially decoupled from denitrification. Our idealized model experiments suggest that factors such as iron limitation or dissolved organic matter cycling can promote such decoupling and allow for negative feedbacks that stabilize the N inventory. Conversely, close spatial co-location of N2 fixation and denitrification could lead to net N loss.

Age model, stable isotopes, major element proportions and endmember unmixing results from four western tropical Atlantic cores

We investigate changes in the delivery and oceanic transport of Amazon sediments related to terrestrial climate variations over the last 250 ka. We present high-resolution geochemical records from four marine sediment cores located between 5 and 12° N along the northern South American margin. The Amazon River is the sole source of terrigenous material for sites at 5 and 9° N, while the core at 12° N receives a mixture of Amazon and Orinoco detrital particles. Using an endmember unmixing model, we estimated the relative proportions of Amazon Andean material ("%-Andes", at 5 and 9° N) and of Amazon material ("%-Amazon", at 12° N) within the terrigenous fraction. The %-Andes and %-Amazon records exhibit significant precessional variations over the last 250 ka that are more pronounced during interglacials in comparison to glacial periods. High %-Andes values observed during periods of high austral summer insolation reflect the increased delivery of suspended sediments by Andean tributaries and enhanced Amazonian precipitation, in agreement with western Amazonian speleothem records. Increased Amazonian rainfall reflects the intensification of the South American monsoon in response to enhanced land-ocean thermal gradient and moisture convergence. However, low %-Amazon values obtained at 12° N during the same periods seem to contradict the increased delivery of Amazon sediments. We propose that reorganizations in surface ocean currents modulate the northwestward transport of Amazon material. In agreement with published records, the seasonal North Brazil Current retroflection is intensified (or prolonged in duration) during cold substages of the last 250 ka (which correspond to intervals of high DJF or low JJA insolation) and deflects eastward the Amazon sediment and freshwater plume.

Age model of ODP Site 162-983

At ODP Site 983, relative geomagnetic paleointensity and planktic and benthic delta18O records have been acquired for the last 350 kyr. The mean sedimentation rate in this interval is 11.3 cm/kyr. Magnetic properties and hysteresis ratios indicate that pseudo-single domain magnetite is the remanence carrier. Volume susceptibility (kappa), anhysteretic (ARM) and isothermal (IRM) remanence values vary by a factor of 3-4, well within the criteria usually cited for paleointensity studies. Natural remanent magnetization (NRM) is normalized by ARM and IRM to acquire the paleointensity proxy. Arithmetic means of NRM/ARM and NRM/IRM, calculated for five demagnetization steps in the 25-45 mT range, constitute the relative paleointensity estimates. Some paleointensity lows (particularly those at ~40, ~120 and ~188 ka) are associated with directional excursions of the field, especially the event at ~188 ka (referred to here as the Iceland Basin Event) that constitutes a short-lived polarity reversal. For the last 200 kyr, the records can be correlated with other high-resolution paleointensity records such as those from the Labrador Sea, Mediterranean/Somali Basin and Sulu Sea, implying that the millennial scale features are globally synchronous. A labeling system for paleointensity features is proposed that ties prominent highs and lows to oxygen isotope stages.