Estimating zero-plane displacement height and aerodynamic roughness length using synthesis of LiDAR and SPOT-5 data

In this study, a combination of low and high density airborne LiDAR and satellite SPOT-5 HRG data were used in conjunction with ground measurements of forest structure to parameterize four models for zero-plane displacement height d(m) and aerodynamic roughness length z_0m(m), over cool-temperate forests in Heihe River basin, an arid region of Northwest China. For the whole study area, forest structural parameters including tree height (Ht) (m), first branch height (FBH) (m), crown width (CW) (m) and stand density (SD)(trees ha^− 1) were derived by stepwise multiple linear regressions of ground-based forest measurements and height quantiles and fractional canopy cover (f_c) derived from the low density LiDAR data. The high density LiDAR data, which covered a much smaller area than the low density LiDAR data, were used to relate SPOT-5's reflectance to the effective plant area index (PAIe) of the forest. This was done by linear spectrum decomposition and Li-Strahler geometric-optical models. The result of the SPOT-5 spectrum decomposition was applied to the whole area to calculate PAIe (and leaf area index LAI). Then, four roughness models were applied to the study area with these vegetation data derived from the LiDAR and SPOT-5 as input. For validation, measurements at an eddy covariance site in the study area were used. Finally, the four models were compared by plotting histograms of the accumulative distribution of modeled d and z_0m in the study area. The results showed that the model using by frontal area index (FAI) produced best d estimate, and the model using both LAI and FAI generated the best z_0m. Furthermore, all models performed much better when the representative tree height was Lorey's mean height instead of using an arithmetic mean.     

Application of high spatial resolution satellite imagery for riparian and forest ecosystem classification

Terrestrial Ecosystem Mapping provides critical information to land and resource managers by incorporating information on climate, physiography, surficial material, soil, and vegetation structure. The main objective of this research was to determine the capacity of high spatial resolution satellite image data to discriminate vegetation structural stages in riparian and adjacent forested ecosystems as defined using the British Columbia Terrestrial Ecosystem Mapping (TEM) scheme. A high spatial resolution QuickBird image, captured in June 2005, and coincident field data covering the riparian area of Lost Shoe Creek and adjacent forests on Vancouver Island, British Columbia, was used in this analysis. Semi-variograms were calculated to assess the separability of vegetation structural stages and assess which spatial scales were most appropriate for calculation of grey-level co-occurrence texture measures to maximize structural class separation. The degree of spatial autocorrelation showed that most vegetation structural types in the TEM scheme could be differentiated and that window sizes of 3 × 3 pixels and 11 × 11 pixels were most appropriate for image texture calculations. Using these window sizes, the texture analysis showed that co-occurrence contrast, dissimilarity, and homogeneity texture measures, based on the bands in the visible part of the spectrum, provided the most significant statistical differentiation between vegetation structural classes. Subsequently, an object-oriented classification algorithm was applied to spectral and textural transformations of the QuickBird image data to map the vegetation structural classes. Using both spectral and textural image bands yielded the highest classification accuracy (overall accuracy = 78.95%). The inclusion of image texture increased the classification accuracies of vegetation structure by 2-19%. The results show that information on vegetation structure can be mapped effectively from high spatial resolution satellite image data, providing an additional tool to ongoing aerial photograph interpretation.     

Quantifying indicators of riparian condition in Australian tropical savannas: integrating high spatial resolution imagery and field survey data

