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: Field-scale crop models are often applied at spatial resolutions coarser than that of the arable field. However, little is known about the response of the models to spatially aggregated climate input data and why these responses can differ across models. Depending on the model, regional yield estimates from large-scale simulations may be biased, compared to simulations with high-resolution input data. We evaluated this so-called aggregation effect for 13 crop models for the region of North Rhine-Westphalia in Germany. The models were supplied with climate data of 1 km resolution and spatial aggregates of up to 100 km resolution raster. The models were used with 2 crops (winter wheat and silage maize) and 3 production situations (potential, water-limited and nitrogen-water-limited growth) to improve the understanding of errors in model simulations related to data aggregation and possible interactions with the model structure. The most important climate variables identified in determining the model-specific input data aggregation on simulated yields were mainly related to changes in radiation (wheat) and temperature (maize). Additionally, aggregation effects were systematic, regardless of the extent of the effect. Climate input data aggregation changed the mean simulated regional yield by up to 0.2 t ha(-1), whereas simulated yields from single years and models differed considerably, depending on the data aggregation. This implies that large-scale crop yield simulations are robust against climate data aggregation. However, large-scale simulations can be systematically biased when being evaluated at higher temporal or spatial resolution depending on the model and its parameterization.

Abstract: Agro-meteorological quantities are often in the form of time series, and knowledge about their temporal scaling properties is fundamental for transferring locally measured fluctuations to larger scales and vice versa. However, the scaling analysis of these quantities is complicated due to the presence of localized trends and nonstationarities. The objective of this study was to characterise scaling properties (i.e. statistical self-similarity) of the chosen agro-meteorological quantities through multifractal detrended fluctuation analysis (MFDFA). For this purpose, MFDFA was performedwith 11 322 measured time series (31 yr) of daily air temperature, wind velocity, relative air humidity, global radiation and precipitation from stations located in Finland, Germany, Poland and Spain. The empirical singularity spectra indicated their multifractal structure. The richness of the studied multifractals was evaluated by the width of their spectrum, indicating considerable differences in dynamics and development. In log-log plots of the cumulative distributions of all meteorological parameters the linear functions prevailed for high values of the response, indicating that these distributions were consistent with power-law asymptotic behaviour. Additionally, we investigated the type of multifractality that underlies the q-dependence of the generalized Hurst exponent by analysing the corresponding shuffled and surrogate time series. For most of the studied meteorological parameters, the multifractality is due to different long-range correlations for small and large fluctuations. Only for precipitation does the multifractality result mainly from broad probability function. This feature may be especially valuable for assessing the effect of change in climate dynamics.