|
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).
|
|
|
Webber, H., White, J. W., Kimball, B. A., Ewert, F., Asseng, S., Rezaei, E. E., et al. (2018). Physical robustness of canopy temperature models for crop heat stress simulation across environments and production conditions. Field Crops Research, 216, 75–88.
Abstract: Despite widespread application in studying climate change impacts, most crop models ignore complex interactions among air temperature, crop and soil water status, CO2 concentration and atmospheric conditions that influence crop canopy temperature. The current study extended previous studies by evaluating Tc simulations from nine crop models at six locations across environmental and production conditions. Each crop model implemented one of an empirical (EMP), an energy balance assuming neutral stability (EBN) or an energy balance correcting for atmospheric stability conditions (EBSC) approach to simulate Tc. Model performance in predicting Tc was evaluated for two experiments in continental North America with various water, nitrogen and CO2 treatments. An empirical model fit to one dataset had the best performance, followed by the EBSC models. Stability conditions explained much of the differences between modeling approaches. More accurate simulation of heat stress will likely require use of energy balance approaches that consider atmospheric stability conditions.
|
|
|
Persson, T., Höglind, M., Gustavsson, A. - M., Halling, M., Jauhiainen, L., Niemeläinen, O., et al. (2014). Evaluation of the LINGRA timothy model under Nordic conditions. Field Crops Research, 161, 87–97.
Abstract: Simulation models are frequently applied to determine the production potential of forage grasses under various scenarios, including climate change. Thorough calibrations and evaluations of forage grass models can help improve their applicability. This study evaluated the ability of the Light Interception and Utilization Simulator-GRAss (LINGRA) model to predict biomass yield of timothy (Phleum pratense L. cv. Grindstad) in the Nordic countries. Variety trial data for the first and second year after establishment were obtained for seven locations: Jokioinen, Finland (60 degrees 48 ‘ N; 23 degrees 29 ‘ E), Maaninka, Finland (63 degrees 09 ‘ N; 27 degrees 18 ‘ E), Korpa, Iceland (64 degrees 09 ‘ N; 21 degrees 45 ‘ W), Srheim, Norway (58 degrees 41 ‘ N; 5 degrees 39 ‘ E), Lillerud, Sweden (59 degrees 24’ N; 13 degrees 16 ‘ E), Ostersund, Sweden (63 degrees 15 ‘ N; 14 degrees 34 ‘ E) and Ulna Sweden (63 degrees 49 ‘ N; 20 degrees 13 ‘ E) from 1992 to 2012. Two calibrations of the LINGRA model were carried out using Bayesian techniques. In the first of these (SRrheim calibration), data on biomass yield and underlying variables obtained from independent field trials at Srheim were used. In the second (Nordic calibration), biomass data from the other locations were used as well. The model was validated against the remaining set of biomass yields from all locations not included in the Nordic calibration. The observed total seasonal yield the first and second year after establishment was 913 and 991 g DM m(-2) respectively on average across the locations. The corresponding average simulated yield after the Srheim calibration was 1044 (root mean square error (RMSE) 258) and 1112 g DM m(-2) (RMSE 312), respectively. After the Nordic calibration, the simulated average total seasonal yield was 863 (RMSE 242) the first year and 927 g DM m(-2) (RMSE 271) the second year after establishment. The differences between the observed and simulated first cut yield followed the same patterns, whereas the prediction accuracy for second cut yield did not differ substantially between the calibration approaches.Using the parameter set from the Nordic region decreased the model predictability at Srheim compared with only using model parameters derived from this location. These results show that using biomass data from several locations, instead of only one specific location, in the calibration of the LINGRA model improved the overall prediction accuracy of first cut dry matter yield and total seasonal dry matter yield across an environmentally heterogeneous region. To further analyse the usefulness of including multi-site data in forage grass model calibrations, other forage grass models could be evaluated against the same dataset.
|
|
|
Grosz, B., Dechow, R., Gebbert, S., Hoffmann, H., Zhao, G., Constantin, J., et al. (2017). The implication of input data aggregation on up-scaling soil organic carbon changes. Env. Model. Softw., 96, 361–377.
Abstract: In up-scaling studies, model input data aggregation is a common method to cope with deficient data availability and limit the computational effort. We analyzed model errors due to soil data aggregation for modeled SOC trends. For a region in North West Germany, gridded soil data of spatial resolutions between 1 km and 100 km has been derived by majority selection. This data was used to simulate changes in SOC for a period of 30 years by 7 biogeochemical models. Soil data aggregation strongly affected modeled SOC trends. Prediction errors of simulated SOC changes decreased with increasing spatial resolution of model output. Output data aggregation only marginally reduced differences of model outputs between models indicating that errors caused by deficient model structure are likely to persist even if requirements on the spatial resolution of model outputs are low. (C)2017 Elsevier Ltd. All rights reserved.
|
|
|
Zhang, S., Tao, F., & Zhang, Z. (2017). Uncertainty from model structure is larger than that from model parameters in simulating rice phenology in China. Europ. J. Agron., 87, 30–39.
Abstract: Rice models have been widely used in simulating and predicting rice phenology in contrasting climate zones, however the uncertainties from model structure (different equations or models) and/or model parameters were rarely investigated. Here, five rice phenological models/modules (Le., CERES-Rice, ORYZA2000, RCM, Beta Model and SIMRIW) were applied to simulate rice phenology at 23 experimental stations from 1992 to 2009 in two major rice cultivation regions of China: the northeastern China and the southwestern China. To investigate the uncertainties from model biophysical parameters, each model was run with randomly perturbed 50 sets of parameters. The results showed that the median of ensemble simulations were better than the simulation by most models. Models couldn’t simulate well in some specific years despite of parameters optimization, suggesting model structure limit model performance in some cases. The models adopting accumulative thermal time function (e.g., CERES-Rice and ORYZA2000) had better performance in the southwestern China, in contrast, those adopting exponential function (e.g., Beta model and RCM model) had better performance in the northeastern China. In northeastern China, the contribution of model structure and model parameters to model total variance was, respectively, about 55.90% and 44.10% in simulating heading date, and about 75.43% and 24.57% in simulating maturity date. In the southwestern China, the contribution of model structure and model parameters to model total variance was, respectively, about 79.97% and 27.03% in simulating heading date, about 92.15% and 7.85% in simulating maturity date. Uncertainty from model structure was the most relevant source. The results highlight that the temperature response functions of rice development rate under extreme climate conditions should be improved based on environment-controlled experimental data.
|
|