Abstract: • ARMOSA model simulated a maize ideotype with drought adaptation under climate change. • The ideotype needs less water for higher yield compared to current hybrids. • Higher production involves more crop residues that enhance soil C sequestration. • Soil organic C may generally decrease and N leaching will increase in sandy soil. The expected climate change will affect the maize yields in view of air temperature increase and scarce water availability. The application of biophysical models offers the chance to design a drought-resistant ideotype and to assist plant breeders and agronomists in the assessment of its suitability in future scenarios. The aim of the present work was to perform a model-based estimation of the yields of two hybrids, current vs ideotype, under future climate scenarios (2030–2060 and 2070–2100) in Lombardy (northern Italy), testing two options of irrigation (small amount at fixed dates vs optimal water supply), nitrogen (N) fertilization (300 vs 400 kg N ha− 1), and crop cycle durations (current vs extended). For the designing of the ideotype we set several parameters of the ARMOSA process-based crop model: the root elongation rate and maximum depth, stomatal resistance, four stage-specific crop coefficients for the actual transpiration estimation, and drought tolerance factor. The work findings indicated that the current hybrid ensures good production only with high irrigation amount (245–565 mm y− 1). With respect to the current hybrid, the ideotype will require less irrigation water (− 13%, p < 0.01) and it resulted in significantly higher yield under water stress condition (+ 15%, p < 0.01) and optimal water supply (+ 2%, p < 0.05). The elongated cycle has a positive effect on yield under any combination of options. Moreover, higher yields projected for the ideotype implicate more crop residues to be incorporated into the soil, which are positively correlated with the SOC sequestration and negatively with N leaching. The crop N uptake is expected to be adequate in view of higher rate of soil mineralization; the N fertilization rate of 400 kg N ha− 1 will involve significant increasing of grain yield, and it is expected to involve a higher rate of SOC sequestration.

Abstract: ► We compared nine crop simulation models for spring barley at seven sites in Europe. ► Applying crop models with restricted calibration leads to high uncertainties. ► Multi-crop model mean yield estimates were in good agreement with observations. ► The degree of uncertainty for simulated grain yield of barley was similar to winter wheat. ► We need more suitable data enabling us to verify different processes in the models. In this study, the performance of nine widely used and accessible crop growth simulation models (APES-ACE, CROPSYST, DAISY, DSSAT-CERES, FASSET, HERMES, MONICA, STICS and WOFOST) was compared during 44 growing seasons of spring barley (Hordeum vulgare L) at seven sites in Northern and Central Europe. The aims of this model comparison were to examine how different process-based crop models perform at multiple sites across Europe when applied with minimal information for model calibration of spring barley at field scale, whether individual models perform better than the multi-model mean, and what the uncertainty ranges are in simulated grain yields. The reasons for differences among the models and how results for barley compare to winter wheat are discussed. Regarding yield estimation, best performing based on the root mean square error (RMSE) were models HERMES, MONICA and WOFOST with lowest values of 1124, 1282 and 1325 (kg ha(-1)), respectively. Applying the index of agreement (IA), models WOFOST, DAISY and HERMES scored best having highest values (0.632, 0.631 and 0.585, respectively). Most models systematically underestimated yields, whereby CROPSYST showed the highest deviation as indicated by the mean bias error (MBE) (-1159 kg ha(-1)). While the wide range of simulated yields across all sites and years shows the high uncertainties in model estimates with only restricted calibration, mean predictions from the nine models agreed well with observations. Results of this paper also show that models that were more accurate in predicting phenology were not necessarily the ones better estimating grain yields. Total above-ground biomass estimates often did not follow the patterns of grain yield estimates and, thus, harvest indices were also different. Estimates of soil moisture dynamics varied greatly. In comparison, even though the growing cycle for winter wheat is several months longer than for spring barley, using RMSE and IA as indicators, models performed slightly, but not significantly, better in predicting wheat yields. Errors in reproducing crop phenology were similar, which in conjunction with the shorter growth cycle of barley has higher effects on accuracy in yield prediction.

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.