Articles | Volume 15, issue 7
Geosci. Model Dev., 15, 2813–2838, 2022
Geosci. Model Dev., 15, 2813–2838, 2022
Model evaluation paper
06 Apr 2022
Model evaluation paper | 06 Apr 2022

Assessing methane emissions for northern peatlands in ORCHIDEE-PEAT revision 7020

Assessing methane emissions for northern peatlands in ORCHIDEE-PEAT revision 7020
Elodie Salmon1, Fabrice Jégou1, Bertrand Guenet2,3, Line Jourdain1, Chunjing Qiu2, Vladislav Bastrikov4, Christophe Guimbaud1, Dan Zhu2,5, Philippe Ciais2, Philippe Peylin2, Sébastien Gogo6, Fatima Laggoun-Défarge6, Mika Aurela7, M. Syndonia Bret-Harte8, Jiquan Chen9, Bogdan H. Chojnicki10, Housen Chu11, Colin W. Edgar8, Eugenie S. Euskirchen8, Lawrence B. Flanagan12, Krzysztof Fortuniak13, David Holl14, Janina Klatt15, Olaf Kolle16, Natalia Kowalska17, Lars Kutzbach14, Annalea Lohila7, Lutz Merbold18, Włodzimierz Pawlak13, Torsten Sachs19, and Klaudia Ziemblińska20 Elodie Salmon et al.
  • 1Laboratoire de Physique et Chimie de l'Environnement et de l'Espace, LPC2E, UMR 7328, Université d'Orléans, CNRS, CNES, 45071 Orléans CEDEX 2, France
  • 2Laboratoire des Sciences du Climat et de l'Environnement, UMR8212, CEA-CNRS-UVSQ, 91191 Gif sur Yvette, France
  • 3Laboratoire de Géologie de l'ENS, IPSL, CNRS, PSL Research University, Laboratoire de Géologie de l'ENS, 24 rue Lhomond, 75231 Paris CEDEX 05, France
  • 4Science Partners, 75010 Paris, France
  • 5Sino-French Institute for Earth System Science, College of Urban and Environmental Sciences, Peking University, Beijing, China
  • 6Institut des Sciences de la Terre d'Orléans, Université d'Orléans, CNRS, BRGM, UMR 7327, 45071 Orléans, France
  • 7Finnish Meteorological Institute, Climate Research Programme, Helsinki, Finland
  • 8Institute of Arctic Biology, University of Alaska Fairbanks, Fairbanks, AK, USA
  • 9Landscape Ecology & Ecosystem Science (LEES) Lab, Department of Geography, Environment, and Spatial Sciences, & Center for Global Change and Earth Observations, Michigan State University, East Lansing, MI 48823, USA
  • 10Laboratory of Bioclimatology, Department of Ecology and Environmental Protection, Faculty of Environmental Engineering and Mechanical Engineering, Poznan University of Life Sciences, Piątkowska 94, 60-649 Poznań, Poland
  • 11Climate and Ecosystem Sciences Division, Lawrence Berkeley National Lab, USA, 1 Cyclotron Rd, Berkeley, CA 94720, USA
  • 12Department of Biological Sciences, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta, Canada
  • 13Department of Meteorology and Climatology, Faculty of Geographical Sciences, University of Lodz, Lodz, Poland
  • 14Institute of Soil Science, Center for Earth System Research and Sustainability (CEN), Universität Hamburg, Hamburg, Germany
  • 15Institute of Ecology and Landscape, Chair of Vegetation Ecology, University of Applied Sciences Weihenstephan-Triesdorf, Am Hofgarten 1, 85354 Freising, Germany
  • 16Field Experiments and Instrumentation, Max Planck Institute for Biogeochemistry, Hans-Knoell-Strasse 10, 07745 Jena, Germany
  • 17Department of Matter and Energy Fluxes, Global Change Research Institute, Czech Academy of Sciences, Bělidla 986/4a, 603 00 Brno, Czech Republic
  • 18Department Agroecology and Environment, Agroscope, Reckenholzstrasse 191, 8046 Zurich, Switzerland
  • 19GFZ German Research Centre for Geosciences, Telegrafenberg, Potsdam, Germany
  • 20Laboratory of Meteorology, Department of Construction and Geoengineering, Faculty of Environmental Engineering and Mechanical Engineering, Poznan University of Life Sciences, Piątkowska 94, 60-649 Poznań, Poland

