Estimating above-ground biomass in young forests with airborne laser scanning

Total above-ground biomass of spruce, pine and birch was estimated in three different field datasets collected in young forests in south-east Norway. The mean heights ranged from 1.77 to 9.66 m. These field data were regressed against metrics derived from canopy height distributions generated from airborne laser scanner (ALS) data with a point density of 0.9–1.2 m^?2. The field data consisted of 79 plots with size 200–232.9 m² and 20 stands with an average size of 3742 m². Total above-ground biomass ranged from 2.27 to 90.42 Mg ha^?1. The influences of (1) regression model form, (2) canopy threshold value and (3) tree species on the relationships between biomass and ALS-derived metrics were assessed. The analysed model forms were multiple linear models, models with logarithmic transformation of the response and explanatory variables, and models with square root transformation of the response. The different canopy thresholds considered were fixed values of 0.5, 1.3 and 2.0 m defining the limit between laser canopy echoes and below-canopy echoes. The proportion of explained variability of the estimated models ranged from 60% to 83%. Tree species had a significant influence on the models. For given values of the ALS-derived metrics related to canopy height and canopy density, spruce tended to have higher above-ground biomass values than pine and deciduous species. There were no clear effects of model form and canopy threshold on the accuracy of predictions produced by cross validation of the various models, but there is a risk of heteroskedasticity with linear models. Cross validation revealed an accuracy of the root mean square error (RMSE) ranging from 3.85 to 13.9 Mg ha^?1, corresponding to 22.6% to 48.1% of mean field-measured biomass. It was concluded that airborne laser scanning has a potential for predicting biomass in young forest stands (> 0.5 ha) with an accuracy of 20–30% of mean ground value.     

Estimating spruce and pine biomass with interferometric X-band SAR

The primary aim of this study was to investigate the suitability of interferometric X-band SAR (InSAR) for inventory of boreal forest biomass. We investigated the relationship between SRTM X-band InSAR height and above-ground biomass in a study area in southern Norway. We generated biomass reference data for each SRTM pixel from a field inventory in combination with airborne laser scanning (ALS). One set of forest inventory plots served for calibrating ALS based biomass models, and another set of field plots was used to validate these models. The biomass values obtained in this way ranged up to 250 t/ha at the stand level. The relationship between biomass and InSAR height was linear, no apparent saturation effect was present, and the accuracy was high (RMSE = 19%). The relationship differed between Norway spruce and Scots pine, where an increase in InSAR height of 1 m corresponded to an increase in biomass of 9.9 and 7.0 t/ha, respectively. Using a high-quality terrain model from ALS enabled biomass to be estimated with a higher accuracy as compared to using a terrain model from topographic maps. Interferometric X-band SAR appears to be a promising method for forest biomass monitoring.     

Mapping U.S. forest biomass using nationwide forest inventory data and moderate resolution information

A spatially explicit dataset of aboveground live forest biomass was made from ground measured inventory plots for the conterminous U.S., Alaska and Puerto Rico. The plot data are from the USDA Forest Service Forest Inventory and Analysis (FIA) program. To scale these plot data to maps, we developed models relating field-measured response variables to plot attributes serving as the predictor variables. The plot attributes came from intersecting plot coordinates with geospatial datasets. Consequently, these models serve as mapping models. The geospatial predictor variables included Moderate Resolution Imaging Spectrometer (MODIS)-derived image composites and percent tree cover; land cover proportions and other data from the National Land Cover Dataset (NLCD); topographic variables; monthly and annual climate parameters; and other ancillary variables. We segmented the mapping models for the U.S. into 65 ecologically similar mapping zones, plus Alaska and Puerto Rico. First, we developed a forest mask by modeling the forest vs. nonforest assignment of field plots as functions of the predictor layers using classification trees in See5©. Secondly, forest biomass models were built within the predicted forest areas using tree-based algorithms in Cubist©. To validate the models, we compared field-measured with model-predicted forest/nonforest classification and biomass from an independent test set, randomly selected from available plot data for each mapping zone. The estimated proportion of correctly classified pixels for the forest mask ranged from 0.79 in Puerto Rico to 0.94 in Alaska. For biomass, model correlation coefficients ranged from a high of 0.73 in the Pacific Northwest, to a low of 0.31 in the Southern region. There was a tendency in all regions for these models to over-predict areas of small biomass and under-predict areas of large biomass, not capturing the full range in variability. Map-based estimates of forest area and forest biomass compared well with traditional plot-based estimates for individual states and for four scales of spatial aggregation. Variable importance analyses revealed that MODIS-derived information could contribute more predictive power than other classes of information when used in isolation. However, the true contribution of each variable is confounded by high correlations. Consequently, excluding any one class of variables resulted in only small effects on overall map accuracy. An estimate of total C pools in live forest biomass of U.S. forests, derived from the nationwide biomass map, also compared well with previously published estimates.     

