<p><span>Wetlands in the boreal zone are a significant source of atmospheric methane, and hence they have been intensively studied with mechanistic models for the assessment of methane dynamics. The arctic-enabled dynamic global vegetation model LPJ-GUESS is one of the models that allow quantification and understanding of the natural methane fluxes at various scales ranging from local to regional and global, but with several uncertainties. Complexity in the underlying environmental processes, warming driven alternative paths of meteorological phenomena and changes in hydrological and vegetation conditions are exigent for a calibrated and optimised LPJ-GUESS. In this study, we used the Markov chain Monte Carlo (using Metropolis-Hastings formula) algorithm to quantify the uncertainties of LPJ-GUESS. Application of this method allows greater search of the posterior distribution, leading to a more complete characterisation of the posterior distribution with reduced risk of sample impoverishment. We will present first results from an assimilation experiment optimising LPJ-GUESS model process parameters using the flux measurement data from 2005 to 2015 from the Siikaneva wetlands in southern Finland. We<span>  </span>analyse the parameter efficiency of LPJ-GUESS by looking into the posterior parameter distributions, parameter correlations, and the interconnections of the processes they control. As a part of this work, knowledge about how the methane data can constrain the parameters and processes is derived. </span></p>

Abstract. Processes behind methane (CH4) emission from boreal wetlands are complex, and hence their model representation is complicated by a large number of parameters and parameter uncertainties. The arctic-enabled dynamic global vegetation model LPJ-GUESS is one such model that allows quantification and understanding of the natural wetland CH4 fluxes at various scales ranging from local to regional and global, but with several uncertainties. The model contains detailed descriptions of CH4 production, oxidation, and transport controlled by several process parameters. Complexities in the underlying environmental processes, warming-driven alternative paths of meteorological phenomena, and changes in hydrological and vegetation conditions are highlight the need for a calibrated and optimised version of LPJ- GUESS. In this study we formulated the parameter calibration as a Bayesian problem, using knowledge of reasonable pa- rameters values as priors. We then used an adaptive Metropolis Hastings (MH) based Markov Chain Monte Carlo (MCMC) algorithm to improve predictions of CH4 emission by LPJ-GUESS and to quantify uncertainties. Application of this method on uncertain parameters allows greater search of their posterior distribution, leading to a more complete characterisation of the posterior distribution with reduced risk of sample impoverishment that can occur when using other optimisation methods. For assimilation, the analysis used flux measurement data gathered during the period 2005 to 2014 from the Siikaneva wetlands in southern Finland with an estimation of measurement uncertainties. The data are used to constrain the processes behind the CH4 dynamics, and the posterior covariance structures are used to explain how the parameters and the processes are related. To further support the conclusions, the CH4 flux and the other component fluxes associated with the flux are examined. The results demonstrate the robustness of MCMC methods to quantitatively assess the interrelationship between objective function choices, parameter identifiability, and data support. As a part of this work, knowledge about how the CH4 data can constrain the parameters and processes is derived. Though the optimisation is performed based on a single site’s flux data from Siikaneva, the algorithm is useful for larger-scale multi-site studies for more robust calibration of LPJ-GUESS and similar models, and the results can highlight where model improvements are needed.

The Adaptive Metropolis algorithm used here contains three key concepts explained in the following sections. The adaptive random walk, allowing the algorithm to learn features of the target distribution; transformed proposals, providing a natural way of including limits on the parameters; and tempering of the target distribution, to reduce the effects of local maxima and allowing better exploration of the target.