Search

Quantifying 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.

European forest monitoring is a central topic nowadays due to the critical role that forests can play in combatting climate change. Crucial information on forests is the number of tree species and the area covered by each of them, as they vary concerning growth rates, wood value, value for biodiversity conservation, and susceptibility to disturbances and global warming. The primary source of forest information is national forest inventories (NFIs). However, they are updated too infrequently to accommodate climate change-related analyses, and their estimates are not based on wall-to-wall information. Remotely sensed data offer new opportunities for up-to-date and large-scale forest monitoring and for enhancing NFI estimates. However, despite the huge scientific efforts, it is still challenging to accurately map forest species through satellite imagery analysis.This study introduces a method for large-scale forest species mapping in the Netherlands using Sentinel-2 (S2) harmonic predictors and demonstrates a scientific procedure for reliably estimating area proportions from remote sensing-based species maps and comparing these estimates with NFI-based estimates.Compared to more standard predictors, harmonic predictors increased the model performance by 8% in terms of overall accuracy and the kappa coefficient by 9% while reducing omission and commission errors by as much as 18% and 13%, respectively. We estimated the area proportion of forest species for each 10-km cell covering the Netherlands first using NFI data and then using the predicted maps. Although the resulting estimates differ by source data and methods, we found an average deviation between NFI and remote sensing-based area proportion estimates of 9%, with deviations approaching 0% when increasing the number of NFI plots per cell.The outcomes of this research play an important role in understanding the relative strengths and limitations of remote sensing-based products and NFI data, as well as be a solid basis for forest species area proportion estimation when (i) no field data are available, (ii) more frequently updated estimates are required, or (iii) wall-to-wall and fine resolution spatially explicit estimates are needed.

A recently proposed modelling strategy predicts the net primary production (NPP) of forest ecosystems by combining the outputs of a NDVI-driven model, Modified C-Fix, and a bio-geochemical model, BIOME-BGC. This combination strategy takes into account the effects of forest disturbances but still assumes the presence of a mixture of differently aged trees. The application of this strategy to even-aged forests, therefore, requires a methodological advancement aimed at properly modifying the modelling of main ecosystem processes. In particular, the adaptation of the method to even-aged forests is based on the use of high-spatial-resolution airborne laser scanning (ALS) datasets, which yields green and woody biomass estimates that regulate the simulation of photosynthetic and respiratory processes, respectively. This approach was experimented in a Mediterranean study area, San Rossore Regional Park (Central Italy), which is covered by even-aged pine stands in different development phases. The modelling strategy is driven by MODIS NDVI images and meteorological data across five years (2011–2015), which are combined with estimates of forest canopy cover and height obtained from ALS data taken in 2015. This allows the production of stand NPP estimates, which, when converted into respective current annual increment (CAI) values, reasonably reproduce the age dependency of the available ground observations. The CAI estimates also show a highly significant correlation with these observations (r = 0.773) and moderate error levels (RMSE = 2.03 m3 ha−1 year−1, MBE = −0.45 m3 ha−1 year−1). These results confirm the potential of the modified simulation method to yield accurate high-spatial-resolution NPP estimates, which can offer valuable insights into C cycling and storage, in even-aged forests.

Remote sensing products are typically assessed using a single accuracy estimate for the entire map, despite significant variations in accuracy across different map areas or classes. Estimating per-pixel uncertainty is a major challenge for enhancing the usability and potential of remote sensing products. This paper introduces the dataDriven open access tool, a novel statistical design-based approach that specifically addresses this issue by estimating per-pixel uncertainty through a bootstrap resampling procedure. Leveraging Sentinel-2 remote sensing data as auxiliary information, the capabilities of the Google Earth Engine cloud computing platform, and the R programming language, dataDriven can be applied in any world region and variables of interest. In this study, the dataDriven tool was tested in the Rincine forest estate study area—eastern Tuscany, Italy—focusing on volume density as the variable of interest. The average volume density was 0.042, corresponding to 420 m3 per hectare. The estimated pixel errors ranged between 93 m3 and 979 m3 per hectare and were 285 m3 per hectare on average. The ability to produce error estimates for each pixel in the map is a novel aspect in the context of the current advances in remote sensing and forest monitoring and assessment. It constitutes a significant support in forest management applications and also a powerful communication tool since it informs users about areas where map estimates are unreliable, at the same time highlighting the areas where the information provided via the map is more trustworthy. In light of this, the dataDriven tool aims to support researchers and practitioners in the spatially exhaustive use of remote sensing-derived products and map validation.

