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.     

Exploring LiDAR-RaDAR synergy—predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR

Scanning Light Detecting and Ranging (LiDAR), Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) were analyzed to determine (1) which of the three sensor systems most accurately predicted forest biomass, and (2) if LiDAR and SAR/InSAR data sets, jointly considered, produced more accurate, precise results relative to those same data sets considered separately. LiDAR ranging measurements, VHF-SAR cross-sectional returns, and X- and P-band cross-sectional returns and interferometric ranges were regressed with ground-estimated (from dbh) forest biomass in ponderosa pine forests in the southwestern United States. All models were cross-validated. Results indicated that the average canopy height measured by the scanning LiDAR produced the best predictive equation. The simple linear LiDAR equation explained 83% of the biomass variability (n = 52 plots) with a cross-validated root mean square error of 26.0 t/ha. Additional LiDAR metrics were not significant to the model. The GeoSAR P-band (λ = 86 cm) cross-sectional return and the GeoSAR/InSAR canopy height (X-P) captured 30% of the forest biomass variation with an average predictive error of 52.5 t/ha. A second RaDAR-FOPEN collected VHF (λ ∼ 7.8 m) and cross-polarized P-band (λ = 88 cm) cross-sectional returns, none of which proved useful for forest biomass estimation (cross-validated R² = 0.09, RMSE = 63.7 t/ha). Joint consideration of LiDAR and RaDAR measurements produced a statistically significant, albeit small improvement in biomass estimation precision. The cross-validated R² increased from 83% to 84% and the prediction error decreased from 26.0 t/ha to 24.9 t/ha when the GeoSAR X-P interferometric height is considered along with the average LiDAR canopy height. Inclusion of a third LiDAR metric, the 60th decile height, further increased the R² to 85% and decreased the RMSE to 24.1 t/ha. On this 11 km² ponderosa pine study area, LiDAR data proved most useful for predicting forest biomass. RaDAR ranging measurements did not improve the LiDAR estimates.     

Empirical relationships between AIRSAR backscatter and LiDAR-derived forest biomass, Queensland, Australia

To evaluate the use of multi-frequency, polarimetric Synthetic Aperture Radar (SAR) data for quantifying the above ground biomass (AGB) of open forests and woodlands, NASA JPL AIRSAR (POLSAR) data were acquired over a 37 × 60 km area west of Injune, central Queensland, Australia. From field measurements recorded within 32 50 × 50 m plots, AGB was estimated by applying species-specific allometric equations to stand measurements. AGB was then scaled-up to the larger area using relationships established with Light Detection and Ranging (LiDAR) data acquired over 150 (10 columns, 15 rows) 500 × 150 m cells (or Primary Sampling Units, PSUs) spaced 4 × 4 km apart in the north- and east-west directions. Large-scale (1 : 4000) stereo aerial photographs were also acquired for each PSU to assess species composition. Based on the LiDAR extrapolations, the median AGB for the PSU grid was 82 Mg ha^− 1 (maximum 164 Mg ha^− 1), with the higher levels associated with forests containing a high proportion of Angophora and Callitris species. Empirical relationships between AGB and SAR backscatter confirmed that C-, L- and P-band saturated at different levels and revealed a greater strength in the relationship at higher incidence angles and a larger dynamic range and consistency of relationships at HV polarizations. A higher level of saturation (above ∼50 Mg ha^− 1) was observed at C-band HV compared to that reported for closed forests which was attributable to a link between foliage projected cover (FPC) and AGB. The study concludes that L-band HV backscatter data acquired at incidence angles approaching or exceeding 45° are best suited for estimating the AGB up to the saturation level of ∼80-85 Mg ha^− 1. For regional mapping of biomass below the level of saturation, the use of the Japanese Space Exploration Agency (JAXA) Advanced Land Observing Satellite (ALOS) Phase Arrayed L-band SAR (PALSAR) is advocated.     

Use of Large-Footprint Scanning Airborne Lidar To Estimate Forest Stand Characteristics in the Western Cascades of Oregon

