Search

Accurately predicting soil respiration (Rs) has received considerable attention recently due to its importance as a significant carbon flux back to the atmosphere. Even small changes in Rs can have a significant impact on the net ecosystem productivity of forests. Variations in Rs have been related to both spatial and temporal variation due to changes in both abiotic and biotic factors. This study focused on soil temperature and moisture and changes in the species composition of the overstory and understory and how these variables impact Rs. Sample plots consisted of four vegetation types: eastern hemlock (Tsuga canadensis L. Carriere)-dominated overstory, mountain laurel (Kalmia latifolia L.)-dominated understory, hardwood-dominated overstory, and cinnamon fern (Osmundastrum cinnamomeum (L.) C. Presl)-dominated understory, with four replications of each. Remotely sensed data collected for each plot, light detection and ranging, and hyperspectral data, were compiled from the National Ecological Observatory Network (NEON) to determine if they could improve predictions of Rs. Soil temperature and soil moisture explained 82% of the variation in Rs. There were no statistically significant differences between the average annual Rs rates among the vegetation types. However, when looking at monthly Rs, cinnamon fern plots had statistically higher rates in the summer when it was abundant and hemlock had significantly higher rates in the dormant months. At the same soil temperature, the vegetation types’ Rs rates were not statistically different. However, the cinnamon fern plots showed the most sensitivity to soil moisture changes and were the wettest sites. Normalized Difference Lignin Index (NDLI) was the only vegetation index (VI) to vary between the vegetation types. It also correlated with Rs for the months of August and September. Photochem

Patterns in foliar nitrogen (N) stable isotope ratios (δ15N) have been shown to reveal trends in terrestrial N cycles, including the identification of ecosystems where N deficiencies limit forest ecosystem productivity. However, there is a gap in our understanding of within-species variation and species-level response to environmental gradients or forest management. Our objective is to examine the relationship between site index, foliar %N, foliar δ15N and spectral reflectance for managed Douglas-fir (Pseudotsuga menziesii) and loblolly pine (Pinus taeda) plantations across their geographic ranges in the Pacific Northwest and the southeastern United States, respectively. Foliage was measured at 28 sites for reflectance using a handheld spectroradiometer, and further analyzed for δ15N and N concentration. Unlike the prior work for grasslands and shrubland species, our results show that foliar δ15N and foliar %N are not well correlated for these tree species. However, multiple linear regression models suggest a strong predictive ability of spectroscopy data to quantify foliar δ15N, with some models explaining more than 65% of the variance in the δ15N. Additionally, moderate to strong explanations of variance were found between site index and foliar δ15N (R2 = 0.49) and reflectance and site index (R2 = 0.84) in the Douglas-fir data set. The development of relationships between foliar spectral reflectance, δ15N and measures of site productivity provides the first step toward mapping canopy δ15N for these managed forests with remote sensing.

Consumer-grade drone-produced digital orthoimagery is a valuable tool for conservation management and enables the low-cost monitoring of remote ecosystems. This study demonstrates the applicability of RGB orthoimagery for the assessment of forest health at the scale of individual trees in a 46-hectare plot of rare southern Appalachian red spruce forest on Whitetop Mountain, Virginia. We used photogrammetric Structure from Motion software Pix4Dmapper with drone-collected imagery to generate a mosaic for point cloud reconstruction and orthoimagery of the plot. Using 3-band RBG digital orthoimagery, we visually classified 9402 red spruce individuals, finding 8700 healthy (92.5%), 251 declining/dying (2.6%), and 451 dead (4.8%). We mapped individual spruce trees in each class and produced kernel density maps of health classes (live, dead, and dying). Our approach provided a nearly gap-free assessment of the red spruce canopy in our study site, versus a much more time-intensive field survey. Our maps provided useful information on stand mortality patterns and canopy gaps that could be used by managers to identify optimal locations for selective thinning to facilitate understory sapling regeneration. This approach, dependent mainly on an off-the-shelf drone system and visual interpretation of orthoimagery, could be applied by land managers to measure forest health in other spruce, or possibly spruce-fir, communities in the Appalachians. Our study highlights the usefulness of drone-produced orthoimagery for conservation monitoring, presenting a valid and accessible protocol for the monitoring and assessment of forest health in remote spruce, and possibly other conifer, populations. Adoption of drone-based monitoring may be especially useful in light of climate change and the possible displacement of southern Appalachian red spruce (and spruce-fir) ecosystems by the upslope migration of deciduous trees.

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.

