Absolute abundances of benthic foraminifers and stable isotopes from the Alvin dives 3709-3712

The presence of gas hydrates on the Blake Ridge diapir, northeastern Atlantic Ocean, offers an opportunity to study the impact of methane seepage on the ecology and geochemistry of benthic foraminifera in the late Holocene. Three push cores, covering a time span of ~ 1000 yrs, were retrieved from three distinct microhabitats at the top of the diapir at a water depth of ~ 2150 m: (i) sediments away from seepage (control core), (ii) sediments overlain by clusters of methanotrophic and thiotrophic bivalves, and (iii) chemoautotrophic microbial mats. The foraminiferal assemblages at the two seep sites are marked by a reduction in benthic foraminiferal species diversity, coupled with a near-absence of agglutinated species. However, an opportunistic population rise in CH4- or H2S-tolerant calcareous species (e.g., Globocassidulina subglobosa and Cassidulina laevigata) that utilize the abundant trophic resources at the seeps has led to an increase in the overall assemblage density there. The delta18O and delta13C values of three species of benthic foraminifera - Gyroidinoides laevigatus, Globocassidulina subglobosa, and Uvigerina peregrina - and the planktonic species Globorotalia menardii were acquired from all three cores. The benthic species from methane seeps yield delta13C values of 0.1 to - 4.2 (per mil VPDB), that are distinctly more 13C-depleted relative to the delta13C of 0.4 to - 1.0 (per mil VPDB) at the control (off seep) site. The species from a mussel-bed site exhibit more negative delta13C values than those from microbial mats, possibly reflecting different food sources and higher rate of anaerobic oxidation of methane. The positive delta13C values in the paired planktonic species suggest that authigenic carbonate precipitation did not overprint the observed 13C depletions. Hence the probable cause of negative delta13C of benthic foraminifera is primary calcification from Dissolved Inorganic Carbon (DIC) containing mixed carbon fractions from (a) highly 13C-depleted, microbially-oxidized methane and (b) a seawater source.

(Table 1) Mean grain sizes and standard deviations of gas hydrate samples

Methane hydrates are present in marine seep systems and occur within the gas hydrate stability zone. Very little is known about their crystallite sizes and size distributions because they are notoriously difficult to measure. Crystal size distributions are usually considered as one of the key petrophysical parameters because they influence mechanical properties and possible compositional changes, which may occur with changing environmental conditions. Variations in grain size are relevant for gas substitution in natural hydrates by replacing CH4 with CO2 for the purpose of carbon dioxide sequestration. Here we show that crystallite sizes of gas hydrates from some locations in the Indian Ocean, Gulf of Mexico and Black Sea are in the range of 200-400 µm; larger values were obtained for deeper-buried samples from ODP Leg 204. The crystallite sizes show generally a log-normal distribution and appear to vary sometimes rapidly with location.

(Table 1) Crystal sizes of measured samples

Due to experimental difficulties grain size distributions of gas hydrate crystallites are largely unknown in natural samples. For the first time, we were able to determine grain size distributions of six natural gas hydrates for samples retrieved from the Gulf of Mexico and from Hydrate Ridge offshore Oregon from varying depths. High-energy synchrotron radiation provides high photon fluxes as well as high penetration depth and thus allows for investigation of bulk sediment samples. The gas hydrate crystallites appear to be (log-) normally distributed in the natural samples and to be of roughly globular shape. The mean grain sizes are in the range from 300-600 µm with a tendency for bigger grains to occur in greater depth, possibly indicating a difference in the formation age. Laboratory produced methane hydrate, starting from ice and aged for 3 weeks, shows half a log-normal curve with a mean value of ~40 µm. This one order-of-magnitude smaller grain sizes suggests that care must be taken when transposing grain-size sensitive (petro-)physical data from laboratory-made gas hydrates to natural settings.

(Table 1) Concentrations and carbon isotopic compositions of CH4 and CO2 at DSDP Site 76-533

