Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Glaciological characteristics and carbon dioxide record of ice cores from Antarctica

We report a decadally resolved record of atmospheric CO2 concentration for the last 1000 years, obtained from the West Antarctic Ice Sheet (WAIS) Divide shallow ice core. The most prominent feature of the pre-industrial period is a rapid ~7 ppm decrease of CO2 in a span of ~20-50 years at ~1600 A.D. This observation confirms the timing of an abrupt atmospheric CO2 decrease of ~10 ppm observed for that time period in the Law Dome ice core CO2 records, but the true magnitude of the decrease remains unclear. Atmospheric CO2 variations over the time period 1000-1800 A.D. are statistically correlated with northern hemispheric climate and tropical Indo-Pacific sea surface temperature. However, the exact relationship between CO2 and climate remains elusive due to regional climate variations and/or uneven geographical data density of paleoclimate records. We observe small differences of 0 ~2% (0 ~ 6 ppm) among the high-precision CO2 records from the Law Dome, EPICA Dronning Maud Land and WAIS Divide Antarctic ice cores. However, those records share common trends of CO2 change on centennial to multicentennial time scales, and clearly show that atmospheric CO2 has been increasing above preindustrial levels since ~1850 A.D.

Strontium isotope ratios of DSDP samples

The isotopic ratio of strontium-87 to strontium-86 shows no detectable variation in present-day ocean water but changes slowly over millions of years. The strontium contained in carbonate shells of marine organisms records the ratio of strontium-87 to strontium-86 of the oceans at the time that the shells form. Sedimentary rocks composed of accumulated fossil carbonate shells can be dated and correlated with the use of high precision measurements of the ratio of strontium-87 to strontium-86 with a resolution that is similar to that of other techniques used in age correlation. This method may prove valuable for many geological, paleontological, paleooceanographic, and geochemical problems.

Chemical and isotopic compositions of sediments from DSDP Hole 21-204

The highly depleted intra-oceanic Tonga-Kermadec island arc forms an endmember of arc systems and a unique location in which to isolate the effects of the slab flux. High precision TIMS uranium, thorium, strontium, neodymium, and lead isotopes, along with complete major and trace element data, have been obtained on an extensive sample set comprising fifty-eight lavas along the arc as well as nineteen samples of the subducting sediments at DSDP site 204 just to the east of the Tonga-Kermadec trench. Ca/Ti and Al/Ti ratios extend from values appropriate to an N-MORB source in the southern Kermadecs to very high ratios in Tonga interpreted to reflect increasing degrees of depletion of the mantle wedge due to backarc basalt extraction. The isotope data emphasize the need for four components in the petrogenesis of the lavas: (1) the mantle wedge; (2) a component with elevated 207Pb/204Pb towards which the Kermadec and southern Tongan lavas extend; (3) a component characterised by high 206Pb/204Pb, Ta/Nd, and low 143Nd/144Nd observed only in the northernmost Tongan islands of Tafahi and Niuatoputapu; (4) a fluid component characterised by strong enrichments of Rb, Ba, U, K, Ph, and Sr, relative to Th, Zr, and the REE and producing large 238U excesses ((230Th/238U) = 0.8-0.5) in the more depleted lavas. The mantle wedge (Component 1) is isotopically similar to the source of the Lau BABB. Component 2 is average pelagic sediment on the downgoing Pacific plate as observed at DSDP sites 595/596 and in the upper sections of the sediment pile at DSDP site 204. Mass balance calculations indicate that less than 0.5% is recycled into the arc lavas; essentially all the subducted sediment is returned to the upper mantle (~0.03 km**3/yr). Exceptionally low concentrations of Ta and Nb relative to Th and the LREE requires that this sediment component is added as a partial melt which was in equilibrium with residual rutile or ilmenite. Component 3 is identified as volcaniclastics from the Louisville Ridge which comprise the lower 44 m of the sediment section intersected at DSDP site 204. These volcaniclastics are spatially restricted to the vicinity of the Louisville Ridge and provide a unique sediment tracer which can be used to show that it takes 4 Myr from the time of subduction to its first appearance in the arc lava signature. Component 4, the fluid contribution to the lava source is inferred to contribute ~1 ppm Rb, 10 ppm Ba, 0.02 ppm U, 600 ppm K, 0.2 ppm Ph, and 30 ppm Sr. It has 87Sr/86Sr = 0.7035 and 206Pb/204Pb = 18.5 and thus it is inferred to have been derived from dehydration of the subducting altered oceanic crust. U-Th isotope disequilibria reflect the time since fluid release from the subducting slab and a reference line through the lowest (230Th/232Th) lavas constrains this to be 30000-50000 yr. The U-Th and Th-Ra isotope systematics are decoupled, and it is suggested that Th-Ra isotope disequilibria record the time since partial melting and thus indicate rapid channelled magma ascent. Olivine gabbro xenoliths from Raoul are interpreted as cumulates to their host lavas with which they form zero age U-Th isochrons indicating that minimal time was spent in magma chambers. The subduction signature is not observed in lavas from the backarc island of Niuafo'ou. These were derived from partial melting of fertile peridotite at 130-160 km depth with melt rates around 0.0002 kg/m**3/yr.