A scanning lidar, a relatively new type of sensor which explicitly measures canopy height, was used to measure structure of conifer forests in the Pacific Northwest. SLICER (Scanning Lidar Imager of Canopies by Echo Recovery), an airborne pulsed laser developed by NASA which scans a swath of five 10-m diameter footprints along the aircraft’s flightpath, captures the power of the reflected laser pulse as a function of height from the top of the canopy to the ground. Ground measurements of forest stand structure were collected on 26 plots with coincident SLICER data. Height, basal area, total biomass, and leaf biomass as estimated from field data could be predicted from SLICER-derived metrics with r² values of 0.95, 0.96, 0.96, and 0.84, respectively. These relationships were strong up to a height of 52 m, basal area of 132 m²/ha and total biomass of 1300 Mg/ha. In light of these strong relationships, large-footprint, airborne scanning lidar shows promise for characterizing stand structure for management and research purposes.     

Height and biomass of mangroves in Africa from ICESat/GLAS and SRTM

The accurate quantification of the three-dimensional (3-D) structure of mangrove forests is of great importance, particularly in Africa where deforestation rates are high and the lack of background data is a major problem. The objectives of this study are to estimate (1) the total area, (2) canopy height distributions, and (3) above-ground biomass (AGB) of mangrove forests in Africa. To derive the 3-D structure and biomass maps of mangroves, we used a combination of mangrove maps derived from Landsat Enhanced Thematic Mapper Plus (ETM+), lidar canopy height estimates from ICESat/GLAS (Ice, Cloud, and land Elevation Satellite/Geoscience Laser Altimeter System), and elevation data from SRTM (Shuttle Radar Topography Mission) for the African continent. The lidar measurements from the large footprint GLAS sensor were used to derive local estimates of canopy height and calibrate the interferometric synthetic aperture radar (InSAR) data from SRTM. We then applied allometric equations relating canopy height to biomass in order to estimate AGB from the canopy height product. The total mangrove area of Africa was estimated to be 25,960 km² with 83% accuracy. The largest mangrove areas and the greatest total biomass were found in Nigeria covering 8573 km² with 132 × 10⁶ Mg AGB. Canopy height across Africa was estimated with an overall root mean square error of 3.55 m. This error includes the impact of using sensors with different resolutions and geolocation error. This study provides the first systematic estimates of mangrove area, height, and biomass in Africa.     

Measuring biomass changes due to woody encroachment and deforestation/degradation in a forest-savanna boundary region of central Africa using multi-temporal L-band radar backscatter

Satellite L-band synthetic aperture radar backscatter data from 1996 and 2007 (from JERS-1 and ALOS PALSAR respectively), were used with field data collected in 2007 and a back-calibration method to produce biomass maps of a 15 000 km² forest-savanna ecotone region of central Cameroon. The relationship between the radar backscatter and aboveground biomass (AGB) was strong (r² = 0.86 for ALOS HV to biomass plots, r² = 0.95 relating ALOS-derived biomass for 40 suspected unchanged regions to JERS-1 HH). The root mean square error (RMSE) associated with AGB estimation varied from ~ 25% for AGB < 100 Mg ha^− 1 to ~ 40% for AGB > 100 Mg ha^− 1 for the ALOS HV data. Change detection showed a significant loss of AGB over high biomass forests, due to suspected deforestation and degradation, and significant biomass gains along the forest-savanna boundary, particularly in areas of low population density. Analysis of the errors involved showed that radar data can detect changes in broad AGB class in forest-savanna transition areas with an accuracy > 95%. However, quantitative assessment of changes in AGB in Mg ha^− 1 at a pixel level will require radar images from sensors with similar characteristics collecting data from the same season over multiple years.     

Regional aboveground live carbon losses due to drought-induced tree dieback in piñon-juniper ecosystems

Recent large-scale dieback of piñon-juniper (P-J) woodlands and forests across the western US occurred as a result of multi-year drought and subsequent insect and disease outbreaks. P-J vegetation is spatially extensive, thus large-scale mortality events such as the one that has occurred over the past several years could significantly alter regional carbon (C) budgets. Our objective was to use a remote sensing technique coupled with field-based data to estimate changes in aboveground live C stocks across a 4100 km² region of Colorado caused by P-J tree mortality. We hypothesized that dieback would amplify the phenological dynamics of P-J vegetation, and these variations would be related to drought-induced losses of live P-J aboveground biomass (AGB) that are discernible using time-series remote sensing vegetation data. Here, we assess live P-J AGB loss using dry season fractional photosynthetic vegetation cover (PV) derived from multi-year Landsat images. Our results showed a strong linear positive relationship between the maximum decline in PV and field-measured losses of live P-J AGB during the period 2000-05 (r² = 0.64, p = 0.002). These results were then used to map AGB losses throughout the study region. Mean live aboveground C loss (± sd) was 10.0 (± 3.4) Mg C ha^− 1. Total aboveground live P-J C loss was 4.6 Tg C, which was approximately 39 times higher than the concurrent C loss attributed to wildfire and management treatments within or near to the national forests of the study region. Our results suggest that spatially extensive mortality events such as the one observed in P-J woodlands across the western US in the past decade may significantly alter the ecosystem C balance for decades to come. Remote sensing techniques to monitor changes in aboveground C stocks, such as the one developed in our study, may support regional and global C monitoring in the future.     