Comparing regression methods in estimation of biophysical properties of forest stands from two different inventories using laser scanner data

Mean tree height, dominant height, mean diameter, stem number, basal area, and timber volume of 233 field sample plots were estimated from various canopy height and canopy density metrics— derived by means of a small-footprint laser scanner over young and mature forest stands— using ordinary least-squares (OLS) regression analysis, seemingly unrelated regression (SUR), and partial least-squares (PLS) regression. The sample plots were distributed systematically throughout two separate inventory areas with size 1000 and 6500 ha, respectively. The plots were divided into three predefined strata. Separate regression models were estimated for each inventory as well as common models utilizing the plots of both inventories simultaneously. In the models estimated by combining data from the two areas, the statistical effect of inventory was found to be significant (p<0.05) in the mean height models only. A total of 115 test stands and plots with size 0.3-11.7 ha were used to validate the estimated regression models. The bias and standard deviations (parenthesized) of the differences between predicted and ground reference values of mean height, dominant height, mean diameter, stem number, basal area, and volume were −5.5% to 4.7% (3.1-7.3%), −6.0% to 0.4% (2.9-8.2%), −0.2% to 7.9% (5.5-15.8%), −21.3% to 12.5% (13.4-29.3%), −7.3% to 8.4% (7.1-13.6%), and −3.9% to 10.1% (8.3-14.9%), respectively. It was revealed that only minor discrepancies occurred between the three investigated estimation techniques. None of the techniques provided predicted values that were superior to the other techniques over all combinations of strata and variables.     

基于光学与SAR因子的森林生物量多元回归估算

基于福建省Landsat-8 OLI影像，利用混合像元分解模型从实测样地数据中筛选出“纯净”的植被像元，并将筛选出的样地分为针叶林、阔叶林和混交林3种植被类型，依次提取3种不同植被类型“纯净”植被像元的树高、林龄、坡度属性信息以及对应的光学NDVI、RVI植被因子和合成孔径雷达（SAR）HH、HV极化后向散射因子，分别构成不同植被类型的“含光学特征多元因子”（NDVI、RVI、树高、林龄、坡度）和“含SAR特征多元因子”（HH、HV、树高、林龄、坡度），开展对比研究。采用含光学特征的多元因子回归模型先估测不同植被类型的森林叶生物量，然后根据叶生物量与地上生物量的关系间接估测森林地上生物量。同时，采用含SAR特征的多元因子回归模型直接估测森林的地上生物量。最后，对比分析这两组多元回归模型的估测精度。结果表明：不同植被类型的含光学特征多元回归模型的验证精度（针叶林：R²为0.483，RMSE为29.522 t/hm²；阔叶林：R²为0.470，RMSE为21.632 t/hm²；混交林：R²为0.351，RSME为25.253 t/hm²）比含SAR特征多元回归模型的验证精度（针叶林：R²为0.319，RMSE为28.352 t/hm²；阔叶林：R²为0.353，RMSE为18.991t/hm²；混交林：R²为0.281，RMSE为26.637 t/hm²）略高，说明在福建省森林生物量估算中采用含光学特征的多元回归模型（先估测叶生物量进而间接估测地上生物量）比利用含SAR特征的多元回归模型（直接估测地上生物量）更具优势。     

