Rb, Sr, Sm and Nd concentrations and isotopic composition of basalts from DSDP Leg 74 (Table 1)

Basement intersected in DSDP holes 525A, 528 and 527 on the Walvis Ridge consists of submarine basalt flows and pillows with minor intercalated sediments. These holes are situated on the crest and mid and lower northwest flank of a NNW-SSE-trending ridge block which would have closely paralleled the paleo mid-ocean ridge (Rabinowitz and LaBrecque, 1979 doi:10.1029/JB084iB11p05973, Moore et al. (1983 doi:10.1130/0016-7606(1983)94<907:TWRTDS>2.0.CO;2). The basalts were erupted approximately 70 m.y. ago, an age equivalent to that of immediately adjacent oceanic crust in the Angola Basin and coraistent with formation at the paleo mid-ocean ridge (Moore et al., 1983). The basalt types vary from aphyric quartz tholeiites on the ridge crest to highly plagioclase phyric olivine tholeiites on the ridge flank. These show systematic differences in incompatible trace element and isotopic composition. Many element and isotope ratio pairs form systematic trends with the ridge crest basalts at one end and the highly phyric ridge flank basalts at the other. The low 143Nd/144Nd (0.51238), 206Pb/204Pb (17.54), 207Pb/204Pb (15.47), 208Pb/204Pb (38.14) and high 87Sr/86Sr (0.70512) ratios of the ridge crest basalts suggest derivation from an old Nd/Sm-, Rb/Sr- and Pb/U-enriched mantle source. This isotopic signature is similar to that of alkaline basalts on Tristan da Cunha but offset to significantly lower Nd and Pb isotopic ratios. The isotopic ratio trends may be extrapolated beyond the ridge flank basalts with higher 143Nd/144Nd (0.51270), 206Pb/204Pb (18.32), 207Pb/204Pb (15.52), 208Pb/204Pb (38.77) and lower 87Sr/86Sr (0.70417) ratios in the direction of increasingly Nd/Sm-, Rb/Sr- and Pb/U-depleted source compositions. These isotopic correlations are equally consistent with mixing of depleted and enriched end member melts or partial melting of an inhomogeneous, variably enriched mantle source. However, observed Zr-Ba-Nb-Y interelement relationships are inconsistent with any simple two-component model of magma mixing, as might result from the rise of a lower mantle plume through the upper mantle. Incompatible element and Pb isotopic systematics also preclude extensive involvement of depleted (N-type) MORB material or its mantle sources. In our preferred petrogenetic model the Walvis Ridge basalts were derived by partial melting of mantle similar to an enriched (E-type) MORB source which had become heterogeneous on a small scale due to the introduction of small-volume melts and metasomatic fluids.

Major and trace elements and Nd and Sr isotope geochemistry of basalts at DSDP Leg 74 Holes

Basement intersected in Holes 525A, 528, and 527 on the Walvis Ridge consists of submarine basalt flows and pillows with minor intercalated sediments. These holes are situated on the crest and mid- and lower NW flank of a NNW-SSE-trending ridge block which would have closely paralleled the paleo mid-ocean ridge. The basalts were erupted approximately 70 Ma, a date consistent with formation at the paleo mid-ocean ridge. The basalt types vary from aphyric quartz tholeiites on the Ridge crest to highly Plagioclase phyric olivine tholeiites on the flank. These show systematic differences in incompatible trace element and isotopic composition, and many element and isotope ratio pairs form systematic trends with the Ridge crest basalts at one end and the highly phyric Ridge flank basalts at the other. The low 143Nd/144Nd (0.51238) and high 87Sr/86Sr (0.70512) ratios of the Ridge crest basalts suggest derivation from an old Nd/Sm and Rb/Sr enriched mantle source. This isotopic signature is similar to that of alkaline basalts on Tristan da Cunha but offset by somewhat lower 143Nd/144Nd values. The isotopic ratio trends may be extrapolated beyond the Ridge flank basalts (which have 143Nd/144Nd of 0.51270 and 87Sr/86Sr of 0.70417) in the direction of typical MORB compositions. These isotopic correlations are equally consistent with mixing of depleted and enriched end-member melts or partial melting of an inhomogeneous, variably enriched mantle source. However, observed Zr-Ba-Nb-Y interelement relationships are inconsistent with any simple two-component model of magma mixing or partial melting. They also preclude extensive involvement of depleted (N-type) MORB material or its mantle sources in the petrogenesis of Walvis Ridge basalts.

Age model and Nd isotope ratios of DSDO Hole 43-386 and ODP Hole 210-1276A

Modern thermohaline circulation plays a role in latitudinal heat transport and in deep-ocean ventilation, yet ocean circulation may have functioned differently during past periods of extreme warmth, such as the Cretaceous. The Late Cretaceous (100-65 Ma) was an important period in the evolution of the North Atlantic Ocean, characterized by opening ocean gateways, long-term climatic cooling and the cessation of intermittent periods of anoxia (oceanic anoxic events, OAEs). However, how these phenomena relate to deep-water circulation is unclear. We use a proxy for deep-water mass composition (neodymium isotopes; e-Nd) to show that, at North Atlantic ODP Site 1276, deep waters shifted in the early Campanian (~78-83 Ma) from e-Nd values of ~-7 to values of ~-9, consistent with a change in the style of deep-ocean circulation but >10 Myr after a change in bottom water oxygenation conditions. A similar, but more poorly dated, trend exists in e-Nd data from DSDP Site 386. The Campanian e-Nd transition observed in the North Atlantic records is also seen in the South Atlantic and proto-Indian Ocean, implying a widespread and synchronous change in deep-ocean circulation. Although a unique explanation does not exist for the change at present, we favor an interpretation that invokes Late Cretaceous climatic cooling as a driver for the formation of Southern Component Water, which flowed northward from the Southern Ocean and into the North Atlantic and proto-Indian Oceans.

CaCO3, foraminifer fragmentation, benthic carbon isotopes, coarse lithic counts and radiogenic isotopes (Sr, Nd, Pb) of the Pliocene and earliest Pleistocene terrigenous component (3.3-2.4 Ma) of IODP Site 306-U1313

We present Plio-Pleistocene records of sediment color, %CaCO3, foraminifer fragmentation, benthic carbon isotopes (d13C) and radiogenic isotopes (Sr, Nd, Pb) of the terrigenous component from IODP Site U1313, a reoccupation of benchmark subtropical North Atlantic Ocean DSDP Site 607. We show that (inter)glacial cycles in sediment color and %CaCO3 pre-date major northern hemisphere glaciation and are unambiguously and consistently correlated to benthic oxygen isotopes back to 3.3 million years ago (Ma) and intermittently so probably back to the Miocene/Pliocene boundary. We show these lithological cycles to be driven by enhanced glacial fluxes of terrigenous material (eolian dust), not carbonate dissolution (the classic interpretation). Our radiogenic isotope data indicate a North American source for this dust (~3.3-2.4 Ma) in keeping with the interpreted source of terrestrial plant wax-derived biomarkers deposited at Site U1313. Yet our data indicate a mid latitude provenance regardless of (inter)glacial state, a finding that is inconsistent with the biomarker-inferred importance of glaciogenic mechanisms of dust production and transport. Moreover, we find that the relation between the biomarker and lithogenic components of dust accumulation is distinctly non-linear. Both records show a jump in glacial rates of accumulation from Marine Isotope Stage, MIS, G6 (2.72 Ma) onwards but the amplitude of this signal is about 3-8 times greater for biomarkers than for dust and particularly extreme during MIS 100 (2.52 Ma). We conclude that North America shifted abruptly to a distinctly more arid glacial regime from MIS G6, but major shifts in glacial North American vegetation biomes and regional wind fields (exacerbated by the growth of a large Laurentide Ice Sheet during MIS 100) likely explain amplification of this signal in the biomarker records. Our findings are consistent with wetter-than-modern reconstructions of North American continental climate under the warm high CO2 conditions of the Early Pliocene but contrast with most model predictions for the response of the hydrological cycle to anthropogenic warming over the coming 50 years (poleward expansion of the subtropical dry zones).

Euphausiid respiration model revamped, link to model results

Euphausiids constitute major biomass component in shelf ecosystems and play a fundamental role in the rapid vertical transport of carbon from the ocean surface to the deeper layers during their daily vertical migration (DVM). DVM depth and migration patterns depend on oceanographic conditions with respect to temperature, light and oxygen availability at depth, factors that are highly dependent on season in most marine regions. Changes in the abiotic conditions also shape Euphausiid metabolism including aerobic and anaerobic energy production. Here we introduce a global krill respiration model which includes the effect of latitude (LAT), the day of the year of interest (DoY), and the number of daylight hours on the day of interest (DLh), in addition to the basal variables that determine ectothermal oxygen consumption (temperature, body mass and depth) in the ANN model (Artificial Neural Networks). The newly implemented parameters link space and time in terms of season and photoperiod to krill respiration. The ANN model showed a better fit (r**2=0.780) when DLh and LAT were included, indicating a decrease in respiration with increasing LAT and decreasing DLh. We therefore propose DLh as a potential variable to consider when building physiological models for both hemispheres. We also tested for seasonality the standard respiration rate of the most common species that were investigated until now in a large range of DLh and DoY with Multiple Linear Regression (MLR) or General Additive model (GAM). GAM successfully integrated DLh (r**2= 0.563) and DoY (r**2= 0.572) effects on respiration rates of the Antarctic krill, Euphausia superba, yielding the minimum metabolic activity in mid-June and the maximum at the end of December. Neither the MLR nor the GAM approach worked for the North Pacific krill Euphausia pacifica, and MLR for the North Atlantic krill Meganyctiphanes norvegica remained inconclusive because of insufficient seasonal data coverage. We strongly encourage comparative respiration measurements of worldwide Euphausiid key species at different seasons to improve accuracy in ecosystem modelling.

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.

Strontium and neodymium isotope data for barite separates and pore waters, and isotope and concentration data for foraminifera and fish teeth from DSDP Leg 85 sediments

Strontium and neodymium isotopic data are reported for barite samples chemically separated from Late Miocene to Pliocene sediments from the eastern equatorial Pacific. At a site within a region of very high productivity close to the equator, 87Sr/86Sr ratios in the barite separates are indistinguishable from those of foraminifera and fish teeth from the same samples. However, at two sites north of the productivity maximum barite separates have slightly, but consistently lower (averaging 0.000062) ratios than the coexisting phases, although values still fall within the total range of published values for the contemporaneous seawater strontium isotope curve. We examine possible causes for this offset including recrystallization of the foraminifera, fish teeth or barite, the presence of non-barite contaminants, or incorporation of older, reworked deep-sea barite; the inclusion of a small amount of hydrothermal barite in the sediments seems most consistent with our data, although there are difficulties associated with adequate production and transportation of this phase. Barite is unlikely to replace calcite as a preferred tracer of seawater strontium isotopes in carbonate-rich sediments, but may prove a useful substitute in cases where calcite is rare or strongly affected by diagenesis. In contrast to the case for strontium, neodymium isotopic ratios in the barite separates are far from expected values for contemporary seawater, and appear to be dominated by an (unobserved) eolian component with high neodymium concentration and low 143Nd/144Nd. These results suggest that the true potential of barite as an indicator of paleocean neodymium isotopic ratios and REE patterns will be realized only when a more selective separation procedure is developed.

Stable isotope record and age model for sapropel stratigraphy of ODP Site 160-967

Stable oxygen isotope data from four holes drilled at the Ocean Drilling Program Site 967, which is located on the lower northern slope of the Eratosthenes Seamount, provide a continuous record of Eastern Mediterranean surface-water conditions during the last 3.2 Ma. A high-resolution stratigraphy for the Pliocene-Pleistocene sequence was established by using a combination of astronomical calibration of sedimentary cycles, nannofossil stratigraphy, and stable oxygen isotope fluctuations. Sapropels and color cycles are present throughout the last 3.2 Ma at Site 967, and their ages, as determined by calibration against the precessional component of the astronomical record, are consistent with those estimated for the sapropels of the classical land-based marine sequences of the Punta Piccola, San Nicola, Singa, and Vrica sections (southern Italy). The Site 967 oxygen isotope record shows large amplitude fluctuations mainly caused by variations in surface water salinity throughout the entire period. Spectral analysis shows that fluctuations in the d18O record were predominantly influenced by orbital obliquity and precessional forcing from 3.2 to 1 Ma, and all main orbital frequencies characterize the d18O record for the last million years. The start of sapropel formation at 3.2 Ma indicates a possible link between sapropel formation and the build up of northern hemisphere ice sheets. The dominance of the obliquity cycle in the interval from 3.2-1 Ma further points to the sensitivity of Eastern Mediterranean climate to the fluctuations in the volume of Arctic ice sheets. An intensification of negative isotope anomalies at Site 967, relative to the open ocean, supports a link between high run-off (during warm periods) and sapropel formation. freshwater input would have inhibited deep-water formation, which led to stagnation of deeper waters. Comparison with the land sections also confirms that differential preservation and diagenesis play a key role in sapropel occurrence.

The North Sea autumn spawning Herring (Clupea harengus L.) Spawning Component Abundance Index (SCAI)

The North Sea autumn-spawning herring (Clupea harengus) stock consists of a set of different spawning components. The dynamics of the entire stock have been well characterized, but although time-series of larval abundance indices are available for the individual components, study of the dynamics at the component level has historically been hampered by missing observations and high sampling noise. A simple state-space statistical model is developed that is robust to these problems, gives a good fit to the data, and proves capable of both handling and predicting missing observations well. Furthermore, the sum of the fitted abundance indices across all components proves an excellent proxy for the biomass of the total stock, even though the model utilizes information at the individual-component level. The Orkney-Shetland component appears to have recovered faster from historic depletion events than the other components, whereas the Downs component has been the slowest. These differences give rise to changes in stock composition, which are shown to vary widely within a relatively short time. The modelling framework provides a valuable tool for studying and monitoring the dynamics of the individual components of the North Sea herring stock.

Sea ice thickness distribution (ITD) Model for the Arctic, links to NetCDF files

Sea ice models contain many different parameterizations of which one of the most commonly used is a subgrid-scale ice thickness distribution (ITD). The effect of this model component and the associated ice strength formulation on the reproduction of observed Arctic sea ice is assessed. To this end the model's performance in reproducing satellite observations of sea ice concentration, thickness and drift is evaluated. For an unbiased comparison, different model configurations with and without an ITD are tuned with an automated parameter optimization. The original combination of ITD and ice strength parameterization does not lead to better results than a simple single category model. Yet changing to a simpler ice strength formulation, which depends linearly on the mean ice thickness across all thickness categories, allows to clearly improve the model-data misfit when using an ITD. In the original formulation, the ice strength depends strongly on the number of thickness categories, so that introducing more categories can lead to thicker albeit weaker ice on average.

(Table 1) Age comparison between the GPTS and ODP Site 181-1123 age model for each chron base within the Middle Miocene section

Today the deep western boundary current (DWBC) east of New Zealand is the most important route for deep water entering the Pacific Ocean. Large-scale changes in deep water circulation patterns are thought to have been associated with the development of the East Antarctic Ice Sheet (EAIS) close to the main source of bottom water for the DWBC. Here we reconstruct the changing speed of the southwest Pacific DWBC during the middle Miocene from ~15.5-12.5 Ma, a period of significant global ice accumulation associated with EAIS growth. Sortable silt mean grain sizes from Ocean Drilling Program Site 1123 reveal variability in the speed of the Pacific inflow on the timescale of the 41 kyr orbital obliquity cycle. Similar orbital period flow changes have recently been demonstrated for the Pleistocene epoch. Collectively, these observations suggest that a strong coupling between changes in the speed of the deep Pacific inflow and high-latitude climate forcing may have been a persistent feature of the global thermohaline circulation system for at least the past 15 Myr. Furthermore, long-term changes in flow speed suggest an intensification of the DWBC under an inferred increase in Southern Component Water production. This occurred at the same time as decreasing Tethyan outflow and major EAIS growth between ~15.5 and 13.5 Ma. These results provide evidence that a major component of the deep thermohaline circulation was associated with the middle Miocene growth of the EAIS and support the view that this time interval represents an important step in the development of the Neogene icehouse climate.

(Table S1) Tie-points used in age model construction for ODP Leg 172 sites

ix Ocean Drilling Program (ODP) sites, in the Northwest Atlantic have been used to investigate kinematic and chemical changes in the "Western Boundary Undercurrent" (WBUC) during the development of full glacial conditions across the Marine Isotope Stage 5a/4 boundary (~70,000 years ago). Sortable silt mean grain size(sort s) measurements are employed to examine changes in near bottom flow speeds, together with carbon isotopes measured in benthic foraminifera and % planktic foraminiferal fragmentation as proxies for changes in water-mass chemistry. A depth transect of cores, spanning 1.8-4.6 km depth, allows changes in both the strength and depth of the WBUC to be constrained across millennial scale events. Sort s measurements reveal that the flow speed structure of the WBUC during warm intervals ("interstadials") was comparable to modern (Holocene) conditions. However, significant differences are observed during cold intervals, with higher relative flow speeds inferred for the shallow component of the WBUC (~2 km depth) during all cold "stadial" intervals (including Heinrich Stadial 6), and a substantial weakening of the deep component (~3-4 km) during full glacial conditions. Our results therefore reveal that the onset of full glacial conditions was associated with a regime shift to a shallower mode of circulation (involving Glacial North Atlantic Intermediate Water) that was quantitatively distinct from preceding cold stadial events. Furthermore, our chemical proxy data show that the physical response of the WBUC during the last glacial inception was probably coupled to basin-wide changes in the water-mass composition of the deep Northwest Atlantic.