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.

National monitoring of forestlands and the processes causing canopy cover loss, be they abrupt or gradual, partial or stand clearing, temporary (disturbance) or persisting (deforestation), are necessary at fine scales to inform management, science and policy. This study utilizes the Landsat archive and an ensemble of disturbance algorithms to produce maps attributing event type and timing to >258 million ha of contiguous Unites States forested ecosystems (1986–2010). Nationally, 75.95 million forest ha (759,531 km²) experienced change, with 80.6% attributed to removals, 12.4% to wildfire, 4.7% to stress and 2.2% to conversion. Between regions, the relative amounts and rates of removals, wildfire, stress and conversion varied substantially. The removal class had 82.3% (0.01 S.E.) user’s and 72.2% (0.02 S.E.) producer’s accuracy. A survey of available national attribution datasets, from the data user’s perspective, of scale, relevant processes and ecological depth suggests knowledge gaps remain.

Accurately monitoring forest canopy height is crucial for sustainable forest management, particularly in southeastern North Carolina, USA, where dense forests and limited accessibility pose substantial challenges. This study presents an explainable machine learning framework that integrates sparse GEDI LiDAR samples with multi-sensor remote sensing data to improve both the accuracy and interpretability of forest canopy height estimation. This framework incorporates multitemporal optical observations from Sentinel-2; C-band backscatter and InSAR coherence from Sentinel-1; quad-polarization L-Band backscatter and polarimetric decompositions from the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR); texture features from the National Agriculture Imagery Program (NAIP) aerial photography; and topographic data derived from an airborne LiDAR-based digital elevation model. We evaluated four machine learning algorithms, K-nearest neighbors (KNN), random forest (RF), support vector machine (SVM), and eXtreme gradient boosting (XGB), and found consistent accuracy across all models. Our evaluation highlights our method’s robustness, evidenced by closely matched R² and RMSE values across models: KNN (R² of 0.496, RMSE of 5.13 m), RF (R² of 0.510, RMSE of 5.06 m), SVM (R² of 0.544, RMSE of 4.88 m), and XGB (R² of 0.548, RMSE of 4.85 m). The integration of comprehensive feature sets, as opposed to subsets, yielded better results, underscoring the value of using multisource remotely sensed data. Crucially, SHapley Additive exPlanations (SHAP) revealed the multi-seasonal red-edge spectral bands of Sentinel-2 as dominant predictors across models, while volume scattering from UAVSAR emerged as a key driver in tree-based algorithms. This study underscores the complementary nature of multi-sensor data and highlights the interpretability of our models. By offering spatially continuous, high-quality canopy height estimates, this cost-effective, data-driven approach advances large-scale forest management and environmental monitoring, paving the way for improved decision-making and conservation strategies.

The U.S. Forest Service, Forest Inventory and Analysis (FIA) program is tasked with making and reporting estimates of various forest attributes using a design-based network of permanent sampling plots. To make its estimates more precise, FIA uses a technique known as post-stratification to group plots into more homogenous classes, which helps lower variance when deriving population means. Currently FIA uses a nationally available map of tree canopy cover for post-stratification, which tends to work well for forest area estimates but less so for structural attributes like volume. Here we explore the use of new statewide digital aerial photogrammetric (DAP) point clouds developed from stereo imagery collected by the National Agricultural Imagery Program (NAIP) to improve these estimates in the southeastern mixed hardwood forests of Tennessee and Virginia, United States (U.S.). Our objectives are to 1. evaluate the relative quality of NAIP DAP point clouds using airborne LiDAR and FIA tree height measurements, and 2. assess the ability of NAIP digital height models (DHMs) to improve operational forest inventory estimates above the gains already achieved from FIA’s current post-stratification approach. Our results show the NAIP point clouds were moderately to strongly correlated with FIA field measured maximum tree heights (average Pearson’s r = 0.74) with a slight negative bias (−1.56 m) and an RMSE error of ~4.0 m. The NAIP point cloud heights were also more accurate for softwoods (R²s = 0.60–0.79) than hardwoods (R²s = 0.33–0.50) with an error structure that was consistent across multiple years of FIA measurements. Several factors served to degrade the relationship between the NAIP point clouds and FIA data, including a lack of 3D points in areas of advanced hardwood senescence, spurious height values in deep shadows and imprecision of FIA plot locations (which were estimated to be off the true locations by +/− 8 m). Using NAIP strata maps for post-stratification yielded forest volume estimates that were 31% more precise on average than estimates stratified with tree canopy cover data. Combining NAIP DHMs with forest type information from national map products helped improve stratification performance, especially for softwoods. The monetary value of using NAIP height maps to post-stratify FIA survey unit total volume estimates was USD 1.8 million vs. the costs of installing more field plots to achieve similar precision gains. Overall, our results show the benefit and growing feasibility of using NAIP point clouds to improve FIA’s operational forest inventory estimates.

Throughout the last three decades, north central Georgia has experienced significant loss in forest land and tree cover. This study revealed the temporal patterns and thematic transitions associated with this loss by augmenting traditional forest inventory data with remotely sensed observations. In the US, there is a network of field plots measured consistently through time from the USDA Forest Service’s Forest Inventory and Analysis (FIA) Program, serial photo-based observations collected through image-based change estimation (ICE) methodology, and historical Landsat-based observations collected through TimeSync. The objective here was to evaluate how these three data sources could be used to best estimate land use and land cover (LULC) change. Using data collected in north central Georgia, we compared agreement between the three data sets, assessed the ability of each to yield adequately precise and temporally coherent estimates of land class status as well as detect net and transitional change, and we evaluated the effectiveness of using remotely sensed data in an auxiliary capacity to improve detection of statistically significant changes. With the exception of land cover from FIA plots, agreement between paired data sets for land use and cover was nearly 85%, and estimates of land class proportion were not significantly different for overlapping time intervals. Only the long time series of TimeSync data revealed significant change when conducting analyses over five-year intervals and aggregated land categories. Using ICE and TimeSync data through a two-phase estimator improved precision in estimates but did not achieve temporal coherence. We also show analytically that using auxiliary remotely sensed data for post-stratification for binary responses must be based on maps that are extremely accurate in order to see gains in precision. We conclude that, in order to report LULC trends in north central Georgia with adequate precision and temporal coherence, we need data collected on all the FIA plots each year over a long time series and broadly collapsed LULC classes.

In this paper, we propose a new estimator for creating expansion factors for survey plots in the US Forest Service (USFS) Forest Inventory and Analysis program. This estimator was previously used in the GIS literature, where it was called Penalized Maximum Entropy Dasymetric Modeling. We show here that the method is a regularized version of the raking estimator widely used in sample surveys. The regularized raking method differs from other predictive modeling methods for integrating survey and ancillary data, in that it produces a single set of expansion factors that can have a general purpose which can be used to produce small-area estimates and wall-to-wall maps of any plot characteristic. This method also differs from other more widely used survey techniques, such as GREG estimation, in that it is guaranteed to produce positive expansion factors. Here, we extend the previous method to include cross-validation, and provide a comparison to expansion factors between the regularized raking and ridge GREG survey calibration.