Deep water formation in the North Pacific and deglacial CO2 rise

Deep water formation in the North Atlantic and Southern Ocean is widely thought to influence deglacial CO2 rise and climate change; here we suggest that deep water formation in the North Pacific may also play an important role. We present paired radiocarbon and boron isotope data from foraminifera from sediment core MD02-2489 at 3640 m in the North East Pacific. These show a pronounced excursion during Heinrich Stadial 1, with benthic-planktic radiocarbon offsets dropping to ~350 years, accompanied by a decrease in benthic d11B. We suggest this is driven by the onset of deep convection in the North Pacific, which mixes young shallow waters to depth, old deep waters to the surface, and low-pH water from intermediate depths into the deep ocean. This deep water formation event was likely driven by an increase in surface salinity, due to subdued atmospheric/monsoonal freshwater flux during Heinrich Stadial 1. The ability of North Pacific Deep Water (NPDW) formation to explain the excursions seen in our data is demonstrated in a series of experiments with an intermediate complexity Earth system model. These experiments also show that breakdown of stratification in the North Pacific leads to a rapid ~30 ppm increase in atmospheric CO2, along with decreases in atmospheric d13C and D14C, consistent with observations of the early deglaciation. Our inference of deep water formation is based mainly on results from a single sediment core, and our boron isotope data are unavoidably sparse in the key HS1 interval, so this hypothesis merits further testing. However we note that there is independent support for breakdown of stratification in shallower waters during this period, including a minimum in d15N, younging in intermediate water 14C, and regional warming. We also re-evaluate deglacial changes in North Pacific productivity and carbonate preservation in light of our new data, and suggest that the regional pulse of export production observed during the Bølling-Allerød is promoted by relatively stratified conditions, with increased light availability and a shallow, potent nutricline. Overall, our work highlights the potential of NPDW formation to play a significant and hitherto unrealized role in deglacial climate change and CO2 rise.

Geochemistry of ODP Leg 127 igneous rocks

Legs 127 and 128 of the Ocean Drilling Program cored basement samples from two sites in the Yamato Basin (Sites 794 and 797) and one site in the Japan Basin (Site 795) of the Japan Sea. These samples represent sills and lava flows erupted or shallowly intruded in a marine environment during backarc extension and spreading in the middle Miocene. In this paper, we describe the geochemical characteristics of these igneous units using 52 new instrumental neutron activation analyses (INAA), 8 new X-ray fluorescence (XRF) analyses, and previous shipboard XRF analyses. The sills intruded into soft sediment at Sites 794 and 797 were subject to extensive hydrothermal activity, estimated at <230° C under subgreenschist facies conditions, which heavily to totally altered the fine-grained unit margins and moderately to heavily altered the coarse-grained unit interiors. Diagenesis further altered the composition of these igneous bodies and lava flows at Sites 794, 795, and 797, most intensely at unit margins. Our study of two well-sampled units shows that Mg, Ca, Sr, and the large-ion lithophile elements (LILE) mobilized during alteration, and that the concentrations of Y, Yb, and Lu decreased and Ce increased in the most severely altered samples. Nevertheless, our study shows that the rare-earth elements (REE) were relatively immobile in the majority of the samples, even where secondary mixed-layer clays comprised the great majority of the rock. Fresher Yamato Basin samples are compositionally heterogenous tholeiitic basalts and dolerites. At Site 794 in the north-central portion of the basin, Units 1 to 5 (upper basement) comprise mildly light rare-earth element (LREE) enriched basalts and dolerites (chondrite-normalized La/Sm of 1.4-1.8), while the stratigraphically lower Units 6 to 9 are less enriched dolerites with (La/Sm)N of 0.7-1.3. All Site 794 samples lack Nb and Ta depletions and LILE enrichments, lacking a strong subduction-related incompatible element geochemical signature. At Site 797 in the western margin of the basin, two stratigraphically-definable unit groups also occur. The upper nine units are incompatible-element depleted tholeiitic sills and flows with strong depletions of Nb and Ta relative to normal mid-ocean ridge basalt (N-MORB). The lower twelve sills represent LREE-enriched tholeiites (normalized La/Sm ranges from 1.1 to 1.8), with distinctly higher LILE and high field-strength element (HFSE) contents. At Site 795 at the northern margin of the Japan Sea, three eruptive units consist of basaltic andesite to calc-alkaline basalt (normalized La/Sm of 1.1 to 1.5) containing moderate depletions of the HFSE relative to N-MORB. The LILE-depleted nature of these samples precludes their origin in a continental arc, indicating that they more likely erupted within a rifting oceanic arc system. The heterogenous nature of the Japan Sea rocks indicate that they were derived at each site from multiple parental magmas generated from a compositionally heterogenous mantle source. Their chemistry is intermediate in character between arc basalts, MORB, and intraplate basalts, and implies little involvement of continental crust at any point in their genesis. Their flat chondrite-normalized, medium-to-heavy rare earth patterns indicate that the primary magmas which produced them last equilibrated with and segregated from spinel lherzolite at shallow depths (<30 kbar). In strong contrast to their isotopic compositional arrays, subduction-related geochemical signatures are usually poorly defined. No basin-wide temporal or geographic systematics of rock chemistry may be confidently detailed; instead, the data show both intimate (site-specific) and widespread backarc mantle heterogeneity over a narrow (2 Ma or so) range in time, with mantle heterogeneity most closely resembling a "plum-pudding" model.

