Abstract. Northern peatlands in permafrost regions contain a large amount of organic carbon (C) in the soil. Climate warming and associated permafrost degradation are expected to have significant impacts on the C balance of these ecosystems, but the magnitude is uncertain. We incorporated a permafrost model, Northern Ecosystem Soil Temperature (NEST), into a biogeochemical model, DeNitrificationDeComposition (DNDC), to model C dynamics in highlatitude peatland ecosystems. The enhanced model was applied to assess effects of permafrost thaw on C fluxes of a subarctic peatland at Stordalen, Sweden. DNDC simulated soil freeze-thaw dynamics, net ecosystem exchange of CO 2 (NEE), and CH 4 fluxes across three typical land cover types, which represent a gradient in the process of ongoing permafrost thaw at Stordalen. Model results were compared with multiyear field measurements, and the validation indicates that DNDC was able to simulate observed differences in seasonal soil thaw, NEE, and CH 4 fluxes across the three land cover types. Consistent with the results from field studies, the modeled C fluxes across the permafrost thaw gradient demonstrate that permafrost thaw and the associated changes in soil hydrology and vegetation not only increase net uptake of C from the atmosphere but also increase the annual to decadal radiative forcing impacts on climate due to increased CH 4 emissions. This study indicates the potential of utilizing biogeochemical models, such as DNDC, to predict the soil thermal regime in permafrost areas and to investigate impacts of permafrost thaw on ecosystem C fluxes after incorporating a permafrost component into the model framework.