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.
•The study aims at constructing wall-to-wall estimates of forest growing stock (GSV).•We combine NFI plot data, remotely sensed and auxiliary variables.•We applied the methodology in Mediterranean ...Forest.•We create a wall-to-wall GSV forest map in a large test area.•The GSV map was used to produce model-based estimates of GSV at small scale.
Spatial predictions of forest variables are required for supporting modern national and sub-national forest planning strategies, especially in the framework of a climate change scenario. Nowadays methods for constructing wall-to-wall maps and calculating small-area estimates of forest parameters are becoming essential components of most advanced National Forest Inventory (NFI) programs. Such methods are based on the assumption of a relationship between the forest variables and predictor variables that are available for the entire forest area. Many commonly used predictors are based on data obtained from active or passive remote sensing technologies. Italy has almost 40% of its land area covered by forests. Because of the great diversity of Italian forests with respect to composition, structure and management and underlying climatic, morphological and soil conditions, a relevant question is whether methods successfully used in less complex temperate and boreal forests may be applied successfully at country level in Italy.
For a study area of more than 48,657 km2 in central Italy of which 43% is covered by forest, the study presents the results of a test regarding wall-to-wall, spatially explicit estimation of forest growing stock volume (GSV) based on field measurement of 1350 plots during the last Italian NFI. For the same area, we used potential predictor variables that are available across the whole of Italy: cloud-free mosaics of multispectral optical satellite imagery (Landsat 5 TM), microwave sensor data (JAXA PALSAR), a canopy height model (CHM) from satellite LiDAR, and auxiliary variables from climate, temperature and precipitation maps, soil maps, and a digital terrain model.
Two non-parametric (random forests and k-NN) and two parametric (multiple linear regression and geographically weighted regression) prediction methods were tested to produce wall-to-wall map of growing stock volume at 23-m resolution. Pixel level predictions were used to produce small-area, province-level model-assisted estimates. The performances of all the methods were compared in terms of percent root mean-square error using a leave-one-out procedure and an independent dataset was used for validation. Results were comparable to those available for other ecological regions using similar predictors, but random forests produced the most accurate results with a pixel level R2 = 0.69 and RMSE% = 37.2% against the independent validation dataset. Model-assisted estimates were more precise than the original design-based estimates provided by the NFI.
Forest canopy gaps are important to ecosystem dynamics. Depending on tree species, small canopy openings may be associated with intra-crown porosity and with space among crowns. Yet, literature on ...the relationships between very fine-scaled patterns of canopy openings and biodiversity features is limited. This research explores the possibility of: (1) mapping forest canopy gaps from a very high spatial resolution orthomosaic (10 cm), processed from a versatile unmanned aerial vehicle (UAV) imaging platform, and (2) deriving patch metrics that can be tested as covariates of variables of interest for forest biodiversity monitoring. The orthomosaic was imaged from a test area of 240 ha of temperate deciduous forest types in Central Italy, containing 50 forest inventory plots each of 529 m2 in size. Correlation and linear regression techniques were used to explore relationships between patch metrics and understory (density, development, and species diversity) or forest habitat biodiversity variables (density of micro-habitat bearing trees, vertical species profile, and tree species diversity). The results revealed that small openings in the canopy cover (75% smaller than 7 m2) can be faithfully extracted from UAV red, green, and blue bands (RGB) imagery, using the red band and contrast split segmentation. The strongest correlations were observed in the mixed forests (beech and turkey oak) followed by intermediate correlations in turkey oak forests, followed by the weakest correlations in beech forests. Moderate to strong linear relationships were found between gap metrics and understory variables in mixed forest types, with adjusted R2 from linear regression ranging from 0.52 to 0.87. Equally strong correlations in the same forest types were observed for forest habitat biodiversity variables (with adjusted R2 ranging from 0.52 to 0.79), with highest values found for density of trees with microhabitats and vertical species profile. In conclusion, this research highlights that UAV remote sensing can potentially provide covariate surfaces of variables of interest for forest biodiversity monitoring, conventionally collected in forest inventory plots. By integrating the two sources of data, these variables can be mapped over small forest areas with satisfactory levels of accuracy, at a much higher spatial resolution than would be possible by field-based forest inventory solely.
The development of laser scanning technologies has gradually modified methods for forest mensuration and inventory. The main objective of this study is to assess the potential of integrating ALS and ...TLS data in a complex mixed Mediterranean forest for assessing a set of five single-tree attributes: tree position (TP), stem diameter at breast height (DBH), tree height (TH), crown base height (CBH) and crown projection area radii (CPAR). Four different point clouds were used: from ZEB1, a hand-held mobile laser scanner (HMLS), and from FARO® FOCUS 3D, a static terrestrial laser scanner (TLS), both alone or in combination with ALS. The precision of single-tree predictions, in terms of bias and root mean square error, was evaluated against data recorded manually in the field with traditional instruments. We found that: (i) TLS and HMLS have excellent comparable performances for the estimation of TP, DBH and CPAR; (ii) TH was correctly assessed by TLS, while the accuracy by HMLS was lower; (iii) CBH was the most difficult attribute to be reliably assessed and (iv) the integration with ALS increased the performance of the assessment of TH and CPAR with both HMLS and TLS.
A forest inventory is often carried out using airborne laser data combined with ground measured reference data. Traditionally, the ground reference data have been collected manually with a caliper ...combined with land surveying equipment. During recent years, studies have shown that the caliper can be replaced by equipment and methods that capture the ground reference data more efficiently. In this study, we compare three different ground based laser measurement methods: terrestrial laser scanner, handheld laser scanner and a backpack laser scanner. All methods are compared with traditional measurements. The study area is located in southeastern Norway and divided into seven different locations with different terrain morphological characteristics and tree density. The main tree species are boreal, dominated by Norway spruce and Scots pine. To compare the different methods, we analyze the estimated tree stem diameter, tree position and data capture efficiency. The backpack laser scanning method captures the data in one operation. For this method, the estimated diameter at breast height has the smallest mean differences of 0.1 cm, the smallest root mean square error of 2.2 cm and the highest number of detected trees with 87.5%, compared to the handheld laser scanner method and the terrestrial laser scanning method. We conclude that the backpack laser scanner method has the most efficient data capture and can detect the largest number of trees.
Nowadays, forest inventories are frequently carried out using a combination of field measurements and remote sensing data, often acquired with light detection and ranging (LiDAR) sensors. Several ...studies have investigated how three-dimensional laser scanning point clouds from different platforms can be used to acquire information traditionally collected with forest instruments, such as hypsometers and callipers to detect single-tree attributes like tree height and diameter at the breast height. The present study has tested the performances of the ZEB1 instrument, a type of hand-held mobile laser scanner, for single-tree attributes estimation in pure Castanea sativa Mill. stands cultivated for fruit production in Central Italy. In particular, the influence of walking scan path density on single-tree attributes estimation (number of trees, tree position, diameter at breast height, tree height, and crown base height) was investigated to test the efficiency of field measures. The point clouds were acquired by walking along straight lines drawn with different spacing: 10 and 15 m apart. A single-tree scan approach, which included walking with the instrument around each tree, was used as reference data. In order to evaluate the efficiency of the survey, the influence of the walking scan path was discussed in relation to the accuracy of single-tree attributes estimation, as well as the time and cost needed for data acquisition, pre-processing, and analysis. Our results show that the 10 m scan path provided the best results, with an omission error of 6%; the assessment of single-tree attributes was successful, with values of the coefficient of determination and the relative root mean square error similar to other studies. The 10 m scan path has also proved to decrease the costs by about €14 for data pre-processing, and a saving of time for data acquisition and data analysis of about 37 min compared to the reference data.
A Landsat time series has been recognized as a viable source of information for monitoring and assessing forest disturbances and for continuous reporting on forest dynamics. This study focused on ...developing automated procedures for detecting disturbances in Mediterranean coppice forests which are characterized by rapid regrowth after a cut. Specifically, new methods specific to Mediterranean coppice forests are needed for mapping clearcut disturbances over time and for estimating related indicators in the context of Sustainable Forest Management and Biodiversity International monitoring frameworks. The aim of this work was to develop a new change detection algorithm for mapping clearcut disturbances in Mediterranean coppice forests with Landsat time series (LTS) using a short time window. Accuracy for the new algorithm, characterized as the Two Thresholds Method (TTM), was evaluated using an independent clearcut reference dataset over a temporal period of the 13 years between 2001 and 2013. TTM was also evaluated against two benchmark approaches: (i) LandTrendr, and (ii) the forest loss category of the Global Forest Change Map. Overall Accuracy for LandTrendr and TTM were greater than 0.94. Meanwhile, smaller accuracies were always obtained for the GFC. In particular, Producer’s Accuracy ranged between 0.45 and 0.84 for TTM and between 0.49 and 0.83 for LT, while for the GFC, PA ranged between 0 and 0.38. User’s Accuracy ranged between 0.86 and 0.96 for TTM and between 0.73 and 0.91 for LT, while for the GFC UA ranged between 0.19 and 1.00. Moreover, to illustrate the utility of TTM for mapping clearcut disturbances in Mediterranean coppice forests, we applied TTM to a Landsat scene that covered almost the entirety of the Tuscany region in Italy.
To combat global deforestation, monitoring forest disturbances at sub-annual scales is a key challenge. For this purpose, the new Planetscope nano-satellite constellation is a game changer, with ...a revisit time of 1 day and a pixel size of 3-m.
We present a near-real time forest disturbance alert system based on PlanetScope imagery: the Thresholding Rewards and Penances algorithm (TRP). It produces a new forest change map as soon as a new PlanetScope image is acquired.
To calibrate and validate TRP, a reference set was constructed as a complete census of five randomly selected study areas in Tuscany, Italy. We processed 572 PlanetScope images acquired between 1 May 2018 and 5 July 2019.
TRP was used to construct forest change maps during the study period for which the final user's accuracy was 86% and the final producer's accuracy was 92%. In addition, we estimated the forest change area using an unbiased stratified estimator that can be used with a small sample of reference data. The 95% confidence interval for the sample-based estimate of 56.89 ha included the census-based area estimate of 56.19 ha.
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.
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).