Search

ABSTRACTQuantifying the carbon-stocking contribution of forest plantations is a crucial but challenging and expensive process, usually performed through field analysis. For this reason, plantations’ carbon storage is often calculated and reported using generic and inaccurate functions relying exclusively on tree species and plantation age. This study introduces a new field-independent (FI) method for forest plantations’ carbon quantification and mapping through automatic analysis of Sentinel-2 data. The study area is a Guatemalan forest plantation of 20 hectares, for which we constructed a reference dataset measuring in the field the diameter and the height of all trees within 20 randomly selected plots (10-meter radius). The CO2 equivalent absorbed by the plantation was first estimated using ground data and a design-based (DB) approach. Then, to obtain CO2 equivalent estimates but also maps, we used both ground and Sentinel-2 data to compare a standard model-assisted (MA) approach relying on Random Forests with the FI approach. Our results demonstrate that the FI method provides carbon stock statistics comparable to those obtained using DB and MA methods and more accurate maps. Accordingly, the RMSE obtained using the FI method was 34% while that obtained by the MA method – exploiting random forest algorithm – was greater (RMSE = 39%). The 95% confidence interval estimates of the CO2 stored in the plantation were 100 ± 18 MgC ha−1 and 102 ± 8 MgC ha−1, for DB and MA respectively. Using the FI method, the CO2 ranged between 89 and 117 Mg C ha−1, all values within the DB confidence interval. In addition, the FI map was surprisingly consistent with the MA-derived map, making our approach a valid alternative for monitoring plantation status and carbon storage when ground data are not available.

In European mountains most beech forest areas have been managed for timber production. This practice has reduced the availability of biomass for the whole forest-dwelling species assemblage and of deadwood for the saproxylic community. Despite most of Italy’s beech stands having a long history of management, its effects on forest species remain poorly understood. To address this gap, we studied beetle abundance and diversity in five beech-dominated forests with increasing management intensity in central Italy’s Apennines (Tuscany). We assessed if forests with similar management intensity exhibited comparable patterns in beetle diversity, abundance, and commonness versus rarity. Three forests were managed with even-aged shelterwood; one was managed with continuous cover forestry; and one was old-growth. We found 25 beetle families and 195 species across all sites with similar total abundance and richness. However, the representation of the most abundant families varied among sampling sites (ANOVA test: always significant for the total abundance of the most abundant families: F ≥ 2.77, d.f. = 4, p ≤ 0.038). The old-growth forest harbored more threatened species than managed sites. Saproxylic assemblages were similar between the recently cut site and the old-growth forest, and between shelterwood and continuous cover sites. While the similarity gradient among the whole species assemblages reflected geographical proximity, the similarity gradient among saproxylic assemblages reflected the successional proximity among forest management systems. Our research underscores the effects of management on beetle diversity, offering insights for sustainable forestry.

Carbon assimilation and wood production are influenced by environmental conditions and endogenous factors, such as species auto-ecology, age, and hierarchical position within the forest structure. Disentangling the intricate relationships between those factors is more pressing than ever due to climate change’s pressure. We employed the 3D-CMCC-FEM model to simulate undisturbed forests of different ages under four climate change (plus one no climate change) Representative Concentration Pathways (RCP) scenarios from five Earth system models. In this context, carbon stocks and increment were simulated via total carbon woody stocks and mean annual increment, which depends mainly on climate trends. We find greater differences among different age cohorts under the same scenario than among different climate scenarios under the same age class. Increasing temperature and changes in precipitation patterns led to a decline in above-ground biomass in spruce stands, especially in the older age classes. On the contrary, the results show that beech forests will maintain and even increase C-storage rates under most RCP scenarios. Scots pine forests show an intermediate behavior with a stable stock capacity over time and in different scenarios but with decreasing mean volume annual increment. These results confirm current observations worldwide that indicate a stronger climate-related decline in conifers forests than in broadleaves.

