Forest canopy height and carbon estimation at Monks Wood National Nature Reserve, UK, using dual-wavelength SAR interferometry

Forest canopy height is a critical parameter in better quantifying the terrestrial carbon cycle. It can be used to estimate aboveground biomass and carbon pools stored in the vegetation, and predict timber yield for forest management. Polarimetric SAR interferometry (PolInSAR) uses polarimetric separation of scattering phase centers derived from interferometry to estimate canopy height. A limitation of PolInSAR is that it relies on sufficient scattering phase center separation at each pixel to be able to derive accurate forest canopy height estimates. The effect of wavelength-dependent penetration depth into the canopy is known to be strong, and could potentially lead to a better height separation than relying on polarization combinations at one wavelength alone. Here we present a new method for canopy height mapping using dual-wavelength SAR interferometry (InSAR) at X- and L-band. The method is based on the scattering phase center separation at different wavelengths. It involves the generation of a smoothed interpolated terrain elevation model underneath the forest canopy from repeat-pass L-band InSAR data. The terrain model is then used to remove the terrain component from the single-pass X-band interferometric surface height to estimate forest canopy height. The ability of L-band to map terrain height under vegetation relies on sufficient spatial heterogeneity of the density of scattering elements that scatter L-band electromagnetic waves within each resolution cell. The method is demonstrated with airborne X-band VV polarized single-pass and L-band HH polarized repeat-pass SAR interferometry using data acquired by the E-SAR sensor over Monks Wood National Nature Reserve, UK. This is one of the first radar studies of a semi-natural deciduous woodland that exhibits considerable spatial heterogeneity of vegetation type and density. The canopy height model is validated using airborne imaging LIDAR data acquired by the Environment Agency. The rmse of the LIDAR canopy height estimates compared to theodolite data is 2.15 m (relative error 17.6%). The rmse of the dual-wavelength InSAR-derived canopy height model compared to LIDAR is 3.49 m (relative error 28.5%). From the canopy height maps carbon pools are estimated using allometric equations. The results are compared to a field survey of carbon pools and rmse values are presented. The dual-wavelength InSAR method could potentially be delivered from a spaceborne constellation similar to the TerraSAR system.     

An airborne lidar sampling strategy to model forest canopy height from Quickbird imagery and GEOBIA

High-resolution digital canopy models derived from airborne lidar data have the ability to provide detailed information on the vertical structure of forests. However, compared to satellite data of similar spatial resolution and extent, the small footprint airborne lidar data required to produce such models remain expensive. In an effort to reduce these costs, the primary objective of this paper is to develop an airborne lidar sampling strategy to model full-scene forest canopy height from optical imagery, lidar transects and Geographic Object-Based Image Analysis (GEOBIA). To achieve this goal, this research focuses on (i) determining appropriate lidar transect features (i.e., location, direction and extent) from an optical scene, (ii) developing a mechanism to model forest canopy height for the full-scene based on a minimum number of lidar transects, and (iii) defining an optimal mean object size (MOS) to accurately model the canopy composition and height distribution. Results show that (i) the transect locations derived from our optimal lidar transect selection algorithm accurately capture the canopy height variability of the entire study area; (ii) our canopy height estimation models have similar performance in two lidar transect directions (i.e., north-south and west-east); (iii) a small lidar extent (17.6% of total size) can achieve similar canopy height estimation accuracies as those modeled from the full lidar scene; and (iv) different MOS can lead to distinctly different canopy height results. By comparing the best canopy height estimate with the full lidar canopy height data, we obtained average estimation errors of 6.0 m and 6.8 m for conifer and deciduous forests at the individual tree crown/small tree cluster level, and an area weighted combined error of 6.2 m, which is lower than the provincial forest inventory height class interval (i.e., ≈ 9.0 m).     

Remotely sensed measurements of forest structure and fuel loads in the Pinelands of New Jersey

