|
Dumont, B., Leemans, V., Mansouri, M., Bodson, B., Destain, J. - P., & Destain, M. - F. (2014). Parameter identification of the STICS crop model, using an accelerated formal MCMC approach. Env. Model. Softw., 52, 121–135.
Abstract: This study presents a Bayesian approach for the parameters’ identification of the STICS crop model based on the recently developed Differential Evolution Adaptive Metropolis (DREAM) algorithm. The posterior distributions of nine specific crop parameters of the STICS model were sampled with the aim to improve the growth simulations of a winter wheat (Triticum aestivum L) culture. The results obtained with the DREAM algorithm were initially compared to those obtained with a Nelder-Mead Simplex algorithm embedded within the OptimiSTICS package. Then, three types of likelihood functions implemented within the DREAM algorithm were compared, namely the standard least square, the weighted least square, and a transformed likelihood function that makes explicit use of the coefficient of variation (CV). The results showed that the proposed CV likelihood function allowed taking into account both noise on measurements and heteroscedasticity which are regularly encountered in crop modelling. (C) 2013 Elsevier Ltd. All rights reserved.
|
|
|
Dumont, B., Basso, B., Bodson, B., Destain, J. - P., & Destain, M. - F. (2015). Climatic risk assessment to improve nitrogen fertilisation recommendations: A strategic crop model-based approach. European Journal of Agronomy, 65, 10–17.
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.
|
|
|
Ma, S., Lardy, R., Graux, A. - I., Ben Touhami, H., Klumpp, K., Martin, R., et al. (2015). Regional-scale analysis of carbon and water cycles on managed grassland systems. Env. Model. Softw., 72, 356–371.
Abstract: Predicting regional and global carbon (C) and water dynamics on grasslands has become of major interest, as grasslands are one of the most widespread vegetation types worldwide, providing a number of ecosystem services (such as forage production and C storage). The present study is a contribution to a regional-scale analysis of the C and water cycles on managed grasslands. The mechanistic biogeochemical model PaSim (Pasture Simulation model) was evaluated at 12 grassland sites in Europe. A new parameterization was obtained on a common set of eco-physiological parameters, which represented an improvement of previous parameterization schemes (essentially obtained via calibration at specific sites). We found that C and water fluxes estimated with the parameter set are in good agreement with observations. The model with the new parameters estimated that European grassland are a sink of C with 213 g C m(-2) yr(-1), which is close to the observed net ecosystem exchange (NEE) flux of the studied sites (185 g C m(-2) yr(-1) on average). The estimated yearly average gross primary productivity (GPP) and ecosystem respiration (RECO) for all of the study sites are 1220 and 1006 g C m(-2) yr(-1), respectively, in agreement with observed average GPP (1230 g C m(-2) yr(-1)) and RECO (1046 g C m(-2) yr(-1)). For both variables aggregated on a weekly basis, the root mean square error (RMSE) was similar to 5-16 g C week(-1) across the study sites, while the goodness of fit (R-2) was similar to 0.4-0.9. For evapotranspiration (ET), the average value of simulated ET (415 mmyr(-1)) for all sites and years is close to the average value of the observed ET (451 mm yr(-1)) by flux towers (on a weekly basis, RMSE similar to 2-8 mm week(-1); R-2 = 0.3-0.9). However, further model development is needed to better represent soil water dynamics under dry conditions and soil temperature in winter. A quantification of the uncertainties introduced by spatially generalized parameter values in C and water exchange estimates is also necessary. In addition, some uncertainties in the input management data call for the need to improve the quality of the observational system.
|
|
|
Francone, C., Cassardo, C., Richiardone, R., & Confalonieri, R. (2012). Sensitivity Analysis and Investigation of the Behaviour of the UTOPIA Land-Surface Process Model: A Case Study for Vineyards in Northern Italy. Boundary-Layer Meteorology, 144(3), 419–430.
Abstract: We used sensitivity-analysis techniques to investigate the behaviour of the land-surface model UTOPIA while simulating the micrometeorology of a typical northern Italy vineyard (Vitis vinifera L.) under average climatic conditions. Sensitivity-analysis experiments were performed by sampling the vegetation parameter hyperspace using the Morris method and quantifying the parameter relevance across a wide range of soil conditions. This method was used since it proved its suitability for models with high computational time or with a large number of parameters, in a variety of studies performed on different types of biophysical models. The impact of input variability was estimated on reference model variables selected among energy (e.g. net radiation, sensible and latent heat fluxes) and hydrological (e.g. soilmoisture, surface runoff, drainage) budget components. Maximum vegetation cover and maximum leaf area index were ranked as the most relevant parameters, with sensitivity indices exceeding the remaining parameters by about one order of magnitude. Soil variability had a high impact on the relevance of most of the vegetation parameters: coefficients of variation calculated on the sensitivity indices estimated for the different soils often exceeded 100 %. The only exceptions were represented by maximum vegetation cover and maximum leaf area index, which showed a low variability in sensitivity indices while changing soil type, and confirmed their key role in affecting model results.
|
|
|
Dumont, B., Basso, B., Leemans, V., Bodson, B., Destain, J. - P., & Destain, M. - F. (2015). A comparison of within-season yield prediction algorithms based on crop model behaviour analysis. Agricultural and Forest Meteorology, 204, 10–21.
Abstract: The development of methodologies for predicting crop yield, in real-time and in response to different agro-climatic conditions, could help to improve the farm management decision process by providing an analysis of expected yields in relation to the costs of investment in particular practices. Based on the use of crop models, this paper compares the ability of two methodologies to predict wheat yield (Triticum aestivum L.), one based on stochastically generated climatic data and the other on mean climate data. It was shown that the numerical experimental yield distribution could be considered as a log-normal distribution. This function is representative of the overall model behaviour. The lack of statistical differences between the numerical realisations and the logistic curve showed in turn that the Generalised Central Limit Theorem (GCLT) was applicable to our case study. In addition, the predictions obtained using both climatic inputs were found to be similar at the inter and intra-annual time-steps, with the root mean square and normalised deviation values below an acceptable level of 10% in 90% of the climatic situations. The predictive observed lead-times were also similar for both approaches. Given (i) the mathematical formulation of crop models, (ii) the applicability of the CLT and GLTC to the climatic inputs and model outputs, respectively, and (iii) the equivalence of the predictive abilities, it could be concluded that the two methodologies were equally valid in terms of yield prediction. These observations indicated that the Convergence in Law Theorem was applicable in this case study. For purely predictive purposes, the findings favoured an algorithm based on a mean climate approach, which needed far less time (by 300-fold) to run and converge on same predictive lead time than the stochastic approach. (C) 2015 Elsevier B.V. All rights reserved.
|
|