A completely data-driven, design-based sampling strategy is proposed for mapping a forest attribute within the spatial units tessellating a survey region. Based on sample data, a model is selected, and model parameters are estimated using least-squares criteria for predicting the attribute of interest within units as a linear function of a set of auxiliary variables. The spatial interpolation of residuals arising from model predictions is performed by inverse distance weighting. The leave-one-out cross validation procedure is adopted for selecting the smoothing parameter used for interpolation. The densities of the attributes of interest within units are estimated by summing predictions and interpolated residuals. Finally, density estimates are rescaled to match the total estimate over the survey region obtained by the traditional regression estimator with the total estimate obtained from the map as the sum of the density estimates within units. A bootstrap procedure accounts for the uncertainty. The consistency of the strategy is proven by incorporating previous results. A simulation study is performed and an application for mapping wood volume densities in the forest estate of Rincine (Central Italy) is described. © 2022 John Wiley & Sons Ltd.

This paper presents a new co-registration procedure of complementary point clouds captured by both Terrestrial (TLS) and Airborne Laser Scanning (ALS) technologies. Starting from the geographic position of the TLS point cloud, a geometric features recognition algorithm, which evaluates digital terrain models obtained from both ALS and TLS, was developed and implemented in a new GIS software (ForeSight®). As a case study, we tested this new approach using point clouds acquired from both hand-held mobile TLS and ALS sensors over 24 test sites located in a protected area in southern Italy, with the ultimate goal of characterizing the different forest stand structures. From each aligned point cloud, a plot-level spatially explicit index (Enhanced Structural Spatial Index, ESCI) was derived to assess the three-dimensional structure of the considered forest stands. Then, we compared structural features derived from the ESCI index with different computed ALS metrics. Finally, the most correlated ALS metrics were used as predictors to produce an ESCI-map of the entire region of interest. © 2021 Elsevier B.V.

Vegetation structure is one of the key factors in forest ecosystems. Especially understory structure has major implications for wildlife habitat selection, reproduction, and survival. Structural indices traditionally used to characterize understory vegetation are based on field vegetation surveys describing general features such as leaf area index (LAI), canopy cover or vegetation height, hiding much of the three-dimensional vegetation structure complexity. The application of terrestrial laser scanning (TLS) in forest ecological and management applications is becoming more effective. In this study, we use TLS data to quantify spatial attributes of forest subcanopy in four different forest strata ranging from 0.5 m to 10 m from the ground. We collected data in 12 plots of mature European beech (Fagus sylvatica L.) forests and 12 plots of mature black pine (Pinus nigra subsp. laricio Maire) forests, located in the Sila National Park, Italy. We propose a TLS-based approach to estimate a fine-scale vegetation density using the Plant Density Index (PDI) and to test the PDI at different height classes. We found a significant relationship between the PDI and the number of trees belonging to the dominant layer, using the Spearman correlation coefficient (r = 0.83, p<inf>val</inf> = 0.001). Basing on PDI values, a cluster analysis of the four subcanopy strata was carried out for deriving clusters of structurally homogeneous forest plots. Results identified three clusters in terms of the vegetation features in the horizontal height classes: the first cluster primarily includes Beech forests characterized by plots with the highest tree densities; the second one includes both Beech and Pine forests characterized by dense ground vegetation and shrubs and an intermediate tree density; the third group is represented by Pine forests with massive presence of vegetation lower strata and moderate tree density. Then, PCA allowed identifying the relationship between the considered subcanopy layers and forest plots. © 2021 Elsevier B.V.

Vegetation structure is a key determinant of species distribution and diversity. Compared to traditional methods, the use of Terrestrial Laser Scanning (TLS) has allowed massive amounts of point cloud data collected for quantifying three-dimensional habitat properties at increasing spatial and temporal scales. We used TLS to characterize the forest plots across a broad range of forest structural diversity, located in the Sila National Park, South Italy. The dataset reports data collected in 24 15-m-radius circular plots, 12 of which were dominated by beech (Fagus sylvatica L.) and 12, by black pine (Pinus nigra subsp. laricio). In detail, this work provides dataset of i) plot-level attributes calculated from raw data, such as the number of trees, ii) tree-level data, comprising a total of 1709 trees, with information related to field-based forest inventory such as the diameter at breast height (DBH), and iii) plot-level information related to the time for conducting both traditional field- and TLS-based forest inventories. Compared to traditional methods, the use of TLS allows a very high-resolution quantification of the 3D forest structural properties, also reducing the time for conducting forest inventories. © 2020

