Evaluating Statewide NAIP Photogrammetric Point Clouds for Operational Improvement of National Forest Inventory Estimates in Mixed Hardwood Forests of the Southeastern U.S.

: 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


Introduction 1.National Forest Inventory
Around the world, National Forest Inventory (NFI) programs collect information on the location, composition, and distribution of forest resources using a combination of field measurements and remote sensing data.This information, which is used to estimate forest characteristics at various points in time, is one of the key sources of data for forecasting the effects of climate change and anthropogenic greenhouse gas emissions [1].To help guide sustainable management, NFI programs often use design-based sampling to obtain unbiased estimates of timber resources in their respective countries [2][3][4].Traditionally, estimates have focused on basic attributes, such as forest volume and area.However, in recent years, NFI programs have evolved to provide information on a wider range of forest characteristics, from damage and disease to ecosystem services such as biodiversity and carbon pools [5,6].To meet the growing needs of local land managers and to improve inputs into international reporting (e.g., United Nations Intergovernmental Panel on Climate Change or IPCC), NFI programs employ techniques such as sample intensification [7], small area methods [8], model-assisted and hybrid-estimation [9], and post-stratification [10] to increase the scale and precision of statistical estimates derived from forest inventory data.

FIA and Nationwide Forest Attribute Estimation
In the U.S., the Forest Inventory and Analysis (FIA) program, managed by United States Forest Service (USFS), has been designated as the design-based NFI to determine the extent, condition, volume, and growth of forests, and monitor changes in forest health.FIA data are collected on a permanent network of fixed sampling plots set approximately every 2400 ha across the U.S. using a random, spatially balanced design [11].The tree (e.g., DBH, height, species) and stand-level attributes (e.g., basal area, trees per acre) collected by FIA have a wide range of applications [12], such as assessment of carbon pools [13], analysis of timber product outputs [14], calculation of forest growth trends [15,16], and study of species' range shifts in response to climate change [17][18][19].Although FIA is the most inclusive and longest running forest monitoring program in the U.S., its statistical estimates still contain uncertainty due to its relatively sparse sampling frame (which is designed to produce strategic-level estimates), number of unsampled plot locations resulting from denied access and hazardous conditions [20], and natural variability in forest composition and structure among sampled plot locations [21].
To increase the precision of forest inventory estimates without growing the number of sampling plots, a statistical technique known as post-stratification can be used.This involves aggregating measurements of forest attributes (response variables) into strata that are more homogeneous than the full dataset, typically through the use of wall-to-wall land cover maps and other sources of classified remote sensing data (e.g., see McRoberts et al., 2002).Reductions in estimate variance are obtained by deriving population means as the weighted mean of each stratum [22,23].The FIA program has long used post-stratification to increase the precision of its forest inventory estimates starting in the 1980s with double sampling with aerial photography [24].Over the years much research has been done on the topic [25][26][27], leading to the use of various auxiliary remote sensing datasets to form strata [28].Currently, all four FIA units use some form of the Landsat-based, National Land Cover Database (NLCD) Tree Canopy Cover (TCC) map [29] for post-stratification.

Use of Point Clouds in Forest Inventory Applications
Although the spectral information from optical satellite imagery provides greater precision of strata for canopy-level attributes such as the proportion of forest area, it is less effective for attributes related to below-canopy structure and quantity, such as aboveground biomass, because the spectral signal tends to saturate, making it ineffective at distinguishing between mature and old-growth forests [30][31][32].In contrast to optical remote sensing, estimating canopy-level forest attributes with airborne laser scanning (ALS) (or airborne LiDAR) data has been more successful [33,34].Recently, many countries have acquired nationwide ALS LiDAR to produce high-resolution topographic elevation data.As a result of this wall-to-wall availability, point cloud data from ALS has become a viable data source for statewide estimation of forest attributes [2,5,35].Despite the growing availability of ALS data in the public domain, there are still significant challenges to using it in operational reporting workflows within FIA.First, since most statewide LiDAR acquisitions are focused on deriving topographic elevation, most are acquired during leaf-off conditions, which limits the data's ability to accurately resolve canopy cover, height, and other structural attributes in hardwood-dominated forests.Second, since the topography is relatively stable and ALS is still expensive to acquire, most states and countries do not have plans for repeat collections, further limiting temporal availability of ALS data for ongoing, broadscale forest monitoring programs.
An affordable alternative for wall-to-wall forest structure data at the state level are point clouds derived from Digital Aerial Photogrammetry (DAP; [36]).Using photogrammetric processing techniques (also known as Structure from Motion, [37]), DAP can be used to produce 3D point clouds from overlapping stereo images.Unlike ALS, which uses an active sensor to penetrate the canopy to obtain information about understory vegetation structure and ground surface topography, DAP point clouds only capture the top surface of 3D objects, such as trees and buildings.Because DAP does not fully resolve the ground surface, a high-resolution Digital Terrain Model (DTM) is needed to obtain normalized surface heights, referred to here as a Digital Height Model (DHM).Here we differentiate a DHM, which represents a rasterized surface of normalized heights for all 3D objects on the landscape including trees, vegetation, and non-forest structures such as houses, buildings, and, in their absence, ground and pavement, from a Digital Surface Model (DSM), a commonly used term for a rasterized surface of non-normalized elevations for all 3D objects; and a Canopy Height Model (CHM), a raster of normalized heights for trees and forests, where heights for non-forest objects have been removed, typically using a mask based on land cover data derived from classified satellite imagery [38].
Although ALS captures more within canopy detail in forests and is the primary source for high-resolution DTMs, DAP offers certain advantages to operational forest inventory programs.For example, DAP point clouds are approximately 1/2 to 1/3 the cost of ALS acquisitions [39,40] and in many cases, the stereo images required to produce them are becoming increasingly available from continuous monitoring programs.For example, in the U.S., the National Agriculture Imagery Program (NAIP) collects new, high-resolution stereo imagery over the conterminous U.S. every 2-3 years, and the acquisition cost of its DSM product is less than $0.01 USD/ha [27].Thus, point clouds derived from NAIP DAP (NAIP point cloud hereafter) are an affordable and repeatable source of high-resolution forest structure data for all states in the conterminous U.S.While shorter update intervals are preferred for tracking changes in areas where forests are disturbed more frequently, the potential 2-3-year return interval for NAIP DAP is less than half the 10-year shelf-life of most ALS datasets used for model-assisted estimation [41].In addition, recent studies performed using small and heterogeneous datasets showed that the precision of height estimates from DAP point clouds is comparable to, albeit slightly lower than, the height estimates from ALS [42].For example, ref. [43] reported that DAP and ALS metrics showed similar levels of correlation in terms of forest structural attributes such as height, volume, and biomass.Slope, canopy cover, and image acquisition date can generate significant differences between the height metrics derived from DAP and ALS point clouds [44].When considering combinations of several biophysical forest characteristics, metrics derived from ALS tend to perform better than those from DAP [45].Other factors which limit the accuracy of DAP DHMs include the selected DTM cell resolution (e.g., 1 m vs. 10 m DEM), the geo-precision of plot coordinate locations (e.g., use of recreational grade vs. survey-grade, HPGNSS receivers, see [27]), forest type, stand density, and topographic complexity (e.g., slope, aspect, and elevation) [44].Despite being slightly less accurate, DAP point clouds have been successfully used in several nationwide applications to estimate forest canopy height [46,47], and are thus a viable and affordable option to consider for use in operational forest inventory programs.

Using NAIP Digital Height Surfaces to Improve FIA Forest Inventory Estimates
Given its ability to enable large-scale, frequent, wall-to-wall assessments of forest structure, there is growing interest in learning more about the quality of NAIP DAP DHMs (NAIP DHM hereafter) and their potential for improving operational forest inventory estimates derived by the FIA program.In the Pacific Northwest region of the U.S., ref. [27] found that using NAIP DHMs for post-stratification reduced the variance of mean forest volume estimates by as much as four times compared to estimates derived solely from simple random sampling.This success, along with other studies indicating the potential for using DAP-derived DHMs to improve estimates of height, density, and volume [48,49], suggest they may also be useful for improving traditional FIA estimates developed for statewide reporting.While the study by [27] directly points to the potential benefits of using NAIP DHMs in the context of FIA estimation, it does so with a limited number of sample plots (only 1/10 of the FIA plots in WA state were used) in a predominantly coniferbased forest system.To date, the effectiveness of NAIP DHMs has not been evaluated in southeastern mixed hardwood forests, where seasonality (e.g., fall leaf senescence), compositional heterogeneity [50], and disturbance rates are higher than in other regions of the country [51,52].

Objectives
Here, our study aims to evaluate NAIP DHMs in an operational context, using thousands of FIA plots to quantify their accuracy in resolving tree heights, as well as to assess their ability to improve FIA estimates in diverse, southern mixed hardwood forests (Figure 1).To achieve these objectives, we use FIA field-measured trees and ALS LiDAR data to validate the accuracy of height estimates derived from NAIP DHMs collected in 2018 for the states of Tennessee (TN) and Virginia (VA).Using a series of scatter plots and linear regression models, we evaluate the effect of forest type (hardwood vs. softwood), resolution of rasterized DHM (1 m vs. 10 m) and the timing of FIA field measurement on the level of variability and bias found in the three sources of height information (FIA field measurements, ALS LiDAR, and NAIP point clouds).Additionally, we also assess the utility of NAIP DHMs when used for post-stratification of FIA forest area, total volume (in m 3 ), and volume per area (in m 3 /acre) estimates.Using 10 m wall-to-wall NAIP DHMs, we derive three different strata maps, which are independently used to post-stratify population totals and species group estimates derived with FIA's standard estimation procedure (described in chapter 4 of [11]).The three strata maps include: 1. the full DHM, which contains height values for both forest and non-forest areas; 2. a CHM, where heights in non-forest areas have been removed using a forest mask derived from a national Landsat-based disturbance map [52]; and 3. a CHM with forest type information (e.g., deciduous, evergreen, and mixed) added from NLCD land cover data.To test the utility of each strata map, we produced two measures of relative efficiency (RE), one comparing the variance from each of the three NAIP stratification maps to the variance from FIA's current post-stratification approach (referred to as RE FIA ), and another comparing the variance from the best strata map to the variance obtained from simple random sampling (referred to as RE SRS ).Furthermore, we also examine the utility of NAIP DHMs to improve general forest estimate totals for specific hardwood and softwood species groups.We conclude our analysis by making general recommendations for future use of NAIP DAP and derivative DHM/CHM products within FIA and the southeast region.

