Modeling sulfate reduction in methane hydrate-bearing continental margin sediments: Does a sulfate-methane transition require anaerobic oxidation of methane?

The sulfate‐methane transition (SMT), a biogeochemical zone where sulfate and methane are metabolized, is commonly observed at shallow depths (1–30 mbsf) in methane‐bearing marine sediments. Two processes consume sulfate at and above the SMT, anaerobic oxidation of methane (AOM) and organoclastic sulfate reduction (OSR). Differentiating the relative contribution of each process is critical to estimate methane flux into the SMT, which, in turn, is necessary to predict deeper occurrences of gas hydrates in continental margin sediments. To evaluate the relative importance of these two sulfate reduction pathways, we developed a diagenetic model to compute the pore water concentrations of sulfate, methane, and dissolved inorganic carbon (DIC). By separately tracking DIC containing ¹²C and ¹³C, the model also computes δ¹³C‐DIC values. The model reproduces common observations from methane‐rich sediments: a well‐defined SMT with no methane above and no sulfate below and a δ¹³C‐DIC minimum at the SMT. The model also highlights the role of upward diffusing ¹³C‐enriched DIC in contributing to the carbon isotope mass balance of DIC. A combination of OSR and AOM, each consuming similar amounts of sulfate, matches observations from Site U1325 (Integrated Ocean Drilling Program Expedition 311, northern Cascadia margin). Without AOM, methane diffuses above the SMT, which contradicts existing field data. The modeling results are generalized with a dimensional analysis to the range of SMT depths and sedimentation rates typical of continental margins. The modeling shows that AOM must be active to establish an SMT wherein methane is quantitatively consumed and the δ¹³C‐DIC minimum occurs. The presence of an SMT generally requires active AOM.