Abstract: Crop models of crop growth are increasingly used to quantify the impact of global changes due to climate or crop management. Therefore, accuracy of simulation results is a major concern. Studies with ensembles of crop models can give valuable information about model accuracy and uncertainty, but such studies are difficult to organize and have only recently begun. We report on the largest ensemble study to date, of 27 wheat models tested in four contrasting locations for their accuracy in simulating multiple crop growth and yield variables. The relative error averaged over models was 24-38% for the different end-of-season variables including grain yield (GY) and grain protein concentration (GPC). There was little relation between error of a model for GY or GPC and error for in-season variables. Thus, most models did not arrive at accurate simulations of GY and GPC by accurately simulating preceding growth dynamics. Ensemble simulations, taking either the mean (e-mean) or median (e-median) of simulated values, gave better estimates than any individual model when all variables were considered. Compared to individual models, e-median ranked first in simulating measured GY and third in GPC. The error of e-mean and e-median declined with an increasing number of ensemble members, with little decrease beyond 10 models. We conclude that multimodel ensembles can be used to create new estimators with improved accuracy and consistency in simulating growth dynamics. We argue that these results are applicable to other crop species, and hypothesize that they apply more generally to ecological system models.

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.