Geochemistry of marine sediments of the world's major subduction zones

This paper presents new major and trace-element data and Lu-Hf and Sm-Nd isotopic compositions for representative suites of marine sediment samples from 14 drill sites outboard of the world's major subduction zones. These suites and samples were chosen to represent the global range in lithology, Lu/Hf ratios, and sediment flux in subducting sediments worldwide. The data reported here represent the most comprehensive data set on subducting sediments and define the Hf-Nd isotopic variations that occur in oceanic sediments and constrain the processes that caused them. Using new marine sediment data presented here, in conjunction with published data, we derive a new Terrestrial Array given by the equation, epsilon-Hf = 1.55 * epsiolon-Nd + 1.21. This array was calculated using >3400 present-day Hf and Nd isotope values. The steeper slope and smaller y-intercept of this array, compared to the original expression (epsilon-Hf = 1.36 * epsilonNd + 2.89; Vervoort et al., 1999, doi:10.1016/S0012-821X(99)00047-3) reflects the use of present day values and the unradiogenic Hf of old continental samples included in the array. In order to examine the Hf-Nd isotopic variations in marine sediments, we have classified our samples into 5 groups based on lithology and major and trace-element geochemical compositions: turbidites, terrigenous clays, and volcaniclastic, hydrothermal and hydrogenetic sediments. Compositions along the Terrestrial Array are largely controlled by terrigenous material derived from the continents and delivered to the ocean basins via turbidites, volcaniclastic sediments, and volcanic inputs from magmatic arcs. Compositions below the Terrestrial Array derive from unradiogenic Hf in zircon-rich turbidites. The anomalous compositions above the Terrestrial Array largely reflect the decoupled behavior of Hf and Nd during continental weathering and delivery to the ocean. Both terrigenous and hydrogenetic clays possess anomalously radiogenic Hf, reflecting terrestrial sedimentary and weathering processes on the one hand and marine inheritance on the other. This probably occurs during complementary processes involving preferential retention of unradiogenic Hf on the continents in the form of zircon and release of radiogenic Hf from the breakdown of easily weathered, high Lu-Hf phases such as apatite.

Geochemistry and age determination of zircons in Cretaceous to Quaternary sediments

In central Antarctica, drainage today and earlier back to the Paleozoic radiates from the Gamburtsev Subglacial Mountains (GSM). Proximal to the GSM past the Permian-Triassic fluvial sandstones in the Prince Charles Mountains (PCM) are Cretaceous, Eocene, and Pleistocene sediment in Prydz Bay (ODP741, 1166, and 1167) and pre-Holocene sediment in AM04 beneath the Amery Ice Shelf. We analysed detrital zircons for U-Pb ages, Hf-isotope compositions, and trace elements to determine the age, rock type, source of the host magma, and "crustal" model age (T(C)DM). These samples, together with others downslope from the GSM and the Vostok Subglacial Highlands (VSH), define major clusters of detrital zircons interpreted as coming from (1) 700 to 460 Ma mafic granitoids and alkaline rock, epsilon-Hf 9 to -28, signifying derivation 2.5 to 1.3 Ga from fertile and recycled crust, and (2) 1200-900 Ma mafic granitoids and alkaline rock, epsilon-Hf 11 to -28, signifying derivation 1.8 to 1.3 Ga from fertile and recycled crust. Minor clusters extend to 3350 Ma. Similar detrital zircons in Permian-Triassic, Ordovician, Cambrian, and Neoproterozoic sandstones located along the PaleoPacific margin of East Antarctica and southeast Australia further downslope from central Antarctica reflect the upslope GSM-VSH nucleus of the central Antarctic provenance as a complex of 1200-900 Ma (Grenville) mafic granitoids and alkaline rocks and older rocks embedded in 700-460 Ma (Pan-Gondwanaland) fold belts. The wider central Antarctic provenance (CAP) is tentatively divided into a central sector with negative ?Hf in its 1200-900 Ma rocks bounded on either side by positive epsilon-Hf. The high ground of the GSM-VSH in the Permian and later to the present day is attributed to crustal shortening by far-field stress during the 320 Ma mid-Carboniferous collision of Gondwanaland and Laurussia. Earlier uplifts in the ~500 Ma Cambrian possibly followed the 700-500 Ma assembly of Gondwanaland, and in the Neoproterozoic the 1000-900 Ma collisional events in the Eastern Ghats-Rayner Province at the end of the 1300-1000 Ma assembly of Rodinia.

Trace element geochemistry of zircons frim in situ ocean lithoshere

