Abstract: Field-scale crop models are increasingly applied at spatio-temporal scales that range from regions to the globe and from decades up to 100 years. Sufficiently detailed data to capture the prevailing spatio-temporal heterogeneity in weather, soil, and management conditions as needed by crop models are rarely available. Effective sampling may overcome the problem of missing data but has rarely been investigated. In this study the effect of sampling weather data has been evaluated for simulating yields of winter wheat in a region in Germany over a 30-year period (1982-2011) using 12 process-based crop models. A stratified sampling was applied to compare the effect of different sizes of spatially sampled weather data (10, 30, 50,100, 500, 1000 and full coverage of 34,078 sampling points) on simulated wheat yields. Stratified sampling was further compared with random sampling. Possible interactions between sample size and crop model were evaluated. The results showed differences in simulated yields among crop models but all models reproduced well the pattern of the stratification. Importantly, the regional mean of simulated yields based on full coverage could already be reproduced by a small sample of 10 points. This was also true for reproducing the temporal variability in simulated yields but more sampling points (about 100) were required to accurately reproduce spatial yield variability. The number of sampling points can be smaller when a stratified sampling is applied as compared to a random sampling. However, differences between crop models were observed including some interaction between the effect of sampling on simulated yields and the model used. We concluded that stratified sampling can considerably reduce the number of required simulations. But, differences between crop models must be considered as the choice for a specific model can have larger effects on simulated yields than the sampling strategy. Assessing the impact of sampling soil and crop management data for regional simulations of crop yields is still needed.

Abstract: Salinization is one of the most serious problems confronting sustainable agriculture in semi-arid and arid regions. Accurate mapping of soil salinization and the associated risk represent a fundamental step in planning agricultural and remediation activities. Geostatistical analysis is very useful for soil quality assessment because it makes it possible to determine the spatial relationships between selected variables and to produce synthetic maps of spatial variation. The main objective of this paper was to map the soil salinization risk in the Delia-Nivolelli alluvial basin (south-western Sicily, southern Italy), using multivariate geostatistical techniques and a set of topographical, physical and soil hydraulic properties. Elevation data were collected from existing topographic maps and analysed preliminarily to improve the estimate precision of sparsely sampled primary variables. For interpolation multi-collocated cokriging was applied to the dataset, including textural and hydraulic properties and electrical conductivity measurements carried out on 128 collected soil samples, using elevation data as auxiliary variable. Spatial dependence among elevation and physical soil properties was explored with factorial kriging analysis (FKA) that could isolate and display the sources of variation acting at different spatial scales. FKA isolated significant regionalised factors which give a concise description of the complex soil physical variability at the different selected spatial scales. These factors mapped, allowed the delineation of zones at different salinisation risk to be managed separately to control and prevent salinization risk. The proposed methodology could be a valid support for land use and soil remediation planning at regional scale.