Accurate structural information about forests, including canopy heights and diameters, is crucial for quantifying tree volume, biomass, and carbon stocks, enabling effective forest ecosystem management, particularly in response to changing environmental conditions. Since late 2018, NASA’s Global Ecosystem Dynamics Investigation (GEDI) mission has monitored global canopy structure using a satellite Light Detection and Ranging (LiDAR) instrument. While GEDI has collected billions of LiDAR shots across a near-global range (between 51.6°N and >51.6°S), their spatial distribution remains dispersed, posing challenges for achieving complete forest coverage. This study proposes and evaluates an approach that generates high-resolution canopy height maps by integrating GEDI data with Sentinel-1, Sentinel-2, and topographical ancillary data through three machine learning (ML) algorithms: random forests (RF), gradient tree boost (GB), and classification and regression trees (CART). To achieve this, the secondary aims included the following: (1) to assess the performance of three ML algorithms, RF, GB, and CART, in predicting canopy heights, (2) to evaluate the performance of our canopy height maps using reference canopy height from canopy height models (CHMs), and (3) to compare our canopy height maps with other two existing canopy height maps. RF and GB were the top-performing algorithms, achieving the best 13.32% and 16% root mean squared error for broadleaf and coniferous forests, respectively. Validation of the proposed approach revealed that the 100th and 98th percentile, followed by the average of the 75th, 90th, 95th, and 100th percentiles (AVG), were the most accurate GEDI metrics for predicting real canopy heights. Comparisons between predicted and reference CHMs demonstrated accurate predictions for coniferous stands (R-squared = 0.45, RMSE = 29.16%).

Future agricultural systems should increase productivity and sustainability of food production and supply. For this, integrated and efficient capture, management, sharing, and use of agricultural and environmental data from multiple sources is essential. However, there are challenges to understand and efficiently use different types of agricultural and environmental data from multiple sources, which differ in format and time interval. In this regard, the role of emerging technologies is considered to be significant for integrated data gathering, analyses and efficient use. In this study, a concept was developed to facilitate the full integration of digital technologies to enhance future smart and sustainable agricultural systems. The concept has been developed based on the results of a literature review and diverse experiences and expertise which enabled the identification of stat-of-the-art smart technologies, challenges and knowledge gaps. The features of the proposed solution include: data collection methodologies using smart digital tools; platforms for data handling and sharing; application of Artificial Intelligent for data integration and analysis; edge and cloud computing; application of Blockchain, decision support system; and a governance and data security system. The study identified the potential positive implications i.e. the implementation of the concept could increase data value, farm productivity, effectiveness in monitoring of farm operations and decision making, and provide innovative farm business models. The concept could contribute to an overall increase in the competitiveness, sustainability, and resilience of the agricultural sector as well as digital transformation in agriculture and rural areas. This study also provided future research direction in relation to the proposed concept. The results will benefit researchers, practitioners, developers of smart tools, and policy makers supporting the transition to smarter and more sustainable agriculture systems.

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.