Study Area
Our study area encompasses the states of TN and VA in the southeastern U.S. (Figure 2A).According to FIA's 2009 state report [53], TN has a land area of 109,150 km 2 , 51.9% of which is covered by forest, with a total growing stock volume of 0.784 billion m 3 .Most of TN falls in the humid subtropical zone (Köppen climate classification Cfa), except for a few areas in the eastern part of the state that have more of an oceanic (Köppen Cfb) climate.For reporting purposes, FIA divides TN into five survey units: West (T1), West Central (T2), Central (T3), Plateau (T4), and East (T5) (Figure 2B).A large amount of the land area in the West survey unit is comprised of cropland; thus, forest cover is lower than in other survey units.West Central is a mixture of cropland, pasture, woodland, and oakhickory-pine forest.The Central survey unit includes Nashville, the state's largest metropolitan area.The East survey unit is located along the Appalachian Mountains.This unit's elevation difference is approximately 1500 m and Appalachian oak and northern hardwood forests are the dominant forest types.On the other hand, VA has an area of 102,280 km 2 , of which 58.7% was covered by forest in 2016, with a total growing stock volume of 1.116 billion m 3 [54].VA's climate is humid subtropical (Cfa).FIA divides VA into five survey units for reporting, which include: Southern Mountains (V1), Northern Mountains (V2), Southern Piedmont (V3), Northern Piedmont (V4), and Coastal Plain (V5) (Figure 2B).The Southern and Northern Mountains are located along the ridge of Appalachian Mountains; thus, the geographical and environmental settings of these two units are similar to the East survey unit in TN.The Northern and Southern Piedmont units represent the transitional area between inland mountains and the coastal plain, while the latter survey unit occupies the eastern border of the state.The natural vegetation cover in the coastal plain of VA is dominated by loblolly-shortleaf pine and oak-hickory-pine forest.A basic overview of the forest characteristics, topography, and disturbance rates (based on data from the USFS Landscape Change Monitoring System or LCMS) in each FIA survey unit can be found in Table 1.

Study Area
Our study area encompasses the states of TN and VA in the southeastern U.S. (Figure 2A).According to FIA's 2009 state report [53], TN has a land area of 109,150 km 2 , 51.9% of which is covered by forest, with a total growing stock volume of 0.784 billion m 3 .Most of TN falls in the humid subtropical zone (Köppen climate classification Cfa), except for a few areas in the eastern part of the state that have more of an oceanic (Köppen Cfb) climate.For reporting purposes, FIA divides TN into five survey units: West (T1), West Central (T2), Central (T3), Plateau (T4), and East (T5) (Figure 2B).A large amount of the land area in the West survey unit is comprised of cropland; thus, forest cover is lower than in other survey units.West Central is a mixture of cropland, pasture, woodland, and oak-hickory-pine forest.The Central survey unit includes Nashville, the state's largest metropolitan area.The East survey unit is located along the Appalachian Mountains.This unit's elevation difference is approximately 1500 m and Appalachian oak and northern hardwood forests are the dominant forest types.On the other hand, VA has an area of 102,280 km 2 , of which 58.7% was covered by forest in 2016, with a total growing stock volume of 1.116 billion m 3 [54].VA's climate is humid subtropical (Cfa).FIA divides VA into five survey units for reporting, which include: Southern Mountains (V1), Northern Mountains (V2), Southern Piedmont (V3), Northern Piedmont (V4), and Coastal Plain (V5) (Figure 2B).The Southern and Northern Mountains are located along the ridge of Appalachian Mountains; thus, the geographical and environmental settings of these two units are similar to the East survey unit in TN.The Northern and Southern Piedmont units represent the transitional area between inland mountains and the coastal plain, while the latter survey unit occupies the eastern border of the state.The natural vegetation cover in the coastal plain of VA is dominated by loblolly-shortleaf pine and oak-hickory-pine forest.A basic overview of the forest characteristics, topography, and disturbance rates (based on data from the USFS Landscape Change Monitoring System or LCMS) in each FIA survey unit can be found in Table 1.

Point Cloud Data
In 2018, FIA and other USFS partners purchased point clouds for the states of TN and VA (among others).In both states the NAIP imagery was acquired using multiple Leica ADS-100 SH100 pushbroom sensors resulting in the collection of 60 cm RGB digital imagery with approximately 30% side lap.Three look angles were collected (25.6 • forward, nadir, 19.4 • backward) of which only the forward lap was used to derive the stereo imagery.Point clouds were derived using the Leica XPro SGM Module and delivered in .laz1.4 format with accompanying color infrared (CIR) and RGB values encoded [55].Contract specifications required the nominal point spacing, or distance between the 3D points, to be double the ground surface distance of the native imagery (or in this case 1.2 m).Both the VA and TN point clouds exceeded these specifications with 0.7 m and 1.0 m point spacings, respectively.Acquisition dates ranged between early August and mid-December (Figure 3D and Table 2A); thus, a portion of the imagery was taken during the leaf-off season for hardwoods.Additionally, because the NAIP point clouds were collected in different seasons and processed by different vendors, there are noticeable differences in point density across the study area (Figure 3A and Table 2B).For example, when TN's NAIP point cloud was rasterized into a 10 m grid product, about 1% of the pixels had fewer than 30 points per 100 m 2 while in VA more than 4% of pixels fell in this range.Areas with few if any points were predominately the result of deep shadows caused by low sun angles and steep topography due to the late fall/early winter acquisition period (Figure 3A-C).
Remote Sens. 2022, 14, 4386 6 of 30   As part of our height validation, ALS LiDAR point clouds covering FIA field plot locations were acquired from the USGS 3D Elevation Program (3DEP, [56]).We acquired LiDAR data categorized as Quality Level 2 by the U.S. Geological Survey.Quality Level 2 products have a minimum nominal pulse spacing of 0.7 m and a vertical error of 10 cm, evaluated as the root mean square error (RMSE) in the elevation dimension; in addition, only the data acquired within the last eight years (2011-2019) were included in the Level 2 product [57].Because the 3DEP 1 m DEMs are derived from the ALS point clouds, their spatial coverage is nearly identical to the coverage shown in red shading in Figure 2A.Note, data comparisons involving ALS are limited to areas with available data.

Point Cloud Data
In 2018, FIA and other USFS partners purchased point clouds for the states of TN and VA (among others).In both states the NAIP imagery was acquired using multiple Leica ADS-100 SH100 pushbroom sensors resulting in the collection of 60 cm RGB digital imagery with approximately 30% side lap.Three look angles were collected (25.6° forward, nadir, 19.4° backward) of which only the forward lap was used to derive the stereo imagery.Point clouds were derived using the Leica XPro SGM Module and delivered in .laz1.4 format with accompanying color infrared (CIR) and RGB values encoded [55].Contract specifications required the nominal point spacing, or distance between the 3D points, to be double the ground surface distance of the native imagery (or in this case 1.2 m).Both the VA and TN point clouds exceeded these specifications with 0.7 m and 1.0 m point spacings, respectively.Acquisition dates ranged between early August and mid-December (Figure 3D and Table 2A); thus, a portion of the imagery was taken during the leaf-off season for hardwoods.Additionally, because the NAIP point clouds were collected in different seasons and processed by different vendors, there are noticeable differences in point density across the study area (Figure 3A and Table 2B).For example, when TN's NAIP point cloud was rasterized into a 10 m grid product, about 1% of the pixels had fewer than 30 points per 100 m 2 while in VA more than 4% of pixels fell in this range.Areas with few if any points were predominately the result of deep shadows caused by low sun angles and steep topography due to the late fall/early winter acquisition period (Figure 3A-C).

Digital Elevation Data
To convert the rasterized NAIP point clouds into normalized heights, we used digital terrain data to represent the ground surface.Here, 1 m and 1/3 arc-second (approximately 10 m) DEMs from USGS 3DEP were used to represent the ground elevation.Both DEM datasets were acquired by several agencies and processed to a common coordinate system and unit of vertical measure.While the 1 m DEM is available for nearly the entire state of TN, large portions of VA are missing 1 m data, including the mountainous highlands and coastal plains located in the west central and eastern parts of the state, respectively (Figure 2A).Unlike the 1 m data, the 10 m DEM is available seamlessly across both states.

