Forest vertical structure from GLAS: An evaluation using LVIS and SRTM data

The Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat) is the first spaceborne lidar instrument for continuous global observation of the Earth. GLAS records a vertical profile of the returned laser energy from its footprint. To help understand the application of the data for forest spatial structure studies in our regional projects, an evaluation of the GLAS data was conducted using NASA's Laser Vegetation Imaging Sensor (LVIS) data in an area near NASA's Goddard Space Flight Center in Greenbelt, Maryland, USA. The tree height indices from airborne large-footprint lidars such as LVIS have been successfully used for estimation of forest structural parameters in many previous studies and served as truth in this study.The location accuracy of the GLAS footprints was evaluated by matching the elevation profile from GLAS with the Shuttle Radar Topography Mission (SRTM) DEM. The results confirmed the location accuracy of the GLAS geolocation, and showed a high correlation between the height of the scattering phase center from SRTM and the top tree height from GLAS data. The comparisons between LVIS and GLAS data showed that the GLAS waveform is similar to the aggregation of the LVIS waveforms within the GLAS footprint, and the tree height indices derived from the GLAS and LVIS waveforms were highly correlated. The best correlations were found between the 75% waveform energy quartiles of LVIS and GLAS (r² = 0.82 for October 2003 GLAS data, and r² = 0.65 for June 2005 GLAS data). The correlations between the 50% waveform energy quartiles of LVIS and GLAS were also high (0.77 and 0.66 respectively). The comparisons of the top tree height and total length of waveform of the GLAS data acquired in fall of 2003 and early summer of 2005 showed a several meter bias. Because the GLAS footprints from these two orbits did not exactly overlap, several other factors may have caused this observed difference, including difference of forest structures, seasonal difference of canopy structures and errors in identifying the ground peak of waveforms.     

Predicting lidar measured forest vertical structure from multi-angle spectral data

A capability to remotely measure the vertical and spatial distribution of forest structure is required for more accurate modeling of energy, carbon, water, and climate over regional, continental, and global scales. We examined the potential of using a multi-angle spectral sensor to predict forest vertical structure as measured by an airborne lidar system. Data were acquired from AirMISR (Airborne Multi-Angle Imaging Spectrometer) and airborne LVIS (Laser Vegetation Imaging Sensor) for a 7000 ha study site near Howland Maine, consisting of small plantations, multi-generation clearings and large natural forest stands. The LVIS data set provided a relatively direct measure of forest vertical structure at a fine scale (20 m diameter footprints). Multivariate linear regression and neural network models were developed to predict the LVIS forest energy height measures from 28 AirMISR multi-angle spectral radiance values. The best model accurately predicted the maximum canopy height (as measured from LVIS) using AirMISR data (rmse = 0.92 m, R² = 0.89). The models developed in this study achieved high accuracies over a study site with an elaborate patchwork of forest communities with exceptional diversity in forest structure. We conclude that models using MISR-like data are capable of accurately predicting the vertical structure of forest canopies.     

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.     

Impact of footprint diameter and off-nadir pointing on the precision of canopy height estimates from spaceborne lidar

