Investigating RaDAR-LiDAR synergy in a North Carolina pine forest

A low frequency (80-120 MHz) VHF RaDAR, BioSAR, specifically designed for forest biomass estimation and a profiling LiDAR, PALS, were flown over loblolly pine plantations in the southeastern United States. LiDAR-only, RaDAR-only, and joint LiDAR-RaDAR linear models were developed to determine if returns from two sensors could be used to estimate pine biomass more accurately and precisely than returns from either sensor alone. The best five-variable RaDAR model explained 81.8% (R²) of the stem green biomass variability, with a regression RMSE of 57.5 t/ha. The best one-variable LiDAR model explained 93.3% of the biomass variation (RMSE = 33.9 t/ha). Combining the RaDAR normalized volumetric returns with the profiling LiDAR ranging measurements did little to improve the best LiDAR-only model. The best LiDAR-RaDAR model explained 93.8% of the biomass variation (RSME = 32.7 t/ha). Cross-validation and training/test validation procedures demonstrated (1) that all models are unbiased and (2) the increased precision of the LiDAR-only and LiDAR-RaDAR models. The results of this investigation and a companion study indicate that there is little to be gained combining VHF-RaDAR volumetric returns and profiling LiDAR ranging measurements in pine forests; a LiDAR ranging system is sufficient for accurate, precise biomass estimation.     

Impact of spatial variability of tropical forest structure on radar estimation of aboveground biomass

Understanding the spatial variability of tropical forest structure and its impact on the radar estimation of aboveground biomass (AGB) is important to assess the scale and accuracy of mapping AGB with future low frequency radar missions. We used forest inventory plots in old growth, secondary succession, and forest plantations at the La Selva Biological Station in Costa Rica to examine the spatial variability of AGB and its impact on the L-band and P-band polarimetric radar estimation of AGB at multiple spatial scales. Field estimation of AGB was determined from tree size measurements and an allometric equation developed for tropical wet forests. The field data showed very high spatial variability of forest structure with no spatial dependence at a scale above 11 m in old-growth forest. Plot sizes of greater than 0.25 ha reduced the coefficients of variation in AGB to below 20% and yielded a stationary and normal distribution of AGB over the landscape. Radar backscatter measurements at all polarization channels were strongly positively correlated with AGB at three scales of 0.25 ha, 0.5 ha, and 1.0 ha. Among these measurements, PHV and LHV showed strong sensitivity to AGB < 300 Mg ha^− 1 and AGB < 150 Mg ha^− 1 respectively at the 1.0 ha scale. The sensitivity varied across forest types because of differences in the effects of forest canopy and gap structure on radar attenuation and scattering. Spatial variability of structure and speckle noise in radar measurements contributed equally to degrading the sensitivity of the radar measurements to AGB at spatial scales less than 1.0 ha. By using algorithms based on polarized radar backscatter, we estimated AGB with RMSE = 22.6 Mg ha^− 1 for AGB < 300 Mg ha^− 1 at P-band and RMSE = 23.8 Mg ha^− 1 for AGB < 150 Mg ha^− 1 at L-band and with the accuracy optimized at 1-ha scale within 95% confidence interval. By adding the forest height, estimated from the C-band Interferometry data as an independent variable to the algorithm, the AGB estimation improved beyond the backscatter sensitivity by 20% at P-band and 40% at L-band. The results suggested the estimation of AGB can be improved substantially from the fusion of lidar or InSAR derived forest height with the polarimetric backscatter.     

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.     

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.     

Model-assisted regional forest biomass estimation using LiDAR and InSAR as auxiliary data: A case study from a boreal forest area

There is a need for accurate inventory methods that produce relevant and timely information on the forest resources and carbon stocks for forest management planning and for implementation of national strategies under the United Nations Collaborative Program on Reduced Emissions from Deforestation and Forest Degradation in Developing Countries (REDD). Such methods should produce information that is consistent across various geographical scales. Airborne scanning Light Detection and Ranging (LiDAR) is among the most promising remote sensing technologies for estimation of forest resource information such as timber volume and biomass, while acquisition of three dimensional data with Interferometric Synthetic Aperture Radar (InSAR) from space is seen as a relevant option for inventory in the tropics because of its ability to “see through the clouds” and its potential for frequent updates at low costs. Based on a stratified probability sample of 201 field survey plots collected in a 960 km² boreal forest area in Norway, we demonstrate how total above-ground biomass (AGB) can be estimated at three distinct geographical levels in such a way that the estimates at a smaller level always sum up to the estimate at a larger level. The three levels are (1) a district (the entire study area), (2) a village, local community or estate level, and (3) a stand or patch level. The LiDAR and InSAR data were treated as auxiliary information in the estimation. At the two largest geographical levels model-assisted estimators were employed. A model-based estimation was conducted at the smallest level. Estimates of AGB and corresponding error estimates based on (1) the field sample survey were compared with estimates obtained by using (2) LiDAR and (3) InSAR data as auxiliary information. For the entire study area, the estimates of AGB were 116.0, 101.2, and 111.3 Mg ha⁻¹, respectively. Corresponding standard error estimates were 3.7, 1.6, and 3.2 Mg ha⁻¹. At the smallest geographical level (stand) an independent validation on 35 large field plots was carried out. RMSE values of 17.1-17.3 Mg ha⁻¹ and 42.6-53.2 Mg ha⁻¹ were found for LiDAR and InSAR, respectively. A time lag of six years between acquisition of InSAR data and field inventory has introduced some errors. Significant differences between estimates and reference values were found, illustrating the risk of using pure model-based methods in the estimation when there is a lack of fit in the models. We conclude that the examined remote sensing techniques can provide biomass estimates with smaller estimated errors than a field-based sample survey. The improvement can be highly significant, especially for LiDAR.     

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.     

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.     

