Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Geoscientific Model Development 11 (2018): 497-519, doi:10.5194/gmd-11-497-2018.

Description: 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 =  0.76; Nash–Sutcliffe modeling efficiency, MEF  =  0.76) and ecosystem respiration (ER, r2 =  0.78, MEF  =  0.75), with lesser accuracy for latent heat fluxes (LE, r2 =  0.42, MEF  =  0.14) and and net ecosystem CO2 exchange (NEE, r2 =  0.38, MEF  =  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 〈 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized Vcmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average Vcmax value.

Description: This study was supported by the European Research Council Synergy grant ERC-2013-SyG- 610028 IMBALANCE-P.

Description: Wetlands are one of the most significant natural sources of methane (CH4) to the atmosphere. They emit CH4 because decomposition of soil organic matter in waterlogged anoxic conditions produces CH4, in addition to carbon dioxide (CO2). Production of CH4 and how much of it escapes to the atmosphere depend on a multitude of environmental drivers. Models simulating the processes leading to CH4 emissions are thus needed for upscaling observations to estimate present CH4 emissions and for producing scenarios of future atmospheric CH4 concentrations. Aiming at a CH4 model that can be added to models describing peatland carbon cycling, we composed a model called HIMMELI that describes CH4 build-up in and emissions from peatland soils. It is not a full peatland carbon cycle model but it requires the rate of anoxic soil respiration as input. Driven by soil temperature, leaf area index (LAI) of aerenchymatous peatland vegetation, and water table depth (WTD), it simulates the concentrations and transport of CH4, CO2, and oxygen (O2) in a layered one-dimensional peat column. Here, we present the HIMMELI model structure and results of tests on the model sensitivity to the input data and to the description of the peat column (peat depth and layer thickness), and demonstrate that HIMMELI outputs realistic fluxes by comparing modeled and measured fluxes at two peatland sites. As HIMMELI describes only the CH4-related processes, not the full carbon cycle, our analysis revealed mechanisms and dependencies that may remain hidden when testing CH4 models connected to complete peatland carbon models, which is usually the case. Our results indicated that (1) the model is flexible and robust and thus suitable for different environments; (2) the simulated CH4 emissions largely depend on the prescribed rate of anoxic respiration; (3) the sensitivity of the total CH4 emission to other input variables is mainly mediated via the concentrations of dissolved gases, in particular, the O2 concentrations that affect the CH4 production and oxidation rates; (4) with given input respiration, the peat column description does not significantly affect the simulated CH4 emissions in this model version.