Abstract: Process-based models (PBM) for simulation of weather dependent grass growth can assist farmers and plant breeders in addressing the challenges of climate change by simulating alternative roads of adaptation. They can also provide management decision support under current conditions. A drawback of existing grass models is that they do not take into account the effect of winter stresses, limiting their use for full-year simulations in areas where winter survival is a key factor for yield security. Here, we present a novel full-year PBM for grassland named BASGRA. It was developed by combining the LINGRA grassland model (Van Oijen et al., 2005a) with models for cold hardening and soil physical winter processes. We present the model and show how it was parameterized for timothy (Phleum pratense L.), the most important forage grass in Scandinavia and parts of North America and Asia. Uniquely, BASGRA simulates the processes taking place in the sward during the transition from summer to winter, including growth cessation and gradual cold hardening, and functions for simulating plant injury due to low temperatures, snow and ice affecting regrowth in spring. For the calibration, we used detailed data from five different locations in Norway, covering a wide range of agroclimatic regions, day lengths (latitudes from 59 degrees to 70 degrees N) and soil conditions. The total dataset included 11 variables, notably above-ground dry matter, leaf area index, tiller density, content of C reserves, and frost tolerance. All data were used in the calibration. When BASGRA was run with the maximum a-posteriori (MAP) parameter vector from the single, Bayesian calibration, nearly all measured variables were simulated to an overall normalized root mean squared error (NRMSE) <0.5. For many site x experiment combinations, NRMSE was <0.3. The temporal dynamics were captured well for most variables, as evaluated by comparing simulated time courses versus data for the individual sites. The results may suggest that BASGRA is a reasonably robust model, allowing for simulation of growth and several important underlying processes with acceptable accuracy for a range of agroclimatic conditions. However, the robustness of the model needs to be tested further using independent data from a wide range of growing conditions. Finally we show an example of application of the model, comparing overwintering risks in two climatically different sites, and discuss future model applications. Further development work should include improved simulation of the dynamics of C reserves, and validation of winter tiller dynamics against independent data. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Predicting regional and global carbon (C) and water dynamics on grasslands has become of major interest, as grasslands are one of the most widespread vegetation types worldwide, providing a number of ecosystem services (such as forage production and C storage). The present study is a contribution to a regional-scale analysis of the C and water cycles on managed grasslands. The mechanistic biogeochemical model PaSim (Pasture Simulation model) was evaluated at 12 grassland sites in Europe. A new parameterization was obtained on a common set of eco-physiological parameters, which represented an improvement of previous parameterization schemes (essentially obtained via calibration at specific sites). We found that C and water fluxes estimated with the parameter set are in good agreement with observations. The model with the new parameters estimated that European grassland are a sink of C with 213 g C m(-2) yr(-1), which is close to the observed net ecosystem exchange (NEE) flux of the studied sites (185 g C m(-2) yr(-1) on average). The estimated yearly average gross primary productivity (GPP) and ecosystem respiration (RECO) for all of the study sites are 1220 and 1006 g C m(-2) yr(-1), respectively, in agreement with observed average GPP (1230 g C m(-2) yr(-1)) and RECO (1046 g C m(-2) yr(-1)). For both variables aggregated on a weekly basis, the root mean square error (RMSE) was similar to 5-16 g C week(-1) across the study sites, while the goodness of fit (R-2) was similar to 0.4-0.9. For evapotranspiration (ET), the average value of simulated ET (415 mmyr(-1)) for all sites and years is close to the average value of the observed ET (451 mm yr(-1)) by flux towers (on a weekly basis, RMSE similar to 2-8 mm week(-1); R-2 = 0.3-0.9). However, further model development is needed to better represent soil water dynamics under dry conditions and soil temperature in winter. A quantification of the uncertainties introduced by spatially generalized parameter values in C and water exchange estimates is also necessary. In addition, some uncertainties in the input management data call for the need to improve the quality of the observational system.