Fractional cover mapping of spruce and pine at 1 ha resolution combining very high and medium spatial resolution satellite imagery

Increases in extreme weather events associated with climate change have the potential to put currently healthy forests at risk. One option to minimize this risk is the application of forest management measures aimed at generating species mixtures predicted to be more resilient to these threats. In order to apply such measures appropriately, forest managers need up-to-date, accurate and consistent forest maps at relatively fine spatial resolutions. Cost efficiency is a major factor when creating such maps. Taking European spruce (Picea abies) and Scots pine (Pinus sylvestris) as an example, this paper describes an innovative approach for mapping two tree species using a combination of commercial very high resolution WorldView-2 (WV2) images and Landsat time series data. As a first step, this study used a supervised object-based classification of WV2 images covering relatively small test sites distributed across the region of interest. Using these classification maps as training data, wall-to-wall mapping of fractional coverages of spruce and pine was achieved using multi-temporal Landsat data and Random Forests (RF) regression. The method was applied for the entire state of Bavaria (Germany), which comprises a total forested area of approximately 26,000km². As applied here, this two-step approach yields consistent and accurate maps of fractional tree cover estimates with a spatial resolution of 1ha. Independent validation of the fractional cover estimates using 3780 reference samples collected through visual interpretation of orthophotos produced root-mean-square errors (RMSE) of 11% (for spruce) and 14% (for pine) with almost no bias, and R² values of 0.74 and 0.79 for spruce and pine, respectively. The majority of the validation samples (75% (spruce) and 84% (pine)) were modeled within the assumed uncertainty of ±15% of the reference sample. Accuracies were significantly better compared to those achieved using a single-step classification of Landsat time series data at the pixel level (30m), because the two-step approach better captures regional variation in the spectral signatures of target classes. Moreover, the increased number of available reference cells mitigates the impact of occasional errors in the reference data set. This two-step approach has great potential for cost-effective operational mapping of dominant forest types over large areas.