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.     

Testing LiDAR models of fractional cover across multiple forest ecozones

Four LiDAR-based models of canopy fractional cover (FC_LiDAR) have been tested against hemispherical photography fractional cover measurements (FC_HP) and compared across five ecozones, eight forest species and multiple LiDAR survey configurations. The four models compared are based on: i) a canopy-to-total first returns ratio (FC_LiDAR(FR)) method; ii) a canopy-to-total returns ratio (FC_LiDAR(RR)); iii) an intensity return ratio (FC_LiDAR(IR)); and iv) a Beer's Law modified (two-way transmission loss) intensity return ratio (FC_LiDAR(BL)). It is found that for the entire dataset, the FC_LiDAR(RR) model demonstrates the lowest overall predictive capability of overhead FC (annulus rings 1-4) (r² = 0.70), with a slight improvement for the FC_LiDAR(FR) model (r² = 0.74). The intensity-based FC_LiDAR(IR) model displays the best results (r² = 0.78). However, the FC_LiDAR(BL) model is considered generally more useful (r² = 0.75) because the associated line of best fit passes through the origin, has a slope near unity and produces a mean estimate of FC_HP within 5%. Therefore, FC_LiDAR(BL) requires the least calibration across a broad range of forest cover types. The FC_LiDAR(FR) and FC_LiDAR(RR) models, on the other hand, were found to be sensitive to variations in both canopy height and sensor pulse repetition frequency (or pulse power); i.e. changing the repetition frequency led to a systematic shift of up to 11% in the mean FC_LiDAR(RR) estimates while it had no effect on the intensity-based FC_LiDAR(IR) or FC_LiDAR(BL) models. While the intensity-based models were generally more robust, all four models displayed at least some sensitivity to variations in canopy structural class, suggesting that some calibration of FC_LiDAR might be necessary regardless of the model used. Short (< 2 m tall) or open canopy forest plots posed the greatest challenge to accurate FC estimation regardless of the model used.     

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.     

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.     

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.     

Modeling forest songbird species richness using LiDAR-derived measures of forest structure

Conservation of biodiversity requires information at many spatial scales in order to detect and preserve habitat for many species, often simultaneously. Vegetation structure information is particularly important for avian habitat models and has largely been unavailable for large areas at the desired resolution. Airborne LiDAR, with its combination of relatively broad coverage and fine resolution provides existing new opportunities to map vegetation structure and hence avian habitat. Our goal was to model the richness of forest songbirds using forest structure information obtained from LiDAR data. In deciduous forests of southern Wisconsin, USA, we used discrete-return airborne LiDAR to derive forest structure metrics related to the height and density of vegetation returns, as well as composite variables that captured major forest structural elements. We conducted point counts to determine total forest songbird richness and the richness of foraging, nesting, and forest edge-related habitat guilds. A suite of 35 LiDAR variables were used to model bird species richness using best-subsets regression and we used hierarchical partitioning analysis to quantify the explanatory power of each variable in the multivariate models. Songbird species richness was correlated most strongly with LiDAR variables related to canopy and midstory height and midstory density (R² = 0.204, p < 0.001). Richness of species that nest in the midstory was best explained by canopy height variables (R² = 0.197, p < 0.001). Species that forage on the ground responded to mean canopy height and the height of the lower canopy (R² = 0.149, p < 0.005) while aerial foragers had higher richness where the canopy was tall and dense and the midstory more sparse (R² = 0.216, p < 0.001). Richness of edge-preferring species was greater where there were fewer vegetation returns but higher density in the understory (R² = 0.153, p < 0.005). Forest interior specialists responded positively to a tall canopy, developed midstory, and a higher proportion of vegetation returns (R² = 0.195, p < 0.001). LiDAR forest structure metrics explained between 15 and 20% of the variability in richness within deciduous forest songbird communities. This variability was associated with vertical structure alone and shows how LiDAR can provide a source of complementary predictive data that can be incorporated in models of wildlife habitat associations across broad geographical extents.     

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.     

Estimating biomass carbon stocks for a Mediterranean forest in central Spain using LiDAR height and intensity data

