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.     

Lidar-based mapping of leaf area index and its use for validating GLOBCARBON satellite LAI product in a temperate forest of the southern USA

Lidar provides enhanced abilities to remotely map leaf area index (LAI) with improved accuracies. We aim to further explore the capability of discrete-return lidar for estimating LAI over a pine-dominated forest in East Texas, with a secondary goal to compare the lidar-derived LAI map and the GLOBCARBON moderate-resolution satellite LAI product. Specific problems we addressed include (1) evaluating the effects of analysts and algorithms on in-situ LAI estimates from hemispherical photographs (hemiphoto), (2) examining the effectiveness of various lidar metrics, including laser penetration, canopy height and foliage density metrics, to predict LAI, (3) assessing the utility of integrating Quickbird multispectral imagery with lidar for improving the LAI estimate accuracy, and (4) developing a scheme to co-register the lidar and satellite LAI maps and evaluating the consistency between them. Results show that the use of different analysts or algorithms in analyzing hemiphotos caused an average uncertainty of 0.35 in in-situ LAI, and that several laser penetration metrics in logarithm models were more effective than other lidar metrics, with the best one explaining 84% of the variation in the in-situ LAI (RMSE = 0.29 LAI). The selection of plot size and height threshold in calculating laser penetration metrics greatly affected the effectiveness of these metrics. The combined use of NDVI and lidar metrics did not significantly improve estimation over the use of lidar alone. We also found that mis-registration could induce a large artificial discrepancy into the pixelwise comparison between the coarse-resolution satellite and fine-resolution lidar-derived LAI maps. By compensating for a systematic sub-pixel shift error, the correlation between two maps increased from 0.08 to 0.85 for pines (n = 24 pixels). However, the absolute differences between the two LAI maps still remained large due to the inaccuracy in accounting for clumping effects. Overall, our findings imply that lidar offers a superior tool for mapping LAI at local to regional scales as compared to optical remote sensing, accuracies of lidar-estimate LAI are affected not only by the choice of models but also by the absolute accuracy of in-situ reference LAI used for model calibration, and lidar-derived LAI maps can serve as reliable references for validating moderate-resolution satellite LAI products over large areas.     

Estimating aboveground carbon in a catchment of the Siberian forest tundra: Combining satellite imagery and field inventory

This study was part of an interdisciplinary research project on soil carbon and phytomass dynamics of boreal and arctic permafrost landscapes. The 45 ha study area was a catchment located in the forest tundra in northern Siberia, approximately 100 km north of the Arctic Circle.The objective of this study was to estimate aboveground carbon (AGC) and assess and model its spatial variability. We combined multi-spectral high resolution remote sensing imagery and sample based field inventory data by means of the k-nearest neighbor (k-NN) technique and linear regression.Field data was collected by stratified systematic sampling in August 2006 with a total sample size of n = 31 circular nested sample plots of 154 m² for trees and shrubs and 1 m² for ground vegetation. Destructive biomass samples were taken on a sub-sample for fresh weight and moisture content. Species-specific allometric biomass models were constructed to predict dry biomass from diameter at breast height (dbh) for trees and from elliptic projection areas for shrubs.Quickbird data (standard imagery product), acquired shortly before the field campaign and archived ASTER data (Level-1B product) of 2001 were geo-referenced, converted to calibrated radiances at sensor and used as carrier data. Spectral information of the pixels which were located in the inventory plots were extracted and analyzed as reference set. Stepwise multiple linear regression was applied to identify suitable predictors from the set of variables of the original satellite bands, vegetation indices and texture metrics. To produce thematic carbon maps, carbon values were predicted for all pixels of the investigated satellite scenes. For this prediction, we compared the kNN distance-weighted classifier and multiple linear regression with respect to their predictions.The estimated mean value of aboveground carbon from stratified sampling in the field is 15.3 t/ha (standard error SE = 1.50 t/ha, SE% = 9.8%). Zonal prediction from the k-NN method for the Quickbird image as carrier is 14.7 t/ha with a root mean square error RMSE = 6.42 t/ha, RMSE_r = 44%) resulting from leave-one-out cross-validation. The k-NN-approach allows mapping and analysis of the spatial variability of AGC. The results show high spatial variability with AGC predictions ranging from 4.3 t/ha to 28.8 t/ha, reflecting the highly heterogeneous conditions in those permafrost-influenced landscapes. The means and totals of linear regression and k-NN predictions revealed only small differences but some regional distinctions were recognized in the maps.     

Multi-Scale Adaptive Sampling with Mobile Agents for Mapping of Forest Fires

