We analysed 10 years (2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017) of continuous eddy covariance (EC) CO 2 flux measurements of net ecosystem exchange (NEE) in a young pedunculate oak forest in Croatia. Measured NEE was gap-filled and partitioned into gross primary productivity (GPP) and ecosystem reparation (R ECO ) using the online tool by Max Planck Institute for Biogeochemistry in Jena, Germany. Annual NEE, GPP, and R ECO were correlated with main environmental drivers. Net primary productivity was estimated from EC (NPP EC ), as a sum of −NEE and R h obtained using a constant R h :R ECO ratio, and from independent periodic biometric measurements (NPP BM ). For comparing the NPP at the seasonal level, we propose a simple model that aimed at accounting for late-summer and autumn carbon storage in the non-structural carbohydrate pool. Over the study period, Jastrebarsko forest acted as a carbon sink, with an average (±std. dev.) annual NEE of −319 (±94) gC m −2 year −1 , GPP of 1594 (±109) gC m −2 year −1 , and R ECO of 1275 (±94) gC m −2 year −1 . Annual NEE showed high inter-annual variability and poor correlation with annual average global radiation, air temperature, and total precipitation, but significant (R 2 = 0.501, p = 0.02) correlation with the change in soil water content between May and September. Comparison of annual NPP EC and NPP BM showed a good overall agreement (R 2 = 0.463, p = 0.03), although in all years NPP BM was lower than NPP EC , with averages of 680 (±88) gC m −2 year −1 and 819 (±89) gC m −2 year −1 , respectively. Lower values of NPP BM indicate that fine roots and grasses contributions to NPP, which were not measured in the study period, could have an important contribution to the overall ecosystem NPP. At a seasonal level, two NPP estimates showed differences in their dynamic, but the application of the proposed model greatly improved the agreement in the second part of the growing season. Further research is needed on the respiration partitioning and mechanisms of carbon allocation.
IntroductionCurrently global forests store approximately 30% of total anthropogenic CO 2 emissions [1], however, the potential of forests to act as a carbon sink in the future is uncertain due to possible saturation effect [2], or negative impact of changed environmental conditions on forest productivity [3]. The eddy covariance (EC) technique is a widely used, state-of-the-art method, and it has become a standard in the estimation and monitoring of high frequency (typically half-hourly) carbon and water fluxes within terrestrial ecosystems [4]. Long-term data series of net ecosystem carbon exchange (NEE) or net ecosystem productivity (NEP, NEP = −NEE), in combination with meteorological and biometric measurements, provide invaluable information on the response of forest ecosystem to environmental conditions and climate change [5][6][7][8].EC data are essential for calibration of process-based models like Biome-BGC and its variants [9,10] and they are also used for calibration and testi...