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.