Abstract: Several process-based models for simulating the growth of perennial grasses have been developed but few include the simulation of regrowth. The model CATIMO simulates the primary growth of timothy (Phleum pratense L), an important perennial forage grass species in northern regions of Europe and North America. Our objective was to further develop the model CATIMO to simulate timothy regrowth using the concept of reserve-dependent growth. The performance of this modified CATIMO model in simulating leaf area index (LAI), biomass dry matter (DM) yield, and N uptake of regrowth was assessed with data from four independent field experiments in Norway, Finland, and western and eastern Canada using an approach that combines graphical comparison and statistical analysis. Biomass DM yield and N uptake of regrowth were predicted at the same accuracy as primary growth with linear regression coefficients of determination between measured and simulated values greater than 0.79, model simulation efficiencies greater than 0.78, and normalized root mean square errors (14-30% for biomass and 24-34% for N uptake) comparable with the coefficients of variation of measured data (1-21% for biomass and 1-25% for N uptake). The model satisfactorily simulated the regrowth LAI but only up to a value of about 4.0. The modified CATIMO model with its capacity to simulate regrowth provides a framework to simulate perennial grasses with multiple harvests, and to explore management options for sustainable grass production under different environmental conditions. Crown Copyright (C) 2012 Published by Elsevier B.V. All rights reserved.

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.