基于光学与SAR因子的森林生物量多元回归估算

基于福建省Landsat-8 OLI影像,利用混合像元分解模型从实测样地数据中筛选出“纯净”的植被像元,并将筛选出的样地分为针叶林、阔叶林和混交林3种植被类型,依次提取3种不同植被类型“纯净”植被像元的树高、林龄、坡度属性信息以及对应的光学NDVI、RVI植被因子和合成孔径雷达（SAR）HH、HV极化后向散射因子,分别构成不同植被类型的“含光学特征多元因子”（NDVI、RVI、树高、林龄、坡度）和“含SAR特征多元因子”（HH、HV、树高、林龄、坡度）,开展对比研究。采用含光学特征的多元因子回归模型先估测不同植被类型的森林叶生物量,然后根据叶生物量与地上生物量的关系间接估测森林地上生物量。同时,采用含SAR特征的多元因子回归模型直接估测森林的地上生物量。最后,对比分析这两组多元回归模型的估测精度。结果表明：不同植被类型的含光学特征多元回归模型的验证精度（针叶林：R²为0.483,RMSE为29.522 t/hm²;阔叶林：R²为0.470,RMSE为21.632 t/hm²;混交林：R²为0.351,RSME为25.253 t/hm²）比含SAR特征多元回归模型的验证精度（针叶林：R²为0.319,RMSE为28.352 t/hm²;阔叶林：R²为0.353,RMSE为18.991t/hm²;混交林：R²为0.281,RMSE为26.637 t/hm²）略高,说明在福建省森林生物量估算中采用含光学特征的多元回归模型（先估测叶生物量进而间接估测地上生物量）比利用含SAR特征的多元回归模型（直接估测地上生物量）更具优势。     

Satellite estimation of tropical secondary forest above-ground biomass: Data from Brazil and Bolivia

This paper reports on a test of the ability to estimate above-ground biomass of tropical secondary forest from canopy spectral reflectance using satellite optical data. Landsat Thematic Mapper data were acquired concurrent with field surveys conducted in secondary forest fallows near Manaus, Brazil and Santa Cruz de la Sierra, Bolivia. Measurements of age and above-ground live biomass were made in 34 regrowth stands. Satellite data were converted to surface reflectances and compared with regrowth stand age, biomass and structural variables. Among the Brazilian stands, significant relationships were observed between middle-infrared reflectance and stand age, height, volume and biomass. The canopy reflectance-biomass relationship saturated at around 15.0 kg m^-2, or over 15 years of age (r > 0.80, p < 0.01). In the Bolivian study area, no significant relationship between canopy spectral reflectance and biomass was observed. These contrasting results are probably caused by a low Sun angle during the satellite measurements from Bolivia. However, regrowth structural and general compositional differences between the two study areas could explain the lack of a significant relationship in Bolivia. The results demonstrate a current potential for biomass estimation of secondary forests with satellite optical data in some, but not all, tropical regions. A discussion of the potential for regional extrapolation of spectral relationships and future satellite imagery is included.     

Lidar remote sensing of forest biomass: A scale-invariant estimation approach using airborne lasers

