|
|
Makowski, D. (2017). A simple Bayesian method for adjusting ensemble of crop model outputs to yield observations. Europ. J. Agron., 88, 76–83.
Abstract: Multi-model forecasting has drawn some attention in crop science for evaluating effect of climate change on crop yields. The principle is to run several individual process-based crop models under several climate scenarios in order to generate ensembles of output values. This paper describes a simple Bayesian method – called Bayes linear method- for updating ensemble of crop model outputs using yield observations. The principle is to summarize the ensemble of crop model outputs by its mean and variance, and then to adjust these two quantities to yield observations in order to reduce uncertainty. The adjusted mean and variance combine two sources of information, i.e., the ensemble of crop model outputs and the observations. Interestingly, with this method, observations collected under a given climate scenario can be used to adjust mean and variance of the model ensemble under a different scenario. Another advantage of the proposed method is that it does not rely on a separate calibration of each individual crop model. The uncertainty reduction resulting from the adjustment of an ensemble of crop models to observations was assessed in a numerical application. The implementation of the Bayes linear method systematically reduced uncertainty, but the results showed the effectiveness of this method varied in function of several factors, especially the accuracy of the yield observation, and the covariance between the crop model output and the observation. (C) 2015 Elsevier B.V. All rights reserved.
|
|
|
|
Sandor, R., Ehrhardt, F., Grace, P., Recous, S., Smith, P., Snow, V., et al. (2020). Ensemble modelling of carbon fluxes in grasslands and croplands. Field Crops Research, 252, 107791.
Abstract: Croplands and grasslands are agricultural systems that contribute to land–atmosphere exchanges of carbon (C). We evaluated and compared gross primary production (GPP), ecosystem respiration (RECO), net ecosystem exchange (NEE) of CO2, and two derived outputs – C use efficiency (CUE=-NEE/GPP) and C emission intensity (IntC= -NEE/Offtake [grazed or harvested biomass]). The outputs came from 23 models (11 crop-specific, eight grassland-specific, and four models covering both systems) at three cropping sites over several rotations with spring and winter cereals, soybean and rapeseed in Canada, France and India, and two temperate permanent grasslands in France and the United Kingdom. The models were run independently over multi-year simulation periods in five stages (S), either blind with no calibration and initialization data (S1), using historical management and climate for initialization (S2), calibrated against plant data (S3), plant and soil data together (S4), or with the addition of C and N fluxes (S5). Here, we provide a framework to address methodological uncertainties and contextualize results. Most of the models overestimated or underestimated the C fluxes observed during the growing seasons (or the whole years for grasslands), with substantial differences between models. For each simulated variable, changes in the multi-model median (MMM) from S1 to S5 was used as a descriptor of the ensemble performance. Overall, the greatest improvements (MMM approaching the mean of observations) were achieved at S3 or higher calibration stages. For instance, grassland GPP MMM was equal to 1632 g C m−2 yr-1 (S5) while the observed mean was equal to 1763 m-2 yr-1 (average for two sites). Nash-Sutcliffe modelling efficiency coefficients indicated that MMM outperformed individual models in 92.3 % of cases. Our study suggests a cautious use of large-scale, multi-model ensembles to estimate C fluxes in agricultural sites if some site-specific plant and soil observations are available for model calibration. The further development of crop/grassland ensemble modelling will hinge upon the interpretation of results in light of the way models represent the processes underlying C fluxes in complex agricultural systems (grassland and crop rotations including fallow periods).
|
|
|
|
Yin, X. G., Kersebaum, K. C., Kollas, C., Manevski, K., Baby, S., Beaudoin, N., et al. (2017). Performance of process-based models for simulation of grain N in crop rotations across Europe. Agric. Syst., 154, 63–77.
Abstract: The accurate estimation of crop grain nitrogen (N; N in grain yield) is crucial for optimizing agricultural N management, especially in crop rotations. In the present study, 12 process-based models were applied to simulate the grain N of i) seven crops in rotations, ii) across various pedo-climatic and agro-management conditions in Europe, under both continuous simulation and single year simulation, and for iv) two calibration levels, namely minimal and detailed calibration. Generally, the results showed that the accuracy of the simulations in predicting grain N increased under detailed calibration. The models performed better in predicting the grain N of winter wheat (Triticum aestivum L.), winter barley (Hordewn vulgare L.) and spring barley (Hordeum vulgare L.) compared to spring oat (Avena saliva L.), winter rye (Secale cereale L.), pea (Piswn sativum L.) and winter oilseed rape (Brassica napus L.). These differences are linked to the intensity of parameterization with better parameterized crops showing lower prediction errors. The model performance was influenced by N fertilization and irrigation treatments, and a majority of the predictions were more accurate under low N and rainfed treatments. Moreover, the multi-model mean provided better predictions of grain N compared to any individual model. In regard to the Individual models, DAISY, FASSET, HERMES, MONICA and STICS are suitable for predicting grain N of the main crops in typical European crop rotations, which all performed well in both continuous simulation and single year simulation. Our results show that both the model initialization and the cover crop effects in crop rotations should be considered in order to achieve good performance of continuous simulation. Furthermore, the choice of either continuous simulation or single year simulation should be guided by the simulation objectives (e.g. grain yield, grain N content or N dynamics), the crop sequence (inclusion of legumes) and treatments (rate and type of N fertilizer) included in crop rotations and the model formalism.
|
|
|
|
Webber, H., Ewert, F., Kimball, B. A., Siebert, S., White, J. W., Wall, G. W., et al. (2016). Simulating canopy temperature for modelling heat stress in cereals. Env. Model. Softw., 77, 143–155.
Abstract: Crop models must be improved to account for the effects of heat stress events on crop yields. To date, most approaches in crop models use air temperature to define heat stress intensity as the cumulative sum of thermal times (TT) above a high temperature threshold during a sensitive period for yield formation. However, observational evidence indicates that crop canopy temperature better explains yield reductions associated with high temperature events than air temperature does. This study presents a canopy level energy balance using Monin ObukhovSimilarity Theory (MOST) with simplifications about the canopy resistance that render it suitable for application in crop models and other models of the plant environment. The model is evaluated for a uniform irrigated wheat canopy in Arizona and rainfed maize in Burkina Faso. No single variable regression relationships for key explanatory variables were found that were consistent across sowing dates to explain the deviation of canopy temperature from air temperature. Finally, thermal times determined with simulated canopy temperatures were able to reproduce thermal times calculated with observed canopy temperature, whereas those determined with air temperatures were not. (C) 2015 Elsevier Ltd. All rights reserved.
|
|
|
|
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.
|
|
