A quantitative assessment of methane cycling in Hikurangi Margin sediments (New Zealand) using geophysical imaging and biogeochemical modeling

Takahe seep, located on the Opouawe Bank, Hikurangi Margin, is characterized by a well-defined subsurface seismic chimney structure ca. 80,500 m2 in area. Sub-seafloor geophysical data based on acoustic anomaly layers indicated the presence of gas hydrate and free gas layers within the chimney structure. Reaction-transport modeling was applied to porewater data from 11 gravity cores to constrain methane turnover rates and benthic methane fluxes in the upper 10 m. Model results show that methane dynamics were highly variable due to transport and dissolution of ascending gas. The dissolution of gas (up to 3761 mmol m−2 yr−1) dwarfed the rate of methanogenesis within the simulated sediment column (2.6 mmol m−2 yr−1). Dissolved methane is mainly consumed by anaerobic oxidation of methane (AOM) at the base of the sulfate reduction zone and trapped by methane hydrate formation below it, with maximum rates in the central part of the chimney (946 and 2420 mmol m−2 yr−1, respectively). A seep-wide methane budget was constrained by combining the biogeochemical model results with geophysical data and led to estimates of AOM rates, gas hydrate formation and benthic dissolved methane fluxes of 3.68 × 104 mol yr−1, 73.85 × 104 mol yr−1and 1.19 × 104 mol yr−1, respectively. A much larger flux of methane probably escapes in gaseous form through focused bubble vents. The approach of linking geochemical model results with spatial geophysical data put forward here can be applied elsewhere to improve benthic methane turnover rates from limited single spot measurements to larger spatial scales