Spatially explicit simulation of hydrologically controlled carbon and nitrogen cycles and associated feedback mechanisms in a boreal ecosystem

Ecosystem models that simulate biogeochemical processes usually ignore hydrological controls that govern them. It is quite possible that topographically driven water fluxes significantly influence the spatial distribution of C sources and sinks because of their large contribution to the local water balance. To investigate this, we simulated biogeochemical processes along with the associated feedback mechanisms in a boreal ecosystem using a spatially explicit hydroecological model, boreal ecosystem productivity simulator (BEPS)–TerrainLab V2.0, that has a tight coupling of ecophysiological, hydrological, and biogeochemical processes. First, the simulated dynamics of snowpack, soil temperature, net ecosystem productivity (NEP), and total ecosystem respiration (TER) were validated with high‐frequency measurements for 2 years. The model was able to explain 80% of the variability in NEP and 84% of the variability in TER. Further, we investigated the influence of topographically driven subsurface base flow on soil C and N cycling and on the spatiotemporal patterns of C sources and sinks using three hydrological modeling scenarios that differed in hydrological conceptualizations. In general, the scenarios that had nonexplicit hydrological representation overestimated NEP, as opposed to the scenario that had an explicit (realistic) representation. The key processes controlling the NEP differences were attributed to the combined effects of variations in photosynthesis (due to changes in stomatal conductance and nitrogen (N) availability), heterotrophic respiration, and autotrophic respiration, all of which occur simultaneously affecting NEP. Feedback relationships were also found to exacerbate the differences. We identified six types of NEP differences (biases), of which the most commonly found was due to an underestimation of the existing C sources, highlighting the vulnerability of regional‐scale ecosystem models that ignore hydrological processes.