We characterize the textural and geochemical features of ocean crustal zircon recovered from plagiogranite, evolved gabbro, and metamorphosed ultramafic host-rocks collected along present-day slow and ultraslow spreading mid-ocean ridges (MORs). The geochemistry of 267 zircon grains was measured by sensitive high-resolution ion microprobe-reverse geometry at the USGS-Stanford Ion Microprobe facility. Three types of zircon are recognized based on texture and geochemistry. Most ocean crustal zircons resemble young magmatic zircon from other crustal settings, occurring as pristine, colorless euhedral (Type 1) or subhedral to anhedral (Type 2) grains. In these grains, Hf and most trace elements vary systematically with Ti, typically becoming enriched with falling Ti-in-zircon temperature. Ti-in-zircon temperatures range from 1,040 to 660°C (corrected for a TiO2 ~ 0.7, a SiO2 ~ 1.0, pressure ~ 2 kbar); intra-sample variation is typically ~60-15°C. Decreasing Ti correlates with enrichment in Hf to ~2 wt%, while additional Hf-enrichment occurs at relatively constant temperature. Trends between Ti and U, Y, REE, and Eu/Eu* exhibit a similar inflection, which may denote the onset of eutectic crystallization; the inflection is well-defined by zircons from plagiogranite and implies solidus temperatures of ~680-740°C. A third type of zircon is defined as being porous and colored with chaotic CL zoning, and occurs in ~25% of rock samples studied. These features, along with high measured La, Cl, S, Ca, and Fe, and low (Sm/La)N ratios are suggestive of interaction with aqueous fluids. Non-porous, luminescent CL overgrowth rims on porous grains record uniform temperatures averaging 615 ± 26°C (2SD, n = 7), implying zircon formation below the wet-granite solidus and under water-saturated conditions. Zircon geochemistry reflects, in part, source region; elevated HREE coupled with low U concentrations allow effective discrimination of ~80% of zircon formed at modern MORs from zircon in continental crust. The geochemistry and textural observations reported here serve as an important database for comparison with detrital, xenocrystic, and metamorphosed mafic rock-hosted zircon populations to evaluate provenance.

Radiogenic Hf isotopic composition of continental eolian dust in Pacific Ocean sediments

The inorganic silicate fraction extracted from bulk pelagic sediments from the North Pacific Ocean is eolian dust. It monitors the composition of continental crust exposed to erosion in Asia. 176Lu/177Hf ratios of modern dust are subchondritic between 0.011 and 0.016 but slightly elevated with respect to immature sediments. Modern dust samples display a large range in Hf isotopic composition (IC), -4.70 < epsilon-Hf < +16.45, which encompasses that observed for the time series of DSDP cores 885/886 and piston core LL44-GPC3 extending back to the late Cretaceous. Hafnium and neodymium isotopic results are consistent with a dominantly binary mixture of dust contributed from island arc volcanic material and dust from central Asia. The Hf-Nd isotopic correlation for all modern dust samples, epsilon-Hf= =0.78 epsilon-Nd = +5.66 (n =22, R**2 =0.79), is flatter than those reported so far for terrestrial reservoirs. Moreover, the variability in epsilon-Hf of Asian dust exceeds that predicted on the basis of corresponding epsilon-Nd values (34.76 epsilon-Hf < +2.5; -10.96< epsilon-Nd <-10.1). This is attributed to: (1) the fixing of an important unradiogenic fraction of Hf in zircons, balanced by radiogenic Hf that is mobile in the erosional cycle, (2) the elevated Lu/Hf ratio in chemical sediments which, given time, results in a Hf signature that is radiogenic compared with Hf expected from its corresponding Nd isotopic components, and (3) the possibility that diagenetic resetting of marine sediments may incorporate a significant radiogenic Hf component into diagenetically grown minerals such as illite. Together, these processes may explain the variability and more radiogenic character of Hf isotopes when compared to the Nd isotopic signatures of Asian dust. The Hf-Nd isotope time series of eolian dust are consistent with the results of modern dust except two samples that have extremely radiogenic Hf for their Nd (epsilon-Hf =+8.6 and +10.3, epsilon-Nd =39.5 and 39.8). These data may point to a source contribution of dust unresolved by Nd and Pb isotopes. The Hf IC of eolian dust input to the oceans may be more variable and more radiogenic than previously anticipated. The Hf signature of Pacific seawater, however, has varied little over the past 20 Myr, especially across the drastic increase of eolian dust flux from Asia around 3.5 Ma. Therefore, continental contributions to seawater Hf appear to be riverine rather than eolian. Current predictions regarding the relative proportions of source components to seawater Hf must account for the presence of a variable and radiogenic continental component. Data on the IC and flux of river-dissolved Hf to the oceans are urgently required to better estimate contributions to seawater Hf. This then would permit the use of Hf isotopes as a monitor of past changes in erosion.

Chemistry of minerals from amphibole-chlorite veins of ODP Site 209-1270

