Integrating Natural Gas Hydrates in the Global Carbon Cycle

We produced a two-dimensional geological time- and basin-scale model of the sedimentary margin in passive and active settings, for the simulation of the deep sedimentary methane cycle including hydrate formation. Simulation of geochemical data required development of parameterizations for bubble transport in the sediment column, and for the impact of the heterogeneity in the sediment pore fluid flow field, which represent new directions in modeling methane hydrates. The model is somewhat less sensitive to changes in ocean temperature than our previous 1-D model, due to the different methane transport mechanisms in the two codes (pore fluid flow vs. bubble migration). The model is very sensitive to reasonable changes in organic carbon deposition through geologic time, and to details of how the bubbles migrate, in particular how efficiently they are trapped as they rise through undersaturated or oxidizing chemical conditions and the hydrate stability zone. The active margin configuration reproduces the elevated hydrate saturations observed in accretionary wedges such as the Cascadia Margin, but predicts a decrease in the methane inventory per meter of coastline relative to a comparable passive margin case, and a decrease in the hydrate inventory with an increase in the plate subduction rate