Abstract: The production of durum wheat in the Mediterranean basin is expected to experience increased variability in yield and quality as a consequence of climate change. To assess how environmental variables and agronomic practices affect grain protein content (GPC), a novel approach based on monthly gridded input data has been implemented to develop empirical model, and validated on historical time series to assess its capability to reproduce observed spatial and inter-annual GPC variability. The model was applied in four Italian regions and at the whole national scale and proved reliable and usable for operational purposes also in a forecast ‘real-time’ mode before harvesting. Precipitable water during autumn to winter and air temperature from anthesis to harvest were extremely important influences on GPC; these and additional variables, included in a linear model, were able to account for 95% of the variability in GPC that has occurred in the last 15 years in Italy. Our results are a unique example of the use of modelling as a predictive real-time platform and are a useful tool to understand better and forecast the impacts of future climate change projections on durum wheat production and quality.

Abstract: Due to the more frequent use of crop models at regional and national scale, the effects of spatial data input resolution have gained increased attention. However, little is known about the influence of variability in crop management on model outputs. A constant and uniform crop management is often considered over the simulated area and period. This study determines the influence of crop management adapted to climatic conditions and input data resolution on regional-scale outputs of crop models. For this purpose, winter wheat and maize were simulated over 30 years with spatially and temporally uniform management or adaptive management for North Rhine-Westphalia ((similar to)34 083 km(2)), Germany. Adaptive management to local climatic conditions was used for 1) sowing date, 2) N fertilization dates, 3) N amounts, and 4) crop cycle length. Therefore, the models were applied with four different management sets for each crop. Input data for climate, soil and management were selected at five resolutions, from 1 x 1 km to 100 x 100 km grid size. Overall, 11 crop models were used to predict regional mean crop yield, actual evapotranspiration, and drainage. Adaptive management had little effect (< 10% difference) on the 30-year mean of the three output variables for most models and did not depend on soil, climate, and management resolution. Nevertheless, the effect was substantial for certain models, up to 31% on yield, 27% on evapotranspiration, and 12% on drainage compared to the uniform management reference. In general, effects were stronger on yield than on evapotranspiration and drainage, which had little sensitivity to changes in management. Scaling effects were generally lower than management effects on yield and evapotranspiration as opposed to drainage. Despite this trend, sensitivity to management and scaling varied greatly among the models. At the annual scale, effects were stronger in certain years, particularly the management effect on yield. These results imply that depending on the model, the representation of management should be carefully chosen, particularly when simulating yields and for predictions on annual scale.

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.