Modelling Nd-isotopes with a coarse resolution ocean circulation model: Sensitivities to model parameters and source/sink distributions

The neodymium (Nd) isotopic composition (Nd) of seawater is a quasi-conservative tracer of water mass mixing and is assumed to hold great potential for paleoceanographic studies. Here we present a comprehensive approach for the simulation of the two neodymium isotopes 143Nd, and 144Nd using the Bern3D model, a low resolution ocean model. The high computational efficiency of the Bern3D model in conjunction with our comprehensive approach allows us to systematically and extensively explore the sensitivity of Nd concentrations and Nd to the parametrisation of sources and sinks. Previous studies have been restricted in doing so either by the chosen approach or by computational costs. Our study thus presents the most comprehensive survey of the marine Nd cycle to date. Our model simulates both Nd concentrations as well as Nd in good agreement with observations. Nd covaries with salinity, thus underlining its potential as a water mass proxy. Results confirm that the continental margins are required as a Nd source to simulate Nd concentrations and Nd consistent with observations. We estimate this source to be slightly smaller than reported in previous studies and find that above a certain magnitude its magnitude affects Nd only to a small extent. On the other hand, the parametrisation of the reversible scavenging considerably affects the ability of the model to simulate both, Nd concentrations and Nd. Furthermore, despite their small contribution, we find dust and rivers to be important components of the Nd cycle. In additional experiments, we systematically varied the diapycnal diffusivity as well as the Atlantic-to-Pacific freshwater flux to explore the sensitivity of Nd concentrations and its isotopic signature to the strength and geometry of the overturning circulation. These experiments reveal that Nd concentrations and Nd are comparatively little affected by variations in diapycnal diffusivity and the Atlantic-to-Pacific freshwater flux. In contrast, an adequate representation of Nd sources and sinks is crucial to simulate Nd concentrations and Nd consistent with observations. The good agreement of our results with observations paves the way for the evaluation of the paleoceanographic potential of Nd in further model studies.

Major and trace element analyses of ODP Leg 120 volcanic rocks (Table 1)

Widespread Lower Cretaceous magmatism occurred along the Indian-Australian/Antarctic margins, and in the juvenile Indian Ocean, during the rifting of eastern Gondwana. The formation of this magmatic province probably began around 120-130 Ma with the eruption of basalts on the Naturaliste Plateau and at Bunbury, western Australia. On the northeast margin of India, activity began around 117 Ma with the Rajmahal continental basalts and associated lamprophyre intrusions. The formation of the Kerguelen Plateau in the Indian Ocean began no later than 114 Ma. Ultramafic lamprophyres (alnoites) were emplaced in the Prince Charles Mountains near the Antarctic continental margin at ~ 110 Ma. These events are considered to be related to a major mantle plume, the remnant of which is situated beneath the region of Kerguelen and Heard islands at the present day. Geochemical data are presented for each of these volcanic suites and are indicative of complex interactions between asthenosphere-derived magmas and the continental lithosphere. Kerguelen Plateau basalts have Sr and Nd isotopic compositions lying outside the field for Indian Ocean mid-ocean ridge basalts (MORB) but, with the exception of Site 738 at the southern end of the plateau, within the range of more recent hotspot basalts from Kerguelen and Heard Islands. However, a number of the plateau tholeiites are characterized by lower 206Pb/204Pb ratios than are basalts from Kerguelen Island, and many also have anomalously high La/Nb ratios. These features suggest that the source of the Kerguelen Plateau basalts suffered contamination by components derived from the Gondwana continental lithosphere. An extreme expression of this lithospheric signature is shown by a tholeiite from Site 738, suggesting that the southernmost part of the Kerguelen Plateau may be underlain by continental crust. The Rajmahal tholeiites mostly fall into two distinct geochemical groups. Some Group I tholeiites have Sr and Nd isotopic compositions and incompatible element abundances, similar to Kerguelen Plateau tholeiites from Sites 749 and 750, indicating that the Kerguelen-Heard mantle plume may have directly furnished Rajmahal volcanism. However, their elevated 207Pb/204Pb ratios indicate that these magmas did not totally escape contamination by continental lithosphere. In contrast to the Group I tholeiites, significant contamination is suggested for Group II Rajmahal tholeiites, on the basis of incompatible element abundances and isotopic compositions. The Naturaliste Plateau and the Bunbury Basalt samples show varying degrees of enrichment in incompatible elements over normal MORB. The Naturaliste Plateau samples (and Bunbury Basalt) have high La/Nb ratios, a feature not inconsistent with the notion that the plateau may consist of stretched continental lithosphere, near the ocean-continent divide.