The objectives of this research were: (1) to quantify indicators of riparian condition; and (2) to assess these indicators for detecting change in riparian condition. Two multi‐spectral QuickBird images were acquired in 2004 and 2005 for a section of the Daly River in north Australia. These data were collected coincidently with vegetation and geomorphic field data. Indicators of riparian condition, including percentage canopy cover, organic litter, canopy continuity, bank stability, flood damage, riparian zone width and vegetation overhang, were then mapped. Field measurements and vegetation indices were empirically related using regression analysis to develop algorithms for mapping organic litter and canopy cover (R ² = 0.59–0.78). Using a standard nearest‐neighbour algorithm, object‐oriented supervised image classification provided thematic information (overall accuracies 81–90%) for mapping riparian zone width and vegetation overhang. Bank stability and flood damage were mapped empirically from a combination of canopy cover information and the image classification products (R ² = 0.70–0.81). Multi‐temporal image analysis of riparian condition indicators (RCIs) demonstrated the advantages of using continuous and discrete data values as opposed to categorical data. This research demonstrates how remote sensing can be used for mapping and monitoring riparian zones in remote tropical savannas and other riparian environments at scales from 1 km to 100s km of stream length.     

Comparative evaluation of airborne LiDAR and ship-based multibeam SoNAR bathymetry and intensity for mapping coral reef ecosystems

Large areas of the world's coastal marine environments remain poorly characterized because they have not been mapped with sufficient accuracy and at spatial resolutions high enough to support a wide range of societal needs. Expediting the rate of seafloor mapping requires the collection of multi-use datasets that concurrently address hydrographic charting needs and support decision-making in ecosystem-based management. While active optical and acoustic sensors have previously been compared for the purpose of hydrographic charting, few studies have evaluated the performance and cost effectiveness of these systems for providing benthic habitat maps. Bathymetric and intensity data were collected in shallow water (< 50 m depth) coral reef ecosystems using two conventional remote sensing technologies: (1) airborne Light Detection and Ranging (LiDAR), and (2) ship-based multibeam (MBES) Sound Navigation and Ranging (SoNAR). A comparative assessment using a suite of twelve metrics demonstrated that LiDAR and MBES were equally capable of discriminating seafloor topography (r = > 0.9), although LiDAR depths were found to be consistently shallower than MBES depths. The intensity datasets were not significantly correlated at a broad 4 × 5 km spatial scale (r = − 0.11), but were moderately correlated in flat areas at a fine 4 × 500 m spatial scale (r = 0.51), indicating that the LiDAR intensity algorithm needs to be improved before LiDAR intensity surfaces can be used for habitat mapping. LiDAR cost 6.6% less than MBES and required 40 fewer hours to map the same study area. MBES provided more detail about the seafloor by fully ensonifying high-relief features, by differentiating between fine and coarse sediments and by collecting data with higher spatial resolutions. Surface fractal dimensions and fast Fourier transformations emerged as useful methods for detecting artifacts in the datasets. Overall, LiDAR provided a more cost effective alternative to MBES for mapping and monitoring shallow water coral reef ecosystems (< 50 m depth), although the unique advantages of MBES may make it a more appropriate choice for answering certain ecological or geological questions requiring very high resolution data.     

Assessment of QuickBird high spatial resolution imagery to detect red attack damage due to mountain pine beetle infestation

High spatial resolution remotely sensed data has the potential to complement existing forest health programs for both strategic planning over large areas, as well as for detailed and precise identification of tree crowns subject to stress and infestation. The area impacted by the current mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in British Columbia, Canada, has increased 40-fold over the previous 5 years, with approximately 8.5 million ha of forest infested in 2005. As a result of the spatial extent and intensity of the outbreak, new technologies are being assessed to help detect, map, and monitor the damage caused by the beetle, and to inform mitigation of future beetle outbreaks. In this paper, we evaluate the capacity of high spatial resolution QuickBird multi-spectral imagery to detect mountain pine beetle red attack damage. ANOVA testing of individual spectral bands, as well as the Normalized Difference Vegetation Index (NDVI) and a ratio of red to green reflectance (Red-Green Index or RGI), indicated that the RGI was the most successful (p < 0.001) at separating non-attack crowns from red attack crowns. Based on this result, the RGI was subsequently used to develop a binary classification of red attack and non-attack pixels. The total number of QuickBird pixels classified as having red attack damage within a 50 m buffer of a known forest health survey point were compared to the number of red attack trees recorded at the time of the forest health survey. The relationship between the number of red attack pixels and observed red attack crowns was assessed using independent validation data and was found to be significant (r² = 0.48, p < 0.001, standard error = 2.8 crowns). A comparison of the number of QuickBird pixels classified as red attack, and a broader scale index of mountain pine beetle red attack damage (Enhanced Wetness Difference Index, calculated from a time series of Landsat imagery), was significant (r² = 0.61, p < 0.001, standard error = 1.3 crowns). These results suggest that high spatial resolution imagery, in particular QuickBird satellite imagery, has a valuable role to play in identifying tree crowns with red attack damage. This information could subsequently be used to augment existing detailed forest health surveys, calibrate synoptic estimates of red attack damage generated from overview surveys and/or coarse scale remotely sensed data, and facilitate the generation of value-added information products, such as estimates of timber volume impacts at the forest stand level.     

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.     

