Abstract: The main aim of the paper is to show how, and how many, simulations carried out using the Land Surface Model UTOPIA (University of TOrino model of land Process Interaction with Atmosphere) are representative of the micro-meteorological conditions and exchange processes at the atmosphere/biosphere interface, with a particular focus on heat and hydrologic transfers, over an area of the Piemonte (Piedmont) region, NW Italy, which is characterized by the presence of many vineyards. Another equally important aim is to understand how much the quality of the simulation outputs was influenced by the input data, whose measurements are often unavailable for long periods over country areas at an hourly basis. Three types of forcing data were used: observations from an experimental campaign carried out during the 2008, 2009, and 2010 vegetative seasons in three vineyards, and values extracted from the freely available Global Land Data Assimilation System (GLDAS, versions 2.0 and 2.1). Since GLDAS also contains the outputs of the simulations performed using the Land Surface Model NOAH, an additional intercomparison between the two models, UTOPIA and NOAH, both driven by the same GLDAS datasets, was performed. The intercomparisons were performed on the following micro-meteorological variables: net radiation, sensible and latent turbulent heat fluxes, and temperature and humidity of soil. The results of this study indicate that the methodology of employing land surface models driven by a gridded database to evaluate variables of micro-meteorological and agronomic interest in the absence of observations is suitable and gives satisfactory results, with uncertainties comparable to measurement errors, thus, allowing us to also evaluate some time trends. The comparison between GLDAS2.0 and GLDAS2.1 indicates that the latter generally produces simulation outputs more similar to the observations than the former, using both UTOPIA and NOAH models.

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.