NAIP DHMs and Strata Map Processing
To produce statewide, wall-to-wall DHMs we used the "gridmetrics" tool, part of the FUSION software package developed and distributed by the USFS Pacific Northwest Research Station [58].To avoid dealing with missing DEMs and the logistical constraints involved with processing and storing numerous high resolution raster datasets, we opted to produce 10 m DHMs for both states.To minimize impacts of missing data and outliers we used the rasterized point density map (Figure 2A) and the North American Forest Dynamics (NAFD) national disturbance map [59,60].For example, we used the NAFD map to identify forest pixels with five or fewer NAIP DAP points and/or maximum heights greater than 50 m (the height of the tallest measured tree in TN and VA) and assigned them the overall average maximum height value for all forest pixels, which was 18 m.Additionally, areas mapped as water in the NAFD map were assigned a maximum height value of 0.
After processing and filtering the DHMs, we combined them with other remote sensing datasets to form three strata maps for use in post-stratifying the FIA estimates.The three strata maps include: 1. the raw DHM, with surface heights for all objects on the landscape, including vegetation, as well as non-forest objects such as buildings, houses, and other manmade structures (referred to as PS DHM ); 2. a CHM where surface heights of non-forest objects have been removed using a forest/non-forest mask derived from the NAFD disturbance map ( [52]; referred to as PS CHM ); and 3. a CHM combined with general forest type information (evergreen, mixed, deciduous) obtained from the 2016 NLCD map ( [61,62] referred to as PS CHM+FT ).To develop the three strata maps we first binned the DHM into 8 roughly equal height classes (0-3 m; 3-7.5 m; 7.5-13.5 m; 13.5-18 m; 18-22.5 m; 22.5-27 m; 27-31.5 m; and 31.5-50m).Then to generate the PS DHM map we added a ninth class representing areas mapped as water in the NAFD map.Similarly, the PS CHM map was generated by adding a ninth and tenth class representing non-forest and water areas in the NAFD map.Lastly, for the PS CHM+FT map we used the NLCD land cover map to assign water and non-forest classes, as well as to further break out the 8 DHM height bins into deciduous (NLCD class 41), evergreen (NLCD class 42), mixed (NLCD class 43) and bottomland hardwood (NLCD class 90) classes.Note that when grouping strata, not all forest types contained pixels from all height bins, thus after iterative recoding, the TN and VA PS CHM+FT maps contained 14 and 16 forest specific height bins (e.g., deciduous 0-7.5 m; evergreen 0-7.5 m; mixed 0-7.5 m, etc.), plus non-forest and water classes, respectively.Moreover, while the NAFD and NLCD maps both have 30 m pixels, they are only used for masking and filtering the DHMs; thus, the final resolution of all three NAIP strata maps is 10 m.
To test the three NAIP strata maps against the precision gained by FIA's current poststratification approach, we also developed a fourth strata map by binning the NLCD TCC map [29] into four separate classes (0-10%; 11-47%; 48-84% and 85-100%).This fourth strata map (referred to as PS TCC ) has a resolution of 30 m and is only used to calculate the relative efficiency gained from using the NAIP height maps instead of tree canopy cover data to post-stratify the FIA estimates (i.e., RE FIA ).

FIA Data
Field data collected by the USFS FIA program in the states of TN and VA were obtained online from the FIA data mart [63].Here we use the FIA data at the plot-level to validate the NAIP point cloud heights and as a sample to evaluate the benefit of using NAIP DHMs for post-stratification.The design of the FIA program intends for each permanent plot in TN and VA to be measured every five years [11]; however, due to logistical constraints and other issues, it can take more than 5 years to collect the full cycle of plots.To derive estimates for the time period closest to 2018 (which is the year the NAIP imagery was acquired) we selected the most recent field measurement for each plot, resulting in 2931 forested plots in TN collected from 2012-2019 (with an average n = 366 plots per year) and 3304 plots in VA collected from 2014-2019 (with an average n = 550 plots per year).Each FIA plot is comprised of four, 7.32 m radius circular subplots which are 36.58m apart and oriented 120 • , 240 • , and 360 • from the plot center.FIA plot coordinates are recorded at the center of subplot 1, which falls in the middle of the 0.042 ha plot area.Here we use the actual FIA plot locations, which are not publicly available; however, since the locations of subplots 2-4 are not precisely known, we only use the FIA tree measurements collected on subplot 1 to validate the point clouds.Since some of the forested plots have no trees on subplot 1 only a subset of measurements from each state are used for the point cloud validation (1936 in TN and 2258 in VA).To determine the tree with the maximum height, we used the FIA tree variable ACTUALHT (m), which is a record of the distance from ground level to the highest remaining part of the tree still present and attached to the bole [64].Hardwood and softwood species groups are defined using the FIA tree variable SPGRPCD.In the eastern U.S. there are 19 hardwood and 9 softwood species groups (see Table S1 in the Supplementary Materials for individual species groups).Because species group is a tree-level attribute, some degree of simplification is necessary to group plots into similar classes.Here we use the species group of the tallest live tree to assign each subplot to the hardwood or softwood class.While this approach may introduce some errors, we believe they will be minor at the hardwood/softwood level of aggregation because in many cases in southern mixed hardwoods the tallest tree on a plot comes from the most dominant species group, either due to competition (in natural stands) or planting (in managed stands).Volume estimates are based on the FIA tree variable VOLCFGRS (m 3 ), which represents the total cubic volume of wood in the central stem of trees with diameter at breast height (DBH) ≥5.0 inches.Finally, for the post-stratification analysis, information from all 4 FIA subplots is used to derive forest area (in ha), total live volume (in m 3 ) and volume/hectare (in m 3 /ha) estimates.

Point Cloud Tree Height Validation
To assess the quality of the NAIP point cloud heights in areas with trees and forest we compared the maximum heights from the ALS and NAIP point clouds to the maximum tree heights taken from the FIA measurements.For this plot-level assessment, we first clip the 3D points falling within a 7.32 m (24 ft) radius circle centered on each FIA plot.Next, using the lidR package in R [65,66] we derive normalized heights by subtracting surface elevation using the corresponding 1 m and 10 m DEMs, then the resulting maximum height values at each resolution are recorded for each plot.This derivation of point cloud maximum height follows the area-based approach [67] and is equivalent to the highest point (in meters) sensed by the point cloud datasets within the 0.02 ha (0.04 acre) area surrounding the center of each FIA plot.Using the FIA tree measurements from this area, which is analogous to the area covered by FIA subplot 1, we present a series of scatter plots to help visualize the impact of increasing the resolution of the DEMs used to derive normalized surface heights, as well as to explore differences in the quality of the NAIP point cloud data by state and species group.The scatter plots include fitted regression lines (R 2 ) and correlation statistics, Pearson r, root mean square error (RMSE), and bias calculated as: (1) where y = point cloud height for the ith observation, x = FIA tree height for the ith observation, y = average point cloud height and x = average FIA tree height.Note for the RMSE and bias estimators x is considered the true value and x the true mean.Since some of the FIA plots were measured up to six years before the NAIP imagery was collected, we used 30 m forest disturbance data from the Landscape Change Monitoring System (LCMS, [68]) to exclude plots that experienced forest cover loss between the year of FIA field measurement and 2018.In addition, a very small number of plots with NAIP maximum heights > 50 m were also excluded, as these were determined to be outliers resulting from noise in the point cloud data.
As 1 m ALS heights are known to be of high quality [69,70], we use them here as part of the validation process to help shed light on the error structure of the NAIP data, as well as to determine if the FIA measurements have more error in certain places or are biased as trees get taller.Given the number of different field crews used to collect FIA data, coupled with challenges associated with measuring tree heights from ground-level, we assume the FIA measurements also contain errors.Therefore, to evaluate the strength of the height relationships we also compare R 2 's using reduced major axis (RMA) regression [71], which makes no assumptions about errors in either the X or Y variables.Lastly, to get a sense of how the height errors propagate through time, we present a series of violin plots showing the distribution, median, and direction of height differences between the ALS and NAIP point clouds and the FIA tree heights collected annually from 2012 to 2019.

Post-Stratified Estimation
To test the utility of the NAIP height maps for post-stratification, we produce a series of FIA estimates using the Forest Inventory ESTimation and Analysis (FIESTA) package available in R [72].Because the effectiveness of post-stratification is directly related to the correlation between the strata maps and the variable being estimated, we expect the benefit of using the NAIP height strata will fluctuate for different variables and geographic areas of interest.Therefore, estimates for three different variables, forest area (in ha), total forest volume (in m 3 ), and forest volume/area (in m 3 /ha), are derived independently for each of the survey units found in TN and VA (see Figure 2B for survey unit locations).Using these estimates, we report RE FIA , which compares the variances obtained from the NAIP strata maps to those from FIA's current estimation approach (using PS TCC for post-stratification) and RE SRS , which compares the variance from the NAIP strata maps to the variance obtained from simple random sampling.For each state and survey unit, RE FIA is reported for all three NAIP strata maps (PS DHM , PS CHM and PS CHM+FT ), whereas RE SRS is only reported for the strata map producing the highest RE FIA .For the best strata map we also report the mean value of the variable of interest, its percent standard error, as well as the number of plots gained or lost [derived as (RE FIA − 1.00) × n], and the apparent sample size after post-stratification (derived as n + the number of plots gained or lost).Note that for certain variables and survey units it is possible the NAIP strata maps will result in higher variance than FIA's current approach.In these cases, the estimates and standard errors are reported for the tree canopy cover strata map (previously denoted as PS TCC ), whereas the plot gains and apparent sample sizes are reported as losses relative to the initial sample.To help determine which of the FIA variables were most improved by the NAIP height strata, and to get a sense of which of the NAIP maps were most useful for post-stratification, we also provide a summary of these results, showing the overall number of improved estimates and the average number of plots gained/lost by state and survey unit.
A number of helpful functions in FIESTA facilitated importing and clipping the survey unit boundaries and strata maps, as well as assigning the strata map values to the forest inventory plots.After processing the boundaries and strata maps, the modGBpop function is used to define the area units, FIA plot, tree and condition tables, and the map to be used for post-stratification.The characteristics defined by modGBpop are then used in the modGBratio function to generate per-acre estimates by domain and estimation unit using calculations based on chapter 4 of [73].Here, we use FIA's base grid plots (i.e., PLOT.INTENSITY = 1) and the ratio of means estimator (described in Section 4.3.4 of [73]) to derive per-hectare estimates of live forest volume using the FIA tree variable VOLCFGRS.To get estimates on a per-hectare basis, the ratio estimator (Equation ( 4) divides the variable of interest (in this case total live volume) by the estimate of forest area in the population of interest (in this case the FIA survey units in TN and VA), resulting in three separate estimates from which to gauge the reduction in variance that occurs from using the various NAIP height maps for post-stratification.
Here Va d is the estimated forest volume per hectare in the dth survey unit, Vt d is the estimated total forest volume in the dth survey unit, Fa d is the estimated forest area in the dth survey unit, Vt hd is the estimated total forest volume in stratum h of the dth survey unit, Fa hd is the estimated forest area in stratum h of the dth survey unit, W h is the weight for stratum h and H is the total number of strata.The variance of the forest volume per hectare estimate v Va d is derived using Equation (5) defined as: Lastly, to get a sense of how well the NAIP strata perform for hardwood and softwood species, we also develop ratio estimates of live forest volume for the individual species groups found in each FIA survey unit.Note that at this level of estimation, several species had very low sample sizes so only the ones with a percent standard error ≤ 20% are reported.Similar to the broader survey-unit results, we summarize these estimates to determine which of the NAIP strata maps is most effective in reducing variance at the species group-level and to quantify the number of plots gained (or lost) on average for hardwood and softwood species.