Radiolarian fauna of sediment core MD96-2089 off Lüderitz, Namibia

The study of radiolarian assemblages from Core MD 962086 provides new information on the variability in the upwelling intensity and origin of upwelled water masses over the past 350 ky in one of the major filamentous regions of the Benguela Upwelling System (BUS), located off Lüderitz, Namibia. The use of key radiolarian species to trace the source of upwelled waters, and the use of a radiolarian-based upwelling index (URI) to reconstruct the upwelling intensity represent the first use of radiolarians for paleoceanographic reconstructions in the BUS. These radiolarian-based proxies indicate strongest upwelling during Marine Isotope Stages (MIS) 3, 5, and 8, which compares well with other studies. While during MIS 3 and 8, the radiolarian-based proxies indicate the influx of waters of Southern Ocean origin, they also point to the increased influence of tropical waters during the lower portion of MIS 5. During MIS 2, 4 and 6 the radiolarian assemblages indicate generally lower upwelling intensities, although this signal is complicated by the increased occurrence of organic carbon in the sediments during these intervals. During MIS 2 there appears to be less of an input of Southern Ocean waters to the BUS, although during the also glacial MIS 4 and 6, there is evidence for an increased influence of cold Antarctic waters. The comparison of the results from Core MD 962086 with other studies in the BUS area indicates a non-uniform pattern of upwelling intensity and advection of cold, southern waters into this system during MIS 2. Weaker upwelling signaled by the radiolarian-based proxy in MIS 4 is in contrast to other studies that indicate higher productivity during this time period. In general, the data show that there is a strong spatiotemporal complexity in upwelling intensity in the BUS and that the advection of water into it is not strongly tied to glacial-interglacial variations in climate.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Major, trace and rare earth element concentrations of sediments from ODP Leg 207 sites

The Cretaceous and Paleogene sediments recovered during Ocean Drilling Program Leg 207 can be divided into three broad modes of deposition: synrift clastics (lithologic Unit V), organic matter-rich, laminated black shales (Unit IV), and open-marine chalk and calcareous claystones (Units III-I). The aim of this study is to provide a quantitative geochemical characterization of sediments representing these five lithologic units. For this work we used the residues (squeeze cakes) obtained from pore water sampling. Samples were analyzed for bulk parameters (total inorganic carbon, total organic carbon, and S) and by X-ray fluorescence for major (Si, Ti, Al, Fe, Mn, Mg, Ca, Na, K, and P) and selected minor (As, Ba, Co, Cr, Cu, Mo, Ni, Pb, Rb, Sr, U, V, Y, Zn, and Zr) elements. Inductively coupled plasma-mass spectrometry analyses for rare earth elements (REEs) were performed on acid digestions of the squeeze cake samples from Site 1258. The major element composition is governed by the mixture of a terrigenous detrital component of roughly average shale (AS) composition with biogenous carbonate and silica. The composition of the terrigenous detritus is close to AS in Units II-IV. For Unit I, a more weathered terrigenous source is suggested. Carbonate contents reach >60 wt% on average in chalks and calcareous claystones of Units II-IV. The SiO2 contribution in excess of the normal terrigenous-detrital background indicates the presence of biogenous silica, with highest amounts in Units II and III. The contents of coarse-grained material (quartz) are enhanced in Unit V, where Ti and Zr contents are also high. This indicates a high-energy depositional environment. REE patterns are generally similar to AS. A more pronounced negative Ce anomaly in Unit IV may indicate low-oxygen conditions in the water column. The Cretaceous black shales of Unit IV are clearly enriched in redox-sensitive and stable sulfide-forming elements (Mo, V, Zn, and As). High phosphate contents point toward enhanced nutrient supply and high bioproductivity. Ba/Al ratios are rather high throughout Unit IV despite the absence of sulfate in the pore water, indicating elevated primary production. Manganese contents are extremely low for most of the interval studied. Such an Mn depletion is only possible in an environment where Mn was mobilized and transported into an expanded oxygen minimum zone ("open system"). The sulfur contents show a complete sulfidation of the reactive iron of Unit IV and a significant excess of sulfur relative to that of iron, which indicates that part of the sulfur was incorporated into organic matter. We suppose extreme paleoenvironmental conditions during black shale deposition: high bioproductivity like in recent coastal upwelling settings together with severe oxygen depletion if not presence of hydrogen sulfide in the water column.

