Molecular and isotopic properties of gases from ODP Holes of ODP Leg 204

We report and discuss molecular and isotopic properties of hydrate-bound gases from 55 samples and void gases from 494 samples collected during Ocean Drilling Program (ODP) Leg 204 at Hydrate Ridge offshore Oregon. Gas hydrates appear to crystallize in sediments from two end-member gas sources (deep allochthonous and in situ) as mixtures of different proportions. In an area of high gas flux at the Southern Summit of the ridge (Sites 1248-1250), shallow (0-40 m below the seafloor [mbsf]) gas hydrates are composed of mainly allochthonous mixed microbial and thermogenic methane and a small portion of thermogenic C2+ gases, which migrated vertically and laterally from as deep as 2- to 2.5-km depths. In contrast, deep (50-105 mbsf) gas hydrates at the Southern Summit (Sites 1248 and 1250) and on the flanks of the ridge (Sites 1244-1247) crystallize mainly from microbial methane and ethane generated dominantly in situ. A small contribution of allochthonous gas may also be present at sites where geologic and tectonic settings favor focused vertical gas migration from greater depth (e.g., Sites 1244 and 1245). Non-hydrocarbon gases such as CO2 and H2S are not abundant in sampled hydrates. The new gas geochemical data are inconsistent with earlier models suggesting that seafloor gas hydrates at Hydrate Ridge formed from gas derived from decomposition of deeper and older gas hydrates. Gas hydrate formation at the Southern Summit is explained by a model in which gas migrated from deep sediments, and perhaps was trapped by a gas hydrate seal at the base of the gas hydrate stability zone (GHSZ). Free gas migrated into the GHSZ when the overpressure in gas column exceeded sealing capacity of overlaying sediments, and precipitated as gas hydrate mainly within shallow sediments. The mushroom-like 3D shape of gas hydrate accumulation at the summit is possibly defined by the gas diffusion aureole surrounding the main migration conduit, the decrease of gas solubility in shallow sediment, and refocusing of gas by carbonate and gas hydrate seals near the seafloor to the crest of the local anticline structure.

Platinum-group, major, and trace element abundances in ODP Leg 183 basalts

Seventeen basalts from Ocean Drilling Program (ODP) Leg 183 to the Kerguelen Plateau (KP) were analyzed for the platinum-group elements (PGEs: Ir, Ru, Rh, Pt, and Pd), and 15 were analyzed for trace elements. Relative concentrations of the PGEs ranged from ~0.1 (Ir, Ru) to ~5 (Pt) times primitive mantle. These relatively high PGE abundances and fractionated patterns are not accounted for by the presence of sulfide minerals; there are only trace sulfides present in thin-section. Sulfur saturation models applied to the KP basalts suggest that the parental magmas may have never reached sulfide saturation, despite large degrees of partial melting (~30%) and fractional crystallization (~45%). First order approximations of the fractionation required to produce the KP basalts from an ~30% partial melt of a spinel peridotite were determined using the PELE program. The model was adapted to better fit the physical and chemical observations from the KP basalts, and requires an initial crystal fractionation stage of at least 30% olivine plus Cr-spinel (49:1), followed by magma replenishment and fractional crystallization (RFC) that included clinopyroxene, plagioclase, and titanomagnetite (15:9:1). The low Pd values ([Pd/Pt]_pm < 1.7) for these samples are not predicted by currently available Kd values. These Pd values are lowest in samples with relatively higher degrees of alteration as indicated by petrographic observations. Positive anomalies are a function of the behavior of the PGEs; they can be reproduced by Cr-spinel, and titanomagnetite crystallization, followed by titanomagnetite resorption during the final stages of crystallization. Our modeling shows that it is difficult to reproduce the PGE abundances by either depleted upper or even primitive mantle sources. Crustal contamination, while indicated at certain sites by the isotopic compositions of the basalts, appears to have had a minimal affect on the PGEs. The PGE abundances measured in the Kerguelen Plateau basalts are best modeled by melting a primitive mantle source to which was added up to 1% of outer core material, followed by fractional crystallization of the melt produced. This reproduces both the abundances and patterns of the PGEs in the Kerguelen Plateau basalts. An alternative model for outer core PGE abundances requires only 0.3% of outer core material to be mixed into the primitive mantle source. While our results are clearly model dependent, they indicate that an outer core component may be present in the Kerguelen plume source.