We used a single-beam, first return profiling LIDAR (Light Detection and Ranging) measurements of canopy height, intensive biometric measurements in plots, and Forest Inventory and Analysis (FIA) data to quantify forest structure and ladder fuels (defined as vertical fuel continuity between the understory and canopy) in the New Jersey Pinelands. The LIDAR data were recorded at 400 Hz over three intensive areas of 1 km² where transects were spaced at 200 m, and along 64 transects spaced 1 km apart (total of ca. 2500 km²). LIDAR and field measurements of canopy height were similar in the three intensive study areas, with the 80th percentile of LIDAR returns explaining the greatest amount of variability (79%). Correlations between LIDAR data and aboveground tree biomass measured in the field were highly significant when all three 1 km² areas were analyzed collectively, with the 80th percentile again explaining the greatest amount of variability (74%). However, when intensive areas were analyzed separately, correlations were poor for Oak/Pine and Pine/Scrub Oak stands. Similar results were obtained using FIA data; at the landscape scale, mean canopy height was positively correlated with aboveground tree biomass, but when forest types were analyzed separately, correlations were significant only for some wetland forests (Pitch Pine lowlands and mixed hardwoods; r² = 0.74 and 0.59, respectively), and correlations were poor for upland forests (Oak/Pine, Pine/Oak and Pine/Scrub Oak, r² = 0.33, 0.11 and 0.21, respectively). When LIDAR data were binned into 1-m height classes, more LIDAR pulses were recorded from the lowest height classes in stands with greater shrub biomass, and significant differences were detected between stands where recent prescribed fire treatments had been conducted and unburned areas. Our research indicates that single-beam LIDAR can be used for regional-scale (forest biomass) estimates, but that relationships between height and biomass can be poorer at finer scales within individual forest types. Binned data are useful for estimating the presence of ladder fuels (vertical continuity of leaves and branches) and horizontal fuel continuity below the canopy.     

Exploring LiDAR-RaDAR synergy—predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR

Scanning Light Detecting and Ranging (LiDAR), Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) were analyzed to determine (1) which of the three sensor systems most accurately predicted forest biomass, and (2) if LiDAR and SAR/InSAR data sets, jointly considered, produced more accurate, precise results relative to those same data sets considered separately. LiDAR ranging measurements, VHF-SAR cross-sectional returns, and X- and P-band cross-sectional returns and interferometric ranges were regressed with ground-estimated (from dbh) forest biomass in ponderosa pine forests in the southwestern United States. All models were cross-validated. Results indicated that the average canopy height measured by the scanning LiDAR produced the best predictive equation. The simple linear LiDAR equation explained 83% of the biomass variability (n = 52 plots) with a cross-validated root mean square error of 26.0 t/ha. Additional LiDAR metrics were not significant to the model. The GeoSAR P-band (λ = 86 cm) cross-sectional return and the GeoSAR/InSAR canopy height (X-P) captured 30% of the forest biomass variation with an average predictive error of 52.5 t/ha. A second RaDAR-FOPEN collected VHF (λ ∼ 7.8 m) and cross-polarized P-band (λ = 88 cm) cross-sectional returns, none of which proved useful for forest biomass estimation (cross-validated R² = 0.09, RMSE = 63.7 t/ha). Joint consideration of LiDAR and RaDAR measurements produced a statistically significant, albeit small improvement in biomass estimation precision. The cross-validated R² increased from 83% to 84% and the prediction error decreased from 26.0 t/ha to 24.9 t/ha when the GeoSAR X-P interferometric height is considered along with the average LiDAR canopy height. Inclusion of a third LiDAR metric, the 60th decile height, further increased the R² to 85% and decreased the RMSE to 24.1 t/ha. On this 11 km² ponderosa pine study area, LiDAR data proved most useful for predicting forest biomass. RaDAR ranging measurements did not improve the LiDAR estimates.     

A LiDAR-derived canopy density model for tree stem and crown mapping in Australian forests

The retrieval of tree and forest structural attributes from Light Detection and Ranging (LiDAR) data has focused largely on utilising canopy height models, but these have proved only partially useful for mapping and attributing stems in complex, multi-layered forests. As a complementary approach, this paper presents a new index, termed the Height-Scaled Crown Openness Index (HSCOI), which provides a quantitative measure of the relative penetration of LiDAR pulses into the canopy. The HSCOI was developed from small footprint discrete return LiDAR data acquired over mixed species woodlands and open forests near Injune, Queensland, Australia, and allowed individual trees to be located (including those in the sub-canopy) and attributed with height using relationships (r² = 0.81, RMSE = 1.85 m, n = 115; 4 outliers removed) established with field data. A threshold contour of the HSCOI surface that encompassed ∼ 90% of LiDAR vegetation returns also facilitated mapping of forest areas, delineation of tree crowns and clusters, and estimation of canopy cover. At a stand level, tree density compared well with field measurements (r² = 0.82, RMSE = 133 stems ha^− 1, n = 30), with the most consistent results observed for stem densities ≤ 700 stems ha^− 1. By combining information extracted from both the HSCOI and the canopy height model, predominant stem height (r² = 0.91, RMSE = 0.77 m, n = 30), crown cover (r² = 0.78, RMSE = 9.25%, n = 30), and Foliage & Branch Projective Cover (FBPC; r² = 0.89, RMSE = 5.49%, n = 30) were estimated to levels sufficient for inventory of woodland and open forest structural types. When the approach was applied to forests in north east Victoria, stem density and crown cover were reliably estimated for forests with a structure similar to those observed in Queensland, but less so for forests of greater height and canopy closure.     

