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.

Abstract: During the past few years, several studies have compared the performance of crop simulation models to assess the uncertainties in model-based climate change impact assessments and other modelling studies. Many of these studies have concentrated on cereal crops, while fewer model comparisons have been conducted for grasses. We compared the predictions for timothy grass (Phleum pratertse L.) yields for first and second cuts along with the dynamics of above-ground biomass for the grass simulation models BASGRA and CATIMO, and the soil -crop model STICS. The models were calibrated and evaluated using field data from seven sites across Northern Europe and Canada with different climates, soil conditions and management practices. Altogether the models were compared using data on timothy grass from 33 combinations of sites, cultivars and management regimes. Model performances with two calibration approaches, cultivar-specific and generic calibrations, were compared. All the models studied estimated the dynamics of above-ground biomass and the leaf area index satisfactorily, but tended to underestimate the first cut yield. Cultivar-specific calibration resulted in more accurate first cut yield predictions than the generic calibration achieving root mean square errors approximately one third lower for the cultivar-specific calibration. For the second cut, the difference between the calibration methods was small. The results indicate that detailed soil process descriptions improved the overall model performance and the model responses to management, such as nitrogen applications. The results also suggest that taking the genetic variability into account between cultivars of timothy grass also improves the yield estimates. Calibrations using both spring and summer growth data simultaneously revealed that processes determining the growth in these two periods require further attention in model development.