Leaf area index (LAI) is an important biophysical parameter used to monitor, model, and manage loblolly pine plantations across the southeastern United States. Landsat provides forest scientists and managers the ability to obtain accurate and timely LAI estimates. The objective of this study was to investigate the relationship between loblolly pine LAI measured in situ (at both leaf area minimum and maximum through two growing seasons at two geographically disparate study areas) and vegetation indices calculated using data from Landsat 7 (ETM+) and Landsat 8 (OLI). Sub-objectives included examination of the impact of georegistration accuracy, comparison of top-of-atmosphere and surface reflectance, development of a new empirical model for the species and region, and comparison of the new empirical model with the current operational standard. Permanent plots for the collection of ground LAI measurements were established at two locations near Appomattox, Virginia and Tuscaloosa, Alabama in 2013 and 2014, respectively. Each plot is thirty by thirty meters in size and is located at least thirty meters from a stand boundary. Plot LAI measurements were collected twice a year using the LI-COR LAI-2200 Plant Canopy Analyzer. Ground measurements were used as dependent variables in ordinary least squares regressions with ETM+ and OLI-derived vegetation indices. We conclude that accurately-located ground LAI estimates at minimum and maximum LAI in loblolly pine stands can be combined and modeled with Landsat-derived vegetation indices using surface reflectance, particularly simple ratio (SR) and normalized difference moisture index (NDMI), across sites and sensors. The best resulting model (LAI = −0.00212 + 0.3329SR) appears not to saturate through an LAI of 5 and is an improvement over the current operational standard for loblolly pine monitoring, modeling, and management in this ecologically and economically important region.

The goal of the current study was to develop methods of estimating the height of vertical components within plantation coniferous forest using airborne discrete multiple return lidar. In the summer of 2008, airborne lidar and field data were acquired for Loblolly pine forest locations in North Carolina and Virginia, USA, which comprised a variety of stand conditions (e.g. stand age, nutrient regime, and stem density). The methods here implement both field plot-scale analysis and an automated approach for the delineation of individual tree crown (ITC) locations and horizontal extents through a marker-based region growing process applied to a lidar derived canopy height model. The estimation of vertical features was accomplished through aggregating lidar return height measurements into vertical height bins, of a given horizontal extent (plot or ITC), creating a vertical ‘stack’ of bins describing the frequency of returns by height. Once height bins were created the resulting vertical distributions were smoothed with a regression curve-line function and canopy layers were identified through the detection of local maxima and minima. Estimates from Lorey’s mean canopy height was estimated from plot-level curve-fitting with an overall accuracy of 5.9% coefficient of variation (CV) and the coefficient of determination (R²) value of 0.93. Estimates of height to the living canopy produced an overall R² value of 0.91 (11.0% CV). The presence of vertical features within the sub-canopy component of the fitted vertical function also corresponded to areas of known understory presence and absence. Estimates from ITC data were averaged to the plot level. Estimates of field Lorey’s mean canopy top height from average ITC data produced an R² value of 0.96 (7.9% CV). Average ITC estimates of height to the living canopy produced the closest correspondence to the field data, producing an R² value of 0.97 (6.2% CV). These results were similar to estimates produced by a statistical regression method, where R² values were 0.99 (2.4% CV) and 0.98 (4.9% CV) for plot average top canopy height and height to the living canopy, respectively. These results indicate that the characteristics of the dominant canopy can be estimated accurately using airborne lidar without the development of regression models, in a variety of intensively managed coniferous stand conditions.

Airborne discrete-return light detection and ranging (lidar) can be used to estimate leaf area index (LAI) with relatively high accuracy. This capacity was explored with regard to assessing the capability of estimating LAI at different heights at the plot level, in the presence of understorey vegetation, within intensively managed Loblolly pine forest in North Carolina, USA. Field measurements utilized the LI-COR LAI-2200 plant canopy analyser for field-based estimates of effective LAI at three elevations within each plot; these were on the ground (0.0 m) and 1.0 m and 2.5 m above the ground within the various understorey heights and densities. A number of new and previously existing lidar metrics and indices were calculated from the distribution of return heights, which have been identified as potentially strong predictors of LAI. A bivariate and stepwise regression approach was then applied to create models for the estimation of LAI from lidar-derived height distribution metrics. The results show that specific logarithm transformed laser penetration indices calculated using a height threshold (e.g. the number of returns below 2.5 m ratioed against all returns) as close to field LAI measurement height (e.g. 2.5 m) was more effective than other lidar metrics. LAI can be estimated for each of the three measurement heights within the understorey component explaining 67 to 76% of the variance (root mean square error 0.42–0.57). The indices that produced the highest correlations and which were selected in stepwise regression analysis were calculated using all returns. The results indicate that LAI can be estimated accurately using lidar data in pine plantation forest over a variety of stand conditions.

