Microbial community analyses inform geochemical reaction network models for predicting pathways of greenhouse gas production

The mechanisms, pathways, and rates of CO2 and CH4 production are central to understanding carbon cycling and greenhouse gas flux in wetlands. Thawing permafrost regions are of particular interest because they are disproportionally affected by climate warming and store large reservoirs of organic C that may be readily converted to CO2 and CH4 upon thaw. This conversion is accomplished by a community of microorganisms interacting in complex ways to transform large organic compounds into fatty acids and ultimately CO2 and CH4. While the central role of microbes in this process is well-known, geochemical rate models rarely integrate microbiological information. Herein, we expanded the geochemical rate model of Neumann et al., (2016, Biogeochemistry 127: 57–87) to incorporate a Bayesian probability analysis and applied the result to quantifying rates of CO2, CH4, and acetate production in closed-system incubations of peat collected from three habitats along a permafrost thaw gradient. The goals of this analysis were twofold. First, we integrated microbial community analyses with geochemical rate modeling by using microbial data to inform the best model choice among equally mathematically feasible model variants. Second, based on model results, we described changes in organic carbon transformation among habitats to understand the changing pathways of greenhouse gas production along the permafrost thaw gradient. We found that acetoclasty, hydrogenotrophy, CO2 production, and homoacetogenesis were the important reactions in this system, with little evidence for anaerobic CH4 oxidation. There was a distinct transition in the reactions across the thaw gradient. The collapsed palsa stage presents an initial disequilibrium where the abrupt (physically and temporally) change in elevation introduces freshly fixed carbon into anoxic conditions then fermentation products build up over time as the system transitions through the acid phase and electron acceptors are depleted. In the bog, fermentation slows, while methanogenesis increases. In the fully thawed fen, most of the terminal electron acceptors are depleted and the system becomes increasingly methanogenic. This suggests that as permafrost regions thaw and dry palsas transition into wet fens, CH4 emissions will rise, increasing the warming potential of these systems and accelerating climate warming feedbacks.