The quantification of tree-related microhabitats (TreMs) and multi-taxon biodiversity is pivotal to the implementation of forest conservation policies, which are crucial under the current climate change scenarios. We assessed the capacity of Airborne Laser Scanning (ALS) data to quantify biodiversity indices related to both forest beetle and bird communities and TreMs, calculating the species richness and types of saproxylic and epixylic TreMs using the Shannon index. As biodiversity predictors, 240 ALS-derived metrics were calculated: 214 were point-cloud based, 14 were pixel-level from the canopy height model, and 12 were RGB spectral statistics. We used the random forests algorithm to predict species richness and the Shannon diversity index, using the field plot measures as dependent variables and the ALS-derived metrics as predictors for each taxon and TreMs type. The final models were used to produce wall-to-wall maps of biodiversity indices. The Shannon index produced the best performance for each group considered, with a mean difference of −6.7%. Likewise, the highest R2 was for the Shannon index (0.17, against 0.14 for richness). Our results confirm the importance of ALS data in assessing forest biodiversity indicators that are relevant for monitoring forest habitats. The proposed method supports the quantification and monitoring of the measures needed to implement better forest stands and multi-taxon biodiversity conservation.

Forest parameter estimation is required to support the sustainable management of forest ecosystems. Currently, forest resource assessment is increasingly linked to auxiliary information obtained from remote sensing (RS) technologies. In forest parameter estimation, airborne laser scanning (ALS) data have been demonstrated to be an invaluable source of information. However, ALS data are often not available for the entire forest area, whereas images from multiple satellite systems offer new opportunities for large-scale forest surveys. This study aims to assess and estimate the contribution of field plot data and ALS data along with Landsat data to the precision of growing stock volume (GSV) estimates. We compared different approaches for model-assisted estimation of mean forest GSV per unit area using different proportions of field sample data, ALS cover data, and wall-to-wall Landsat data. Model-assisted estimators were used with NFI sample data in an Italian study area using 10 RS predictors, specifically the seven Landsat 7 ETM+ bands and three fine-resolution metrics based on ALS-derived canopy height. We found that relative to the standard simple expansion estimator, the model-assisted estimators produced relative efficiency of 1.16 when using only Landsat data and relative efficiencies as great as 1.33 when using increasing levels of ALS coverage.

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.

In the context of the potential future use of unmanned ground vehicles for forest inventories, we present the first experiences with SPOT, a legged robot equipped with a LiDAR instrument and several cameras that have been used with a teleoperation approach for single-tree detection and measurements. This first test was carried out using the default LiDAR system (the so-called enhanced autonomy payload-EAP, installed on the board of SPOT to guide autonomous movements) to understand advantages and limitations of this platform to support forest inventory activities. The test was carried out in the Vallombrosa forest (Italy) by assessing different data acquisition methods. The first results showed that EAP LiDAR generated noisy point clouds where only large trees (DBH ≥ 20 cm) could be identified. The results showed that the accuracy in tree identification and DBH measurements were strongly influenced by the path used for data acquisition, with average errors in tree positioning no less than 1.9 m. Despite this, the best methods allowed the correct identification of 97% of large trees.

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.

Several political initiatives aim to achieve net-zero emissions by the middle of the twenty-first century. In this context, forests are crucial as a carbon sink to store unavoidable emissions. Assessing the carbon sequestration potential of forest ecosystems is pivotal to the availability of accurate forest variable estimates for supporting international reporting and appropriate forest management strategies. Spatially explicit estimates are even more important for Mediterranean countries such as Italy, where the capacity of forests to act as sinks is decreasing due to climate change. This study aimed to develop a spatial approach to obtain high-resolution maps of Italian forest above-ground biomass (ITA-BIO) and current annual volume increment (ITA-CAI), based on remotely sensed and meteorological data. The ITA-BIO estimates were compared with those obtained with two available biomass maps developed in the framework of two international projects (i.e., the Joint Research Center and the European Space Agency biomass maps, namely, JRC-BIO and ESA-BIO). The estimates from ITA-BIO, JRC-BIO, ESA-BIO, and ITA-CAI were compared with the 2nd Italian NFI (INFC) official estimates at regional level (NUT2). The estimates from ITA-BIO are in good agreement with the INFC estimates (R2 = 0.95, mean difference = 3.8 t ha−1), while for JRC-BIO and ESA-BIO, the estimates show R2 of 0.90 and 0.70, respectively, and mean differences of 13.5 and of 21.8 t ha−1 with respect to the INFC estimates. ITA-CAI estimates are also in good agreement with the INFC estimates (R2 = 0.93), even if they tend to be slightly biased. The produced maps are hosted on a web-based forest resources management Decision Support System developed under the project AGRIDIGIT (ForestView) and represent a key element in supporting the new Green Deal in Italy, the European Forest Strategy 2030 and the Italian Forest Strategy.