The use of robotics in distributed monitoring applications requires wireless sensors that are deployed efficiently. A very important aspect of sensor deployment includes positioning them for sampling at locations most likely to yield information about the spatio-temporal field of interest, for instance, the spread of a forest fire. In this paper, we use mobile robots (agents) that estimate the time-varying spread of wildfires using a distributed multi-scale adaptive sampling strategy. The proposed parametric sampling algorithm, “EKF-NN-GAS” is based on neural networks, the extended Kalman filter (EKF), and greedy heuristics. It combines measurements arriving at different times, taken at different scale lengths, such as from ground, airborne, and spaceborne observation platforms. One of the advantages of our algorithm is the ability to incorporate robot localization uncertainty in addition to sensor measurement and field parameter uncertainty into the same EKF model. We employ potential fields, generated naturally from the estimated fire field distribution, in order to generate fire-safe trajectories that could be used to rescue vehicles and personnel. The covariance of the EKF is used as a quantitative information measure for sampling locations most likely to yield optimal information about the sampled field distribution. Neural net training is used infrequently to generate initial low resolution estimates of the fire spread parameters. We present simulation and experimental results for reconstructing complex spatio-temporal forest fire fields “truth models”, approximated by radial basis function (RBF) parameterizations. When compared to a conventional raster scan approach, our algorithm shows a significant reduction in the time necessary to map the fire field.     

Methods for evaluating volunteers’ contributions in a deforestation detection citizen science project

Today, due to the availability of free remote sensing data, efficient algorithms for image classification and increased connectivity and computing power, together with international policy initiatives, such as the United Nations Programme on Reducing Emissions from Deforestation and Forest Degradation (UN-REDD), more and more countries are investing in their own national forest monitoring schemes. However, tropical forests remain under threat worldwide. Recently, a citizen science project that enables citizens around the globe to be involved in forest monitoring tasks has been proposed, called “ForestWatchers” (www.forestwatchers.net). Its main goal is to allow volunteers (many of them with no scientific training) around the globe, with their own smartphones, tablets and notebooks, review satellite images of forested regions and confirm whether automatic assignments of forested and deforested regions are correct. Inspected images are then sent to a central database where the results are integrated to generate up-to-date deforestation maps. This approach offers a low-cost way to both strengthen the scientific infrastructure and engage members of the public in science. Here, we describe the methods developed within the scope of the ForestWatchers project to assess the volunteers’ performance. These tools have been evaluated with data of two of the project’s preliminary tasks. The first, called “BestTile”, asks volunteers to select which of several images of the same area has the least cloud cover, while in the second, called “Deforestation”, volunteers draw polygons on satellite images delimiting areas they believe have been deforested. The results from more than 500 volunteers show that using simple statistical tests, it is possible to achieve a triple goal: to increase the overall efficiency of the data collecting tasks by reducing the required number of volunteers per task, to identify malicious behavior and outliers, and to motivate volunteers to continue their contributions.     

Using ICESat's Geoscience Laser Altimeter System (GLAS) to assess large-scale forest disturbance caused by hurricane Katrina

In 2005, hurricane Katrina resulted in a large disturbance to U.S. forests. Recent estimates of damage from hurricane Katrina have relied primarily on optical remote sensing and field data. This paper is the first large-scale study to use satellite-based lidar data to quantify changes in forest structure from that event. GLAS data for the years prior to and following hurricane Katrina were compared to wind speed, forest cover, and damage data to assess the adequacy of sensor sampling, and to estimate changes in Mean Canopy Height (MCH) over all areas that experienced tropical force winds and greater. Statistically significant decreases in MCH post-Katrina were found to increase with wind intensity: Tropical Storm ?MCH = − 0.5 m, Category 1 ?MCH = − 2 m, and Category 2 ?MCH = − 4 m. A strong relationship was also found between changes in non-photosynthetic vegetation (?NPV), a metric previously shown to be related to storm damage, and post-storm MCH. The season of data acquisition was shown to influence calculations of MCH and MCH loss, but did not preclude the detection of major large-scale patterns of damage. Results from this study show promise for using space-borne lidar for large-scale assessments of forest disturbance, and highlight the need for future data on vegetation structure from space.     

Sensitivity of SAR data to post-fire forest regrowth in Mediterranean and boreal forests

