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.

This study examines the effectiveness of using the Normalized Difference Vegetation Index (NDVI) derived from 1326 different Landsat Thematic Mapper and Enhanced Thematic Mapper images in finding low density development within the Commonwealth of Virginia's forests. Individual NDVI images were stacked by year for the years 1995-2011 and the yearly maximum for each pixel was extracted, resulting in a 17-year image stack of all yearly maxima (a 98.7% data reduction). Using location data from housing starts and well permits, known previously forested housing starts were isolated from all other forest disturbance types. Samples from development disturbances and other forest disturbances, as well as from undisturbed forest, were used to derive vegetation index thresholds enabling separation of disturbed forest from undisturbed forest. Disturbances, once identified, could be separated into Development Disturbances and Non-Development Disturbances using a classification tree and only two variables from the Disturbance Detection and Diagnostics (D 3 ) algorithm: the maximum NDVI in the available recovery period and the slope between the NDVI value at the time of the disturbance and the maximum NDVI in the available recovery period. Low density development disturbances of previous forest land cover had an F-measure, combining precision and recall into a single class-specific accuracy (β = 1), of 0.663. We compared our results to the NLCD 2001-2011 land cover changes from any forest (classes 41, 42, 43, and 90) to any developed (classes 21, 22, 23, and 24), resulting in an F-measure of 0.00 for the same validation points. Landsat time series stacks thus show promise for identifying even the small changes associated with low density development that have been historically overlooked/underestimated by prior mapping efforts. However, further research is needed to ensure that (1) the approach will work in other forest biomes and (2) enabling detection of these important, but spatially and spectrally subtle, disturbances still ensures accurate detection of other forest disturbances.