The principal gaseous carbon-containing components identified in the first 400 m of sediment at Deep Sea Drilling Project Site 533, Leg 76, are methane (CH4) and carbon dioxide (CO2). Below a sub-bottom depth of about 25 m, sediment cores commonly contained pockets caused by the expansion of gas upon core recovery. The carbon isotopic composition (d13C per mil relative to PDB standard) of CH4 and CO2 in these gas pockets has been measured, resulting in the following observations: (1) d13C-CH4 values increase with depth from approximately -94 per mil in the uppermost sediment to about -66 per mil in the deepest sediment, reflecting a systematic but nonlinear depletion of 12C with depth. (2) d13C-CO2 values also increase with depth of sediment from about -25 per mil to about -4 per mil, snowing a depletion of 12C that closely parallels the trend of the isotopic composition of CH4. The magnitude and parallel distribution of d13C values for both CH4 and CO2 are consistent with the concept that the formation of the CH4 resulted from the microbiological reduction of CO2 from organic substances. These results imply that CH4 and CO2 incorporated in gas hydrates at this site are biogenic.

(Table 1) Crystal sizes of measured samples from the Bush Hill region in the Gulf of Mexico and from ODP Leg 204 at the Hydrate Ridge offshore Oregon

The grain sizes of gas hydrate crystallites are largely unknown in natural samples. Single grains are hardly detectable with electron or optical microscopy. For the first time, we have used high-energy synchrotron diffraction to determine grain sizes of six natural gas hydrates retrieved from the Bush Hill region in the Gulf of Mexico and from ODP Leg 204 at the Hydrate Ridge offshore Oregon from varying depth between 1 and 101 metres below seafloor. High-energy synchrotron radiation provides high photon fluxes as well as high penetration depth and thus allows for investigation of bulk sediment samples. Gas hydrate grain sizes were measured at the Beam Line BW 5 at the HASYLAB/Hamburg. A 'moving area detector method', originally developed for material science applications, was used to obtain both spatial and orientation information about gas hydrate grains within the sample. The gas hydrate crystal sizes appeared to be (log-)normally distributed in the natural samples. All mean grain sizes lay in the range from 300 to 600 µm with a tendency for bigger grains to occur in greater depth. Laboratory-produced methane hydrate, aged for 3 weeks, showed half a log-normal curve with a mean grain size value of c. 40 µm. The grains appeared to be globular shaped.

(Table 1) Carbon, hydrocarbons and observations at DSDP Holes 67-497 and 67-498A

Gas hydrates are icelike materials that form when specific conditions of temperature, pressure, and gas composition are simultaneously satisfied. Among the first descriptions of gas hydrates under natural conditions was that of Hammerschmidt (1940), who found them in pipelines used to transport natural gas. Milton (1976) indicates that conditions are suitable for the presence of gas hydrates in areas affected by permafrost and cites studies suggesting that large quantities of gas exist in hydrate form.

(Table 1) Deuterium ratios, chloride content and hydrate volume of pore water and hydrate meltwater from ODP Site 164-996

Site 996 is located above the Blake Diapir where numerous indications of vertical fluid migration and the presence of hydrate existed prior to Ocean Drilling Program (ODP) Leg 164. Direct sampling of hydrates and visual observations of hydrate-filled veins that could be traced 30-40 cm along cores suggest a connection between fluid migration and hydrate formation. The composition of pore water squeezed from sediment cores showed large variations due to melting of hydrate during core recovery and influence of saline water from the evaporitic diapir below. Analysis of water released during hydrate decomposition experiments showed that the recovered hydrates contained significant amounts of pore water. Solutions of the transport equations for deuterium (d2H) and chloride (Cl-) were used to determine maximum (d2H) and minimum (Cl-) in situ concentrations of these species. Minimum in situ concentrations of hydrate were estimated by combining these results with Cl- and d2H values measured on hydrate meltwaters and pore waters obtained by squeezing of sediments, by the means of a method based on analysis of distances in the two-dimensional Cl- d2H space. The computed Cl- and d2H distribution indicates that the minimum hydrate amount solutions are representative of the actual hydrate amount. The highest and mean hydrate concentrations estimates from our model are 31% and 10% of the pore space, respectively. These concentrations agree well with visual core observations, supporting the validity of the model assumptions. The minimum in situ Cl- concentrations were used to constrain the rates of upward fluid migration. Simulation of all available data gave a mean flow rate of 0.35 m/k.y. (range: 0.125-0.5 m/k.y.).

(Table 1) Ratio of hydroxy-GDGT vs. the total core GDGT was calculated based on the detection of hydroxy-GDGT