Whole-rock and mineral geochemistry of ODP Hole 176-735B gabbros

We report mineral chemistry, whole-rock major element compositions, and trace element analyses on Hole 735B samples drilled and selected during Leg 176. We discuss these data, together with Leg 176 shipboard data and Leg 118 sample data from the literature, in terms of primary igneous petrogenesis. Despite mineral compositional variation in a given sample, major constituent minerals in Hole 735B gabbroic rocks display good chemical equilibrium as shown by significant correlations among Mg# (= Mg/[Mg + Fe2+]) of olivine, clinopyroxene, and orthopyroxene and An (=Ca/[Ca + Na]) of plagioclase. This indicates that the mineral assemblages olivine + plagioclase in troctolite, plagioclase + clinopyroxene in gabbro, plagioclases + clinopyroxene + olivine in olivine gabbro, and plagioclase + clinopyroxene + olivine + orthopyroxene in gabbronorite, and so on, have all coprecipitated from their respective parental melts. Fe-Ti oxides (ilmenite and titanomagnetite), which are ubiquitous in most of these rocks, are not in chemical equilibrium with olivine, clinopyroxene, and plagioclase, but precipitated later at lower temperatures. Disseminated oxides in some samples may have precipitated from trapped Fe-Ti-rich melts. Oxides that concentrate along shear bands/zones may mark zones of melt coalescence/transport expelled from the cumulate sequence as a result of compaction or filter pressing. Bulk Hole 735B is of cumulate composition. The most primitive olivine, with Fo = 0.842, in Hole 735B suggests that the most primitive melt parental to Hole 735B lithologies must have Mg# 0.637, which is significantly less than Mg# = 0.714 of bulk Hole 735B. This suggests that a significant mass fraction of more evolved products is needed to balance the high Mg# of the bulk hole. Calculations show that 25%-45% of average Eastern Atlantis II Fracture Zone basalt is needed to combine with 55%-75% of bulk Hole 735B rocks to give a melt of Mg# 0.637, parental to the most primitive Hole 735B cumulate. On the other hand, the parental melt with Mg# 0.637 is far too evolved to be in equilibrium with residual mantle olivine of Fo > 0.89. Therefore, a significant mass fraction of more primitive cumulate (e.g., high Mg# dunite and troctolite) is yet to be sampled. This hidden cumulate could well be deep in the lower crust or simply in the mantle section. We favor the latter because of the thickened cold thermal boundary layer atop the mantle beneath slow-spreading ridges, where cooling and crystallization of ascending mantle melts is inevitable. These observations and data interpretation require reconsideration of the popular concept of primary mantle melts and relationships among the extent of mantle melting, melt production, and the composition and thickness of igneous crust.

Geochemistry of Central American Volcanic Arc basement samples