Large-scale forest monitoring benefits greatly from change detection analysis based on remote sensing data because it enables characterizing forest dynamics of disturbance and recovery by detecting both gradual and abrupt changes on Earth’s surface. In this study, two of the main disturbances occurring in Mediterranean forests, harvesting operations and forest fires, were analyzed through the analysis of Landsat Times Series images in a case study in Central Italy (Tuscany region). Disturbances were characterized based on their distinct temporal behaviors before and after the event: a period of 20 years (1999–2018) was used to extract and analyze at pixel level spectral trajectories for each disturbance and produce descriptive temporal trends of the phenomena. Recovery metrics were used to characterize both short- (5 years) and long-term aspects of recovery for harvested and burned areas. Spectral, recovery, and trend analysis metrics were then used with the Random Forest classifier to differentiate between the two disturbance classes and to investigate their potential as predictors. Among spectral bands, the Landsat SWIR 1 band proved the best to detect areas interested by harvesting, while forest fires were better detected by the SWIR 2 band; among spectral indices, the NBR scored as the best for both classes. On average, harvested areas recovered faster in both short- and long-term aspects and showed less variability in the magnitude of the disturbance event and recovery rate over time. This tendency is confirmed by the results of the classifier, which obtained an overall accuracy of 98.6%, and identified the mean of the post-disturbance values of the trend as the best predictor to differentiate between disturbances.

Forests play a prominent role in the battle against climate change, as they absorb a relevant part of human carbon emissions. However, precisely because of climate change, forest disturbances are expected to increase and alter forests’ capacity to absorb carbon. In this context, forest monitoring using all available sources of information is crucial. We combined optical (Landsat) and photonic (GEDI) data to monitor four decades (1985–2019) of disturbances in Italian forests (11 Mha). Landsat data were confirmed as a relevant source of information for forest disturbance mapping, as forest harvestings in Tuscany were predicted with omission errors estimated between 29% (in 2012) and 65% (in 2001). GEDI was assessed using Airborne Laser Scanning (ALS) data available for about 6 Mha of Italian forests. A good correlation (r2 = 0.75) between Above Ground Biomass Density GEDI estimates (AGBD) and canopy height ALS estimates was reported. GEDI data provided complementary information to Landsat. The Landsat mission is capable of mapping disturbances, but not retrieving the three-dimensional structure of forests, while our results indicate that GEDI is capable of capturing forest biomass changes due to disturbances. GEDI acquires useful information not only for biomass trend quantification in disturbance regimes but also for forest disturbance discrimination and characterization, which is crucial to further understanding the effect of climate change on forest ecosystems.

