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: Two independently developed simulation models – the grassland-specific PaSim and the biome-generic Biome-BGC MuSo (BBGC MuSo) – linking climate, soil, vegetation and management to ecosystem biogeochemical cycles were compared in a simulation of carbon (C) and water fluxes. The results were assessed against eddy-covariance flux data from five observational grassland sites representing a range of conditions in Europe: Grillenburg in Germany, Laqueuille in France with both extensive and intensive management, Monte Bondone in Italy and Oensingen in Switzerland. Model comparison (after calibration) gave substantial agreement, the performances being marginal to acceptable for weekly-aggregated gross primary production and ecosystem respiration (R-2 similar to 0.66 – 0.91), weekly evapotranspiration (R-2 similar to 0.78 – 0.94), soil water content in the topsoil (R-2 similar to 0.1 -0.7) and soil temperature (R-2 similar to 0.88 – 0.96). The bias was limited to the range -13 to 9 g C m(-2) week(-1) for C fluxes (-11 to 8 g C m(-2) week(-1) in case of BBGC MuSo, and -13 to 9 g C m(-2) week(-1) in case of PaSim) and -4 to 6 mm week for water fluxes (with BBGC MuSo providing somewhat higher estimates than PaSim), but some higher relative root mean square errors indicate low accuracy for prediction, especially for net ecosystem exchange The sensitivity of simulated outputs to changes in atmospheric carbon dioxide concentration ([CO2]), temperature and precipitation indicate, with certain agreement between the two models, that C outcomes are dominated by [CO2] and temperature gradients, and are less due to precipitation. ET rates decrease with increasing [CO2] in PaSim (consistent with experimental knowledge), while lack of appropriate stomatal response could be a limit in BBGC MuSo responsiveness. Results of the study indicate that some of the errors might be related to the improper representation of soil water content and soil temperature. Improvement is needed in the model representations of soil processes (especially soil water balance) that strongly influence the biogeochemical cycles of managed and unmanaged grasslands. (C) 2015 Elsevier B.V. All rights reserved.