Search

Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r2 Combining double low line 0.76; Nash-Sutcliffe modeling efficiency, MEF Combining double low line 0.76) and ecosystem respiration (ER, r2 Combining double low line 0.78, MEF Combining double low line 0.75), with lesser accuracy for latent heat fluxes (LE, r2 Combining double low line 0.42, MEF Combining double low line 0.14) and and net ecosystem CO2 exchange (NEE, r2 Combining double low line 0.38, MEF Combining double low line 0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r2 values (0.57-0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r2 values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r2

The processes responsible for methane (CH4) emissions from boreal wetlands are complex; hence, their model representation is complicated by a large number of parameters and parameter uncertainties. The arctic-enabled dynamic global vegetation model LPJ-GUESS (Lund–Potsdam–Jena General Ecosystem Simulator) 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 the 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 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 parameters 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 for a greater search of their posterior distribution, leading to a more complete characterisation of the posterior distribution with a reduced risk of the sample impoverishment that can occur when using other optimisation methods. For assimilation, the analysis used flux measurement data gathered during the period from 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. The experiment using real observations from Siikaneva resulted in a reduction in the root-mean-square error (RMSE), from 0.044 to 0.023 gC m -2 d -1 , and a 93.89 % reduction in the cost function value. As a part of this work, knowledge about how CH4 data can constrain the parameters and processes is derived. Although 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 a more robust calibration of LPJ-GUESS and similar models, and the results can highlight where model improvements are needed. [ABSTRACT FROM AUTHOR]