Sustainable forest management is hugely dependent on high-quality estimates of forest site productivity, but it is challenging to generate productivity maps over large areas. We present a method for generating site index (a measure of such forest productivity) maps for plantation loblolly pine (Pinus taeda L.) forests over large areas in the southeastern United States by combining airborne laser scanning (ALS) data from disparate acquisitions and Landsat-based estimates of forest age. For predicting canopy heights, a linear regression model was developed using ALS data and field measurements from the Forest Inventory and Analysis (FIA) program of the US Forest Service (n = 211 plots). The model was strong (R2 = 0.84, RMSE = 1.85 m), and applicable over a large area (~208,000 sq. km). To estimate the site index, we combined the ALS estimated heights with Landsat-derived maps of stand age and planted pine area. The estimated bias was low (−0.28 m) and the RMSE (3.8 m, relative RMSE: 19.7%, base age 25 years) was consistent with other similar approaches. Due to Landsat-related constraints, our methodology is valid only for relatively young pine plantations established after 1984. We generated 30 m resolution site index maps over a large area (~832 sq. km). The site index distribution had a median value of 19.4 m, the 5th percentile value of 13.0 m and the 95th percentile value of 23.3 m. Further, using a watershed level analysis, we ranked these regions by their estimated productivity. These results demonstrate the potential and value of remote sensing based large-area site index maps.

Sustainable forest management is hugely dependent on high-quality estimates of forest site productivity, but it is challenging to generate productivity maps over large areas. We present a method for generating site index (a measure of such forest productivity) maps for plantation loblolly pine (Pinus taeda L.) forests over large areas in the southeastern United States by combining airborne laser scanning (ALS) data from disparate acquisitions and Landsat-based estimates of forest age. For predicting canopy heights, a linear regression model was developed using ALS data and field measurements from the Forest Inventory and Analysis (FIA) program of the US Forest Service (n = 211 plots). The model was strong (R2 = 0.84, RMSE = 1.85 m), and applicable over a large area (~208,000 sq. km). To estimate the site index, we combined the ALS estimated heights with Landsat-derived maps of stand age and planted pine area. The estimated bias was low (−0.28 m) and the RMSE (3.8 m, relative RMSE: 19.7%, base age 25 years) was consistent with other similar approaches. Due to Landsat-related constraints, our methodology is valid only for relatively young pine plantations established after 1984. We generated 30 m resolution site index maps over a large area (~832 sq. km). The site index distribution had a median value of 19.4 m, the 5th percentile value of 13.0 m and the 95th percentile value of 23.3 m. Further, using a watershed level analysis, we ranked these regions by their estimated productivity. These results demonstrate the potential and value of remote sensing based large-area site index maps.

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 continued development of algorithms using multitemporal Landsat data creates opportunities to develop and adapt imputation algorithms to improve the quality of that data as part of preprocessing. One example is de-striping Enhanced Thematic Mapper Plus (ETM+, Landsat 7) images acquired after the Scan Line Corrector failure in 2003. In this study, we apply window regression, an algorithm that was originally designed to impute low-quality Moderate Resolution Imaging Spectroradiometer (MODIS) data, to Landsat Analysis Ready Data from 2014–2016. We mask Operational Land Imager (OLI; Landsat 8) image stacks from five study areas with corresponding ETM+ missing data layers, using these modified OLI stacks as inputs. We explored the algorithm’s parameter space, particularly window size in the spatial and temporal dimensions. Window regression yielded the best accuracy (and moderately long computation time) with a large spatial radius (a 7 × 7 pixel window) and a moderate temporal radius (here, five layers). In this case, root mean square error for deviations from the observed reflectance ranged from 3.7–7.6% over all study areas, depending on the band. Second-order response surface analysis suggested that a 15 × 15 pixel window, in conjunction with a 9-layer temporal window, may produce the best accuracy. Compared to the neighborhood similar pixel interpolator gap-filling algorithm, window regression yielded slightly better accuracy on average. Because it relies on no ancillary data, window regression may be used to conveniently preprocess stacks for other data-intensive algorithms.