We examined small-scale shear zones in drillcore samples of abyssal peridotites from the Mid-Atlantic Ridge. These shear zones are associated with veins consisting of chlorite + actinolite/tremolite assemblages, with accessory phases zircon and apatite, and they are interpreted as altered plagiogranite melt impregnations, which originate from hydrous partial melting of gabbroic intrusion in an oceanic detachment fault. Ti-in-zircon thermometry yields temperatures around 820°C for the crystallization of the evolved melt. Reaction path modeling indicates that the alteration assemblage includes serpentine of the adjacent altered peridotites. Based on the model results, we propose that formation of chlorite occurred at higher temperatures than serpentinization, thus leading to strain localization around former plagiogranites during alteration. The detachment fault represents a major pathway for fluids through the oceanic crust, as evidenced by extremely low d18O of altered plagiogranite veins (+3.0-4.2 per mil) and adjacent serpentinites (+ 2.6-3.7 per mil). The uniform oxygen isotope data indicate that fluid flow in the detachment fault system affected veins and adjacent host serpentinites likewise.

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.

Stable oxygen isotope ratios in zircons from oceanic crust

Lower ocean crust is primarily gabbroic, although 1-2% felsic igneous rocks that are referred to collectively as plagiogranites occur locally. Recent experimental evidence suggests that plagiogranite magmas can form by hydrous partial melting of gabbro triggered by seawater-derived fluids, and thus they may indicate early, high-temperature hydrothermal fluid circulation. To explore seawater-rock interaction prior to and during the genesis of plagiogranite and other late-stage magmas, oxygen-isotope ratios preserved in igneous zircon have been measured by ion microprobe. A total of 197 zircons from 43 plagiogranite, evolved gabbro, and hydrothermally altered fault rock samples have been analyzed. Samples originate primarily from drill core acquired during Ocean Drilling Program and Integrated Ocean Drilling Program operations near the Mid-Atlantic and Southwest Indian Ridges. With the exception of rare, distinctively luminescent rims, all zircons from ocean crust record remarkably uniform d18O with an average value of 5.2 ± 0.5 per mil (2SD). The average d18O(Zrc) would be in magmatic equilibrium with unaltered MORB [d18O(WR) ~5.6-5.7 per mil], and is consistent with the previously determined value for equilibrium with the mantle. The narrow range of measured d18O values is predicted for zircon crystallization from variable parent melt compositions and temperatures in a closed system, and provides no indication of any interactions between altered rocks or seawater and the evolved parent melts. If plagiogranite forms by hydrous partial melting, the uniform mantle-like d18O(Zrc) requires melting and zircon crystallization prior to significant amounts of water-rock interactions that alter the protolith d18O. Zircons from ocean crust have been proposed as a tectonic analog for >3.9 Ga detrital zircons from the earliest (Hadean) Earth by multiple workers. However, zircons from ocean crust are readily distinguished geochemically from zircons formed in continental crustal environments. Many of the >3.9 Ga zircons have mildly elevated d18O (6.0-7.5 per mil), but such values have not been identified in any zircons from the large sample suite examined here. The difference in d18O, in combination with newly acquired lithium concentrations and published trace element data, clearly shows that the >3.9 Ga detrital zircons did not originate by processes analogous to those in modern mid-ocean ridge settings.

(Table 1, page 641) Hf isotope time series for the uppermost part of a ferromanganese crust from the San Pablo ferro-manganese pavement

The Hf isotope composition of seawater does not match that expected from dissolution of bulk continental crust. This mismatch is generally considered to be due to retention of unradiogenic Hf in resistant zircons during incomplete weathering of continental crust. During periods of intense glacial weathering, zircons should break down more efficiently, resulting in the release of highly unradiogenic Hf to the oceans. We test this hypothesis by comparing Nd and Hf isotope time series obtained from NW Atlantic ferromanganese crusts. Both isotope systems show a decrease associated with the onset of northern hemisphere glaciation. The observed changes display distinct trajectories in epsilon Nd- epsilon Hf space, which differ from previously reported arrays of bulk terrestrial material and seawater. Such patterns are consistent with the release of highly unradiogenic Hf from very old zircons, facilitated by enhanced mechanical weathering.

Composition of melt inclusions and age of zircons of plagiogneisses from the Kola Superdeep Borehole

A comprehensive study of melt inclusions and SHRIMP dating of zircons from trondhjemite gneisses of the sequence VIII from the Kola Superdeep Borehole has revealed presence of old primary magmatic crystals with age up to 2887+/-15 Ma. This is not consistent with the previous view, according to which the oldest zircons from the Archean Complex in SG-3 are products of granulite metamorphism. Primary magmatic zircons of early generation (from 2887 to 2842 Ma) formed in deep-seated magma chambers during partial crystallization of CO2-saturated trondhjemite estimates on duration of generation of tonalite-trondhjemite-granite melts through partial melting of mafic rocks.