|
|
Minet, J., Laloy, E., Tychon, B., & François, L. (2015). Bayesian inversions of a dynamic vegetation model at four European grassland sites. Biogeosciences, 12(9), 2809–2829.
Abstract: Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM((ZS)) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m(-2) day(-1) and 0.50 to 1.28 mm day(-1), respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic mortality time parameters can be explained by site-specific ecophysiological characteristics.
|
|
|
|
Constantin, J., Raynal, H., Casellas, E., Hoffman, H., Bindi, M., Doro, L., et al. (2019). Management and spatial resolution effects on yield and water balance at regional scale in crop models. Agricultural and Forest Meteorology, 275, 184–195.
Abstract: Due to the more frequent use of crop models at regional and national scale, the effects of spatial data input resolution have gained increased attention. However, little is known about the influence of variability in crop management on model outputs. A constant and uniform crop management is often considered over the simulated area and period. This study determines the influence of crop management adapted to climatic conditions and input data resolution on regional-scale outputs of crop models. For this purpose, winter wheat and maize were simulated over 30 years with spatially and temporally uniform management or adaptive management for North Rhine-Westphalia ((similar to)34 083 km(2)), Germany. Adaptive management to local climatic conditions was used for 1) sowing date, 2) N fertilization dates, 3) N amounts, and 4) crop cycle length. Therefore, the models were applied with four different management sets for each crop. Input data for climate, soil and management were selected at five resolutions, from 1 x 1 km to 100 x 100 km grid size. Overall, 11 crop models were used to predict regional mean crop yield, actual evapotranspiration, and drainage. Adaptive management had little effect (< 10% difference) on the 30-year mean of the three output variables for most models and did not depend on soil, climate, and management resolution. Nevertheless, the effect was substantial for certain models, up to 31% on yield, 27% on evapotranspiration, and 12% on drainage compared to the uniform management reference. In general, effects were stronger on yield than on evapotranspiration and drainage, which had little sensitivity to changes in management. Scaling effects were generally lower than management effects on yield and evapotranspiration as opposed to drainage. Despite this trend, sensitivity to management and scaling varied greatly among the models. At the annual scale, effects were stronger in certain years, particularly the management effect on yield. These results imply that depending on the model, the representation of management should be carefully chosen, particularly when simulating yields and for predictions on annual scale.
|
|
|
|
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).
|
|
|
|
Angulo, C., Rötter, R., Trnka, M., Pirttioja, N., Gaiser, T., Hlavinka, P., et al. (2013). Characteristic ‘fingerprints’ of crop model responses to weather input data at different spatial resolutions. European Journal of Agronomy, 49, 104–114.
Abstract: Crop growth simulation models are increasingly used for regionally assessing the effects of climate change and variability on crop yields. These models require spatially and temporally detailed, location-specific, environmental (weather and soil) and management data as inputs, which are often difficult to obtain consistently for larger regions. Aggregating the resolution of input data for crop model applications may increase the uncertainty of simulations to an extent that is not well understood. The present study aims to systematically analyse the effect of changes in the spatial resolution of weather input data on yields simulated by four crop models (LINTUL-SLIM, DSSAT-CSM, EPIC and WOFOST) which were utilized to test possible interactions between weather input data resolution and specific modelling approaches representing different degrees of complexity. The models were applied to simulate grain yield of spring barley in Finland for 12 years between 1994 and 2005 considering five spatial resolutions of daily weather data: weather station (point) and grid-based interpolated data at resolutions of 10 km x 10 km; 20 km x 20 km; 50 km x 50 km and 100 km x 100 km. Our results show that the differences between models were larger than the effect of the chosen spatial resolution of weather data for the considered years and region. When displaying model results graphically, each model exhibits a characteristic ‘fingerprint’ of simulated yield frequency distributions. These characteristic distributions in response to the inter-annual weather variability were independent of the spatial resolution of weather input data. Using one model (LINTUL-SLIM), we analysed how the aggregation strategy, i.e. aggregating model input versus model output data, influences the simulated yield frequency distribution. Results show that aggregating weather data has a smaller effect on the yield distribution than aggregating simulated yields which causes a deformation of the model fingerprint. We conclude that changes in the spatial resolution of weather input data introduce less uncertainty to the simulations than the use of different crop models but that more evaluation is required for other regions with a higher spatial heterogeneity in weather conditions, and for other input data related to soil and crop management to substantiate our findings. Our results provide further evidence to support other studies stressing the importance of using not just one, but different crop models in climate assessment studies. (C) 2013 Elsevier B.V. All rights reserved.
|
|
|
|
Kim, Y., Seo, Y., Kraus, D., Klatt, S., Haas, E., Tenhunen, J., et al. (2015). Estimation and mitigation of N₂O emission and nitrate leaching from intensive crop cultivation in the Haean catchment, South Korea. Science of the Total Environment, 529, 40–53.
Abstract: Considering intensive agricultural management practices and environmental conditions, the LandscapeDNDC model was applied for simulation of yields, N2O emission and nitrate leaching from major upland crops and temperate deciduous forest of the Haean catchment, South Korea. Fertilization rates were high (up to 314 kg N ha(-1) year(-1)) and resulted in simulated direct N2O emissions from potato, radish, soybean and cabbage fields of 1.9 and 2.1 kg N ha(-1) year(-1) in 2009 and 2010, respectively. Nitrate leaching was identified as the dominant pathway of N losses in the Haean catchment with mean annual rates of 112.2 and 125.4 kg N ha(-1) year(-1), causing threats to water quality and leading to substantial indirect N2O emissions of 0.84 and 0.94 kg N ha(-1) year(-1) in 2009 and 2010 as estimates by applying the IPCC EF5. Simulated N2O emissions from temperate deciduous forest were low (approx. 0.50 kg N ha(-1) year(-1)) and predicted nitrate leaching rates were even negligible (≤0.01 kg N ha(-1) year(-1)). On catchment scale more than 50% of the total N2O emissions and up to 75% of nitrate leaching originated from fertilized upland fields, only covering 24% of the catchment area. Taking into account area coverage of simulated upland crops and other land uses these numbers agree well with nitrate loads calculated from discharge and concentration measurements at the catchment outlet. The change of current agricultural management practices showed a high potential of reducing N2O emission and nitrate leaching while maintaining current crop yields. Reducing (39%) and splitting N fertilizer application into 3 times was most effective and lead to about 54% and 77% reducing of N2O emission and nitrate leaching from the Haean catchment, the latter potentially contributing to improved water quality in the Soyang River Dam, which is the major source of drinking water for metropolitan residents.
|
|