Development of a simulation model to predict LiDAR interception in forested environments

Airborne scanning LiDAR systems are used to predict a range of forest attributes. However, the accuracy with which this can be achieved is highly dependent on the sensor configuration and the structural characteristics of the forest examined. As a result, there is a need to understand laser light interactions with forest canopies so that LiDAR sensor configurations can be optimised to assess particular forest types. Such optimisation will not only ensure the targeted forest attributes can be accurately and consistently quantified, but may also minimise the cost of data acquisition and indicate when a survey configuration will not deliver information needs.In this paper, we detail the development and application of a model to simulate laser interactions within forested environments. The developed model, known as the LiDAR Interception and Tree Environment (LITE) model, utilises a range of structural configurations to simulate trees with variable heights, crown dimensions and foliage clumping. We developed and validated the LITE model using field data obtained from three forested sites covering a range of structural classes. Model simulations were then compared to coincident airborne LiDAR data collected over the same sites. Results indicate that the LITE model can be used to produce comparable estimates of maximum height of trees within plots (differences < 2.42 m), mean heights of first return data (differences < 2.27 m), and canopy height percentiles (r² = 0.94, p < 0.001) when compared to airborne LiDAR data. In addition, the distribution of airborne LiDAR hits through the canopy profile was closely matched by model predictions across the range of sites. Importantly, this demonstrates that the structural differences between forest stands can be characterised by LITE. Models that are capable of interpreting the response of small-footprint LiDAR waveforms can facilitate algorithm development, the generation of corrections for actual LiDAR data, and the optimisation of sensor configurations for differing forest types, benefiting a range of experimental and commercial LiDAR applications. As a result, we also performed a scenario analysis to demonstrate how differences in forest structure, terrain, and sensor configuration can influence the interception of LiDAR beams.     

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.     

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.     

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.     

Assessing the extent of agriculture/pasture and secondary succession forest in the Brazilian Legal Amazon using SPOT VEGETATION data

