Does growing atmospheric CO2 explain increasing carbon sink in a boreal coniferous forest?

Abstract

The terrestrial net ecosystem productivity (NEP) has increased during the past three decades, but the mechanisms responsible are still unclear. We analyzed 17 years (2001–2017) of eddy-covariance measurements of NEP, evapotranspiration (ET) and light and water use efficiency from a boreal coniferous forest in Southern Finland for trends and inter-annual variability (IAV). The forest was a mean annual carbon sink (252 [(Formula presented.) 42] gC (Formula presented.) (Formula presented.)), and NEP increased at rate +6.4–7.0 gC (Formula presented.) (Formula presented.) (or ca. +2.5% (Formula presented.)) during the period. This was attributed to the increasing gross-primary productivity GPP and occurred without detectable change in ET. The start of annual carbon uptake period was advanced by 0.7 d (Formula presented.), and increase in GPP and NEP outside the main growing season contributed ca. one-third and one-fourth of the annual trend, respectively. Meteorological factors were responsible for the IAV of fluxes but did not explain the long-term trends. The growing season GPP trend was strongest in ample light during the peak growing season. Using a multi-layer ecosystem model, we showed that direct (Formula presented.) fertilization effect diminishes when moving from leaf to ecosystem, and only 30–40% of the observed ecosystem GPP increase could be attributed to (Formula presented.). The increasing trend in leaf-area index (LAI), stimulated by forest thinning in 2002, was the main driver of the enhanced GPP and NEP of the mid-rotation managed forest. It also compensated for the decrease of mean leaf stomatal conductance with increasing (Formula presented.) and LAI, explaining the apparent proportionality between observed GPP and (Formula presented.) trends. The results emphasize that attributing trends to their physical and physiological drivers is challenged by strong IAV, and uncertainty of LAI and species composition changes due to the dynamic flux footprint. The results enlighten the underlying mechanisms responsible for the increasing terrestrial carbon uptake in the boreal zone.