Point Cloud Validation
The plot-level validation results comparing the ALS and NAIP point cloud heights and the FIA field-measured tree heights (shown in Table 3) indicate that, regardless of the source or DEM resolution, the point cloud maximum heights underestimated the maximum heights obtained from the field measured trees.Across both states the 1 m ALS point clouds had an average height difference of −0.36 m, which was significantly less than the −1.85 m and −1.27 m average difference observed for the 1 m and 10 m NAIP products, respectively.The results also showed that the average amount of underestimation for both ALS and NAIP is lower for hardwoods (ALS 0.24 m, NAIP −1.34 m) than for softwoods (ALS −0.96 m, NAIP −1.78 m).Despite the overall underestimation, the level of correlation between the point cloud maximum heights and FIA maximum tree heights was moderately to fairly strong, with average across state Pearson r values of 0.82 for the 1 m ALS and 0.76 and 0.72 for the 1 m and 10 m NAIP data, respectively (Table 3).We also found correlations between the point cloud heights and FIA maximum tree heights were much lower for hardwoods (ALS 0.77, NAIP 1 m 0.67, and NAIP 10 m 0.64) than softwoods (ALS 0.88, NAIP 1 m 0.86, and NAIP 10 m 0.79).The patterns of underestimation of maximum height by the point clouds and the generally better performance for hardwood species revealed in the maximum height comparison are also reflected in the scatter plots and fitted RMA regression lines shown in Figure 4.For example, the scatter plots reveal a much tighter fit and less noise for the point cloud heights derived with the 1 m vs. 10 m DEM data.This is true for softwoods more so than hardwoods.Overall, the 1 m ALS maximum heights have the best model fits, followed closely by the 1 m NAIP (Table 3 and Figure 4).The scatter plots also show a much higher level of noise in the NAIP point cloud data collected in VA.This is apparent by the larger number of significant errors along the bottom part of the x-axis, where NAIP maximum height is at or near zero but FIA maximum heights range from 5-36 m.While the count shading shown in Figure 4 reveals the majority of NAIP heights do fall along the 1:1 line, the fitted relationships show that regardless of DEM resolution, the NAIP point cloud heights are underestimated at the high end (i.e., maximum heights ≥ 30 m) and overestimated at the low end (i.e., maximum heights 0-10 m) compared to the FIA tree measurements.Because ALS LiDAR itself has proven to be an accurate method for measuring the height and structure of trees and vegetation, we also directly compared the 1 m ALS and NAIP maximum heights to help inform the interpretation of their relationships with FIA data.In general, the ALS and NAIP maximum heights are moderately (R 2 = 0.53) to strongly (R 2 = 0.87) related, as can be seen by the fitted RMA regression lines shown in Figure 5.These scatter plots show the relationship between NAIP and ALS is stronger in TN than in VA (e.g., R 2 s in TN are 0.22 to 0.32 higher) and that the relationship improves when the time between acquisition dates is kept within 2 years (Figure 5, right panel).This is especially true in TN where R 2 improves from 0.75 to 0.87 within this narrower acquisition window.Aside from a few isolated errors at the high end, most of the noise in the NAIP data seems to occur in maximum heights ≤ 25 m.In this range, the NAIP data often significantly underestimates the ALS heights, especially in VA (see lower row, Figure 5).Model fits notwithstanding, the count shading and alignment of the individual errors in relation to the 1:1 line show that above 30 m, the NAIP and ALS maximum heights are generally in good agreement, with little to no systematic bias present (especially in TN, upper right, Figure 5).
Aside from a few isolated errors at the high end, most of the noise in the NAIP data seems to occur in maximum heights ≤ 25 m.In this range, the NAIP data often significantly underestimates the ALS heights, especially in VA (see lower row, Figure 5).Model fits notwithstanding, the count shading and alignment of the individual errors in relation to the 1:1 line show that above 30 m, the NAIP and ALS maximum heights are generally in good agreement, with little to no systematic bias present (especially in TN, upper right, Figure 5).Finally, we plotted the differences between the FIA and point cloud maximum heights using a series violin/box plots broken out by FIA measurement year (Figure 6).These plots indicate that, regardless of the DEM resolution or year of FIA measurement, the medians of the point cloud maximum heights are consistently 1 to 5 m less than the FIA measured tree heights.Overall, the box plots in Figure 6 show that across time, the median height differences for the 1 m ALS are lower, and the residual errors are better balanced above and below zero, as opposed to the 1 m NAIP, which has most of its errors consistently below zero.Although the 10 m NAIP contains slightly less canopy detail (Figure 7) its height errors display slightly less systematic underestimation than the 1 m NAIP.Starting in 2016 through the end of 2019 a slight trend of increasing underestimation is observed for both the NAIP and ALS point clouds (Figure 6).
Prior to running wall-to-wall, statewide DHMs and other heights metrics, we used FUSION to process the NAIP data for a number of smaller test areas to help determine the most appropriate grid cell resolution to best meet the needs of our study.Outputs from these initial tests showed that although the 2 m DHMs captured more canopy detail than the 10 m DHMs (e.g., see Figure 7A vs. Figure 7B), the overall differences were not deemed significant enough to justify the additional processing time and storage space required to develop the higher resolution product.Therefore, 10 m statewide DHMs were used as the basis for developing the NAIP strata maps (PS DHM , PS CHM and PS CHM+FT ) for TN and VA (Figure 8).In addition, Figure 8 also highlights the difference between the PS DSM , which includes heights for buildings and other non-forest structures, and PS CHM , which removes heights in non-forest areas using a forest mask.Finally, we plotted the differences between the FIA and point cloud maximum heights using a series violin/box plots broken out by FIA measurement year (Figure 6).These plots indicate that, regardless of the DEM resolution or year of FIA measurement, the medians of the point cloud maximum heights are consistently 1 to 5 m less than the FIA measured tree heights.Overall, the box plots in Figure 6 show that across time, the median height differences for the 1 m ALS are lower, and the residual errors are better balanced above and below zero, as opposed to the 1 m NAIP, which has most of its errors consistently below zero.Although the 10 m NAIP contains slightly less canopy detail (Figure 7) its height errors display slightly less systematic underestimation than the 1 m NAIP.
Starting in 2016 through the end of 2019 a slight trend of increasing underestimation is observed for both the NAIP and ALS point clouds (Figure 6).
Prior to running wall-to-wall, statewide DHMs and other heights metrics, we used FUSION to process the NAIP data for a number of smaller test areas to help determine the most appropriate grid cell resolution to best meet the needs of our study.Outputs from these initial tests showed that although the 2 m DHMs captured more canopy detail than the 10 m DHMs (e.g., see Figure 7A vs. Figure 7B), the overall differences were not deemed significant enough to justify the additional processing time and storage space required to develop the higher resolution product.Therefore, 10 m statewide DHMs were used as the basis for developing the NAIP strata maps (PSDHM, PSCHM and PSCHM+FT) for TN and VA (Figure 8).In addition, Figure 8 also highlights the difference between the PSDSM, which includes heights for buildings and other non-forest structures, and PSCHM, which removes heights in non-forest areas using a forest mask.

Forest Volume
The NAIP strata maps were used to produce post-stratified volume per area (m 3 /ha), total volume (m 3 ), and forest area (ha) estimates for the 10 survey units in TN and VA (Table 4).These results show that, for both total forest volume and volume/area estimates, the NAIP strata maps yielded lower variances than the NLCD tree canopy cover map (PSTCC) currently used by FIA.Although all 20 individual survey unit volume estimates were improved by the use of the NAIP strata maps, the level of improvement varied by survey unit.For example, based on the bold shading in Table 4 (which indicates the best overall performing NAIP strata map per variable in terms of REFIA), it is apparent that for

