Geochemical relation between gas hydrates and water venting in the seismic blanking zone on the northern Cascadia continental margin offshore Vancouver Island, Canada

The isotopic composition (δD and δ¹⁸O) and chloride concentration (Cl⁻) of pore waters from the northern Cascadia continental margin offshore Vancouver Island were measured to characterize the relations between the water flow regime and the distribution, formation and dissociation of gas hydrates. The δD values of pore waters in gas hydrate-bearing sediments are slightly higher ( 1‰) than those of seawater as the result of gas hydrate dissociation during core recovery and handling. Within the seismic blanking zone, the δD values were slightly lower (− 1‰) than values measured from sites outside the blanking area (0‰). We attribute these differences to 1) distillation of D-rich water during hydrate formation in the center of the blanking zone and 2) limited migration of pore water between inside and outside of the blanking zone due to different fluid fluxes. In contrast, the δ¹⁸O values and Cl⁻ concentrations do not show significant spatial variation due to decreased isotopic fractionation of oxygen and small fraction of chloride relative to hydrogen isotope during gas hydrate formation. The δD value of pore water, therefore, appears to be a sensitive indicator of gas hydrate occurrence. We estimate that gas hydrate occupied at least 2.0 to 6.3% of sediment pore space using δD distribution in this area.     

Fluid flow, methane fluxes, carbonate precipitation and biogeochemical turnover in gas hydrate-bearing sediments at Hydrate Ridge, Cascadia Margin: numerical modeling and mass balances

A numerical model was applied to investigate and to quantify biogeochemical processes and methane turnover in gas hydrate-bearing surface sediments from a cold vent site situated at Hydrate Ridge, an accretionary structure located in the Cascadia Margin subduction zone. Steady state simulations were carried out to obtain a comprehensive overview on the activity in these sediments which are covered with bacterial mats and are affected by strong fluid flow from below. The model results underline the dominance of advective fluid flow that forces a large inflow of methane from below (869 μmol cm⁻² a⁻¹) inducing high oxidation rates in the surface layers. Anaerobic methane oxidation is the major process, proceeding at a depth-integrated rate of 870 μmol cm⁻² a⁻¹. A significant fraction (14%) of bicarbonate produced by anaerobic methane oxidation is removed from the fluids by precipitation of authigenic aragonite and calcite. The total rate of carbonate precipitation (120 μmol cm⁻² a⁻¹) allows for the build-up of a massive carbonate layer with a thickness of 1 m over a period of 20,000 years. Aragonite is the major carbonate mineral formed by anaerobic methane oxidation if the flow velocity of methane-charge fluids is high enough (≥10 cm a⁻¹) to maintain super-saturation with respect to this highly soluble carbonate phase. It precipitates much faster within the studied surface sediments than previously observed in abiotic laboratory experiments, suggesting microbial catalysis. The investigated station is characterized by high carbon and oxygen turnover rates (≈1000 μmol cm⁻² a⁻¹) that are well beyond the rates observed at other continental slope sites not affected by fluid venting. This underlines the strong impact of fluid venting on the benthic system, even though the flow velocity of 10 cm a⁻¹ derived by the model is relative low compared to fluid flow rates found at other cold vent sites. Non-steady state simulations using measured fluid flow velocities as forcing demonstrate a rapid respond of the sediments within a few days to changes in advective flow. Moreover, they reveal that efficient methane oxidation in these sediments prevents methane outflow into the bottom water over a wide range of fluid flow velocities (<80 cm a⁻¹). Only at flow rates exceeding approximately 100 cm a⁻¹, does dissolved methane break through the sediment surface to induce large fluxes of up to 5000 μmol CH₄ cm² a⁻¹ into the overlying bottom water.     

The methane hydrate formation and the resource estimate resulting from free gas migration in seeping seafloor hydrate stability zone

It is a typical multiphase flow process for hydrate formation in seeping seafloor sediments. Free gas can not only be present but also take part in formation of hydrate. The volume fraction of free gas in local pore of hydrate stable zone (HSZ) influences the formation of hydrate in seeping seafloor area, and methane flux determines the abundance and resource of hydrate-bearing reservoirs. In this paper, a multiphase flow model including water (dissolved methane and salt)-free gas hydrate has been established to describe this kind of flow-transfer-reaction process where there exists a large scale of free gas migration and transform in seafloor pore. In the order of three different scenarios, the conversions among permeability, capillary pressure, phase saturations and salinity along with the formation of hydrate have been deducted. Furthermore, the influence of four sorts of free gas saturations and three classes of methane fluxes on hydrate formation and the resource has also been analyzed and compared. Based on the rules drawn from the simulation, and combined information gotten from drills in field, the methane hydrate(MH) formation in Shenhu area of South China Sea has been forecasted. It has been speculated that there may breed a moderate methane flux below this seafloor HSZ. If the flux is about 0.5 kg m⁻² a⁻¹, then it will go on to evolve about 2700 ka until the hydrate saturation in pore will arrive its peak (about 75%). Approximately 1.47 × 10⁹ m³ MH has been reckoned in this marine basin finally, is about 13 times over preliminary estimate.     