Mapping biomass and stress in the Sierra Nevada using lidar and hyperspectral data fusion

In this paper, we explored fusion of structural metrics from the Laser Vegetation Imaging Sensor (LVIS) and spectral characteristics from the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) for biomass estimation in the Sierra Nevada. In addition, we combined the two sensors to map species-specific biomass and stress at landscape scale. Multiple endmember spectral mixture analysis (MESMA) was used to classify vegetation from AVIRIS images and obtain sub-pixel fractions of green vegetation, non-photosynthetic vegetation, soil, and shade. LVIS metrics, AVIRIS spectral indices, and MESMA fractions were compared with field measures of biomass using linear and stepwise regressions at stand (1 ha) level. AVIRIS metrics such as water band indices and shade fractions showed strong correlation with LVIS canopy height (r² = 0.69, RMSE = 5.2 m) and explained around 60% variability in biomass. LVIS variables were found to be consistently good predictors of total and species specific biomass (r² = 0.77, RMSE = 70.12 Mg/ha). Prediction by LVIS after species stratification of field data reduced errors by 12% (r² = 0.84, RMSE = 58.78 Mg/ha) over using LVIS metrics alone. Species-specific biomass maps and associated errors created from fusion were different from those produced without fusion, particularly for hardwoods and pines, although mean biomass differences between the two techniques were not statistically significant. A combined analysis of spatial maps from LVIS and AVIRIS showed increased water and chlorophyll stress in several high biomass stands in the study area. This study provides further evidence that lidar is better suited for biomass estimation, per se, while the best use of hyperspectral data may be to refine biomass predictions through a priori species stratification, while also providing information on canopy state, such as stress. Together, the two sensors have many potential applications in carbon dynamics, ecological and habitat studies.     

Geographic variability in lidar predictions of forest stand structure in the Pacific Northwest

Estimation of the amount of carbon stored in forests is a key challenge for understanding the global carbon cycle, one which remote sensing is expected to help address. However, carbon storage in moderate to high biomass forests is difficult to estimate with conventional optical or radar sensors. Lidar (light detection and ranging) instruments measure the vertical structure of forests and thus hold great promise for remotely sensing the quantity and spatial organization of forest biomass. In this study, we compare the relationships between lidar-measured canopy structure and coincident field measurements of forest stand structure at five locations in the Pacific Northwest of the U.S.A. with contrasting composition. Coefficient of determination values (r²) ranged between 41% and 96%. Correlations for two important variables, LAI (81%) and aboveground biomass (92%), were noteworthy, as was the fact that neither variable showed an asymptotic response.Of the 17 stand structure variables considered in this study, we were able to develop eight equations that were valid for all sites, including equations for two variables generally considered to be highly important (aboveground biomass and leaf area index). The other six equations that were valid for all sites were either related to height (which is most directly measured by lidar) or diameter at breast height (which should be closely related to height). Four additional equations (a total of 12) were applicable to all sites where either Douglas-fir (Pseudotsuga menziesii), western hemlock (Tsuga heterophylla) or Sitka spruce (Picea sitchensi) were dominant. Stand structure variables in sites dominated by true firs (Abies sp.) or ponderosa pine (Pinus ponderosa) had biases when predicted by these four additional equations. Productivity-related variables describing the edaphic, climatic and topographic environment of the sites where available for every regression, but only two of the 17 equations (maximum diameter at breast height, stem density) incorporated them. Given the wide range of these environmental conditions sampled, we conclude that the prediction of stand structure is largely independent of environmental conditions in this study area.Most studies of lidar remote sensing for predicting stand structure have depended on intensive data collections within a relatively small study area. This study indicates that the relationships between many stand structure indices and lidar measured canopy structure have generality at the regional scale. This finding, if replicated in other regions, would suggest that mapping of stand structure using lidar may be accomplished by distributing field sites extensively over a region, thus reducing the overall inventory effort required.     

Generation of crown bulk density for Pinus sylvestris L. from lidar