There has been growing concern about land use/land cover change in tropical regions, as there is evidence of its influence on the observed increase in atmospheric carbon dioxide concentration and consequent climatic changes. Mapping of deforestation by the Brazil's National Space Research Institute (INPE) in areas of primary tropical forest using satellite data indicates a value of 587,727 km² up to the year 2000. Although most of the efforts have been concentrated in mapping primary tropical forest deforestation, there is also evidence of large-scale deforestation in the cerrado savanna, the second most important biome in the region.The main purpose of this work was to assess the extent of agriculture/pasture and secondary succession forest in the Brazilian Legal Amazon (BLA) in 2000, using a set of multitemporal images from the 1-km SPOT-4 VEGETATION (VGT) sensor. Additionally, we discriminated primary tropical forest, cerrado savanna, and natural/artificial waterbodies. Four classification algorithms were tested: quadratic discriminant analysis (QDA), simple classification trees (SCT), probability-bagging classification trees (PBCT), and k-nearest neighbors (K-NN). The agriculture/pasture class is a surrogate for those areas cleared of its original vegetation cover in the past, acting as a source of carbon. On the contrary, the secondary succession forest class behaves as a sink of carbon.We used a time series of 12 monthly composite images of the year 2000, derived from the SPOT-4 VGT sensor. A set of 19 Landsat scenes was used to select training and testing data. A 10-fold cross validation procedure rated PBCT as the best classification algorithm, with an overall sample accuracy of 0.92. High omission and commission errors occurred in the secondary succession forest class, due to confusion with agriculture/pasture and primary tropical forest classes. However, the PBCT algorithm generated the lower misclassification error in this class. Besides, this algorithm yields information about class membership probability, with ∼80% of the pixels with class membership probability greater or equal than 0.8. The estimated total area of agriculture/pasture and secondary succession forest in 2000 in the BLA was 966 × 10³ and 140 × 10³ km², respectively. Comparison with an existing land cover map indicates that agriculture/pasture occurred primarily in areas previously occupied by primary tropical forest (46%) and cerrado savanna (33%), and also in transition forest (19%), and other vegetation types (2%). This further confirms the existing evidence of extensive cerrado savanna conversion.This study also concludes that SPOT-4 VGT data are adequate for discriminating several major land cover types in tropical regions. Agriculture/pasture was mapped with errors of about 5%. Very high classification errors were associated with secondary succession forest, suggesting that a different methodology/sensor has to be used to address this difficult land cover class (namely with the inclusion of ancillary data). For the other classes, we consider that accurate maps can be derived from SPOT-4 VGT data with errors lower than 20% for the cerrado savanna, and errors lower than 10% for the other land cover classes. These estimates may be useful to evaluate impacts of land use/land cover change on the carbon and water cycles, biotic diversity, and soil degradation.     

Mapping seagrass species, cover and biomass in shallow waters: An assessment of satellite multi-spectral and airborne hyper-spectral imaging systems in Moreton Bay (Australia)

Understanding, monitoring and modelling attributes of seagrass biodiversity, such as species composition, richness, abundance, spatial patterns, and disturbance dynamics, requires spatial information. This work assessed the accuracy of commonly available airborne hyper-spectral and satellite multi-spectral image data sets for mapping seagrass species composition, horizontal horizontal-projected foliage cover and above-ground dry-weight biomass. The work was carried out on the Eastern Banks in Moreton Bay, Australia, an area of shallow and clear coastal waters, containing a range of seagrass species, cover and biomass levels. Two types of satellite image data were used: Quickbird-2 multi-spectral and Landsat-5 Thematic Mapper multi-spectral. Airborne hyper-spectral image data were acquired from a CASI-2 sensor using a pixel size of 4.0 m. The mapping was constrained to depths shallower than 3.0 m, based on past modelling of the separability of seagrass reflectance signatures at increasing water depths. Our results demonstrated that mapping of seagrass cover, species and biomass to high accuracy levels (> 80%) was not possible across all image types. For each parameter mapped, airborne hyper-spectral data produced the highest overall accuracies (46%), followed by Quickbird-2 and then Landsat-5 Thematic Mapper. The low accuracy levels were attributed to the mapping methods and difficulties in matching locations on image and field data sets. Accurate mapping of seagrass cover, species composition and biomass, using simple approaches, requires further work using high-spatial resolution (< 5 m) and/or hyper-spectral image data. Further work is required to determine if and how the seagrass maps produced in this work are suitable for measuring attributes of seagrass biodiversity, and using these data for modelling floral and fauna biodiversity properties of seagrass environments, and for scaling-up seagrass ecosystem models.     

Characterizing forest ecological structure using pulse types and heights of airborne laser scanning