Background: Gradients in local environmental characteristics may favour the abundance of species with particular traits, while other species decline, or favour species with different traits at the same time, without an increase in average species abundances. Therefore, we asked: do variations in species and traits differ along gradients of deadwood variables? Do species abundance and trait occurrence change with species richness within or between functional groups? Thus, we analysed the beetle assemblages of five forest sites located in Italy, along the Apennines mountains. Methods: From 2012 to 2018 we sampled beetles and five deadwood types in 193 plots to characterise the deadwood gradient: standing dead trees, snags, dead downed trees, coarse woody debris, and stumps. We modelled beetle species relative abundances and trophic traits occurrences against the deadwood variables using joint species distribution models. Results: Out of 462 species, only 77 showed significant responses to at least one deadwood type, with a weak mean response across species. Trophic groups showed mostly negative responses to deadwood variables. Species abundance increased with species richness among sites only for phytophagous and saproxylophagous. Trait occurrence did not increase with species richness among sites, except for phytophagous and saproxylophagous. However, trait occurrence changed significantly with species richness of several trophic groups within some sites. We found that increases in species richness do not result in decreases in species abundance of a given trophic group, but rather null or positive relationships were found suggesting low interspecific competition. Conclusions: Our findings suggest that in Mediterranean mountain forests there is still room for increasing the level of naturalness, at least for what concerns deadwood management. On one side, our findings suggest that competition for deadwood substrates is still low, on the other side they indicate that increasing deadwood volume and types to improve overall beetle richness may increase also beetle abundances.

Forests cover 30% of the Earth's landmass, host 80% of the biodiversity on land, and represent one of the main sinks of carbon. Studying forest ecosystems and dynamics is more crucial than ever now that the climate is changing. On the other hand, forest structural attributes and microhabitats data acquisition is challenging, and require huge efforts.Here we provide a georeferenced dataset of living trees, deadwood, and microhabitats referring to 199 plots (13 m radius), collected between 2012 and 2018, and located over six Apennine mountainous forest types across Italy. The dataset we provide promotes collaboration among researchers and improves the possibilities to analyze the evolution of forest ecosystems.

Abstract Worldwide, tropospheric ozone (O3) is a potential threat to wood production, but our understanding of O3 economic impacts on forests is still limited. To overcome this issue, we developed an approach for integrating O3 risk modelling and economic estimates, by using the Italian forests as a case study. Results suggested a significant impact of O3 expressed in terms of stomatal flux with an hourly threshold of uptake (Y = 1 nmol O3 m−2 leaf area s−1 to represent the detoxification capacity of trees), i.e. POD1. In 2005, the annual POD1 averaged over Italy was 20.4 mmol m−2 and the consequent potential damage ranged from 790.90 M€ to 2.85 B€ of capital value (i.e. 255–869 € ha−1, on average) depending on the interest rate. The annual damage ranged from 31.6 to 57.1 M€ (i.e. 10–17 € ha−1 per year, on average). There was also a 1.1% reduction in the profitable forest areas, i.e. with a positive Forest Expectation Value (FEV), with significant declines of the annual national wood production of firewood (− 7.5%), timber pole (− 7.4%), roundwood (− 5.0%) and paper mill (− 4.8%). Results were significantly different in the different Italian regions. We recommend our combined approach for further studies under different economic and phytoclimatic conditions.

Forests absorb 30% of human emissions associated with fossil fuel burning. For this reason, forest disturbances monitoring is needed for assessing greenhouse gas balance. However, in several countries, the information regarding the spatio-temporal distribution of forest disturbances is missing. Remote sensing data and the new Sentinel-2 satellite missions, in particular, represent a game-changer in this topic.Here we provide a spatially explicit dataset (10-meters resolution) of Italian forest disturbances and magnitude from 2017 to 2020 constructed using Sentinel-2 level-1C imagery and exploiting the Google Earth Engine GEE implementation of the 3I3D algorithm. For each year between 2017 and 2020, we provide three datasets: (i) a magnitude of the change map (between 0 and 255), (ii) a categorical map of forest disturbances, and (iii) a categorical map obtained by stratification of the previous maps that can be used to estimate the areas of several different forest disturbances. The data we provide represent the state-of-the-art for Mediterranean ecosystems in terms of omission and commission errors, they support greenhouse gas balance, forest sustainability assessment, and decision-makers forest managing, they help forest companies to monitor forest harvestings activity over space and time, and, supported by reference data, can be used to obtain the national estimates of forest harvestings and disturbances that Italy is called upon to provide.

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.

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.

