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, marieke.scheel@gmail.com

## 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)