Biomass fractions (total aboveground, branches and foliage) were estimated from a small footprint discrete-return LiDAR system in an unmanaged Mediterranean forest in central Spain. Several biomass estimation models based on LiDAR height, intensity or height combined with intensity data were explored. Raw intensity data were normalized to a standard range in order to remove the range dependence of the intensity signal. In general terms, intensity-based models provided more accurate predictions of the biomass fractions. Height models selected were mainly based on a percentile of the height distribution. Intensity models selected included variables that consider the percentage of the intensity accumulated at different height percentiles, which implicitly take into account the height distribution. The general models derived considering all species together were based on height combined with intensity data. These models yielded R² values greater than 0.58 for the different biomass fractions considered and RMSE values of 28.89, 18.28 and 1.51 Mg ha⁻¹ for aboveground, branch and foliage biomass, respectively. Results greatly improved for species-specific models using the main species present in each plot, with R² values greater than 0.85, 0.70 and 0.90 for black pine, Spanish juniper and Holm oak, respectively, and with lower RMSE for the biomass fractions. Reductions in LiDAR point density had only a small effect on the results obtained, except for those models based on a variation of the Canopy Reflection Sum, which was weighted by the mean point density. Based on the species-specific equations derived, Holm oak dominated plots showed the highest average carbon contained by aboveground biomass and branch biomass 44.66 and 31.42 Mg ha^− 1 respectively, while for foliage biomass carbon, Spanish juniper showed the highest average value (3.04 Mg ha^− 1).     

Red tide detection and tracing using MODIS fluorescence data: A regional example in SW Florida coastal waters

Near real-time data from the MODIS satellite sensor was used to detect and trace a harmful algal bloom (HAB), or red tide, in SW Florida coastal waters from October to December 2004. MODIS fluorescence line height (FLH in W m^− 2 μm^− 1 sr^− 1) data showed the highest correlation with near-concurrent in situ chlorophyll-a concentration (Chl in mg m^− 3). For Chl ranging between 0.4 to 4 mg m^− 3 the ratio between MODIS FLH and in situ Chl is about 0.1 W m^− 2 μm^− 1 sr^− 1 per mg m^− 3 chlorophyll (Chl = 1.255 (FLH × 10)^0.86, r = 0.92, n = 77). In contrast, the band-ratio chlorophyll product of either MODIS or SeaWiFS in this complex coastal environment provided false information. Errors in the satellite Chl data can be both negative and positive (3-15 times higher than in situ Chl) and these data are often inconsistent either spatially or temporally, due to interferences of other water constituents. The red tide that formed from November to December 2004 off SW Florida was revealed by MODIS FLH imagery, and was confirmed by field sampling to contain medium (10⁴ to 10⁵ cells L^− 1) to high (> 10⁵ cells L^− 1) concentrations of the toxic dinoflagellate Karenia brevis. The FLH imagery also showed that the bloom started in mid-October south of Charlotte Harbor, and that it developed and moved to the south and southwest in the subsequent weeks. Despite some artifacts in the data and uncertainty caused by factors such as unknown fluorescence efficiency, our results show that the MODIS FLH data provide an unprecedented tool for research and managers to study and monitor algal blooms in coastal environments.     

Remote estimation of chl-a concentration in turbid productive waters — Return to a simple two-band NIR-red model?

Today the water quality of many inland and coastal waters is compromised by cultural eutrophication in consequence of increased human agricultural and industrial activities. Remote sensing is widely applied to monitor the trophic state of these waters. This study investigates the performance of near infrared-red models for the remote estimation of chlorophyll-a concentrations in turbid productive waters and evaluates several near infrared-red models developed within the last 34 years. Three models were calibrated for a dataset with chlorophyll-a concentrations from 0 to 100 mg m⁻³ and validated for independent and statistically different datasets with chlorophyll-a concentrations from 0 to 100 mg m⁻³ and 0 to 25 mg m⁻³ for the spectral bands of the MEdium Resolution Imaging Spectrometer (MERIS) and MODerate resolution Imaging Spectroradiometer (MODIS). The MERIS two-band model estimated chlorophyll-a concentrations slightly more accurately than the more complex models, with mean absolute errors of 2.3 mg m⁻³ for chlorophyll-a concentrations from 0 to 100 mg m⁻³ and 1.2 mg m⁻³ for chlorophyll-a concentrations from 0 to 25 mg m⁻³. Comparable results from several near infrared-red models with different levels of complexity, calibrated for inland and coastal waters around the world, indicate a high potential for the development of a simple universally applicable near infrared-red algorithm.     