Assessing the utility of airborne hyperspectral and LiDAR data for species distribution mapping in the coastal Pacific Northwest, Canada

To effectively manage forested ecosystems an accurate characterization of species distribution is required. In this study we assess the utility of hyperspectral Airborne Imaging Spectrometer for Applications (AISA) imagery and small footprint discrete return Light Detection and Ranging (LiDAR) data for mapping 11 tree species in and around the Gulf Islands National Park Reserve, in coastal South-western Canada. Using hyperspectral imagery yielded producer's and user's accuracies for most species ranging from > 52-95.4 and > 63-87.8%, respectively. For species dominated by definable growth stages, pixel-level fusion of hyperspectral imagery with LiDAR-derived height and volumetric canopy profile data increased both producer's (+ 5.1-11.6%) and user's (+ 8.4-18.8%) accuracies. McNemar's tests confirmed that improvements in overall accuracies associated with the inclusion of LiDAR-derived structural information were statistically significant (p < 0.05). This methodology establishes a specific framework for mapping key species with greater detail and accuracy then is possible using conventional approaches (i.e., aerial photograph interpretation), or either technology on its own. Furthermore, in the study area, acquisition and processing costs were lower than a conventional aerial photograph interpretation campaign, making hyperspectral/LiDAR fusion a viable replacement technology.     

A pseudo-waveform technique to assess forest structure using discrete lidar data

The use of airborne laser scanning systems (lidar) to describe forest structure has increased dramatically since height profiling experiments nearly 30 years ago. The analyses in most studies employ a suite of frequency-based metrics calculated from the lidar height data, which are systematically eliminated from a full model using stepwise multiple linear regression. The resulting models often include highly correlated predictors with little physical justification for model formulation. We propose a method to aggregate discrete lidar height and intensity measurements into larger footprints to create “pseudo-waves”. Specifically, the returns are first sorted into height bins, sliced into narrow discrete elements, and finally smoothed using a spline function. The resulting “pseudo-waves” have many of the same characteristics of traditional waveform lidar data. We compared our method to a traditional frequency-based method to estimate tree height, canopy structure, stem density, and stand biomass in coniferous and deciduous stands in northern Wisconsin (USA). We found that the pseudo-wave approach had strong correlations for nearly all tree measurements including height (cross validated adjusted R² (R²_cv) = 0.82, RMSE_cv = 2.09 m), mean stem diameter (R²_cv = 0.64, RMSE_cv = 6.15 cm), total aboveground biomass (R²_cv = 0.74, RMSE_cv = 74.03 kg ha^− 1), and canopy coverage (R²_cv = 0.79, RMSE_cv = 5%). Moreover, the type of wave (derived from height and intensity or from height alone) had little effect on model formulation and fit. When wave-based and frequency-based models were compared, fit and mean square error were comparable, leading us to conclude that the pseudo-wave approach is a viable alternative because it has 1) an increased breadth of available metrics; 2) the potential to establish new meaningful metrics that capture unique patterns within the waves; 3) the ability to explain metric selection based on the physical structure of forests; and 4) lower correlation among independent variables.     

Forest structure and aboveground biomass in the southwestern United States from MODIS and MISR

