Abstract: Rice models have been widely used in simulating and predicting rice phenology in contrasting climate zones, however the uncertainties from model structure (different equations or models) and/or model parameters were rarely investigated. Here, five rice phenological models/modules (Le., CERES-Rice, ORYZA2000, RCM, Beta Model and SIMRIW) were applied to simulate rice phenology at 23 experimental stations from 1992 to 2009 in two major rice cultivation regions of China: the northeastern China and the southwestern China. To investigate the uncertainties from model biophysical parameters, each model was run with randomly perturbed 50 sets of parameters. The results showed that the median of ensemble simulations were better than the simulation by most models. Models couldn’t simulate well in some specific years despite of parameters optimization, suggesting model structure limit model performance in some cases. The models adopting accumulative thermal time function (e.g., CERES-Rice and ORYZA2000) had better performance in the southwestern China, in contrast, those adopting exponential function (e.g., Beta model and RCM model) had better performance in the northeastern China. In northeastern China, the contribution of model structure and model parameters to model total variance was, respectively, about 55.90% and 44.10% in simulating heading date, and about 75.43% and 24.57% in simulating maturity date. In the southwestern China, the contribution of model structure and model parameters to model total variance was, respectively, about 79.97% and 27.03% in simulating heading date, about 92.15% and 7.85% in simulating maturity date. Uncertainty from model structure was the most relevant source. The results highlight that the temperature response functions of rice development rate under extreme climate conditions should be improved based on environment-controlled experimental data.

Abstract: The accurate estimation of crop grain nitrogen (N; N in grain yield) is crucial for optimizing agricultural N management, especially in crop rotations. In the present study, 12 process-based models were applied to simulate the grain N of i) seven crops in rotations, ii) across various pedo-climatic and agro-management conditions in Europe, under both continuous simulation and single year simulation, and for iv) two calibration levels, namely minimal and detailed calibration. Generally, the results showed that the accuracy of the simulations in predicting grain N increased under detailed calibration. The models performed better in predicting the grain N of winter wheat (Triticum aestivum L.), winter barley (Hordewn vulgare L.) and spring barley (Hordeum vulgare L.) compared to spring oat (Avena saliva L.), winter rye (Secale cereale L.), pea (Piswn sativum L.) and winter oilseed rape (Brassica napus L.). These differences are linked to the intensity of parameterization with better parameterized crops showing lower prediction errors. The model performance was influenced by N fertilization and irrigation treatments, and a majority of the predictions were more accurate under low N and rainfed treatments. Moreover, the multi-model mean provided better predictions of grain N compared to any individual model. In regard to the Individual models, DAISY, FASSET, HERMES, MONICA and STICS are suitable for predicting grain N of the main crops in typical European crop rotations, which all performed well in both continuous simulation and single year simulation. Our results show that both the model initialization and the cover crop effects in crop rotations should be considered in order to achieve good performance of continuous simulation. Furthermore, the choice of either continuous simulation or single year simulation should be guided by the simulation objectives (e.g. grain yield, grain N content or N dynamics), the crop sequence (inclusion of legumes) and treatments (rate and type of N fertilizer) included in crop rotations and the model formalism.

Abstract: Variability of crop yields is detrimental for food security. Under climate change its amplitude is likely to increase, thus it is essential to understand the underlying causes and mechanisms. Crop models are the primary tool to project future changes in crop yields under climate change. Asystematic overview of drivers and mechanisms of crop yield variability (YV) can thus inform crop model development and facilitate improved understanding of climate change impacts on crop yields. Yet there is a vast body of literature on crop physiology and YV, which makes a prioritization of mechanisms for implementation in models challenging. Therefore this paper takes on a novel approach to systematically mine and organize existing knowledge from the literature. The aim is to identify important mechanisms lacking in models, which can help to set priorities in model improvement. We structure knowledge from the literature in a semi-quantitative network. This network consists of complex interactions between growing conditions, plant physiology and crop yield. We utilize the resulting network structure to assign relative importance to causes of YV and related plant physiological processes. As expected, our findings confirm existing knowledge, in particular on the dominant role of temperature and precipitation, but also highlight other important drivers of YV. More importantly, our method allows for identifying the relevant physiological processes that transmit variability in growing conditions to variability in yield. We can identify explicit targets for the improvement of crop models. The network can additionally guide model development by outlining complex interactions between processes and by easily retrieving quantitative information for each of the 350 interactions. We show the validity of our network method as a structured, consistent and scalable dictionary of literature. The method can easily be applied to many other research fields.