Accurate estimates of crown cover (CC) are central for a wide range of forestry studies. As direct measurements do not exist to retrieve this variable in the field, CC is conventionally determined from optical measurements as the complement of gap fraction close to the zenith. As an alternative to passive optical measurements, active sensors like terrestrial Light Detection And Ranging (LiDAR) allows for characterizing the 3D canopy structure with unprecedented detail. We evaluated the reliability of terrestrial LiDAR (TLS) to estimate CC using a voxel-based approach. Specifically, we tested how different voxel sizes (ranging from 5-20 cm) and voxel densities (1-9 points/dm³) influenced the retrieval of CC. Results were compared against benchmark values obtained from digital cover photography (DCP). The trial was performed in hybrid poplar plantations in Northern Italy. Results indicate that TLS can be used for obtaining accurate estimates of CC, but the choice of voxel size and point density is critical for achieving such accuracy. In hybrid poplars, the best performance was obtained using voxel size of 10 cm and point density of 8 points/dm³. The combined ability of measuring and mapping CC also holds great potential to use TLS for calibrating and upscaling results using coarser-scale remotely sensed products. © 2021 Centro di Ricerca per la Selvicoltura, Consiglio per la Ricerca in Agricoltura e l'Analisi dell'Economia Agraria. All rights reserved.

Poplars are one of the most widespread fast-growing tree species used for forest plantations. Owing to their distinct features (fast growth and short rotation) and the dependency on the timber price market, poplar plantations are characterized by large inter-annual fluctuations in their extent and distribution. Therefore, monitoring poplar plantations requires a frequent update of information–not feasible by National Forest Inventories due to their periodicity–achievable by remote sensing systems applications. In particular, the new Sentinel-2 mission, with a revisiting period of 5 days, represents a potentially efficient tool for meeting this need. In this paper, we present a deep learning approach for mapping poplar plantations using Sentinel-2 time series. A reference dataset of poplar plantations was available for a large study area of more than 46,000 km² in Northern Italy and served as training and testing data. Two classification methods were compared: (1) a fully connected neural network (also called multilayer perceptron), and (2) a traditional logistic regression. The performance of the two approaches was estimated through bootstrapping procedure with a confidence interval of 99%. Results indicated for deep learning an omission error rate of 2.77%±2.76%, showing improvements compared to logistic regression, omission error rate = 8.91%±4.79%. © 2021 Informa UK Limited, trading as Taylor & Francis Group.

Accurate estimates of canopy cover (CC) are central for a wide range of forestry studies. As direct measurements are impractical, indirect optical methods have often been used to estimate CC from the complement of gap fraction measurements obtained with restricted-view sensors. In this short note we evaluated the influence of the image pixel resolution (ground sampling distance; GSD) on CC estimation in poplar plantations obtained from field (cover photography; GSD < 1 cm), unmanned aerial (UAV; GSD <10 cm) and satellite (Sentinel-2; GSD = 10 m) imagery. The trial was conducted in poplar tree plantations in Northern Italy, with varying age and canopy cover. Results indicated that the coarser resolution available from satellite data is suitable to obtain estimates of canopy cover, as compared with field measurements obtained from cover photography; therefore, S2 is recommended for larger scale monitoring and routine assessment of canopy cover in poplar plantations. The higher resolution of UAV compared with Sentinel-2 allows finer assessment of canopy structure, which could also be used for calibrating metrics obtained from coarser-scale remote sensing products, avoiding the need of ground measurements. © 2021 Centro di Ricerca per la Selvicoltura, Consiglio per la Ricerca in Agricoltura e l'Analisi dell'Economia Agraria. All rights reserved.

Accurate and frequently updated tree volume estimates are required for poplar plantations, which are characterized by fast growth rate and short rotation. In this study, we tested the potential of terrestrial laser scanning (TLS) as a reliable method for developing nondestructive tree volume allometries in poplar plantations. The trial was conducted in Italy, where 4- to 10-year-old hybrid plantations were sampled to develop tree crown volume allometry in leaf-on conditions, tree stem volume, and height-diameter allometries in leaf-off conditions. We tested one-entry models based on diameter and two-entry models based on both diameter and height. Model performance was assessed by residual analysis. Results indicate that TLS can provide accurate models of tree stem and crown volume, with percentage of root-mean-square error of about 20 percent and 15 percent, respectively. The inclusion of height does not bring relevant improvement in the models, so that only diameter can be used to predict tree stem and crown volume. The TLS-measured stem volume estimates agreed with an available formula derived from harvesting. We concluded that TLS is a reliable method for developing nondestructive volume allometries in poplar plantations and holds great potential to enhance conventional tree inventory and monitoring. © The Author(s) 2020. Published by Oxford University Press on behalf of the Society of American Foresters. All rights reserved.