Afforestation processes, natural and anthropogenic, involve the conversion of other land uses to forest, and they represent one of the most important land use transformations, influencing numerous ecosystem services. Although remotely sensed data are commonly used to monitor forest disturbance, only a few reported studies have used these data to monitor afforestation. The objectives of this study were two fold: (1) to develop and illustrate a method that exploits the 1985–2019 Landsat time series for predicting afforestation areas at 30 m resolution at the national scale, and (2) to estimate afforestation areas statistically rigorously within Italian administrative regions and land elevation classes. We used a Landsat best-available-pixel time series (1985–2019) to calculate a set of temporal predictors that, together with the random forests prediction technique, facilitated construction of a map of afforested areas in Italy. Then, the map was used to guide selection of an estimation sample dataset which, after a complex photointerpretation phase, was used to estimate afforestation areas and associated confidence intervals. The classification approach achieved an accuracy of 87%. At the national level, the afforestation area between 1985 and 2019 covered 2.8 ± 0.2 million ha, corresponding to a potential C-sequestration of 200 million t. The administrative region with the largest afforested area was Sardinia, with 260,670 ± 58,522 ha, while the smallest area of 28,644 ± 12,114 ha was in Valle d’Aosta. Considering elevation classes of 200 m, the greatest afforestation area was between 400 and 600 m above sea level, where it was 549,497 ± 84,979 ha. Our results help to understand the afforestation process in Italy between 1985 and 2019 in relation to geographical location and altitude, and they could be the basis of further studies on the species composition of afforestation areas and land management conditions.

Forest condition, biodiversity, and ecosystem services are strongly interlinked. The biodiversity levels depend to a large extent on the integrity, health, and vitality of forests at the same time as losses of forest biodiversity lead to decreased forest productivity and sustainability. Under this conceptual framework, this study presents a methodology for mapping forest condition at European scale supporting the attainment of the 2020 Aichi Biodiversity Target 5 “the rate of loss of all natural habitats, including forests, is at least halved and where feasible brought close to zero, and degradation and fragmentation is significantly reduced” and the implementation of Sustainable Development Goals (SDG), as well as the EU forest strategy since the sustainable forest management is oriented to support the provision of forest services and to enhance the condition of biodiversity forests’ host.The work presents the developments of an operational indicator at European scale. This spatially explicit information on forest condition can be the baseline map with a 1 km resolution to monitor the state and changes of condition by exposition to pressures and threats. This condition indicator considers structural, functional, and compositional aspects of forest with relevance for health and vitality of species and habitats hosted by forest ecosystems.The methodology implemented used harmonized, published and open datasets. It provided confident results for the assessment of the condition within hemiboreal, temperate and alpine forests, showing the Carpathian, Dinaric Alps and Alps, among others, as hotspots with pre-dominantly good condition. The results were validated with data derived from the reporting for the EU Habitat Directive and explicit dataset on known primary forests in Europe. However, this method underestimated the forest condition in the Mediterranean and Boreal forest types due to data gaps, regional specific characteristics, and design limitations.This study illustrates an operational and transferable approach for addressing the assessment of ecosystem forest condition at European scale being considered as a support tool for European countries when mapping and assessing their national territory, as potential common approach to map forest ecosystems that allows for consistent aggregation and comparisons across scales.

The relationships between structural complexity, deadwood abundance, microhabitat type and species-diversity indicators are excellent tools to monitor biodiversity in forest ecosystems.In spite of their importance, correlations between structural traits and Coleoptera communities in Mediterranean mountain forests have only rarely been investigated. Consequently, the magnitude and direction of the relationships between forest traits and biodiversity indicators remain poorly understood. In this study, we analyzed whether biodiversity indices of saproxylic and non-saproxylic beetle communities could be influenced by stand structure, microhabitat type, and deadwood abundance in two protected beech forests located in the central and southern Apennines (namely Gran Sasso e Monti della Laga National Park, GSML, and Cilento, Vallo di Diano e Alburni National Park, CVDA). Standard measurements of forest structural traits and quantitative assessment of tree microhabitats and deadwood were carried out. Adult beetles were collected using window flight traps and emergence traps on decaying deadwood. The two beech forests were different in terms of both beetle communities and structural traits. A two-block partial least squares analysis 2B-PLS highlighted differences in biodiversity indices and structural traits between the two forest ecosystems. In GSML, we observed that biodiversity indices were positively correlated with the volume of coarse woody debris and the presence fungal infections, clefts into the sapwood, and woodpecker cavities, while more dominant beetle communities were found under denser canopy cover. In CVDA, Coleoptera abundance was positively correlated with the basal area and crown broken microhabitats. Our results point toward the relevance of ecological attributes in tracking changes in beetle biodiversity in specific forest contexts. In these protected Mediterranean mountain beech stands, in which the main forest management strategies have the primary objective of biodiversity conservation, we suggest to progressively increase the structural diversity and canopy dynamics, as well as the volume of coarse woody debris.

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.

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