Geochemical characteristics of sediment pore water from Site XS-01 in the Xisha trough of South China Sea and their significance for gas hydrate occurrence

Gas hydrate is a recently-found new source of energy that mostly exists in marine sediments. In recent years, we have conducted gas hydrate exploration in the South China Sea. The Xisha trough, one of the promising target areas for gas hydrate, is located in the northern margin of the South China Sea, adjacent to several large oil and gas fields. The Xisha trough extends 420 km long with the water depth of 1 500 m in the west part and 3 400 m in the east part and deposits thick sediments with organic matter content of 0.41%–1.02%. Previous studies on topographical features, geological P-T conditions, structural geology, sedimentary geology and geophysical bottom simulating reflectors (BSR) in the Xisha trough suggest that this area is favorable for the formation and accumulation of gas hydrate. In this paper, we present geochemical analyses for the sediment and pore water from a piston core at Site XS-01 in the Xisha trough. Seven pore water samples were analyzed for their anion (Cl⁻, SO₄ ²⁻, Br⁻, I⁻) contents, cation (Na, K, Ca, Mg) contents and trace element (Li, B, Sr, Ba, Rb, Mn) contents. Eight sediment samples were analyzed for stable carbon and oxygen isotopic compositions. A number of geochemical anomalies such as anions (e.g. Cl⁻, SO₄ ²⁻), cations (e.g. Ca, Mg) and trace elements (e.g. Sr, Ba, B) were found in this study. For example, the concentrations of Cl⁻ and SO₄ ²⁻ in pore water show a decreasing trend with depth. The estimated sulfate/methane interface (SMI) is only 18 m, which is quite similar to the SMI value of 23 m in the ODP164 Leg 997 at Blake Ridge. The Ca, Mg and Sr concentrations of pore water also decrease with depth, but concentrations of Ba, and Mg/Ca and Sr/Ca ratios increase with depth. These geochemical anomalies are quite similar to those found in gas hydrate locations in the world such as the Blake Ridge and may be related to the formation and dissociation of gas hydrates. The salt exclusion effect during the gas hydrate formation will cause an increase in major ion concentrations in the pore waters that diffused upward such as Cl. The anaerobic methane oxidation (AMO) may lead to the change of SO₄ ²⁻ and other cations such as Ca, Mg, Sr and Ba in pore water. Low δ ¹³C value of authigenic carbonates is a good indicator for gas hydrate occurrence. However, the bulk sediment samples we analyzed all show normal δ ¹³C values similar to biogenic marine carbonates, and this may also suggest that no gas hydrate-related authigenic carbonates exist or their amount is so small that they are not detectable by using this bulk analytical method. In conclusion, we suggest that the Site XS-01 in the Xisha trough of the northern margin of the South China Sea is a potential target for further gas hydrate exploration. Translated from Quaternary Sciences, 2006, 26(3): 442–448 [译自: 第四纪研究]     

Thermodynamic and pore water halogen constraints on gas hydrate distribution at ODP Site 997 (Blake Ridge)

