Abstract:Over the past years, Bayesian calibration methods have been successfully applied to calibrate ecosystem models. Bayesian methods combine prior probability distributions of model parameters, based on assumptions about their magnitude and uncertainty, with estimates of the likelihood of the simulation results by comparison with observed values. Bayesian methods also quantify the uncertainty in the updated posterior parameters, which can be used to perform an analysis of model output uncertainty. In this paper, we applied Bayesian techniques to calibrate the VSD soil acidification model using data from 182 intensively monitored forest sites in Europe. Out of these 182 plots, 122 plots were used to calibrate VSD and the remaining 60 plots to validate the calibrated model. Prior distributions for the model parameters were based on available literature. Since the available literature shows a strong dependence of some VSD parameters on, for example, soil texture, prior distributions were allowed to depend on soil group (i.e. soils with similar texture or C/N ratio). The likelihood was computed by comparing modelled soil solution concentrations with observed concentrations for the period 1996–2001. Markov Chain Monte Carlo (MCMC) was used to sample the posterior parameter space. Two calibration approaches were applied. In the single-site calibration, the plots were calibrated separately to obtain plot-specific posterior distributions. In the multi-site approach priors were assumed constant in space for each soil group, and all plots were calibrated simultaneously yielding one posterior probability distribution for each soil group. Results from the single-site calibrations show that the model performed much better after calibration compared to a run with standard input parameters when validated on the 60 validation plots. Posterior distributions for H-Al equilibrium constants narrowed down, thus decreasing parameter uncertainty. For base cation weathering of coarse textured soils the posterior distribution shifted to larger values, indicating an initial underestimation of the weathering rate for these soils. Results for the parameters related to nitrogen modelling showed that the nitrogen processes model formulations in VSD may have to be reconsidered as the relationship between nitrogen immobilization and the C/N ratio of the soil, as assumed in VSD, was not substantiated by the validation. The multi-site calibration also strongly decreased model error for most model output parameters, but model error was somewhat larger than the median model error from the single-site calibration except for nitrate. Because the large number of plots calibrated at the same time provided very many observations, the Markov Chain converged to a very narrow parameter space, leaving little room for posterior parameter uncertainty. For an uncertainty analysis with VSD on the European scale, this study provides promising results, but more work is needed to investigate how the results can be used on a European scale by looking at regional patterns in calibrated parameters from the site calibration or by calibrating for regions instead of all of Europe.

A model to calculate effects of atmospheric deposition on soil acidification, eutrophication and carbon sequestration

A very simple dynamic soil acidification model for scenario analyses and target load calculations

Modelling the Impact of Acid Deposition on Forest Soils in North Bohemian Mountains with Two Dynamic Models: the Very Simple Dynamic Model (VSD) and the Model of Acidification of Groundwater in Catchments (MAGIC)

Feasibility of coupled empirical and dynamic modeling to assess climate change and air pollution impacts on temperate forest vegetation of the eastern United States

Sensitivity of soil acidification model to deposition and forest growth

Can models adequately reflect how long-term nitrogen enrichment alters the forest soil carbon cycle?

MAGIC, SAFE and SMART Model Applications at Integrated Monitoring Sites: Effects of Emission Reduction Scenarios

Comparison of measured and simulated changes in base cation amounts using a one-layer and a multi-layer soil acidification model

Atmospheric deposition, CO 2 , and change in the land carbon sink

Modelling carbon sources and sinks in terrestrial vegetation

Ecosystem C and N cycle interactions – diverse model representations and divergent model predictions versus collective empirical constraints

Bayesian calibration of the VSD soil acidification model using European forest monitoring data

Atmospheric deposition, CO2, and change in the land carbon sink

Modelling and Mapping of Critical Loads of Acidification and Eutrophication on Forest Ecosystems in Latvia

Simulating Global Terrestrial Carbon and Nitrogen Biogeochemical Cycles With Implicit and Explicit Representations of Soil Microbial Activity

The contribution of nitrogen deposition to the eutrophication signal in understorey plant communities of European forests

Nitrogen and Sulfur Deposition Reductions Projected to Partially Restore Forest Soil Conditions in the US Northeast, While Understory Composition Continues to Shift with Future Climate Change

Simulating ectomycorrhizal fungi and their role in carbon and nitrogen cycling in forest ecosystems

Does forest tree species composition impact modelled soil recovery from acidic deposition?

Representing nitrogen, phosphorus, and carbon interactions in the E3SM Land Model: Development and global benchmarking

Where does the carbon go? A model-data intercomparison of vegetation carbon allocation and turnover processes at two temperate forest free-air CO2 enrichment sites