A spaceborne lidar mission could serve multiple scientific purposes including remote sensing of ecosystem structure, carbon storage, terrestrial topography and ice sheet monitoring. The measurement requirements of these different goals will require compromises in sensor design. Footprint diameters that would be larger than optimal for vegetation studies have been proposed. Some spaceborne lidar mission designs include the possibility that a lidar sensor would share a platform with another sensor, which might require off-nadir pointing at angles of up to 16°. To resolve multiple mission goals and sensor requirements, detailed knowledge of the sensitivity of sensor performance to these aspects of mission design is required.This research used a radiative transfer model to investigate the sensitivity of forest height estimates to footprint diameter, off-nadir pointing and their interaction over a range of forest canopy properties. An individual-based forest model was used to simulate stands of mixed conifer forest in the Tahoe National Forest (Northern California, USA) and stands of deciduous forests in the Bartlett Experimental Forest (New Hampshire, USA). Waveforms were simulated for stands generated by a forest succession model using footprint diameters of 20 m to 70 m. Off-nadir angles of 0 to 16° were considered for a 25 m diameter footprint diameter.Footprint diameters in the range of 25 m to 30 m were optimal for estimates of maximum forest height (R² of 0.95 and RMSE of 3 m). As expected, the contribution of vegetation height to the vertical extent of the waveform decreased with larger footprints, while the contribution of terrain slope increased. Precision of estimates decreased with an increasing off-nadir pointing angle, but off-nadir pointing had less impact on height estimates in deciduous forests than in coniferous forests. When pointing off-nadir, the decrease in precision was dependent on local incidence angle (the angle between the off-nadir beam and a line normal to the terrain surface) which is dependent on the off-nadir pointing angle, terrain slope, and the difference between the laser pointing azimuth and terrain aspect; the effect was larger when the sensor was aligned with the terrain azimuth but when aspect and azimuth are opposed, there was virtually no effect on R² or RMSE. A second effect of off-nadir pointing is that the laser beam will intersect individual crowns and the canopy as a whole from a different angle which had a distinct effect on the precision of lidar estimates of height, decreasing R² and increasing RMSE, although the effect was most pronounced for coniferous crowns.     

Combining lidar estimates of aboveground biomass and Landsat estimates of stand age for spatially extensive validation of modeled forest productivity

Extensive estimates of forest productivity are required to understand the relationships between shifting land use, changing climate and carbon storage and fluxes. Aboveground net primary production of wood (NPP_Aw) is a major component of total NPP and of net ecosystem production (NEP). Remote sensing of NPP and NPP_Aw is generally based on light use efficiency or process-based biogeochemistry models. However, validating these large area flux estimates remains a major challenge. In this study we develop an independent approach to estimating NPP_Aw, based on stand age and biomass, that could be implemented over a large area and used in validation efforts. Stand age is first mapped by iterative unsupervised classification of a multi-temporal sequence of images from a passive optical sensor (e.g. Landsat TM). Stand age is then cross-tabulated with estimates of stand height and aboveground biomass from lidar remote sensing. NPP_Aw is then calculated as the average increment in lidar-estimated biomass over the time period determined using change detection. In western Oregon, productivity estimates made using this method compared well with forest inventory estimates and were significantly different than estimates from a spatially distributed biogeochemistry model. The generality of the relationship between lidar-based canopy characteristics and stand biomass means that this approach could potentially be widely applicable to landscapes with stand replacing disturbance regimes, notably in regions where forest inventories are not routinely maintained.     

Validation of surface height from shuttle radar topography mission using shuttle laser altimeter

Spaceborne Interferometric SAR (InSAR) technology used in the Shuttle Radar Topography Mission (SRTM) and spaceborne lidar such as Shuttle Laser Altimeter-02 (SLA-02) are two promising technologies for providing global scale digital elevation models (DEMs). Each type of these systems has limitations that affect the accuracy or extent of coverage. These systems are complementary in developing DEM data. In this study, surface height measured independently by SRTM and SLA-02 was cross-validated. SLA data was first verified by field observations, and examinations of individual lidar waveforms. The geolocation accuracy of the SLA height data sets was examined by checking the correlation between the SLA surface height with SRTM height at 90 m resolution, while shifting the SLA ground track within its specified horizontal errors. It was found that the heights from the two instruments were highly correlated along the SLA ground track, and shifting the positions did not improve the correlation significantly. Absolute surface heights from SRTM and SLA referenced to the same horizontal and vertical datum (World Geodetic System (WGS) 84 Ellipsoid) were compared. The effects of forest cover and surface slope on the height difference were also examined. After removing the forest effect on SRTM height, the mean height difference with SLA-02 was near zero. It can be further inferred from the standard deviation of the height differences that the absolute accuracy of SRTM height at low vegetation area is better than the SRTM mission specifications (16 m). The SRTM height bias caused by forest cover needs to be further examined using future spaceborne lidar (e.g. GLAS) data.     

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