Marine sediment sequences with CH₄ hydrate are characterized by an atypical depth profile in dissolved Cl⁻ squeezed from pore space: a shallow subsurface Cl⁻ maximum overlies a lengthy and pronounced Cl⁻ minimum. This pore water Cl⁻ profile represents a combination of multiple processes including glacial–interglacial variations in ocean salinity, advection and diffusion of ions that are excluded during gas hydrate formation at depth, and release of fresh water from dissociation of hydrate during core recovery. In situ quantities of gas hydrate can be determined from a measured pore water Cl⁻ profile provided the in situ pore water signature prior to core recovery can be separated. Ocean Drilling Program (ODP) Site 997 was drilled into a large CH₄ hydrate reservoir on the Blake Ridge in the western Atlantic Ocean. Previously we have constructed a high-resolution pore water Cl⁻ profile at this location; here we present a `coupled chloride-hydrate' numerical model to explain basic trends in the Cl⁻ profile and to isolate in situ Cl⁻ concentrations. The model is based on thermodynamic and ecological considerations, and uses established equations for describing chemical behavior in marine sediment–pore water systems. The model incorporates four key concepts: (1) most gas hydrate is formed immediately below the SO₄²⁻ reduction zone; (2) fluid, dissolved ions and gas advect upward through the sediment column; (3) CH₄ hydrate dissociates at the base of hydrate stability conditions; and (4) seawater salinity fluctuates during glacial–interglacial cycles of the late Pliocene and Quaternary. Rates of upward advection in the model are sufficient to account for measured Br⁻ and I⁻ concentrations as well as CH₄ oxidation at the base of the SO₄²⁻ reduction zone. In situ pore water Cl⁻ inferred from the model is similar to that determined by limited direct sampling; in situ CH₄ hydrate amounts inferred from the model (an average of about 4% of porosity) are broadly consistent with those determined by direct gas sampling and indirect geophysical techniques. The model also predicts production of substantial quantities of free CH₄ gas bubbles (>2.5% of porosity) at a depth immediately below the lowest accumulation of CH₄ hydrate in the sediment column. Our explanation for the pore water Cl⁻ profile at Site 997 is important because it provides a theoretical mechanism for understanding the distribution of interstitial water Cl⁻, gas hydrate, and free gas in a marine sediment column.     

甲烷水合物生成过程中海水常量离子浓度的变化规律

本文自行研制了一套甲烷水合物合成装置,模拟海洋环境甲烷水合物的生成过程,对该过程水合物生成位置、形态、反应时间、环境温压条件进行观测,同时连续测试体系海水中常量离子K+、Na+、Ca2+、Mg2+、C1-、SO42-的浓度及海水盐度,探讨水合物生成过程的温压变化及离子浓度变化之间的关系和离子浓度的变化规律.结果表明,海水中甲烷水合物生成具有很大的随机性,在相同的初始条件下可能有不同的水合物成核、聚集过程;甲烷水合物在生成过程中,耗气量不断增加,孔隙水的盐度和海水中常量阴阳离子的浓度也在不断增加,这种变化具有较高的线性相关性(相关系数为0.9848～0.9950),且不受甲烷水合物生成位置及状态的影响;在水合物生成过程的微环境下耗气量相同时,离子浓度存在细微的差异.这些特征为通过测定海底水合物周围孔隙水中常量离子的浓度初步推算水合物的甲烷耗气量提供了依据.     

盐类对甲烷水合物稳定性影响研究进展

甲烷水合物稳定性主要控制着甲烷水合物稳定带的厚度,温度、压力、孔隙水盐度和气体组分等因素影响着水合物稳定带的厚度。甲烷水合物的形成与地层水关系密切,而地层水中的各种盐离子(Cl~-、Na~+、Mg~(2+)、SO_4(~2-)、Ca~(2+))以及过渡金属(Fe、Mn、Cu、Co、Ni等)会影响天然气水合物的形成和分解条件。因此,研究盐类对甲烷水合物的稳定性认识有助于更加深入了解天然气水合物的成藏条件。本文分析了氯化物、硫酸盐、碳酸盐三大盐类对甲烷水合物稳定性的影响:同一盐类不同盐度条件下,随着盐度的增加,甲烷水合物相平衡曲线向低温高压偏移。总结了不同盐类和阴阳离子对甲烷水合物的抑制作用大小:在相同浓度、不同盐类条件下,盐类浓度在1.0~1.5 mol/L时盐类对甲烷水合物的抑制作用大小为MgCl_2CaCl_2Na ClKCl,盐类浓度大于1.5 mol/L时CaCl_2的抑制作用较强;阴离子对甲烷水合物的抑制作用大小争议较大,阳离子中Mg~(2+)对甲烷水合物的抑制作用最强。从目前的研究成果来看,已有数据与实际地质条件还存在一定差距,需要在真实实验条件下加强氯化物-硫酸盐-碳酸盐-甲烷-水体系的详细研究。本文提出,将高压可视反应腔与显微激光拉曼技术相结合,有望准确获取天然气水合物稳定形成时的温压条件,明确盐类和阴阳离子的抑制作用大小,以及盐类和离子特性如何影响水合物的形成和稳定,以便为未来的水合物勘探开发提供参考。     

Geochemical studies of pore fluid in surface sediment on the Daini Atsumi Knoll

We collected sediment samples and pore water samples from the surface sediment on the Daini Atsumi Knoll, and analyzed the sediments for CH₄, C₂H₆, and δ¹³C_CH4, and the pore fluids for CH₄, C₂H₆, δ¹³C_CH4, Cl⁻, SO₄²⁻, δ¹⁸O_H2O, and δD_H2O, respectively. A comparison of the measured concentration and isotopic composition of methane in pore water samples with those in sediment samples revealed that methane was present in the sediment samples at a higher concentration and was isotopically heavier than those in the pore water samples. It suggests that the effect of the release of a sorbed gas bound to organic particles when heated prior to analysis of hydrocarbons was larger than that of the degassing process. A large amount of a sorbed gas would be a significant source of natural gas. Two striking features are the chemical and isotopic composition of the pore water samples taken from the different sites around the Daini Atsumi Knoll. In the KL09, KL10, and KP07 samples, Cl⁻ concentrations in the pore water samples showed depletion to a minimum of 460 mmol/kg, correspond to  17% dilution of seawater, however the latter was not enriched in CH₄. The isotopic compositions of pore water samples suggested the low-Cl⁻ fluids in the pore water were not derived from dissociation of methane hydrate, but were derived from input of meteoric water. In contrast, in the KP05 samples from the north flank of the Daini Atsumi Knoll, pore water were characterized by CH₄ enrichment more than 370 μmol/kg, but not depleted in Cl⁻ concentrations. The observed methane concentration in the KP05 samples is not sufficient for methane hydrate to form in situ, indicating that the existence of methane hydrate in the surface sediment is negligible, as supported by Cl⁻ concentration. Based on the stable carbon isotope ratio of methane in the pore fluid from the KP05 site (δ¹³C_CH4 < − 50‰PDB), methane is thought to be of microbial origin. The pore waters in the surface sediments in the north flank of the Daini Atsumi Knoll were not directly influenced by upward fluid bearing methane of thermogenic origin from a deeper part of the sedimentary layer. However, extremely high methane concentration in the north flank site as compared with the concentration of pore water taken from the normal seafloor suggests that the north flank site is not the normal seafloor. We hypothesize that upward migration of chemically-reduced fluids from a deeper zone of the sedimentary layer reduces chemically-oxidized solutes in the surface sediment. As a consequence methane production replaced sulfate reduction as the microbial metabolism in the reduced environment of the surface sediment.     

Spatial variation in fluid flow and geochemical fluxes across the sediment-seawater interface at the Carlos Ribeiro mud volcano (Gulf of Cadiz)

Submarine mud volcanism is an important pathway for transfer of deep-sourced fluids enriched in hydrocarbons and other elements into the ocean. Numerous mud volcanoes (MVs) have been discovered along oceanic plate margins, and integrated elemental fluxes are potentially significant for oceanic chemical budgets. Here, we present the first detailed study of the spatial variation in fluid and chemical fluxes at the Carlos Ribeiro MV in the Gulf of Cadiz. To this end, we combine analyses of the chemical composition of pore fluids with a 1-D transport-reaction model to quantify fluid fluxes, and fluxes of boron, lithium and methane, across the sediment-seawater interface. The pore fluids are significantly depleted in chloride, but enriched in lithium, boron and hydrocarbons, relative to seawater. Pore water profiles of sulphate, hydrogen sulphide and total alkalinity indicate that anaerobic oxidation of methane occurs at 34-180 cm depth below seafloor. Clay mineral dehydration, and in particular the transformation of smectite to illite, produces pore fluids that are depleted in chloride and potassium. Profiles of boron, lithium and potassium are closely related, which suggests that lithium and boron are released from the sediments during this transformation. Pore fluids are expelled into the water column by advection; fluid flow velocities are 4 cm yr⁻¹ at the apex of the MV but they rapidly decrease to 0.4 cm yr⁻¹ at the periphery. The associated fluxes of boron, lithium and methane vary between 7-301, 0.5-6 and 0-806 mmol m⁻² yr⁻¹, respectively. We demonstrate that fluxes of Li and B due to mud volcanism may be important on a global scale, however, release of methane into the overlying water column is suppressed by microbial methanotrophy.     

The effect of dissolved barium on biogeochemical processes at cold seeps

A numerical model was applied to investigate and quantify the biogeochemical processes fueled by the expulsion of barium and methane-rich fluids in the sediments of a giant cold-seep area in the Derugin Basin (Sea of Okhotsk). Geochemical profiles of dissolved Ba²⁺, Sr²⁺, Ca²⁺, SO₄²⁻, HS⁻, DIC, I⁻ and of calcium carbonate (CaCO₃) were fitted numerically to constrain the transport processes and the kinetics of biogeochemical reactions. The model results indicate that the anaerobic oxidation of methane (AOM) is the major process proceeding at a depth-integrated rate of 4.9 μmol cm⁻² a⁻¹, followed by calcium carbonate and strontian barite precipitation/dissolution processes having a total depth-integrated rate of 2.1 μmol cm⁻² a⁻¹. At the low seepage rate prevailing at our study site (0.14 cm a⁻¹) all of the rising barium is consumed by precipitation of barite in the sedimentary column and no benthic barium flux is produced. Numerical experiments were run to investigate the response of this diagenetic environment to variations of hydrological and biogeochemical conditions. Our results show that relatively low rates of fluid flow (<∼5 cm a⁻¹) promote the dispersed precipitation of up to 26 wt% of barite and calcium carbonate throughout the uppermost few meters of the sedimentary column. Distinct and persistent events (several hundreds of years long) of more vigorous fluid flow (from 20-110 cm a⁻¹), instead, result in the formation of barite-carbonate crusts near the sediment surface. Competition between barium and methane for sulfate controls the mineralogy of these sediment precipitates such that at low dissolved methane/barium ratios (<4-11) barite precipitation dominates, while at higher methane/barium ratios sulfate availability is limited by AOM and calcium carbonate prevails. When seepage rates exceed 110 cm a⁻¹, barite precipitation occurs at the seafloor and is so rapid that barite chimneys form in the water column. In the Derugin Basin, spectacular barite constructions up to 20 m high, which cover an area of roughly 22 km² and contain in excess of 5 million tons of barite, are built through this process. In these conditions, our model calculates a flux of barium to the water column of at least 20 μmol cm⁻² a⁻¹. We estimate that a minimum of 0.44 × 10⁶ mol a⁻¹ are added to the bottom waters of the Derugin Basin by cold seep processes, likely affecting the barium cycle in the Sea of Okhotsk.     

Pleistocene recharge to midcontinent basins: effects on salinity structure and microbial gas generation

The hydrogeochemistry of saline-meteoric water interface zones in sedimentary basins is important in constraining the fluid migration history, chemical evolution of basinal brines, and physical stability of saline formation waters during episodes of freshwater recharge. This is especially germane for interior cratonic basins, such as the Michigan and Illinois basins. Although there are large differences in formation water salinity and hydrostratigraphy in these basins, both are relatively quiescent tectonically and have experienced repeated cycles of glaciation during the Pleistocene. Exploration for unconventional microbial gas deposits, which began in the upper Devonian-age Antrim Shale at the northern margin of the Michigan Basin, has recently extended into the age-equivalent New Albany Shale of the neighboring Illinois Basin, providing access to heretofore unavailable fluid samples. These reveal an extensive regional recharge system that has profoundly changed the salinity structure and induced significant biogeochemical modification of formation water elemental and isotope geochemistry.New-formation water and gas samples were obtained from Devonian-Mississippian strata in the Illinois Basin. These included exploration wells in the New Albany Shale, an organic-rich black shale of upper Devonian age, and formation waters from over- and underlying regional aquifer systems (Siluro-Devonian and Mississippian age). The hydrostratigraphic relations of major aquifers and aquitards along the eastern margin of the Illinois Basin critically influenced fluid migration into the New Albany Shale. The New Albany Shale formation water chemistry indicates significant invasion of meteoric water, with δD values as low as −46.05‰, into the shale. The carbon stable isotope system (δ¹³C values as high as 29.4‰), coupled with δ¹⁸O, δD, and alkalinity of formation waters (alkalinity ≤24.08 meq/kg), identifies the presence of microbial gas associated with meteoric recharge. Regional geochemical patterns identify the underlying Siluro-Devonian carbonate aquifer system as the major conduit for freshwater recharge into the fractured New Albany Shale reservoirs. Recharge from overlying Mississippian carbonates is only significant in the southernmost portion of the basin margin where carbonates directly overlie the New Albany Shale.Recharge of dilute waters (Cl⁻ <1000 mM) into the Siluro-Devonian section has suppressed formation water salinity to depths as great as 1 km across the entire eastern Illinois Basin margin. Taken together with salinity and stable isotope patterns in age-equivalent Michigan Basin formation waters, they suggest a regional impact of recharge of δ¹⁸O- and δD-depleted fluids related to Pleistocene glaciation. Devonian black shales at both basin margins have been affected by recharge and produced significant volumes of microbial methane. This recharge is also manifested in different salinity gradients in the two basins because of their large differences in original formation water salinity. Given the relatively quiet tectonic history and subdued current topography in the midcontinent region, it is likely that repeated cycles of glacial meltwater invasion across this region have induced a strong disequilibrium pattern in fluid salinity and produced a unique class of unconventional shale-hosted gas deposits.     

南海北部陆坡西沙海槽XS-01站位沉积物孔隙水的 地球化学特征及其对天然气水合物的指示意义*

天然气水合物是近年来国际上发现的一种新型能源,大量赋存在海底沉积物中。西沙海槽位于南海北部陆坡区,周边有多个大型深水油气田区。对该区地形地貌、地质构造和沉积条件分析以及地球物理BSR分布表明,西沙海槽是我国海洋天然气水合物资源勘查的一个有利远景区。文章主要研究了位于西沙海槽最大BSR区内的XS-01站位沉积物孔隙水的地球化学特征,发现该站位孔隙水阴阳离子浓度和微量元素组成特征变化显示出可能与天然气水合物有关的明显地球化学异常,与国际上己发现有天然气水合物地区的异常相类似。因此,认为该站位是西沙海槽区最有利的天然气水合物赋存区,值得进一步的勘查工作。     

A review of the geochemistry of methane in natural gas hydrate

The largest accumulations on Earth of natural gas are in the form of gas hydrate, found mainly offshore in outer continental margin sediment and, to a lesser extent, in polar regions commonly associated with permafrost. Measurements of hydrocarbon gas compositions and of carbon-isotopic compositions of methane from natural gas hydrate samples, collected in subaquatic settings from around the world, suggest that methane guest molecules in the water clathrate structures are mainly derived by the microbial reduction of CO₂ from sedimentary organic matter. Typically, these hydrocarbon gases are composed of > 99% methane, with carbon-isotopic compositions (δ¹³C_PDB) ranging from − 57 to − 73‰. In only two regions, the Gulf of Mexico and the Caspian Sea, has mainly thermogenic methane been found in gas hydrate. There, hydrocarbon gases have methane contents ranging from 21 to 97%, with δ¹³C values ranging from − 29 to − 57‰. At a few locations, where the gas hydrate contains a mixture of microbial and thermal methane, microbial methane is always dominant. Continental gas hydrate, identified in Alaska and Russia, also has hydrocarbon gases composed of > 99% methane, with carbon-isotopic compositions ranging from − 41 to − 49‰. These gas hydrate deposits also contain a mixture of microbial and thermal methane, with thermal methane likely to be dominant. Published by Elsevier Science Ltd     

A cryptic sulfur cycle driven by iron in the methane zone of marine sediment (Aarhus Bay, Denmark)

Sulfate reduction and sulfur-iron geochemistry were studied in 5-6 m deep gravity cores of Holocene mud from Aarhus Bay (Denmark). A goal was to understand whether sulfate is generated by re-oxidation of sulfide throughout the sulfate and methane zones, which might explain the abundance of active sulfate reducers deep below the main sulfate zone. Sulfate penetrated down to 130 cm where methane started to build up and where the concentration of free sulfide peaked at 5.5 mM. Below this sulfate-methane transition, sulfide diffused downwards to a sulfidization front at 520 cm depth, below which dissolved iron, Fe²⁺, accumulated in the pore water. Sulfate reduction rates measured by ³⁵S-tracer incubations in the sulfate zone were high due to high concentrations of reactive organic matter. Within the sulfate-methane transition, sulfate reduction was distinctly stimulated by the anaerobic oxidation of methane. In the methane zone below, sulfate remained at positive “background” concentrations of <0.5 mM down to the sulfidization front. Sulfate reduction decreased steeply to rates which at 300-500 cm depth were 0.2-1 pmol SO₄²⁻ cm⁻³ d⁻¹, i.e., 4-5 orders of magnitude lower than rates measured near the sediment surface. The turn-over time of sulfate increased from 3 years at 12 cm depth to 100-1000 years down in the methane zone. Sulfate reduction in the methane zone accounted for only 0.1% of sulfate reduction in the entire sediment column and was apparently limited by the low pore water concentration of sulfate and the low availability of organic substrates. Amendment of the sediment with both sulfate and organic substrates immediately caused a 10- to 40-fold higher, “potential sulfate reduction” which showed that a physiologically intact community of sulfate reducing bacteria was present. The “background” sulfate concentration appears to be generated from the reaction of downwards diffusing sulfide with deeply buried Fe(III) species, such as poorly-reactive iron oxides or iron bound in reactive silicates. The oxidation of sulfide to sulfate in the sulfidic sediment may involve the formation of elemental sulfur and thiosulfate and their further disproportionation to sulfide and sulfate. The net reaction of sulfide and Fe(III) to form pyrite requires an additional oxidant, irrespective of the formation of sulfate. This could be CO₂ which is reduced with H₂ to methane. The methane subsequently diffuses upwards to become re-oxidized at the sulfate-methane transition and thereby removes excess reducing power and enables the formation of excess sulfate. We show here how the combination of these well-established sulfur-iron-carbon reactions may lead to the deep formation of sulfate and drive a cryptic sulfur cycle. The iron-rich post-glacial sediments underlying Holocene marine mud stimulate the strong sub-surface sulfide reoxidation observed in Aarhus Bay and are a result of the glacial to interglacial history of the Baltic Sea area. Yet, processes similar to the ones described here probably occur widespread in marine sediments, in particular along the ocean margins.     

Origin and age of pore waters in an actively venting gas hydrate field near Sado Island,Japan Sea: Interpretation of halogen and 129I distributions

A gas hydrate field with highly active venting of methane was recently found near Sado Island in the eastern Japan Sea. Piston cores were collected from active venting sites and nearby locations in the Umitaka Spur–Joetsu Knoll area during two cruises in 2004 (UT04) and 2005 (KY05-08). We report here halogen concentrations and ¹²⁹I/I ratios in pore waters associated with gas hydrates from these expeditions. The strongly biophilic behavior of I and, to a lesser degree, of Br together with the presence of the long-lived iodine radioisotope (¹²⁹I) allow evaluation of potential source materials for methane in gas hydrate systems. Depth profiles of all three halogens, particularly the very rapid downward increases of Br and I concentrations, strongly suggest input of deep fluids enriched in Br and I, but the profiles also display the effects of gas hydrate formation and dissociation. Although the ¹²⁹I/I ratios are modified by ¹²⁹I from seawater and sediments at shallow depth, likely ratios of the deep fluids are estimated to be between 400 × 10^− 15 and 600 × 10^− 15, equivalent to a Late Oligocene to Early Miocene age. Ages in the active methane venting sites typically are closer to the old end of this range than those in the reference sites. This age range suggests that the methane associated with venting and gas hydrate formation in this area is derived from organic materials accumulated during the initial opening of the Japan Sea. The Umitaka Spur–Joetsu Knoll gas hydrate field demonstrates the movement of deep fluids associated with the release of significant amounts of methane from the seafloor, processes which might be important components of mass transfer and carbon cycle in the shallow geosphere.     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(II) and Fe(III)

Equilibrium and kinetic Fe isotope fractionation between aqueous ferrous and ferric species measured over a range of chloride concentrations (0, 11, 110 mM Cl⁻) and at two temperatures (0 and 22°C) indicate that Fe isotope fractionation is a function of temperature, but independent of chloride contents over the range studied. Using ⁵⁷Fe-enriched tracer experiments the kinetics of isotopic exchange can be fit by a second-order rate equation, or a first-order equation with respect to both ferrous and ferric iron. The exchange is rapid at 22°C, ∼60-80% complete within 5 seconds, whereas at 0°C, exchange rates are about an order of magnitude slower. Isotopic exchange rates vary with chloride contents, where ferrous-ferric isotope exchange rates were ∼25 to 40% slower in the 11 mM HCl solution compared to the 0 mM Cl⁻ (∼10 mM HNO₃) solutions; isotope exchange rates are comparable in the 0 and 110 mM Cl⁻ solutions.The average measured equilibrium isotope fractionations, Δ_{Fe(III)-Fe(II)}, in 0, 11, and 111 mM Cl⁻ solutions at 22°C are identical within experimental error at +2.76±0.09, +2.87±0.22, and +2.76±0.06 ‰, respectively. This is very similar to the value measured by Johnson et al. (2002a) in dilute HCl solutions. At 0°C, the average measured Δ_{Fe(III)-Fe(II)} fractionations are +3.25±0.38, +3.51±0.14 and +3.56±0.16 ‰ for 0, 11, and 111 mM Cl⁻ solutions. Assessment of the effects of partial re-equilibration on isotope fractionation during species separation suggests that the measured isotope fractionations are on average too low by ∼0.20 ‰ and ∼0.13 ‰ for the 22°C and 0°C experiments, respectively. Using corrected fractionation factors, we can define the temperature dependence of the isotope fractionation from 0°C to 22°C as: where the isotopic fractionation is independent of Cl⁻ contents over the range used in these experiments. These results confirm that the Fe(III)-Fe(II) fractionation is approximately half that predicted from spectroscopic data, and suggests that, at least in moderate Cl⁻ contents, the isotopic fractionation is relatively insensitive to Fe-Cl speciation.     

Microbial methane turnover at mud volcanoes of the Gulf of Cadiz

The Gulf of Cadiz is a tectonically active area of the European continental margin and characterised by a high abundance of mud volcanoes, diapirs, pockmarks and carbonate chimneys. During the R/V SONNE expedition “GAP-Gibraltar Arc Processes (SO-175)” in December 2003, several mud volcanoes were surveyed for gas seepage and associated microbial methane turnover. Pore water analyses and methane oxidation measurements on sediment cores recovered from the centres of the mud volcanoes Captain Arutyunov, Bonjardim, Ginsburg, Gemini and a newly discovered, mud volcano-like structure called “No Name” show that thermogenic methane and associated higher hydrocarbons rising from deeper sediment strata are completely consumed within the seabed. The presence of a distinct sulphate-methane transition zone (SMT) overlapping with high sulphide concentrations suggests that methane oxidation is mediated under anaerobic conditions with sulphate as the electron acceptor. Anaerobic oxidation of methane (AOM) and sulphate reduction (SR) rates show maxima at the SMT, which was found between 20 and 200 cm below seafloor at the different mud volcanoes. In comparison to other methane seeps, AOM activity (<383 mmol m⁻² year⁻¹) and diffusive methane fluxes (<321 mmol m⁻² year⁻¹) in mud volcano sediments of the Gulf of Cadiz are low to mid range. Corresponding lipid biomarker and 16S rDNA clone library analysis give evidence that AOM is mediated by a mixed community of anaerobic methanotrophic archaea and associated sulphate reducing bacteria (SRB) in the studied mud volcanoes. Little is known about the variability of methane fluxes in this environment. Carbonate crusts littering the seafloor of mud volcanoes in the northern part of the Gulf of Cadiz had strongly ¹³C-depleted lipid signatures indicative of higher seepage activities in the past. However, actual seafloor video observations showed only scarce traces of methane seepage and associated biological processes at the seafloor. No active fluid or free gas escape to the hydrosphere was observed visually at any of the surveyed mud volcanoes, and biogeochemical measurements indicate a complete methane consumption in the seafloor. Our observations suggest that the emission of methane to the hydrosphere from the mud volcano structures studied here may be insignificant at present.     

Chemical, biological and hydrological controls on the C content of cold seep carbonate crusts: numerical modeling and implications for convection at cold seeps

Understanding the hydrology of cold seep environments is crucial to perform accurate estimates of fluid and chemical fluxes at sedimentary wedges. Shallow convection processes may affect fluid flux estimates and could favor the destabilization of gas hydrate accumulations, increasing the sediment-ocean methane flux. Evidence for the occurrence of convection at cold seeps, however, is still limited. We use the concentration of ¹⁴C (D¹⁴C) in carbonate crusts formed at cold seeps of the eastern Mediterranean Sea as a tracer for convective recirculation of seawater-derived fluids. A numerical model is applied to investigate the controls on ¹⁴C incorporation in cold seep carbonates. Our simulations show that increased amounts of CH₄ in the expelled fluids result in elevated crust D¹⁴C, while high Ca²⁺ and HCO₃⁻ concentrations produce the opposite effect. Convection is the only transport process that can significantly increase crust D¹⁴C. Advection, bioirrigation, eddy diffusion and bioturbation instead, have little effect on, or produce a decrease of, crust D¹⁴C. In addition, the presence of old or modern carbon (MC) in host sediments prior to cementation and the ¹⁴C-decay associated to the time needed to form the crust contribute in defining the D¹⁴C of carbonate crusts. We then use the model to reproduce the ¹⁴C content of the eastern Mediterranean Sea crusts to constrain the chemical and hydrological conditions that led to their formation. Some crusts contain relatively low amounts of ¹⁴C (−945.0<D¹⁴C ‰<−930.2) which, assuming no ageing after crust formation, can be reproduced without considering convection. Other crusts from two sites (the Amsterdam and Napoli mud volcanoes), instead, have a very high ¹⁴C-content (−899.0<D¹⁴C ‰<−838.4) which can only be reproduced by the model if convection mixes deep fluids with seawater. Order-of-magnitude calculations using the Rayleigh criterion for convection suggest that the slow seepage (about 10 cm year⁻¹) of low salinity (20‰) fluids at the Amsterdam sites could trigger haline convection there. On the Napoli mud volcano, where high-density brines are expelled, density-driven convection cannot take place and other processes, possibly involving the rapid movement of free gas in the sediment, could be important.     

Consequences of Water Level Drops for Soft Sediment Deformation and Vertical Fluid Leakage

Both water level drops and erosion have previously been suggested as causes of fluid overpressures in the subsurface. Quantification of the relevance of these processes to supra-lithostatic fluid pressure formation with a wide selection of input parameters is lacking, and thus desired. The magnitudes and drop times that are required for water level drops to result in supra-lithostatic pore pressures in a variety of situations are calculated. Situations with pore fluids consisting of water, water with dissolved methane, water with a gas hydrate layer and dissolved methane in the underlying sediments, and water with dissolved methane, a gas hydrate layer, and free gas accumulation below the hydrate layer are separately addressed. The overpressure formation from reservoir gas expansion is also simulated. The simulation results demonstrate that high fluid overpressures can develop in a rock as a response to a water level drop without the presence of gas, provided that the rock has a sufficiently low compressibility. The contribution to fluid overpressuring is however dramatically increased if the pore water is saturated with methane prior to the water level drop, and is further amplified by the presence of gas hydrates and free gas accumulations beneath such hydrates. Gas expansion in reservoirs should be expected to significantly increase the fluid overpressures in shallow, sealed pressure compartments that experience erosion or water level drops, even if the water level drop duration exceeds one million years. Enough relationships between the calculated overpressure formation and the main controlling factors are provided in order to enable readers to make inferences about a variety of geological settings. Analyses of simulation results prompt us to suggest that pockmarks are likely to be triggered by gas expansion in vertical fluid migration pathways, that the giant craters at the seabed west of Albatross South in the Barents Sea result from hydrate dissociation, and that overpressure build-up due to gas expansion has contributed to reservoir overpressuring in many eroded basins, including the Hammerfest Basin in the Barents Sea.     

Weak coupling between sulfate reduction and the anaerobic oxidation of methane in methane-rich seafloor sediments during ex situ incubation

We investigated coupling between sulfate reduction (SR) and anaerobic oxidation of methane (AOM) by quantifying pore water geochemical profiles, determining rates of microbial processes, and examining microbial community structure at two sites within Mississippi Canyon lease block 118 (MC118) in the Northern Gulf of Mexico. Sediments from the northwest seep contained high concentrations of methane while sediments from the southwest seep contained methane, gaseous n-alkanes and liquid hydrocarbons and had abundant surficial accumulations of gas hydrate. Volumetric (21.5 μmol cm⁻³ day⁻¹) and integrated (1429 mmol m⁻² day⁻¹) rates of SR at MC118 in ex situ incubations are the highest reported thus far for seafloor environments. AOM rates were small in comparison, with volumetric rates ranging from 0.1 to 12.6 nmol cm⁻³ day⁻¹. Diffusion cannot adequately supply the sulfate required to support these high SR rates so additional mechanisms, possibly biological sulfide oxidation and/or downward advection, play important roles in supplying sulfate at these sites. The microbial communities at MC118 included sulfate-reducing bacteria phylogenetically associated with Desulfobacterium anilini, which is capable of complex hydrocarbon degradation. Despite low AOM rates, the majority of archaea identified were phylogenetically related to previously described methane oxidizing archaea. To evaluate whether weak coupling between SR and AOM occurs in habitats lacking the complex hydrocarbon milieu present at MC118, we compiled available SR and AOM rates and found that the global median ratio of SR to AOM was 10.7:1 rather than the expected 1:1. The global median integrated AOM rate was used to refine global estimates for AOM rates at cold seeps; these new estimates are only 5% of the previous estimate.