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: The size of arable crop yields depends on many weather factors, such as precipitation and air temperature during the vegetation period. When studying the relation between yields and precipitation, not only the total amount of precipitation, but also the occurrence of long periods without precipitation must be taken into account. The paper [Bojar et al., 2014] demonstrated that barley yield significantly statistically depends on the length of the series of days without precipitation. This paper attempts to analyse the statistical data on daily precipitation totals recorded during the January – December periods in the years 1971 – 2013 at the weather station of the University of Science and Technology in Bydgoszcz, Faculty of Agriculture and Biotechnology, in the Research Centre located in an agricultural area in the Mochle township, situated 17 kilometres from Bydgoszcz. The primary statistical operation in the study is an attempt to estimate the Markov chain order. To this end, two criteria of chain order determination are applied: BIC (Bayesian information criterion, Schwarz 1978) and AIC (Akaike information criterion, Akaike 1974). Both are based on the log-likelihood functions for transition probability of the Markov chain constructed on certain data series. Statistical analysis of precipitation totals data leads to the conclusion that both AIC and BIC indicate the 2nd order for the studied Markov chain. The proposed method of estimating the variability of precipitation occurrence in the future will be utilised to improve region-related bio-physical and economical models, and to assess the risk of extreme events in the context of growing climate hazards. It will serve as basis for a search in agriculture for solutions mitigating those hazards.