Abstract: Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM((ZS)) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m(-2) day(-1) and 0.50 to 1.28 mm day(-1), respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic mortality time parameters can be explained by site-specific ecophysiological characteristics.

Abstract: CARAIB (Dury et al., 2011) is a state-of-the-art dynamic vegetation model with various modules dealing with (i) soil hydrology, (ii) photosynthesis/stomatal regulation, (iii) carbon allocation and biomass growth, (iv) litter/soil carbon dynamics, (v) vegetation cover dynamics, (vi) seed dispersal, and (vii) vegetation fires. Climate and atmospheric CO2 are the primary inputs. The model calculates all major water and CO2/carbon fluxes and pools. It can be run with plant functional types or species (up to 100 different species) at various spatial scales, from the municipality to country or continental levels.   Within the VOTES project (Fontaine et al., 2013), the model has been improved to include crops and meadows, and some modules have been written to translate model outputs into quantitative indicators of ecosystem services (e.g., evaluate crop yield from net primary productivity or calculate soil erosion from runoff, slope, grown species and various soil attributes). The model was run over an area covering four municipalities in central Belgium, where land-use is dominated by crops, meadows, housing and some forests and was introduced in the model at the land parcel level. Simulations were also performed for the future. In these simulations, CARAIB was combined with the Aporia Agent-Based Model,  to project land-use changes up to 2050. This approach is currently extended within the MASC project (funded by Belgian Science Policy, BELSPO) to the whole Belgian territory (at 1 km2) and to Western Europe (at 20  km x 20 km).