Researchers in lidar (Light Detection And Ranging) strive to search for the most appropriate laser-based metrics as predictors in regression models for estimating forest structural variables. Many previously developed models are scale-dependent that need to be fitted and then applied both at the same scale or pixel size. The objective of this paper is to develop methods for scale-invariant estimation of forest biomass using lidar data. We proposed two scale-invariant models for biomass: a linear functional model and an equivalent nonlinear model that use lidar-derived canopy height distributions (CHD) and canopy height quantile functions (CHQ) as predictors, respectively. The two models are called functional regression models because the predictors CHD and CHQ are themselves functions or functional data. The model formulation was justified mathematically under moderate assumptions. We also created a fine-resolution biomass map by mapping individual tree component biomass in a temperate forest of eastern Texas with a lidar tree-delineation approach. The map was used as reference data to synthesize training and test datasets at multiple scales for validating the two scale-invariant models. Results suggest that the models can accurately predict biomass and yield consistent predictive performances across a variety of scales with an R² ranging from 0.80 to 0.95 (RMSE: from 14. 3 Mg/ha to 33.7 Mg/ha) among all the fitted models. Results also show that a training data size of around 50 plots or less was enough to guarantee a good fitting of the linear functional model. Our findings demonstrate the effectiveness of CHD and CHQ as lidar metrics for estimating biomass as well as the capability of lidar for mapping biomass at a range of scales. The functional regression models of this study are useful for lidar-based forest inventory tasks where the analysis units vary in size and shape. They also hold promise for estimating other forest characteristics such as below-ground biomass, timber volume, crown fuel weight, and Leaf Area Index.     

Estimation of above-ground forest biomass using metrics based on Gaussian decomposition of waveform lidar data