The dataset reports data collected in 38 square (50 x 50m) 0.25 ha plots representative of poplar plantations in Lombardy Region (Northern Italy), which were used to calibrate optical information derived from unmanned aerial vehicle (UAV) and satellite (Sentinel-2) sensors. In each plot, the diameter at breast height was measured using a caliper; height, stem and crown volume of each tree were then derived from diameter using allometric equations developed in an independent study. Additional canopy attributes (foliage and crown cover, crown porosity, leaf area index) were derived in each plot from 12-20 optical images collected using digital cover photography (DCP). The collected data allows characterizing the assessment of structure of these plantations, along with their variation over the rotation time. Canopy and crown data also enable the evaluation of optimal rotation and tree spacing, as well as the relationship between stand and canopy structure. The raw datasets consist of 2,591 records (trees) associated with inventory measurements and 616 records (images) associated with optical canopy measurements. An R code was also provided to calculate plot-level attributes from raw data. Dataset and associated metadata are freely available at http://dx.doi.org/10.17632/ycr7w5pvkt.1.

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 (R2 = 0.87), compared to Pines (R2 = 0.61), and to Spruce-dominated plots (R2 = 0.37).Lidar data were employed to estimate above ground biomass, using an exponential regression model; overall the R2 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.

The worldwide increase in urban and industrial areas represents a challenge for urban green management, pollutant mitigation and environmental monitoring. We propose an analysis approach for the spatial and spatial-temporal distribution of pollutants in the environment through dendrochemistry techniques, in order to gauge the value of this discipline in urban ecosystem. The proposed analysis models can be useful to evaluate significant changes in space and time related to economic activities and to describe the “impacts” of adopted strategies, as demonstrated in the case study of the opening or closure of factories, and therefore to describe the cause-effect relation connected with human activities. Trees represent a key factor for urban planning, providing a wide variety of ecosystem services including in-depth environmental monitoring, which is one of the main elements to be included in a high quality urban design. The proposed approach aims at suggesting the dendrochemistry as a novel and feasible tool definable as a cost-saving ecosystem service in the urban context.

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.

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.

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.

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.

National forest inventories are a primary source of data for the assessment of forest resources and lastly more often biodiversity at national scales. The diversity of adopted sampling designs and measurements reduces the prospect for a reliable comparison of generated estimates. The ICP Forest dataset represents a unique opportunity for a standardized approach of forest estimates through Europe. This work aims to provide a distribution map of the mean deadwood volume in European forest. A total of 3243 ICP Forests plots were analysed and presented. The study area extends over 3,664,576 km2 interesting 19 countries. We observed that the highest percentage of plots show a deadwood volume lower than 50 m3 ha−1, with a few of forests attaining around the maximum of 300 m3 ha−1. Forests with more than 100 m3 ha−1 are concentrated in mountainous regions, central Europe and other regions, linked to high-forest management types, while coppices-derived forest systems (part of the Great Britain, Mediterranean region) show lower deadwood content. The map of deadwood volume on European Forests is of interests for scientists, land planners, forest managers and decision-makers, as a reference for further evaluation of changes, stratified sampling, ground reference for model validation, restoration and conservation purposes.

