Increased Central European forest mortality explained by higher harvest rates driven by enhanced productivity

# Increased Central European forest mortality explained by higher harvest rates driven by enhanced productivity ## Author Marieke Scheel, Lund University, Sweden, [email protected] ## Description Data underlying analysis in: Marieke Scheel, Mats Lindeskog, Benjamin Smith, Susanne Suvanto, Thomas A. M. Pugh Increased Central European forest mortality explained by higher harvest rates driven by enhanced productivity Scripts underlying analysis: https://github.com/mariekesche/harvest_driven_canopy_mortality Folder: harvest_checks - NFI_Germany.txt, National Forest Index data from Germany as difference between inventories in 200-2003 and 2011-2013, values are given as fraction, n: number of NFI plots in a grid cell, HARVEST_ALL: clear-cut harvest, HARVEST_PARTIAL: thinning harvest, NATDEAD: natural dead (not harvested) Folder: manag_climfix (S_man,clim) - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4) - crownloss.txt, m^2 of crown/m^2 of ground lost per year - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: manag_co2fix (S_man,CO2) - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4) - crownloss.txt, m^2 of crown/m^2 of ground lost per year - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: manag_ndepfix (S_man,N) - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4) - crownloss.txt, m^2 of crown/m^2 of ground lost per year - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: manag_nofix (S_man) - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4) - closs.txt, biomass loss in kg [C]/m^2 year split into DBH classes - closs_harv.txt, biomass loss due to harvest in kg [C]/m^2 year split into DBH classes - cpool_forest_rel.txt, biomass of forests in kg [C]/m^2 year - crownloss.txt, m^2 of crown/m^2 of ground lost per year - crownloss_age.txt, m^2 of crown/m^2 of ground lost per year due to age - crownloss_dist.txt, m^2 of crown/m^2 of ground lost per year due to disturbance - crownloss_fire.txt, m^2 of crown/m^2 of ground lost per year due to fire disturbance - crownloss_greff.txt, m^2 of crown/m^2 of ground lost per year due to growth efficiency - crownloss_harv.txt, m^2 of crown/m^2 of ground lost per year due to harvest - crownloss_other.txt, m^2 of crown/m^2 of ground lost per year due to other reasons - crownloss_thin.txt, m^2 of crown/m^2 of ground lost per year due to natural thinning - diam_cmass_wood.txt, wooden biomass in kg [C]/m^2 year split into DBH classes - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes - diam_dens.txt, total number of trees/m^2 year split into DBH classes - stemloss.txt, number of trees/m^2 year lost in respective cells - stemloss_harv.txt, number of trees/m^2 year lost in respective cells due to harvest Folder: manag_nothin (S_nothin) - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4) - crownloss.txt, m^2 of crown/m^2 of ground lost per year - diam_cmass_wood.txt, wooden biomass in kg [C]/m^2 year split into DBH classes - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year - lai.txt, leaf area index (LAI) for simulated species and plant functional types (PFTs) Folder: PNV (S_PNV) - cflux.txt, Net Primary Production (NPP) values in kg [C]/ m^2 year (column 4) - crownloss.txt, m^2 of crown/m^2 of ground lost per year - crownloss_age.txt, m^2 of crown/m^2 of ground lost per year due to age - crownloss_dist.txt, m^2 of crown/m^2 of ground lost per year due to disturbance - crownloss_fire.txt, m^2 of crown/m^2 of ground lost per year due to fire disturbance - crownloss_greff.txt, m^2 of crown/m^2 of ground lost per year due to growth efficiency - crownloss_other.txt, m^2 of crown/m^2 of ground lost per year due to other reasons - crownloss_thin.txt, m^2 of crown/m^2 of ground lost per year due to natural thinning - diam_cmass_wood.txt, wooden biomass in kg [C]/m^2 year split into DBH classes - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: dependencies - gridlist.txt, coordinates of 0.5°x0.5° grid cells that simulations were run on; longitude, latitude, FAO number, size in m^2 - landcover_eu.txt, input file LPJ-GUESS model showing changes from natural to forest (harvest); longitude, latitude, year, natural, forest, barren ## Dependencies - canopy mortality rates published in "Senf C Pflugmacher D Zhiqiang Y Sebald J Knorn J Neumann M Hostert P and Seidl R 2018 Canopy mortality has doubled in Europe’s temperate forests over the last three decades Nature Communications 9 4978 10.1038/s41467-018-07539-6" - harvest removal rates published in "Ceccherini G Duveiller G Grassi G Lemoine G Avitabile V Pilli R and Cescatti A 2020 Abrupt increase in harvested forest area over Europe after 2015 Nature 583 72-77 10.1038/s41586-020-2438-y" - FAO forest removal area data data retrieved from https://www.fao.org/faostat/en/#data/GF (24.07.2021), modified for overview (1st column: Area Code, 2nd column: Year, 3rd column: Area in 1000 ha)