Airborne laser scanning (ALS) is increasingly used to estimate various forest characteristics. Technological improvements in unmanned aerial vehicles (UAVs) and drone laser scanning (DLS) sensors have permitted the acquisition of high pulse density datasets. There is an assumption that higher pulse densities yield greater accuracies in estimating forest characteristics. In this study, we investigated the effect of pulse density (.25, .5, 1, 5, 10, 50, 100 and 300 pulses m⁻²) on the ability to delineate individual tree crowns (ITCs) and estimate ITC height and crown horizontal diameter, in addition to plot-level leaf area index (LAI). The current study took place in an experimentally varied Pinus taeda L. forest, which included three stem densities: (i) 618; (ii) 1236; and (iii) 1853 trees per hectare (TPH). ITCs were classified directly from the DLS point cloud for each of the pulse densities. The correct delineation of ITCs relative to field tree-coordinates was relatively consistent (±5%) for pulse densities of 5 to 300 pulses m⁻². ITC delineation accuracy decreased with lower pulse densities. Planting stem density did impact ITC delineation accuracy. Higher pulse densities, plots with 618 TPH correctly classified ~88% of ITCs, and plots with the 1853 TPH correctly classified ~50% of ITCs. Estimates of tree height were largely unaffected by changes in tree density. Root mean square error (RMSE) for tree height varied from .5 to 2.5 m at pulse densities of 300 to .25 pulses m⁻², respectively. Estimates of crown horizontal diameter varied with regard to both pulse and stem density from 1.2 (300 ppm⁻² and 1853 TPH) to 4.2 m (.25 ppm⁻² and 618 TPH). RMSE varied among stem densities from .6 to 1.2 m as pulse density decreased. There was significant difference in ITC delineation accuracy, particularly when considering stem density, and the estimates of tree height and crown horizontal diameter among the DLS pulse densities used. The accuracy of predicted LAI was largely unaffected by changes in pulse density, when pulse density was above .5 pulses m⁻². There was little or no difference in estimates of LAI at these pulse densities. Our results suggest that low-density DLS data may be capable of estimating plot-level forest metrics reliably in some situations, however once the analysis scale is reduced to the individual-tree-level, the influence of pulse density is more substantial. The results here provide guidance to forest managers who must balance metric estimation accuracy and price when planning new ALS or DLS acquisitions.

The competitive neighbourhood surrounding an individual tree can have a significant influence on its diameter at breast height (DBH) and individual tree stem volume (SV). Distance dependent competition index metrics are rarely recorded in traditional field campaigns because they are laborious to collect and are spatially limited. Remote sensing data could overcome these limitations while providing estimation of forest attributes over a large area. We used unoccupied aerial vehicle laser scanning data to delineate individual tree crowns (ITCs) and calculated crown size and distance-dependent competition indices to estimate DBH and SV. We contrasted two methods: (i) Random Forest (RF) and (ii) backwards-stepwise, linear multiple regression (LMR). We utilized an existing experiment in Pinus taeda L. plantations including multiple planting densities, genotypes and silvicultural levels. While the tree planting density did affect the correct delineation of ITCs, between 61% and 99% (mean 86%) were correctly linked to the planting location. The most accurate RF and LMR models all included metrics related to ITC size and competitive neighbourhood. The DBH estimates from RF and LMR were similar: RMSE 3.05 and 3.13 cm (R² 0.64 and 0.62), respectively. Estimates of SV from RF were slightly better than for LMR: RMSE 0.06 and 0.07 m³ (R² 0.77 and 0.70), respectively. Our results provide evidence that ITC size and competition index metrics may improve DBH and SV estimation accuracy when analysing laser-scanning data. The ability to provide accurate, and near-complete, forest inventories holds a great deal of potential for forest management planning.

The tree crown’s directionality of growth may be an indicator of how aggressive the tree is in terms of foraging for light. Airborne drone laser scanning (DLS) has been used to accurately classify individual tree crowns (ITCs) and derive size metrics related to the crown. We compare ITCs among 6 genotypes exhibiting different crown architectures in managed loblolly pine (Pinus taeda L.) in the United States. DLS data are classified into ITC objects, and we present novel methods to calculate ITC shape metrics. Tree stems are located using (a) model-based clustering and (b) weighting cluster-based size. We generated ITC shape metrics using 3-dimensional (3D) alphashapes in 2 DLS acquisitions of the same location, 4 years apart. Crown horizontal distance from the stem was estimated at multiple heights, in addition to calculating 3D volume in specific azimuths. Crown morphologies varied significantly (P < 0.05) spatially, temporally, and among the 6 genotypes. Most genotypes exhibited larger crown volumes facing south (150° to 173°). We found that crown asymmetries were consistent with (a) the direction of solar radiation, (b) the spatial arrangement and proximity of the neighboring crowns, and (c) genotype. Larger crowns were consistent with larger increases in stem volume, but that increases in the southern portions of crown volume were consistent with larger stem volume increases, than in the north. This finding suggests that row orientation could influence stem growth rates in plantations, particularly impacting earlier development. These differences can potentially reduce over time, especially if stands are not thinned in a timely manner once canopy growing space has diminished.