Chemical composition of DSDP and ODP cherts

Rare earth element (REE), major, and trace element abundances and relative fractionations in forty nodular cherts sampled by the Deep Sea Drilling Project (DSDP) and Ocean Drilling Program (ODP) indicate that the REE composition of chert records the interplay between terrigenous sources and scavenging from the local seawater. Major and (non-REE) trace element ratios indicate that the aluminosilicate fraction within the chert is similar to NASC (North American Shale Composite), with average Pacific chert including ~7% NASC-like particles, Indian chert ~11% NASC, Atlantic chert ~17% NASC, and southern high latitude (SHL) chert 53% NASC. Using La as a proxy for sum REE, approximations of excessive La (the amount of La in excess of that supplied by the detrital aluminosilicate fraction) indicate that Pacific chert contains the greatest excessive La (85% of total La) and SHL chert the least (38% of total La). As shown by interelement associations, this excessive La is most likely an adsorbed component onto aluminosilicate and phosphatic phases. Accordingly, chert from the large Pacific Ocean, where deposition occurs relatively removed from significant terrigenous input, records a depositional REE signal dominated by adsorption of dissolved REEs from seawater. Pacific chert Ce/Ce* <<1 and normative La/Yb ~ 0.8-1, resulting from adsorption of local Ce-depleted seawater and preferential adsorption of LREEs from seawater (e.g., normative La/Yb ~0.4), which increases the normative La/Yb ratio recorded in chert. Chert from the Atlantic basin, a moderately sized ocean basin lined by passive margins and with more terrigenous input than the Pacific, records a mix of adsorptive and terrigenous REE signals, with moderately negative Ce anomalies and normative La/Yb ratios intermediate to those of the Pacific and those of terrigenous input. Chert from the SHL region is dominated by the large terrigenous input on the Antarctic passive margin, with inherited Ce/Ce* ~1 and inherited normative La/Yb values of ~1.2-1.4. Ce/Ce* does not vary with age, either throughout the entire data base or within a particular basin. Overall, Ce/Ce* does not correlate with P2O5 concentrations, even though phosphatic phases may be an important REE carrier.

Uranium contents and isotope ratios of uranium and strontium in sediments, leachates and pore fluids from ODP Hole 162-984A

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the 234U/238U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%). The pore water 234U/238U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment 234U/238U activity ratios are close to 1.0. The 234U/238U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess 234U to the pore fluid by a-recoil injection of 234Th. The fraction of 238U decays that result in a-recoil injection of 234U to pore fluid is estimated to be 0.10 to 0.20 based on the 234U/238U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 0.0000004 to 0.000002 1/yr. There is significant down-hole variability in pore water 234U/238U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 1000 times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials. The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ~ 0.1 1/age for ages spanning 1000 to 500,000,000 yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.

Revised calibration of the MBT-CBT paleotemperature proxy

