Numerical modeling study of mineralization induced by methane cold seep at the sea bottom

In order to investigate the response of authigenic minerals to gas hydrate geo-systems, the biogeochemical processes and its induced mineralization were predicted by employing the comprehensive reactive transport modeling approach. Based on the available data extracted from the northern continental slope area of the South China Sea, a 1-D vertical column model was developed. Three cases with different upward methane flux rates and three cases with different mineral compositions, i.e., a total of six cases were designed to investigate the effects of variations in the depth of sulfate methane transition zone (SMTZ) and in the mineral composition on the formation of authigenic minerals. The simulation results indicate that the SMTZ depth influenced by both the upward methane flux rate and the initial composition played an important role in the formation of authigenic minerals. The AOM reaction is intensive at the interface, and the precipitation amount of calcite is large, which is mainly controlled by AOM. When the methane leakage rate is 20 times higher than the base case, aragonite starts to precipitate. During the simulation, oligoclase, k-feldspar, smectite-Na, smectite-Ca, chlorite dissolved. Our study specific to this area as a starting point may provide a quantitative approach for investigating carbonate and pyrite formation in hydrate-bearing sediments accounting for methane oxidation and sulfate reduction. The method presented here and the model built in this study can be used for other sites with similar conditions. In addition, this study may serve as an indication for the potential natural gas hydrate reservoir in depth, and is also significant for marine carbon and sulfur cycle. •TOUGHREACT was used to simulate the AOM biogeochemical processes and their induced mineralization.•Controlling factors such as upward methane flux rates on the SMI and AOM rate are discussed.•The relationship of biogeochemical change and mineralization was analyzed.•Precipitation of carbonates and sulfide minerals was numerically obtained and compared with field observations.