Characterizing forest structure is an important part of any comprehensive biodiversity assessment. However, current methods for measuring structural complexity require a laborious process that involves many logistically expensive point based measurements. An automated or semi-automated method would be ideal. In this study, the utility of airborne laser scanning (LiDAR; Light Detection and Ranging) for characterizing the ecological structure of a forest landscape is examined. The innovation of this paper is to use different laser pulse return properties from a full waveform LiDAR to characterize forest ecological structure. First the LiDAR dataset is stratified into four vertical layers: ground, low vegetation (0-1 m from the ground), medium vegetation (1-5 m from the ground) and high vegetation (> 5 m). Subsequently the “Type” of LiDAR return is analysed: Type 1 (singular returns); Type 2 (first of many returns); Type 3 (intermediate returns); and Type 4 (last of many returns). A forest characterization scheme derived from LiDAR point clouds is proposed. A validation of the scheme is then presented using a network of field sites that recorded commonly used metrics of biodiversity. The proposed forest characterization categories allow for quantification of gaps (above bare ground, low vegetation and medium vegetation), canopy cover and its vertical density as well as the presence of various canopy strata (low, medium and high). Regression analysis showed that LiDAR derived variables were good predictors of field recorded variables (R² = 0.82, P < 0.05 between LiDAR derived presence of low vegetation and field derived LAI for low vegetation). The proposed scheme clearly shows the potential of full waveform LiDAR to provide information on the complexity of habitat structure.     

Segment-constrained regression tree estimation of forest stand height from very high spatial resolution panchromatic imagery over a boreal environment

Mean stand height is an important parameter for forest volume and biomass estimation in support of monitoring and management activities. Information on mean stand height is typically obtained through the manual interpretation of aerial photography, often supplemented by the collection of field calibration data. In remote areas where forest management practices may not be spatially exhaustive or where it is difficult to acquire aerial photography, alternate approaches for estimating stand height are required. One approach is to use very high spatial resolution (VHSR) satellite imagery (pixels sided less than 1 m) as a surrogate for air photos. In this research we demonstrate an approach for modelling mean stand height at four sites in the Yukon Territory, Canada, from QuickBird panchromatic imagery. An object-based approach was used to generate homogenous segments from the imagery (analogous to manually delineated forest stands) and an algorithm was used to automatically delineate individual tree crowns within the segments. A regression tree was used to predict mean stand height from stand-level metrics generated from the image grey-levels and within-stand objects relating individual tree crown characteristics. Heights were manually interpreted from the QuickBird imagery and divided into separate sets of calibration and validation data. The effects of calibration data set size and the input metrics used on the regression tree results were also assessed. The approach resulted in a model with a significant R² of 0.53 and an RMSE of 2.84 m. In addition, 84.6% of the stand height estimates were within the acceptable error for photo interpreted heights, as specified by the forest inventory standards of British Columbia. Furthermore, residual errors from the model were smallest for the stands that had larger mean heights (i.e., > 20 m), which aids in reducing error in subsequent estimates of biomass or volume (since stands with larger trees contribute more to overall estimates of volume or biomass). Estimated and manually interpreted heights were reclassified into 5-metre height classes (a schema frequently used for forest analysis and modelling applications) and compared; classes corresponded in 54% of stands assessed, and all stands had an estimated height class that was within ± 1 class of their actual class. This study demonstrates the capacity of VHSR panchromatic imagery (in this case QuickBird) for generating useful estimates of mean stand heights in unmonitored, remote, or inaccessible forest areas.     

Wetland monitoring using classification trees and SPOT-5 seasonal time series

Multiseason reflectance data from radiometrically and geometrically corrected multispectral SPOT-5 images of 10-m resolution were combined with thorough field campaigns and land cover digitizing using a binary classification tree algorithm to estimate the area of marshes covered with common reeds (Phragmites australis) and submerged macrophytes (Potamogeton pectinatus, P. pusillus, Myriophyllum spicatum, Ruppia maritima, Chara sp.) over an area of 145,000 ha. Accuracy of these models was estimated by cross-validation and by calculating the percentage of correctly classified pixels on the resulting maps. Robustness of this approach was assessed by applying these models to an independent set of images using independent field data for validation. Biophysical parameters of both habitat types were used to interpret the misclassifications. The resulting trees provided a cross-validation accuracy of 98.7% for common reed and 97.4% for submerged macrophytes. Variables discriminating reed marshes from other land covers were the difference in the near-infrared band between March and June, the Optimized Soil Adjusted Vegetation Index of December, and the Normalized Difference Water Index (NDWI) of September. Submerged macrophyte beds were discriminated with the shortwave-infrared band of December, the NDWI of September, the red band of September and the Simple Ratio index of March. Mapping validations provided accuracies of 98.6% (2005) and 98.1% (2006) for common reed, and 86.7% (2005) and 85.9% (2006) for submerged macrophytes. The combination of multispectral and multiseasonal satellite data thus discriminated these wetland vegetation types efficiently. Misclassifications were partly explained by digitizing inaccuracies, and were not related to biophysical parameters for reedbeds. The classification accuracy of submerged macrophytes was influenced by the proportion of plants showing on the water surface, percent cover of submerged species, water turbidity, and salinity. Classification trees applied to time series of SPOT-5 images appear as a powerful and reliable tool for monitoring wetland vegetation experiencing different hydrological regimes even with a small training sample (N = 25) when initially combined with thorough field measurements.     

