Ages and transit times as important diagnostics of model performance for predicting carbon dynamics in terrestrial vegetation models

The global carbon cycle is strongly controlled by the source/sink strength of vegetation as well as the capacity of terrestrial ecosystems to retain this carbon. These dynamics, as well as processes such as the mixing of old and newly fixed carbon, have been studied using ecosystem models, but different assumptions regarding the carbon allocation strategies and other model structures may result in highly divergent model predictions. We modeled three systems of vegetation compartments and assessed their performance by calculating the age of the carbon in vegetation system and within each compartment, and the overall transit time of C in the system. We used these diagnostics to assess the influence of three different carbon allocation schemes on the rates of C cycling in vegetation. First, we used published measurements of ecosystem C compartments from the Harvard Forest Environmental Measurement Site to find the best set of parameters for the different model structures. Second, we calculated C stocks, release fluxes, radiocarbon values based on the bomb spike, ages, and transit times. We found a good fit of the three model structures to the available data, but the time series of C in foliage and wood need to be complemented with other ecosystem compartments in order to reduce the high parameter collinearity that we observed, and reduce model equifinality. Differences in model structures had a small impact on predicting C stocks in ecosystem compartments, but overall they resulted in very different predictions of age and transit time distributions. In particular, the inclusion of two storage compartments resulted in the prediction of a system mean age that was 10–20 years older than in the models with one or no storage compartments. The age of carbon in the wood compartment of this model was also distributed towards older ages, whereas fast cycling compartments had an age distribution that did not exceed 5 years. As expected, models with C distributed towards older ages also had longer transit times. These results suggest that ages and transit times, which can be indirectly measured using isotope tracers, serve as important diagnostics of model structure and could largely help to reduce uncertainties in model predictions. Furthermore, by considering age and transit times of C in vegetation compartments as distributions, not only their mean values, we obtain additional insights on the temporal dynamics of carbon use, storage, and allocation to plant parts, which not only depends on the rate at which this C is transferred in and out of the compartments, but also on the stochastic nature of the process itself