In the recent decades, growing demand for wood products, combined with efforts to conserve natural forests, has supported a steady increase in the global extent of planted forests. In this paper, a two-phase sampling strategy for large-scale assessment of hybrid poplar plantations in Northern Italy was implemented. The first phase was performed by means of tessellation stratified sampling on high-resolution remotely sensed imagery, covering the survey area by a grid of regular polygons of equal size and randomly and independently selecting one point per quadrat. All the plantations spotted by at least one sample point were selected. In the second phase, we randomly chosen a subset of plantations by stratified sampling that were visited on the ground to collect qualitative and quantitative attributes. The resulting estimates were reliable, and the survey demonstrated relatively easy to be implemented and replicated. These considerations support the use of the proposed sampling strategy to frequently update information on fast-growing forest plantations within agricultural farms, like hybrid poplar crops. Moreover, the results of the case study here presented highlight the relevance of hybrid poplar plantations in Italy, in the context of sustainable development strategies under a green economy perspective. © 2020, Springer-Verlag GmbH Germany, part of Springer Nature.

The objective of this research is to test and evaluate hyperspectral and lidar data to derive information on tree species dominance and above ground biomass in the Białowieża Forest in Poland. This forest is threatened by climate change, fire, bark beetles attacks, and logging, with changes in species composition and dominance. In this conservation valuable area, the monitoring of forest resources is thus critical. Results indicate that vegetation indices from hyperspectral data can support species dominance detection: using a Classification and Regression Trees algorithm the three main plot types (dominated by Deciduous, Spruce, and Pines species) were classified with an Overall Accuracy > 0.9. The accuracy decreased when a ‘Mixed’ group was added to account for very heterogeneous plots, and plots dominated by Spruce were not correctly detected. Hyperspectral vegetation indices were also used to estimate the level of species dominance in the forest plots, using a Multivariate Multiple Linear Regression model; the obtained accuracy varied according to groups, being higher for Deciduous (R² = 0.87), compared to Pines (R² = 0.61), and to Spruce-dominated plots (R² = 0.37). Lidar data were employed to estimate above ground biomass, using an exponential regression model; overall the R² resulted equal to 0.66 but ranged from 0.57 to 0.78 when considering subgroups according to species dominance; the addition of hyperspectral vegetation indices improved the result only for Pines. The illustrated methods provide a reliable description of important forest characteristics and simplify resource monitoring, supporting local authorities to address the challenges imposed by climate change and other forest threats. © 2020 The Authors

Light detection and ranging (Lidar) is considered the most advanced technology to assess forest aboveground biomass (AGB). Currently, this technology is shared by different sensors ranging from ground [terrestrial laser scanning (TLS)], airborne [aerial laser scanning (ALS)] up to spaceborne ones, which entail different spatial scales. However, few studies tested the simultaneous and combined use of Lidar to estimate AGB, linking ground measurements up to satellite observations. To fill this gap, we performed a study in two Mediterranean forest types [i.e., mountainous beech (Fagus sylvatica) and black pine (Pinus nigra subsp. laricio)] with contrasting structures (i.e., broadleaf versus needleleaf forests), where field inventory, TLS, ALS, and the recent spaceborne Global Ecosystem Dynamics Investigation (GEDI) data were simultaneously acquired. A three-step procedure was followed, which involved (i) the validation of AGB estimates obtained from TLS against reference values obtained from conventional field inventory; (ii) the calibration and validation of AGB estimates derived from ALS against TLS measurements, and (iii) the calibration and validation of AGB estimates derived from GEDI against mapped AGB values obtained from ALS. Our main results indicated that TLS provides consistent measurements of AGB as compared with field measurements (R2 ranged between 0.6 and 0.9 and root-mean-square error ranged between 29% and 49%), indicating its potential as ground reference for airborne Lidar observations. The combined availability of ground, airborne, and spaceborne observations is suitable to link ground measurements up to satellite observations. Differences in Lidar performance between needleleaf and broadleaf forests are also considered and discussed. © 2020 Society of Photo-Optical Instrumentation Engineers (SPIE).

