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: 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.