Forest Volume
The NAIP strata maps were used to produce post-stratified volume per area (m 3 /ha), total volume (m 3 ), and forest area (ha) estimates for the 10 survey units in TN and VA (Table 4).These results show that, for both total forest volume and volume/area estimates, the NAIP strata maps yielded lower variances than the NLCD tree canopy cover map (PS TCC ) currently used by FIA.Although all 20 individual survey unit volume estimates were improved by the use of the NAIP strata maps, the level of improvement varied by survey unit.For example, based on the bold shading in Table 4 (which indicates the best overall performing NAIP strata map per variable in terms of RE FIA ), it is apparent that for forest volume/area and total volume the NAIP strata maps performed best in survey units T4 and T2 in TN and V4 and V2 in VA.The overall highest relative efficiency gains relative to the FIA approach (RE FIA ) were in TN, where both forest volume/area and total volume estimates resulted in maximum RE FIA 's = 1.46.This level of variance reduction is analogous to a 46% increase in the sample size used to make the estimates, which is analogous to adding 251 plots to the volume/area estimate in T4, and 190 plots to the total volume estimate in T2 (Table 4).In VA, the gains for forest volume were slightly lower than in TN, with maximum RE FIA's of 1.38 for total forest volume and 1.33 for forest volume/area.The 38% increase in sample size in V2 is like adding 265 plots to the total forest volume estimate, while the 33% increase in V4 theoretically adds 167 plots to the V4 forest volume/area estimate (Table 4).While the NAIP strata maps resulted in similar levels of RE FIA for both types of volume estimates, the gains relative to simple random sampling (RE SRS ) were two times higher for total forest volume (which ranged from 2.28 to 2.87) than for forest volume/area (which ranged from 1.24 to 1.58) (Table 4).
To get a sense of the overall performance of the NAIP strata maps across units and to help identify which of the strata maps performed best, we summarized the individual survey unit results in Table 4 to get the average percentage of plots gained/lost resulting from the use of the NAIP strata maps, as well as to sum the number of times each individual map resulted in the highest RE FIA (denoted by the gray shading in Table 4).This summary (Table 5) indicates that the average percent of plots gained across all the survey units for forest volume ranged between 24% and 36%.The lowest percentage gain for volume, 15% for V3 volume/area, is well above zero indicating there was no penalty for using the NAIP strata maps to post-stratify the forest volume estimates.Table 5 also shows that for forest volume, the PS CHM strata map produced the lowest variances in 7 of the 10 survey units.This was especially true for total volume, where it had the best results in 4 out of the 5 survey units in both TN and VA (Table 5).Note that because the results in Table 4 allow for ties (with gray shading used for multiple maps) there were two survey units (T4 and V3) where the PS DHM strata map equaled the performance of the PS CHM map.Improved survey unit represents the number of survey units where RE FIA is greater than 1.Best PS map indicates the number of times each NAIP strata map was found to be the best method (note in the event of a tie units can have more than 1 best map per survey unit).Plot gain/loss column shows the average percentage of plots gained or lost relative to the initial sample size, as well as the minimum and maximum number of plots gained per survey unit in each state (note losses relative to the NAIP maps are shown as negatives).

Forest Area
In contrast to the consistently positive efficiency gains observed for the forest volume estimates, the NAIP strata maps were much less effective at improving the estimates of forest area (Table 4).For example, in TN the NAIP DHM strata maps did produce lower variances than the FIA approach in 4 of the 5 survey units; however, in 2 of the improved survey units the RE FIA 's were barely above 1.00 (see T3 and T5, Table 4).The other two survey units did see slightly higher gains with RE FIA 's of 1.24 (T2) and 1.11 (T4), which resulted in 100 and 61 additional plots, respectively.In the fifth survey unit (T1), the NAIP strata map generated estimates with higher variance than the current FIA approach, resulting in a theoretical reduction of 21% of its sample plots (Table 4).Overall, on average there was a 4% gain in the number of plots added to the forest area estimates in TN.
On the other hand, in VA, only one of the 5-survey unit forest area estimates was improved by the NAIP strata maps (Table 4).The 8% maximum gain observed in V2 (Table 4) was the maximum achieved in VA, resulting in the theoretical addition of 55 plots to the sample (Table 4).Use of the NAIP height strata in the other 4 survey units yielded forest area estimates with higher variances than the current FIA approach, resulting in a theoretical loss of 22 (V4) to 119 (V1) sample plots (Table 4).Across all the survey units in VA, there was a −8% average reduction in sample size compared to estimates from FIA's current post-stratification approach (Table 5).
Although the NAIP strata maps were not particularly effective in reducing the variance of the forest area estimates relative to FIA's current approach, the gains relative to simple random sampling were the highest of all the variables tested (ranging from 2.48 (V2) to 4.05 (T2), Table 4).This level of improvement, which yields a 2.5 to 4× increase in sample size, far exceeds the gains achieved for forest volume (see RE SRS 's, Table 4).

Species Group Volume
Post-stratified forest volume estimates were also run for the individual species groups found in the 10 survey units in TN and VA.Each of the survey units contained a unique combination of the 9 softwood and 28 hardwood species groups used by FIA in the eastern U.S (Table S1).In total there were 84 (TN) and 80 (VA) species group estimates which produced significant results (i.e., standard errors ≤ 20%).Similar to the general forest volume estimates reported in Table 4, the results for the individual species groups (Table S2) show that in 4 out of 5 survey units in TN and 3 out of 5 survey units in VA, a softwood species group achieved the maximum gain from using the NAIP strata maps.In TN, the volume/area estimate for the other eastern softwood group in unit T3 had the maximum RE FIA of 1.72, while in VA the maximum RE FIA of 1.60 was observed for the eastern white pine group in unit V2 (Table S2).To get a sense of performance of the NAIP strata maps across units, and to determine which map was most effective at lowering the variance of the species group volume estimates, we summarized the individual species results by grouping them into hardwood and softwood classes.These results (Table 6) show that across both states, forest volume/area was improved the most by the NAIP height strata, with 85% (TN) and 86% (VA) of its species group estimates improved versus 75% (TN) and 66% (VA) for total forest volume (see the number of improved species, Table 6).Overall, 78% of the significant species group volume estimates were improved by the NAIP strata maps.Table 6 also shows that on average softwood species groups tended to gain the most, with average gains of 10-20% versus 4-7% for hardwoods.Note there were some penalties for using the NAIP height strata, as some of the hardwood and softwood species groups had minimum gains well below zero, indicating a theoretical loss of plots compared to the use of the tree canopy cover strata (PS TCC ).Total species: the number of species groups observed in the FIA sampling plots.Improved species: the number of species groups which have RE FIA greater than 1.00.Softwood plot gain/loss: the number of plots gained (or lost) for softwood species groups by using the NAIP strata maps instead of tree canopy cover data.Hardwood plot gain/loss: the number of plots gained (or lost) for hardwood species groups by using the NAIP strata maps instead of tree canopy cover data.
To determine which NAIP strata map performed best for the species group volume estimates, we tabulated the number of times each map recorded the highest gain (Table 6).Overall, across both states the NAIP PS CHM+FT map had the highest RE FIA in 59% of the improved species group volume estimates.The results also show that the PS CHM+FT map had the highest RE FIA in 65% of the improved species group models in TN and 52% in VA.The other two NAIP strata maps, improved ~20% (PS CHM ) and ~30% (PS DHM ) of the species group volume estimates in TN and VA, respectively.

Validation of NAIP DAP and ALS Point Clouds
To evaluate the quality of the remote sensing point clouds we used a 3-way validation approach to compare ALS and NAIP DAP maximum surface heights with FIA field measured tree heights.This approach not only allowed a robust evaluation of the point clouds ability to accurately resolve the maximum height of trees and forest canopies, but given the high accuracy of ALS LiDAR (e.g., [33,74,75]), it also allowed us to evaluate potential bias in the FIA field measurements.As expected, these comparisons (presented in Figures 4 and 5 and Table 3) show the ALS point cloud maximum heights are in better agreement with the FIA tree heights than the 1 m NAIP data, as evidenced by their higher correlations (r ranging from 0.74-0.88for ALS vs. 0.63-0.89for NAIP) and levels of explained variance (R 2 s ranging from 0.55-0.78for ALS vs. 0.39-0.79for NAIP).
Perhaps a more important criterion of data quality is the level of systematic bias in the point cloud datasets.Overall, regardless of the grid cell size of the underlying DEM, the species group, or the state of data collection, the NAIP point clouds underestimated the FIA measured tree heights by an average of −1.56 m (Table 3).Similar levels of underestimation, particularly at the high end of the height range (>15 m), have been reported in other studies (e.g., [76,77]).We also found that the NAIP point clouds had higher levels of bias for softwoods (average bias = −1.78m) than hardwoods (average bias = −1.34m), even though in most cases the softwood heights were more highly correlated with the FIA measurements.This occurs because the flaws and artifacts introduced by disturbance and phenology have a lower impact on the NAIP softwood heights (particularly in TN, see Figure 4), so the bias calculations are more impacted by the underestimation that occurs at the high end of the height spectrum (i.e., maximum heights > 30 m).On the other hand, the NAIP hardwood heights displayed a lot more noise that resulted in distinct overestimation at the low end and underestimation at the high end of the height range (see Figure 4), which is much more problematic to account for when making wall-to-wall maps.This phenomenon also explains why the 10 m NAIP heights have less bias than the 1 m NAIP heights (−1.27 m vs. −1.85m, Table 3), although the 10 m data clearly have more error and residual scatter in their fitted relationships with the FIA measurements (Figure 4).These examples illuminate the challenge of interpreting bias simply by the calculated values for the different point cloud datasets.