An assessment of chlorophyll-a algorithms available for SeaWiFS in coastal and open areas of the Bay of Bengal and Arabian Sea

Three ocean colour algorithms, OC4v6, Carder and OC5 were tested for retrieving Chlorophyll-a (Chla) in coastal areas of the Bay of Bengal and open ocean areas of the Arabian Sea. Firstly, the algorithms were run using ~ 80 in situ Remote Sensing Reflectance, (R_rs(λ)) data collected from coastal areas during eight cruises from January 2000 to March 2002 and the output was compared to in situ Chla. Secondly, the algorithms were run with ~ 20 SeaWiFS R_rs(λ) and the results were compared with coincident in situ Chla. In both cases, OC5 exhibited the lowest log₁₀-RMS, bias, had a slope close to 1 and this algorithm appears to be the most accurate for both coastal and open ocean areas. Thirdly the error in the algorithms was regressed against Total Suspended Material (TSM) and Coloured Dissolved Organic Material (CDOM) data to assess the co-variance with these parameters. The OC5 error did not co-vary with TSM and CDOM. OC4v6 tended to over-estimate Chla > 2 mg m⁻³ and the error in OC4v6 co-varied with TSM. OC4v6 was more accurate than the Carder algorithm, which over-estimated Chla at concentrations > 1 mg m⁻³ and under-estimated Chla at values < 0.5 mg m⁻³. The error in Carder Chla also co-varied with TSM. The algorithms were inter-compared using > 5500 SeaWiFS R_rs(λ) data from coastal to offshore transects in the Northern Bay of Bengal. There was good agreement between OC4v6 and OC5 in open ocean waters and in coastal areas up to 2 mg m⁻³. There was a strong divergence between Carder and OC5 in open ocean and coastal waters. OC4v6 and Carder tended to over-estimate Chla in coastal areas by a factor of 2 to 3 when TSM > 25 g m⁻³. We strongly recommend the use of OC5 for coastal and open ocean waters of the Bay of Bengal and Arabian Sea. A Chla time series was generated using OC5 from 2000 to 2003, which showed that concentrations at the mouths of the Ganges reach a maxima (~ 5 mg m⁻³) in October and November and were 0.08 mg m⁻³ further offshore increasing to 0.2 mg m⁻³ during December. Similarly in early spring from February to March, Chla was 0.08 to 0.2 mg m⁻³ on the east coast of the Bay.     

Fast, prime factor, discrete Fourier transform algorithms over GF(2) for 8 ? m ? 10

In this paper it is shown that Winograd’s algorithm for computing convolutions and a fast, prime factor, discrete Fourier transform (DFT) algorithm can be modified to compute Fourier-like transforms of long sequences of 2^m − 1 points over GF(2^m), for 8 ? m ? 10. These new transform techniques can be used to decode Reed-Solomon (RS) codes of block length 2^m − 1. The complexity of this new transform algorithm is reduced substantially from more conventional methods. A computer simulation verifies these new results.     

Substrate age and precipitation effects on Hawaiian forest canopies from spaceborne imaging spectroscopy

We used spaceborne imaging spectroscopy provided by the Earth Observing-1 Hyperion sensor to quantify the relative importance of precipitation and substrate age that control ecosystem development and functioning in Metrosideros polymorpha rainforests of Hawaii. Four hyperspectral vegetation indices provided metrics of forest canopy structure, biochemistry and physiology to compare along gradients of annual rainfall (750 to > 6000 mm year^− 1) and substrate age (0 to 250,000 years). The canopy greenness index NDVI increased with annual precipitation and substrate age, but saturated in forests with rainfall of 3000 mm year^− 1. Precipitation and substrate age were roughly equal contributors to the observed greenness of the forests. A canopy water content index (NDWI) also increased with precipitation and substrate age, but did not reach a maximum until very wet (> 5000 mm year^− 1) forest conditions were encountered on the oldest substrates. The water index appears superior to the NDVI in capturing spatial and climate-substrate driven variations in canopy structure. The photochemical reflectance index (PRI) indicated highest light-use efficiency levels in canopies on the most developed substrates and at annual precipitation levels of 3-4500 mm year^− 1. A leaf carotenoid index (CRI) suggested a maximum canopy photosynthetic capacity at ∼ 4000 mm rainfall year^− 1 on the oldest substrates. These results quantify the sensitivity of rainforest canopies to changing precipitation and soil conditions, and they corroborate plot-scale analyses in native Hawaiian forests ecosystems. Structural and functional studies of remote rainforest regions are possible with spaceborne imaging spectroscopy, and could be used to understand the dynamics of rainforests with climate change.     

