DOI: 10.5281/zenodo.8354397
title: Data and R-scripts for estimating carbon dioxide emissions from drained peatland forest soils for the greenhouse gas inventory of Finland

description: <strong> Introduction</strong> A new method for estimating carbon dioxide emissions from rained peatland forest soils was developed for the Greenhouse Gas Inventory of Finland (GHG inventory). The method is based on a set of models (Ojanen et al. 2014, Tuomi et al., 2009) that dynamically compile all relevant carbon inputs and outputs into a time series of soil CO<sub>2</sub> emission. A complete description of the method is described in Alm et al. (2023). Here we present the input data and R-scripts (R Core Team, 2020) for computing the time series from year 1990 to 2022 of CO<sub>2</sub> emission from soil in forest land on drained organic soil, like it was reported by the Finnish GHG inventory (Statistics Finland, 2023). <strong>Time series data </strong> The source of forest and area data is the Finnish National Forest Inventory (NFI) as a part of Luke Statutory Services. The NFI standing forest data in the data files includes annual country-wide estimates of mean basal area and standing biomass of Scots pine (<em>Pinus sylvestris</em> L.), Norway spruce (Picea abies (L.) H. Karst) and all the broadleaved forest trees combined. The data concerns forest land on drained organic soil only (class FRA 1 according to the FAO forest land definition). The NFI data for each year has been averaged by different drained peatland forest site types (FTYPE) and by inventory regions of southern and northern Finland. The areas and proportions of FTYPEs of all drained peatland “forests remaining forests” (i.e., forests that have not undergone another change in land use in the past 20 years) in southern and northern Finland (Alm et al., 2023), derived from NFI12 (2014–2018). Annual litter input from harvest residues was estimated using statistics of harvested stem volumes by species, collected and published by Luke (Luke statistics). The stem volumes were converted to whole trees and further to litter fractions and further to The share of residues remaining in forest is estimated by subtracting the amount of the logging residues collected for energy use, the data obtained from Luke statistics/energy. The biomass of live trees, annual litterfall from live trees aboveground and root litter belowground are derived from the National Forest Inventory of Finland (inventory rounds NFI8 to NFI13). The R-code also includes calculation of annual litter production from the harvesting residues. The regression-based transfer models, implemented in the R-code, also need meteorological time series inputs: The soil organic matter decomposition model (Ojanen et al. 2014) uses May-October mean temperature. Decomposition model yasso07 (Tuomi et al., 2009), applied for estimating the CO<sub>2</sub> release by decomposition of harvesting residues and above ground litter from natural mortality, is constrained by annual temperature, annual temperature amplitude and annual precipitation. Starting from the original country-wide grid produced by the Finnish Meteorological Institute (FMI) the weather time series were spatially averaged so that the FMI weather grid values were collected from those locations where peatlands representing each FTYPE in southern and northern Finland were observed by the NFI, respectively. The pre-prepared input data are given in files, see Table 1 for descriptions. Table 1. Description of input data files. <strong>File</strong> <strong>Description of data</strong> basal.areas.csv Time series of years 1990-2022 for annual average basal area (m<sup>2</sup> ha<sup>-1</sup>) by year, by peatland forest site type (peat_type) and by tree species or group (tree_type). Values of peat_type correspond to FTYPE: 1 Herb-rich type 2 <em>Vaccinium myrtillus</em> type 4 <em>Vaccinium vitis-idaea</em> type 6 Dwarf shrub type 7 <em>Cladina</em> type Values of tree species or group correspond to: 1 Scots pine 2 Norway spruce 3 Broadleaved species biomass.csv Time series of years 1990-2022 for annual biomass (biomass, t ha<sup>-1</sup> of dry mass) by year, by biomass component, by tree species and by peatland forest site type (tkg). Values of peat_type correspond to FTYPE: 1 Herb-rich type 2 <em>Vaccinium myrtillus</em> type 4 <em>Vaccinium vitis-idaea</em> type 6 Dwarf shrub type 7 <em>Cladina</em> type dead_litter.csv Time series of years 1990-2022 of annual aboveground litter from dead wood: Harvesting residues and natural mortality combined (C, t ha<sup>-1</sup> of dry mass; lognat_litter). Values of region correspond to GHG inventory region: south South Finland north North Finland ghgi_litter.csv Time series of years 1990-2022 for litter AWEN-fractions (A=acid soluble, W=water soluble, E=ethanol soluble, N=non-soluble; C, t ha<sup>-1</sup>) by different litter types: Above-ground coarse woody litter (coarse_woody_litter), fine woody litter (fine_woody_litter), non-woody litter (non_woody_litter) by litter source and deposition type by region. “org” denotes organic soil. Values of region correspond to GHG inventory region: south South Finland north North Finland Values of ground correspond to litter deposition environment: above Above-ground litter below Below-ground litter lognat_decomp.csv Time series of years 1990-2022 for C, t ha<sup>-1</sup> of dry mass, decomposed from logging residues and natural mortality by region. Values of variable “region” correspond to GHG inventory region: south South Finland north North Finland logyasso_weather_data.csv Time series of years 1990-2022 for regional (region) precipitation sum (mm, sum_P), average annual temperature (°C, mean_T) and amplitude of the annual temperature (°C , ampli_T). Values of region correspond to GHG inventory region: south South Finland north North Finland total_area.csv Areas (ha) of drained peatland forests remaining forest land by region and peat_type. Values of variable “region” correspond to GHG inventory region: south South Finland north North Finland Values of peat_type correspond to FTYPE: 1 Herb-rich type 2 <em>Vaccinium myrtillus</em> type 4 <em>Vaccinium vitis-idaea</em> type 6 Dwarf shrub type 7 <em>Cladina</em> type weather_data.csv Time series of years 1990-2022 for 30-year rolling mean temperature for the May-October period (roll_T) used by the soil decomposition models. The values are calculated for each FTYPE (peat_type) using their spatial distributions (see details in Alm et al., 2023). Values of variable “region” correspond to GHG inventory region: south South Finland north North Finland Values of peat_type correspond to FTYPE: 1 Herb-rich type 2 <em>Vaccinium myrtillus</em> type 4 <em>Vaccinium vitis-idaea</em> type 6 Dwarf shrub type 7 <em>Cladina</em> type <strong>The R-scripts</strong> The scripts are an excerpt from the Finnish greenhouse gas inventory code set, applying the necessary pre-processed input data and producing the soil CO<sub>2</sub> emissions for each FTYPE separately. The necessary R-packages (R Core Team, 2020) are managed in the script LIBRARIES.R. Guidance for running the R-scripts is given in the README.txt. <strong>References</strong> Alm, J., Wall, A., Myllykangas, J-P., Ojanen, P., Heikkinen, J., Henttonen, H. M., Laiho, R., Minkkinen, K., Tuomainen, T. and Mikola, J. A new method for estimating carbon dioxide emissions from drained peatland forest soils for the greenhouse gas inventory of Finland. Biogeosciences https://doi.org/10.5194/bg-20-1-2023, 2023. LUKE Statistics https://www.luke.fi/en/statistics/total-roundwood-removals-and-drain, last access 8.12.2022. https://www.luke.fi/en/statistics/commercial-fellings/commercial-fellings-72023. last access 8.12.2022. Statistics Finland 2023. URL: https://unfccc.int/documents/627718 (last access 13.9.2023). Ojanen, P., Lehtonen, A., Heikkinen, J., Penttilä, T., and Minkkinen, K.: Soil CO2 balance and its uncertainty in forestry drained peatlands in Finland, Forest Ecol. Manage., 325, 60–73, 2014. R Core Team: R: A language and environment for statistical computing. R Foundation forStatistical Computing, Vienna, Austria, URL https://www.R-project.org, 2020. Tuomi, M., Thum, T., Järvinen, H., Fronzek, S., Berg, B., Harmon, M., Trofymow, J.A., Sevanto, S. and Liski, J.: Leaf litter decomposition - Estimates of global variability based on Yasso07 model, Ecol. Modell. 220 (23):3362-3371, 2009.
description: Includes both input dataset (time series 1990-2022) and software (R-scripts) for producing the respective result time series of soil CO2 balance in Finland for greenhouse gas inventory. The data directory structure for outputs is also given.