The application of Monteith’s models to estimate the gross primary production (GPP) of forest ecosystems is a common and efficient practice. A modified version of one of these parametric models, C-Fix, has been widely tested in Italy at daily to annual temporal scales. The current paper is a step forward demonstrating the capability of this model to reproduce also inter-annual forest GPP variations. To this aim, dendrochronological measurements collected in two study areas were utilized to validate the model. The first study area is a coastal Mediterranean pine forest and the second is a more temperate ecosystem dominated by the presence of European beech. The experiment was organized into the following main steps. Modified C-Fix was applied to simulate annual GPP over a 15-year period for both sites. Both the main driving factors used by the model and the estimated GPP were then analysed and assessed against the detrended tree-ring width series of the collected cores. The study yields indications on the factors most influential on GPP for the two forest ecosystems and supports the capability of the proposed modelling approach to reproduce relevant inter-annual GPP variations, which is useful for predicting the impacts of ongoing climatic changes.

The study of land cover and land use dynamics are fundamental to understanding the radical changes that human activity is causing locally and globally and to analyse the continuous metamorphosis of landscape. In Europe, the Copernicus Program offers numerous territorial monitoring tools to users and decision makers, such as Sentinel data. This research aims at developing and implementing a land cover mapping and change detection methodology through the classification of Copernicus Sentinel-1 and Sentinel-2 satellite data. The goal is to create a versatile and economically sustainable algorithm capable of rapidly processing large amounts of data, allowing the creation of national-scale products with high spatial resolution and update frequency for operational purposes. Great attention was paid to compatibility with the main activities planned in the near future at the national and European level. In this sense, a land cover classification system consistent with the European specifications of the EAGLE group has been adopted. The methodology involves the definition of distinct sets of decision rules for each of the land cover macro-classes and for the land cover change classes. The classification refers to pixels’ spectral and backscatter characteristics, exploiting the main multi-temporal indices while proposing two new ones: the NDCI to distinguish between broad-leaved and needle-leaved trees, and the Burned Index (BI) to identify burned areas. This activity allowed for the production of a land cover map for 2018 and the change detection related to forest disturbances and land consumption for 2017–2018, reaching an overall accuracy of 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.

Since 2010, the Italian Ministry of University and Research issued new evaluation protocols to select candidates for University professorships and assess the bibliometric productivity of Universities and Research Institutes based on bibliometric indicators, i.e. scientific paper and citation numbers and the h-index. Under this framework, the objective of this study was to quantify the bibliometric productivity of the Italian forest research community during the 2002-2012 period. We examined the following productivity parameters: (i) the bibliometric productivity under the Forestry subject category at the global level; (ii) compared the aggregated bibliometric productivity of Italian forest scientists with scientists from other countries; (iii) analyzed publication and citation temporal trends of Italian forest scientists and their international collaborations; and (iv) characterized productivity distribution among Italian forest scientists at different career levels. Results indicated the following: (i) the UK is the most efficient country based on the ratio between Gross Domestic Spending (GDS) on Research and Development (R&D) and bibliometric productivity under the Forestry subject category, followed by Italy; (ii) Italian forest scientist productivity exhibited a significant positive time trend, but was characterized by high inequality across authors; (iii) one-half of the Italian forest scientist publications were written in collaboration with foreign scientists; (iv) a strong relationship exists between bibliometric indicators calculated by WOS and SCOPUS, suggesting these two databases have the same potential to evaluate the forestry research community; and (v) self-citations did not significantly affect the rank of Italian forest scientists.

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.

Forest growing stock volume (GSV) is an important parameter in the context of forest resource management. National Forest Inventories (NFIs) are routinely used to estimate forest parameters, including GSV, for national or international reporting. Remotely sensed data are increasingly used as a source of auxiliary information for NFI data to improve the spatial precision of forest parameter estimates. In this study, we combine data from the NFI in Poland with satellite images of Landsat 7 and 3D point clouds collected with airborne laser scanning (ALS) technology to develop predictive models of GSV. We applied an area-based approach using 13,323 sample plots measured within the second cycle of the NFI in Poland (2010–2014) with poor positional accuracy from several to 15 m. Four different predictive approaches were evaluated: multiple linear regression, k-Nearest Neighbours, Random Forest and Deep Learning fully connected neural network. For each of these predictive methods, three sets of predictors were tested: ALS-derived, Landsat-derived and a combination of both. The developed models were validated at the stand level using field measurements from 360 reference forest stands. The best accuracy (RMSE% = 24.2%) and lowest systematic error (bias% = −2.2%) were obtained with a deep learning approach when both ALS- and Landsat-derived predictors were used. However, the differences between the evaluated predictive approaches were marginal when using the same set of predictor variables. Only a slight increase in model performance was observed when adding the Landsat-derived predictors to the ALS-derived ones. The obtained results showed that GSV can be predicted at the stand level with relatively low bias and reasonable accuracy for coniferous species, even using field sample plots with poor positional accuracy for model development. Our findings are especially important in the context of GSV prediction in areas where NFI data are available but the collection of accurate positions of field plots is not possible or justified because of economic reasons.

