Abstract: This paper describes the environmental consequences of agricultural adaptation on eutrophication of the nearby ecological network for a study area in the Netherlands. More specifically, we explored (i) likely responses of farmers to changes in climate, technology, policy, and markets; (ii) subsequent changes in nitrogen (N) emissions in responses to farmer adaptations; and (iii) to what extent the emitted N was deposited in nearby nature reserves, in view of the potential impacts on plant species diversity and desired nature targets. For this purpose, a spatially-explicit study at landscape level was performed by integrating the environmental model INITIATOR, the farm model FSSIM, and the land-use model RULEX. We evaluated two alternative scenarios of change in climate, technology, policy, and markets for 2050: one in line with a ‘global economy’ (GE) storyline and the other in line with a ‘regional communities’ (RC) storyline. Results show that the GE storyline resulted in a relatively strong increase in agricultural production compared to the RC storyline. Despite the projected conversions of agricultural land to nature (as part of the implementation of the National Ecological Network), we project an increase in N losses and N deposition due to N emissions in the study area of about 20 %. Even in the RC storyline, with a relatively modest increase in agricultural production and a larger expansion of the nature reserve, the N losses and deposition remain at the current level, whereas a reduction is required. We conclude that more ambitious green policies are needed in view of nature protection.

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.