Major and trace element compositions of the Buendnerschiefer samples

The Bündnerschiefer of the Swiss-Italian Alps is a large sedimentary complex deposited on the Piemonte-Liguria and Valais oceans and associated continental margins from the upper Jurassic to Eocene. It is made of a large variety of sequences associated or not with an ophiolitic basement. The Bündnerschiefer makes an accretionary prism that developed syn-tectonically from the onset of alpine subduction, and it records orogenic metamorphism following episodes of HP metamorphism. The Bündnerschiefer shares important similarities with the Otago schists of New Zealand and with the Wepawaug schists of Connecticut, both of which form accretionary prisms and have an orogenic metamorphic imprint. With the aim of testing the hypothesis of mobility of chemical components as a function of metamorphic grade, in this work I present fifty-five bulk chemical analyses of various lithological facies of the Bündnerschiefer collected along the well-studied field gradient of the Lepontine dome of Central Switzerland, in the Prättigau half window of East Switzerland, and in the Tsaté Nappe of Valle d'Aosta (Italy). The dataset includes the concentration of major components, large ion lithophile elements (Rb, Sr, Ba, Cs), high field strength elements (Zr, Ti, Nb, Th, U, Ta, Hf), fluid-mobile light elements (B, Li), volatiles (CO2, S), REEs, and Y, V, Cr, Co, Sn, Pb, Cu, Zn, Tl, Sb, Be, and Au. These data are compared against the compositions of the global marine sediment reservoir, typical crustal reservoirs, and against the previously measured compositions of Otago and Wepawaug schists. Results reveal that, irrespective of their metamorphic evolution, the bulk chemical compositions of orogenic metasediments are characterized by mostly constant compositional ratios (e.g., K2O/Al2O3, Ba/Al2O3, Sr/CaO, etc.), whose values in most cases are undistinguishable from those of actual marine sediments and other crustal reservoirs. For these rocks, only volatile concentrations decrease dramatically as a function of metamorphic temperature, and significant deviations from the reservoir signatures are evident for SiO2, B, and Li. These results are interpreted as an indication of residual enrichment in the sediments, a process taking place during syn-metamorphic dehydration from the onset of metamorphism in a regime of chemical immobility. Residual enrichment increased the absolute concentrations of the chemical components of these rocks, but did not modify significantly their fundamental ratios. This poor compositional modification of the sediments indicates that orogenic metamorphism in general does not promote significant mass transfer from accretionary prisms. In contrast, mass transfer calculations carried out in a shear zone crosscutting the Bündnerschiefer shows that significant mass transfer occurs within these narrow zones, resulting in gains of H2O, SiO2, Al2O3, K2O, Ba, Y, Rb, Cu, V, Tl, Mo, and Ce during deformation and loss of Na2O, CO2, S, Ni, B, U, and Pb from the rock. These components were presumably transported by an aquo-carbonic fluid along the shear zone. These distinct attitudes to mobilize chemical elements from orogenic sediments may have implications for a potentially large number of geochemical processes in active continental margins, from the recycling of chemical components at plate margins to the genesis of hydrothermal ore deposits.

Radionuclides measured on water bottle samples during SONNE cruise SO202 (INOPEX)