NAIP Point Cloud Anomalies
While ALS performed better overall, a closer look at the scatter plots in Figure 4 reveals that most of the improvement occurred in hardwoods, where ALS had far fewer errors along the lower part of the x-axis (e.g., where FIA maximum height is high and point cloud maximum height is near zero).These errors, which are evident in all the NAIP plots except TN softwoods (top two rows of Figure 4), have two causes.First, some of these errors are caused by gaps of missing data where the NAIP imagery was acquired in leaf-off conditions.In areas with missing leaves the photogrammetric processing algorithms were unable to successfully triangulate enough 3D points, resulting in zero heights above ground.Although most of the NAIP images covering the study area were acquired between October and December (i.e., late fall/early winter, see Figure 3D), the impact of missing leaves was much more pronounced in VA, where a notable increase in the frequency and magnitude of errors along the x-axis can be seen especially for hardwoods (Figure 4).Overall, the effect of missing leaves was mostly confined to eastern TN and western VA (i.e., survey units T5, V1 and V2, Figure 2) along the higher elevations of the Blue Ridge mountains, where cooler temperatures resulted in earlier periods of senescence and leaf-drop compared to lower elevations.The high spatial variability in leaf-off conditions across the entire study area precluded us from using acquisition date alone to identify areas where the quality of the NAIP point cloud was poor.In fact, in many cases (particularly in TN) the late season imagery captured active periods of senescence where leaves were in the process of changing colors but had not yet fallen off the trees.Our analysis found that in these situations the NAIP height estimates still tended to be of high quality.
The late season acquisition window and steep topography also caused isolated patches of deep shadow, which resulted in extremely low point returns.Like the areas of missing leaves, shadowed areas (highlighted in red, Figure 3A-C) tended to occur in isolated patches containing 1-30 points per 10 m grid cell.Unlike areas of missing data, which had 0 points and no height estimate, areas of deep shadow often resulted in erroneously high height values, a finding that has been previously reported in other studies [76].For the point cloud height validation, we avoided these areas by eliminating plots with estimated height values > 50 m, which is the height of the tallest tree measured by FIA in TN and VA.Overall, fewer than 10 plots were omitted for this reason.Although areas of missing data (due to missing leaves, deep shadows and other anomalies) represented a small percentage of the full NAIP data set, we opted to leave the impacted plots in the validation since these areas are inherent flaws in the data that ultimately limit its effectiveness when used for post-stratification of the FIA estimates (discussed below).

Forest Disturbance
Although the inherent flaws in the NAIP imagery did limit its overall quality, the biggest factor contributing to the large errors in the NAIP point cloud heights is forest management and disturbance.Although we attempted to remove impacts of forest harvesting using LCMS data, there were still many instances where this screening process did not remove all artifacts related to disturbance, resulting in substantial deviations in observed and predicted heights along both the x-and y-axes.This is particularly evident in VA where partial harvesting and thinning operations, which have been shown in other studies to cause errors in NAIP CHMs [77], resulted in additional noise.
Overall, there were far more errors in the VA NAIP data (Figure 4, right two columns) than in TN (Figure 4, left two columns), likely due to better image acquisition dates in TN resulting in higher point densities compared to VA (Figure 3A-D), as well as lower rates of disturbance (10.8% in TN vs. 14.7% in VA based on LCMS data).We also found there was very little impact of using the 1 m vs. the 10 m resolution DEM for height normalization of the NAIP point cloud data, as evidenced by the fact that the amount of residual scatter around the fitted RMA regression models in TN (Figure 4) did not significantly worsen with coarser DEM resolution.This bodes well for the use of 10 m DEMs when running wall-to-wall height maps for use in post-stratification.

Temporal Offset
One aspect that contributes to the noise added by disturbance is the amount of time between the FIA measurements and the date the NAIP images were collected.Figure 6 shows that neither the DEM resolution (1 m vs. 10 m) nor the number of years between field measurement and image acquisition (ranging from 6 years before to 1 year after) significantly impacts the error structure or bias of the NAIP maximum height estimates.For example, the median errors observed across time range between slightly less than 0 and −3.5 m, with the majority (25-75%) of data each year falling within +/− 3.0 m (Figure 6).While Figure 6 does indicate a slight increase in the amount of underestimation as the FIA year approaches the NAIP acquisition year (i.e., 2018), the effect is minimal, with very little of the data exceeding +/− 5.0 m.Overall, the lack of temporal degradation in the NAIP heights is encouraging, as FIA estimates are always made with plots that range backwards in time across a number of years, thus our results suggest NAIP point clouds should continue to provide gains in precision for a number of years after its acquisition date.

GPS Imprecision
Another factor that is difficult to quantify but certainly impacts the level of agreement between the point clouds and field plots is the inherent imprecision of the FIA plot locations.In the southeastern U.S., FIA currently uses recreational grade GPS receivers to collect plot coordinates in the field.Although a broad assessment of FIA plot coordinate precision is beyond the scope of this paper, initial results from a pilot study aimed at collecting High-Precision Global Navigation Satellite System (HPGNSS) coordinates for forested FIA plots in VA reveals that typical FIA plot coordinates collected with recreational grade GPS receivers have an RMSE of ~8.0 m (n = 31) when compared to HPGNSS locations.Although [41] concluded that global positioning error was not a significant factor in determining the level of LiDAR-assisted post-stratification improvement of aboveground biomass estimates in Minnesota, their reasoning was based on the fact that plots measured with HPGNSS and recreational receivers resulted in similar mean estimates of biomass, despite finding that the ratio of standard errors indicated a relative gain in efficiency of around 50% when using the more precise coordinates.This level of gain is equivalent to adding 50% more plots to the FIA sample, which is a significant improvement in both precision and potential cost savings that should not be overlooked.While we did not explicitly examine the effects of plot location accuracy, other studies have shown offsets from imprecise GPS locations can significantly degrade the relationship between point cloud heights and FIA field measurements, leading to a reduced gain in estimated precision compared to plots with HPGNSS locations (e.g., see Table 1 in [27]).
The amount of noise added due to GPS inaccuracy will undoubtedly vary, depending on the resolution of the underlying DEM and the heterogeneous nature of the landscape in which the plot is located.For example, if a 10 m DEM (which covers 100 m 2 ) is used to normalize the point cloud heights, the currently observed offset in VA (+/− 8 m) could cause a plot to be linked to a pixel that mostly lands outside the field measured subplot, which covers 168.16 m 2 (or roughly a 13 m pixel).On the other hand, if a 1 m DEM is used, and an area-based calculation is performed summarizing all the pixels that are purported to cover the FIA subplot area, the GPS plot locations will result in pixels outside the plot area contributing to the calculation; however, even if a portion of the field area is covered, the inclusion of some of the plot pixels helps smooth the noise and ultimately limits the effect of the spatial location error.This is especially true in dense forests where shifting the spatial location of plot center has little effect due to homogenous canopy conditions.However, in highly managed forests where frequent harvesting and thinning result in a patchier, more structurally diverse landscape, the effect of GPS imprecision may be more significant.For instance, in TN where rates of forest management are lower and the quality of the NAIP imagery is higher, upscaling the NAIP point cloud maximum heights from 1 m to 10 m did not result in an expansion of errors or addition of significant outliers compared to the FIA measurements.On the other hand, in VA, where disturbance rates from harvesting are much higher and the quality of the NAIP imagery is lower, the upscaling from 1 m to 10 m resulted in a proliferation of new errors along the x-axis.While this evidence is far from conclusive, it does warrant further investigation, especially as FIA seeks ways to maximize the efficiency gained from using point clouds in its operational estimation system.

FIA Tree Measurement Uncertainty
Because nearly all the point clouds displayed some level of underestimation for maximum heights > 30 m and given other studies have shown errors in field measured tree heights as trees get taller [78,79], we compared the 1 m NAIP and 1 m ALS maximum heights directly, with the idea of using the proven accuracy of ALS to check for potential bias in the FIA measurements.As with the FIA height validations (presented in Figure 4), we also attempted to remove the effects of disturbance using the LCMS data.Despite our attempts, Figure 5 shows a good amount of residual scatter remained after filtering for disturbance, particularly when the time between measurements exceeded two years (left column, Figure 5).The upper right panel of Figure 5 shows the best overall relationship between ALS and NAIP maximum height occurred in TN, where the lower rate of disturbance and reduced measurement gap of two years yielded a strong linear relationship (R 2 = 0.87) between the two remote sensing-based height estimates.The fact that this relationship also displays very little residual scatter and more importantly, next to no bias at either end of the height spectrum, suggests the remote sensing point clouds may not underestimate the FIA height measurements as much as Figure 4 might suggest.This is especially true as other studies have also shown better, less biased relationships between LiDAR and NAIP point cloud heights than between NAIP and field measured tree heights (e.g., see [76]).Since FIA tree heights are only recorded to the nearest foot from ground level and considering the challenges associated with locating the highest part of upper canopy trees in dense deciduous forests, there are no doubt errors in the FIA tree heights.Although the ALS and NAIP maximum heights in VA are much noisier (and thus less valuable for evaluating agreement) there is still a noticeable lack of systematic bias evident in these graphs (especially in the lower right panel, Figure 5), which further supports the fact that, despite the effects of disturbance and other noise, the NAIP point cloud heights are in excellent agreement with the best available remotely sensed measures of height from ALS LiDAR.