Mapping forest disturbances is an essential component of forest monitoring systems both to support local decisions and for international reporting. Between the 28 and 29 October 2018, the VAIA storm hit the Northeast regions of Italy with wind gusts exceeding 200 km h−1. The forests in these regions have been seriously damaged. Over 490 Municipalities in six administrative Regions in Northern Italy registered forest damages caused by VAIA, that destroyed or intensely damaged forest stands spread over an area of 67,000 km2. The present work tested the use of two continuous change detection algorithms, i.e., the Bayesian estimator of abrupt change, seasonal change, and trend (BEAST) and the continuous change detection and classification (CCDC) to map and estimate forest windstorm damage area using a normalized burned ration (NBR) time series calculated on three years Sentinel-2 (S2) images collection (i.e., January 2017–October 2019). We analyzed the accuracy of the maps and the damaged forest area using a probability-based stratified estimation within 12 months after the storm with an independent validation dataset. The results showed that close to the storm (i.e., 1 to 6 months November 2018–March 2019) it is not possible to obtain accurate results independently of the algorithm used, while accurate results were observed between 7 and 12 months from the storm (i.e., May 2019–October 2019) in terms of Standard Error (SE), percentage SE (SE%), overall accuracy (OA), producer accuracy (PA), user accuracy (UA), and gmean for both BEAST and CCDC (SE < 3725.3 ha, SE% < 9.69, OA > 89.7, PA and UA > 0.87, gmean > 0.83).

Information about forest cover and its characteristics are essential in national and international forest inventories, monitoring programs, and reporting activities. Two of the most common forest variables needed to support sustainable forest management practices are forest cover area and growing stock volume (GSV m3 ha−1). Nowadays, national forest inventories (NFI) are complemented by wall-to-wall maps of forest variables which rely on models and auxiliary data. The spatially explicit prediction of GSV is useful for small-scale estimation by aggregating individual pixel predictions in a model-assisted framework. Spatial knowledge of the area of forest land is an essential prerequisite. This information is contained in a forest mask (FM). The number of FMs is increasing exponentially thanks to the wide availability of free auxiliary data, creating doubts about which is best-suited for specific purposes such as forest area and GSV estimation. We compared five FMs available for the entire area of Italy to examine their effects on the estimation of GSV and to clarify which product is best-suited for this purpose. The FMs considered were a mosaic of local forest maps produced by the Italian regional forest authorities; the FM produced from the Copernicus Land Monitoring System; the JAXA global FM; the hybrid global FM produced by Schepaschencko et al., and the FM estimated from the Corine Land Cover 2006. We used the five FMs to mask out non-forest pixels from a national wall-to-wall GSV map constructed using inventory and remotely sensed data. The accuracies of the FMs were first evaluated against an independent dataset of 1,202,818 NFI plots using four accuracy metrics. For each of the five masked GSV maps, the pixel-level predictions for the masked GSV map were used to calculate national and regional-level model-assisted estimates. The masked GSV maps were compared with respect to the coefficient of correlation (ρ) between the estimates of GSV they produced (both in terms of mean and total of GSV predictions within the national and regional boundaries) and the official NFI estimates. At the national and regional levels, the model-assisted GSV estimates based on the GSV map masked by the FM constructed as a mosaic of local forest maps were closest to the official NFI estimates with ρ = 0.986 and ρ = 0.972, for total and mean GSV, respectively. We found a negative correlation between the accuracies of the FMs and the differences between the model-assisted GSV estimates and the NFI estimate, demonstrating that the choice of the FM plays an important role in GSV estimation when using the model-assisted estimator.

Different forest types based on different tree species composition may have similar spectral signatures if observed with traditional multispectral satellite sensors. Hyperspectral imagery, with a more continuous representation of their spectral behavior may instead be used for their classification. The new hyperspectral Precursore IperSpettrale della Missione Applicativa (PRISMA) sensor, developed by the Italian Space Agency, is able to capture images in a continuum of 240 spectral bands ranging between 400 and 2500 nm, with a spectral resolution smaller than 12 nm. The new sensor can be employed for a large number of remote sensing applications, including forest types discrimination. In this study, we compared the capabilities of the new PRISMA sensor against the well-known Sentinel-2 Multi-Spectral Instrument (MSI) in recognition of different forest types through a pairwise separability analysis carried out in two study areas in Italy, using two different nomenclature systems and four separability metrics. The PRISMA hyperspectral sensor, compared to Sentinel-2 MSI, allowed for a better discrimination in all forest types, increasing the performance when the complexity of the nomenclature system also increased. PRISMA achieved an average improvement of 40% for the discrimination between two forest categories (coniferous vs. broadleaves) and of 102% in the discrimination between five forest types based on main tree species groups.