The MBT-CBT proxy for the reconstruction of paleotemperatures and past soil pH is based on the distribution of branched glycerol dialkyl glycerol tetraether (brGDGT) membrane lipids. The Methylation of Branched Tetraether (MBT) and the Cyclisation of Branched Tetraether (CBT) indices were developed to quantify these distributions, and significant empirical relations between these indices and annual mean air temperature (MAT) and/or soil pH were found in a large data set of soils. In this study, we extended this soil dataset to 278 globally distributed surface soils. Of these soils, 26% contains all nine brGDGTs, while in 63% of the soils the seven most common brGDGTs were detected, and the latter were selected for calibration purposes. This resulted in new transfer functions for the reconstruction of pH based on the CBT index: pH = 7.90-1.97 × CBT (r**2 = 0.70; RMSE = 0.8; n = 176), as well as for MAT based on the CBT index and methylation index based on the seven most abundant GDGTs (defined as MBT'): MAT = 0.81-5.67 × CBT + 31.0 × MBT' (r**2 = 0.59; RMSE = 5.0 °C; n = 176). The new transfer function for MAT has a substantially lower correlation coefficient than the original equation (r**2 = 0.77). To investigate possible improvement of the correlation, we used our extended global surface soil dataset to statistically derive the indices that best describe the relations of brGDGT composition with MAT and soil pH. These new indices, however, resulted in only a relatively minor increase in correlation coefficients, while they cannot be explained straightforwardly by physiological mechanisms. The large scatter in the calibration cannot be fully explained by local factors or by seasonality, but MAT for soils from arid regions are generally substantially (up to 20 °C) underestimated, suggesting that absolute brGDGT-based temperature records for these areas should be interpreted with caution. The applicability of the new MBT'-CBT calibration function was tested using previously published MBT-CBT-derived paleotemperature records covering the last deglaciation in Central Africa and East Asia, the Eocene-Oligocene boundary and the Paleocene-Eocene thermal maximum. The results show that trends remain similar in all records, but that absolute temperature estimates and the amplitude of temperature changes are lower for most records, and generally in better agreement with independent proxy data.

Helium isotope concentrations and ratios in sediments from ODP Hole 130-806B

We have determined the helium abundance and isotopic composition of seafloor carbonate sediments from the flanks of the Ontong Java Plateau, western equatorial Pacific Ocean (ODP Site 806). These results provide a two million year record of the burial flux of extraterrestrial 3He, which we believe is a proxy for the terrestrial accretion rate of interplanetary dust particles. The 3He burial flux prior to ~700 ka was relatively low, ~0.5 pcc/cm**2/kyr, but from 700 ka to the present, the burial flux gradually increased to a value of ~1.0 pcc/cm**2/kyr. 100 kyr periodicity in the 3He burial flux is apparent over the last 700 kyr and correlates with the oxygen isotope record of global climate, with high 3He burial fluxes associated with interglacial periods. This periodicity and phase are consistent with previous 3He measurements in North Atlantic sediments. Although 100 kyr periodicity in 3He burial flux is in agreement with recent predictions of the accretion rate of interplanetary dust based on a model of the orbital evolution of asteroidal debris, the measurements and predictions differ by one half cycle in phase. Nevertheless, our observations suggest the terrestrial accretion rate of interplanetary dust is controlled by orbital eccentricity and/or inclination relative to the solar-system invariable plane. Such control is a necessary but not sufficient condition for the hypothesis of that variations in extraterrestrial dust accretion modulates terrestrial climate with a 100 kyr period. We also identify several brief (<25 kyr) intervals of strongly enhanced 3He burial, possibly related to random and transient fluctuations in the accretion rate of asteroidal or cometary dust particles.

Concentrations and accumulation rates of platinum-group elements and rhenium in Pacific and Atlantic sediments, DSDP and ODP data

The nature of Re-platinum-group element (PGE; Pt, Pd, Ir, Os, Ru) transport in the marine environment was investigated by means of marine sediments at and across the Cretaceous-Tertiary boundary (KTB) at two hemipelagic sites in Europe and two pelagic sites in the North and South Pacific. A traverse across the KTB in the South Pacific pelagic clay core found elevated levels of Re, Pt, Ir, Os, and Ru, each of which is approximately symmetrically distributed over a distance of ~1.8 m across the KTB. The Re-PGE abundance patterns are fractionated from chondritic relative abundances: Ru, Pt, Pd, and Re contents are slightly subchondritic relative to Ir, and Os is depleted by ~95% relative to chondritic Ir proportions. A similar depletion in Os (~90%) was found in a sample of the pelagic KTB in the North Pacific, but it is enriched in Ru, Pt, Pd, and Re relative to Ir. The two hemipelagic KTB clays have near-chondritic abundance patterns. The ~1.8-m-wide Re-PGE peak in the pelagic South Pacific section cannot be reconciled with the fallout of a single impactor, indicating that postdepositional redistribution has occurred. The elemental profiles appear to fit diffusion profiles, although bioturbation could have also played a role. If diffusion had occurred over ~65 Ma, the effective diffusivities are ~10**?13 cm**2/s, much smaller than that of soluble cations in pore waters (~10**?6 cm**2/s). The coupling of Re and the PGEs during redistribution indicates that postdepositional processes did not significantly fractionate their relative abundances. If redistribution was caused by diffusion, then the effective diffusivities are the same. Fractionation of Os from Ir during the KTB interval must therefore have occurred during aqueous transport in the marine environment. Distinctly subchondritic Os/Ir ratios throughout the Cenozoic in the South Pacific core further suggest that fractionation of Os from Ir in the marine environment is a general process throughout geologic time because most of the inputs of Os and Ir into the ocean have Os/Ir ratios >/=1. Mass balance calculations show that Os and Re burial fluxes in pelagic sediments account for only a small fraction of the riverine Os (<10%) and Re (<0.1%) inputs into the oceans. In contrast, burial of Ir in pelagic sediments is similar to the riverine Ir input, indicating that pelagic sediments are a much larger repository for Ir than for Os and Re. If all of the missing Os and Re is assumed to reside in anoxic sediments in oceanic margins, the calculated burial fluxes in anoxic sediments are similar to observed burial fluxes. However, putting all of the missing Os and Re into estuarine sediments would require high concentrations to balance the riverine input and would also fail to explain the depletion of Os at pelagic KTB sites, where at most ~25% of the K-T impactor's Os could have passed through estuaries. If Os is preferentially sequestered in anoxic marine environments, it follows that the Os/Ir ratio of pelagic sediments should be sensitive to changes in the rates of anoxic sediment deposition. There is thus a clear fractionation of Os and Re from Ir in precipitation out of sea water in pelagic sections. Accordingly, it is inferred here that Re and Os are removed from sea water in anoxic marine depositional regimes.

Contents, partition, and isotopic composition of sulphur in sediments from ODP Sites 680 and 686

We have determined (1) the abundance and isotopic composition of pyrite, monosulphide, elemental sulphur, organically bound sulphur, and dissolved sulphide; (2) the partition of ferric and ferrous iron; (3) the organic carbon contents of sediments recovered at two sites drilled on the Peru Margin during Leg 112 of the Ocean Drilling Program. Sediments at both sites are characterised by high levels of organically bound sulphur (OBS). OBS comprises up to 50% of total sedimentary sulphur and up to 1% of bulk sediment. The weight ratio of S to C in organic matter varies from 0.03 to 0.15 (mean = 0.10). Such ratios are like those measured in lithologically similar, but more deeply buried petroleum source rocks of the Monterey and Sisquoc formations in California. The sulphur content of organic matter is not limited by the availability of porewater sulphide. Isotopic data suggest that sulphur is incorporated into organic matter within a metre of the sediment surface, at least partly by reaction with polysulphides. Most inorganic Sulphur occurs as pyrite. Pyrite formation occurred within surface sediments and was limited by the availability of reactive iron. But despite highly reducing sulphidic conditions, only 35-65% of the total iron was converted to sulphide; 10-30% of the total iron still occurs as Fe(III). In surface sediments, the isotopic composition of pyrite is similar to that of both iron monosulphide and dissolved sulphide. Either pyrite, like monosulphide, formed by direct reaction between dissolved sulphide and detrital iron, and/or the sulphur species responsible for converting FeS to FeS2 is isotopically similar to dissolved sulphide. Likely stoichiometries for the reaction between ferric iron and excess sulphide imply a maximum resulting FeS2:FeS ratio of 1:1. Where pyrite dominates the pool of iron sulphides, at least some pyrite must have formed by reaction between monosulphide and elemental sulphur and/or polysulphide. Elemental sulphur (S°) is most abundant in surface sediments and probably formed by oxidation of sulphide diffusing across the sediment-water interface. In surface sediments, S° is isotopically heavier than dissolved sulphide, FeS and FeS2 and is unlikely to have been involved in the conversion of FeS to FeS2. Polysulphides are thus implicated as the link between FeS and FeS2.

Chemical composition of sediments from DSDP Hole 149

The diagenesis and geochemical evolution of deep-sea sediments are controlled by the interaction between sediments and their associated pore waters. With increasing depth, the pore water of Hole 149 (DSDP) exhibits a strong depletion in Mg and a corresponding enrichment in Ca, while the alkalinity remains relatively constant. Dissolved SiO2 is nearly constant in the upper 100 m of sediment, but is highly enriched in the deepest pore waters. The pore waters exhibit a depletion in K with increasing depth, and O18/O16 pore water ratios also decrease. The sediment section has three zones of sedimentary regimes with increasing depth in the drill hole: an upper 100 m section of detrital clays, a middle section enriched in calc-akalic volcanics which have undergone submarine weathering to a smectite phase, and a lower section of siliceous ooze which still has a diagenetic smectite phase. The quartz-feldspar ratios and O18/O16 composition of the silicate phases are in agreement with these interpretations. The submarine weathering of volcanics to a smectite can account for the observed pore water gradients. Volcanics release Ca and Mg to the pore waters causing the alkalinity values to increase. Smectite is formed, depletes the pore waters in Mg and O18 and causes the alkalinity to decrease. The net reaction allows for the observed relationship between pore water Ca and Mg gradients with little net change in alkalinity. Given the abundance of volcanics in many deep-sea sediments, especially in lower sections which often form near ridge crests, the submarine formation of smectite may be an additional oceanic Mg sink which has not yet been fully considered.

Geochemistry of Atlantis Massif serpentinites

The Lost City hydrothermal system at the southern Atlantis Massif (Mid-Atlantic Ridge, 30°N) provides a natural laboratory for studying serpentinization processes, the temporal evolution of ultramafic-hosted hydrothermal systems, and alteration conditions during formation and emplacement of an oceanic core complex. Here we present B, O, and Sr isotope data to investigate fluid/rock interaction and mass transfer during detachment faulting and exhumation of lithospheric sequences within the Atlantis Massif. Our data indicate that extensive serpentinization was a seawater-dominated process that occurred predominately at temperatures of 150-250 °C and at high integrated W/R ratios that led to a marked boron enrichment (34-91 ppm). Boron removal from seawater during serpentinization is positively correlated with changes in d11B (11-16 per mil) but shows no correlation with O-isotope composition. Modeling indicates that B concentrations and isotope values of the serpentinites are controlled by transient temperature-pH conditions. In contrast to prior studies, we conclude that low-temperature marine weathering processes are insignificant for boron geochemistry of the Atlantis Massif serpentinites. Talc- and amphibole-rich fault rocks formed within a zone of detachment faulting at temperatures of approximately 270-350 °C and at low W/R ratios. Talc formation in ultramafic domains in the massif was subsequent to an early stage of serpentinization and was controlled by the access of Si-rich fluids derived through seawater-gabbro interactions. Replacement of serpentine by talc resulted in boron loss and significant lowering of d11B values (9-10 per mil), which we model as the product of progressive extraction of boron. Our study provides new constraints on the boron geochemical cycle at oceanic spreading ridges and suggests that serpentinization associated with ultramafic-hosted hydrothermal systems may have important implications for the behavior of boron in subduction zone settings.

Organic matter and sediment properties of samples from site GeoB15105 from the Black Sea

Dissolved organic matter (DOM) in marine sediments is a complex mixture of thousands of individual constituents that participate in biogeochemical reactions and serve as substrates for benthic microbes. Knowledge of the molecular composition of DOM is a prerequisite for a comprehensive understanding of the biogeochemical processes in sediments. In this study, interstitial water DOM was extracted with Rhizon samplers from a sediment core from the Black Sea and compared to the corresponding water-extractable organic matter fraction (<0.4 µm) obtained by Soxhlet extraction, which mobilizes labile particulate organic matter and DOM. After solid phase extraction (SPE) of DOM, samples were analyzed for the molecular composition by Fourier Transform Ion-Cyclotron Resonance Mass Spectrometry (FT-ICR MS) with electrospray ionization in negative ion mode. The average SPE extraction yield of the dissolved organic carbon (DOC) in interstitial water was 63%, whereas less than 30% of the DOC in Soxhlet-extracted organic matter was recovered. Nevertheless, Soxhlet extraction yielded up to 4.35% of the total sedimentary organic carbon, which is more than 30-times the organic carbon content of the interstitial water. While interstitial water DOM consisted primarily of carbon-, hydrogen- and oxygen-bearing compounds, Soxhlet extracts yielded more complex FT-ICR mass spectra with more peaks and higher abundances of nitrogen- and sulfur-bearing compounds. The molecular composition of both sample types was affected by the geochemical conditions in the sediment; elevated concentrations of HS- promoted the early diagenetic sulfurization of organic matter. The Soxhlet extracts from shallow sediment contained specific three- and four-nitrogen-bearing molecular formulas that were also detected in bacterial cell extracts and presumably represent proteinaceous molecules. These compounds decreased with increasing sediment depth while one- and two-nitrogen-bearing molecules increased, resulting in a higher similarity of both sample types in the deep sediment. In summary, Soxhlet extraction of sediments accessed a larger and more complex pool of organic matter than present in interstitial water DOM.

Magnesium isotope ratios of pore-fluids and dolomites

Magnesium concentrations in deep-sea sediment pore-fluids typically decrease down core due to net precipitation of dolomite or clay minerals in the sediments or underlying crust. To better characterize and differentiate these processes, we have measured magnesium isotopes in pore-fluids and sediment samples from Ocean Drilling Program sites (1082, 1086, 1012, 984, 1219, and 925) that span a range of oceanographic settings. At all sites, magnesium concentrations decrease with depth. At sites where diagenetic reactions are dominated by the respiration of organic carbon, pore-fluid d26Mg values increase with depth by as much as 2 per mil. Because carbonates preferentially incorporate 24Mg (low d26Mg), the increase in pore-fluid d26Mg values at these sites is consistent with the removal of magnesium in Mg-carbonate (dolomite). In contrast, at sites where the respiration of organic carbon is not important and/or weatherable minerals are abundant, pore-fluid d26Mg values decrease with depth by up to 2 per mil. The decline in pore-fluid d26Mg at these sites is consistent with a magnesium sink that is isotopically enriched relative to the pore-fluid. The identity of this enriched magnesium sink is likely clay minerals. Using a simple 1D diffusion-advection-reaction model of pore-fluid magnesium, we estimate rates of net magnesium uptake/removal and associated net magnesium isotope fractionation factors for sources and sinks at all sites. Independent estimates of magnesium isotope fractionation during dolomite precipitation from measured d26Mg values of dolomite samples from sites 1082 and 1012 are very similar to modeled net fractionation factors at these sites, suggesting that local exchange of magnesium between sediment and pore-fluid at these sites can be neglected. Our results indicate that the magnesium incorporated in dolomite is 2.0-2.7 per mil depleted in d26Mg relative to the precipitating fluid. Assuming local exchange of magnesium is minor at the rest of the studied sites, our results suggest that magnesium incorporated into clay minerals is enriched in d26Mg by 0 per mil to +1.25 per mil relative to the precipitating fluid. This work demonstrates the utility of magnesium isotopes as a tracer for magnesium sources/sinks in low-temperature aqueous systems.

Biogeochemistry and microbiology of groundwater during acetate and bicarbonate amendment to an alluvial aquifer

These data are from a field experiment conducted in a shallow alluvial aquifer along the Colorado River in Rifle, Colorado, USA. In this experiment, bicarbonate-promoted uranium desorption and acetate amendment were combined and compared to an acetate amendment-only experiment in the same experimental plot. Data include names and location data for boreholes, geochemical data for all the boreholes between June 1, 2010 and January 1, 2011, microarray data provided as signal to noise ratio (SNR) for individual microarray probes, microarray data provided as signal to noise ratio (SNR) by Genus.

Sulphate concentrations and sulphur isotopic composition in pore fluids from ODP Leg 168 sites

A numerical model of sulfate reduction and isotopic fractionation has been applied to pore fluid SO4**2- and d34S data from four sites drilled during Ocean Drilling Program (ODP) Leg 168 in the Cascadia Basin at 48°N, where basement temperatures reach up to 62°C. There is a source of sulfate both at the top and the bottom of the sediment column due to the presence of basement fluid flow, which promotes bacterial sulfate reduction below the sulfate minimum zone at elevated temperatures. Pore fluid d34S data show the highest values (135 per mil) yet found in the marine environment. The bacterial sulfur isotopic fractionation factor, a, is severely underestimated if the pore fluids of anoxic marine sediments are assumed to be closed systems and Rayleigh fractionation plots yield erroneous values for a by as much as 15 per mil in diffusive and advective pore fluid regimes. Model results are consistent with a = 1.077+/-0.007 with no temperature effect over the range 1.8 to 62°C and no effect of sulfate reduction rate over the range 2 to 10 pmol/ccm/day. The reason for this large isotopic fractionation is unknown, but one difference with previous studies is the very low sulfate reduction rates recorded, about two orders of magnitude lower than literature values that are in the range of µmol/ccm/day to tens of nmol/ccm/day. In general, the greatest 34S depletions are associated with the lowest sulfate reduction rates and vice versa, and it is possible that such extreme fractionation is a characteristic of open systems with low sulfate reduction rates.