Post-Stratifying FIA Estimates with Digital Height Maps vs. Tree Canopy Cover Data
Our other main objective, aside from evaluating the quality of the NAIP point cloud heights, was to quantify the level of precision gained from using NAIP statewide height maps to post-stratify FIA forest area and volume estimates.The results in Table 4 clearly show that for forest volume, the height based NAIP strata maps did a better job of reducing variance than the strata maps based on tree canopy cover data.In fact, although the level of relative efficiency gained varied by state (avg RE FIA in TN = 1.35 vs. VA = 1.28) and survey unit (min/max RE FIA across TN = 1.28-1.46 vs. 1.15-1.38 in VA), the NAIP height maps outperformed canopy cover for all 20 volume estimates in TN and VA (Table 4).Overall, for forest volume we found that on average, the best NAIP strata map (shaded in gray, Table 4) resulted in a 31% gain in sample size compared to the tree canopy cover maps.This level of improvement resulted in the theoretical addition of 175 sample plots per survey unit on average.Although the least number of plots gained for any one survey unit was 71 (V3 in VA Northern Piedmont, Figure 2), the other survey units gained between 126-265 plots depending on the initial sample size and strength of the NAIP heights relationship with forest volume.The quality of this relationship hinges on the level of noise introduced by GPS imprecision (resulting in spatial mismatch between plots and pixels) and our assignment of plots to strata using a single 10 m pixel (which introduces error because only a portion of the area sampled on the ground is being captured by the remote sensing data).
Because the amount of efficiency gained in the sample is based on the strength of the relationship between the map and the variable of interest (i.e., R 2 , see [27]) substantial errors in the Y direction will serve to flatten the fitted linear relationship, causing the predictions to regress to the sample mean, which ultimately negates the benefit of having wall-to-wall height maps.These errors, which are readily apparent in scatter plots between NAIP maximum height and FIA forest volume (m 3 /ha) (Figure 9), ultimately reduce the strength of the fitted regression relationships, especially for softwoods, where the outliers along the y-axis capped R 2 below 0.5.Minus these outliers, the R 2 's for softwoods would likely have approached 0.7 or higher given their similarity to the relationships presented in Figure 7 of [27].Due to the exponential relationship between R 2 and RE shown in Figure 5 of [27], there is likely a twofold level of improvement that could still be obtained for softwood volume just by eliminating the errors in the NAIP height maps.In light of these potential improvements, Figure 9 also highlights the enormous potential of NAIP height maps to achieve even further gains in efficiency vs. tree canopy cover maps, which in reality have a poor relationship with volume and only serve to reduce estimate variance because of a slightly positive, but mostly spurious relationship with volume (Figure 9, bottom row).However, when estimates were made for specific species groups (Table S2), the added forest type information in the PSCHM+FT map began to take on added importance.For example, of the 78% of the species group estimates that were improved by the NAIP height maps, just over half (53%) were attributed to the NAIP CHM with added forest type information (Tables 6 and S2).While the percentage improvement tended to be lower at this scale of estimation, we note that in many cases the difference varied between realizing slight positive gains in efficiency vs. small losses in sample size compared to the tree canopy cover strata maps (Table S2).Overall, our results show that strata maps developed from NAIP height information are clearly better at reducing the variance of forest volume estimates than strata derived from tree canopy cover data, and when the NAIP strata are better, information from NAFD (e.g., PSCHM) and NLCD (PSCHM+FT) can lead to even bigger gains in efficiency than just using the wall-to-wall NAIP DHM (PSDHM), which only achieved the highest gains for 20% and 23% of the survey unit and species group estimates, respectively.4.3.6.Potential for Operational Use of NAIP Point Clouds within FIA Because the net effect of post-stratification can be quantified in terms of the theoretical number of plots gained (or lost) relative to the initial sample size, perhaps the best way to determine the level of added benefit provided by the statewide NAIP point clouds is to compare the amount it would take to field measure the number of additional plots gained from post-stratification vs. the total cost of the NAIP point cloud data.For example, at a cost of USD 1000 per field visit the 71 plots added to the forest volume/area estimate in unit V3 (Table 4) equates to USD 71,000, which is slightly more than it cost to purchase both the TN and VA NAIP point clouds.In other words, for reasonably sized states, the potential savings achieved in just one survey unit will more than cover the cost of purchasing the NAIP point cloud data.At the high end, stratifying the survey unit total forest What is also apparent from Table 4 is that, unlike the improvements realized for volume, the NAIP height maps were not a clear improvement over the canopy cover maps when it came to stratifying the FIA forest area estimates.Although 4 of the 5 survey units in TN were improved by the NAIP height maps, the level of relative efficiency gained (avg RE FIA = 1.10) was 20% lower than was observed for volume.This average level of improvement is somewhat misleading because survey units T3 and T5 only saw 1-2% improvements in sample size, leading to only 6 and 15 plots gained, respectively (Table 4).The fact that only 1 of 5 survey units in VA was improved further draws into question the usefulness of the NAIP height maps for stratifying forest area estimates.Overall, the NAIP strata maps mostly resulted in less precise forest area estimates which at times caused a theoretical loss in sample size in VA compared to the estimates stratified with the tree canopy cover map.This is not terribly problematic, nor is it terribly surprising.Canopy cover is a better indicator of where trees and forests are located spatially (due to the satellite data's ability to capture information about leaf area), and because the tree canopy cover map is in the public domain and is a product co-produced by the USFS, it can continue to be used by FIA for stratification at no additional cost to the program.
Although not as widely successful as the survey unit estimates, the NAIP strata maps also did a nice job of improving the species group volume estimates.With 66-86% of the significant species group estimates gaining some benefit from the NAIP strata maps it is clear they still add value, even though on average only 5 and 15 plots were gained for hardwoods and softwood models, respectively (Table 6).Although softwood species groups typically had fewer plots, they clearly benefited more from the NAIP height maps than hardwoods, with maximum gains adding between 52-72 plots to their overall sample size (Table 6).Despite this success there were still times when the NAIP strata maps were not an improvement over tree canopy cover, resulting in a small loss of precision (and sample size) for both hardwood and softwood species groups.The same sources of error discussed above (and highlighted in Figure 9) impact the strength of the species group poststratifications and, because the sample size is much smaller, the presence of any outliers in the NAIP height maps is particularly problematic.For example, our results showed only a few extreme outliers (primarily in the Y direction) were responsible for limiting efficiency gains for the Loblolly/shortleaf pine and Other yellow pine volume estimates (FIA species codes 2 and 3-Table S1).
The last question pertaining to post-stratification is which of the NAIP strata maps is best?While our results are far from conclusive, they do show that, for traditional FIA estimates made at the survey unit level, combining the NAIP DHM with a secondary forest mask (in this case from NAFD) to produce a CHM (or PS CHM ) led to higher efficiency gains for many of the forest area and volume estimates (Tables 4 and S2).In fact, of the 83% of survey unit estimates that were improved by the NAIP height data across TN and VA, 67% of them involved use of the NAIP CHM (PS CHM ).Although there were times when the other NAIP strata maps (PS DHM and PS CHM+FT ) yielded better results, we note that in most of these instances the efficiency reported for PS CHM was only 1-2% lower than these other maps.It was also evident that, at this scale of estimation, adding the extra information about forest type to the NAIP height data (e.g., using NLCD) was not particularly useful.
However, when estimates were made for specific species groups (Table S2), the added forest type information in the PS CHM+FT map began to take on added importance.For example, of the 78% of the species group estimates that were improved by the NAIP height maps, just over half (53%) were attributed to the NAIP CHM with added forest type information (Table 6 and Table S2).While the percentage improvement tended to be lower at this scale of estimation, we note that in many cases the difference varied between realizing slight positive gains in efficiency vs. small losses in sample size compared to the tree canopy cover strata maps (Table S2).Overall, our results show that strata maps developed from NAIP height information are clearly better at reducing the variance of forest volume estimates than strata derived from tree canopy cover data, and when the NAIP strata are better, information from NAFD (e.g., PS CHM ) and NLCD (PS CHM+FT ) can lead to even bigger gains in efficiency than just using the wall-to-wall NAIP DHM (PS DHM ), which only achieved the highest gains for 20% and 23% of the survey unit and species group estimates, respectively.4.3.6.Potential for Operational Use of NAIP Point Clouds within FIA Because the net effect of post-stratification can be quantified in terms of the theoretical number of plots gained (or lost) relative to the initial sample size, perhaps the best way to determine the level of added benefit provided by the statewide NAIP point clouds is to compare the amount it would take to field measure the number of additional plots gained from post-stratification vs. the total cost of the NAIP point cloud data.For example, at a cost of USD 1000 per field visit the 71 plots added to the forest volume/area estimate in unit V3 (Table 4) equates to USD 71,000, which is slightly more than it cost to purchase both the TN and VA NAIP point clouds.In other words, for reasonably sized states, the potential savings achieved in just one survey unit will more than cover the cost of purchasing the NAIP point cloud data.At the high end, stratifying the survey unit total forest volume estimates with the NAIP height maps added 938 and 951 plots to the FIA sample in TN and VA, respectively, which in total nets nearly USD 2 million in potential cost savings from not having to measure 1889 additional field plots.As other studies have reported similar precision gains when using NAIP point clouds with area-level small area estimation approaches [80], there is growing evidence to suggest NAIP DAP is a highly cost-effective option for increasing the precision of NFI forest volume estimates.While these numbers are impressive, they omit the costs associated with processing and storing the NAIP point clouds, which are not insignificant.Based on our experience, the level of effort to produce 10 m statewide DHMs in FUSION can be a bit time consuming but is otherwise a straightforward process that should not impose significant operational or financial burdens; therefore, FIA should consider wider adoption and use of NAIP point clouds in its operational estimation system.We do note that computational limitations and storage requirements still make it difficult to operationally produce wall-to-wall 2 m DHMs, although we envision with continued improvements in computing technology that this barrier will be lessened in the future.Once wall-to-wall DHMs are produced it is not difficult to incorporate the NAIP height data into different stratification routines as we did here with FIESTA.Our results show use of different stratification maps will result in different levels of precision gain, so further testing is needed to fully understand the best way to maximize the use of NAIP point clouds in FIA's estimation system.As illuminated here, future FIA estimates may well benefit from use of both height and tree canopy cover maps, especially if more advanced model-assisted estimation approaches are used.