An airborne lidar sampling strategy to model forest canopy height from Quickbird imagery and GEOBIA

High-resolution digital canopy models derived from airborne lidar data have the ability to provide detailed information on the vertical structure of forests. However, compared to satellite data of similar spatial resolution and extent, the small footprint airborne lidar data required to produce such models remain expensive. In an effort to reduce these costs, the primary objective of this paper is to develop an airborne lidar sampling strategy to model full-scene forest canopy height from optical imagery, lidar transects and Geographic Object-Based Image Analysis (GEOBIA). To achieve this goal, this research focuses on (i) determining appropriate lidar transect features (i.e., location, direction and extent) from an optical scene, (ii) developing a mechanism to model forest canopy height for the full-scene based on a minimum number of lidar transects, and (iii) defining an optimal mean object size (MOS) to accurately model the canopy composition and height distribution. Results show that (i) the transect locations derived from our optimal lidar transect selection algorithm accurately capture the canopy height variability of the entire study area; (ii) our canopy height estimation models have similar performance in two lidar transect directions (i.e., north-south and west-east); (iii) a small lidar extent (17.6% of total size) can achieve similar canopy height estimation accuracies as those modeled from the full lidar scene; and (iv) different MOS can lead to distinctly different canopy height results. By comparing the best canopy height estimate with the full lidar canopy height data, we obtained average estimation errors of 6.0 m and 6.8 m for conifer and deciduous forests at the individual tree crown/small tree cluster level, and an area weighted combined error of 6.2 m, which is lower than the provincial forest inventory height class interval (i.e., ≈ 9.0 m).     

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.     

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.     

Measuring effective leaf area index, foliage profile, and stand height in New England forest stands using a full-waveform ground-based lidar

Effective leaf area index (LAI) retrievals from a scanning, ground-based, near-infrared (1064 nm) lidar that digitizes the full return waveform, the Echidna Validation Instrument (EVI), are in good agreement with those obtained from both hemispherical photography and the Li-Cor LAI-2000 Plant Canopy Analyzer. We conducted trials at 28 plots within six stands of hardwoods and conifers of varying height and stocking densities at Harvard Forest, Massachusetts, Bartlett Experimental Forest, New Hampshire, and Howland Experimental Forest, Maine, in July 2007. Effective LAI values retrieved by four methods, which ranged from 3.42 to 5.25 depending on the site and method, were not significantly different (β < 0.1 among four methods). The LAI values also matched published values well. Foliage profiles (leaf area with height) retrieved from the lidar scans, although not independently validated, were consistent with stand structure as observed and as measured by conventional methods. Canopy mean top height, as determined from the foliage profiles, deviated from mean RH100 values obtained from the Lidar Vegetation Imaging Sensor (LVIS) airborne large-footprint lidar system at 27 plots by − 0.91 m with RMSE = 2.04 m, documenting the ability of the EVI to retrieve stand height. The Echidna Validation Instrument is the first realization of the Echidna® lidar concept, devised by Australia's Commonwealth Scientific and Industrial Research Organization (CSIRO), for measuring forest structure using full-waveform, ground-based, scanning lidar.     

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.     

Integrating waveform lidar with hyperspectral imagery for inventory of a northern temperate forest