Estimating forest canopy height and terrain relief from GLAS waveform metrics

Quantifying aboveground biomass in forest ecosystems is required for carbon stock estimation, aspects of forest management, and further developing a capacity for monitoring carbon stocks over time. Airborne Light Detection And Ranging (LiDAR) systems, of all remote sensing technologies, have been demonstrated to yield the most accurate estimates of aboveground biomass for forested areas over a wide range of biomass values. However, these systems are limited by considerations including large data volumes and high costs. Within the constraints imposed by the nature of the satellite mission, the GeoScience Laser Altimeter System (GLAS) aboard ICESat has provided data conferring information regarding forest vertical structure for large areas at a low end user cost. GLAS data have been demonstrated to accurately estimate forest height and aboveground biomass especially well in topographically smooth areas with homogeneous forested conditions. However in areas with dense forests, high relief, or heterogeneous vegetation cover, GLAS waveforms are more complex and difficult to consistently characterize. We use airborne discrete return LiDAR data to simulate GLAS waveforms and to subsequently deconstruct coregistered GLAS waveforms into vegetation and ground returns. A series of waveform metrics was calculated and compared to topography and vegetation information gleaned from the airborne data. A model to estimate maximum relief directly from waveform metrics was developed with an R² of 0.76 (n = 110), and used for the classification of the maximum relief of the areas sensed by GLAS. Discriminant analysis was also conducted as an alternative classification technique. A model was also developed estimating forest canopy height from waveform metrics for all of the data (R² = 0.81, n = 110) and for the three separate relief classes; maximum relief 0-7 m (R² = 0.83, n = 44), maximum relief 7-15 m (R² = 0.88, n = 41) and maximum relief > 15 m (R² = 0.75, n = 25). The moderate relief class model yielded better predictions of forest height than the low relief class model which is attributed to the increasing variability of waveform metrics with terrain relief. The moderate relief class model also yielded better predictions than the high relief class model because of the mixing of vegetation and terrain signals in waveforms from high relief footprints. This research demonstrates that terrain can be accurately modeled directly from GLAS waveforms enabling the inclusion of terrain relief, on a waveform specific basis, as supplemental model input to improve estimates of canopy height.     

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.     

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.     

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).     

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.     

Measuring forest structure and biomass in New England forest stands using Echidna ground-based lidar

A ground-based, upward-scanning, near-infrared lidar, the Echidna® validation instrument (EVI), built by CSIRO Australia, retrieves forest stand structural parameters, including mean diameter at breast height (DBH), stem count density (stems/area), basal area, and above-ground woody biomass with very good accuracy in six New England hardwood and conifer forest stands. Comparing forest structural parameters retrieved using EVI data with extensive ground measurements, we found excellent agreement at the site level using five EVI scans (plots) per site (R² = 0.94-0.99); very good agreement at the plot level for stem count density and biomass (R² = 0.90-0.85); and good agreement at the plot level for mean DBH and basal area (R² = 0.48-0.66). The observed variance at site and plot levels suggest that a sample area of at least 1 ha (10⁴ m²) is required to estimate these parameters accurately at the stand level using either lidar-based or conventional methods. The algorithms and procedures used to retrieve these structural parameters are dependent on the unique ability of the Echidna® lidar to digitize the full waveform of the scattered lidar pulse as it returns to the instrument, which allows consistent separation of scattering by trunks and large branches from scattering by leaves. This successful application of ground-based lidar technology opens the door to rapid and accurate measurement of biomass and timber volume in areal sampling scenarios and as a calibration and validation tool for mapping biomass using airborne or spaceborne remotely sensed data.     

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.     

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.     

The BIOMASS mission: Mapping global forest biomass to better understand the terrestrial carbon cycle