Protection and recovery of natural resource and biodiversity requires accurate monitoring at multiple scales. Airborne Laser Scanning (ALS) provides high-resolution imagery that is valuable for monitoring structural changes to vegetation, providing a reliable reference for ecological analyses and comparison purposes, especially if used in conjunction with other remote-sensing and field products. However, the potential of ALS data has not been fully exploited, due to limits in data availability and validation. To bridge this gap, the global network for airborne laser scanner data (GlobALS) has been established as a worldwide network of ALS data providers that aims at linking those interested in research and applications related to natural resources and biodiversity monitoring. The network does not collect data itself but collects metadata and facilitates networking and collaborative research amongst the end-users and data providers. This letter describes this facility, with the aim of broadening participation in GlobALS.

This paper aims to provide general considerations, in the form of a scientific review, with reference to selected experiences of ALS applications under alpine, temperate and Mediterranean environments in Italy as case studies. In Italy, the use of ALS data have been mainly focused on the stratification of forest stands and the estimation of their timber volume and biomass at local scale. Potential for ALS data exploitation concerns their integration in forest inventories on large territories, their usage for silvicultural systems detection and their use for the estimation of fuel load in forest and pre-forest stands. Multitemporal ALS may even be suitable to support the assessment of current annual volume increment and the harvesting rates.

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.

A ten-year data-set descriptive of Italian forest gross primary production (GPP) has been recently constructed by the application of Modified C-Fix, a parametric model driven by remote sensing and ancillary data. That data-set is currently being used to develop multivariate regression models which link the inter-year GPP variations of five forest types (white fir, beech, chestnut, deciduous and evergreen oaks) to seasonal values of temperature and precipitation. The five models obtained, which explain from 52% to 88% of the inter-year GPP variability, are then applied to predict the effects of expected environmental changes (+2 °C and increased CO2 concentration). The results show a variable response of forest GPP to the simulated climate change, depending on the main ecosystem features. In contrast, the effects of increasing CO2 concentration are always positive and similar to those given by a combination of the two environmental factors. These findings are analyzed with reference to previous studies on the subject, particularly concerning Mediterranean environments. The analysis confirms the plausibility of the scenarios obtained, which can cast light on the important issue of forest carbon pool variations under expected global changes.

A large part of arid areas in tropical and sub-tropical regions are dominated by sparse xerophytic vegetation, which are essential for providing products and services for local populations. While a large number of researches already exist for the derivation of wall-to-wall estimations of above ground biomass (AGB) with remotely sensed data, only a few of them are based on the direct use of non-photogrammetric aerial photography. In this contribution we present an experiment carried out in a study area located in the Santiago Island in the Cape Verde archipelago where a National Forest Inventory (NFI) was recently carried out together with a new acquisition of a visible high-resolution aerial orthophotography. We contrasted two approaches: single-tree, based on the automatic delineation of tree canopies; and area-based, on the basis of an automatic image classification. Using 184 field plots collected for the NFI we created parametric models to predict AGB on the basis of the crown projection area (CPA) estimated from the two approaches. Both the methods produced similar root mean square errors (RMSE) at pixel level 45% for the single-tree and 42% for the area-based. However, the latest was able to better predict the AGB along all the variable range, limiting the saturation problem which is evident when the CPA tends to reach the full coverage of the field plots. These findings demonstrate that in regions dominated by sparse vegetation, a simple aerial orthophoto can be used to successfully create AGB wall-to-wall predictions. The level of these estimations’ uncertainty permits the derivation of small area estimations useful for supporting a more correct implementation of sustainable management practices of wood resources.