Boundary scavenging, or the enhanced removal of adsorption-prone elements from the ocean in areas of high particle flux, is an often cited, though not well-quantified, concept used to understand the oceanic distribution of many trace metals. Because 230Th and 231Pa are produced uniformly from uranium decay and removed differentially by scavenging, the process of boundary scavenging can be elucidated by a more detailed knowledge of their water column distributions. To this end, filtered seawater was collected across the gradients in particle flux which span the subarctic Pacific: in the west during the Innovative North Pacific Experiment (INOPEX) and in the east along Line P. Lateral concentration gradients of dissolved 230Th are small throughout the subarctic Pacific at 12 sites of variable particle flux. This contradicts the prediction of the traditional boundary scavenging model. A compilation of water column data from throughout the North Pacific reveals much larger lateral concentration gradients for 230Th between the subarctic North Pacific and subtropical gyre, over lateral gradients in scavenging intensity similar to those found within the subarctic. This reflects a biogeochemical-province aspect to scavenging. Upper water column distributions of 231Pa and 231Pa/230Th ratio are consistent with the influence of scavenging by biogenic opal, while deep waters (>2.5 km) reveal an additional 231Pa sink possibly related to manganese oxides produced at continental margins or ridge crests.

Age models and stable isotope ratios from the Mid-Miocene

Mid-Miocene pelagic sedimentary sections can be correlated using intermediate and high resolution oxygen and carbon isotopic records of benthic foraminifera. Precision of a few tens of thousands of years is readily achievable at sites with high sedimentation rates, for example, Deep Sea Drilling Project sites 289 and 574. The mid-Miocene carbon isotope records are characterized by an interval of high d13C values between 17 and 13.5 Ma (the Monterey Excursion of Vincent and Berger 1985) upon which are superimposed a series of periodic or quasi-periodic fluctuations in d13C values. These fluctuations have a period of approximately 440 kyr, suggestive of the 413 kyr cycle predicted by Milankovitch theory. Vincent and Berger proposed that the Monterey Excursion was the result of increased organic carbon burial in continental margins sediments. The increased d13C values (called 13C maxima) superimposed on the generally high mid-Miocene signal coincide with increases in d18O values suggesting that periods of cooling and/or ice buildup were associated with exceptionally rapid burial of organic carbon and lowered atmospheric CO2 levels. It is likely that during the Monterey Excursion the ocean/atmosphere system became progressively more sensitive to small changes in insolation, ultimately leading to major cooling of deep water and expansion of continental ice. We have assigned an absolute chronology, based on biostratigraphic and magneto-biostratigraphic datum levels, to the isotope stratigraphy and have used that chronology to correlate unconformities, seismic reflectors, carbonate minima, and dissolution intervals. Intervals of sediment containing 13C maxima are usually better preserved than the overlying and underlying sediments, indicating that the d13C values of TCO2 in deep water and the corrosiveness of seawater are inversely correlated. This again suggests that the 13C maxima were associated with rapid burial of organic carbon and reduced levels of atmospheric CO2. The absolute chronology we have assigned to the isotopic record indicates that the major mid-Miocene deepwater cooling/ice volume expansion took 2 m.y. and was not abrupt as had been reported previously. The cooling appears abrupt at many sites because the interval is characterized by a number of dissolution intervals. The cooling was not monotonic, and the 2 m.y. interval included an episode of especially rapid cooling as well as a brief return to warmer conditions before the final phase of the cooling period. The increase in d18O values of benthic foraminifera between 14.9 and 12.9 Ma was greatest at deeper water sites and at sites closest to Antarctica. The data suggest that the d18O value of seawater increased by no more than about 1.1 per mil during this interval and that the remainder of the change in benthic d18O values resulted from cooling in Antarctic regions of deepwater formation. Equatorial planktonic foraminifera from sites 237 and 289 exhibit a series of 0.4 per mil steplike increases in d13C values. Only one of these increases in planktonic d13C is correlated with any of the features in the mid-Miocene benthic carbon isotope record.

Geochemistry at DSDP Legs 36 and 71 Holes