Conclusions
Our results revealed several interesting findings about the quality of NAIP DAP point clouds and their potential utility for improving forest inventory estimates in southeastern mixed hardwood forests.Overall, regardless of the underlying DEM resolution (1 m or 10 m), the NAIP point clouds underestimated FIA maximum tree heights by an average of −1.56 m but otherwise were in strong agreement with the field measured trees (R 2 's = 0.58-0.89).The biggest difference between the NAIP point cloud and field measured maximum tree heights occurred at the state level, where the underestimation in VA was 1 m more than in TN (bias in TN = −0.99m vs. −2.13m in VA), likely due to the use of different vendors and processing routines.In terms of post-stratification, using the NAIP DHMs instead of tree canopy cover maps resulted in forest volume estimates which were 15 to 46% more precise than FIA's current approach.The additional precision gained from using the NAIP height data added 175 plots per survey unit on average to the FIA sample, offering huge financial savings compared to the costs of installing more field plots.Although more testing is needed to determine the full benefit of using NAIP DAP point clouds in other forest types across the U.S., our results (and those of others, e.g., [27,76,77,80]) highlight the benefits, potential, and growing feasibility of using NAIP DAP point clouds to improve the precision of FIA's operational forest inventory estimates.Moving forward, future work will focus on using NAIP DHMs and tree canopy cover data with more sophisticated model-assisted approaches, such as the generalized regression estimator (GREG), which uses linear regression and calibration techniques such as raking and lasso to improve estimate precision.More work is also needed to improve the quality and consistency of the photogrammetric models used to derive 3D points from the NAIP stereo images so that point clouds from different states and acquisition dates can be used more consistently over time to estimate changes in forest growth.Collecting high precision GPS data at FIA subplot locations would help reduce misalignment errors, leading to cleaner, less noisy models and improved stratification performance.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10.3390/rs14174386/s1, Table S1: The full list of species groups in eastern U.S., Table S2: Total volume and Volume/area estimates by each species group and survey unit.

Figure 1 .
Figure 1.The workflow used to validate the NAIP point clouds and assess the NAIP digital height model (DHM)/canopy height models (CHMs) ability to post-stratify the FIA estimates.

Figure 1 .
Figure 1.The workflow used to validate the NAIP point clouds and assess the NAIP digital height model (DHM)/canopy height models (CHMs) ability to post-stratify the FIA estimates.

Figure 2 .
Figure 2. Study area showing (A) the states of Tennessee (TN) and Virginia (VA) with red outline indicating the 1 m DEM coverage and (B) the elevation and survey units in each state.Combinations of letters and numbers represent survey units.T1: West, T2: West Central, T3: Central, T4: Plateau, T5: East, V1: Coastal Plain, V2: Southern Piedmont, V3: Northern Piedmont, V4: Northern Mountains, V5: Southern Mountains.

Figure 2 .
Figure 2. Study area showing (A) the states of Tennessee (TN) and Virginia (VA) with red outline indicating the 1 m DEM coverage and (B) the elevation and survey units in each state.Combinations of letters and numbers represent survey units.T1: West, T2: West Central, T3: Central, T4: Plateau, T5: East, V1: Coastal Plain, V2: Southern Piedmont, V3: Northern Piedmont, V4: Northern Mountains, V5: Southern Mountains.

Figure 3 .
Figure 3. NAIP point cloud density and acquisition month.(A) Point density raster showing the number of 3D NAIP points in each 10-m grid cell.(B) A zoomed view of the point density raster.(C) A zoomed view of the NAIP imagery natural color composite.(D) NAIP imagery acquisition month.

Figure 3 .
Figure 3. NAIP point cloud density and acquisition month.(A) Point density raster showing the number of 3D NAIP points in each 10-m grid cell.(B) A zoomed view of the point density raster.(C) A zoomed view of the NAIP imagery natural color composite.(D) NAIP imagery acquisition month.

Figure 4 .
Figure 4. Scatter plot of FIA maximum tree heights (x-axis) vs. maximum point cloud heights (yaxis).Columns are broken out by state and species group (i.e., hardwood and softwood), while rows show the point cloud heights obtained for NAIP and ALS using different DEM resolutions.Plots are shaded from blue (low) to yellow (high) to indicate the density of observations that occur in places where multiple plots overlap.Blue lines are the reduced major axis regression (RMA) fits and the red dashed lines are the 1:1 lines.

Figure 4 .
Figure 4. Scatter plot of FIA maximum tree heights (x-axis) vs. maximum point cloud heights (y-axis).Columns are broken out by state and species group (i.e., hardwood and softwood), while rows show the point cloud heights obtained for NAIP and ALS using different DEM resolutions.Plots are shaded from blue (low) to yellow (high) to indicate the density of observations that occur in places where multiple plots overlap.Blue lines are the reduced major axis regression (RMA) fits and the red dashed lines are the 1:1 lines.

Figure 5 .
Figure 5.Comparison of the 1 m ALS and 1 m NAIP point cloud heights at the FIA plot locations.Columns split the dataset into the plots where ALS was collected more (left) or less than 2 years from the NAIP acquisitions, which were collected in 2018; rows break out the plots collected in TN (top) and VA (bottom).Shading from blue (low) to yellow (high) indicates the density of observations in places with multiple overlapping plots.Blue lines are the reduced major axis regression (RMA) fits and the red dashed lines are the 1:1 lines.

Figure 5 .
Figure 5.Comparison of the 1 m ALS and 1 m NAIP point cloud heights at the FIA plot locations.Columns split the dataset into the plots where ALS was collected more (left) or less than 2 years from the NAIP acquisitions, which were collected in 2018; rows break out the plots collected in TN (top) and VA (bottom).Shading from blue (low) to yellow (high) indicates the density of observations in places with multiple overlapping plots.Blue lines are the reduced major axis regression (RMA) fits and the red dashed lines are the 1:1 lines.Remote Sens. 2022, 14, 15 of 30

Figure 6 .
Figure 6.Differences between FIA maximum tree height and point cloud maximum height by FIA measurement year.Dashed red lines are drawn at +/− 5 m and dashed black lines at zero.For the boxplots, the bottom and top of the rectangles span the 25% to 75% height range, while the thick black bars are the median values.Columns break out the plots by point cloud source (ALS or NAIP) and DEM resolution (1 m and 10 m), while rows highlight the data collected in different states (TN top, VA bottom).

Figure 6 .
Figure 6.Differences between FIA maximum tree height and point cloud maximum height by FIA measurement year.Dashed red lines are drawn at +/− 5 m and dashed black lines at zero.For the boxplots, the bottom and top of the rectangles span the 25% to 75% height range, while the thick black bars are the median values.Columns break out the plots by point cloud source (ALS or NAIP) and DEM resolution (1 m and 10 m), while rows highlight the data collected in different states (TN top, VA bottom).

Figure 6 .
Figure 6.Differences between FIA maximum tree height and point cloud maximum height by FIA measurement year.Dashed red lines are drawn at +/− 5 m and dashed black lines at zero.For the boxplots, the bottom and top of the rectangles span the 25% to 75% height range, while the thick black bars are the median values.Columns break out the plots by point cloud source (ALS or NAIP) and DEM resolution (1 m and 10 m), while rows highlight the data collected in different states (TN top, VA bottom).

Figure 7 .
Figure 7.Comparison between 1 m NAIP natural color composite (left), 10 m NAIP DHM (center), and 2 m NAIP DHM (right) for softwood (A) and hardwood (B) areas in the state of TN.Figure 7. Comparison between 1 m NAIP natural color composite (left), 10 m NAIP DHM (center), and 2 m NAIP DHM (right) for softwood (A) and hardwood (B) areas in the state of TN.

Figure 9 .
Figure 9. Scatter plots showing the relationship between NAIP maximum height (in meters, top row) and NLCD tree canopy cover (in %, bottom row) vs. FIA volume per area (m 3 /ha) for hardwood (left column) and softwood (right column) species groups.

Figure 9 .
Figure 9. Scatter plots showing the relationship between NAIP maximum height (in meters, top row) and NLCD tree canopy cover (in %, bottom row) vs. FIA volume per area (m 3 /ha) for hardwood (left column) and softwood (right column) species groups.

Unit Plots Total Volume Forest Area VPA Max Height Mean Height Forest Loss Mean Elevation Major Species Group
Summary of forest attributes in each FIA survey unit.Total volume, Volume per Area (VPA), and Forest area are estimated from the FIA sample using simple random sampling.Max height represents the tallest measured tree while Mean height is the average height of all the trees measured in each survey unit.Forest loss was estimated using disturbance data from the U.S. Forest Service (USFS) Landscape Change Monitoring System (LCMS).Major species group represents the most frequently observed species group in each survey unit.See TableS1for the full list of FIA species groups.

Table 2 .
NAIP point density by acquisition month.
(A) Percentage of NAIP imagery in each acquisition month (%).(B) Percentage of 10 m pixels in each point density class (%) and the Min, Mean and Maximum point density in each state.

Table 2 .
NAIP point density by acquisition month.(A) Percentage of NAIP imagery in each acquisition month (%).Percentage of 10 m pixels in each point density class (%) and the Min, Mean and Maximum point density in each state.

Table 3 .
Summary of maximum height comparison between the NAIP point cloud data and FIA field measured trees.
FIA indicates the strata map with the highest variance ratio.Columns with asterisks (*) have values taken from the best post-stratification method while double asterisks (**) indicate survey units where the tree canopy cover map (PS TCC ) resulted in a lower variance than the NAIP strata maps.

Table 5 .
Summary of post-stratification results by survey unit.

Table 6 .
Post-stratified estimation results by FIA species group.