Geochemical provenance analysis of fine-grained sediment of sediment core AND1-1

Bulk chemical fine-grained sediment compositions from southern Victoria Land glacimarine sediments provide significant constraints on the reconstruction of sediment provenance models in the McMurdo Sound during Late Cenozoic time. High-resolution (~ 1 ka) geochemical data were obtained with a non-destructive AVAATECH XRF Core Scanner (XRF-CS) on the 1285 m long ANDRILL McMurdo Ice Shelf Project (MIS) sediment core AND-1B. This data set is complemented by high-precision chemical analyses (XRF and ICP-OES) on discrete samples. Statistical analyses reveal three geochemical facies which are interpreted to represent the following sources for the sediments recovered in the AND-1B core: 1) local McMurdo Volcanic Group (MVG) rocks, 2) Transantarctic Mountain rocks west of Ross Island (W TAM), and 3) Transantarctic Mountain rocks from more southerly areas (S TAM). Data indicate in combination with other sediment facies analyses (McKay et al., 2009, doi:10.1130/B26540.1) and provenance scenarios (Talarico and Sandroni, 2009, doi:10.1016/j.gloplacha.2009.04.007) that diamictites at the drill site are largely dominated by local sources (MVG) and are interpreted to indicate cold polar conditions with dry-based ice. MVG is interpreted to indicate cold polar condition with dry-based ice. A mixture of MVG and W TAM is interpreted to represent polar conditions and the S TAM facies is interpreted to represent open-marine conditions. Down-core variations in geochemical facies in the AND-1B core are interpreted to represent five major paleoclimate phases over the past 14 Ma. Cold polar conditions with major MVG influence occur below 1045 mbsf and above 120 mbsf. A section of warmer climate conditions with extensive peaks of S TAM influence characterizes the rest of the core, which is interrupted by a section from 525 to 855 mbsf of alternating influences of MVG and W TAM.

Geochemistry of North Atlantic Igneous Province basalts