For Middle Jurassic to Pleistocene times, clay mineralogical and geochemical data provide information on the evolution of continental and marine paleoenvironments. They are a source of information on marginal instability, on the continental and shallow marine environments related to the development of the Southern Ocean during the Middle and Late Jurassic, and on tectonic relaxation of the continental margins at the end of the Late Jurassic. They also provide evidence for the influences of the South Atlantic opening and the movement of the Falkland Plateau in a reduced marine environment until Aptian-Albian times, and the transition to an open marine environment during Albian time; the influences of the Albian-Turonian and Coniacian-Santonian Andean deformations in an open marine environment; the limited tectonic effects and strong influence of marine currents at the Cretaceous/Tertiary boundary; the influences of the global climatic cooling and inferred bottom water circulation during the late Eocene and Oligocene; the widening of the South Atlantic Ocean during Oligocene time, which was accompanied by an increased influence of the biogenic components on sedimentation; increased carbonate dissolution from late Oligocene to early Miocene, related to the deepening of the ocean; limited mineralogical and important geochemical modifications when the Drake Passage opened in the early Miocene; the influence of the late Miocene development of the Antarctic ice-sheet; the major Antarctic cooling and Patagonian glaciation during Pliocene time; and the change in the Antarctic Bottom Water circulation at the Pliocene/Pleistocene boundary.

Uk'37, derived sea surface temperatures and coccolith abundances in North Atlantic sediment cores

The UK37' index has proven to be a robust proxy to estimate past sea surface temperatures (SSTs) over a range of time scales, but like any other proxy, it has uncertainties. For instance, in reconstructions of the Last Glacial Maximum (LGM) in the northern North Atlantic, UK37' indicates higher temperatures than those derived from foraminiferal proxies. Here we evaluate whether such warm glacial estimates are caused by the advection of reworked alkenones in ice-rafted debris (IRD) to deep-sea sediments. We have quantified both coccolith assemblages and alkenones in sediments from glaciogenic debris flows in the continental margins of the northern North Atlantic, and from a deep-sea core from the Reykjanes Ridge. Certain debris flow deposits in the North Atlantic were generated by the presence of massive ice-sheets in the past, and their associated ice streams. Such deposits are composed of the same materials that were present in the IRD at the time they were generated. We conclude that ice rafting from some locations was a transport pathway to the deep sea floor of reworked alkenones and pre-Quaternary coccolith species during glacial stages, but that not all of the IRD contained alkenones, even when reworked coccoliths were present. We speculate that the ratio of reworked coccoliths to alkenone concentration might be useful to infer whether significant reworked alkenone inputs from IRD did occur at a particular site in the glacial North Atlantic. We also observe that alkenones in some of the debris flows contain a colder signal than estimated for LGM sediments in the northern North Atlantic. This is also clear in the deep-sea core studied where the warmest intervals do not correspond to the intervals with large inputs of reworked coccoliths or IRD. We conclude that any possible bias to UK37' estimates associated with reworked alkenones is not necessarily towards higher values, and that the high SST anomalies for the LGM are unlikely to be the result of a bias caused by IRD inputs.

(Table S1) Molecular analyses and morphological comparisons of deep-sea clams

Large vesicomyid clams are common inhabitants of sulphidic deep-sea habitats such as hydrothermal vents, hydrocarbon seeps and whale-falls. Yet, the species- and genus-level taxonomy of these diverse clams has been unstable due to insufficiencies in sampling and absence of detailed taxonomic studies that would consistently compare molecular and morphological characters. To clarify uncertainties about species-level assignments, we examined DNA sequences from mitochondrial cytochrome-c-oxidase subunit I (COI) in conjunction with morphological characters. New and published COI sequences were used to create a molecular database for 44 unique evolutionary lineages corresponding to species. Overall, the congruence between molecular and morphological characters was good. Several discrepancies due to synonymous species designations were recognized, and acceptable species names were rectified with published COI sequences in cases where morphological specimens were available. We identified seven species with trans-Pacific distributions, and two species with Indo-Pacific distributions. Presently, 27 species have only been documented from one region, which might reflect limited ranges, or insufficient geographical sampling. Vesicomyids exhibit the greatest species diversity along the northwest Pacific ridge systems and in the eastern Pacific, along the western America margin, where depth zonation typically results in segregation of closely related species. The broad distributions of several vesicomyid species suggest that their required chemosynthetic habitats might be more common than previously recognized and occur along most continental margins.