In response to the urgent need for improved mapping of global biomass and the lack of any current space systems capable of addressing this need, the BIOMASS mission was proposed to the European Space Agency for the third cycle of Earth Explorer Core missions and was selected for Feasibility Study (Phase A) in March 2009. The objectives of the mission are 1) to quantify the magnitude and distribution of forest biomass globally to improve resource assessment, carbon accounting and carbon models, and 2) to monitor and quantify changes in terrestrial forest biomass globally, on an annual basis or better, leading to improved estimates of terrestrial carbon sources (primarily from deforestation); and terrestrial carbon sinks due to forest regrowth and afforestation. These science objectives require the mission to measure above-ground forest biomass from 70° N to 56° S at spatial scale of 100-200 m, with error not exceeding ± 20% or ± 10 t ha^− 1 and forest height with error of ± 4 m. To meet the measurement requirements, the mission will carry a P-Band polarimetric SAR (centre frequency 435 MHz with 6 MHz bandwidth) with interferometric capability, operating in a dawn-dusk orbit with a constant incidence angle (in the range of 25°-35°) and a 25-45 day repeat cycle. During its 5-year lifetime, the mission will be capable of providing both direct measurements of biomass derived from intensity data and measurements of forest height derived from polarimetric interferometry. The design of the BIOMASS mission spins together two main observational strands: (1) the long heritage of airborne observations in tropical, temperate and boreal forest that have demonstrated the capabilities of P-band SAR for measuring forest biomass; (2) new developments in recovery of forest structure including forest height from Pol-InSAR, and, crucially, the resistance of P-band to temporal decorrelation, which makes this frequency uniquely suitable for biomass measurements with a single repeat-pass satellite. These two complementary measurement approaches are combined in the single BIOMASS sensor, and have the satisfying property that increasing biomass reduces the sensitivity of the former approach while increasing the sensitivity of the latter. This paper surveys the body of evidence built up over the last decade, from a wide range of airborne experiments, which illustrates the ability of such a sensor to provide the required measurements.At present, the BIOMASS P-band radar appears to be the only sensor capable of providing the necessary global knowledge about the world's forest biomass and its changes. In addition, this first chance to explore the Earth's environment with a long wavelength satellite SAR is expected to make yield new information in a range of geoscience areas, including subsurface structure in arid lands and polar ice, and forest inundation dynamics.     

The use of waveform lidar to measure northern temperate mixed conifer and deciduous forest structure in New Hampshire

The direct retrieval of canopy height and the estimation of aboveground biomass are two important measures of forest structure that can be quantified by airborne laser scanning at landscape scales. These and other metrics are central to studies attempting to quantify global carbon cycles and to improve understanding of the spatial variation in forest structure evident within differing biomes. Data acquired using NASA's Laser Vegetation Imaging Sensor (LVIS) over the Bartlett Experimental Forest (BEF) in central New Hampshire (USA) was used to assess the performance of waveform lidar in a northern temperate mixed conifer and deciduous forest.Using coincident plots established for this study, we found strong agreement between field and lidar measurements of height (r² = 0.80, p < 0.000) at the footprint level. Allometric calculations of aboveground biomass (AGBM) and LVIS metrics (AGBM: r² = 0.61, PRESS RMSE = 58.0 Mg ha^− 1, p < 0.000) and quadratic mean stem diameter (QMSD) and LVIS metrics (r² = 0.54, p = 0.002) also showed good agreement at the footprint level. Application of a generalized equation for determining AGBM proposed by Lefsky et al. (2002a) to footprint-level field data from Bartlett resulted in a coefficient of determination of 0.55; RMSE = 64.4 Mg ha^− 1; p = 0.002. This is slightly weaker than the strongest relationship found with the best-fit single term regression model.Relationships between a permanent grid of USDA Forest Service inventory plots and the mean values of aggregated LVIS metrics, however, were not as strong. This discrepancy suggests that validation efforts must be cautious in using pre-existing field data networks as a sole means of calibrating and verifying such remote sensing data. Stratification based on land-use or species composition, however, did provide the means to improve regression relationships at this scale. Regression models established at the footprint level for AGBM and QMSD were applied to LVIS data to generate predicted values for the whole of Bartlett. The accuracy of these models was assessed using varying subsets of the USFS NERS plot data. Coefficient of determinations ranged from fair to strong with aspects of land-use history and species composition influencing both the fit and the level of error seen in the predicted relationships.     

L- and P-band backscatter intensity for biomass retrieval in hemiboreal forest

At present, the greatest source of uncertainty in the global carbon cycle is in the terrestrial ecosystems. In order to reduce these uncertainties it is necessary to provide consistent and accurate global estimates of the world forest biomass. One of the most promising methods for obtaining such estimates is through polarimetric SAR backscatter measurements at low frequencies. In this paper, the relation between polarimetric SAR backscatter at L- and P-bands and forest biomass is investigated using data acquired within the BioSAR-I campaign in southern Sweden during 2007. Methods for estimating biomass on stand level using these data are developed and evaluated, and the results for the two frequency bands are compared. For L-band data, the best results were obtained using HV-polarized backscatter only, giving estimation errors in terms of root mean square errors (RMSE) between 31% and 46% of the mean biomass for stands with biomass ranging from 10 to 290 t/ha, and an (adjusted) coefficient of determination (R²) between 0.4 and 0.6. For P-band data, the results are better than for L-band. Models using HV- or HH-polarized P-band backscatter give similar results, as does a model including both HV and HH. The RMSEs were between 18 and 27%, and the R² values were between 0.7 and 0.8.     

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.