Abstract: A spatial resolution needs to be determined prior to using models to simulate crop yields at a regional scale, but a dilemma exists in compromising between different demands. A fine spatial resolution demands extensive computation load for input data assembly, model runs, and output analysis. A coarse spatial resolution could result in loss of spatial detail in variability. This paper studied the impact of spatial resolution, data aggregation and spatial heterogeneity of weather data on simulations of crop yields, thus providing guidelines for choosing a proper spatial resolution for simulations of crop yields at regional scale. Using a process-based crop model SIMPLACE (LINTUL2) and daily weather data at 1 km resolution we simulated a continuous rainfed winter wheat cropping system at the national scale of Germany. Then we aggregated the weather data to four resolutions from 10 to 100 km, repeated the simulation, compared them with the 1 km results, and correlated the difference with the intra-pixel heterogeneity quantified by an ensemble of four semivariogram models. Aggregation of weather data had small effects over regions with a flat terrain located in northern Germany, but large effects over southern regions with a complex topography. The spatial distribution of yield bias at different spatial resolutions was consistent with the intra-pixel spatial heterogeneity of the terrain and a log-log linear relationship between them was established. By using this relationship we demonstrated the way to optimize the model resolution to minimize both the number of simulation runs and the expected loss of spatial detail in variability due to aggregation effects. We concluded that a high spatial resolution is desired for regions with high spatial environmental heterogeneity, and vice versa. This calls for the development of multi-scale approaches in regional and global crop modeling. The obtained results require substantiation for other production situations, crops, output variables and for different crop models. (C) 2014 Elsevier B.V. All rights reserved.

Abstract: Modelling N transformations within cropping systems is crucial for N management optimization in order to increase N use efficiency and reduce N losses. Such modelling remains challenging because of the complexity of N cycling in soil–plant systems. In the current study, the uncertainties of six widely used process-based models (PBMs), including APSIM, CROPSYST, DAISY, FASSET, HERMES and STICS, were tested in simulating different N managements (catch crops (CC) and different N fertilizer rates) in 12-year rotations in Western Europe. Winter wheat, sugar beet and pea were the main crops, and radish was the main CC in the tested systems. Our results showed that PBMs simulated yield, aboveground biomass, N export and N uptake well with low RMSE values, except for sugar beet, which was generally less well parameterized. Moreover, PBMs provided more accurate crop simulations (i.e. N export and N uptake) compared to simulations of soil (N mineralization and soil mineral N (SMN)) and environmental variables (N leaching). The use of multi-model ensemble mean or median of four PBMs significantly reduced the mean absolute percentage error (MAPE) between simulations and observations to less than 15% for yield, aboveground biomass, N export and N uptake. Multi-model ensemble also significantly reduced the MAPE for net N mineralization and annual N leaching to around 15%, while it was larger than 20% for SMN. Generally, PBMs well simulated the CC effects on N fluxes, i.e. increasing N mineralization and reducing N leaching in both short-term and long-term, and all PBMs correctly predicted the effects of the reduced N rate on all measured variables in the study. The uncertainties of multi-model ensemble for N mineralization, SMN and N leaching were larger, mainly because these variables are influenced by plant-soil interactions and subject to cumulative long-term effects in crop rotations, which makes them more difficult to simulate. Large differences existed between individual PBMs due to the differences in formalisms for describing N processes in soil–plant systems, the skills of modelers and the model calibration level. In addition, the model performance also depended on the simulated variables, for instance, HERMES and FASSET performed better for yield and crop biomass, APSIM, DAISY and STICS performed better for N export and N uptake, STICS provided best simulation for SMN and N leaching among the six individual PBMs in the study, but all PBMs met difficulties to well predict either average or variance of soil N mineralization. Our results showed that better calibration for soil N variables is needed to improve model predictions of N cycling in order to optimize N management in crop rotations.

Abstract: The success of using agro-ecosystem models for the high-resolution simulation of agricultural yields for larger areas is often hampered by a lack of input data. We investigated the effect of different spatially resolved soil and weather data used as input for the MONICA model on its ability to reproduce winter wheat yields in the Federal State of Thuringia, Germany (16,172 km(2)). The combination of one representative soil and one weather station was insufficient to reproduce the observed mean yield of 6.66 +/- 0.87 t ha(-1) for the federal state. Use of a 100 m x 100 m grid of soil and relief information combined with just one representative weather station yielded a good estimator (7.01 +/- 1.47 t ha(-1)). The soil and relief data grid used in combination with weather information from 14 weather stations in a nearest neighbour approach produced even better results (6.60 +/- 1.37 t ha(-1)); the same grid used with 39 additional rain gauges and an interpolation algorithm that included an altitude correction of temperature data slightly overpredicted the observed mean (7.36 +/- 1.17 t ha(-1)). It was concluded that the apparent success of the first two high-resolution approaches over the latter was based on two effects that cancelled each other out: the calibration of MONICA to match high-yield experimental data and the growth-defining and -limiting effect of weather data that is not representative for large parts of the region. At the county and farm level the MONICA model failed to reproduce the 1992-2010 time series of yields, which is partly explained by the fact that many growth-reducing factors were not considered in the model. (C) 2013 Elsevier B.V. All rights reserved.