The use of lidar data to estimate critical variables needed for modeling wildfire behavior was tested on a Scots pine forest (Pinus sylvestris L.) in central Spain. Lidar data accurately estimated crown bulk density at the plot level (r²=0.80). Lidar data could be used to directly estimate crown volume (r²=0.92) and foliage biomass (r²=0.84), which together produced better results than directly fitting the lidar data to crown bulk density. Incorporating equations that relate tree diameter at breast height and other forest parameters improved estimates of foliage biomass. Individual tree level analyses were not completely successful due to difficulty in accurately assigning laser pulses to the correct tree (r²=0.14).     

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.     

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.     

Aboveground forest biomass derived using multiple dates of WorldView-2 stereo-imagery: quantifying the improvement in estimation accuracy

ABSTRACT

The aim of this study was to investigate the capabilities of two date satellite-derived image-based point clouds (IPCs) to estimate forest aboveground biomass (AGB). The data sets used include panchromatic WorldView-2 stereo-imagery with 0.46 m spatial resolution representing 2014 and 2016 and a detailed digital elevation model derived from airborne laser scanning data. Altogether, 332 field sample plots with an area of 256 m² were used for model development and validation. Predictors describing forest height, density, and variation in height were extracted from the IPC 2014 and 2016 and used in k-nearest neighbour imputation models developed with sample plot data for predicting AGB. AGB predictions for 2014 (AGB₂₀₁₄) were projected to 2016 using growth models (AGB_{Projected_2016}) and combined with the AGB estimates derived from the 2016 data (AGB₂₀₁₆). AGB prediction model developed with 2014 data was also applied to 2016 data (AGB_{2016_pred2014}). Based on our results, the change in the 90^th percentile of height derived from the WorldView-2 IPC was able to characterize forest height growth between 2014 and 2016 with an average growth of 0.9 m. Features describing canopy cover and variation in height derived from the IPC were not as consistent. The AGB₂₀₁₆ had a bias of ?7.5% (?10.6 Mg ha^?1) and root mean square error (RMSE) of 26.0% (36.7 Mg ha^?1) as the respective values for AGB_{Projected_2016} were 7.0% (9.9 Mg ha^?1) and 21.5% (30.8 Mg ha^?1). AGB_{2016_pred2014} had a bias of ?19.6% (?27.7 Mg ha^?1) and RMSE of 33.2% (46.9 Mg ha^?1). By combining predictions of AGB₂₀₁₆ and AGB_{Projected_2016} at sample plot level as a weighted average, we were able to decrease the bias notably compared to estimates made on any single date. The lowest bias of ?0.25% (?0.4 Mg ha^?1) was obtained when equal weights of 0.5 were given to the AGB_{Projected_2016} and AGB₂₀₁₆ estimates. Respectively, RMSE of 20.9% (29.5 Mg ha^?1) was obtained using equal weights. Thus, we conclude that combination of two date WorldView-2 stereo-imagery improved the reliability of AGB estimates on sample plots where forest growth was the only change between the two dates.     

Distinguishing between live and dead standing tree biomass on the North Rim of Grand Canyon National Park, USA using small-footprint lidar data

Accurate estimation of live and dead biomass in forested ecosystems is important for studies of carbon dynamics, biodiversity, wildfire behavior, and for forest management. Lidar remote sensing has been used successfully to estimate live biomass, but studies focusing on dead biomass are rare. We used lidar data, in conjunction with field measurements from 58 plots to distinguish between and map standing live and dead tree biomass in the mixed coniferous forest of the North Rim of Grand Canyon National Park, USA. Lidar intensity and canopy volume were key variables for estimating live biomass, whereas for dead biomass, lidar intensity alone was critical for accurate estimation. Regression estimates of both live and dead biomass ranged between 0 and 600 Mg ha^− 1, with means of 195.08 Mg ha^− 1 and 65.73 Mg ha^− 1, respectively. Cross validation with field data resulted in correlation coefficients for predicted vs. observed of 0.85 for live biomass (RMSE = 50 Mg ha^− 1 and %RMSE (RMSE as a percent of the mean) = 26). For dead biomass, correlation was 0.79, RMSE was 42 Mg ha^− 1, and %RMSE was 63. Biomass maps revealed interesting patterns of live and dead standing tree biomass. Live biomass was highest in the ponderosa pine zone, and decreased from south to north through the mixed conifer and spruce-fir forest zones. Dead biomass exhibited a background range of values in these mature forests from zero to 100 Mg ha^− 1, with lower values in locations having higher live biomass. In areas with high dead biomass values, live biomass was near zero. These areas were associated with recent wildfires, as indicated by fire maps derived from the Monitoring Trends in Burn Severity Project (MTBS). Combining our dead biomass maps with the MTBS maps, we demonstrated the complementary power of these two datasets, revealing that MTBS burn intensity class can be described quantitatively in terms of dead biomass. Assuming a background range of dead biomass up to 100 Mg ha^− 1, it is possible to estimate and map the contribution to the standing dead tree biomass pool associated with recent wildfire.     

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.     

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.     

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

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

