Drivers of dissolved organic carbon export in a subarctic catchment: Importance of microbial decomposition, sorption-desorption, peatland and lateral flow

Tundra soils account for 50% of global stocks of soil organic carbon (SOC), and it is expected that the amplified climate warming in high latitude could cause loss of this SOC through decomposition. Decomposed SOC could become hydrologically accessible, which increase downstream dissolved organic carbon (DOC) export and subsequent carbon release to the atmosphere, constituting a positive feedback to climate warming. However, DOC export is often neglected in ecosystem models. In this paper, we incorporate processes related to DOC production, mineralization, diffusion, sorption-desorption, and leaching into a customized arctic version of the dynamic ecosystem model LPJ-GUESS in order to mechanistically model catchment DOC export, and to link this flux to other ecosystem processes. The extended LPJ-GUESS is compared to observed DOC export at Stordalen catchment in northern Sweden. Vegetation communities include flood-tolerant graminoids (Eriophorum) and Sphagnum moss, birch forest and dwarf shrub communities. The processes, sorption-desorption and microbial decomposition (DOC production and mineralization) are found to contribute most to the variance in DOC export based on a detailed variance-based Sobol sensitivity analysis (SA) at grid cell-level. Catchment-level SA shows that the highest mean DOC exports come from the Eriophorum peatland (fen). A comparison with observations shows that the model captures the seasonality of DOC fluxes. Two catchment simulations, one without water lateral routing and one without peatland processes, were compared with the catchment simulations with all processes. The comparison showed that the current implementation of catchment lateral flow and peatland processes in LPJ-GUESS are essential to capture catchment-level DOC dynamics and indicate the model is at an appropriate level of complexity to represent the main mechanism of DOC dynamics in soils. The extended model provides a new tool to investigate potential interactions among climate change, vegetation dynamics, soil hydrology and DOC dynamics at both stand-alone to catchment scales.