It has been suggested that attempts to use remote sensing to map the spatial and structural patterns of individual tree species abundances in heterogeneous forests, such as those found in northeastern North America, may benefit from the integration of hyperspectral or multi-spectral information with other active sensor data such as lidar. Towards this end, we describe the integrated and individual capabilities of waveform lidar and hyperspectral data to estimate three common forest measurements - basal area (BA), above-ground biomass (AGBM) and quadratic mean stem diameter (QMSD) - in a northern temperate mixed conifer and deciduous forest. The use of this data to discriminate distribution and abundance patterns of five common and often, dominant tree species was also explored. Waveform lidar imagery was acquired in July 2003 over the 1000 ha. Bartlett Experimental Forest (BEF) in central New Hampshire (USA) using NASA's airborne Laser Vegetation Imaging Sensor (LVIS). High spectral resolution imagery was likewise acquired in August 2003 using NASA's Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). Field data (2001-2003) from over 400 US Forest Service Northern Research Station (USFS NRS) plots were used to determine actual site conditions.Results suggest that the integrated data sets of hyperspectral and waveform lidar provide improved outcomes over use of either data set alone in evaluating common forest metrics. Across all forest conditions, 8-9% more of the variation in AGBM, BA, and QMSD was explained by use of the integrated sensor data in comparison to either AVIRIS or LVIS metrics applied singly, with estimated error 5-8% lower for these variables. Notably, in an analysis using integrated data limited to unmanaged forest tracts, AGBM coefficients of determination improved by 25% or more, while corresponding error levels decreased by over 25%. When data were restricted based on the presence of individual tree species within plots, AVIRIS data alone best predicted species-specific patterns of abundance as determined by species fraction of biomass. Nonetheless, use of LVIS and AVIRIS data - in tandem - produced complementary maps of estimated abundance and structure for individual tree species, providing a promising adjunct to traditional forest inventory and conservation biology planning efforts.     

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.     

A systematic method for 3D mapping of mangrove forests based on Shuttle Radar Topography Mission elevation data, ICEsat/GLAS waveforms and field data: Application to Ciénaga Grande de Santa Marta, Colombia

Mangrove forests are found within the intertropical zone and are one of the most biodiverse and productive wetlands on Earth. We focus on the Ciénaga Grande de Santa Marta (CGSM) in Colombia, the largest coastal lagoon–delta ecosystem in the Caribbean area with an extension of 1280 km², where one of the largest mangrove rehabilitation projects in Latin America is currently underway. Extensive man-made hydrological modifications in the region caused hypersaline soil (> 90 g kg^− 1) conditions since the 1960s triggering a large dieback of mangrove wetlands (~ 247 km²). In this paper, we describe a new systematic methodology to measure mangrove height and aboveground biomass by remote sensing. The method is based on SRTM (Shuttle Radar Topography Mission) elevation data, ICEsat/GLAS waveforms (Ice, Cloud, and Land Elevation Satellite/Geoscience Laser Altimeter System) and field data. Since the locations of the ICEsat and field datasets do not coincide, they are used independently to calibrate SRTM elevation and produce a map of mangrove canopy height. We compared height estimation methods based on waveform centroids and the canopy height profile (CHP). Linear relationships between ICEsat height estimates and SRTM elevation were derived. We found the centroid of the canopy waveform contribution (CWC) to be the best height estimator. The field data was used to estimate a SRTM canopy height bias (− 1.3 m) and estimation error (rms = 1.9 m). The relationship was applied to the SRTM elevation data to produce a mangrove canopy height map. Finally, we used field data and published allometric equations to derive an empirical relationship between canopy height and biomass. This relationship was used to scale the mangrove height map and estimate aboveground biomass distribution for the entire CGSM. The mean mangrove canopy height in CGSM is 7.7 m and most of the biomass is concentrated in forests around 9 m in height. Our biomass maps will enable estimation of regeneration rates of mangrove forests under hydrological rehabilitation at large spatial scales over the next decades. They will also be used to assess how highly disturbed mangrove forests respond to increasing sea level rise under current global climate change scenarios.     

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.     

Characterizing forest canopy structure with lidar composite metrics and machine learning

