Information on land use and land cover (LULC) including forest cover is important for the development of strategies for land planning and management. Satellite remotely sensed data of varying resolutions have been an unmatched source of such information that can be used to produce estimates with a greater degree of confidence than traditional inventory estimates. However, use of these data has always been a challenge in tropical regions owing to the complexity of the biophysical environment, clouds, and haze, and atmospheric moisture content, all of which impede accurate LULC classification. We tested a parametric classifier (logistic regression) and three non-parametric machine learning classifiers (improved k-nearest neighbors, random forests, and support vector machine) for classification of multi-temporal Sentinel 2 satellite imagery into LULC categories in Dak Nong province, Vietnam. A total of 446 images, 235 from the year 2017 and 211 from the year 2018, were pre-processed to gain high quality images for mapping LULC in the 6516 km² study area. The Sentinel 2 images were tested and classified separately for four temporal periods: (i) dry season, (ii) rainy season, (iii) the entirety of the year 2017, and (iv) the combination of dry and rainy seasons. Eleven different LULC classes were discriminated of which five were forest classes. For each combination of temporal image set and classifier, a confusion matrix was constructed using independent reference data and pixel classifications, and the area on the ground of each class was estimated. For overall temporal periods and classifiers, overall accuracy ranged from 63.9% to 80.3%, and the Kappa coefficient ranged from 0.611 to 0.813. Area estimates for individual classes ranged from 70 km² (1% of the study area) to 2200 km² (34% of the study area) with greater uncertainties for smaller classes.

Globally, forests are a crucial natural resource, and their sound management is critical for human and ecosystem health and well-being. Efforts to manage forests depend upon reliable data on the status of and trends in forest resources. When these data come from well-designed natural resource monitoring (NRM) systems, decision makers can make science-informed decisions. National forest inventories (NFIs) are a cornerstone of NRM systems, but require capacity and skills to implement. Efficiencies can be gained by incorporating auxiliary information derived from remote sensing (RS) into ground-based forest inventories. However, it can be difficult for countries embarking on NFI development to choose among the various RS integration options, and to develop a harmonized vision of how NFI and RS data can work together to meet monitoring needs. The NFI of the United States, which has been conducted by the USDA Forest Service’s (USFS) Forest Inventory and Analysis (FIA) program for nearly a century, uses RS technology extensively. Here we review the history of the use of RS in FIA, beginning with general background on NFI, FIA, and sampling statistics, followed by a description of the evolution of RS technology usage, beginning with paper aerial photography and ending with present day applications and future directions. The goal of this review is to offer FIA’s experience with NFI-RS integration as a case study for other countries wishing to improve the efficiency of their NFI programs.

Background: Large area forest inventories often use regular grids (with a single random start) of sample locations to ensure a uniform sampling intensity across the space of the surveyed populations. A design-unbiased estimator of variance does not exist for this design. Oftentimes, a quasi-default estimator applicable to simple random sampling (SRS) is used, even if it carries with it the likely risk of overestimating the variance by a practically important margin. To better exploit the precision of systematic sampling we assess the performance of five estimators of variance, including the quasi default. In this study, simulated systematic sampling was applied to artificial populations with contrasting covariance structures and with or without linear trends. We compared the results obtained with theSRS, Matérn’s, successive difference replication, Ripley’s, and D’Orazio’s variance estimators.Results: The variances obtained with the four alternatives to theSRSestimator of variance were strongly correlated, and in all study settings consistently closer to the target design variance than the estimator forSRS. The latter always produced the greatest overestimation. In populations with a near zero spatial autocorrelation, all estimators, performed equally, and delivered estimates close to the actual design variance.Conclusion: Without a linear trend, theSDRandDORestimators were best with variance estimates more narrowly distributed around the benchmark; yet in terms of the least average absolute deviation, Matérn’s estimator held a narrow lead. With a strong or moderate linear trend, Matérn’s estimator is choice. In large populations, and a low sampling intensity, the performance of the investigated estimators becomes more similar.

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.

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 m³ ha⁻¹). 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.

Multisource forest inventory methods were developed to improve the precision of national forest inventory estimates. These methods rely on the combination of inventory data and auxiliary information correlated with forest attributes of interest. As these methods have been predominantly tested over coniferous forests, the present study used this approach for heterogeneous and complex deciduous forests in the center of France. The auxiliary data considered included a forest type map, Landsat 8 spectral bands and derived vegetation indexes, and 3D variables derived from photogrammetric canopy height models. On a subset area, changes in canopy height estimated from two successive photogrammetric models were also used. A model-assisted inference framework, using a k nearest-neighbors approach, was used to predict 11 field inventory variables simultaneously. The results showed that among the auxiliary variables tested, 3D metrics improved the precision of dendrometric estimates more than other auxiliary variables. Relative efficiencies (RE) varying from 2.15 for volume to 1.04 for stand density were obtained using all auxiliary variables. Canopy height changes also increased RE from 3% to 26%. Our results confirmed the importance of 3D metrics as auxiliary variables and demonstrated the value of canopy change variables for increasing the precision of estimates of forest structural attributes such as density and quadratic mean diameter.