The m-pancycle-connectivity of a WK-Recursive network

In this paper, we study the m-pancycle-connectivity of a WK-Recursive network. We show that a WK-Recursive network with amplitude W and level L is strictly (5 × 2^L−1 − 2)-pancycle-connected for W ? 3. That is, each pair of vertices in a WK-recursive network with amplitude greater than or equal to 3 resides in a common cycle of every length ranging from 5 × 2^L−1 − 2 to N, where N is the size of the interconnection network; and the value 5 × 2^L−1 − 2 reaches the lower bound of the problem.     

Estimating Siberian timber volume using MODIS and ICESat/GLAS

Geosciences Laser Altimeter System (GLAS) space LiDAR data are used to attribute a MODerate resolution Imaging Spectrometer (MODIS) 500 m land cover classification of a 10° latitude by 12° longitude study area in south-central Siberia. Timber volume estimates are generated for 16 forest classes, i.e., four forest cover types × four canopy density classes, across this 811,414 km² area and compared with a ground-based regional volume estimate. Two regional GLAS/MODIS timber volume products, one considering only those pulses falling on slopes ≤ 10° and one utilizing all GLAS pulses regardless of slope, are generated. Using a two-phase(GLAS-ground plot) sampling design, GLAS/MODIS volumes average 163.4 ± 11.8 m³/ha across all 16 forest classes based on GLAS pulses on slopes ≤ 10° and 171.9 ± 12.4 m³/ha considering GLAS shots on all slopes. The increase in regional GLAS volume per-hectare estimates as a function of increasing slope most likely illustrate the effects of vertical waveform expansion due to the convolution of topography with the forest canopy response. A comparable, independent, ground-based estimate is 146 m³/ha [Shepashenko, D., Shvidenko, A., and Nilsson, S. (1998). Phytomass (live biomass) and carbon of Siberian forests. Biomass and Bioenergy, 14, 21-31], a difference of 11.9% and 17.7% for GLAS shots on slopes ≤ 10° and all GLAS shots regardless of slope, respectively. A ground-based estimate of total volume for the entire study area, 7.46 × 10⁹ m³, is derived using Shepashenko et al.'s per-hectare volume estimate in conjunction with forest area derived from a 1990 forest map [Grasia, M.G. (ed.). (1990). Forest Map of USSR. Soyuzgiproleskhoz, Moscow, RU. Scale: 1:2,500,000]. The comparable GLAS/MODIS estimate is 7.38 × 10⁹ m³, a difference of less than 1.1%. Results indicate that GLAS data can be used to attribute digital land cover maps to estimate forest resources over subcontinental areas encompassing hundreds of thousands of square kilometers.     

On the reproducibility of field-measured reflectance factors in the context of vegetation studies

This paper describes a study aimed at quantifying uncertainty in field measurements of vegetation canopy hemispherical conical reflectance factors (HCRF). The use of field spectroradiometers is common for this purpose, but the reliability of such measurements is still in question. In this paper we demonstrate the impact of various measurement uncertainties on vegetation canopy HCRF, using a combined laboratory and field experiment employing three spectroradiometers of the same broad specification (GER 1500). The results show that all three instruments performed similarly in the laboratory when a stable radiance source was measured (NEΔL < 1 mW m⁻² sr⁻¹ nm⁻¹ in the range of 400-1000 nm). In contrast, field-derived standard uncertainties (u = SD of 10 consecutive measurements of the same surface measured in ideal atmospheric conditions) significantly differed from the lab-based uncertainty characterisation for two targets: a control (75% Spectralon panel) and a cropped grassland surface. Results indicated that field measurements made by a single instrument of the vegetation surface were reproducible to within ± 0.015 HCRF and of the control surface to within ± 0.006 HCRF (400-1000 nm (± 1σ)). Field measurements made by all instruments of the vegetation surface were reproducible to within ± 0.019 HCRF and of the control surface to within ± 0.008 HCRF (400-1000 nm (± 1σ)). Statistical analysis revealed that even though the field conditions were carefully controlled and the absolute values of u were small, different instruments yielded significantly different reflectance values for the same target. The results also show that laboratory-derived uncertainty quantities do not present a useful means of quantifying all uncertainties in the field. The paper demonstrates a simple method for u characterisation, using internationally accepted terms, in field scenarios. This provides an experiment-specific measure of u that helps to put measurements in context and forms the basis for comparison with other studies.     

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.     