A lack of reliable observations for canopy science research is being partly overcome by the gradual use of lidar remote sensing. This study aims to improve lidar-based canopy characterization with airborne laser scanners through the combined use of lidar composite metrics and machine learning models. Our so-called composite metrics comprise a relatively large number of lidar predictors that tend to retain as much information as possible when reducing raw lidar point clouds into a format suitable as inputs to predictive models of canopy structural variables. The information-rich property of such composite metrics is further complemented by machine learning, which offers an array of supervised learning models capable of relating canopy characteristics to high-dimensional lidar metrics via complex, potentially nonlinear functional relationships. Using coincident lidar and field data over an Eastern Texas forest in USA, we conducted a case study to demonstrate the ubiquitous power of the lidar composite metrics in predicting multiple forest attributes and also illustrated the use of two kernel machines, namely, support vector machine and Gaussian processes (GP). Results show that the two machine learning models in conjunction with the lidar composite metrics outperformed traditional approaches such as the maximum likelihood classifier and linear regression models. For example, the five-fold cross validation for GP regression models (vs. linear/log-linear models) yielded a root mean squared error of 1.06 (2.36) m for Lorey's height, 0.95 (3.43) m for dominant height, 5.34 (8.51) m²/ha for basal area, 21.4 (40.5) Mg/ha for aboveground biomass, 6.54 (9.88) Mg/ha for belowground biomass, 0.75 (2.76) m for canopy base height, 2.2 (2.76) m for canopy ceiling height, 0.015 (0.02) kg/m³ for canopy bulk density, 0.068 (0.133) kg/m² for available canopy fuel, and 0.33 (0.39) m²/m² for leaf area index. Moreover, uncertainty estimates from the GP regression were more indicative of the true errors in the predicted canopy variables than those from their linear counterparts. With the ever-increasing accessibility of multisource remote sensing data, we envision a concomitant expansion in the use of advanced statistical methods, such as machine learning, to explore the potentially complex relationships between canopy characteristics and remotely-sensed predictors, accompanied by a desideratum for improved error analysis.     

Forest carbon densities and uncertainties from Lidar, QuickBird, and field measurements in California

Greenhouse gas inventories and emissions reduction programs require robust methods to quantify carbon sequestration in forests. We compare forest carbon estimates from Light Detection and Ranging (Lidar) data and QuickBird high-resolution satellite images, calibrated and validated by field measurements of individual trees. We conducted the tests at two sites in California: (1) 59 km² of secondary and old-growth coast redwood (Sequoia sempervirens) forest (Garcia-Mailliard area) and (2) 58 km² of old-growth Sierra Nevada forest (North Yuba area). Regression of aboveground live tree carbon density, calculated from field measurements, against Lidar height metrics and against QuickBird-derived tree crown diameter generated equations of carbon density as a function of the remote sensing parameters. Employing Monte Carlo methods, we quantified uncertainties of forest carbon estimates from uncertainties in field measurements, remote sensing accuracy, biomass regression equations, and spatial autocorrelation. Validation of QuickBird crown diameters against field measurements of the same trees showed significant correlation (r = 0.82, P < 0.05). Comparison of stand-level Lidar height metrics with field-derived Lorey's mean height showed significant correlation (Garcia-Mailliard r = 0.94, P < 0.0001; North Yuba R = 0.89, P < 0.0001). Field measurements of five aboveground carbon pools (live trees, dead trees, shrubs, coarse woody debris, and litter) yielded aboveground carbon densities (mean ± standard error without Monte Carlo) as high as 320 ± 35 Mg ha^− 1 (old-growth coast redwood) and 510 ± 120 Mg ha^− 1 (red fir [Abies magnifica] forest), as great or greater than tropical rainforest. Lidar and QuickBird detected aboveground carbon in live trees, 70-97% of the total. Large sample sizes in the Monte Carlo analyses of remote sensing data generated low estimates of uncertainty. Lidar showed lower uncertainty and higher accuracy than QuickBird, due to high correlation of biomass to height and undercounting of trees by the crown detection algorithm. Lidar achieved uncertainties of < 1%, providing estimates of aboveground live tree carbon density (mean ± 95% confidence interval with Monte Carlo) of 82 ± 0.7 Mg ha^− 1 in Garcia-Mailliard and 140 ± 0.9 Mg ha^− 1 in North Yuba. The method that we tested, combining field measurements, Lidar, and Monte Carlo, can produce robust wall-to-wall spatial data on forest carbon.     

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.     

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.