The Central American Volcanic Arc (CAVA) has been the subject of intensive research over the past few years, leading to a variety of distinct models for the origin of CAVA lavas with various source components. We present a new model for the NW Central American Volcanic Arc based on a comprehensive new geochemical data set (major and trace element and Sr-Nd-Pb-Hf-O isotope ratios) of mafic volcanic front (VF), behind the volcanic front (BVF) and back-arc (BA) lava and tephra samples from NW Nicaragua, Honduras, El Salvador and Guatemala. Additionally we present data on subducting Cocos Plate sediments (from DSDP Leg 67 Sites 495 and 499) and igneous oceanic crust (from DSDP Leg 67 Site 495), and Guatemalan (Chortis Block) granitic and metamorphic continental basement. We observe systematic variations in trace element and isotopic compositions both along and across the arc. The data require at least three different endmembers for the volcanism in NW Central America. (1) The NW Nicaragua VF lavas require an endmember with very high Ba/(La, Th) and U/Th, relatively radiogenic Sr, Nd and Hf but unradiogenic Pb and low d18O, reflecting a largely serpentinite-derived fluid/hydrous melt flux from the subducting slab into a depleted N-MORB type of mantle wedge. (2) The Guatemala VF and BVF mafic lavas require an enriched endmember with low Ba/(La, Th), U/Th, high d18O and radiogenic Sr and Pb but unradiogenic Nd and Hf isotope ratios. Correlations of Hf with both Nd and Pb isotopic compositions are not consistent with this endmember being subducted sediments. Granitic samples from the Chiquimula Plutonic Complex in Guatemala have the appropriate isotopic composition to serve as this endmember, but the large amounts of assimilation required to explain the isotope data are not consistent with the basaltic compositions of the volcanic rocks. In addition, mixing regressions on Nd vs. Hf and the Sr and O isotope plots do not go through the data. Therefore, we propose that this endmember could represent pyroxenites in the lithosphere (mantle and possibly lower crust), derived from parental magmas for the plutonic rocks. (3) The Honduras and Caribbean BA lavas define an isotopically depleted endmember (with unradiogenic Sr but radiogenic Nd, Hf and Pb isotope ratios), having OIB-like major and trace element compositions (e.g. low Ba/(La, Th) and U/Th, high La/Yb). This endmember is possibly derived from melting of young, recycled oceanic crust in the asthenosphere upwelling in the back-arc. Mixing between these three endmember types of magmas can explain the observed systematic geochemical variations along and across the NW Central American Arc.

Geochemical and mineralogical analyses of clay-rich intervals in ODP Sites 166-1006 and 166-1007

The western flank of the Great Bahama Bank, drilled during ODP Leg 166 at seven sites, represents a prograding carbonate sequence from late Oligocene to Holocene [Eberli et al., Proc. ODP Init. Reports 166 (1997)]. The signatures of the detrital input and of diagenetic alteration are evident in clay enriched intervals from the most distal Sites 1006 and 1007 in the Straits of Florida. Mineralogical and chemical investigations (XRD, TEM, SEM, ICP-MS) run on bulk rocks and on the clay fractions enable the origin and evolution of silicate parageneses to be characterized. Plio-Pleistocene silt and clay interbeds contain detrital clay assemblages comprising chlorite, illite, interstratified illite smectite, smectite, kaolinite and palygorskite. The greater smectite input within late Pliocene units than in Pleistocene oozes may relate either varying source areas or change in paleoclimatic conditions and weathering intensity. The clay intervals from Miocene-upper Oligocene wackestone sections are fairly different, with prevalent smectite in the fine fraction, whose high crystallinity and Mg contents that point towards an authigenic origin. The lower Miocene section, below 1104 mbsf, at depths where compaction features are well developed, is particularly characterized by abundant authigenic Na-K-clinoptilolite filling foraminifer tests. The authigenic smectite and clinoptilolite paragenesis is recorded by the chemical trends, both of the sediment and the interstitial fluid. This diagenetic evolution implies Si- and Mg rich fluids circulating in deeper and older sequences. For lack of any local volcaniclastic input, the genesis of zeolite and the terms of water rock interaction are discussed. The location of the diagenetic front correlates with that of the seismic sequence boundary P2 dated as 23.2 Ma. This correspondence may allow the chronostratigraphic significance of some specific seismic reflections to be reassessed.

Geochemistry of abyssal peridotites from ODP Hole 209-1274A