Understanding irrigator responses to changes in water availability is critical for building strategies to support effective management of water resources. Using remote sensing data, we examine farmer responses to seasonal changes in water availability in Idaho’s Snake River Plain for the time series 1984–2016. We apply a binary threshold based on the seasonal maximum of the Normalized Difference Moisture Index (NDMI) using Landsat 5–8 images to distinguish irrigated from non-irrigated lands. We find that the NDMI of irrigated lands increased over time, consistent with trends in irrigation technology adoption and increased crop productivity. By combining remote sensing data with geospatial data describing water rights for irrigation, we show that the trend in NDMI is not universal, but differs by farm size and water source. Farmers with small farms that rely on surface water are more likely than average to have a large contraction (over −25%) in irrigated area over the 33-year period of record. In contrast, those with large farms and access to groundwater are more likely than average to have a large expansion (over +25%) in irrigated area over the same period.

The continued development of algorithms using multitemporal Landsat data creates opportunities to develop and adapt imputation algorithms to improve the quality of that data as part of preprocessing. One example is de-striping Enhanced Thematic Mapper Plus (ETM+, Landsat 7) images acquired after the Scan Line Corrector failure in 2003. In this study, we apply window regression, an algorithm that was originally designed to impute low-quality Moderate Resolution Imaging Spectroradiometer (MODIS) data, to Landsat Analysis Ready Data from 2014–2016. We mask Operational Land Imager (OLI; Landsat 8) image stacks from five study areas with corresponding ETM+ missing data layers, using these modified OLI stacks as inputs. We explored the algorithm’s parameter space, particularly window size in the spatial and temporal dimensions. Window regression yielded the best accuracy (and moderately long computation time) with a large spatial radius (a 7 × 7 pixel window) and a moderate temporal radius (here, five layers). In this case, root mean square error for deviations from the observed reflectance ranged from 3.7–7.6% over all study areas, depending on the band. Second-order response surface analysis suggested that a 15 × 15 pixel window, in conjunction with a 9-layer temporal window, may produce the best accuracy. Compared to the neighborhood similar pixel interpolator gap-filling algorithm, window regression yielded slightly better accuracy on average. Because it relies on no ancillary data, window regression may be used to conveniently preprocess stacks for other data-intensive algorithms.

Understanding irrigator responses to changes in water availability is critical for building strategies to support effective management of water resources. Using remote sensing data, we examine farmer responses to seasonal changes in water availability in Idaho’s Snake River Plain for the time series 1984–2016. We apply a binary threshold based on the seasonal maximum of the Normalized Difference Moisture Index (NDMI) using Landsat 5–8 images to distinguish irrigated from non-irrigated lands. We find that the NDMI of irrigated lands increased over time, consistent with trends in irrigation technology adoption and increased crop productivity. By combining remote sensing data with geospatial data describing water rights for irrigation, we show that the trend in NDMI is not universal, but differs by farm size and water source. Farmers with small farms that rely on surface water are more likely than average to have a large contraction (over −25%) in irrigated area over the 33-year period of record. In contrast, those with large farms and access to groundwater are more likely than average to have a large expansion (over +25%) in irrigated area over the same period.

In the published paper [1], the title and Appendix Tables A4, A5, A7, and A8 contain typographical errors. The correct title and table captions are as follows: [...]

