Search

Biodiversity monitoring represents a major challenge to supporting proper forest ecosystem management and biodiversity conservation. The latter is indeed shifting in recent years from single-species to multi-taxon approaches. However, multi-taxonomic studies are quite rare due to the effort required for performing field surveys. In this context, remote sensing is a powerful tool, continuously providing consistent and open access data at a different range of spatial and temporal scales. In particular, the Sentinel-2 (S2) mission has great potential to produce reliable proxies for biological diversity. In beech forests of two Italian National Parks, we sampled the beetle fauna, breeding birds, and epiphytic lichens. First, we calculated Shannon’s entropy and Simpson’s diversity. Then, to produce variables for biodiversity assessment, we exploited S2 data acquired in the 4 years 2017–2021. S2 images were used to construct spectral bands and photosynthetic indices time series, from which 91 harmonic metrics were derived. For each taxon and multi-taxon community, we assessed the correlation with S2 harmonic metrics, biodiversity indices, and forest structural variables. Then, to assess the potential of the harmonic metrics in predicting species diversity in terms of Shannon’s and Simpson’s biodiversity indices, we also fit a random forests model between each diversity index and the best 10 harmonic metrics (in terms of absolute correlation, that is, the magnitude of the correlation) for each taxon. The models’ performance was evaluated via the relative root mean squared error (RMSE%). Overall, 241 beetle, 27 bird, and 59 lichen species were recorded. The diversity indices were higher for the multi-taxon community than for the single taxa. They were generally higher in the CVDA site than in GSML, except for the bird community. The highest correlation values between S2 data and biodiversity indices were recorded in CVDA for multi-taxon and beetle communities (| r| = 0.52 and 0.38, respectively), and in GSML for lichen and beetle communities (| r| = 0.34 and 0.26, respectively). RMSE% ranged between 2.53 and 9.99, and between 8.1 and 16.8 for the Simpson and Shannon index, respectively. The most important variables are phase and RMSE of red-Edge bands for bird and lichen communities, while RMSE and time of tassel cap and from EVI indices for beetles and multi-taxon diversity. Our results demonstrate that S2 data can be used for identifying potential biodiversity hotspots, showing that the herein presented harmonic metrics are informative for several taxa inhabiting wood, giving concrete support to cost-effective biodiversity monitoring and nature-based forest management in complex mountain systems.

Afforestation is one of the most effective processes for removing carbon dioxide from the atmosphere and combating global warming. Landsat data and machine learning approaches can be used to map afforestation (i) indirectly, by constructing two maps of the same area over different periods and then predicting changes, or (ii) directly, by constructing a single map and analyzing observations of change in both the response and remotely sensed variables. Of crucial importance, no comprehensive comparisons of direct and indirect approaches for afforestation monitoring are known to have been conducted to date. Afforestation maps estimated through the analysis of remotely sensed data may serve as intermediate products for guiding the selection of samples and the production of statistics. In this and similar studies, a huge effort is dedicated to collecting validation data. In turn, those validation datasets have varying sampling intensities in different areas, which complicates their use for assessing the accuracies of new maps. As a result, the work done to collect data is often not sufficiently exploited, with some validation datasets being used just once. In this study, we addressed two main aims. First, we implemented a methodology to reuse validation data acquired via stratified sampling with strata constructed from remote sensing maps. Second, we used this method for acquiring data for comparing map accuracy estimates and the precision of estimates for direct and indirect approaches for country-wide mapping of afforestation that occurred in Italy between 1985 and 2019. To facilitate these comparisons, we used Landsat imagery, random forest classification, and Google Earth Engine. The herein-presented method produced different accuracy estimates with 95% confidence interval and for different map classes. Afforestation accuracies ranged between 53 ± 5.9% for the indirect map class inside the buffer—defined as a stratum within 120 m of the forest/non-forest mask boundaries—and 26 ± 3.4% for the direct map outside the buffer. The accuracy in non-afforestation map classes was much greater, ranging from 87 ± 1.9% for the indirect map inside the buffer to 99 ± 1.3% for the direct map outside the buffer. Additionally, overall accuracies (with 95% CI) were estimated with large precision for both direct and indirect maps (87 ± 1.3% and 89 ± 1.6%, respectively), confirming (i) the effectiveness of the method we introduced for reusing samples and (ii) the relevance of remotely sensed data and machine learning for monitoring afforestation.