DOI: 10.5281/zenodo.6393789
title: Increased Central European forest mortality explained by higher harvest rates driven by enhanced productivity

description: # Increased Central European forest mortality explained by higher harvest rates driven by enhanced productivity ## Author<br> Marieke Scheel, Lund University, Sweden, marieke.scheel@gmail.com ## Description<br> Data underlying analysis in:<br> Marieke Scheel, Mats Lindeskog, Benjamin Smith, Susanne Suvanto, Thomas A. M. Pugh<br> Increased Central European forest mortality explained by higher harvest rates driven by enhanced productivity<br> Scripts underlying analysis:<br> https://github.com/mariekesche/harvest_driven_canopy_mortality Folder: harvest_checks<br> - 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)<br> - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4)<br> - crownloss.txt, m^2 of crown/m^2 of ground lost per year<br> - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: manag_co2fix (S_man,CO2)<br> - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4)<br> - crownloss.txt, m^2 of crown/m^2 of ground lost per year<br> - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: manag_ndepfix (S_man,N)<br> - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4)<br> - crownloss.txt, m^2 of crown/m^2 of ground lost per year<br> - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: manag_nofix (S_man)<br> - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4)<br> - closs.txt, biomass loss in kg [C]/m^2 year split into DBH classes<br> - closs_harv.txt, biomass loss due to harvest in kg [C]/m^2 year split into DBH classes<br> - cpool_forest_rel.txt, biomass of forests in kg [C]/m^2 year<br> - crownloss.txt, m^2 of crown/m^2 of ground lost per year<br> - crownloss_age.txt, m^2 of crown/m^2 of ground lost per year due to age<br> - crownloss_dist.txt, m^2 of crown/m^2 of ground lost per year due to disturbance<br> - crownloss_fire.txt, m^2 of crown/m^2 of ground lost per year due to fire disturbance<br> - crownloss_greff.txt, m^2 of crown/m^2 of ground lost per year due to growth efficiency<br> - crownloss_harv.txt, m^2 of crown/m^2 of ground lost per year due to harvest<br> - crownloss_other.txt, m^2 of crown/m^2 of ground lost per year due to other reasons<br> - crownloss_thin.txt, m^2 of crown/m^2 of ground lost per year due to natural thinning<br> - diam_cmass_wood.txt, wooden biomass in kg [C]/m^2 year split into DBH classes<br> - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes<br> - diam_dens.txt, total number of trees/m^2 year split into DBH classes<br> - stemloss.txt, number of trees/m^2 year lost in respective cells<br> - stemloss_harv.txt, number of trees/m^2 year lost in respective cells due to harvest Folder: manag_nothin (S_nothin)<br> - cflux_forest_rel.txt, Net Primary Production (NPP) values of forests in kg [C]/ m^2 year (column 4)<br> - crownloss.txt, m^2 of crown/m^2 of ground lost per year<br> - diam_cmass_wood.txt, wooden biomass in kg [C]/m^2 year split into DBH classes<br> - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year<br> - lai.txt, leaf area index (LAI) for simulated species and plant functional types (PFTs) Folder: PNV (S_PNV)<br> - cflux.txt, Net Primary Production (NPP) values in kg [C]/ m^2 year (column 4)<br> - crownloss.txt, m^2 of crown/m^2 of ground lost per year<br> - crownloss_age.txt, m^2 of crown/m^2 of ground lost per year due to age<br> - crownloss_dist.txt, m^2 of crown/m^2 of ground lost per year due to disturbance<br> - crownloss_fire.txt, m^2 of crown/m^2 of ground lost per year due to fire disturbance<br> - crownloss_greff.txt, m^2 of crown/m^2 of ground lost per year due to growth efficiency<br> - crownloss_other.txt, m^2 of crown/m^2 of ground lost per year due to other reasons<br> - crownloss_thin.txt, m^2 of crown/m^2 of ground lost per year due to natural thinning<br> - diam_cmass_wood.txt, wooden biomass in kg [C]/m^2 year split into DBH classes<br> - diam_crownarea.txt, total m^2 of crown/m^2 of ground per year split into DBH classes Folder: dependencies<br> - gridlist.txt, coordinates of 0.5°x0.5° grid cells that simulations were run on; longitude, latitude, FAO number, size in m^2<br> - landcover_eu.txt, input file LPJ-GUESS model showing changes from natural to forest (harvest); longitude, latitude, year, natural, forest, barren ## Dependencies<br> - 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"<br> - harvest removal rates published in "Ceccherini G Duveiller G Grassi G Lemoine G Avitabile V Pilli R and Cescatti A 2020<br> Abrupt increase in harvested forest area over Europe after 2015 Nature 583 72-77 10.1038/s41586-020-2438-y"<br> - 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)