There are pressing concerns about the interplay between agricultural productivity, water demand, and water availability in semi-arid to arid regions of the world. Currently, irrigated agriculture is the dominant water user in these regions and is estimated to consume approximately 80% of the world’s diverted freshwater resources. We develop an improved irrigated land-use mapping algorithm that uses the seasonal maximum value of a spectral index to distinguish between irrigated and non-irrigated parcels in Idaho’s Snake River Plain. We compare this approach to two alternative algorithms that differentiate between irrigated and non-irrigated parcels using spectral index values at a single date or the area beneath spectral index trajectories for the duration of the agricultural growing season. Using six different pixel and county-scale error metrics, we evaluate the performance of these three algorithms across all possible combinations of two growing seasons (2002 and 2007), two datasets (MODIS and Landsat 5), and three spectral indices, the Normalized Difference Vegetation Index, Enhanced Vegetation Index and Normalized Difference Moisture Index (NDVI, EVI, and NDMI). We demonstrate that, on average, the seasonal-maximum algorithm yields an improvement in classification accuracy over the accepted single-date approach, and that the average improvement under this approach is a 60% reduction in county scale root mean square error (RMSE), and modest improvements of overall accuracy in the pixel scale validation. The greater accuracy of the seasonal-maximum algorithm is primarily due to its ability to correctly classify non-irrigated lands in riparian and developed areas of the study region.

Remote detection of forest disturbance remains a key area of interest for scientists and land managers. Subtle disturbances such as drought, disease, insect activity, and thinning harvests have a significant impact on carbon budgeting and forest productivity, but current change detection algorithms struggle to accurately identify them, especially over decadal timeframes. We introduce an algorithm called Edyn, which inputs a time series of residuals from harmonic regression into a control chart to signal low-magnitude, consistent deviations from the curve as disturbances. After signaling, Edyn retrains a new baseline curve. We compared Edyn with its parent algorithm (EWMACD—Exponentially Weighted Moving Average Change Detection) on over 3500 visually interpreted Landsat pixels from across the contiguous USA, with reference data for timing and type of disturbance. For disturbed forested pixels, Edyn had a mean per-pixel commission error of 31.1% and omission error of 70.0%, while commission and omission errors for EWMACD were 39.9% and 65.2%, respectively. Edyn had significantly less overall error than EWMACD (F1 = 0.19 versus F1 = 0.13). These patterns generally held for all of the reference data, including a direct comparison to other contemporary change detection algorithms, wherein Edyn and EWMACD were found to have lower omission error rates for a category of subtle changes over long periods.

In the published paper [1], the title and Appendix Tables A4, A5, A7, and A8 contain typographical errors. The correct title and table captions are as follows: [...]

Remote detection of forest disturbance remains a key area of interest for scientists and land managers. Subtle disturbances such as drought, disease, insect activity, and thinning harvests have a significant impact on carbon budgeting and forest productivity, but current change detection algorithms struggle to accurately identify them, especially over decadal timeframes. We introduce an algorithm called Edyn, which inputs a time series of residuals from harmonic regression into a control chart to signal low-magnitude, consistent deviations from the curve as disturbances. After signaling, Edyn retrains a new baseline curve. We compared Edyn with its parent algorithm (EWMACD—Exponentially Weighted Moving Average Change Detection) on over 3500 visually interpreted Landsat pixels from across the contiguous USA, with reference data for timing and type of disturbance. For disturbed forested pixels, Edyn had a mean per-pixel commission error of 31.1% and omission error of 70.0%, while commission and omission errors for EWMACD were 39.9% and 65.2%, respectively. Edyn had significantly less overall error than EWMACD (F1 = 0.19 versus F1 = 0.13). These patterns generally held for all of the reference data, including a direct comparison to other contemporary change detection algorithms, wherein Edyn and EWMACD were found to have lower omission error rates for a category of subtle changes over long periods.

