|
|
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.
|
|
|
|
Minet, J., Laloy, E., Tychon, B., & François, L. (2014). Bayesian inference of a dynamic vegetation model for grassland. FACCE MACSUR Mid-term Scientific Conference, 3(S) Sassari, Italy.
Abstract: As a part of the MACSUR task L2.4, we probabilistically calibrated the CARAIB dynamic vegetation model by Markov chain Monte Carlo (MCMC) simulation with the DREAMZS sampler. CARAIB is a mechanistic model that calculates the carbon assimilation of the vegetation as a function of the soil and climatic conditions, and can thus be used for simulating grassland production under cutting or grazing management. Bayesian model inversion was performed at 4 grassland sites across Europe: Oensingen, CH; Grillenburg, DE; Laqueuille, FR and Monte-Bodone, IT. Four daily measured variables from these sites: the Gross Primary Productivity (GPP), Net Ecosystem Exchange (NEE), Evapotranspiration (ET) and Soil Water Content (SWC) were used to sample 10 parameters related to rooting depth, stomatal conductance, specific leaf area, carbon-nitrogen ratio and water stresses. The maximized likelihood function therefore involved four objectives, whereas the applied Bayesian framework allowed for assessing the so called parameter posterior probability density function (pdf), which quantifies model parameter uncertainty caused by measurement and model errors. Sampling trials were performed using merged data from all sites (all-sites-sampling) and for each site (site-specific sampling) separately. The derived posterior parameter pdfs from the all-sites sampling and site-specific sampling runs showed differences in relation with the specificities of each site. Analysis of these distributions also revealed model sensitivity to parameters conditioned on the measured data, as well as parameter correlations.
|
|
|
|
Minet, J., Tychon, B., Jacquemin, I., & François, L. (2014). Can a global dynamic vegetation model be used for both grassland and crop modeling at the local scale..
|
|
|
|
Pirttioja, N., Carter, T. R., Fronzek, S., Bindi, M., Hoffmann, H., Palosuo, T., et al. (2015). Temperature and precipitation effects on wheat yield across a European transect: a crop model ensemble analysis using impact response surfaces. Clim. Res., 65, 87–105.
Abstract: This study explored the utility of the impact response surface (IRS) approach for investigating model ensemble crop yield responses under a large range of changes in climate. IRSs of spring and winter wheat Triticum aestivum yields were constructed from a 26-member ensemble of process-based crop simulation models for sites in Finland, Germany and Spain across a latitudinal transect. The sensitivity of modelled yield to systematic increments of changes in temperature (-2 to +9°C) and precipitation (-50 to +50%) was tested by modifying values of baseline (1981 to 2010) daily weather, with CO2 concentration fixed at 360 ppm. The IRS approach offers an effective method of portraying model behaviour under changing climate as well as advantages for analysing, comparing and presenting results from multi-model ensemble simulations. Though individual model behaviour occasionally departed markedly from the average, ensemble median responses across sites and crop varieties indicated that yields decline with higher temperatures and decreased precipitation and increase with higher precipitation. Across the uncertainty ranges defined for the IRSs, yields were more sensitive to temperature than precipitation changes at the Finnish site while sensitivities were mixed at the German and Spanish sites. Precipitation effects diminished under higher temperature changes. While the bivariate and multi-model characteristics of the analysis impose some limits to interpretation, the IRS approach nonetheless provides additional insights into sensitivities to inter-model and inter-annual variability. Taken together, these sensitivities may help to pinpoint processes such as heat stress, vernalisation or drought effects requiring refinement in future model development.
|
|
|
|
Rodriguez, A., Ruiz-Ramos, M., Palosuo, T., Carter, T. R., Fronzek, S., Lorite, I. J., et al. (2019). Implications of crop model ensemble size and composition for estimates of adaptation effects and agreement of recommendations. Agricultural and Forest Meteorology, 264, 351–362.
Abstract: unless local adaptation can ameliorate these impacts. Ensembles of crop simulation models can be useful tools for assessing if proposed adaptation options are capable of achieving target yields, whilst also quantifying the share of uncertainty in the simulated crop impact resulting from the crop models themselves. Although some studies have analysed the influence of ensemble size on model outcomes, the effect of ensemble composition has not yet been properly appraised. Moreover, results and derived recommendations typically rely on averaged ensemble simulation results without accounting sufficiently for the spread of model outcomes. Therefore, we developed an Ensemble Outcome Agreement (EOA) index, which analyses the effect of changes in composition and size of a multi-model ensemble (MME) to evaluate the level of agreement between MME outcomes with respect to a given hypothesis (e.g. that adaptation measures result in positive crop responses). We analysed the recommendations of a previous study performed with an ensemble of 17 crop models and testing 54 adaptation options for rainfed winter wheat (Triticum aestivwn L.) at Lleida (NE Spain) under perturbed conditions of temperature, precipitation and atmospheric CO2 concentration. Our results confirmed that most adaptations recommended in the previous study have a positive effect. However, we also showed that some options did not remain recommendable in specific conditions if different ensembles were considered. Using EOA, we were able to identify the adaptation options for which there is high confidence in their effectiveness at enhancing yields, even under severe climate perturbations. These include substituting spring wheat for winter wheat combined with earlier sowing dates and standard or longer duration cultivars, or introducing supplementary irrigation, the latter increasing EOA values in all cases. There is low confidence in recovering yields to baseline levels, although this target could be attained for some adaptation options under moderate climate perturbations. Recommendations derived from such robust results may provide crucial information for stakeholders seeking to implement adaptation measures.
|
|