New Sr-Nd-Pb-Hf data require the existence of at least four mantle components in the genesis of basalts from the the North Atlantic Igneous Province (NAIP): (1) one (or more likely a small range of) enriched component(s) within the Iceland plume, (2) a depleted component within the Iceland plume (distinct from the shallow N-MORB source), (3) a depleted sheath surrounding the plume and (4) shallow N-MORB source mantle. These components have been available since the major phase of igneous activity associated with plume head impact during Paleogene times. In Hf-Nd isotope space, samples from Iceland, DSDP Leg 49 (Sites 407, 408 and 409), ODP Legs 152 and 163 (southeast Greenland margin), the Reykjanes Ridge, Kolbeinsey Ridge and DSDP Leg 38 (Site 348) define fields that are oblique to the main ocean island basalt array and extend toward a component with higher 176Hf/177Hf than the N-MORB source available prior to arrival of the plume, as indicated by the compositions of Cretaceous basalts from Goban Spur (~95 Ma). Aside from Goban Spur, only basalts from Hatton Bank on the oceanward side of the Rockall Plateau (DSDP Leg 81) lie consistently within the field of N-MORB, which indicates that the compositional influence of the plume did not reach this far south and east ~55 Ma ago. Thus, Hf-Nd isotope systematics are consistent with previous studies which indicate that shallow MORB-source mantle does not represent the depleted component within the Iceland plume (Thirlwall, J. Geol. Soc. London 152 (1995) 991-996; Hards et al., J. Geol. Soc. London 152 (1995) 1003-1009; Fitton et al., 1997 doi:10.1016/S0012-821X(97)00170-2). They also indicate that the depleted component is a long-lived and intrinsic feature of the Iceland plume, generated during an ancient melting event in which a mineral (such as garnet) with a high Lu/Hf was a residual phase. Collectively, these data suggest a model for the Iceland plume in which a heterogeneous core, derived from the lower mantle, consists of 'enriched' streaks or blobs dispersed in a more depleted matrix. A distinguishing feature of both the enriched and depleted components is high Nb/Y for a given Zr/Y (i.e. positive DeltaNb), but the enriched component has higher Sr and Pb isotope ratios, combined with lower epsilon-Nd and epsilon-Hf. This heterogeneous core is surrounded by a sheath of depleted material, similar to the depleted component of the Iceland plume in its epsilon-Nd and epsilon-Hf, but with lower 87Sr/86Sr, 208Pb/204Pb and negative DeltaNb; this material was probably entrained from near the 670 km discontinuity when the plume stalled at the boundary between the upper and lower mantle. The plume sheath displaced more normal MORB asthenosphere (distinguished by its lower epsilon-Hf for a given epsilon-Nd or Zr/Nb ratio), which existed in the North Atlantic prior to plume impact. Preliminary data on MORBs from near the Azores plume suggest that much of the North Atlantic may be 'polluted' not only by enriched plume material but also by depleted material similar to the Iceland plume sheath. If this hypothesis is correct, it may provide a general explanation for some of the compositional diversity and variations in inferred depth of melting (Klein and Langmuir, 1987 doi:10.1029/JB092iB08p08089) along the MAR in the North Atlantic.

Lead isotopes inthe Eastern Equatorial Pacific

The influence of atmospheric dust on climate and biogeochemical cycles in the oceans is well understood but poorly quantified. Glacial atmospheric dust loads were generally greater than those during the Holocene, as shown, for example, by the covariation of dust fluxes in the Equatorial Pacific and Antarctic ice cores. Nevertheless, it remains unclear whether these increases in dust flux were associated with changes in sources of dust, which would in turn suggest variations in wind patterns, climate or paleo-environment. Such questions can be answered using radiogenic isotope tracers of dust provenance. Here, we present a 160-kyr high-precision lead isotope time-series of dust input to the Eastern Equatorial Pacific (EEP) from core ODP Leg 138, Site 849 (0°11.59' N, 110°31.18' W). The Pb isotope record, combined with Nd isotope data, rules out contributions from Northern Hemisphere dust sources, north of the Intertropical Convergence Zone, such as Asia or North Africa/Sahara; similarly, eolian sources in Australia, Central America, the Northern Andes and Patagonia appear insignificant based upon the radiogenic isotope data. Fluctuations in Pb isotope ratios throughout the last 160 kyr show, instead, that South America remained the prevailing source of dusts to the EEP. There are two distinct South American Pb isotope end-members, constrained to be located in the south Central Volcanic Zone (CVZ, 22° S - 27.5° S) and the South Volcanic Zone (SVZ, 33° S - 43° S), with the former most likely originating in the Atacama Desert. Dust availability in the SVZ appears to be related to the weathering of volcanic deposits and the development of ash-derived Andosols, and influenced by local factors that might include vegetation cover. Variations in the dust fluxes from the two sources are in phase with both the dust flux and temperature records from Antarctican ice cores. We show that the forcing of dust provenance over time in the EEP overall is influenced by high-southerly-latitude climate conditions, leading to changes in the latitudinal position and strength of the South Westerlies as well as the coastal winds that blow northward along the Chilean margin. The net result is a modulation of dust emission from the Atacama Desert and the SVZ via a northward migration of the South Westerlies during cold periods and southward retreat during glacial terminations.