Estimating aboveground biomass of a mangrove plantation on the Northern coast of Vietnam using machine learning techniques with an integration of ALOS-2 PALSAR-2 and Sentinel-2A data

ABSTRACT

Aboveground biomass (AGB) of mangrove forest plays a crucial role in global carbon cycle by reducing greenhouse gas emissions and mitigating climate change impacts. Monitoring mangrove forests biomass accurately still remains challenging compared to other forest ecosystems. We investigated the usability of machine learning techniques for the estimation of AGB of mangrove plantation at a coastal area of Hai Phong city (Vietnam). The study employed a GIS database and support vector regression (SVR) to build and verify a model of AGB, drawing upon data from a survey in 25 sampling plots and an integration of Advanced Land Observing Satellite-2 Phased Array Type L-band Synthetic Aperture Radar-2 (ALOS-2 PALSAR-2) dual-polarization horizontal transmitting and horizontal receiving (HH) and horizontal transmitting and vertical receiving (HV) and Sentinel-2A multispectral data. The performance of the model was assessed using root mean square error (RMSE), mean absolute error (MAE), coefficient of determination (R²), and leave-one-out cross-validation. Usability of the SVR model was assessed by comparing with four state-of-the-art machine learning techniques, i.e. radial basis function neural networks, multi-layer perceptron neural networks, Gaussian process, and random forest. The SVR model shows a satisfactory result (R² = 0.596, RMSE = 0.187, MAE = 0.123) and outperforms the four machine learning models. The SVR model-estimated AGB ranged between 36.22 and 230.14 Mg ha^?1 (average = 87.67 Mg ha^?1). We conclude that an integration of ALOS-2 PALSAR-2 and Sentinel-2A data used with SVR model can improve the AGB accuracy estimation of mangrove plantations in tropical areas.     

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.     

Regional aboveground forest biomass using airborne and spaceborne LiDAR in Québec

Aboveground dry biomass was estimated for the 1.3 M km² forested area south of the treeline in the eastern Canadian province of Québec by combining data from an airborne and spaceborne LiDAR, a Landsat ETM+ land cover map, a Shuttle Radar Topographic Mission (SRTM) digital elevation model, ground inventory plots, and vegetation zone maps. Plot-level biomass was calculated using allometric relationships between tree attributes and biomass. A small footprint portable laser profiler then flew over these inventory plots to develop a generic airborne LiDAR-based biomass equation (R² = 0.65, n = 207). The same airborne LiDAR system flew along four portions of orbits of the ICESat Geoscience Laser Altimeter System (GLAS). A square-root transformed equation was developed to predict airborne profiling LiDAR estimates of aboveground dry biomass from GLAS waveform parameters combined with an SRTM slope index (R² = 0.59, n = 1325).Using the 104,044 quality-filtered GLAS pulses obtained during autumn 2003 from 97 orbits over the study area, we then predicted aboveground dry biomass for the main vegetation areas of Québec as well as for the entire Province south of the treeline. Including cover type covariances both within and between GLAS orbits increased standard errors of the estimates by two to five times at the vegetation zone level and as much as threefold at the provincial level. Aboveground biomass for the whole study area averaged 39.0 ± 2.2 (standard error) Mg ha^? 1 and totalled 4.9 ± 0.3 Pg. Biomass distributions were 12.6% northern hardwoods, 12.6% northern mixedwood, 38.4% commercial boreal, 13% non-commercial boreal, 14.2% taiga, and 9.2% treed tundra. Non-commercial forests represented 36% of the estimated aboveground biomass, thus highlighting the importance of remote northern forests to C sequestration. This study has shown that space-based forest inventories of northern forests could be an efficient way of estimating the amount, distribution, and uncertainty of aboveground biomass and carbon stocks at large spatial scales.