ORCHIDEE-PEAT (revision 4596), a model for northern peatland
CO2, water and energy fluxes on daily to annual scales
Authors: Qiu Chunjing , Zhu Dan , Ciais Philippe , Guenet Bertrand , Krinner Gerhard , Peng Shushi , Aurela Mika , Bernhofer Christian , Brümmer Christian , Bret-Harte Syndonia , Chu Housen , Chen Jiquan , Desai Ankur R , Dušek Jiří , Euskirchen Eugénie S. , Fortuniak Krzysztof , Flanagan Lawrence B. , Friborg Thomas , Grygoruk Mateusz , Gogo Sébastien , Grünwald Thomas , Hansen Birger U. , Holl David , Humphreys Elyn , Hurkuck Miriam , Kiely Gerard , Klatt Janina , Kutzbach Lars , Largeron Chloé , Laggoun-Défarge Fatima , Lund Magnus , Lafleur Peter M. , Li Xuefei , Mammarella Ivan , Merbold Lutz , Nilsson Mats B. , Olejnik Janusz , Ottosson-Löfvenius Mikaell , Oechel Walter , Parmentier Frans-Jan W. , Peichl Matthias , Pirk Norbert , Peltola Olli , Pawlak Włodzimierz , Rebmann Corinna , Rasse Daniel , Rinne Janne , Shaver Gaius , Schmid Hans Peter , Sottocornola Matteo , Steinbrecher Rainer , Sachs Torsten , Urbaniak Marek , Zona Donatella , Ziemblinska Klaudia 
Published: September 2017

Abstract. Peatlands store substantial amount of carbon, are vulnerable to climate change. To predict the fate of carbon stored in peatlands, the complex interactions between water, peat and vegetations need more attention. This study describes a modified version of the ORCHIDEE land surface model for simulating the hydrology, surface energy and CO2 fluxes of peatlands on daily to annual time scales. The model, referred to as ORCHIDEE-PEAT, includes a separate soil tile in each 0.5° grid-cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation with a grid-cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model is evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site to match the peak of growing season gross primary productivity (GPP), derived from direct EC measurements. Regarding short-term variations from day to day, the model performance was good for the variations in GPP (r2 = 0.76, Nash-Sutcliff modeling efficiency, MEF = 0.76), with lesser accuracy for latent heat fluxes (LE, r2 = 0.42, MEF = 0.14) and Net ecosystem CO2 exchange (NEE, r2 = 0.38, MEF = 0.26). Seasonal variations in GPP, NEE and energy fluxes on monthly scales showed moderate to high r2 values ranging from 0.57 to 0.86. For spatial across-sites gradients of annual mean GPP, NEE and LE, r2 of 0.93, 0.27, and 0.71, respectively, were achieved. The water table variations are not well predicted (r2 < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, when using the observed water table in the carbon module to define the fraction of oxic and anoxic decomposition instead of the modeled water table, ORCHIDEE-PEAT shows a small improvement in reproducing NEE. Moreover, we found a significant relationship between optimized Vcmax and the latitude (temperature), which can better reflect the spatial gradients of annual NEE than using an average Vcmax value. In a future version of ORCHIDEE-PEAT, the influences of water table on photosynthesis and depth-dependent influences of soil temperature on respiration may be included.