Red band bidirectional reflectance factor data from the NASA MODerate resolution Imaging Spectroradiometer (MODIS) acquired over the southwestern United States were interpreted through a simple geometric-optical (GO) canopy reflectance model to provide maps of fractional crown cover (dimensionless), mean canopy height (m), and aboveground woody biomass (Mg ha^− 1) on a 250 m grid. Model adjustment was performed after dynamic injection of a background contribution predicted via the kernel weights of a bidirectional reflectance distribution function (BRDF) model. Accuracy was assessed with respect to similar maps obtained with data from the NASA Multiangle Imaging Spectroradiometer (MISR) and to contemporaneous US Forest Service (USFS) maps based partly on Forest Inventory and Analysis (FIA) data. MODIS and MISR retrievals of forest fractional cover and mean height both showed compatibility with the USFS maps, with MODIS mean absolute errors (MAE) of 0.09 and 8.4 m respectively, compared with MISR MAE of 0.10 and 2.2 m, respectively. The respective MAE for aboveground woody biomass was ~ 10 Mg ha^− 1, the same as that from MISR, although the MODIS retrievals showed a much weaker correlation, noting that these statistics do not represent evaluation with respect to ground survey data. Good height retrieval accuracies with respect to averages from high resolution discrete return lidar data and matches between mean crown aspect ratio and mean crown radius maps and known vegetation type distributions both support the contention that the GO model results are not spurious when adjusted against MISR bidirectional reflectance factor data. These results highlight an alternative to empirical methods for the exploitation of moderate resolution remote sensing data in the mapping of woody plant canopies and assessment of woody biomass loss and recovery from disturbance in the southwestern United States and in parts of the world where similar environmental conditions prevail.     

Effects of different sensors, flying altitudes, and pulse repetition frequencies on forest canopy metrics and biophysical stand properties derived from small-footprint airborne laser data

Canopy height distributions were created from small-footprint airborne laser scanner (ALS) data collected over 40 field sample plots with size 1000 m² located in mature conifer forest. ALS data were collected with two different instruments, i.e., the ALTM 1233 and ALTM 3100 laser scanners (Optech Inc.). The ALTM 1233 data were acquired at a flying altitude of 1200 m and a pulse repetition frequency (PRF) of 33 kHz. Three different acquisitions were carried out with ALTM 3100, i.e., (1) a flying altitude of 1100 m and a PRF of 50 kHz, (2) a flying altitude of 1100 m and a PRF of 100 kHz, and (3) a flying altitude of 2000 m and a PRF of 50 kHz. Height percentiles, mean and maximum height values, coefficients of variation of the heights, and canopy density at different height intervals above the ground were derived from the four different ALS datasets and for single + first and last echoes of the ALS data separately. The ALS-derived height- and density variables were assessed in pair-wise comparisons to evaluate the effects of (a) instrument, (b) flying altitude, and (c) PRF. A systematic shift in height values of up to 0.3 m between sensors when the first echoes were compared was demonstrated. Also the density-related variables differed significantly between the two instruments. Comparisons of flying altitudes and PRFs revealed upwards shifted canopy height distributions for the highest flying altitude (2000 m) and the lowest PRF (50 kHz). The distribution of echoes on different echo categories, i.e., single and multiple (first and last) echoes, differed significantly between acquisitions. The proportion of multiple echoes decreased with increasing flying altitude and PRF. Different echo categories have different properties since it is likely that single echoes tend to occur in the densest parts of the tree crowns, i.e., near the apex where the concentration of biological matter is highest and distance to the ground is largest. To assess the influence of instrument, flying altitude, and PRF on biophysical properties derived from ALS data, regression analysis was carried out to relate ALS-derived metrics to mean tree height (h_L) and timber volume (V). Cross validation revealed only minor differences in precision for the different ALS acquisitions, but systematic differences between acquisitions of up to 2.5% for h_L and 10.7% for V were found when comparing data from different acquisitions.     

Height growth reconstruction of a boreal forest canopy over a period of 58 years using a combination of photogrammetric and lidar models

Field data describing the height growth of trees or stands over several decades are very scarce. Consequently, our capacity of analyzing forest dynamics over large areas and long periods of time is somewhat limited. This study proposes a new method for retrospectively reconstructing plot-wise average dominant tree height based on a time series of high-resolution canopy height maps, termed canopy height models (CHMs). The absolute elevation of the canopy surface, or digital surface model (DSM), was first reconstructed by applying image-matching techniques to stereo-pairs of aerial photographs acquired in 1945, 1965, 1983, and 2003. The historical CHMs were then created by subtracting the bare earth elevation provided from a recent lidar survey from the DSMs. A method for estimating average dominant tree height from these historical CHMs was developed and calibrated for each photographic year. The accuracy of the resulting remote sensing height estimates was compared to age-height data reconstructed based on dendrometric measurements. The height bias of the remote sensing estimates relative to the verification data ranged from 0.52 m to 1.55 m (1.16 m on average). The corresponding root-mean-square errors varied between 1.49 m and 2.88 m (2.03 m average). Despite being slightly less accurate than historical field data, the quality of the remote sensing estimates is sufficient for many types of forest dynamics studies. The procedures for implementing this method, with the exception of the calibration phase, are entirely automated such that forest height growth curves can be reconstructed and mapped over large areas for which recent lidar data and historical photographs exist.     

