Abstract
Methane is a potent greenhouse gas, which traps significantly more heat than carbon dioxide over short timescales, but remains underrepresented in models used for projecting future climate scenarios, particularly in permafrost landscapes. Understanding methane dynamics in northern wetlands is critical for improving projections of high-latitude carbon feedbacks to the global climate system. To address this, we integrated ecological processes driving methane production, oxidation, and transport pathways into the DVM-DOS-TEM terrestrial biosphere model. Using long-term site observations from a thermokarst bog in boreal Alaska, we calibrated and benchmarked the model against observations of methane and carbon dioxide fluxes, carbon and nitrogen stocks, and soil temperature and moisture. Parameter sensitivity analyses identified opposing correlations for methanogenesis and heterotrophic respiration, revealing a dependence on carbon availability. We also found sensitivity to the distribution of carbon relative to the water table position. Investigation of the dominance of different methane transport pathways demonstrated the need for observations of methane flux partitioning, and the utility in comparing simulated to observed seasonality. The cumulative methane efflux projected to the year 2100 had a range of 211-3470 g C m(-2) between low and high warming scenarios. Methane emissions dominated by diffusion showed the greatest variability across projections. Increased soil temperature and carbon availability from permafrost thaw resulted in increasing methane emissions between 2030 and 2060. But methane emission was then limited by a deepening of the water table. Nevertheless, we estimated a positive radiative forcing (i.e., warming) from these mid-century methane emissions that persisted until 2100. An ebullition-dominant parameterization led to lower variability but a net negative radiative forcing (i.e., cooling) on average. Our study highlights the importance of representing methane emission pathways and the uncertainty associated with partitioning them on predicting the carbon budget and radiative forcing of wetland ecosystems in high latitudes.