Large-footprint waveform light detection and ranging (lidar) data have been widely used in above-ground forest biomass estimation. Waveform metrics derived from basic statistics (e.g. percentile of energy) of the lidar waveform, such as canopy height and height of median energy, have been applied to biomass estimation in numerous studies. In this study, a set of metrics based on Gaussian decomposition (GD) results were developed and evaluated for forest above-ground biomass estimation using NASA’s laser vegetation imaging sensor (LVIS) data. The GD metrics were designed to explicitly incorporate lidar intensity and vertical structures of canopy layers for biomass estimation. The proposed GD metrics used information related to the above-ground height of each Gaussian centroid and the Gaussian area index (GAI), where GAI is the area covered by a Gaussian function. Two types of novel GD metrics were developed: (1) percentile-height GAI metrics expressing the GAI summation of a subset of Gaussian centroids located within a certain percentile height range; and (2) height-weighted GAI metrics, a summation of GAIs of a waveform weighted by the corresponding heights of their Gaussian centroids. A biomass regression model was built by eight newly developed GD metrics using GAI information and five pre-existing GD-derived metrics that have not previously been used for biomass estimation. The performance of the regression model was then compared to another regression model using 12 previously published metrics (non-GD metrics). The Random Forests (RF) regression algorithm was employed for predicting biomass. The RF out-of-bag results indicated that above-ground biomass estimations using GD metrics achieved significantly better results than those derived from non-GD metrics for deciduous plots (19% lower root mean square error (RMSE), 25% higher coefficient of determination (R²), and marginally better results in coniferous plots (4% lower RSME, 6% higher R²). The combination of GD and non-GD metrics achieved comparable biomass estimation results to the model using exclusively GD metrics. GD metrics also showed strong correlation with forest attributes such as mean diameter at breast height (DBH) and stem density. This study contributes to the usage of GD results for accurate estimation of forest above-ground biomass in large-footprint lidar waveform data in temperate deciduous forests, because temperate deciduous forests have been proved challenging in regard to lidar-derived biomass estimations.     

Forest mapping by geoinformatics for landscape fire behaviour modelling in coastal forests,Greece

This study aims at quantifying and mapping fire-related characteristics of forest structure through field inventories, statistics, remote sensing, and geographical information systems in the island of Lesvos, northeast Aegean Sea, Greece. Simulation of fire behaviour requires forest biomass inputs that describe surface fuel types/models along with canopy fuel properties, such as canopy cover, stand height, crown base height, and crown bulk density, to accurately predict surface and crown fire spread and spotting potential. Forest canopy characteristics and other vegetation attributes were sampled and derived in over 100 field plots, the majority of which were located in coastal pine forest stands. Regression models involving four dependent forest stand variables (stand height, canopy cover, crown base height, and crown bulk density) were developed using generalized additive models. The values of adjusted R ² were 0.72 for stand height, 0.68 for canopy cover, 0.51 for crown base height, and 0.33 for crown bulk density. These regression models were used to create forest fuel characteristics layers, which can be used as inputs to fire management applications and state-of-the-art landscape-scale fire behaviour models.     

Relationships between polarimetric SAR backscatter and forest canopy and sub-canopy biophysical properties

This paper highlights the potential of multiwaveband polarimetric SAR data for the estimation of both canopy (percentage canopy closure) and sub-canopy (stem biomass) biophysical variables of a Sitka spruce forest in upland Wales. Stand stem biomass was estimated using forest survey data on diameter at breast height (DBH) and tree height from 0.01 ha plots. Photographs of the forest canopy were taken using a camera fitted with a wide-angle fisheye lens from a number of locations within a stand. The photographs were later digitized and estimates of stand percentage canopy closure were derived using image processing software. It was found that C-band HV and VV, and L-band HV and VV polarization backscatter were significantly related to stem biomass. There was no sensitivity to percentage canopy closure using single polarization backscatter but highly significant relationships were obtained using ratios of single polarization backscatter and variables derived from the polarization signatures. The strong correlations between C-band backscatter and stem biomass indicated a relationship between the structure of the top crown layer and sub-canopy biomass.     

Estimation of tropical rain forest aboveground biomass with small-footprint lidar and hyperspectral sensors

Tropical forests are an important component of the global carbon balance, yet there is considerable uncertainty in estimates of their carbon stocks and fluxes, which are typically estimated through analysis of aboveground biomass in field plots. Remote sensing technology is critical for assessing fine-scale spatial variability of tropical forest biomass over broad spatial extents. The goal of our study was to evaluate relatively new technology, small-footprint, discrete-return lidar and hyperspectral sensors, for the estimation of aboveground biomass in a Costa Rican tropical rain forest landscape. We derived a suite of predictive metrics for field plots: lidar metrics were calculated from plot vertical height profiles and hyperspectral metrics included fraction of spectral mixing endmembers and narrowband indices that respond to photosynthetic vegetation, structure, senescence, health and water and lignin content. We used single- and two-variable linear regression analyses to relate lidar and hyperspectral metrics to aboveground biomass of plantation, managed parkland and old-growth forest plots. The best model using all 83 biomass plots included two lidar metrics, plot-level mean height and maximum height, with an r² of 0.90 and root-mean-square error (RMSE) of 38.3 Mg/ha. When the analysis was constrained to plantation plots, which had the most accurate field data, the r² of the model increased to 0.96, with RMSE of 10.8 Mg/ha (n = 32). Hyperspectral metrics provided lower accuracy in estimating biomass than lidar metrics, and models with a single lidar and hyperspectral metric were no better than the best model using two lidar metrics. These results should be viewed as an initial assessment of using these combined sensors to estimate tropical forest biomass; hyperspectral data were reduced to nine indices and three spectral mixture fractions, lidar data were limited to first-return canopy height, sensors were flown only once at different seasons, and we explored only linear regression for modeling. However, this study does support conclusions from studies at this and other climate zones that lidar is a premier instrument for mapping biomass (i.e., carbon stocks) across broad spatial scales.     

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.     

Technical note: Canopy height models and airborne lasers to estimate forest biomass: Two problems

Three ground datasets were used to simulate the canopy height characteristics of tropical forests in Costa Rica for the purposes of forest biomass estimation. The canopy height models (CHMs) were used in conjunction with airborne laser data. Gross biomass estimation errors on the order of 50-90% arose in two of the three analyses. The characteristics of the datasets and the biomass estimation procedure are reviewed to identify sources of error. In one dataset, the width of the fixed-area ground plots were small enough (5 m) that significant portions of the overstory canopy above the plots were not accounted for in the ground samples. The use of mapped stand data from thin ground plots resulted in inaccurate CHMs, which in turn lead to gross overestimates of forest biomass (about 90% larger than the ground reference value). CHMs generated using mensuration data collected on thin, fixed-area plots may significantly underestimate the average canopy height and crown closure actually found on that plot. This underestimation problem is directly related to plot width. Below a critical threshold, the thinner the plot, the greater the underestimation bias. In the second dataset, it is believed that lower-than-normal rainfalls at the beginning and end of the wet season may have produced a forest canopy with reduced leaf area. The airborne laser pulses penetrated further into the canopy, resulting in airborne laser estimates of forest biomass which grossly underestimated reference values by about 50%. Changing canopy conditions (e.g. leaf loss due to drought, insect defoliation, storm damage) can affect the accuracy of a CHM. Sources of error are reported in order to forewarn those researchers who produce and utilize canopy height models.     

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.     

Forest biomass mapping from lidar and radar synergies

The use of lidar and radar instruments to measure forest structure attributes such as height and biomass at global scales is being considered for a future Earth Observation satellite mission, DESDynI (Deformation, Ecosystem Structure, and Dynamics of Ice). Large footprint lidar makes a direct measurement of the heights of scatterers in the illuminated footprint and can yield accurate information about the vertical profile of the canopy within lidar footprint samples. Synthetic Aperture Radar (SAR) is known to sense the canopy volume, especially at longer wavelengths and provides image data. Methods for biomass mapping by a combination of lidar sampling and radar mapping need to be developed.In this study, several issues in this respect were investigated using aircraft borne lidar and SAR data in Howland, Maine, USA. The stepwise regression selected the height indices rh50 and rh75 of the Laser Vegetation Imaging Sensor (LVIS) data for predicting field measured biomass with a R² of 0.71 and RMSE of 31.33 Mg/ha. The above-ground biomass map generated from this regression model was considered to represent the true biomass of the area and was used as a reference map since no better biomass map exists for the area. Random samples were taken from the biomass map and the correlation between the sampled biomass and co-located SAR signature was studied. The best models were used to extend the biomass from lidar samples into all forested areas in the study area, which mimics a procedure that could be used for the future DESDYnI mission. It was found that depending on the data types used (quad-pol or dual-pol) the SAR data can predict the lidar biomass samples with R² of 0.63-0.71, RMSE of 32.0-28.2 Mg/ha up to biomass levels of 200-250 Mg/ha. The mean biomass of the study area calculated from the biomass maps generated by lidar-SAR synergy was within 10% of the reference biomass map derived from LVIS data. The results from this study are preliminary, but do show the potential of the combined use of lidar samples and radar imagery for forest biomass mapping. Various issues regarding lidar/radar data synergies for biomass mapping are discussed in the paper.     

Estimating forest canopy fuel parameters using LIDAR data

Fire researchers and resource managers are dependent upon accurate, spatially-explicit forest structure information to support the application of forest fire behavior models. In particular, reliable estimates of several critical forest canopy structure metrics, including canopy bulk density, canopy height, canopy fuel weight, and canopy base height, are required to accurately map the spatial distribution of canopy fuels and model fire behavior over the landscape. The use of airborne laser scanning (LIDAR), a high-resolution active remote sensing technology, provides for accurate and efficient measurement of three-dimensional forest structure over extensive areas. In this study, regression analysis was used to develop predictive models relating a variety of LIDAR-based metrics to the canopy fuel parameters estimated from inventory data collected at plots established within stands of varying condition within Capitol State Forest, in western Washington State. Strong relationships between LIDAR-derived metrics and field-based fuel estimates were found for all parameters [sqrt(crown fuel weight): R²=0.86; ln(crown bulk density): R²=0.84; canopy base height: R²=0.77; canopy height: R²=0.98]. A cross-validation procedure was used to assess the reliability of these models. LIDAR-based fuel prediction models can be used to develop maps of critical canopy fuel parameters over forest areas in the Pacific Northwest.     

Estimating aboveground forest biomass from canopy reflectance model inversion in mountainous terrain

The amount and spatial distribution of aboveground forest biomass (AGB) are required inputs to forest carbon budgets and ecosystem productivity models. Satellite remote sensing offers distinct advantages for large area and multi-temporal applications, however, conventional empirical methods for estimating forest canopy structure and AGB can be difficult in areas of high relief and variable terrain. This paper introduces a new method for obtaining AGB from forest structure estimates using a physically-based canopy reflectance (CR) model inversion approach. A geometric-optical CR model was run in multiple forward mode (MFM) using SPOT-5 imagery to derive forest structure and biomass at Kananaskis, Alberta in the Canadian Rocky Mountains. The approach first estimates tree crown dimensions and stem density for satellite image pixels which are then related to tree biomass and AGB using a crown spheroid surface area approach. MFM estimates of AGB were evaluated for 36 deciduous (trembling aspen) and conifer (lodgepole pine) field validation sites and compared against spectral mixture analysis (SMA) and normalised difference vegetation index (NDVI) biomass predictions from atmospherically and topographically corrected (SCS+C) imagery. MFM provided the lowest error for all validation plots of 31.7 tonnes/hectare (t/ha) versus SMA (32.6 t/ha error) and NDVI (34.7 t/ha) as well as for conifer plots (MFM: 23.0 t/ha; SMA 27.9 t/ha; NDVI 29.7 t/ha) but had higher error than SMA and NDVI for deciduous plots (by 4.5 t/ha and 2.1 t/ha, respectively). The MFM approach was considerably more stable over the full range of biomass values (67 to 243 t/ha) measured in the field. Field plots with biomass > 1 standard deviation from the field mean (over 30% of plots) had biomass estimation errors of 37.9 t/ha using MFM compared with 65.5 t/ha and 67.5 t/ha error from SMA and NDVI, respectively. In addition to providing more accurate overall results and greater stability over the range of biomass values, the MFM approach also provides a suite of other biophysical structural outputs such as density, crown dimensions, LAI, height and sub-pixel scale fractions. Its explicit physical-basis and minimal ground data requirements are also more appropriate for larger area, multi-scene, multi-date applications with variable scene geometry and in high relief terrain. MFM thus warrants consideration for applications in mountainous and other, less complex terrain for purposes such as forest inventory updates, ecological modeling and terrestrial biomass and carbon monitoring studies.     

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.     

Spatial analysis of radiometric fractions from high-resolution multispectral imagery for modelling individual tree crown and forest canopy structure and health

Research was conducted in a forest adjacent to an abandoned acid mine tailings site to assess forest structural health using high spatial and spectral resolution digital camera imagery. Conventional approaches to this problem involve the use image spectral information, basic spectral transformations, or occasionally spatial transformations of image brightness. This research introduces fractional textures and semivariance analysis of image fractions. They were integrated with conventional image measures in stepwise multiple regression modelling of forest structure (canopy and crown closure, stem density, tree height, crown size) and health (a visual stress index). The goal was to conduct a relative comparison of the potential of the various image variable types in modelling of forest structure and health. Analysis was conducted for both canopy (crowns and shadows) and individual tree crown sample data sets extracted from 10 nm bandwidth spectral bands at three resolutions (0.25, 0.5, 1.0 m). Spatial transformations (texture, semivariogram range) of image brightness (DN) and image fractions (IF) were consistently the most significant and first entered variables in the best models of the forest parameters. At the canopy-scale, despite a limited number of available plots (6), stable models were produced that demonstrated the potential for spatially transformed variables. Semivariogram range explained 88% of the total variation of 9 of the 18 models and represented 56% of the variables used in all models while texture variables explained 51% of model variance in 8 of the 18 models and represented 40% of the variables used. At the tree crown scale (n=31), 88% of the total variation of six of eight models was explained by texture variables and 6% by semivariogram variables. DN and IF variables that were not spatially transformed contributed little to the models at both scales. They represented 4% and 6%, respectively, of the variables used in all models. Spatial information in image fractions and image brightness has proven to be more significant than spectral information in these analyses. Of the spatial resolutions evaluated, 0.5 m consistently produced similar or better models than those using the 0.25 or 1.0 m resolutions. These results demonstrate the potential for integration of spatial transforms of image fractions and raw brightness in high-resolution modelling of forest structure and health.