Satellite lidar vs. small footprint airborne lidar: Comparing the accuracy of aboveground biomass estimates and forest structure metrics at footprint level

The use of lidar remote sensing for mapping the spatial distribution of canopy characteristics has the potential to allow an accurate and efficient estimation of tree dimensions and canopy structural properties from local to regional and continental scales. The overall goal of this paper was to compare biomass estimates and height metrics obtained by processing GLAS waveform data and spatially coincident discrete-return airborne lidar data over forest conditions in east Texas. Since biomass estimates are derived from waveform height metrics, we also compared ground elevation measurements and canopy parameters. More specific objectives were to compare the following parameters derived from GLAS and airborne lidar: (1) ground elevations; (2) maximum canopy height; (3) average canopy height; (4) percentiles of canopy height; and (5) above ground biomass. We used the elliptical shape of GLAS footprints to extract canopy height metrics and biomass estimates derived from airborne lidar. Results indicated a very strong correlation for terrain elevations between GLAS and airborne lidar, with an r value of 0.98 and a root mean square error of 0.78 m. GLAS height variables were able to explain 80% of the variance associated with the reference biomass derived from airborne lidar, with an RMSE of 37.7 Mg/ha. Most of the models comparing GLAS and airborne lidar height metrics had R-square values above 0.9.     

Physically based vertical vegetation structure retrieval from ICESat data: Validation using LVIS in White Mountain National Forest, New Hampshire, USA

Vegetation structure retrieval accuracies from spaceborne Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat) data are affected by surface topography, background noise and sensor saturation. This study uses a physical approach to remove surface topography effect from lidar returns to retrieve vegetation height from ICESat/GLAS data over slope terrains. Slope-corrected vegetation heights from ICESat/GLAS data were compared to airborne Laser Vegetation Imaging Sensor (LVIS) (20 m footprint size) and small-footprint lidar data collected in White Mountain National Forest, NH. Impact of slope on LVIS vegetation height estimates was assessed by comparing LVIS height before and after slope correction with small-footprint discrete-return lidar and field data.Slope-corrected GLAS vegetation heights match well with 98 percentile heights from small-footprint lidar (R² = 0.77, RMSE = 2.2 m) and top three LVIS mean (slope-corrected) heights (R² = 0.64, RMSE = 3.7 m). Impact of slope on LVIS heights is small, however, comparison of LVIS heights (without slope correction) with either small footprint lidar or field data indicates that our scheme improves the overall LVIS height accuracy by 0.4-0.7 m in this region. Vegetation height can be overestimated by 3 m over a 15° slope without slope correction. More importantly, both slope-corrected GLAS and LVIS height differences are independent of slope. Our results demonstrate the effectiveness of the physical approach to remove surface topography from large footprint lidar data to improve accuracy of maximum vegetation height estimates.GLAS waveforms were compared to aggregated LVIS waveforms in Bartlett Experimental Forest, NH, to evaluate the impact of background noise and sensor saturation on vegetation structure retrievals from ICESat/GLAS. We found that GLAS waveforms with sensor saturation and low background noise match well with aggregated LVIS waveforms, indicating these waveforms capture vertical vegetation structure well. However, waveforms with large noise often lead to mismatched waveforms with LVIS and underestimation of waveform extent and vegetation height. These results demonstrate the quality of ICESat/GLAS vegetation structure estimates.     

美国高分辨率土地覆盖信息提取技术研究进展