Correspondence: Elodie Salmon ( and Fabrice Jégou (


In the global methane budget, the largest natural source is attributed to wetlands, which encompass all ecosystems composed of waterlogged or inundated ground, capable of methane production. Among them, northern peatlands that store large amounts of soil organic carbon have been functioning, since the end of the last glaciation period, as long-term sources of methane (CH4) and are one of the most significant methane sources among wetlands. To reduce uncertainty of quantifying methane flux in the global methane budget, it is of significance to understand the underlying processes for methane production and fluxes in northern peatlands. A methane model that features methane production and transport by plants, ebullition process and diffusion in soil, oxidation to CO2, and CH4 fluxes to the atmosphere has been embedded in the ORCHIDEE-PEAT land surface model that includes an explicit representation of northern peatlands. ORCHIDEE-PCH4 was calibrated and evaluated on 14 peatland sites distributed on both the Eurasian and American continents in the northern boreal and temperate regions. Data assimilation approaches were employed to optimized parameters at each site and at all sites simultaneously. Results show that methanogenesis is sensitive to temperature and substrate availability over the top 75 cm of soil depth. Methane emissions estimated using single site optimization (SSO) of model parameters are underestimated by 9 g CH4 m−2 yr−1 on average (i.e., 50 % higher than the site average of yearly methane emissions). While using the multi-site optimization (MSO), methane emissions are overestimated by 5 g CH4 m−2 yr−1 on average across all investigated sites (i.e., 37 % lower than the site average of yearly methane emissions).

1 Introduction

The atmospheric methane level estimated from ice cores analysis (Etheridge et al., 1998) and in situ measurements (Blake et al., 1982; Dlugokencky, 2021; Prinn et al., 2018) has nearly tripled since the preindustrial equilibrium value, i.e., from 680 ppb to reach a value of 1892 ppb in December 2020 (Dlugokencky, 2021; Saunois et al., 2020). This increase is consistent with the world population increase and industrialization, such as the increase in fossil fuel extraction and use, organic waste generation, and livestock numbers (Raynaud et al., 2003).

Methane is the second most important anthropogenic greenhouse gas (GHG) after CO2 and accounts for about 23 % of the cumulative total radiative forcing (Etminan et al., 2016). In the troposphere methane is an ozone precursor, and in the stratosphere, methane interacts with hydroxyl radicals and carbon monoxide to produce water vapor. About 90 % of CH4 is oxidized by the hydroxyl radical in the troposphere (Smith et al., 2003) and reactions with chlorine in the stratosphere or in the marine boundary layer (Allan et al., 2007; Thornton et al., 2010), leading to a residence time of about 9 years (Prather et al., 2012). At the continental surface, 5 % to 10 % of all methane sources is removed from the atmosphere by diffusion in soils and oxidation by soil microorganisms (Krüger et al., 2002; Prather et al., 1995; Smith et al., 2003, 1991; Tokida et al., 2007a, b). Among natural sources, natural wetlands are the largest contributor and the most uncertain one in the global budget (Kirschke et al., 2013; Saunois et al., 2016). They contribute 25 %–30 % of total methane emissions estimated by Saunois et al. (2020) and encompass anaerobic ecosystems composed of waterlogged or inundated ground that are capable of methane production, which include peatlands, mineral soil wetlands, and floodplains. Peatlands are of particular interests because peat is composed of organic detritus and has an average carbon content of 52 % dry mass (Gorham, 1991). Consequently, peatlands are large soil organic carbon reservoirs that could be functioning as a source of CH4 and a source or sink of CO2 to the atmosphere. They cover around 3 % of surface continental lands but store approximately one-third of the global soil carbon (Gorham, 1991). They are located in boreal and sub-arctic regions (80 %, Strack et al., 2008), although some smaller areas are found in temperate and tropical regions (10 %–12 %). Since the end of the last glaciation period (around 16 500 years ago), northern peatlands have been functioning as long-term carbon sinks. This storage results from a delicate balance between carbon inputs (CO2 absorbed by photosynthesis) and carbon outputs (CO2 and CH4 production, dissolved and particulate carbon). Clearly, in these ecosystems, processes controlling methane production, fluxes between the land surface and the atmosphere, and feedback on climate are intimately connected.

The major pathway for methane production is via microbial processes, which is limited by the availability of substrates (polymeric and monomeric compounds derived from carbohydrates, fatty acids, amino acids, acetate, and hydrogen; (Blodau, 2002; Le Mer and Roger, 2001), the low oxygen content that is directly correlated with soil water content, and soil temperature. After its production, CH4 migrates to the soil surface and is emitted to the atmosphere through three main processes (Bridgham et al., 2013): (1) diffusion through porous soil media; (2) ebullition, whereby bubbles form in pores filled with water then quickly migrate to the surface; and (3) plant-mediated fluxes via some vascular plant adapted to live in flooded environments. These plants have developed aerenchyma to channel gas fluxes; oxygen is transported to roots and cells, and CH4 is transported from roots to the atmosphere (Bridgham et al., 2013; Smith et al., 2003).

Since the late 1980s, many CH4 cycling processes have been mathematically described and included in terrestrial ecosystem models (Xu et al., 2016). These terrestrial ecosystem models have been outlined in two broad categories by the Xu et al. (2016) review: (1) empirical models employed to evaluate observed processes of the CH4 cycling and (2) process-based models used for budget quantification and to study sensitivity of CH4 processes to environmental drivers. Unfortunately, so far only a few global-scale models have featured peatland ecosystems, permafrost dynamics, and CH4 fluxes, which are essential features to evaluate future climate changes and interactions between the land surface and the atmosphere (Anav et al., 2013). Recent developments of the ORCHIDEE land surface model have led to simulations of soil hydrology, permafrost thermodynamics, and the carbon cycle at northern latitudes (Guimberteau et al., 2018) and in northern peatlands specifically (Qiu et al., 2018), including peat carbon decomposition controlled by soil water content and temperature as well as CO2 production and consumption processes (Largeron et al., 2018; Qiu et al., 2018). In the present study we adapt the Khvorostyanov et al. (2008a, b) methane model to ORCHIDEE-PEAT (Sect. 2.1) and calibrate and evaluate simulated emissions at northern peatland sites. To achieve model calibration, parameters were optimized with a data assimilation approach described in Sect. 2.3. Parameters were optimized against methane fluxes at each site and from multiple sites simultaneously (Sect. 3) in order to highlight parameter uncertainties while scaling up simulations from site scale to larger scale. The model evaluation is performed by discussing both optimization methods.

2 Model description

A general presentation of ORCHIDEE-PCH4 and associated processes is provided in Sect. 2.1. Implementations of methane production and oxidation as well as transport are respectively specified in Sect. 2.1.1 and 2.1.2, whereas parameter values established for the site simulation conditions before observation periods are given in Sect. 2.2. Section 2.3 describes the parameter optimization approaches.


The ORCHIDEE land surface model is a dynamic global process-oriented model that simulates carbon, water, and energy fluxes between the biosphere, land surface geosphere, and atmosphere. The carbon scheme describes photosynthesis, respiration, soil carbon cycle, and CO2 production and emissions. One of the branches of the ORCHIDEE land model aimed to improve the implementation of high-latitude physical, hydrological, and biogeochemical processes such as soil thermal processes, hydraulic processes, snowpack properties, and plant and soil carbon fluxes (ORCHIDEE-MICT, Guimberteau et al., 2018).

A northern peatlands scheme has been recently integrated into the model (ORCHIDEE-PEAT, Largeron et al., 2018; Qiu et al., 2018), which includes a peatland PFT (plant functional type) with adapted biological parameters created to allow a separate calculation of the water balance. This PFT is defined as a flood-tolerant C3 grass with reduced productivity due to the lack of nutrients and with a reduced rooting depth. For the present study, ORCHIDEE-PEAT v2.0 (Qiu et al., 2019) has been further enriched with a module simulating methane production, oxidation, and transport in northern peatlands; it is named ORCHIDEE-PCH4. To achieve this, the methane scheme described by Khvorostyanov et al. (2008a, b) was revised according to high-latitude processes and peatland ecosystem features. This early version was an idealized 1D soil model that accounted for heat and gas transport as well as soil organic carbon decomposition and production of CO2 and CH4 driven by soil water content and temperature in the soil column. In that early version, only a moss layer that serves as a thermal insulator was considered for the vegetation above ground (Khvorostyanov et al., 2008a). Soil moisture and carbon dynamics were treated as a single-layer bucket scheme of 1 m depth containing a fixed amount of soil carbon content. In contrast, ORCHIDEE-PCH4 is integrated into the peatland soil hydrological diffusion model (Largeron et al., 2018; Qiu et al., 2018) that incorporates water supply by precipitation and runoff collected from other soils surrounding the peatland in the same grid cell. The deep drainage is blocked to maintain soil water content at saturation in the bottom part of the peat soil. At the top of the water column, a dynamic water reservoir was added to represent standing water above the soil surface when water inputs exceed outputs and when soil is fully saturated. ORCHIDEE-PEAT simulates peat accumulation and decomposition to CO2 of the three soil carbon pools (active, slow, and passive) that are vertically discretized in 32 layers, accounting for a total maximum depth of 38 m (Qiu et al., 2019).

The methane scheme in Fig. 1 delineates (1) methanogenesis of the three carbon pools, (2) methane and oxygen transport in the soil and snow layers, (3) transport of methane to the atmosphere by ebullition, (4) plant-mediated transport, and (5) methanotrophy by soil oxic conditions and root exudates.

Figure 1Model diagram of methane cycling processes in ORCHIDEE-PCH4. Carbon fluxes are indicated by thin black arrows. Other variables that influence each carbon flux are displayed on white arrows.


Each of these processes is constrained by soil temperature, soil water content (θsoil), soil O2 concentration, atmospheric CH4 concentration, leaf area, and snow cover. The temporal variation of CH4 in the soil layer (z) is assessed by

(1) [ CH 4 ] ( z , t ) t = f MG a + f MG s + f MG p - f Diff - f Ebu - f PMT - f MT ,

where each term varies in time (t) and with depth (z). The equation expresses methane production (fMG, MG: methanogenesis, a: active pool, s: slow pool, p: passive pool), transport by diffusion, ebullition, and plant (fDiff, fEbu, fPMT) and oxidation (fMT, MT: methanotrophy) processes. Net methane fluxes to the atmosphere are the sum of methane transport processes fEbu (Ebu: ebullition) and fPMT (PMT: plant-mediated transport) as well as the amount of CH4 that diffuses from the topsoil layer at z=0 to the atmosphere. Prognostic variables are defined per air volume, i.e., the volume of gas in the air-filled pores (ν) and gas dissolved in the water-filled pores (Khvorostyanov et al., 2008a; Tans, 1998; Tang et al., 2013; Tang and Riley, 2014), assuming a constant equilibrium between gas concentrations in the air-filled and the water-filled part of pores. This gas volume is linked to the soil volume by the total CH4 and O2 in pores (εgas, gas = O2, CH4) defined as

(2) ε gas = ν + θ soil π soil B gas ,

where θsoil is the volumetric water content of the soil, πsoil is the soil porosity, and Bgas is the Bunsen gas solubility coefficient defined for CH4 and O2, respectively, with BCH4=0.043 and BO2=0.038 (Hodgman, 1936; Wiesenburg and Guinasso, 1979).

2.1.1 Methane production and oxidation

Methanogenesis in soil occurs when oxygen concentration is limited for microorganisms and is considered for each type of soil carbon pool ([C]i, i=a,s,p; in g C m−3 of soil), active, slow, and passive:

(3) f MG i = [ C ] i k i q MG e - O 2 p / O 2 anoxia f clay ,

where the rate of methanogenesis (ki in s−1) depends on soil temperature and moisture according to the same function as for the heterotrophic respiration (Qiu et al., 2019). This rate has been defined by Khvorostyanov et al. (2008a) to be 10 times lower than the rate of heterotrophic respiration. Here, qMG determines the ratio between the rate of soil oxic and anoxic decomposition, [O2]p is the oxygen concentration in the soil [O2]soil (in g O2 m−3 of soil) per unit porous volume (εO2πsoil , πsoil is the soil porosity), and [O2]anoxia is the soil oxygen concentration at which anoxic conditions are reached and enable methane production. This oxygen concentration threshold is assumed to be 2 g m−3 (Duval and Goodwin, 2000). Soil clay content affects the decomposition of the active soil carbon pool (Parton et al., 1988):

(4) f clay = 1 - 0.75 clay ,

where clay is the clay fraction and has a value of 0.2, the neither the slow nor the passive pools are modified by fclay. Methane is oxidized to CO2 in aerated soil layers. The amount of methane consumed by methanotrophy is limited by the soil oxygen concentration, [O2]soil, following a 1:2 CH4:O2 molar ratio:

(5) f MT = k MT 1 2 O 2 soil Mw CH 4 Mw O 2 ε O 2 ε CH 4 ,

where kMT is the rate of methanotrophy, the value of which ranges from 0.06 to 5 d−1 (Morel et al., 2019). The conversion of oxygen to methane content is provided by methane and oxygen molecular weights MwCH4 and MwO2 and their respective total gas porosities εCH4 and εO2.

2.1.2 Methane transport

The formation of methane bubbles in water-filled pores is determined by

(6) f Ebu = k Ebu CH 4 soil - CH 4 ET p Ebu ,

where kEbu is a rate constant of 1 h−1. Methane ebullition occurs when methane concentration exceeds a concentration threshold that depends on soil temperature (Tsoil) and pressure (Psoil in Pa). Above 0.75 m depth it is calculated as follows:

(7) CH 4 ET = mxr CH 4 P soil Mw CH 4 R T soil B CH 4 ,

where mxrCH4 is the methane mixing ratio in the bubbles. Walter and Heimann (2000) determined this mixing ratio to range between 27 % and 53 % for totally vegetated and unvegetated soil, and Riley et al. (2011) calculated it at 15 %. It is converted to g CH4 per unit porous volume by an ideal gas constant (R), MwCH4, and the Bunsen methane solubility coefficients (BCH4). It has been suggested that ebullition in soil occurs when the partial pressure of dissolved gases exceeds the hydrostatic pressure (Chanton and Whiting, 1995). We estimated that in our model below the layer corresponding to 0.75 m the hydrostatic pressure is always higher than the partial pressure of dissolved gases. Therefore, we considered the methane ebullition threshold to be constant below 0.75 m and equal to the value defined at 0.75 m in order to avoid methane accumulation in the deeper layers. The methane flux provided by ebullition (fEbu) is modulated by the probability of methane bubbles reaching the soil surface. Indeed, in the soil column the water table level fluctuates, modifying the connectivity between water-filled pores involving variation of the surface methane flux. Therefore, the probability that methane bubbles will escape to the atmosphere is expressed as

(8) p Ebu = θ soil z Δ z / ( wsize × η ) ,

where θsoil(z) is the soil water content, Δz is the soil layer thickness, and the tortuosity η that depicts the sinuous path of bubbles is defined to be 2/3 (Hillel, 1982). The term wsize sizes the extent of the connected network of water-filled pores envisioned that can be depicted as droplets dispersed in the pores. Khvorostyanov et al. (2008a, b) defined wsize = 1 cm for a carbon-rich loess deposit of the Yedoma.

In wetlands, some vascular plants have developed a strategy to carry oxygen down to their root tips by employing aerenchyma tissue. These tissues are air channels in which gas exchange depends on the gradient of gas concentrations between the soil and the atmosphere. Oxygen is transferred from the atmosphere to the roots and creates an aerobic zone around them in which methane will be oxidized. The proportion of methane oxidized (Mrox) in the root zone is emitted as CO2 to the atmosphere. Walter and Heimann (2000) estimated Mrox to range between 39 % and 98 % of methane located in the root zone. Conversely, the methane concentration gradient results in a flux to the atmosphere through plants that is expressed by

(9) f PMT = k PMT × T veg × f root × LAI × CH 4 soil - CH 4 atm × ( 1 - M rox ) ,

where kPMT is a rate constant of the unit 0.01 h−1, and Tveg has been defined by Walter and Heimann (2000) as a factor that describes the efficiency of plants in methane transport depending on the type and density of these plants. Its value ranges between 0 and 15, with shrubs and trees being poorly efficient and grasses and sedges being very efficient in gas transport. The methane concentration gradient is also modified by the vertical distribution of roots in the soil as

(10) f root = 2 × z root - z soil z root .

This function describes the vertical distribution of roots in the soil in which zroot is the rooting depth and zsoil the soil depth. The leaf area index (LAI) influences the methane flux, which varies by growing stage of the plants.

The gas diffusion scheme features the diffusion of CH4 and O2 in the three top layers of snow when snow cover is formed and in the 32 soil layers that correspond to 38 m depth. This scheme considered (1) the diffusion of oxygen from the topsoil to the soil layer, (2) the diffusion of methane produced and remaining in the soil, and (3) methane exchange between the soil and the atmosphere at z=0:

(11) f Diff = D gas ( z ) gas soil t , z z .

Diffusion coefficients, Dgas, are based on the diffusivity of each gas in air (Dgas, air) and in water (Dgas, water):

(12) D gas = D gas , air ν + D gas , water θ soil π soil B gas η ,

where ν is the volume of gas in the air-filled pores, θsoil is the volumetric water content of soil, πsoil is the soil porosity, and Bgas is the Bunsen coefficient of the gas; the tortuosity η is defined to be 2/3 (Hillel, 1982). Diffusivities of O2 in air and in water are respectively defined to 1.6 × 10−5 and 1.6 × 10−9 m2 s−1 and for methane 1.7 × 10−5 and 2.0  × 10−9 m2 s−1 (Khvorostyanov et al., 2008a). The diffusion is discretized using a forward time-centered space method (Press et al., 1993) and converted into a tridiagonal system of equations before being solved using a forward then backward substitution method. A time-splitting option is also implemented for the diffusion of large concentrations of gas per time step.

The only source of oxygen considered is from the atmosphere and is determined using atmospheric surface pressure, temperature, and an atmospheric O2 mixing ratio of 20.9 %. Atmospheric methane content is also defined in the same way by employing a methane mixing ratio of 1.7 ppm and is used as a boundary condition when the topsoil layer is in contact with the atmosphere. In winter, when snow accumulates above the topsoil, these atmospheric boundary conditions are applied to the top snow layer, and then gases diffuse from and to the atmosphere through the snow layers, then soil layers. Methane and oxygen diffusivity in the snow are defined by

(13) D gas = D gas , air 1 - ρ snow ρ ice η snow ,

where Dgas, air the diffusion coefficient of each gas in free air, the snow porosity is defined by the ratio of density of snow ρsnow and ice ρice, and the tortuosity (ηsnow) is equal to 31-ρsnowρice. Snow density is determined by the snowpack scheme (Wang et al., 2013), with the density of the ice being 920.0 kg m−3.

2.2 Site description and simulation setup

The model was evaluated on 14 peatland sites distributed on the Eurasian and American continents in boreal and temperate northern regions (from 41 to 69 N). These sites are a subset of the 30 peatland sites collected for the calibration of ORCHIDEE-PEAT (Qiu et al., 2018), for which, in addition to eddy covariance data and physical variables (water table, snow depth, soil temperature), methane emissions were measured by eddy covariance at a daily timescale at US-Los, hourly timescale at DK-Nuf, and otherwise at a half-hourly timescale or chamber measurements at a monthly timescale for FR-Lag and RU-Che. All methane emissions data were monthly averages. At DE-Sfn, DE-Hmm, FI-Lom, PL-Kpt, PL-Wet, and US-Wpt, year-round data were available, and zero values were filled for the first and the last month of years at the beginning and the end of the observation period. Otherwise, winter months were filled with zero, and during spring, summer, and fall months missing data were gap-filled using a linear regression. Descriptions of the sites were provided in Qiu et al. (2018). In Table 1, sites are assembled by increasing extreme values of mean monthly measurements of methane emission, then by locations and ecological characteristics. The extreme values of mean monthly measurements are the most reliable quantity of methane fluxes since periods of observation and monitoring frequency differ. Among the 14 peatlands, 9 sites are located in temperate regions, 3 in boreal regions, and 2 in arctic permafrost regions. The majority of the sites are fen (9 sites) and the others are three bogs (DE-Sfn, US-Bog, DE-Hmm), a marsh (US-Wpt), and a tundra (RU-Che). It is worth noticing that there is no obvious correlation between the magnitude of the monthly mean fluxes and types of ecosystems. Indeed, US-Los and DE-Spw are temperate fens that release less than 10 mg CH4 m−2 d−1. Sites emitting 10–150 mg m−2 d−1 are located in Germany, northwestern America, and France, among which half are fens and the other half are bogs. Half of them, including DE-Sfn, US-Bog, and CA-Wp1, are forested peatlands that release less than 55 mg CH4 m−2 d−1. The others, including DE-Zrk, DE-Hmm, and FR-Lag, experienced a temporary drainage event because of anthropogenic activities during years earlier than the observed period. Sites located in Finland, Denmark, and Poland are fens emitting between 150 and 400 mg m−2 d−1. The largest methane emitters are the arctic tundra RU-Che and the marsh US-Wpt, which released more than 500 mg m−2 d−1. All sites are covered with some snow during winter, and US-Bog and RU-Che are underlaid with permafrost located below 0.5 m.

Table 1Site characteristics. Site identification includes the country initials and a three-letter name for each site; locations of the sites are provided by the country, latitude (Lat), and longitude (Lon) values. Hydrological characteristics are distinguished by the type of ecosystem: fen, bog, tundra, and marsh. Y and N indicate the presence and absence of snow cover in winter, permafrost soil, and forest above the peat. Temporary drawdown of the water table level is specified by the presence and absence indicators Y or N.

Download Print Version | Download XLSX

Each peatland site is a sub-grid area embedded in the 0.5× 0.5 grid cells whose extent is determined by a fraction of grid area as defined in Table 2. These sub-grid areas enable the representation of ecosystem variability in which a specific scheme simulates soil hydrology, vegetation characteristics, and soil carbon cycling for northern peatlands. The fraction of peatlands per grid cell was defined by modifying the prescribed values employed by Qiu et al. (2018) in order to collect enough water to fill the peatland by runoff from the other soil fractions and elevate the water table level for northern peatlands. We employed vegetation phenotype properties and peatland fractions described in Qiu et al. (2019) as well as peatland hydrology and a carbon model as described in Qiu et al. (2019). Site simulations were then constrained at the grid cell scale with a half-hourly time series of meteorological conditions, e.g., air temperature, wind speed, wind direction, longwave incoming radiation, shortwave incoming radiation, specific humidity, atmospheric pressure, and precipitation. These time series are flux tower measurements that were gap-filled by the 6-hourly CRU-NCEP 0.5 global climate forcing dataset (Qiu et al., 2018). Other variables measured on a half-hourly time step at sites, e.g., CO2 and energy (latent heat: LE; sensible heat: H) fluxes, water table position, soil temperature, and snow depth, served for the calibration of peatland soil and vegetation phenotype characteristics such as the maximum rate of carboxylation (Vcmax). Optimized Vcmax values (Qiu et al., 2018) are utilized to capture spatial carbon flux gradients (gross primary production, ecosystem respiration, and net ecosystem exchange) at each peatland site. The peat model (Qiu et al., 2019) enables a vertical buildup of peat by simulating a downward movement of C when the discretized organic layers reach a threshold defined from a regression relationship between the carbon fraction and measured bulk density. This scheme in ORCHIDEE-PCH4 serves to constrain the vertical distribution of the soil carbon stock to the observed maximum peat depth. Simulations with ORCHIDEE-PCH4 driven by repeated site-specific meteorological conditions were performed for various periods of time to reach the observed soil carbon content and maximum peat depth (Table 2). During the first part of those simulations, atmospheric CO2 concentration was set to the preindustrial value at 285 ppm, and then from 1860 until the beginning of the respective observation period of methane emissions listed in Table 1, the CO2 concentration had risen. During soil carbon accumulation simulations, methane model parameters were defined as the default values defined in Table 3. Then during the site-specific measurement periods (Table 1), methane variables are calibrated against observed monthly average methane flux time series. A site-specific simulation over the observed period is run again using the optimized parameters.

Table 2Simulation conditions and framework to constrain peatland soil carbon stocks.

Download Print Version | Download XLSX

Table 3List of parameters driving the methane production, oxidation, and transport scheme in ORCHIDEE-PCH4.

Download Print Version | Download XLSX

2.3 Optimization of methane parameters

The methane scheme revisited in ORCHIDEE-PCH4 (described in Sect. 2.1) is driven by seven parameters (Table 3) that constrain methane production (qMG), oxidation (kMT, Mrox), and transport (mxrCH4, wsize, Tveg, zroot). In order to optimize these parameters, we employed the ORCHIDEE data assimilation system (Bastrikov et al., 2018) that relies on the minimization of a cost function employing a Bayesian statistical formalism that expresses the discrepancy between observations and simulated methane emissions as well as the difference between the optimized parameter values and the prior information on them, weighted by the uncertainties assigned to both observations and parameters. A random search algorithm based on the genetic algorithm (GA) serves to randomly iterate the set of seven parameters following the principles of genetics and natural selection similar to chromosome genetic sequencing (Goldberg, 1989; Haupt and Haupt, 2004). At each iteration, eight sets of parameters are defined from the previous iteration following crossover and mutation rules (Bastrikov et al., 2018). The frequency at which these rules are used is governed by the crossover-to-mutation ratio fixed to 4:1, the number of parameter blocks exchanged during crossover, which is 2, and the number of parameters perturbed during mutation, which is equal to 1. In addition, a ranking in ascending order of the corresponding cost function values of all sets of parameters serves to selectively preserve the set of parameters that reduces the gap between observations and simulation data.

Two types of simulations are performed over the site-specific observation period defined in Table 1: a single site (SS) experiment for which parameters are optimized for each site and a multi-site (MS) experiment that aims at refining one set of parameters considering all sites together. The single site experiments are performed for 100 iterations and aim at finding the lowest cost function employing the model–data root mean square difference (RMSD). Prior conditions for the single site experiment are described and listed in Table 3. Initial parameter values and ranges were derived from the literature and expert knowledge, and parameter uncertainties are defined as 40 % of the prescribed ranges. Across sites, mean values of each parameter serve as prior conditions for the multi-site experiment. The latter was performed for 50 iterations and aims to evaluate methane emission uncertainties at hemispheric scale when only one set of parameters is employed.

3 Results

3.1 Single site optimization (SSO)

For each site, to minimize the discrepancy between observed and simulated methane emissions, iterative single site simulations were performed. Successive runs serve to ensure that the minimum reached is not a local minimum. Results from the last minimization experience are reported in Table 4 (uncertainties in parameters at sites are in Table S1). As expected, most optimized parameters fit within the initial range defined in Table 3 except for four of the sites. One of these four sites, DE-Spw, is among the sites that emits the lowest amount of methane (up to 7 mg m−2 d−1) and features a larger stock of carbon of 84 kg C m−2 than at US-Los that features 27 kg C m−2 and emits up to 4 mg m−2 d−1. This explains, at the DE-Spw site, that the optimized value of wsize was reduced to 0.5 mm to maintain low methane emissions. The other three sites, for which some of the optimized parameters are out of the initial range (DK-Nuf, PL-Wet, and US-Wpt), are among the sites that emit more than 150 mg CH4 m−2 d−1. The carbon stocks at DK-Nuf and PL-Wet are respectively 55 and 38 kg C m−2, which is lower than at FI-Lom and PL-Kpt that accumulated more than 200 kg C m−2. Three parameter ranges were modified for DK-Nuf; the minimum value of qMG was lowered to 7.0, zroot maximum is increased to the maximum peat depth at 0.75 m in order to consider plant-mediated transport in all the peat layers, the maximum value of Tveg was increased to 40.0, and the maximum rate of methanotrophy kMT was enlarged up to 8 d−1 to decrease the methane oxidation and to obtain in the simulation methane emissions higher than 150 mg CH4 m−2 d−1. PL-Wet also required modifying range values of qMG to 1.0–11.0, leading to the lowest optimized qMG value of 4.0, which significantly reduced the RMSD from 227.4 to 80.5 (Fig. S1 and Table S2). For the US-Wpt site, qMG, kMT, and Tveg were adjusted to increase methane production and fluxes in order to balance the carbon stock of 5 kg C m−2, which is lower than the one at RU-Che.

Table 4Single site optimized values of methane scheme parameters for each peatland site. In parentheses are the prior parameter ranges which differ from the values in Table 3. Uncertainties for these ranges are specified in parentheses.

Download Print Version | Download XLSX

Table 5Discrepancies between observed and simulated methane emissions are quantified by the root mean square difference (RMSD) approach. Minimization efficiency of each test is indicated by the relationship between the prior using default values and posterior RMSD as (1  RMSDpost/ RMSDprior) × 100. Normalized root mean square difference (NRMSD) is defined by the RMSD posterior normalized by the annual mean of observed methane emissions.

Download Print Version | Download XLSX

Across sites, qMG values extend between 4.0 and 10.7, and optimized kMT values vary between 1 and 5.25 d−1. The fraction of methane that is oxidized at the root (Mrox) level fluctuates between 0.004 and 0.99, with the lowest values obtained at US-Wpt and RU-Che sites that emitted up to 500 mg CH4 m−2 d−1 and the largest values at US-Los that released the lowest amount of methane. The optimization of the maximum root depth (zroot) results in values ranging between 0.057 and 0.68 with a maximum value at the DK-Nuf site, which is an arctic fen in Greenland. Optimized values for plant-mediated transport efficiency (Tveg) fell between 0.003 and 23.6. The largest Tveg values of 23.6 and 22.3 were obtained for DK-Nuf and US-Wpt, respectively, and the lowest value of 0.003 at DE-Spw. The dimension of water droplets dispersed in the soil depicts the probability of methane-rich bubbles being released to the atmosphere (wsize). The optimized wsize values vary within the range 0.005 and 0.032. And the optimized mixed ratios of methane involved in the ebullition process (mxrCH4) range between 0.06 and 0.53.

Differences between observed and simulated methane fluxes employing initial and optimized parameters are quantified by the RMSD prior (RMSDprior) and posterior (RMSDpost), respectively (Table 5). At sites where methane fluxes were small, such as US-Los and DE-Spw, RMSDpost values are respectively 1.1 and 9.5, whereas at US-Wpt, and RU-Che where monthly mean methane emissions reached up to 550 mg CH4 m−2 d−1, RMSDpost are larger, i.e., respectively 249 and 140. At sites that emitted between 10 and 150 mg CH4 m−2 d−1, RMSD values fluctuate between 4 and 26, and when methane fluxes were between 150 and 400 mg CH4 m−2 d−1, RMSD was 38–80. Performances of the optimization at each site are also evaluated utilizing the relationship (1− RMSDpost/ RMSDprior) × 100, which compares the RMSDprior defined by using the prior values and ranges with the RMSDpost obtained after parameter optimization. It might seem that optimizations are more efficient at sites with low methane emissions than at sites that emitted the most, whereas NRMSD values, which are the RMSDpost normalized by the annual mean of the observed emissions, are close to 1 at each site except for US-Los and DE-Spw for which NRMSDs are 10 and 19, respectively. This suggests that the optimizations are less efficient for sites that emitted the least amount of methane. Direct comparisons during the period of observation between observed and simulated methane emissions are displayed for each site in Figs. 2b, 3b, 4b, and 5b. The temporal and average magnitude are equivalent as in measurements except for US-Wpt and RU-Che for which simulated emissions are much lower than observed emissions.

Figure 2Temporal distribution of methane at sites emitting less than 10 mg CH4 m−2 d−1. (a) Simulated water table position estimated from the soil water content. (b) Simulated (dark line) and observed (gray line) methane emissions released to the atmosphere. (c) Cumulative amount of simulated methane emitted by diffusion, plant-mediated transport, and ebullition. (d) Methane concentration in the soil layers (dark line) and in the snow layers of the model (gray line). (e) On the left is the depth at which simulated methane production is the highest in the soil, scaled to the maximum peat depth. On the right is the amount of simulated methane produced at these depths.


In addition to the mismatch between observed and simulated methane emissions during the observed period, Figs. 2, 3, 4, and 5 show the simulated water table position, the amount of methane that is emitted by diffusion, plant transport, and ebullition, the temporal methane concentration in the soil and in the snow, and the depth at which the largest amount of methane is produced together with the rate of production at that depth. These variables show the consistency of the model regarding peatland functioning. US-Los and DE-Spw emitted less than 10 mg CH4 m−2 d−1, and their simulated water table positions fluctuate below the surface between 10 and 60 cm, while showing a clear seasonal pattern, and are lower in summer than in winter. In winter, simulated emissions are the result of methane diffusion between the soil and the atmosphere, while in spring and summer methane mainly diffuses through aerenchyma of vascular plants. At DE-Spw, the simulated methane concentration in the soil that ranges between 40 and 140 g m−2 is more than 10 times higher than at US-Los, for which the observed concentration barely reaches 5 g CH4 m−2 in the fall. The model simulates methane accumulation in the soil at DE-Spw that stimulates a small release of methane to the atmosphere by ebullition. In the model, the largest production of methane occurs consistently around 20 cm for US-Los and 40 cm for DE-Spw, which is above the simulated water table position. It is commonly expected for methanogenesis to take place below the observed water table position. However, here the simulated water table position is a prognostic variable defined by the cumulative amount of soil water content over the soil column (Figs. S2 and S3). Indeed, in these simulations above the water table position soil moisture is still higher than 80 % (Figs. S4 and S5). At those depths the simulated methane productions reach up to 0.2 and 1.0 g CH4 m−2, respectively, in the summer. In the winter, simulated methane productions are very small, and some methane is diffused in the simulated snowpack covering the peatlands: up to 0.025 g CH4 m−2 at US-Los and 0.17 g CH4 m−2 at DE-Spw. This explains the negative methane flux (Fig. 2c) produced in winter by the model via simulated diffusion of atmospheric methane in the snow cover (Fig. 2d). Then the positive flux that appears in the spring occurs simultaneously with snow melting.

Other sites that emitted less than 150 mg CH4 m−2 d−1 are shown in Fig. 3. Except for CA-Wp1 and US-Bog, during winter these peatlands are nearly inundated in the simulations with a simulated water table position near 10 cm above ground level. CA-Wp1 and US-Bog are respectively fen and bog boreal peatlands, and their simulated water table position is lower than at the other sites. US-Bog is affected by permafrost, which might explain the unexpectedly low position of the simulated water table. At DE-Sfn, methane is mainly transported in the model via vascular plants and by ebullition, whereas at the other sites, simulated methane is predominantly carried via vascular plants only. As for US-Los and DE-Spw, at CA-Wp1, during the winter the simulations show that in the topsoil layers some methane is transferred by diffusion (Fig. 3c) to the snow cover (Fig. 3d). Then a small part of the simulated methane is temporarily stored in the snow (Fig. 3d) and the other part is released to the atmosphere via diffusion (Fig. 3c). More simulated snow accumulated at DE-Sfn, DE-Zrk, CA-Wp1, and US-Bog where up to 0.8–0.04 g CH4 m−2 is temporarily stored in the snow (Fig. 3d). At FR-Lag and DE-Hmm, less methane, with values less than 0.005 g CH4 m−2, is contained in the simulated snow cover (Fig. 3d). As for DE-Spw, at DE-Sfn, simulation results show that up to 140 g CH4 m−2 accumulates in the soil layers of the model during winter and provides sufficient methane to be expelled to the surface by ebullition. In contrast, methane accumulated up to 80 g CH4 m−2 in the soil layers of the model at CA-Wp1 is not sufficient to trigger the methane ebullition process. In all the other sites, methane concentrations in the soil layers of the model are smaller: between 5 and 35 g CH4 m−2. The maximum of simulated methanogenesis takes place steadily at around 20 cm depth at DE-Sfn, FR-Lag, and DE-Hmm, which in winter is about 30 cm under the simulated water table position. At this depth simulated methane production fluctuated at 0.01–0.12 g CH4 m−2. At DE-Sfn, CA-Wp1, and US-Bog, simulations show that in the winter most of the methane is produced at around 75 cm depth, and then in spring and summer the depth of maximum simulated production becomes shallower to reach 20 cm. In early spring at US-Bog, the maximum simulated production is temporarily near the surface at 1 cm depth, which correlates with an increase in methane that accumulates in the simulated snow. At DE-Sfn, the depth at which the maximum simulated production occurred fluctuates more than at both other sites of CA-Wp1 and US-Bog. Unlike CA-Wp1 and US-Bog, during the first 2 years the maximum simulated production deepens at 75 cm when the maximum value of simulated production is reached.

Figure 3Temporal distribution of methane for sites emitting between 10 and 150 mg CH4 m−2 d−1. (a) Simulated water table position estimated from the soil water content. (b) Simulated (dark line) and observed (gray line) methane emissions released to the atmosphere. (c) Cumulative amount of simulated methane emitted by diffusion, plant-mediated transport, and ebullition. (d) Methane concentration in the soil layers (dark line) and in the snow layers (gray line) of the model. (d) On the left is the depth at which simulated methane production is the highest in the soil, scaled to the maximum peat depth. On the right is the amount of simulated methane produced at these depths.


Sites that emitted between 150 and 400 mg CH4 m−2 d−1 are temperate, sub-arctic, and arctic fens (Fig. 4). Simulated water table positions at FI-Lom, DK-Nuf, and PL-Wet are lower in winter than in summer. During the observed period of 3 years, the simulated water table position at PL-Kpt is lower in summer the first and the last year of observations and higher in summer during the second year. In the winter the methane fluxes are stored in the simulated snow cover at FI-Lom (Fig. 4d); therefore, the simulated surface fluxes above the snow are driven by diffusion (Fig. 4c). However, during summer simulated methane fluxes essentially originate from plant-mediated transport. At DK-Nuf, PL-Kpt, and PL-Wet, simulation results show that less methane, with values less than 0.4 g CH4 m−2 d−1, accumulates in the simulated snow during winter (Fig. 4d). Methane is transported by vascular plants in summer at DK-Nuf and PL-Wet, whereas at PL-Kpt simulated methane fluxes are provided by both vascular plants and ebullition. This is consistent with high soil methane concentrations at PL-Kpt during summer that are near 70 g CH4 m−2 the first year and near 90 g CH4 m−2 the last 2 years of observation. In contrast, at FI-Lom simulated soil methane concentrations are near 50 g CH4 m−2 during summer, whereas the winter concentrations are near 80 g CH4 m−2 (Fig. 4d), which is not sufficient to cause methane ebullition (Fig. 4c). Indeed, the ebullition in Eqs. (7) and (8) results from the balance of soil temperature, pressure, gas content, and porosity, which explain the large diversity of methane flux responses by ebullition at each site. At DK-Nuf and PL-Wet simulated soil methane concentrations are less than 10 g CH4 m−2, and therefore ebullition is not produced. At FI-Lom, PL-Kpt, and PL-Wet, the highest simulated methane production rates are maximum at 0.3 g CH4 m−2 d−1, steadily near 20 cm at PL-Wet, at about 20 cm depth in summer, and deepen down to 75 cm depth in winter for the two other sites. While at DK-Nuf the highest simulated methane production rates are lower with values up to 0.08 g CH4 m−2 d−1 and take place around 20 cm in the summer and 40 cm in winter.

Figure 4Temporal distribution of methane for sites emitting between 150 and 400 mg CH4 m−2 d−1. (a) Simulated water table position estimated from the soil water content. (b) Simulated (dark line) and observed (gray line) methane emissions released to the atmosphere. (c) Cumulative amount of simulated methane emitted by diffusion, plant-mediated transport, and ebullition. (d) Methane concentration in the soil layers (dark line) and in the snow layers (gray line) of the model. (e) On the left is the depth at which simulated methane production is the highest in the soil, scaled to the maximum peat depth. On the right is the amount of simulated methane produced at these depths.


The highest simulated methane fluxes of 600 mg CH4 m−2 d−1 were observed at US-Wpt and RU-Che that are respectively a temperate marsh and an arctic tundra site. The simulated water table positions at both sites are lower in the summer than in the winter and vary for US-Wpt between 10 cm above ground and 40 cm below ground level. At RU-Che the prognostic water table depth is very low, i.e., 60 to 90 cm below the soil surface as for US-Bog. Indeed, both sites are underlaid with permafrost, which limits water infiltration to the deepest soil layers and can explain these deeper simulated water table positions. At US-Wpt and RU-Che, site simulations could only provide methane fluxes up to 100 mg CH4 m−2 d−1 despite the expansion of ranges for the optimization of the parameters. These simulated fluxes are entirely transported via vascular plant tissues. During the year of highest fluxes at both sites, simulated methane concentrations are around 0.2 g CH4 m−2 of soil; however, simulated methane concentrations in snow are 10 times lower at the marsh site at 0.3 mg CH4 m−2 than at the tundra site at 3.0–4.0 mg CH4 m−2. At US-Wpt, simulations show that methane is primarily produced around 20 cm depth at a rate of 40–60 mg CH4 m−2 d−1. However, at RU-Che, the simulated methane production rate is higher around 100 mg CH4 m−2 d−1 and occurs at 20 cm depth during summer and a few centimeters below the surface during winter.

Figure 5Temporal distribution of methane for sites emitting more than 400 mg CH4 m−2 d−1. (a) Simulated water table position estimated from the soil water content. (b) Simulated (dark line) and observed (gray line) methane emissions released to the atmosphere. (c) Cumulative amount of simulated methane emitted by diffusion, plant-mediated transport, and ebullition. (d) Methane concentration in the soil layers (dark line) and in the snow layers (gray line) of the model. (e) On the left is the depth at which simulated methane production is the highest in the soil, scaled to the maximum peat depth. On the right is the amount of simulated methane produced at these depths.


3.2 Multi-site optimization (MSO)

For large-scale simulations only one set of parameters is needed for the simulation of methane emissions to achieve the average of each parameter value optimized on-site being commonly employed. Here, a multi-site optimization has been performed for which prior values correspond to the average values of each parameter obtained from the single site optimizations described in Sect. 3.1. This multi-site optimization serves to assess to what extent a multi-site optimization is more efficient than using average values of parameters optimized on-site independently. Multi-site optimized parameter values acquired by using average values of parameters defined at each site and the initial ranges (Table 3) are shown in Table 6. Compared to the prior values, qMG stayed about the same, optimized kMT shifted to values that promote lower oxidation of methane, and near the root area the proportion of methane oxidation Mrox is increased. The plant-mediated transport rate is stimulated by the increase in Tveg to a value of 9 and the rooting depth is about the same at 0.27 for the prior and 0.26 for the posterior. Then the capability of methane ebullition in the model is decreased by the increase in the ebullition threshold deriving from mxrCH4 and the decrease in the probability of bubbles reaching the surface (wsize).

Table 6Multi-site prior and optimized values of methane scheme parameters. Parameter prior values are the average value of the parameters optimized at each site. Parameter descriptions and references are in Table 3.

Download Print Version | Download XLSX

In Table 7, RMSDMS prior constitutes the difference between observed and simulated emissions resulting from average single site optimized parameter values. RMSDMS post is generated from the multi-site optimization of the parameters. For eight sites, posterior values of the RMSDMS are smaller than prior values (RMSDMS prior), thereby reducing the deviation of simulated emissions from the observation. The RMSDMS post values of the six other sites are larger than the RMSDMS prior. Among those RMSDMS values, posterior and prior values are very similar by less than 1 unit for FI-Lom and DK-Nuf. At DE-Hmm, PL-Wet, and US-Bog the differences are lower than 16 units, whereas at RU-Che RMSDMS post is larger by more than 100 units than the RMSDMS prior. NRMSDMS values are larger at US-Los, DE-Spw, and DE-Sfn where methane emissions are lower. At the other sites, the differences of NRMSDMS and NRMSDSS are lower than 1.7 units. These results suggest that for global-scale simulation parameters defined by the multi-site optimization should provide methane emissions estimation with lower uncertainties than when parameters are defined from the average of single site optimization values. Indeed, differences using single site and multi-site optimized parameters, displayed in Fig. 6, are of the same order of magnitude for most sites except for the three sites that emitted the largest amount of methane (PL-Wet, RU-Che, and US-Wpt) and the lowest amount of methane (US-Los, DE-Spw, and DE-Sfn). However, for those six sites methane emission differences between observations and simulations are lower when using multi-site optimized parameters.

A multi-site optimization has also been performed employing extended ranges of parameter values that are enlarged to the maximum and minimum values obtained for the single site optimizations (Tables S4 to S6 and Fig. S9). Despite a different set of parameters being defined (Table S3), discrepancies between observed and simulated emissions (Tables S5 and S6 and Fig. S10) are similar to the ones obtained using default parameter ranges.

Table 7Discrepancies between observed and simulated methane emissions are quantified by the root mean square difference (RMSD) approach. Minimization efficiency of the multi-site optimization is indicated by the relationship between the prior using average values of parameters optimized by the single site optimization and posterior RMSDMS as (1− RMSDMSpost/RMSDMS prior) × 100. Normalized root mean square difference (NRMSDMS) is defined by the RMSDMS posterior normalized by the annual mean of observed methane emissions in Table 5.

Download Print Version | Download XLSX

Figure 6Simulated and observed (gray line) methane emissions using single site (dashed dark line) and multi-site (solid dark line) optimized parameters.


4 Discussion

4.1 Parameterization sensitivity

Sensitivity analyses were previously performed to assess methane emission model responsiveness to parameter values (Meng et al., 2012; Riley et al., 2011; Spahni et al., 2011; Wania et al., 2009; Zhu et al., 2014). These studies (van Huissteden et al., 2009; Riley et al., 2011) suggested that temperature dependency of methanogenesis is the most influential parameter affecting methane production, whereas methane emissions are mostly sensitive to oxidation and plant transport. Indeed, in large-scale models such as CLM4Me, LPJ-GUESS, LPX-Bern, CNRM, and ORCHIDEE (Potter, 1997; Riley et al., 2011; Khvorostyanov et al., 2008b; Wania et al., 2009, 2010; Zhu et al., 2014; Morel et al., 2019) methane production results from anoxic decomposition of soil organic matter, the rate of which is constrained by the soil oxic and anoxic decomposition ratio (qMG). Therefore, the methanogenesis rate is driven by the same variables as the oxic decomposition that depends on soil temperature and primary production. This ratio was first established from experimental studies that determine the microbial production ratio CO2 to CH4 (Potter et al., 1996; Segers, 1998) for various water table positions. These ratio values were found to be between 0.58 and 10 000. Because of this wide range of values, process-based models employed this CO2-to-CH4 ratio as an adjustable parameter that is weighted by environmental factors such as soil moisture and temperature. Wania et al. (2009) performed a sensitivity analysis study of the LPJ-WHyMe model using seven sites in which the multi-site optimization value of the CO2/CH4 ratio was defined at 10, while other models such as CLM4Me use a value of 5. Khvorostyanov et al. (2008a) and Morel et al. (2019) respectively used qMG values of 9 and 10 to simulate methane emissions from arctic peatlands. Therefore, in the present study at first qMG was optimized in the range of 9–11, and then this range was enlarged only for sites that underestimate methane emissions. Results show that for 13 sites out of 14, qMG values range 8.0–10.7 for the single site optimization approach, and using the multi-site approach a value of 9.6 was found. As in the previous sensitivity analysis studies (Riley et al., 2011) lower qMG values were obtained at sites located at the highest latitudes.

After methanogenesis, methane is mobilized in pores and ultimately emitted to the atmosphere or is oxidized by methanotrophs depending on whether methane travels along the anoxic or the oxic parts of the soil. In large-scale models, methanotrophy is formulated employing a Michaelis–Menten or a first-order kinetic framework based on soil methane and oxygen content (Morel et al., 2019). These formulations are then driven by the oxidation rate, the values of which vary from a few hours to days. In the present work, we employed the first-order kinetic formulation of Khvorostyanov et al. (2008a) that is driven by methane and oxygen content. Optimization of the oxidation rate leads to values that are spread over the full range of 1 to 5 per day. This is consistent with the review paper of Smith et al. (2003), highlighting the fact that methanotrophy is more sensitive to soil moisture than soil temperature and that there is a direct link between methane oxidation rate and gas diffusivity. Thus, the optimization of the oxidation rate results from the balance between model inputs and outputs that are respectively available methane and oxygen substrates as well as methane fluxes, which explain this large variability in oxidation rate. In addition, in our model, snow is considered in the diffusion scheme, which in part controls diffusivity of oxygen from the atmosphere to the ground in winter (e.g., Fig. 2c).

Methane emissions mediated by vascular plants result from series of processes that include (1) the diffusion and advective transport of methane and oxygen in aerenchyma tissues, (2) autotrophic respiration of a fraction of oxygen transiting in aerenchyma of vascular plants (Colmer, 2003; Nielsen et al., 2017), (3) methane production by microbial decomposition of plant exudates, and (4) methane oxidation by exudates and by remaining oxygen at the root level brought through aerenchyma that increase methanotroph activities. Modeling these processes requires (1) understanding and quantifying them (Kaiser et al., 2017; Raivonen et al., 2017; Riley et al., 2011; Wania, 2007) as well as (2) evaluating the average density of vascular plants that are capable of significant gas transport across ecosystems. While a significant number of studies provide insight on gas exchanges through vascular plants, densities of vascular plants with aerenchyma in peatlands are poorly characterized. In the most recent models, formulations of various complexity were used to simulate vegetation-mediated gas transport considering mainly CH4 and O2 (Kaiser et al., 2017; Morel et al., 2019; Raivonen et al., 2017; Riley et al., 2011; Wania, 2007). These schemes consider plant transport at the scale of the plant and are based on gas concentration gradients between the atmosphere and the soil as well as some plant traits and properties such as plant height, root diameters, aerenchyma porosity, and permeability. Because of the biodiversity of peatlands, calibration of parameters accounting for plant traits and properties of each plant species or family is a cumbersome achievement, and the lack of quantification of aerenchymatous plants at the scale of the ecosystem reduces the benefit in considering these characteristics. In the present scheme, vegetation transport of methane is simulated by employing the rather simple scheme of Walter and Heimann (2001) that is driven by the rooting depth (zroot) of vascular plants with aerenchyma and by the proportion of methane that is oxidized by the rhizosphere (Mrox). Optimized zroot values at sites ranges between 6 and 68 cm depth with the average depth defined at 26 cm, which is also the value obtained using the multi-site approach. These values are consistent with values utilized by Walter and Heimann (2001) that ranged between 0 and 74 cm. It could be expected for zroot to be set near the depth of maximum methanogenesis as is the case at DE-Sfn where zroot is defined at 40 cm. Half of the sites have a zroot defined between 10 and 60 cm above the depth of maximum methanogenesis, and the other remaining values are established between 10 and 50 cm below the depth of maximum methanogenesis. In the rhizosphere methane can also be oxidized at a rate (Mrox) that is independent of the rate of methanotrophy. Results of the optimization at site level provided Mrox values that are scattered over the range of 0 to 1, with the highest values of 0.99 at site US-Los, which emitted the least methane. The lowest value of 0.003 was found at RU-Che; the site emitted the largest amount of methane. Two trends can be distinguished; for sites that emitted less than 150 mg CH4 m−2 d−1 an average of 60 % of methane is oxidized by the rhizosphere against 22 % at sites emitting more. Across all sites the average proportion of methane oxidized is 44 %, whereas the optimized value obtained with the multi-site approach is 70 %. In previous models, Zhuang et al. (2004) and Wania et al. (2010) employed a fixed value of 40 % and 50 %, respectively, at the global scale. With a more realistic and complex formulation in CLM4Me, Riley et al. (2011) estimated that 60 % of methane that would have been transferred to the atmosphere by aerenchyma tissues is instead oxidized by the rhizosphere. Tveg was introduced by Walter et al. (1996) to describe the density of plants and their efficiency in methane transport for site estimation. It is an adjustable parameter that was scaled to be between 0 and 15, with lower values for ecosystems dominated by trees and shrubs and the highest values for ecosystems dominated by grasses and sedges. For our 14 sites, optimization at sites established Tveg values between 0.003 and 24 with an average value of 7 and an optimized value at 8.6 for the multi-site approach. Only two values have been defined above 10 at US-Wpt and DK-Nuf, which are two sites that are limited in methane substrates in the model; this explains these high values of Tveg.

When methane is significantly produced in the soil, the accumulation of methane in the water-saturated pores involves the formation of methane-rich bubbles that will migrate in the soil layers and eventually deliver methane to the atmosphere. This flux of methane is commonly prompted in land surface models by the amount of methane that is no longer soluble in saturated water-filled pores. This excess amount is defined here from the mixing ratio (mrxCH4) of methane in bubbles. Then this volumetric content of methane is converted to methane concentration per soil volume in each layer depending on soil temperature and pressure. The optimization of mxrCH4 at each site leads to values ranging between 3 % and 53 % with a mean value at 24 %, whereas the multi-site optimization evaluates mxrCH4 at 57 %. It has been suggested in the literature that the methane partial pressure is sensitive to fluctuations of the hydrostatic and atmospheric pressure (Tokida et al., 2007b) and of the water table position (Fechner-Levy and Hemond, 1996). Vegetation also impacts the ebullition flux by increasing substrate availability and by indefinitely stabilizing bubbles around roots (Klapstein et al., 2014). Migration of methane-rich bubbles to the soil surface can be modeled as an instantaneous transport to the atmosphere or to upper layers or by an advective layer-by-layer transport. Here we considered the probability of a methane-rich bubble reaching the surface depending on the connectivity between water-filled pores (wsize). Khvorostyanov et al. (2008a) defined wsize at 1 cm, which establishes a probability of 1 at the surface that decreases to zero at 1.5 m depth when soil is saturated. Probability increases when wsize increases and quickly decreases when soil moisture decreases. In the present study, at each site wsize is optimized to values of 0.05–3 cm. At most sites, optimized wsize values are near or below 1 cm except for US-Los, DK-Nuf, and RU-Che. This might be explained by the low methane concentration in the model soil layers at these sites, which annihilates possible emissions by ebullition in the model. The average value across sites corresponds to the same value determined by Khvorostyanov et al. (2008a) at 0.9 cm. A lower value is obtained for the multi-site optimization of 0.2 cm, which reduces occurrence of methane flux by ebullition in our model.

4.2 Methane sources

Soil and litter organic carbon and plant exudates are recognized to be the main substrates for methanogenesis (Chang et al., 2019; Riley et al., 2011; Whalen, 2005). Recent work of Hopple et al. (2019) demonstrates that dissolved organic carbon (DOC) also significantly contributes to anoxic decomposition in peatlands. Some field studies suggested that high-latitude methanogenesis can be substrate-limited (Chang et al., 2019; Riley et al., 2011; Whalen, 2005). In large-scale models, soil organic carbon (SOC) is considered to be the primary source of methane; however, in order to increase the rate of methanogenesis, labile organic matter, such as litter carbon and plant exudates, is directly combined with soil carbon, bypassing oxic decomposition processes to account for them as substrates for the methane production scheme (Morel et al., 2019; Khvonostyanov et al., 2008b). In the present study, SOC is the only substrate for methanogenesis for which total soil carbon stock and maximum peat depth have been adjusted to observation data at each site (Table 2). Simulation results show that at sites that emitted more than 400 mg CH4 m−2 d−1, which are US-WPT and RU-Che, methane emissions are lower than expected, reflecting the lack of substrate for methanogenesis. Indeed, in land surface models, soil carbon is distributed in three types: the active, slow, and passive pool. The active pool features labile SOC, whereas the slow and passive pools exert more stable SOC with slower decomposition rates. Figures 2e to 5e display the depth of maximum methane production and reveal that the deepest methane production depth is 0.75 m in all the simulation results. Integrated SOC accumulated up to 0.75 m by our model for each site is reported in Table 8. These carbon stocks correspond to available substrate for methanogenesis occurring at a depth lower than 0.75 m. The lowest carbon stocks were obtained at US-Los, CA-Wp1, PL-Wet, US-Wpt, and RU-Che with a total SOC lower than 50 kg m−2. Unlike the other sites, the active SOC contents at US-Wpt and RU-Che are very small at 4 and 3.5 kg m−2, respectively, which limits methane production in the model. At both sites, simulated vertical carbon contents were constrained using observed soil bulk density and the carbon accumulation model described in Qiu et al. (2019). Khvorostyanov et al. (2008b) previously performed site simulation at RU-Che in which they prescribed 15 g C m−2 yr−1 of root exudates that was added to the active SOC, leading to emissions up to 300 mg m−2 d−1. As US-Wpt is a marsh it is expected to have a lower total SOC than the other peatland sites. It is also expected that root exudates and DOC in pore water as well as in aboveground reservoirs significantly contribute to methanogenesis, which is not explicitly considered in the present version of the model.

Table 8Integrated simulated soil organic carbon content of peatland sites up to 0.75 m depth.

Download Print Version | Download XLSX

4.3 Methane fluxes

Sensitivity of methane fluxes to model parameters was evaluated by comparing annual methane emissions obtained by employing single site (SS) and multi-site (MS) optimized parameters. Table 9 reports annual observed and simulated methane fluxes as well as the contributions among the three types of methane transport, i.e., diffusion, ebullition, and plant-mediated. Considering all 14 sites, average annual methane emissions for the observed values are 18 ± 18 g m−2 yr−1 and 9 ± 6 as well as 25 ± 38 g m−2 yr−1 for simulations using SS and MS optimized parameters, respectively. Diffusion of methane in the topsoil layers of the model was minor compared to the other emissions and appeared to act as a sink of methane rather than a source. Plant-mediated transport (PMT) was the largest simulated flux during the plant's growth period. For SSO simulations these PMT fluxes represent between 52 % and 74 % of the total fluxes at US-Los, DE-Spw, DE-Sfn, and PL-Kpt and more than 97 % at all the other sites, whereas for MSO simulations PMT fluxes are all higher than 98 %. Given that diffusion released small amounts of methane to the atmosphere, remaining fluxes are emitted by ebullition. The largest ebullition fluxes were obtained in SSO simulations, whereas less methane was released by ebullition in MSO simulations. For about half of the sites, 3 %–11 % of fluxes were furnished via ebullition and less than 1 % at the other sites using SSO parameter values. In simulations employing MSO parameter values, ebullition contributed to less than 2 % of the total fluxes at each site.

Table 9Yearly methane emissions defined from the observed data (Obs) as well as simulations employing optimized parameters obtained by the single site optimization (SSO) and multi-site optimization (MSO). The methane fluxes combine methane emitted by diffusion, plant-mediated transport, and ebullition.

Download Print Version | Download XLSX

Discrepancies between the observation data and the SSO and MSO simulations are displayed in Fig. 6. At sites that emitted the largest amount of methane e.i. PL-Wet, RU-Che, and US-Wpt, SSO and MSO simulations were underestimated up to 46 and 53 g CH4 m−2 yr−1, respectively (Figs. S6 to S8). At the other sites when using SSO parameters methane emissions were still underestimated even though this was only about 7 g CH4 m−2 yr−1. In MSO simulations only the three sites of DE-Hmm, FI-Lom, and DK-Nuf underestimated methane emissions of 11 g CH4 m−2 yr−1 compared to observation data. Simulations that display, in Fig. 7, an overestimation of methane emissions were all performed using MSO parameters. At DE-Spw and DE-Sfn methane emissions were overestimated by 118 and 95 g CH4 m−2 yr−1. This large excess of methane emissions results from a significant increase in the parameter Tveg between the SSO and MSO. Indeed, optimized Tveg values at these sites are 0.003 and 0.1 when optimized at site level, whereas it was defined at 8.6 with the multi-site approach. In the model, Tveg established the magnitude of plant-mediated fluxes, which are constrained by soil methane content, plant growth, and root expansion in the soil. This shows that for peatlands where methanogenesis is not substrate-limited, Tveg is a key parameter to evaluate methane fluxes. Other sites that display an overestimation of methane emissions using MSO parameters are US-Los, CA-Wp1, and PL-Kpt. For these sites the excess of emissions compared to the observations only extends up to 12 g CH4 m−2 yr−1. Across sites, differences between observed emissions and simulated emissions employing SSO parameters average around 9 g CH4 m−2 yr−1 of methane deficiency. On the contrary, emissions obtained with MSO parameters are in excess of about 5 g CH4 m−2 yr−1 on average compared to observations. Average differences between observations and simulation results significantly decrease to −1.2 and 0.5 g CH4 m−2 yr−1 for SSO and MSO simulations when excluding sites that emitted more than 300 and less than 20 mg CH4 m−2 d−1, i.e., PL-Wet, RU-Che, and US-Wpt for the SSO simulations and DE-Spw, DE-Sfn, PL-Wet, RU-Che, and US-Wpt for the MSO simulations. This shows that the model is better constrained at sites emitting between 20 and 300 mg CH4 m−2 d−1.

Figure 7Difference in annual methane emissions defined between the observed data (Obs) and simulations employing optimized parameters obtained by the single site optimization (SSO) and by multi-site optimization (MSO).


Average methane emissions estimated from these 14 sites can be utilized to roughly calculate emissions from peatlands located northern of 30 N. In Qiu et al. (2019), northern peatland extent has been estimated using ORCHIDEE_PEAT v2.0 and compared with three other peatland inventories and soil data (Batjes, 2016; Joosten, 2009; Xu et al., 2018). All four estimates of northern peatland areas range between 2823 and 3896 × 103 km2. Assessment of methane emissions for these northern peatland areas estimated using the average fluxes from measurements yields annual methane fluxes of 51–71 Tg CH4 yr−1 (Table 9). These annual fluxes are in good agreement with annual methane emissions determined from upscaling of flux measurements of 44–54 Tg CH4 yr−1 by Zhu et al. (2013). Estimates of annual methane fluxes obtained from the SSO and MSO simulations lead to values of 25–35 and 70–96 Tg CH4 yr−1, respectively. Estimates from SSO simulations are consistent with annual methane emissions calculated from inversion models (Bruhwiler et al., 2014; Spahni et al., 2011) and other process-based models (Chen et al., 2015; Peltola et al., 2019; Treat et al., 2018; Zhang et al., 2016). Annual methane emissions assessed from MSO simulations are above the upper range of annual methane fluxes provided by the global methane budget for natural wetlands located north of 30 N of 12–61 Tg CH4 yr−1 for a bottom-up approach and 31–64 Tg CH4 yr−1 for a top-down approach (Saunois et al., 2020).

5 Conclusion

The methane model developed by Khvorostyanov et al. (2008a) has been modified to encompass northern peatlands and permafrost features embedded in the most recent version of ORCHIDEE-PEAT v2.0. This modified version, ORCHIDEE-PCH4, which was used in this study, integrates a vertical discretization of oxic and anoxic decomposition of soil organic carbon of northern peatlands and subsequent methane production, oxidation and transport by vascular plants, and ebullition and diffusion in soil and snow layers. A sensitivity analysis of methane emissions was performed on changes of seven model parameters optimized with site-level measurements of 14 sites located north of 41 N on the Eurasian and American continents. The ORCHIDEE data assimilation system (Bastrikov et al., 2018) with a genetic algorithm for random search has been successfully employed to optimize these seven parameters at each site and consider methane emissions from all sites simultaneously. Our results show that, as in previous methane emissions models (Meng et al., 2012; Riley et al., 2011; Spahni et al., 2011; Wania et al., 2009; Zhu et al., 2014), simulated methanogenesis is strongly correlated with simulated soil temperature and moisture content, whereas methane emissions are more strongly correlated with plant-mediated fluxes and soil methane oxidation proportion. Surprisingly, a weak correlation has been established between the observed water table positions and the prognostic water table positions established from simulated soil moisture content. A correlation between soil moisture content and water table position in the field is needed to improve representation of the water table position in models.

Single site optimization results highlighted the fact that the depth of the highest methane production fluctuates between 20 cm during the warmer season and 75 cm during the cold season. This demonstrates the sensitivity of methanogenesis to soil temperature and provides insight on the extent to which methanogenesis takes place in the soil layers. This also serves to identify sites that are substrate-limited and to emphasize the need for global-scale models to consider dissolved organic matter as a source of methane substrate. Indeed, in some site simulation studies prescribed methane substrate originating from litter decomposition or plant exudates was added to soil organic content in order to balance out the lack of labile substrate. In the scheme of ORCHIDEE-PCH4, the addition of methane diffusion in the snow layers during winter exposes the potential of snow to delay methane emissions coming from the soil.

Optimization of parameters simultaneously employing methane emissions from all 14 sites produce a reduction in the rate of methanotrophy and in methane transport in the soil by ebullition, promoting methane oxidation at the root level and transport of methane by vascular plants. These involve a large overestimation of sites emitting small amounts of methane. Nonetheless, on average methane emissions simulated employing the multi-site optimization approach are only overestimated by about 5 g CH4 m−2 yr−1 because the overestimation of low emitting sites is counterbalanced by the high emitting sites that are limited in methane substrates. In contrast, average methane emissions obtained from the simulations using parameters from the single site optimization underestimate the average observed fluxes by 9 g CH4 m−2 yr−1. Nevertheless, extrapolation of these average methane emissions to northern peatland areas reveals that emissions estimated from the multi-site simulations are much larger than emissions estimated from other peatland process-based models and inventories, whereas emissions calculated from the single site optimizations are in good agreement with other estimates. This demonstrates the complexity of the interactions of the methane cycle with environmental conditions considered at various scales and the need for more detailed on-site studies.

Code availability

The source code (, Salmon, 2021) is available online via (last access: 24 March 2022). Readers interested in running the model should follow the guidelines at (last access: 24 March 2022).

The optimization tool is available through a dedicated website for data assimilation with ORCHIDEE (, Bastrikov, 2018).

Data availability

Measured eddy covariance fluxes and related meteorological data can be obtained from the SNO-T (, Doc SNO-T, 2022) for FR-Lag; from the European flux database cluster (, European Fluxes Database Cluster, 2022) for DE-Hmm (, last access: 25 March 2022), DE-Spw (, last access: 25 March 2022), PL-Kpt (, last access: 25 March 2022), PL-Wet (, last access: 25 March 2022), and RU-Che (, last access: 25 March 2022); from the FLUXNET database (, ORNL DAAC, 2022) for DE-Sfn, DE-Zrk, FI-Lom, and DK-Nuf; from the AmeriFlux database (, AmeriFlux, 2022) for CA-Wp1 (, last access: 25 March 2022), US-Bog (, last access: 25 March 2022), US-Los (, last access: 25 March 2022), and US-Wpt (, last access: 25 March 2022); and from the investigators upon request. Model outputs are available upon request.


The supplement related to this article is available online at:

Author contributions

ES revised and modified the implementation of the methane module in ORCHIDEE-PEAT, performed model optimization simulations in ORCHIDEE-PCH4 employing the ORCHIDEE data assimilation system, investigated simulation results, and prepared the paper with contributions from all co-authors. FJ, BG, CG, PC, SG, and FLD conceptualized, secured funding for, and supervised the project. CQ, DZ, and BG provided ORCHIDEE-PEAT and assisted with model setup and development. LJ, FJ, and BG assisted in implementing ORCHIDEE-PCH4 and investigating simulation results. VB and PP provided the model and the expertise on the ORCHIDEE data assimilation system and contributed to the interpretation of simulation results. CG, SG, FLD, MA, MSBH, JC, BHC, HC, CWE, ESE, LBF, KF, DH, JK, OK, NK, LK, AL, LM, WP, TS, and KZ produced and provided the quality field dataset employed to constrain and validate ORCHIDEE-PEAT and ORCHIDEE-PCH4.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The modeling work was supported by the European Union's Horizon 2020 Project CRESCENDO under contract 641816 and Labex VOLTAIRE ANR-10-LABX-100-01. The authors acknowledge the support of staff at each site. Research at US-Los was supported by the AmeriFlux Network Management Project under contract no. 7544821 to the ChEAS core site cluster. Funding for the measurements in Biebrza National Park was provided by the Polish National Science Centre under projects UMO-2015/17/B/ST10/02187 and UMO-2020/37/B/ST10/01219. US-Bog was supported by the United States National Science Foundation: NSF OPP 1107892, 1503912, 1936712. Lawrence B. Flanagan acknowledges funding from the Natural Sciences and Engineering Research Council of Canada (NSERC), the FLUXNET-Canada Network (NSERC, the Canadian Foundation for Climate and Atmospheric Sciences – CFCAS, and BIOCAP Canada), and the Canadian Carbon Program (CFCAS). Fieldwork at FR-Lag was funded as part of the Labex VOLTAIRE and the PIVOTS project of the Région Centre – Val de Loire (ARD 2020 program and CPER 2015–2020) in the framework of the French Peatland Observatory, SNO Tourbières, endorsed by CNRS-INSU. Research work at DE-Hmm was funded by the Deutsche Forschungsgemeinschaft under Germany's Excellence Strategy – EXC 177 “CliSAP – Integrated Climate System Analysis and Prediction” – contributing to the Center for Earth System Research and Sustainability (CEN) of Universität Hamburg. The work at PL-Wet is based on use of Large Research Infrastructure CzeCOS supported by the Ministry of Education, Youth and Sports of CR within the CzeCOS program under grant number LM2018123. Natalia Kowalska acknowledges the support by SustES – Adaptation strategies for sustainable ecosystem services and food security under adverse environmental conditions (CZ.02.1.01/0.0/0.0/16_019/0000797).

Financial support

This research has been supported by Horizon 2020 (CRESCENDO (grant no. 641816)) and Labex VOLTAIRE (grant no. ANR-10-LABX-100-01).​​​​​​​

Review statement

This paper was edited by Carlos Sierra and reviewed by two anonymous referees.


AmeriFlux: Homepage,, last access: 25 March 2022. 

Allan, W., Struthers, H., and Lowe, D. C.: Methane carbon isotope effects caused by atomic chlorine in the marine boundary layer: Global model results compared with Southern Hemisphere measurements, J. Geophys. Res.-Atmos., 112, D04306,, 2007. 

Anav, A., Friedlingstein, P., Kidston, M., Bopp, L., Ciais, P., Cox, P., Jones, C., Jung, M., Myneni, R., and Zhu, Z.: Evaluating the land and ocean components of the global carbon cycle in the CMIP5 earth system models, J. Climate, 26, 6801–6843,, 2013. 

Bastrikov, V.: ORCHIDEE Data Assimilation Systems, ORCHIDEE Land Surface Model [code], (last access: 25 March 2022), 2018. 

Bastrikov, V., MacBean, N., Bacour, C., Santaren, D., Kuppel, S., and Peylin, P.: Land surface model parameter optimisation using in situ flux data: comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2), Geosci. Model Dev., 11, 4739–4754,, 2018. 

Batjes, N. H.: Harmonized soil property values for broad-scale modelling (WISE30sec) with estimates of global soil carbon stocks, Geoderma, 269, 61–68,, 2016. 

Benscoter, B. W., Thompson, D. K., Waddington, J. M., Flannigan, M. D., Wotton, B. M., de Groot, W. J., and Turetsky, M. R.: Interactive effects of vegetation, soil moisture and bulk density on depth of burning of thick organic soils, Int. J. Wildl. Fire, 20, 418–429,, 2011. 

Blake, D. R., Mayer, E. W., Tyler, S. C., Makide, Y., Montague, D. C., and Rowland, F. S.: Global increase in atmospheric methane concentrations between 1978 and 1980, Geophys. Res. Lett., 9, 477–480,, 1982.​​​​​​​ 

Blodau, C.: Carbon cycling in peatlands – A review of processes and controls, Environ. Rev., 10, 111–134,, 2002. 

Bradley-Cook, J. I. and Virginia, R. A.: Soil carbon storage, respiration potential, and organic matter quality across an age and climate gradient in southwestern Greenland, Polar Biol., 39, 1283–1295,, 2016. 

Bridgham, S. D., Cadillo-Quiroz, H., Keller, J. K., and Zhuang, Q.: Methane emissions from wetlands: Biogeochemical, microbial, and modeling perspectives from local to global scales, Glob. Change Biol., 19, 1325–1346,, 2013. 

Bruhwiler, L., Dlugokencky, E., Masarie, K., Ishizawa, M., Andrews, A., Miller, J., Sweeney, C., Tans, P., and Worthy, D.: CarbonTracker-CH4: an assimilation system for estimating emissions of atmospheric methane, Atmos. Chem. Phys., 14, 8269–8293,, 2014. 

Chang, K. Y., Riley, W. J., Brodie, E. L., McCalley, C. K., Crill, P. M., and Grant, R. F.: Methane Production Pathway Regulated Proximally by Substrate Availability and Distally by Temperature in a High-Latitude Mire Complex, J. Geophys. Res.-Biogeo., 124, 3057–3074,, 2019. 

Chanton, J. P. and Whiting, G. J.: Trace gas exchange in freshwater and coastal marine environments: ebullition and trans-port by plants, in: Biogenic trace gases: measuring emissions from soil and water, edited by: Matson, P. A. and Harriss, R. C., Blackwell Science Ltd, ISBN 0-632-03641-9, 1995. 

Chason, D. B. and Siegel, D. I.: Hydraulic conductivity and related physical properties of peat, Lost River Peatland, northern Minnesota., Soil Sci., 142, 91–99, 1986. 

Chen, X., Bohn, T. J., and Lettenmaier, D. P.: Model estimates of climate controls on pan-Arctic wetland methane emissions, Biogeosciences, 12, 6259–6277,, 2015. 

Chu, H., Chen, J., Gottgens, J. F., Ouyang, Z., John, R., Czajkowski, K., and Becker, R.: Net ecosystem methane and carbon dioxide exchanges in a Lake Erie coastal marsh and a nearby cropland, J. Geophys. Res.-Biogeo., 119, 722–740,, 2014. 

Colmer, T. D.: Long-distance transport of gases in plants: A perspective on internal aeration and radial oxygen loss from roots, Plant Cell Environ., 26, 17–36,, 2003. 

Dettmann, U., Bechtold, M., Frahm, E., and Tiemeyer, B.: On the applicability of unimodal and bimodal van Genuchten–Mualem based models to peat and other organic soils under evaporation conditions, J. Hydrol., 515, 103–115,, 2014. 

Dlugokencky, E.: Trends in atmospheric methane, NOAA/GML​​​​​​​,, last access: January 2021. 

Doc SNO-T: Documentation du SI SNO-Tourbières (SNO-T), Doc SNO-T [data set],, last access: 25 March 2022. 

Dutta, K., Schuur, E. A. G., Neff, J. C., and Zimov, S. A.: Potential carbon release from permafrost soils of Northeastern Siberia, Glob. Change Biol., 12, 2336–2351,, 2006. 

Duval, B. and Goodwin, S.: Methane production and release from two New England peatlands, Int. Microbiol., 3, 89–96,, 2000. 

Etheridge, D. M., Steele, L. P., Francey, R. J., and Langenfelds, R. L.: Atmospheric methane between 1000 A.D. and present: Evidence of anthropogenic emissions and climatic variability, J. Geophys. Res.-Atmos., 103, 15979–15993,, 1998. 

Etminan, M., Myhre, G., Highwood, E. J., and Shine, K. P.: Radiative forcing of carbon dioxide, methane, and nitrous oxide: A significant revision of the methane radiative forcing, Geophys. Res. Lett., 43, 12614–12623,, 2016. 

European Fluxes Database Cluster: Sites list, European Fluxes Database Cluster [data set],, last access: 25 March 2022. 

Fechner-Levy, E. J. and Hemond, H. F.: Trapped methane volume and potential effects on methane ebullition in a northern peatland, Limnol. Oceanogr., 41, 1375–1383,, 1996. 

Gogo, S., Laggoun-Défarge, F., Delarue, F., and Lottier, N.: Invasion of a Sphagnum-peatland by Betula spp and Molinia caerulea impacts organic matter biochemistry. Implications for carbon and nutrient cycling, Biogeochemistry, 106, 53–69,, 2011. 

Goldberg, D.: Genetic Algorithms, in: Search, Optimization and Ma- chine Learning, Addison-Wesley, 432 pp., ISBN 978-0-201-15767-3, 1989. 

Gorham, E.: Northern peatlands: role in the carbon cycle and probable responses to climatic warming, Ecol. Appl., 1, 182–195,, 1991. 

Guimberteau, M., Zhu, D., Maignan, F., Huang, Y., Yue, C., Dantec-Nédélec, S., Ottlé, C., Jornet-Puig, A., Bastos, A., Laurent, P., Goll, D., Bowring, S., Chang, J., Guenet, B., Tifafi, M., Peng, S., Krinner, G., Ducharne, A., Wang, F., Wang, T., Wang, X., Wang, Y., Yin, Z., Lauerwald, R., Joetzjer, E., Qiu, C., Kim, H., and Ciais, P.: ORCHIDEE-MICT (v8.4.1), a land surface model for the high latitudes: model description and validation, Geosci. Model Dev., 11, 121–163,, 2018. 

Haupt, R. L. and Haupt, S. E.: Practical Genetic Algorithms, John Wiley & Sons, Inc., ISBN 9780471671749,, 2004. 

Hillel, D.: Introduction to soil physics, Academic Press Inc., New york, 392 pp., ISBN 9780123485205, 1982. 

Hodgman, C. D.: Handbook of Chemistry and Physics, Chemical Rubber Publishing Co., Cleveland, 1951 pp.​​​​​​​, ISSN 0363-3055,, 1936. 

Hommeltenberg, J., Schmid, H. P., Drösler, M., and Werle, P.: Can a bog drained for forestry be a stronger carbon sink than a natural bog forest?, Biogeosciences, 11, 3477–3493,, 2014. 

Hopple, A. M., Pfeifer-Meister, L., Zalman, C. A., Keller, J. K., Tfaily, M. M., Wilson, R. M., Chanton, J. P., and Bridgham, S. D.: Does dissolved organic matter or solid peat fuel anaerobic respiration in peatlands?, Geoderma, 349, 79–87,, 2019. 

Jaszczynski, J.: The relationship between dissolved organic carbon and hydro-climatic factors in peat-muck soil, J. Water L. Dev., 24, 27–33,​​​​​​​​​​​​​​, 2015. 

Joosten, H.: The Global Peatland CO2 Picture. Peatland status and drainage related emissions in all countries of the world, Wetlands International,, (last access: 25 March 2022), 2009. 

Kaiser, S., Göckede, M., Castro-Morales, K., Knoblauch, C., Ekici, A., Kleinen, T., Zubrzycki, S., Sachs, T., Wille, C., and Beer, C.: Process-based modelling of the methane balance in periglacial landscapes (JSBACH-methane), Geosci. Model Dev., 10, 333–358,, 2017. 

Khvorostyanov, D. V., Krinner, G., Ciais, P., Heimann, M., and Zimov, S. A.: Vulnerability of permafrost carbon to global warming. Part I: Model description and role of heat generated by organic matter decomposition, Tellus B, 60, 250–264,, 2008a. 

Khvorostyanov, D. V., Ciais, P., Krinner, G., Zimov, S. A., Corradi, C., and Guggenberger, G.: Vulnerability of permafrost carbon to global warming. Part II: Sensitivity of permafrost carbon stock to global warming, Tellus B, 60, 265–275,, 2008b. 

Kirschke, S., Bousquet, P., Ciais, P., Saunois, M., Canadell, J. G., Dlugokencky, E. J., Bergamaschi, P., Bergmann, D., Blake, D. R., Bruhwiler, L., Cameron-Smith, P., Castaldi, S., Chevallier, F., Feng, L., Fraser, A., Heimann, M., Hodson, E. L., Houweling, S., Josse, B., Fraser, P. J., Krummel, P. B., Lamarque, J. F., Langenfelds, R. L., Le Quéré, C., Naik, V., O'doherty, S., Palmer, P. I., Pison, I., Plummer, D., Poulter, B., Prinn, R. G., Rigby, M., Ringeval, B., Santini, M., Schmidt, M., Shindell, D. T., Simpson, I. J., Spahni, R., Steele, L. P., Strode, S. A., Sudo, K., Szopa, S., Van Der Werf, G. R., Voulgarakis, A., Van Weele, M., Weiss, R. F., Williams, J. E., and Zeng, G.: Three decades of global methane sources and sinks, Nat. Geosci., 6, 813–823,, 2013. 

Klapstein, S. J., Turetsky, M. R., McGuire, A. D., Harden, J. W., Czimczik, C. I., Xu, X., Chanton, J. P., and Waddington, J. M.: Controls on methane released through ebullition in peatlands affected by permafrost degradation, J. Geophys. Res.-Biogeo., 119, 418–431,, 2014. 

Krüger, M., Eller, G., Conrad, R., and Frenzel, P.: Seasonal variation in pathways of CH4 production and in CH4 oxidation in rice fields determined by stable carbon isotopes and specific inhibitors, Glob. Change Biol.,, 2002. 

Largeron, C., Krinner, G., Ciais, P., and Brutel-Vuilmet, C.: Implementing northern peatlands in a global land surface model: description and evaluation in the ORCHIDEE high-latitude version model (ORC-HL-PEAT), Geosci. Model Dev., 11, 3279–3297,, 2018. 

Le Mer, J. and Roger, P.: Production, oxidation, emission and consumption of methane by soils: A review, Eur. J. Soil Biol., 37, 25–50,, 2001. 

Leroy, F., Gogo, S., Guimbaud, C., Francez, A. J., Zocatelli, R., Défarge, C., Bernard-Jannin, L., Hu, Z., and Laggoun-Défarge, F.: Response of C and N cycles to N fertilization in Sphagnum and Molinia-dominated peat mesocosms, J. Environ. Sci., 77, 264–272,, 2019. 

Lohila, A., Aurela, M., Hatakka, J., Pihlatie, M., Minkkinen, K., Penttilä, T., and Laurila, T.: Responses of N2O fluxes to temperature, water table and N deposition in a northern boreal fen, Eur. J. Soil Sci., 61, 651–661,, 2010. 

Long, K. D., Flanagan, L. B., and Cai, T.: Diurnal and seasonal variation in methane emissions in a northern Canadian peatland measured by eddy covariance, Glob. Change Biol., 16, 2420–2435,, 2010. 

Manies, K. L., Fuller, C. C., Jones, M. C., Waldrop, M. P., and McGeehin, J. P.: Soil data for a thermokarst bog and the surrounding permafrost plateau forest, located at Bonanza Creek Long Term Ecological Research Site, Interior Alaska, U.S. Geological Survey Open-File Report 2016–1173​​​​​​​, 11 p.​​​​​​​,, 2017. 

Meng, L., Hess, P. G. M., Mahowald, N. M., Yavitt, J. B., Riley, W. J., Subin, Z. M., Lawrence, D. M., Swenson, S. C., Jauhiainen, J., and Fuka, D. R.: Sensitivity of wetland methane emissions to model assumptions: application and model testing against site observations, Biogeosciences, 9, 2793–2819,, 2012. 

Milecka, K., Kowalewski, G., Fiałkiewicz-Kozieł, B., Gałka, M., Lamentowicz, M., Chojnicki, B. H., Goslar, T., and Barabach, J.: Hydrological changes in the Rzecin peatland (Puszcza Notecka, Poland) induced by anthropogenic factors: Implications for mire development and carbon sequestration, Holocene, 27, 651–664,, 2016. 

Morel, X., Decharme, B., Delire, C., Krinner, G., Lund, M., Hansen, B. U., and Mastepanov, M.: A New Process-Based Soil Methane Scheme: Evaluation Over Arctic Field Sites With the ISBA Land Surface Model, J. Adv. Model. Earth Syst., 11, 293–326,, 2019. 

ORNL DAAC: FLUXNET, Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC) for Biogeochemical Dynamics [data set],, last access: 25 March​​​​​​​ 2022. 

Nielsen, C. S., Michelsen, A., Ambus, P., Chamindu Deepagoda, T. K. K.​​​​​​​, and Elberling, B.: Linking rhizospheric CH4 oxidation and net CH4 emissions in an arctic wetland based on 13CH4 labeling of mesocosms, Plant Soil, 412, 201–213,, 2017. 

Parton, W. J., Stewart, J. W. B., and Cole, C. V.: Dynamics of C, N, P and S in grassland soils: a model, Biogeochemistry, 5, 109–131,, 1988. 

Peltola, O., Vesala, T., Gao, Y., Räty, O., Alekseychik, P., Aurela, M., Chojnicki, B., Desai, A. R., Dolman, A. J., Euskirchen, E. S., Friborg, T., Göckede, M., Helbig, M., Humphreys, E., Jackson, R. B., Jocher, G., Joos, F., Klatt, J., Knox, S. H., Kowalska, N., Kutzbach, L., Lienert, S., Lohila, A., Mammarella, I., Nadeau, D. F., Nilsson, M. B., Oechel, W. C., Peichl, M., Pypker, T., Quinton, W., Rinne, J., Sachs, T., Samson, M., Schmid, H. P., Sonnentag, O., Wille, C., Zona, D., and Aalto, T.: Monthly gridded data product of northern wetland methane emissions based on upscaling eddy covariance observations, Earth Syst. Sci. Data, 11, 1263–1289,, 2019. 

Potter, C. S.: An ecosystem simulation model for methane production and emission from wetlands, Global Biogeochem. Cy., 11, 495–506,, 1997. 

Potter, C. S., Davidson, E. A., and Verchot, L. V.: Estimation of global biogeochemical controls and seasonality in soil methane consumption, Chemosphere, 32, 2219–2246,, 1996. 

Prather, M. J., Derwent, R., Ehhalt, D., Fraser, P. J., Sanhueza, E., and Zhou, X.: Other trace gases and atmospheric chemistry, in: Climate change 1994: radiative forcing of climate change an evaluation of the IPCC IS92 emission scenarios, edited by: Houghton, J. T., Cambridge University Press, Cambridge, 73–126, (last access: 25 March 2022), 1995. 

Prather, M. J., Holmes, C. D., and Hsu, J.: Reactive greenhouse gas scenarios: Systematic exploration of uncertainties and the role of atmospheric chemistry, Geophys. Res. Lett., 39, L09803,, 2012. 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P.: Numerical Recipes in Fortran 77: the art of scientific computing, in: Numerical Recipies in Fortran 77, 2nd edn., Cambridge University Press, ISBN 9780521430647, 1993. 

Prinn, R. G., Weiss, R. F., Arduini, J., Arnold, T., DeWitt, H. L., Fraser, P. J., Ganesan, A. L., Gasore, J., Harth, C. M., Hermansen, O., Kim, J., Krummel, P. B., Li, S., Loh, Z. M., Lunder, C. R., Maione, M., Manning, A. J., Miller, B. R., Mitrevski, B., Mühle, J., O'Doherty, S., Park, S., Reimann, S., Rigby, M., Saito, T., Salameh, P. K., Schmidt, R., Simmonds, P. G., Steele, L. P., Vollmer, M. K., Wang, R. H., Yao, B., Yokouchi, Y., Young, D., and Zhou, L.: History of chemically and radiatively important atmospheric gases from the Advanced Global Atmospheric Gases Experiment (AGAGE), Earth Syst. Sci. Data, 10, 985–1018,, 2018. 

Qiu, C.: ORCHIDEE_PEAT_V2 revision 5488, Forge IPSL Jussieu [code],, 2019. 

Qiu, C., Zhu, D., Ciais, P., Guenet, B., Krinner, G., Peng, S., Aurela, M., Bernhofer, C., Brümmer, C., Bret-Harte, S., Chu, H., Chen, J., Desai, A. R., Dušek, J., Euskirchen, E. S., Fortuniak, K., Flanagan, L. B., Friborg, T., Grygoruk, M., Gogo, S., Grünwald, T., Hansen, B. U., Holl, D., Humphreys, E., Hurkuck, M., Kiely, G., Klatt, J., Kutzbach, L., Largeron, C., Laggoun-Défarge, F., Lund, M., Lafleur, P. M., Li, X., Mammarella, I., Merbold, L., Nilsson, M. B., Olejnik, J., Ottosson-Löfvenius, M., Oechel, W., Parmentier, F.-J. W., Peichl, M., Pirk, N., Peltola, O., Pawlak, W., Rasse, D., Rinne, J., Shaver, G., Schmid, H. P., Sottocornola, M., Steinbrecher, R., Sachs, T., Urbaniak, M., Zona, D., and Ziemblinska, K.: ORCHIDEE-PEAT (revision 4596), a model for northern peatland CO2, water, and energy fluxes on daily to annual scales, Geosci. Model Dev., 11, 497–519,, 2018. 

Qiu, C., Zhu, D., Ciais, P., Guenet, B., Peng, S., Krinner, G., Tootchi, A., Ducharne, A., and Hastie, A.: Modelling northern peatland area and carbon dynamics since the Holocene with the ORCHIDEE-PEAT land surface model (SVN r5488), Geosci. Model Dev., 12, 2961–2982,, 2019. 

Raivonen, M., Smolander, S., Backman, L., Susiluoto, J., Aalto, T., Markkanen, T., Mäkelä, J., Rinne, J., Peltola, O., Aurela, M., Lohila, A., Tomasic, M., Li, X., Larmola, T., Juutinen, S., Tuittila, E.-S., Heimann, M., Sevanto, S., Kleinen, T., Brovkin, V., and Vesala, T.: HIMMELI v1.0: HelsinkI Model of MEthane buiLd-up and emIssion for peatlands, Geosci. Model Dev., 10, 4665–4691,, 2017. 

Raynaud, D., Blunier, T., Ono, Y., and Delmas, R. J.: The Late Quaternary History of Atmospheric Trace Gases and Aerosols: Interactions Between Climate and Biogeochemical Cycles, in: Paleoclimate, Global Change and the Future. Global Change – The IGBP Series, edited by: Alverson, K. D., Pedersen, T. F., and Bradley, R. S., Springer, Berlin, Heidelberg,, 2003. 

Riley, W. J., Subin, Z. M., Lawrence, D. M., Swenson, S. C., Torn, M. S., Meng, L., Mahowald, N. M., and Hess, P.: Barriers to predicting changes in global terrestrial methane fluxes: analyses using CLM4Me, a methane biogeochemistry model integrated in CESM, Biogeosciences, 8, 1925–1953,, 2011. 

Salmon, E.: ORCHIDEE_MICT-PEAT-CH4 r 7020, IPSL Data Catalog [code],, 2021. 

Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Brovkin, V., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Curry, C., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., McDonald, K. C., Marshall, J., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prigent, C., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Steele, P., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., van Weele, M., van der Werf, G. R., Weiss, R., Wiedinmyer, C., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: The global methane budget 2000–2012, Earth Syst. Sci. Data, 8, 697–751,, 2016. 

Saunois, M., Stavert, A. R., Poulter, B., Bousquet, P., Canadell, J. G., Jackson, R. B., Raymond, P. A., Dlugokencky, E. J., Houweling, S., Patra, P. K., Ciais, P., Arora, V. K., Bastviken, D., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Carlson, K. M., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., Crill, P. M., Covey, K., Curry, C. L., Etiope, G., Frankenberg, C., Gedney, N., Hegglin, M. I., Höglund-Isaksson, L., Hugelius, G., Ishizawa, M., Ito, A., Janssens-Maenhout, G., Jensen, K. M., Joos, F., Kleinen, T., Krummel, P. B., Langenfelds, R. L., Laruelle, G. G., Liu, L., Machida, T., Maksyutov, S., McDonald, K. C., McNorton, J., Miller, P. A., Melton, J. R., Morino, I., Müller, J., Murguia-Flores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., Parker, R. J., Peng, C., Peng, S., Peters, G. P., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., Riley, W. J., Rosentreter, J. A., Segers, A., Simpson, I. J., Shi, H., Smith, S. J., Steele, L. P., Thornton, B. F., Tian, H., Tohjima, Y., Tubiello, F. N., Tsuruta, A., Viovy, N., Voulgarakis, A., Weber, T. S., van Weele, M., van der Werf, G. R., Weiss, R. F., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The Global Methane Budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623,, 2020. 

Segers, R.: Methane production and methane consumption: A review of processes underlying wetland methane fluxes, Biogeochemistry, 41, 23–51,, 1998. 

Smith, K. A., Ball, T., Conen, F., Dobbie, K. E., Massheder, J., and Rey, A.: Exchange of greenhouse gases between soil and atmosphere: Interactions of soil physical factors and biological processes, Eur. J. Soil Sci., 69, 10–20,, 2003. 

Smith, R. L., Howes, B. L., and Garabedian, S. P.: In situ measurement of methane oxidation in groundwater by using natural-gradient tracer tests, Appl. Environ. Microbiol., 57, 1997–2004,, 1991. 

Spahni, R., Wania, R., Neef, L., van Weele, M., Pison, I., Bousquet, P., Frankenberg, C., Foster, P. N., Joos, F., Prentice, I. C., and van Velthoven, P.: Constraining global methane emissions and uptake by ecosystems, Biogeosciences, 8, 1643–1665,, 2011. 

Strack, M., Waddington, J. M., Turetsky, M., Roulet, N. T., and Byrne, K. A.: Northern peatlands, greenhouse gas exchange and climate change, in: Peatlands and Climate Change, edited by: Strack, M., International Peat Society, 44–69, ISBN 978-952-99401-1-0, 2008. 

Sulman, B. N., Desai, A. R., Cook, B. D., Saliendra, N., and Mackay, D. S.: Contrasting carbon dioxide fluxes between a drying shrub wetland in Northern Wisconsin, USA, and nearby forests, Biogeosciences, 6, 1115–1126,, 2009. 

Tang, J. Y., Riley, W. J., Koven, C. D., and Subin, Z. M.: CLM4-BeTR, a generic biogeochemical transport and reaction module for CLM4: model development, evaluation, and application, Geosci. Model Dev., 6, 127–140,, 2013. 

Tang, J. Y. and Riley, W. J.: Technical Note: Simple formulations and solutions of the dual-phase diffusive transport for biogeochemical modeling, Biogeosciences, 11, 3721–3728,, 2014. 

Tans, P. P.: Oxygen isotopic equilibrium between carbon dioxide and water in soils, Tellus B, 50, 163–178,, 1998. 

Thornton, J. A., Kercher, J. P., Riedel, T. P., Wagner, N. L., Cozic, J., Holloway, J. S., Dubì, W. P., Wolfe, G. M., Quinn, P. K., Middlebrook, A. M., Alexander, B., and Brown, S. S.: A large atomic chlorine source inferred from mid-continental reactive nitrogen chemistry, Nature, 464, 271–274,, 2010. 

Tokida, T., Mizoguchi, M., Miyazaki, T., Kagemoto, A., Nagata, O., and Hatano, R.: Episodic release of methane bubbles from peatland during spring thaw, Chemosphere, 70, 165–171,, 2007a. 

Tokida, T., Miyazaki, T., Mizoguchi, M., Nagata, O., Takakai, F., Kagemoto, A., and Hatano, R.: Falling atmospheric pressure as a trigger for methane ebullition from peatland, Global Biogeochem. Cy., 21, GB2003,, 2007b. 

Treat, C. C., Bloom, A. A., and Marushchak, M. E.: Nongrowing season methane emissions – a significant component of annual emissions across northern ecosystems​​​​​​​, Glob. Change Biol., 24, 3331–3343,, 2018. 

van Huissteden, J., Petrescu, A. M. R., Hendriks, D. M. D., and Rebel, K. T.: Sensitivity analysis of a wetland methane emission model based on temperate and arctic wetland sites, Biogeosciences, 6, 3035–3051,, 2009. 

Vybornova, O.: Effect of re-wetting on greenhouse gas emissions from different microtopes in a cut-over bog in Northern Germany, PhD thesis, Hamburg University, (last access: 25 March 2022), 2017. 

Walter, B. P., Heimann, M., Shannon, R. D., and White, J. R.: A process-based model to derive methane emissions from natural wetlands, Geophys. Res. Lett., 23, 3731–3734,, 1996. 

Walter, P. and Heimann, M.: A process-based, climate-sensitive model to derive methane emissions from natural wetlands: Application to five wetland sites, sensitivity to model parameters, and climate, New York, Global Biogeochem. Cy., 14, 745–765,, 2000. 

Wang, T., Ottlé, C., Boone, A., Ciais, P., Brun, E., Morin, S., Krinner, G., Piao, S., and Peng, S.: Evaluation of an improved intermediate complexity snow scheme in the ORCHIDEE land surface model, J. Geophys. Res.-Atmos., 118, 6064–6079,, 2013. 

Wania, R.: Modelling northern peatland land surface processes, vegetation dynamics and methane emissions, PhD thesis, Bristol University,​​​​​​​, (last access: 25 March 2022), 2007. 

Wania, R., Ross, I., and Prentice, I. C.: Integrating peatlands and permafrost into a dynamic global vegetation model: 1. Evaluation and sensitivity of physical land surface processes, Global Biogeochem. Cy., 23, GB3014,, 2009. 

Wania, R., Ross, I., and Prentice, I. C.: Implementation and evaluation of a new methane model within a dynamic global vegetation model: LPJ-WHyMe v1.3.1, Geosci. Model Dev., 3, 565–584,, 2010. 

Whalen, S. C.: Biogeochemistry of methane exchange between natural wetlands and the atmosphere, Environ. Eng. Sci., 22, 73–94,, 2005. 

Wiesenburg, D. A. and Guinasso, N. L.: Equilibrium Solubilities of Methane, Carbon Monoxide, and Hydrogen in Water and Sea Water, J. Chem. Eng. Data, 24, 356–360,, 1979. 

Xu, J., Morris, P. J., Liu, J., and Holden, J.: PEATMAP: Refining estimates of global peatland distribution based on a meta-analysis, Catena, 160, 134–140,, 2018. 

Xu, X., Yuan, F., Hanson, P. J., Wullschleger, S. D., Thornton, P. E., Riley, W. J., Song, X., Graham, D. E., Song, C., and Tian, H.: Reviews and syntheses: Four decades of modeling methane cycling in terrestrial ecosystems, Biogeosciences, 13, 3735–3755,, 2016.  

Zak, D., Gelbrecht, J., Wagner, C., and Steinberg, C. E. W.: Evaluation of phosphorus mobilization potential in rewetted fens by an improved sequential chemical extraction procedure, Eur. J. Soil Sci., 59, 1191–1201,, 2008. 

Zhang, Z., Zimmermann, N. E., Kaplan, J. O., and Poulter, B.: Modeling spatiotemporal dynamics of global wetlands: comprehensive evaluation of a new sub-grid TOPMODEL parameterization and uncertainties, Biogeosciences, 13, 1387–1408,, 2016. 

Zhu, Q., Liu, J., Peng, C., Chen, H., Fang, X., Jiang, H., Yang, G., Zhu, D., Wang, W., and Zhou, X.: Modelling methane emissions from natural wetlands by development and application of the TRIPLEX-GHG model, Geosci. Model Dev., 7, 981–999,, 2014. 

Zhu, X., Zhuang, Q., Qin, Z., Glagolev, M., and Song, L.: Estimating wetland methane emissions from the northern high latitudes from 1990 to 2009 using artificial neural networks, Global Biogeochem. Cy., 27, 592–604,, 2013. 

Zhuang, Q., Melillo, J. M., Kicklighter, D. W., Prinn, R. G., McGuire, A. D., Steudler, P. A., Felzer, B. S., and Hu, S.: Methane fluxes between terrestrial ecosystems and the atmosphere at northern high latitudes during the past century: A retrospective analysis with a process-based biogeochemistry model, Global Biogeochem. Cy., 18, 18, GB3010,, 2004. 

Short summary
A methane model that features methane production and transport by plants, the ebullition process and diffusion in soil, oxidation to CO2, and CH4 fluxes to the atmosphere has been embedded in the ORCHIDEE-PEAT land surface model, which includes an explicit representation of northern peatlands. This model, ORCHIDEE-PCH4, was calibrated and evaluated on 14 peatland sites. Results show that the model is sensitive to temperature and substrate availability over the top 75 cm of soil depth.