Retrieval of forest canopy attributes based on a geometric-optical model using airborne LiDAR and optical remote-sensing data

In the retrieval of forest canopy attributes using a geometric-optical model, the spectral scene reflectance of each component should be known as prior knowledge. Generally, these reflectances were acquired by a foregone survey using an analytical spectral device. This article purposed to retrieve the forest structure parameters using light detection and ranging (LiDAR) data, and used a linear spectrum decomposition model to determine the reflectances of the spectral scene components, which are regarded as prior knowledge in the retrieval of forest canopy cover and effective plant area index (PAI_e) using a simplified Li–Strahler geometric-optical model based on a Satellites Pour l'Observation de la Terre 5 (SPOT-5) high-resolution geometry (HRG) image. The airborne LiDAR data are first used to retrieve the forest structure parameters and then the proportion of the SPOT pixel not covered by crown or shadow K_g of each pixel in the sample was calculated, which was used to extract the reflectances of the spectral scene components by a linear spectrum decomposition model. Finally, the forest canopy cover and PAI_e are retrieved by the geometric-optical model. As the acquired time of SPOT-5 image and measured data has a discrepancy of about 2 months, the retrieved result of forest canopy cover needs a further validation. The relatively high value of R ² between the retrieval result of PAI_e and the measurements indicates the efficiency of our methods.     

Optimising the use of hyperspectral and LiDAR data for mapping reedbed habitats

Reedbeds are important habitats for supporting biodiversity and delivering a range of ecosystem services, yet reedbeds in the UK are under threat from intensified agriculture, changing land use and pollution. To develop appropriate conservation strategies, information on the distribution of reedbeds is required. Field surveys of these wetland environments are difficult, time consuming and expensive to execute for large areas. Remote sensing has the potential to replace or complement such field surveys, yet the specific application to reedbed habitats has not been fully investigated. In the present study, airborne hyperspectral and LiDAR imagery were acquired for two sites in Cumbria, UK. The research aimed to determine the most effective means of analysing hyperspectral data covering the visible, near infrared (NIR) and shortwave infrared (SWIR) regions for mapping reedbeds and to investigate the effects of incorporating image textural information and LiDAR-derived measures of canopy structure on the accuracy of reedbed delineation. Due to the high dimensionality of the hyperspectral data, three image compression algorithms were evaluated: principal component analysis (PCA), spectrally segmented PCA (SSPCA) and minimum noise fraction (MNF). The LiDAR-derived measures tested were the canopy height model (CHM), digital surface model (DSM) and the DSM-derived slope map. The SSPCA-compressed data produced the highest reedbed accuracy and processing efficiency. The optimal SSPCA dataset incorporated 12 PCs comprised of the first 3 PCs derived from each of the spectral segments: visible (392-700 nm), NIR (701-972 nm), SWIR-1 (973-1366 nm) and SWIR-2 (1530-2240 nm). Incorporating image textural measures produced a significant improvement in the classification accuracy when using MNF-compressed data, but had no impact when using the SSPCA-compressed imagery. A significant improvement (+ 11%) in the accuracy of reedbed delineation was achieved when a mask generated by applying a 3 m threshold to the LiDAR-derived CHM was used to filter the reedbed map derived from the optimal SSPCA dataset. This paper demonstrates the value in combining appropriately compressed hyperspectral imagery with LiDAR data for the effective mapping of reedbed habitats.     