土地覆盖数据是进行全球变化研究的基础。美国地质调查局组织的土地覆盖特征研究项目成功的开展了两次工作, 建设了土地覆盖数据库, 得到了多方面的认可, 代表了美国相关领域的前沿技术。项目中涉及一系列技术, 包括土地覆盖分类分区、遥感影像选取、数据预处理过程、图像转换变化分析、土地覆盖分类技术、数据产品验证、土地覆盖数据库建设等。另外还有关键技术研究, 包括: 土地覆盖多边形分析、纹理特征分析、树冠密度信息提取、城市不透性表面估算等。就美国地质调查局项目中土地覆盖遥感影像数据处理方法、技术流程、数据库建设, 以及有关土地覆盖度等相关关键技术进行了介绍, 期望能够反映其进展情况, 对国内的相关工作起到启示作用。     

Updating the 2001 National Land Cover Database land cover classification to 2006 by using Landsat imagery change detection methods

The recent release of the U.S. Geological Survey (USGS) National Land Cover Database (NLCD) 2001, which represents the nation's land cover status based on a nominal date of 2001, is widely used as a baseline for national land cover conditions. To enable the updating of this land cover information in a consistent and continuous manner, a prototype method was developed to update land cover by an individual Landsat path and row. This method updates NLCD 2001 to a nominal date of 2006 by using both Landsat imagery and data from NLCD 2001 as the baseline. Pairs of Landsat scenes in the same season in 2001 and 2006 were acquired according to satellite paths and rows and normalized to allow calculation of change vectors between the two dates. Conservative thresholds based on Anderson Level I land cover classes were used to segregate the change vectors and determine areas of change and no-change. Once change areas had been identified, land cover classifications at the full NLCD resolution for 2006 areas of change were completed by sampling from NLCD 2001 in unchanged areas. Methods were developed and tested across five Landsat path/row study sites that contain several metropolitan areas including Seattle, Washington; San Diego, California; Sioux Falls, South Dakota; Jackson, Mississippi; and Manchester, New Hampshire. Results from the five study areas show that the vast majority of land cover change was captured and updated with overall land cover classification accuracies of 78.32%, 87.5%, 88.57%, 78.36%, and 83.33% for these areas. The method optimizes mapping efficiency and has the potential to provide users a flexible method to generate updated land cover at national and regional scales by using NLCD 2001 as the baseline.     

Forest canopy height from the Multiangle Imaging SpectroRadiometer (MISR) assessed with high resolution discrete return lidar

In this study retrievals of forest canopy height were obtained through adjustment of a simple geometric-optical (GO) model against red band surface bidirectional reflectance estimates from NASA's Multiangle Imaging SpectroRadiometer (MISR), mapped to a 250 m grid. The soil-understory background contribution was partly isolated prior to inversion using regression relationships with the isotropic, geometric, and volume scattering kernel weights of a Li-Ross kernel-driven bidirectional reflectance distribution function (BRDF) model. The height retrievals were assessed using discrete return lidar data acquired over sites in Colorado as part of the Cold Land Processes Experiment (CLPX) and used with fractional crown cover retrievals to obtain aboveground woody biomass estimates. For all model runs with reasonable backgrounds and initial b/r (vertical to horizontal crown radii) values < 2.0, root mean square error (RMSE) distributions were centered between 2.5 and 3.7 m while R² distributions were centered between 0.4 and 0.7. The MISR/GO aboveground biomass estimates predicted via regression on fractional cover and mean canopy height for the CLPX sites showed good agreement with U.S. Forest Service Interior West map data (adjusted R² = 0.84). The implication is that multiangle sensors such as MISR can provide spatially contiguous retrievals of forest canopy height, cover, and aboveground woody biomass that are potentially useful in mapping distributions of aboveground carbon stocks, tracking disturbance, and in initializing, constraining, and validating ecosystem models. This is important because the MISR record is spatially comprehensive and extends back to the year 2000 and the launch of the NASA Earth Observing System (EOS) Terra satellite; it might thus provide a ~ 10-year baseline record that would enhance exploitation of data from the NASA Deformation, Ecosystem Structure and Dynamics of Ice (DESDynI) mission, as well as furthering realization of synergies with active instruments.     

Retrieval of canopy height using moderate-resolution imaging spectroradiometer (MODIS) data