There are pressing concerns about the interplay between agricultural productivity, water demand, and water availability in semi-arid to arid regions of the world. Currently, irrigated agriculture is the dominant water user in these regions and is estimated to consume approximately 80% of the world’s diverted freshwater resources. We develop an improved irrigated land-use mapping algorithm that uses the seasonal maximum value of a spectral index to distinguish between irrigated and non-irrigated parcels in Idaho’s Snake River Plain. We compare this approach to two alternative algorithms that differentiate between irrigated and non-irrigated parcels using spectral index values at a single date or the area beneath spectral index trajectories for the duration of the agricultural growing season. Using six different pixel and county-scale error metrics, we evaluate the performance of these three algorithms across all possible combinations of two growing seasons (2002 and 2007), two datasets (MODIS and Landsat 5), and three spectral indices, the Normalized Difference Vegetation Index, Enhanced Vegetation Index and Normalized Difference Moisture Index (NDVI, EVI, and NDMI). We demonstrate that, on average, the seasonal-maximum algorithm yields an improvement in classification accuracy over the accepted single-date approach, and that the average improvement under this approach is a 60% reduction in county scale root mean square error (RMSE), and modest improvements of overall accuracy in the pixel scale validation. The greater accuracy of the seasonal-maximum algorithm is primarily due to its ability to correctly classify non-irrigated lands in riparian and developed areas of the study region.

We recruit an online labor force through Amazon.com’s Mechanical Turk platform to identify clouds and cloud shadows in Landsat satellite images. We find that a large group of workers can be mobilized quickly and relatively inexpensively. Our results indicate that workers’ accuracy is insensitive to wage, but deteriorates with the complexity of images and with time-on-task. In most instances, human interpretation of cloud impacted area using a majority rule was more accurate than an automated algorithm (Fmask) commonly used to identify clouds and cloud shadows. However, cirrus-impacted pixels were better identified by Fmask than by human interpreters. Crowd-sourced interpretation of cloud impacted pixels appears to be a promising means by which to augment or potentially validate fully automated algorithms.

We recruit an online labor force through Amazon.com’s Mechanical Turk platform to identify clouds and cloud shadows in Landsat satellite images. We find that a large group of workers can be mobilized quickly and relatively inexpensively. Our results indicate that workers’ accuracy is insensitive to wage, but deteriorates with the complexity of images and with time-on-task. In most instances, human interpretation of cloud impacted area using a majority rule was more accurate than an automated algorithm (Fmask) commonly used to identify clouds and cloud shadows. However, cirrus-impacted pixels were better identified by Fmask than by human interpreters. Crowd-sourced interpretation of cloud impacted pixels appears to be a promising means by which to augment or potentially validate fully automated algorithms.

Information on foliar macronutrients is required in order to understand plant physiological and ecosystem processes such as photosynthesis, nutrient cycling, respiration and cell wall formation. The ability to measure, model and map foliar macronutrients (nitrogen (N), phosphorus (P), potassium (K), calcium (Ca) and magnesium (Mg)) at the forest canopy level provides information on the spatial patterns of ecosystem processes (e.g., carbon exchange) and provides insight on forest condition and stress. Imaging spectroscopy (IS) has been used particularly for modeling N, using airborne and satellite imagery mostly in temperate and tropical forests. However, there has been very little research conducted at these scales to model P, K, Ca, and Mg and few studies have focused on boreal forests. We report results of a study of macronutrient modeling using spaceborne IS and airborne light detection and ranging (LiDAR) data for a mixedwood boreal forest canopy in northern Ontario, Canada. Models incorporating Hyperion data explained approximately 90% of the variation in canopy concentrations of N, P, and Mg; whereas the inclusion of LiDAR data significantly improved the prediction of canopy concentration of Ca (R2 = 0.80). The combined used of IS and LiDAR data significantly improved the prediction accuracy of canopy Ca and K concentration but decreased the prediction accuracy of canopy P concentration. The results indicate that the variability of macronutrient concentration due to interspecific and functional type differences at the site provides the basis for the relationship observed between the remote sensing measurements (i.e., IS and LiDAR) and macronutrient concentration. Crown closure and canopy height are the structural metrics that establish the connection between macronutrient concentration and IS and LiDAR data, respectively. The spatial distribution of macronutrient concentration at the canopy scale mimics functional type distribution at the sit

