Europe's forest harvesting regimes dataset

The data contains quantification of forest harvesting regimes in 11 European countries based on permanent plots of forest inventories. The results are aggregated on a resolution of 1 degree latitude/longitude. The results and the methods for obtaining them are reported in detail in the publication below:

Suvanto, S., Esquivel-Muelbert, A., Schelhaas, M.J., Astigarraga, J., Astrup, R., Cienciala, E., Fridman, J., Henttonen, H.M., Kunstler, G., Kändler, G., König, L.A., Ruiz-Benito, P., Senf, C., Stadelmann, G., Starcevic, A., Talarczyk, A., Zavala, M.A., Pugh, T.A.M. (2025). Understanding Europe’s Forest Harvesting Regimes. Earth's Future 13 (2), e2024EF005225, https://doi.org/10.1029/2024EF005225.

The file variables_harvest_regimes_zenodo1.csv contains the results for each 1-degree grid cell. The columns in the file are described below (please see more details in the paper and it's supplementary material):

Colum name
Description
Unit
In Suvanto et al.

lat
Latitude of the grid cell center point
degree

lon
Longitude of the grid cell center point
degree

country
Country (in case of grid cells covering area from multiple countries, each country has a separate record in the data)

n_plots
Number of forest inventory plots in the grid cell
count
Supplement, Fig. S3

n_harvest
Number of forest inventory plots with harvest in the grid cell
count
Supplement, Fig. S3

harvest_rate
Total harvest rate, calculated as the percentage of tree basal area harvested in the grid cell per year.
percentage (%)
Fig. 2, Fig. 3

harvest_freq
Frequency of harvest events, calculated as the percentage of plots harvested per year.
percentage (%)
Fig. 2, Fig. 3

harvest_intensity_mean
Mean intensity of harvest events, calculated as the average percentage of tree basal area removed in a harvest event.
percentage (%)
Fig. 2, Fig. 3

harvest_intensity_sd
Standard deviation of the intensity of harvest events, calculated as the standard deviation of the percentage of tree basal area removed in a harvest event.
percentage (%)
Supplement, Fig. S3

p_harvests_in_0_25
Fraction of harvest events with intensity of 0-25%
fraction (between 0-1)
Fig. 4

p_harvests_in_25_50
Fraction of harvest events with intensity of 25-50%
fraction (0-1)
Fig. 4

p_harvests_in_50_75
Fraction of harvest events with intensity of 50-75%
fraction (0-1)
Fig. 4

p_harvests_in_75_100
Fraction of harvest events with intensity of 75-100%
fraction (0-1)
Fig. 4