Ages and stable isotope ratios from speleothem calcite from Cold Water Cave, Iowa

High-precision uranium-thorium mass spectrometric chronology and 18O-13C isotopic analysis of speleothem calcite from Cold Water Cave in northeast Iowa have been used to chart mid-Holocene climate change. Significant shifts in d18O and d13C isotopic values coincide with well-documented Holocene vegetation changes. Temperature estimates based on 18O/16O ratios suggest that the climate warmed rapidly by about 3°C at 5900 years before present and then cooled by 4°C at 3600 years before present. Initiation of a gradual increase in ?d13C at 5900 years before present suggests that turnover of the forest soil biomass was slow and that equilibrium with prairie vegetation was not attained by 3600 years before present.

Diatom biostratigraphy, Argon ages and chronology of ODP Hole 183-1138A

A thick Neogene section was recovered in the upper ~300 m of Ocean Drilling Program Hole 1138A, drilled on the Central Kerguelen Plateau in the Indian sector of the Southern Ocean. Sediment lithologies consist primarily of mixed carbonate and biosiliceous clays and oozes, with several thin (1-3 cm) tephra horizons. The tephras are glass rich, well sorted, and dominantly trachytic to rhyolitic in composition. Volcaniclastic material in these horizons is interpreted to have originated from Heard Island, 180 km northwest of Site 1138, and was likely emplaced through both primary ash fall and turbiditic, submarine flows. A Neogene age-depth model for Hole 1138A is constructed primarily from 36 diatom biostratigraphic datums. Nannofossil and planktonic foraminifer biostratigraphy provides supporting age information. Additionally, four high-precision 40Ar-39Ar ages are derived from ash and tephra horizons, and these radiometric ages are in close agreement with the biostratigraphic ages. The integrated age-depth model reveals a reasonably complete lower Miocene to upper Pleistocene section in Hole 1138A, with the exception of a ~1-m.y. hiatus at the Miocene/Pliocene boundary. Another possible hiatus is also identified at the Oligocene/Miocene boundary. High Neogene sedimentation rates and the presence of both calcareous and siliceous microfossils, combined with datable tephra horizons, establish Site 1138 as a suitable target for future drilling legs with paleoceanographic objectives. This report also proposes two new diatom species, Fragilariopsis heardensis and Azpeitia harwoodii, from Pliocene strata of Hole 1138A.

Grain size and XRD analysis on Site 162-983

The grain size of deep-sea sediments provides an apparently simple proxy for current speed. However, grain size-based proxies may be ambiguous when the size distribution reflects a combination of processes, with current sorting only one of them. In particular, such sediment mixing hinders reconstruction of deep circulation changes associated with ice-rafting events in the glacial North Atlantic because variable ice-rafted detritus (IRD) input may falsely suggest current speed changes. Inverse modeling has been suggested as a way to overcome this problem. However, this approach requires high-precision size measurements that register small changes in the size distribution. Here we show that such data can be obtained using electrosensing and laser diffraction techniques, despite issues previously raised on the low precision of electrosensing methods and potential grain shape effects on laser diffraction. Down-core size patterns obtained from a sediment core from the North Atlantic are similar for both techniques, reinforcing the conclusion that both techniques yield comparable results. However, IRD input leads to a coarsening that spuriously suggests faster current speed. We show that this IRD influence can be accounted for using inverse modeling as long as wide size spectra are taken into account. This yields current speed variations that are in agreement with other proxies. Our experiments thus show that for current speed reconstruction, the choice of instrument is subordinate to a proper recognition of the various processes that determine the size distribution and that by using inverse modeling meaningful current speed reconstructions can be obtained from mixed sediments.