In this study we use the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF) product to develop multivariate linear regression models that estimate canopy heights over study sites at Howland Forest, Maine, Harvard Forest, Massachusetts and La Selva Forest, Costa Rica using (1) directional escape probabilities that are spectrally independent and (2) the directional spectral reflectances used to derive the directional escape probabilities. These measures of canopy architecture are compared with canopy height information retrieved from the airborne Laser Vegetation Imaging Sensor (LVIS). Both the escape probability and the directional reflectance approaches achieve good results, with correlation coefficients in the range 0.54-0.82, although escape probability results are usually slightly better. This suggests that MODIS 500 m BRDF data can be used to extrapolate canopy heights observed by widely-spaced satellite LIDAR swaths to larger areas, thus providing wide-area coverage of canopy height.     

Vegetation height estimation from Shuttle Radar Topography Mission and National Elevation Datasets

A study was conducted to determine the feasibility of obtaining estimates of vegetation canopy height from digital elevation data collected during the 2000 Shuttle Radar Topography Mission (SRTM). The SRTM sensor mapped 80% of the Earth's land mass with a C-band Interferometric Synthetic Aperture Radar (InSAR) instrument, producing the most complete digital surface map of Earth. Due to the relatively short wavelength (5.6 cm) of the SRTM instrument, the majority of incoming electromagnetic energy is reflected by scatterers located within the vegetation canopy at heights well above the “bald-Earth” surface. Interferometric SAR theory provides a basis for properly identifying and accounting for the dependence of this scattering phase center height on both instrument and target characteristics, including relative and absolute vertical error and vegetation structural attributes.An investigation to quantify the magnitude of the vertical error component was conducted using SRTM data from two vegetation-free areas in Iowa and North Dakota, revealing absolute errors of −4.0 and −1.1 m, respectively. It was also shown that the relative vertical error due to phase noise can be reduced significantly through sample averaging. The relative error range for the Iowa site was reduced from 13 to 4 m and for the North Dakota site from 7 to 3 m after averaging of 50 samples. Following error reduction, it was demonstrated that the SRTM elevation data can be successfully correlated via linear regression models with ground-measured canopy heights acquired during the general mission timeframe from test sites located in Georgia and California. Prior to outlier removal and phase noise reduction, initial adjusted r² values for the Georgia and California sites were 0.15 and 0.20, respectively. Following outlier analysis and averaging of at least 20 SRTM pixels per observation, adjusted r² values for the Georgia and California sites improved to 0.79 (rmse=1.1 m) and 0.75 (rmse=4.5 m), respectively. An independent validation of a novel bin-based modeling strategy designed for reducing phase noise in sample plot data confirmed both the robustness of the California model (adjusted r²=0.74) as well as the capacity of the binning strategy to produce stable models suitable for inversion (validated rmse=4.1 m). The results suggest that a minimum mapping unit of approximately 1.8 ha is appropriate for SRTM-based vegetation canopy height mapping.     

The uncertainty in conifer plantation growth prediction from multi-temporal lidar datasets

An evaluation of the use of airborne lidar for multi-temporal forest height growth assessment in a temperate mature red pine (Pinus resinosa Ait.) plantation over a five-year period is presented. The objective was to evaluate the level of uncertainty in lidar-based growth estimates through time so that the optimal repeat interval necessary for statistically meaningful growth measurements could be evaluated. Four airborne lidar datasets displaying similar survey configuration parameters were collected between 2000 and 2005. Coincident with the 2002 and 2005 acquisitions, field mensuration for 126 trees within 19 plots was carried out. Field measurements of stem height were compared to both coincident plot-level laser pulse return (LPR) height percentile metrics and stand level raster canopy height models (CHM).The average plot-level field heights were found to be 23.8 m (standard deviation (σ) = 0.4 m) for 2002 and 25.0 m (σ = 0.6 m) for 2005, with an approximate annual growth rate of 0.4 m/yr (σ = 0.5 m). The standard deviation uncertainty for field height growth estimates over the three year period was 41% at the plot-level (n = 19) and 92% at the individual tree level (n = 126). Of the lidar height percentile metrics tested, the 90th (L90), 95th (L95) and maximum (Lmax) LPR distribution heights demonstrated the highest overall correlations with field-measured tree height. While all lidar-based methods, including raster CHM comparison, tended to underestimate the field estimate of growth, Lmax provided the most robust overall direct estimate (0.32 m/yr, σ = 0.37 m). A single factor analysis of variance demonstrated that there was no statistically significant difference between all plot-level field and Lmax mean growth rate estimates (P = 0.38) and, further, that there was no difference in Lmax growth rate estimates across the examined time intervals (P = 0.59). A power function relationship between time interval and the standard deviation of error in growth estimate demonstrated that over a one-year period, the growth uncertainty was in the range of 0.3 m (∼ 100% of total growth) reducing to less than 0.1 m (∼ 6% of total growth) after 5 years. Assuming a 10% uncertainty is acceptable for operational or research-based conifer plantation growth estimates, this can be achieved at a three-year time interval.     