Disturbed forests may need decades to reach a mature stage and optically-based vegetation indices are usually poorly suited for monitoring purposes due to the rapid saturation of the signal with increasing canopy cover. Spaceborne synthetic aperture radar (SAR) data provide an alternate monitoring approach since the backscattered microwave energy is sensitive to the vegetation structure. Images from two regions in Spain and Alaska were used to analyze SAR metrics (cross-polarized backscatter and co-polarized interferometric coherence) from regrowing forests previously affected by fire. TerraSAR-X X-band backscatter showed the lowest sensitivity to forest regrowth, with the average backscatter increasing by 1-2 dB between the most recent fire scar and the unburned forest. Increased sensitivity (around 3-4 dB) was observed for C-band Envisat Advanced Synthetic Aperture (ASAR) backscatter. The Advanced Land Observing Satellite (ALOS) Phased Array-type L-band Synthetic Aperture Radar (PALSAR) L-band backscatter presented the highest dynamic range from unburned to recently burned forests (approximately 8 dB). The interferometric coherence showed low sensitivity to forest regrowth at all SAR frequencies. For Mediterranean forests, five phases of forest regrowth were discerned whereas for boreal forest, up to four different regrowth phases could be discerned with L-band SAR data. In comparison, the Normalized Difference Vegetation Index (NDVI) provided reliable differentiation only for the most recent development stages. The results obtained were consistent in both environments.     

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

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

Broadband, red-edge information from satellites improves early stress detection in a New Mexico conifer woodland

Multiple plant stresses can affect the health, esthetic condition, and timber harvest value of conifer forests. To monitor spatial and temporal dynamic forest stress conditions, timely, accurate, and cost-effective information is needed that could be provided by remote sensing. Recently, satellite imagery has become available via the RapidEye satellite constellation to provide spectral information in five broad bands, including the red-edge region (690-730 nm) of the electromagnetic spectrum. We tested the hypothesis that broadband, red-edge satellite information improves early detection of stress (as manifest by shifts in foliar chlorophyll a + b) in a woodland ecosystem relative to other more commonly utilized band combinations of red, green, blue, and near infrared band reflectance spectra. We analyzed a temporally dense time series of 22 RapidEye scenes of a piñon-juniper woodland in central New Mexico acquired before and after stress was induced by girdling. We found that the Normalized Difference Red-Edge index (NDRE) allowed stress to be detected 13 days after girdling — between and 16 days earlier than broadband spectral indices such as the Normalized Difference Vegetation Index (NDVI) and Green NDVI traditionally used for satellite based forest health monitoring. We conclude that red-edge information has the potential to considerably improve forest stress monitoring from satellites and warrants further investigation in other forested ecosystems.     

Aboveground biomass retrieval in tropical forests — The potential of combined X- and L-band SAR data use

In the context of reducing emissions from deforestation and forest degradation (REDD) and the international effort to reduce anthropogenic greenhouse gas emissions, a reliable assessment of aboveground forest biomass is a major requirement. Especially in tropical forests which store huge amounts of carbon, a precise quantification of aboveground biomass is of high relevance for REDD activities. This study investigates the potential of X- and L-band SAR data to estimate aboveground biomass (AGB) in intact and degraded tropical forests in Central Kalimantan, Borneo, Indonesia. Based on forest inventory data, aboveground biomass was first estimated using LiDAR data. These results were then used to calibrate SAR backscatter images and to upscale the biomass estimates across large areas and ecosystems. This upscaling approach not only provided aboveground biomass estimates over the whole biomass range from woody regrowth to mature pristine forest but also revealed a spatial variation due to varying growth condition within specific forest types. Single and combined frequencies, as well as mono- and multi-temporal TerraSAR-X and ALOS PALSAR biomass estimation models were analyzed for the development of accurate biomass estimations. Regarding the single frequency analysis overall ALOS PALSAR backscatter is more sensitive to AGB than TerraSAR-X, especially in the higher biomass range (> 100 t/ha). However, ALOS PALSAR results were less accurate in low biomass ranges due to a higher variance. The multi-temporal L- and X-band combined model achieved the best result and was therefore tested for its temporal and spatial transferability. The achieved accuracy for this model using nearly 400 independent validation points was r² = 0.53 with an RMSE of 79 t/ha. The model is valid up to 307 t/ha with an accuracy requirement of 50 t/ha and up to 614 t/ha with an accuracy requirement of 100 t/ha in flat terrain. The results demonstrate that direct biomass measurements based on the synergistic use of L- and X-band SAR can provide large-scale AGB estimations for tropical forests. In the context of REDD monitoring the results can be used for the assessment of the spatial distribution of the biomass, also indicating trends in high biomass ranges and the characterization of the spatial patterns in different forest types.