Robust Inverse Modeling of Growing Season Net Ecosystem Exchange in a Mountainous Peatland : Influence of Distributional Assumptions on Estimated Parameters and Total Carbon Fluxes

While boreal lowland bogs have been extensively studied using the eddy‐covariance (EC) technique, less knowledge exists on mountainous peatlands. Hence, half‐hourly CO2 fluxes of an ombrotrophic peat bog in the Harz Mountains, Germany, were measured with the EC technique during a growing season with exceptionally dry weather spells. A common biophysical process model for net ecosystem exchange was used to describe measured CO2 fluxes and to fill data gaps. Model parameters and uncertainties were estimated by robust inverse modelling in a Bayesian framework using a population‐based Markov Chain Monte Carlo sampler. The focus of this study was on the correct statistical description of error, i.e. the differences between the measured and simulated carbon fluxes, and the influence of distributional assumptions on parameter estimates, cumulative carbon fluxes, and uncertainties. We tested the Gaussian, Laplace, and Student's t distribution as error models. The t‐distribution was identified as best error model by the deviance information criterion. Its use led to markedly different parameter estimates, a reduction of parameter uncertainty by about 40%, and, most importantly, to a 5% higher estimated cumulative CO2 uptake as compared to the commonly assumed Gaussian error distribution. As open‐path measurement systems have larger measurement error at high humidity, the standard deviation of the error was modeled as a function of measured vapor pressure deficit. Overall, this paper demonstrates the importance of critically assessing the influence of distributional assumptions on estimated model parameters and cumulative carbon fluxes between the land surface and the atmosphere.