Abstract: At the plot level, crop simulation models such as STICS have the potential to evaluate risk associated with management practices. In nitrogen (N) management, however, the decision-making process is complex because the decision has to be taken without any knowledge of future weather conditions. The objective of this paper is to present a general methodology for assessing yield variability linked to climatic uncertainty and variable N rate strategies. The STICS model was coupled with the LARS-Weather Generator. The Pearson system and coefficients were used to characterise the shape of yield distribution. Alternatives to classical statistical tests were proposed for assessing the normality of distributions and conducting comparisons (namely, the Jarque-Bera and Wilcoxon tests, respectively). Finally, the focus was put on the probability risk assessment, which remains a key point within the decision process. The simulation results showed that, based on current N application practice among Belgian farmers (60-60-60 kgN ha(-1)), yield distribution was very highly significantly non-normal, with the highest degree of asymmetry characterised by a skewness value of -1.02. They showed that this strategy gave the greatest probability (60 %) of achieving yields that were superior to the mean (10.5 t ha(-1)) of the distribution.

Abstract: ARMOSA is a dynamic simulation model which was developed to simulate crop growth and development, water and nitrogen dynamics under different pedoclimatic conditions and cropping systems in the arable land. The model is meant to be a tool for the evaluation of the impact of different crop management practices on soil nitrogen and carbon cycles and groundwater nitrate pollution. A large data set collected over three to six years from six monitoring sites in Lombardia plain was used to calibrate and validate the model parameters. Measured meteorological data, soil chemical and physical characterizations, crop-related data of different cropping systems allowed for a proper parameterization. Fit indexes showed the reliability of the model in adequately predicting crop-related variables, such as above ground biomass (RRMSE=11.18, EF=0.94, r=0.97), Leaf Area Index maximum value (RRMSE=8.24, EF=0.37, r=0.72), harvest index (RRMSE=19.4, EF=0.32, r=0.74), and crop N uptake (RRMSE=20.25, EF=0.69, r=0.85). Using two different one-year data set from each monitoring site, the model was calibrated and validated, getting to encouraging results: RRMSE=6.28, EF=0.52, r=0.68 for soil water content at different depths, and RRMSE=34.89, EF=0.59, r=0.75 for soil NO3-N content along soil profile. The simulated N leaching was in full agreement with measured data (RRMSE=26.62, EF=0.88, r=0.98).

Abstract: Within the context of nitrogen (N) management, since 1950, with the rapid intensification of agriculture, farmers have often applied much larger fertiliser quantities than what was required to reach the yield potential. However, to prevent pollution of surface and groundwater induced by nitrates, The European Community launched The European Nitrates Directive 91/6/76/EEC. In 2002, in Wallonia (Belgium), the Nitrates Directive has been transposed under the Sustainable Nitrogen Management in Agriculture Program (PGDA), with the aim of maintaining productivity and revenue for the country’s farmers, while reducing the environmental impact of excessive N application. A feasible approach for addressing climatic uncertainty lies in the use of crop models such as the one commonly known as STICS (simulateur multidisciplinaire pour les cultures standard). These models allow the impact on crops of the interaction between cropping systems and climatic records to be assessed. Comprehensive historical climatic records are rare, however, and therefore the yield distribution values obtained using such an approach can be discontinuous. In order to obtain better and more detailed yield distribution information, the use of a high number of stochastically generated climate time series was proposed, relying on the LARS-Weather Generator. The study focused on the interactions between varying N practices and climatic conditions. Historically and currently, Belgian farmers apply 180 kg N ha(-1), split into three equal fractions applied at the tillering, stem elongation and flag-leaf stages. This study analysed the effectiveness of this treatment in detail, comparing it to similar practices where only the N rates applied at the flag-leaf stage were modified. Three types of farmer decision-making were analysed. The first related to the choice of N strategy for maximising yield, the second to obtaining the highest net revenue, and the third to reduce the environmental impact of potential N leaching, which carries the likelihood of taxation if inappropriate N rates are applied. The results showed reduced discontinuity in the yield distribution values thus obtained. In general, the modulation of N levels to accord with current farmer practices showed considerable asymmetry. In other words, these practices maximised the probability of achieving yields that were at least superior to the mean of the distribution values, thus reducing risk for the farmers. The practice based on applying the highest amounts (60-60-100 kg N ha(-1)) produced the best yield distribution results. When simple economical criteria were computed, the 60-60-80 kg N ha(-1) protocol was found to be optimal for 80-90% of the time. There were no statistical differences, however, between this practice and Belgian farmers’ current practice. When the taxation linked to a high level of potentially leachable N remaining in the soil after harvest was considered, this methodology clearly showed that, in 3 years out of 4,30 kg N ha(-1) could systematically be saved in comparison with the usual practice.