Reliability of LiDAR derived predictors of forest inventory attributes: A case study with Norway spruce

To increase the application domain (re-use) of LiDAR-based models the random replication effects in the predictor(s) must be considered. We quantify these effects in a linear predictor (X) of four forest inventory attributes (Lorey's height HT, basal area BA, volume VOL, and stem density TPH) with LiDAR data acquired over 40 spruce-dominated large plots in southeastern Norway. A grid-based random thinning of the raw multi-echo LiDAR data, to five target densities between 0.25 m^− 2 and 2.0 m^− 2, generated 100 replications with each density. A DTM was estimated for each replicate and target pulse density. The four linear predictors were constructed from two indicators of canopy density and a posited average effect of a power-transform of echoes classified as canopy returns. Replication variance varied significantly among plots but the reliability ratio of X was high (≥ 0.92) for HT, BA and VOL but lower for TPH, especially at low pulse densities. Reliability ratios increased with pulse density. Replication variance attenuated the linear regression coefficients by about 10% and inflated the residual variance by 3-6%. A proposed calibration was effective in reducing the impact of replication effects. A proposed bootstrap procedure can be used in practice to obtain good approximations of the replication variance. With echo-densities of approximately 1 m^− 2 or higher the replication effects do not warrant the effort of a calibration.     

Forest carbon densities and uncertainties from Lidar, QuickBird, and field measurements in California

Greenhouse gas inventories and emissions reduction programs require robust methods to quantify carbon sequestration in forests. We compare forest carbon estimates from Light Detection and Ranging (Lidar) data and QuickBird high-resolution satellite images, calibrated and validated by field measurements of individual trees. We conducted the tests at two sites in California: (1) 59 km² of secondary and old-growth coast redwood (Sequoia sempervirens) forest (Garcia-Mailliard area) and (2) 58 km² of old-growth Sierra Nevada forest (North Yuba area). Regression of aboveground live tree carbon density, calculated from field measurements, against Lidar height metrics and against QuickBird-derived tree crown diameter generated equations of carbon density as a function of the remote sensing parameters. Employing Monte Carlo methods, we quantified uncertainties of forest carbon estimates from uncertainties in field measurements, remote sensing accuracy, biomass regression equations, and spatial autocorrelation. Validation of QuickBird crown diameters against field measurements of the same trees showed significant correlation (r = 0.82, P < 0.05). Comparison of stand-level Lidar height metrics with field-derived Lorey's mean height showed significant correlation (Garcia-Mailliard r = 0.94, P < 0.0001; North Yuba R = 0.89, P < 0.0001). Field measurements of five aboveground carbon pools (live trees, dead trees, shrubs, coarse woody debris, and litter) yielded aboveground carbon densities (mean ± standard error without Monte Carlo) as high as 320 ± 35 Mg ha^− 1 (old-growth coast redwood) and 510 ± 120 Mg ha^− 1 (red fir [Abies magnifica] forest), as great or greater than tropical rainforest. Lidar and QuickBird detected aboveground carbon in live trees, 70-97% of the total. Large sample sizes in the Monte Carlo analyses of remote sensing data generated low estimates of uncertainty. Lidar showed lower uncertainty and higher accuracy than QuickBird, due to high correlation of biomass to height and undercounting of trees by the crown detection algorithm. Lidar achieved uncertainties of < 1%, providing estimates of aboveground live tree carbon density (mean ± 95% confidence interval with Monte Carlo) of 82 ± 0.7 Mg ha^− 1 in Garcia-Mailliard and 140 ± 0.9 Mg ha^− 1 in North Yuba. The method that we tested, combining field measurements, Lidar, and Monte Carlo, can produce robust wall-to-wall spatial data on forest carbon.