Aqueous dihydrogen (H2,aq) is produced in copious amounts when seawater interacts with peridotite and H2O oxidizes ferrous iron in olivine to ferric iron in secondary magnetite and serpentine. Poorly understood in this process is the partitioning of iron and its oxidation state in serpentine, although both impose an important control on dihydrogen production. We present results of detailed petrographic, mineral chemical, magnetic and Mößbauer analyses of partially to fully serpentinized peridotites from the Ocean Drilling Program (ODP) Leg 209, Mid-Atlantic Ridge (MAR) 15°N area. These results are used to constrain the fate of iron during serpentinization and are compared with phase equilibria considerations and peridotite-seawater reaction path models. In samples from Hole 1274A, mesh-rims reveal a distinct in-to-out zoning from brucite at the interface with primary olivine, followed by a zone of serpentine + brucite ± magnetite and finally serpentine + magnetite in the outermost mesh-rim. The compositions of coexisting serpentine (Mg# 95) and brucite (Mg# 80) vary little throughout the core. About 30-50% of the iron in serpentine/brucite mesh-rims is trivalent, irrespective of subbasement depth and protolith (harzburgite versus dunite). Model calculations suggest that both partitioning and oxidation state of iron are very sensitive to temperature and water-to-rock ratio during serpentinization. At temperatures above 330 °C the dissolution of olivine and coeval formation of serpentine, magnetite and dihydrogen depends on the availability of an external silica source. At these temperatures the extent of olivine serpentinization is insufficient to produce much hydrogen, hence conditions are not reducing enough to form awaruite. At T < 330 °C, hydrogen generation is facilitated by the formation of brucite, as dissolution of olivine to form serpentine, magnetite and brucite requires no addition of silica. The model calculations suggest that the iron distribution observed in serpentine and brucite is consistent with formation temperatures ranging from <150 to 250 °C and bulk water-to-rock ratios between 0.1 and 5. These conditions coincide with peak hydrogen fugacities during serpentinization and are conducive to awaruite formation during main stage serpentinization. The development of the common brucite rims around olivine is either due to an arrested reaction olivine -> brucite -> serpentine + brucite, or reflects metastable olivine-brucite equilibria developing in the strong gradient in silica activity between orthopyroxene (talc-serpentine) and olivine (serpentine-brucite).

Radiocarbon ages and stable isotope measurements on sediment cores of the subarctic Pacific and its marginal seas

Under modern conditions only North Pacific Intermediate Water is formed in the northwest Pacific Ocean. This situation might have changed in the past. Recent studies with general circulation models indicate a switch to deep-water formation in the northwest Pacific during Heinrich Stadial 1 (17.5-15.0 ka) of the last glacial termination. Reconstructions of past ventilation changes based on paleoceanographic proxy records are still insufficient to test whether a deglacial mode of deep-water formation in the North Pacific Ocean existed. Here we present deglacial ventilation records based on radiocarbon-derived ventilation ages in combination with epibenthic stable carbon isotopes from the northwest Pacific including the Okhotsk Sea and Bering Sea, the two potential source regions for past North Pacific ventilation changes. Evidence for most rigorous ventilation of the intermediate-depth North Pacific occurred during Heinrich Stadial 1 and the Younger Dryas, simultaneous to significant reductions in Atlantic Meridional Overturning Circulation. Concurrent changes in d13C and ventilation ages point to the Okhotsk Sea as driver of millennial-scale changes in North Pacific Intermediate Water ventilation during the last deglaciation. Our records additionally indicate that changes in the d13C intermediate-water (700-1750 m water depth) signature and radiocarbon-derived ventilation ages are in antiphase to those of the deep North Pacific Ocean (>2100 m water depth) during the last glacial termination. Thus, intermediate- and deep-water masses of the northwest Pacific have a differing ventilation history during the last deglaciation.

40Ar-39Ar dating and isotopic composition of groundmass and felspar samples from core AND-2A

The AND-2A drillcore (Antarctic Drilling Program-ANDRILL) was successfully completed in late 2007 on the Antarctic continental margin (Southern McMurdo Sound, Ross Sea) with the aim of tracking ice proximal to shallow marine environmental fluctuations and to document the 20-Ma evolution of the Erebus Volcanic Province. Lava clasts and tephra layers from the AND-2A drillcore were investigated from a petrographic and stratigraphic point of view and analyzed by the 40Ar-39Ar laser technique in order to constrain the age model of the core and to gain information on the style and nature of sediment deposition in the Victoria Land Basin since Early Miocene. Ten out of 17 samples yielded statistically robust 40Ar-39Ar ages, indicating that the AND-2A drillcore recovered <230 m of Middle Miocene (~128-358 m below sea floor, ~11.5-16.0 Ma) and >780 m of Early Miocene (~358-1093 m below sea floor, ~16.0-20.1 Ma). Results also highlight a nearly continuous stratigraphic record from at least 358 m below sea floor down hole, characterized by a mean sedimentation rate of ~19 cm/ka, possible oscillations of no more than a few hundreds of ka and a break within ~17.5-18.1 Ma. Comparison with available data from volcanic deposits on land, suggests that volcanic rocks within the AND-2A core were supplied from the south, possibly with source areas closer to the drill site for the upper core levels, and from 358 m below sea floor down hole, with the 'proto-Mount Morning' as the main source.