Hydroxylated glycerol dialkyl glycerol tetraethers (hydroxy-GDGTs) were detected in marine sediments of diverse depositional regimes and ages. Mass spectrometric evidence, complemented by information gleaned from two-dimensional (2D) 1H-13C nuclear magnetic resonance (NMR) spectroscopy on minute quantities of target analyte isolated from marine sediment, allowed us to identify one major compound as a monohydroxy-GDGT with acyclic biphytanyl moieties (OH-GDGT-0). NMR spectroscopic and mass spectrometric data indicate the presence of a tertiary hydroxyl group suggesting the compounds are the tetraether analogues of the widespread hydroxylated archaeol derivatives that have received great attention in geochemical studies of the last two decades. Three other related compounds were assigned as acyclic dihydroxy-GDGT (2OH-GDGT-0) and monohydroxy-GDGT with one (OH-GDGT-1) and two cyclopentane rings (OH-GDGT-2). Based on the identification of hydroxy-GDGT core lipids, a group of previously reported unknown intact polar lipids (IPLs), including the ubiquitously distributed H341-GDGT (Lipp J. S. and Hinrichs K. -U. (2009) Structural diversity and fate of intact polar lipids in marine sediments. Geochim. Cosmochim. Acta 73, 6816-6833), and its analogues were tentatively identified as glycosidic hydroxy-GDGTs. In addition to marine sediments, we also detected hydroxy-GDGTs in a culture of Methanothermococcus thermolithotrophicus. Given the previous finding of the putative polar precursor H341-GDGT in the planktonic marine crenarchaeon Nitrosopumilus maritimus, these compounds are synthesized by representatives of both cren- and euryarchaeota. The ubiquitous distribution and apparent substantial abundance of hydroxy-GDGT core lipids in marine sediments (up to 8% of total isoprenoid core GDGTs) point to their potential as proxies.

Concentrations of major and minor constituents of samples from the lower end of the cores from the eastern Mediterranean Sea (Table 2)

The pore water chemistry of mud volcanoes from the Olimpi Mud Volcano Field and the Anaximander Mountains in the eastern Mediterranean Sea have been studied for three major purposes: (1) modes and velocities of fluid transport were derived to assess the role of (upward) advection, and bioirrigation for benthic fluxes. (2) Differences in the fluid chemistry at sites of Milano mud volcano (Olimpi area) were compiled in a map to illustrate the spatial heterogeneity reflecting differences in fluid origin and transport in discrete conduits in near proximity. (3) Formation water temperatures of seeping fluids were calculated from theoretical geothermometers to predict the depth of fluid origin and geochemical reactions in the deeper subsurface. No indications for downward advection as required for convection cells have been found. Instead, measured pore water profiles have been simulated successfully by accounting for upward advection and bioirrigation. Advective flow velocities are found to be generally moderate (3-50 cm/y) compared to other cold seep areas. Depth-integrated rates of bioirrigation are 1-2 orders of magnitude higher than advective flow velocities documenting the importance of bioirrigation for flux considerations in surface sediments. Calculated formation water temperatures from the Anaximander Mountains are in the range of 80 to 145 °C suggesting a fluid origin from a depth zone associated with the seismic decollement. It is proposed that at that depth clay mineral dehydration leads to the formation and advection of fluids reduced in salinity relative to sea water. This explains the ubiquitous pore water freshening observed in surface sediments of the Anaximander Mountain area. Multiple fluid sources and formation water temperatures of 55 to 80 °C were derived for expelled fluids of the Olimpi area.

Gas analysis of sediment cores from the Batumi seep area

