Statistical mapping of tree species over Europe

In order to map the spatial distribution of twenty tree species groups over Europe at 1 km x 1 km resolution, the ICP-Forest Level-I plot data were extended with the National Forest Inventory (NFI) plot data of eighteen countries. The NFI grids have a much smaller spacing than the ICP grid. In areas with NFI plot data, the proportions of the land area covered by the tree species were mapped by compositional kriging. Outside these areas, these proportions were mapped with a multinomial multiple logistic regression model. A soil map, a biogeographical map and bioindicators derived from temperature and precipitation data were used as predictors. Both methods ensure that the predicted proportions are in the interval [0,1] and sum to 1. The regression predictions were iteratively scaled to the National Forest Inventory statistics and the Forest map of Europe. The predicted proportions for the twenty tree species were validated by the Bhattacharryya distance between predicted and observed proportions at 230 plot data separated from the calibration data. Besides, the map with the predicted dominant species was validated by computing the error matrix. The median Bhattacharryya distance in the subarea with NFI plot data was 1.712, whereas in the subarea with ICP-Level-I data, this was 2.131. The scaling did not significantly decrease the Bhattacharryya distance. The estimated overall accuracy of this map was 43%. In areas with NFI plot data, overall accuracy was 57%, outside these areas 33%. This gain was mainly attributable to the much denser plot data, less to the prediction method.