Accurate in situ estimates of leaf area index (LAI) are essential for a wide range of ecological studies and applications. Due to the destructiveness and impracticality of direct measurements, indirect optical methods have mostly been used in the field to derive estimates of LAI from gap fraction measurements. Terrestrial laser scanning (TLS) is strongly supporting use of this active technology, which possesses several advantages compared to passive sensors. However, edge effects and partial beam interceptions are significantly challenges for the accurate retrieval of gap fraction from 3D point cloud data available from TLS, particularly in phase-shift instruments, which in turns require point cloud filtering to correct erroneous point measurements. As the limitations above influences the point cloud, we proposed a new method which is based only on the laser return intensity (LRI) information derived from raw TLS data, which are used to generate 2D intensity images. The intensity image contains all the unfiltered LRI information captured by TLS, which is used to separate gap from non-gap pixels, using a procedure comparable to the standard image analysis processing of digital hemispherical images. This allows a theoretically consistent comparison between active and passive optical measurements of gap fraction across all the zenith angle range. The method was tested in real and simulated forests. Gap fraction, canopy openness and effective leaf area index derived from real and simulated intensity TLS images were compared with those obtained using digital hemispherical photography (DHP). Results indicated that the intensity, image-based method outperformed DHP, as the higher pixel resolution of the intensity images and the larger distance covered by TLS allowed detection of many small canopy elements, particularly at higher zenith angles (longer optical distance), which are not detected in DHP. The main findings support the reliability of the intensity, image-based method to standardize protocols for TLS phase-shift scan data processing and use of the produced canopy estimates as a benchmark for passive optical measurements. © 2019 Elsevier B.V.

Growing demand for wood products, combined with efforts to conserve natural forests, have supported a steady increase in the global extent of planted forests. Here, a two-phase sampling strategy for large-scale assessment of the total area and the total wood volume of fast-growing forest tree crops within agricultural land is presented. The first phase is performed using tessellation stratified sampling on high-resolution remotely sensed imagery and is sufficient for estimating the total area of plantations by means of a Monte Carlo integration estimator. The second phase is performed using stratified sampling of the plantations selected in the first phase and is aimed at estimating total wood volume by means of an approximation of the first-phase Horvitz-Thompson estimator. Vegetation indices from Sentinel-2 are exploited as freely available auxiliary information in a linear regression estimator to improve the design-based precision of the estimator based on the sole sample data. Estimators of the totals and of the design-based variances of total estimators are presented. A simulation study is developed in order to check the design-based performance of the two alternative estimators under several artificial distributions supposed for poplar plantations (random, clustered, spatially trended). An application in Northern Italy is also reported. The regression estimator turns out to be invariably better than that based on the sole sample information. Possible integrations of the proposed sampling scheme with conventional national forest inventories adopting tessellation stratified sampling in the first phase are discussed. © 2020, Finnish Society of Forest Science. All rights reserved.

Mixed forests are particularly interesting for forest structure and diversity analyses, as higher complexity and diversity can be expected in these forests compared to pure ones. Integrating different approaches in the analyses of structure and diversity in these forests can provide complementary information on non-spatial, spatial and functional diversity patterns. The study aimed at evaluating the spatio-temporal dynamics in forest structure and diversity in a semi-natural mixed oak-hornbeam floodplain forest. All standing trees were mapped and inventoried in 1995, 2005 and 2016 in three 1-ha mixed forest stands, with different soil moisture regime (xeric, mesic, moist conditions). Traditional, non-spatial structure and diversity measures were coupled with spatially-explicit and functional diversity measures. Results indicated that the three stands showed limited variation in stand structure and similar non-spatial diversity attributes, despite the different species composition. Only the extension to spatial and functional analyses was able to reveal more pronounced differences of diversity patterns, as higher complexity, species mingling, and functional tree complementarity was observed in the moister stand. These findings support use of spatially-explicit measurements in traditional inventory measurement protocols to allow more refined analysis of diversity patterns. On the other hand, functional diversity can be easily implemented in diversity analyses, as it requires species abundance information (which is traditionally collected in forest inventory) and species-specific tree traits which can be inferred from literature. © 2019 Elsevier Ltd

Accurate inventory and mapping of olive (Olea europaea L.) tree attributes represents a central issue to support the olive production system. With reference to the cultivation, there is a high heterogeneity and complexity in the cultivation of olive trees, which is reflected in the large variability in olive grove surfaces. This poses some challenge in accurately estimating olive tree attributes via traditional inventory approaches, as commonly adopted in national forest inventory. From a methodological point of view, the complexity and heterogeneity of olive tree groves can be comparable to the problem of accurately estimating tree outside forests (TOF) attributes. In this study, we tested whether a plot sampling approach formerly developed for TOF is suitable for estimating olive tree attributes at large scale. We tested this approach in a case study where the census of the olive crop area and the number of olive groves was conducted from photo-interpretation of high resolution aerial orthoimagery, used as benchmark to test the effectiveness of the plot sampling approach. The main result of this study is that the plot sampling method can be applied for estimating olive tree attributes. Our obtained RSEs were below 20%, with a limited sampling effort of about 6% of the studied population; the obtained RSEs were below 6% when increasing sampling up to about 21% the studied population. Using robust statistical procedures among countries, should allow obtaining harmonized and comparable information, which can increase the knowledge of olive geographical distribution and structure at its relevant Mediterranean scale. © 2019, Italian Society of Remote Sensing. All rights reserved.