Simulated impact of sample plot size and co-registration error on the accuracy and uncertainty of LiDAR-derived estimates of forest stand biomass

Regression has been widely applied in Light Detection And Ranging (LiDAR) remote sensing to spatially extend predictions of total aboveground biomass (TAGB) and other biophysical properties over large forested areas. Sample (field) plot size has long been considered a key sampling design parameter and focal point for optimization in forest surveys, because of its impact on sampling effort and the estimation accuracy of forest inventory attributes. In this study, we demonstrate how plot size and co-registration error interact to influence the estimation of LiDAR canopy height and density metrics, regression model coefficients, and the prediction accuracy of least-squares estimators of TAGB. We made use of simulated forest canopies and synthetic LiDAR point clouds, so that we could maintain strict control over the spatial scale and complexity of forest scenes, as well as the magnitude and type of planimetric error inherent in ground-reference and LiDAR datasets. Our results showed that predictions of TAGB improved markedly as plot size increased from 314 (10 m radius) to 1964 m² (25 m radius). The co-registration error (spatial overlap) between ground-reference and LiDAR samples negatively impacted the estimation of LiDAR metrics, regression model fit, and the prediction accuracy of TAGB. We found that larger plots maintained a higher degree of spatial overlap between ground-reference and LiDAR datasets for any given GPS error, and were therefore more resilient to the ill effects of co-registration error compared to small plots. The impact of co-registration error was more pronounced in tall, spatially heterogeneous stands than short, homogeneous stands. We identify and briefly discuss three possible ways that LiDAR data could be used to optimize plot size, sample selection, and the deployment of GPS resources in forest biomass surveys.     

Remote sensing of forest biophysical variables using HyMap imaging spectrometer data

This study systematically evaluated linear predictive models between vegetation indices (VI) derived from radiometrically corrected airborne imaging spectrometer (HyMap) data and field measurements of biophysical forest stand variables (n=40). Ratio-based and soil-line-related broadband VI were calculated after HyMap reflectance had been spectrally resampled to Landsat TM channels. Hyperspectral VI involved all possible types of two-band combinations of ratio VI (RVI) and perpendicular VI (PVI) and the red edge inflection point (REIP) computed from two techniques, inverted Gaussian Model and Lagrange Interpolation. Cross-validation procedure was used to assess the prediction power of the regression models. Analyses were performed on the entire data set or on subsets stratified according to stand age. A PVI based on wavebands at 1088 nm and 1148 nm was linearly related to leaf area index (LAI) (R²=0.67, RMSE=0.69 m² m⁻² (21% of the mean); after removal of one forest stand subjected to clearing measures: R²=0.77, RMSE=0.54 m² m⁻² (17% of the mean). A PVI based on wavebands at 885 nm and 948 nm was linearly related to the crown volume (VOL) (R²=0.79, RMSE=0.52). VOL was derived from measured biophysical variables through factor analysis (varimax rotation). The study demonstrates that for hyperspectral image data, linear regression models can be applied to quantify LAI and VOL with good accuracy. For broadband multispectral data, the accuracy was generally lower. It can be stated that the hyperspectral data set contains more information relevant to the estimation of the forest stand variables LAI and VOL than multispectral data. When the pooled data set was analysed, soil-line-related VI performed better than ratio-based VI. When age classes were analysed separately, hyperspectral VI performed considerably better than broadband VI. Best hyperspectral VI in relation with LAI were typically based on wavebands related to prominent water absorption features. Such VI are related to the total amount of canopy water; as the leaf water content is considered to be relatively constant in the study area, variations of LAI are retrieved.