Detailed knowledge of the extent of post-genetic modifications affecting shallow submarine hydrocarbons fueled from the deep subsurface is fundamental for evaluating source and reservoir properties. We investigated gases from a submarine high-flux seepage site in the anoxic Eastern Black Sea in order to elucidate molecular and isotopic alterations of low-molecular-weight hydrocarbons (LMWHC) associated with upward migration through the sediment and precipitation of shallow gas hydrates. For this, near-surface sediment pressure cores and free gas venting from the seafloor were collected using autoclave technology at the Batumi seep area at 845 m water depth within the gas hydrate stability zone. Vent gas, gas from pressure core degassing, and from hydrate dissociation were strongly dominated by methane (>99.85 mol.% of Sum[C1-C4, CO2]). Molecular ratios of LMWHC (C1/[C2 + C3] > 1000) and stable isotopic compositions of methane (d13C = -53.5 per mill V-PDB; D/H around -175 per mill SMOW) indicated predominant microbial methane formation. C1/C2+ ratios and stable isotopic compositions of LMWHC distinguished three gas types prevailing in the seepage area. Vent gas discharged into bottom waters was depleted in methane by >0.03 mol.% (Sum[C1-C4, CO2]) relative to the other gas types and the virtual lack of 14C-CH4 indicated a negligible input of methane from degradation of fresh organic matter. Of all gas types analyzed, vent gas was least affected by molecular fractionation, thus, its origin from the deep subsurface rather than from decomposing hydrates in near-surface sediments is likely. As a result of the anaerobic oxidation of methane, LMWHC in pressure cores in top sediments included smaller methane fractions [0.03 mol.% Sum(C1-C4, CO2)] than gas released from pressure cores of more deeply buried sediments, where the fraction of methane was maximal due to its preferential incorporation in hydrate lattices. No indications for stable carbon isotopic fractionations of methane during hydrate crystallization from vent gas were found. Enrichments of 14C-CH4 (1.4 pMC) in short cores relative to lower abundances (max. 0.6 pMC) in gas from long cores and gas hydrates substantiates recent methanogenesis utilizing modern organic matter deposited in top sediments of this high-flux hydrocarbon seep area.

Pore-water concentration of chloride in sediment cores and equilibrium temperatures and gas bubble analysis at ROV stations from the Vodyanitskii mud volcano

Vodyanitskii mud volcano is located at a depth of about 2070 m in the Sorokin Trough, Black sea. It is a 500-m wide and 20-m high cone surrounded by a depression, which is typical of many mud volcanoes in the Black Sea. 75 kHz sidescan sonar show different generations of mud flows that include mud breccia, authigenic carbonates, and gas hydrates that were sampled by gravity coring. The fluids that flow through or erupt with the mud are enriched in chloride (up to 650 mmol L**-1 at 150-cm sediment depth) suggesting a deep source, which is similar to the fluids of the close-by Dvurechenskii mud volcano. Direct observation with the remotely operated vehicle Quest revealed gas bubbles emanating at two distinct sites at the crest of the mud volcano, which confirms earlier observations of bubble-induced hydroacoustic anomalies in echosounder records. The sediments at the main bubble emission site show a thermal anomaly with temperatures at 60 cm sediment depth that were 0.9 °C warmer than the bottom water. Chemical and isotopic analyses of the emanated gas revealed that it consisted primarily of methane (99.8%) and was of microbial origin (dD-CH4 = -170.8 per mil (SMOW), d13C-CH4 = -61.0 per mil (V-PDB), d13C-C2H6 = -44.0 per mil (V-PDB)). The gas flux was estimated using the video observations of the ROV. Assuming that the flux is constant with time, about 0.9 ± 0.5 x 10**6 mol of methane is released every year. This value is of the same order-of-magnitude as reported fluxes of dissolved methane released with pore water at other mud volcanoes. This suggests that bubble emanation is a significant pathway transporting methane from the sediments into the water column.

Ex situ sulphate reduction and methane oxidation rates of sediment at the rim of a vesicomyd clam colony in the Japan Deep Sea Trench (dive 956)

Sediment samples were collected from the rim of a large vesicomyid clam colony in the Japan Deep Sea Trench. Immediately after sample recovery onboard, the sediment core was sub-sampled for ex situ rate measurements. Sulfate reduction and anaerobic oxidation of methane were measured ex situ by the whole core injection method with three replicate measurements for each method. We incubated the samples at in situ temperature (1.5°C) for 48 hours with either 14C-methane (dissolved in water, 2.5 kBq) or carrier-free 35S-sulfate (dissolved in water, 50 kBq). Sediment was fixed in 25 ml sodium hydroxide (NaOH) solution (2.5%, w/v) or 20 ml ZnAc solution (20%, w/v) for AOM or SR, respectively. Turnover rates were measured as previously described (Kallmeyer et al., 2004; Treude et al., 2003).