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.

Abstract: We assessed the weather data aggregation effect (DAE) on the simulation of cropping systems for different crops, response variables, and production conditions. Using 13 process-based crop models and the ensemble mean, we simulated 30 yr continuous cropping systems for 2 crops (winter wheat and silage maize) under 3 production conditions for the state of North Rhine-Westphalia, Germany. The DAE was evaluated for 5 weather data resolutions (i.e. 1, 10, 25, 50, and 100 km) for 3 response variables including yield, growing season evapotranspiration, and water use efficiency. Five metrics, viz. the spatial bias (Delta), average absolute deviation (AAD), relative AAD, root mean squared error (RMSE), and relative RMSE, were used to evaluate the DAE on both the input weather data and simulated results. For weather data, we found that data aggregation narrowed the spatial variability but widened the., especially across mountainous areas. The DAE on loss of spatial heterogeneity and hotspots was stronger than on the average changes over the region. The DAE increased when coarsening the spatial resolution of the input weather data. The DAE varied considerably across different models, but changed only slightly for different production conditions and crops. We conclude that if spatially detailed information is essential for local management decision, higher resolution is desirable to adequately capture the spatial variability for heterogeneous regions. The required resolution depends on the choice of the model as well as the environmental condition of the study area.

Abstract: Continuous crop rotation modeling is a prospective trend that, compared to 1-crop or discrete year-by-year calculations, can provide more accurate results that are closer to real conditions. The goal of this study was to compare the water balance and yields estimated by the HERMES crop rotation model for present and future climatic conditions in the Czech Republic. Three locations were selected, representing important agricultural regions with different climatic conditions. Crop rotation (spring barley, silage maize, winter wheat, winter rape) was simulated from 1981-2080. The 1981-2010 period was covered by measured meteorological data, while 2011-2080 was represented by a transient synthetic weather series from the weather generator M& Rfi. The data were based on 5 circulation models, representing an ensemble of 18 CMIP3 global circulation models, to preserve much of the uncertainty of the original ensemble. Two types of crop management were compared, and the influences of soil quality, increasing atmospheric CO2 and adaptation measures (i. e. sowing date changes) were also considered. Results suggest that under a ‘dry’ scenario (such as GFCM21), C-3 crops in drier regions will be devastated for a significant number of seasons. Negative impacts are likely even on premium-quality soils regardless of flexible sowing dates and accounting for increasing CO2 concentrations. Moreover, in dry conditions, the use of crop rotations with catch crops may have negative impacts, exacerbating the soil water deficit for subsequent crops. This approach is a promising method for determining how various management strategies and crop rotations can affect yields as well as water, carbon and nitrogen cycling.