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: Highlights • Regularization techniques for regression outperformed the classical regression techniques in predicting crop yields. • Different regression techniques with similar prediction accuracy showed different responses of major climatic variables to crop yield. • The regression models showed some responses of crop yield to climatic conditions that is mostly absent in process based crop models. Abstract Statistical regression models represent alternatives to process-based dynamic models for predicting the response of crop yields to variation in climatic conditions. Regression models can be used to quantify the effect of change in temperature and precipitation on yields. However, it is difficult to identify the most relevant input variables that should be included in regression models due to the high number of candidate variables and to their correlations. This paper compares several regression techniques for modeling response of winter oilseed rape yield to a high number of correlated input variables. Several statistical regression methods were fitted to a dataset including 689 observations of winter oilseed rape yield from replicated field experiments conducted in 239 sites in Denmark, covering nearly all regions of the country from 1992 to 2013. Regression methods were compared by cross-validation. The regression methods leading to the most accurate yield predictions were Lasso and Elastic Net, and the least accurate methods were ordinary least squares and stepwise regression. Partial least squares and ridge regression methods gave intermediate results. The estimated relative yield change for a +1°C temperature increase during flowering was estimated to range between 0 and +6 %, depending on choice of regression method. Precipitation was found to have an adverse effect on yield during autumn and winter. It was estimated that an increase in precipitation of +1 mm/day would result in a relative yield change ranging from 0 to −4 %. Soil type was also important for crop yields with lower yields on sandy soils compared to loamy soils. Later sowing was found to result in increased crop yield. The estimated effect of climate on yield was highly sensitive to the chosen regression method. Regression models showing similar performance led in some cases to different conclusions with respect to effect of temperature and precipitation. Hence, it is recommended to apply an ensemble of regression models, in order to account for the sensitivity of the data driven models for projecting crop yield under climate change.