Abstract: Two independently developed simulation models – the grassland-specific PaSim and the biome-generic Biome-BGC MuSo (BBGC MuSo) – linking climate, soil, vegetation and management to ecosystem biogeochemical cycles were compared in a simulation of carbon (C) and water fluxes. The results were assessed against eddy-covariance flux data from five observational grassland sites representing a range of conditions in Europe: Grillenburg in Germany, Laqueuille in France with both extensive and intensive management, Monte Bondone in Italy and Oensingen in Switzerland. Model comparison (after calibration) gave substantial agreement, the performances being marginal to acceptable for weekly-aggregated gross primary production and ecosystem respiration (R-2 similar to 0.66 – 0.91), weekly evapotranspiration (R-2 similar to 0.78 – 0.94), soil water content in the topsoil (R-2 similar to 0.1 -0.7) and soil temperature (R-2 similar to 0.88 – 0.96). The bias was limited to the range -13 to 9 g C m(-2) week(-1) for C fluxes (-11 to 8 g C m(-2) week(-1) in case of BBGC MuSo, and -13 to 9 g C m(-2) week(-1) in case of PaSim) and -4 to 6 mm week for water fluxes (with BBGC MuSo providing somewhat higher estimates than PaSim), but some higher relative root mean square errors indicate low accuracy for prediction, especially for net ecosystem exchange The sensitivity of simulated outputs to changes in atmospheric carbon dioxide concentration ([CO2]), temperature and precipitation indicate, with certain agreement between the two models, that C outcomes are dominated by [CO2] and temperature gradients, and are less due to precipitation. ET rates decrease with increasing [CO2] in PaSim (consistent with experimental knowledge), while lack of appropriate stomatal response could be a limit in BBGC MuSo responsiveness. Results of the study indicate that some of the errors might be related to the improper representation of soil water content and soil temperature. Improvement is needed in the model representations of soil processes (especially soil water balance) that strongly influence the biogeochemical cycles of managed and unmanaged grasslands. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: Highlights • The parameters to be fixed were consistent across sites. • Model calibration must be performed separately for each specific case. • Possible to reduce model parameters from 66 to 45. • Strong model reductions must be avoided. • The error term for the training data were characterised by timing (phase shift). Abstract Proper parameterisation and quantification of model uncertainty are two essential tasks in improvement and assessment of model performance. Bayesian calibration is a method that combines both tasks by quantifying probability distributions for model parameters and outputs. However, the method is rarely applied to complex models because of its high computational demand when used with high-dimensional parameter spaces. We therefore combined Bayesian calibration with sensitivity analysis, using the screening method by Morris (1991), in order to reduce model complexity by fixing parameters to which model output was only weakly sensitive to a nominal value. Further, the robustness of the model with respect to reduction in the number of free parameters were examined according to model discrepancy and output uncertainty. The process-based grassland model BASGRA was examined in the present study on two sites in Norway and in Germany, for two grass species (Phleum pratense and Arrhenatherum elatius). According to this study, a reduction of free model parameters from 66 to 45 was possible. The sensitivity analysis showed that the parameters to be fixed were consistent across sites (which differed in climate and soil conditions), while model calibration had to be performed separately for each combination of site and species. The output uncertainty decreased slightly, but still covered the field observations of aboveground biomass. Considering the training data, the mean square error for both the 66 and the 45 parameter model was dominated by errors in timing (phase shift), whereas no general pattern was found in errors when using the validation data. Stronger model reduction should be avoided, as the error term increased and output uncertainty was underestimated.

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.