The value of wood for different timber assortments can vary by a factor of ten, optimization of stems assortment is hence a key element in the wood products supply chain, particularly for plantations. ‘Taper functions’ are commonly used in other countries to tackle this issue. In Italy, this approach has not yet entered operational use. These functions are developed based on measures of stem diameters taken at different distances from the base. Such measurements are commonly taken felling the tree and using a tape meter and the tree calliper, clearly assuming some approximations. This research assesses the advantages, in terms of assortments evaluation, that can be obtained if the diameters at different heights are extracted adequately processing Terrestrial Laser Scanning (TLS) output. TLS data have been collected, in a poplar plantation, on 36 trees distributed on three stands with different plantation densities in Padana Plane, Italy. The estimated profiles display a very high variability with an average of 1.8 cm of lateral compression. The results from this study demonstrate the potential and feasibility of estimating bole eccentricity by TLS, providing preliminary tools that will hopefully favour the diffusion of taper functions in operational environments. © 2019, Italian Society of Remote Sensing. All rights reserved.

The value of wood for different timber assortments can vary by a factor of ten. Optimization of stem assortments is, hence, a key element in the wood products supply chain, particularly for plantations. 'Taper functions' are commonly used in other countries to tackle this issue. In Italy, this approach has not yet entered operational use. These functions are developed based on measures of stem diameters taken at different distances from the base. Such measurements are commonly taken felling the tree and using a tape meter and tree caliper, clearly assuming some approximations. This research assesses the advantages, in terms of assortments evaluation, that can be obtained if the diameters at different heights are extracted adequately to process terrestrial laser scanning (TLS) output. TLS data have been collected, in a poplar plantation, on 36 trees distributed on three stands with different plantation densities in Padana Plane, Italy. The estimated profiles display high variability with an average of 1.6 cm of lateral compression. The results from this study demonstrate the potential and feasibility of estimating bole eccentricity by TLS, providing preliminary tools that will hopefully favor the diffusion of taper functions in operational environments. © 2019 by the authors.

Cities are responsible for more than 80% of global greenhouse gas emissions. Sequestration of air pollutants is one of the main ecosystem services that urban forests provide to the citizens. The atmospheric concentration of several pollutants such as carbon dioxide (CO2), tropospheric ozone (O3), and particulate matter (PM) can be reduced by urban trees through processes of adsorption and deposition. We predict the quantity of CO2, O3, and PM removed by urban tree species with the multilayer canopy model AIRTREE in two representative urban parks in Italy: Park of Castel di Guido, a 3673 ha reforested area located northwest of Rome, and Park of Valentino, a 42 ha urban park in downtown Turin. We estimated a total annual removal of 1005 and 500 kg of carbon per hectare, 8.1 and 1.42 kg of ozone per hectare, and 8.4 and 8 kg of PM10 per hectare. We highlighted differences in pollutant sequestration between urban areas and between species, shedding light on the importance to perform extensive in situ measurements and modeling analysis of tree characteristics to provide realistic estimates of urban parks to deliver ecosystem services. ©

Monitoring vegetation structure and functioning is critical for modelling terrestrial ecosystems and energy cycles. Leaf area index (LAI) is an important structural property of vegetation used in many land-surface, climate, and forest monitoring applications. Remote sensing provides a unique way to obtain estimates of leaf area index at spatially extensive areas. However, the analysis and extraction of quantitative information from remotely-sensed data require accurate cross-calibration with in situ forest measurements, which are generally spatially-and temporally-limited, thereby limiting the ability to compare the seasonal dynamic patterns between field and remotely-sensed time series. This is particularly relevant in temperate broadleaved forests, which are characterized by high level of complexity, which can complicate the retrieval of vegetation attributes from remotely-sensed data. In this study, we performed a long-term comparison of MODIS LAI products with continuous in situ leaf area index measurements collected monthly in temperate and Mediterranean forests from 2000 to 2016. Results indicated that LAI showed a good correlation between satellite and ground data for most of the stands, and the pattern in seasonal changes were highly overlapping between the time-series. We conclude that MODIS LAI data are suitable for phenological application and for up-scaling LAI from the stand level to larger scales. © 2019, Italian Society of Remote Sensing. All rights reserved.