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.

Calcium carbonate global LGM Burial rate data

Global databases of calcium carbonate concentrations and mass accumulation rates in Holocene and last glacial maximum sediments were used to estimate the deep-sea sedimentary calcium carbonate burial rate during these two time intervals. Sparse calcite mass accumulation rate data were extrapolated across regions of varying calcium carbonate concentration using a gridded map of calcium carbonate concentrations and the assumption that accumulation of noncarbonate material is uncorrelated with calcite concentration within some geographical region. Mean noncarbonate accumulation rates were estimated within each of nine regions, determined by the distribution and nature of the accumulation rate data. For core-top sediments the regions of reasonable data coverage encompass 67% of the high-calcite (>75%) sediments globally, and within these regions we estimate an accumulation rate of 55.9 ± 3.6 x 10**11 mol/yr. The same regions cover 48% of glacial high-CaCO3 sediments (the smaller fraction is due to a shift of calcite deposition to the poorly sampled South Pacific) and total 44.1 ± 6.0 x 10**11 mol/yr. Projecting both estimates to 100 % coverage yields accumulation estimates of 8.3 x 10**12 mol/yr today and 9.2 x 10**12 mol/yr during glacial time. This is little better than a guess given the incomplete data coverage, but it suggests that glacial deep sea calcite burial rate was probably not considerably faster than today in spite of a presumed decrease in shallow water burial during glacial time.