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: Apart from measurements, numerical models are the most convenient instruments to analyze the carbon and water balance of terrestrial ecosystems and their interactions with changing environmental conditions. The process-based Biome-BGC model is widely used to simulate the storage and flux of water, carbon, and nitrogen within the vegetation, litter, and soil of unmanaged terrestrial ecosystems. Considering herbaceous vegetation related simulations with Biome-BGC, soil moisture and growing season control on ecosystem functioning is inaccurate due to the simple soil hydrology and plant phenology representation within the model. Consequently, Biome-BGC has limited applicability in herbaceous ecosystems because (1) they are usually managed; (2) they are sensitive to soil processes, most of all hydrology; and (3) their carbon balance is closely connected with the growing season length. Our aim was to improve the applicability of Biome-BGC for managed herbaceous ecosystems by implementing several new modules, including management. A new index (heatsum growing season index) was defined to accurately estimate the first and the final days of the growing season. Instead of a simple bucket soil sub-model, a multilayer soil sub-model was implemented, which can handle the processes of runoff, diffusion and percolation. A new module was implemented to simulate the ecophysiological effect of drought stress on plant mortality. Mowing and grazing modules were integrated in order to quantify the functioning of managed ecosystems. After modifications, the Biome-BGC model was calibrated and validated using eddy covariance-based measurement data collected in Hungarian managed grassland ecosystems. Model calibration was performed based on the Bayes theorem. As a result of these developments and calibration, the performance of the model was substantially improved. Comparison with measurement-based estimate showed that the start and the end of the growing season are now predicted with an average accuracy of 5 and 4 days instead of 46 and 85 days as in the original model. Regarding the different sites and modeled fluxes (gross primary production, total ecosystem respiration, evapotranspiration), relative errors were between 18-60% using the original model and 10-18% using the developed model; squares of the correlation coefficients were between 0.02-0.49 using the original model and 0.50-0.81 using the developed model. (c) 2011 Elsevier B.V. All rights reserved.