Generating accurate and unbiased wall-to-wall canopy height maps from airborne lidar data for large regions is useful to forest scientists and natural resource managers. However, mapping large areas often involves using lidar data from different projects, with varying acquisition parameters. In this work, we address the important question of whether one can accurately model canopy heights over large areas of the Southeastern US using a very heterogeneous dataset of small-footprint, discrete-return airborne lidar data (with 76 separate lidar projects). A unique aspect of this effort is the use of nationally uniform and extensive field data (~1800 forested plots) from the Forest Inventory and Analysis (FIA) program of the US Forest Service. Preliminary results are quite promising: Over all lidar projects, we observe a good correlation between the 85th percentile of lidar heights and field-measured height (r = 0.85). We construct a linear regression model to predict subplot-level dominant tree heights from distributional lidar metrics (R2 = 0.74, RMSE = 3.0 m, n = 1755). We also identify and quantify the importance of several factors (like heterogeneity of vegetation, point density, the predominance of hardwoods or softwoods, the average height of the forest stand, slope of the plot, and average scan angle of lidar acquisition) that influence the efficacy of predicting canopy heights from lidar data. For example, a subset of plots (coefficient of variation of vegetation heights <0.2) significantly reduces the RMSE of our model from 3.0–2.4 m (~20% reduction). We conclude that when all these elements are factored into consideration, combining data from disparate lidar projects does not preclude robust estimation of canopy heights.

The need for an efficient and standard technique for optimal spectral sampling of hyperspectral data during the inversion of canopy reflectance models has been the subject of many studies. The objective of this study was to investigate the utility of the discrete wavelet transform (DWT) for extracting useful features from hyperspectral data with which forest LAI can be estimated through inversion of a three dimensional radiative transfer model, the Discrete Anisotropy Radiative Transfer (DART) model. DART, coupled with the leaf optical properties model PROSPECT, was inverted with AVIRIS data using a look-up-table (LUT)-based inversion approach. We used AVIRIS data and in situ LAI measurements from two different hardwood forested sites in Wisconsin, USA. Prior to inversion, model-simulated and AVIRIS hyperspectral data were transformed into discrete wavelet coefficients using Haar wavelets. The LUT inversion was performed with three different datasets, the original reflectance bands, the full set of wavelet extracted features, and two wavelet subsets containing 99.99% and 99.0% of the cumulative energy of the original signal. The energy subset containing 99.99% of the cumulative signal energy provided better estimates of LAI (RMSE = 0.46, R2 = 0.77) than the original spectral bands (RMSE = 0.60, R2 = 0.47). The results indicate that the discrete wavelet transform can increase the accuracy of LAI estimates by improving the LUT-based inversion of DART (and, potentially, by implication, other terrestrial radiative transfer models) using hyperspectral data. The improvement in accuracy of LAI estimates is potentially due to different properties of wavelet analysis such as multi-scale representation, dimensionality reduction, and noise removal.

The objective of this study was to determine whether leaf area index (LAI) in temperate mixed forests is best estimated using multiple-return airborne laser scanning (lidar) data or dual-band, single-pass interferometric synthetic aperture radar data (from GeoSAR) alone, or both in combination. In situ measurements of LAI were made using the LiCor LAI-2000 Plant Canopy Analyzer on 61 plots (21 hardwood, 36 pine, 4 mixed pine hardwood; stand age ranging from 12-164 years; mean height ranging from 0.4 to 41.2 m) in the Appomattox-Buckingham State Forest, Virginia, USA. Lidar distributional metrics were calculated for all returns and for ten one meter deep crown density slices (a new metric), five above and five below the mode of the vegetation returns for each plot. GeoSAR metrics were calculated from the X-band backscatter coefficients (four looks) as well as both X- and P-band interferometric heights and magnitudes for each plot. Lidar metrics alone explained 69% of the variability in LAI, while GeoSAR metrics alone explained 52%. However, combining the lidar and GeoSAR metrics increased the R2 to 0.77 with a CV-RMSE of 0.42. This study indicates the clear potential for X-band backscatter and interferometric height (both now available from spaceborne sensors), when combined with small-footprint lidar data, to improve LAI estimation in temperate mixed forests.