Integrating profiling LIDAR with Landsat data for regional boreal forest canopy attribute estimation and change characterization

Forest dynamics are characterized by both continuous (i.e., growth) and discontinuous (i.e., disturbance) changes. Change detection techniques that use optical remotely sensed data to capture disturbance related changes are established and commonly applied; however, approaches for the capture of continuous forest changes are less mature. Optical remotely sensed imagery is well suited for capturing horizontally distributed conditions, structures, and changes, while Light Detection And Ranging (LIDAR) data are more appropriate for capturing vertically distributed elements of forest structure and change. The integration of optical remotely sensed imagery and LIDAR data provides improved opportunities to fully characterize forest canopy attributes and dynamics.The study described in this paper captures forest conditions along a corridor approximately 600 km long through the boreal forest of Canada. Two coincident LIDAR transects, representing 1997 and 2002 forest conditions respectively, are compared using image segments generated from Landsat ETM+ imagery. The image segments are used to provide a spatial framework within which the attributes and temporal dynamics of the forest canopy are estimated and compared. Segmented and classified Landsat imagery provides a context for the comparison of sufficiently spatially related LIDAR profiles and for the provision of categories to aid in the application of empirical models requiring knowledge of land cover.Global and local approaches were employed for characterizing changes in forest attributes over time. The global approach, emphasized the overall trend in forest change along the length of the entire transect, and indicated that key canopy attributes were stable, and transect characteristics, including forest canopy height, did not change significantly over the five-year period of this study (two sample t-test, p = 0.08). The local approach analyzed segment-based changes in canopy attributes, providing spatially explicit indications of forest growth and depletion. The local approach identified that 84% of the Landsat segments intercepted by both LIDAR transects either have no change, or have a small average increase in canopy height (0.7 m), while the other 16% of segments have an average decrease in canopy height of 1.6 m. As expected, the difference in the magnitude of the changes was markedly greater for depletions than it was for growth, but was less spatially extensive. Growth tends to occur incrementally over broad areas; whereas, depletions are dramatic and spatially constrained. The approach presented holds potential for investigating the impacts of climate change across a latitudinal gradient of boreal forest.     

Three-dimensional canopy fuel loading predicted using upward and downward sensing LiDAR systems

We calibrated upward sensing profiling and downward sensing scanning LiDAR systems to estimates of canopy fuel loading developed from field plots and allometric equations, and then used the LiDAR datasets to predict canopy bulk density (CBD) and crown fuel weight (CFW) in wildfire prone stands in the New Jersey Pinelands. LiDAR-derived height profiles were also generated in 1-m layers, and regressed on CBD estimates calculated for 1-m layers from field plots to predict three-dimensional canopy fuel loading. We then produced maps of canopy fuel metrics for three 9 km² forested areas in the Pinelands.Correlations for standard LiDAR-derived parameters between the two LiDAR systems were all highly significant, with correlation coefficients ranging between 0.82 and 0.98. Stepwise linear regression models developed from the profiling LiDAR data predicted maximum CBD and CFW (r² = 0.94 and 0.92) better than those developed from the scanning LiDAR data (r² = 0.83 and 0.71, respectively). A single regression for the prediction of CBD at all canopy layers had r² values of 0.93 and 0.82 for the profiling and scanning datasets, respectively. Individual bin regressions for predicting CBD at each canopy height layer were also highly significant at most canopy heights, with r² values for each layer ranging between 0.36 and 0.89, and 0.44 and 0.99 for the profiling and scanning datasets, respectively. Relationships were poorest mid-canopy, where highest average values and highest variability in fuel loading occurred. Fit of data to Gaussian distributions of canopy height profiles facilitated a simpler expression of these parameters for analysis and mapping purposes, with overall r² values of 0.86 and 0.92 for the profiling and scanning LiDAR datasets, respectively. Our research demonstrates that LiDAR data can be used to generate accurate, three-dimensional representations of canopy structure and fuel loading at high spatial resolution by linking 1-m return height profiles to biometric estimates from field plots.