Geoscience Laser Altimeter System Footprints Research Articles

New Metrics and the Combinations for Estimating Forest Biomass From GLAS Data

Geoscience laser altimeter system (GLAS) data have been widely used for forest aboveground biomass (AGB) estimation, but there is no consensus on the optimal metrics. To explore whether a few optimal GLAS metrics could generate accurate AGB estimates, we proposed five metrics and explored their combinations with ten existing ones. The importance of these metrics was measured according to their contributions to changes in the cross-validated mean-squared error. The two to eight most important metrics were then selected to develop AGB models, and their performances were evaluated using field AGB. The optimal combination of GLAS metrics was finally used for AGB estimation at GLAS footprints from 2004 to 2007 within a 2°×2° spatial extent in Tahe and Changbai Mountain, China. The results showed that four GLAS metrics, including our proposed metric CRH25 (25th percentile of canopy reflection heights) combined with Lead, quadratic mean canopy height, and H75, yield the best AGB estimates, with an R 2 of 0.61±0.15 and RMSE of 52.20±23.50 Mg/ha, and the inclusion of more GLAS metrics did not improve the results. The estimated forest AGB in Tahe was 89.03±19.16 Mg/ha and 103.07±23.42 Mg/ha in Changbai Mountain. In both regions, the annual average forest AGB estimates for 2005 were higher than the AGB estimates for 2004, 2006, and 2007. The results of this study suggested that a few waveform parameters could enable the accurate estimation of forest AGB. Moreover, this study indicated that GLAS data might be able to monitor forest AGB changes, which require further investigation.

Retrieval of Vertical Foliage Profile and Leaf Area Index Using Transmitted Energy Information Derived from ICESat GLAS Data

The vertical foliage profile (VFP) and leaf area index (LAI) are critical descriptors in terrestrial ecosystem modeling. Although light detection and ranging (lidar) observations have been proven to have potential for deriving the VFP and LAI, existing methods depend only on the received waveform information and are sensitive to additional input parameters, such as the ratio of canopy to ground reflectance. In this study, we proposed a new method for retrieving forest VFP and LAI from Ice, Cloud and land Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS) data over two sites similar in their biophysical parameters. Our method utilized the information from not only the interaction between the laser and the forest but also the sensor configuration, which brought the benefit that our method was free from an empirical input parameter. Specifically, we first derived the transmitted energy profile (TEP) through the lidar 1-D radiative transfer model. Then, the obtained TEP was utilized to calculate the vertical gap distribution. Finally, the vertical gap distribution was taken as input to derive the VFP based on the Beer–Lambert law, and the LAI was calculated by integrating the VFP. Extensive validations of our method were carried out based on the discrete anisotropic radiative transfer (DART) simulation data, ground-based measurements, and the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI product. The validation based on the DART simulation data showed that our method could effectively characterize the VFP and LAI under various canopy architecture scenarios, including homogeneous turbid and discrete individual-tree scenes. The ground-based validation also proved the feasibility of our method: the VFP retrieved from the GLAS data showed a similar trend with the foliage distribution density in the GLAS footprints; the GLAS LAI had a high correlation with the field measurements, with a determination coefficient (R2) of 0.79, root mean square error (RMSE) of 0.49, and bias of 0.17. Once the outliers caused by low data quality and large slope were identified and removed, the accuracy was further improved, with R2 = 0.85, RMSE = 0.35, and bias = 0.10. However, the MODIS LAI product did not present a good relationship with the GLAS LAI. Relative to the GLAS LAI, the MODIS LAI showed an overestimation in the low and middle ranges of the LAI and a saturation at high values of approximately LAI = 5.5. Overall, this method has the potential to produce continental- and global-scale VFP and LAI datasets from the spaceborne lidar system.

Geometric Accuracy Improvement Method for High-Resolution Optical Satellite Remote Sensing Imagery Combining Multi-Temporal SAR Imagery and GLAS Data

With the widespread availability of satellite data, a single region can be described using multi-source and multi-temporal remote sensing data, such as high-resolution (HR) optical imagery, synthetic aperture radar (SAR) imagery, and space-borne laser altimetry data. These have become the main source of data for geopositioning. However, due to the limitation of the direct geometric accuracy of HR optical imagery and the effect of the small intersection angle of HR optical imagery in stereo pair orientation, the geometric accuracy of HR optical imagery cannot meet the requirements for geopositioning without ground control points (GCPs), especially in uninhabited areas, such as forests, plateaus, or deserts. Without satellite attitude error, SAR usually provides higher geometric accuracy than optical satellites. Space-borne laser altimetry technology can collect global laser footprints with high altitude accuracy. Therefore, this paper presents a geometric accuracy improvement method for HR optical satellite remote sensing imagery combining multi-temporal SAR Imagery and GLAS data without GCPs. Based on the imaging mechanism, the differences in the weight matrix determination of the HR optical imagery and SAR imagery were analyzed. The laser altimetry data with high altitude accuracy were selected and applied as height control point in combined geopositioning. To validate the combined geopositioning approach, GaoFen2 (GF2) optical imagery, GaoFen6 (GF6) optical imagery, GaoFen3 (GF3) SAR imagery, and the Geoscience Laser Altimeter System (GLAS) footprint were tested. The experimental results show that the proposed model can be effectively applied to combined geopositioning to improve the geometric accuracy of HR optical imagery. Moreover, we found that the distribution and weight matrix determination of SAR images and the distribution of GLAS footprints are the crucial factors influencing geometric accuracy. Combined geopositioning using multi-source remote sensing data can achieve a plane accuracy of 1.587 m and an altitude accuracy of 1.985 m, which is similar to the geometric accuracy of geopositioning of GF2 with GCPs.

Canopy Height Layering Biomass Estimation Model (CHL-BEM) with Full-Waveform LiDAR

Forest biomass is an important descriptor for studying carbon storage, carbon cycles, and global change science. The full-waveform spaceborne Light Detection And Ranging (LiDAR) Geoscience Laser Altimeter System (GLAS) provides great possibilities for large-scale and long-term biomass estimation. To the best of our knowledge, most of the existing research has utilized average tree height (or height metrics) within a GLAS footprint as the key parameter for biomass estimation. However, the vertical distribution of tree height is usually not as homogeneous as we would expect within such a large footprint of more than 2000 m2, which would limit the biomass estimation accuracy vastly. Therefore, we aim to develop a novel canopy height layering biomass estimation model (CHL-BEM) with GLAS data in this study. First, all the trees with similar height were regarded as one canopy layer within each GLAS footprint. Second, the canopy height and canopy cover of each layer were derived from GLAS waveform parameters. These parameters were extracted using a waveform decomposition algorithm (refined Levenberg–Marquardt—RLM), which assumed that each decomposed vegetation signal corresponded to a particular canopy height layer. Third, the biomass estimation model (CHL-BEM) was established by using the canopy height and canopy cover of each height layer. Finally, the CHL-BEM was compared with two typical biomass estimation models of GLAS in the study site located in Ejina, China, where the dominant species was Populus euphratica. The results showed that the CHL-BEM presented good agreement with the field measurement biomass (R2 = 0.741, RMSE = 0.487, %RMSE = 24.192) and achieved a significantly higher accuracy than the other two models. As a whole, we expect our method to advance all the full-waveform LiDAR development and applications, e.g., the newly launched Global Ecosystem Dynamics Investigation (GEDI).

Mapping of aboveground biomass in Gabon

The aim of this paper is to map the aboveground biomass (AGB) in Gabon. First, a random forest (RF) model that relates reference AGB values to remote sensing (RS)-derived variables (mainly radar and optical images) was built, and the significant predictive variables were determined. Second, the built RF model was applied to the significant RS-derived variables to predict AGB across Gabon. The results showed that the overall RMSE (Root Mean Square Error) on the RS-derived AGB map with a spatial resolution of 50 m was 63.3 t/ha (R2 = 0.53).To improve the accuracy of the RS-derived AGB map, the integration of LiDAR data provided by the Geoscience Laser Altimeter System (GLAS) onboard the Ice Cloud and Land Elevation Satellite (ICESat) was investigated. First, an RF model that relates reference AGB values to GLAS-derived metrics and a DEM (Digital Elevation Model) was built. Second, the calibrated RF model was applied to obtain a spatially distributed estimation of AGB (GLAS footprints geolocation) covering forested areas in Gabon, with a density of 0.13 footprints/km2. Third, the semivariogram of residuals (RS-derived AGB map – GLAS-derived AGB “surrogate AGB”) was computed. Later, a regression kriging interpolation was performed by taking into account the spatial structure of residuals to provide a continuous residual map. Finally, the RS-derived AGB map and the residual map were summed, and a final AGB map was obtained. The results showed that the integration of GLAS surrogate AGB data slightly improves the accuracy of the RS-derived AGB map only for AGB values lower than 100 t/ha (bias and RMSE reduced by 13.9 and 10 t/ha, respectively).

Integration of multi-resource remotely sensed data and allometric models for forest aboveground biomass estimation in China

Quantification of forest aboveground biomass density (AGB) is useful in forest carbon cycle studies, biodiversity protection and climate-change mitigation actions. However, a finer resolution and spatially continuous forest AGB map is inaccessible at national level in China. In this study, we developed forest type- and ecozone-specific allometric models based on 1607 field plots. The allometric models were applied to Geoscience Laser Altimeter System (GLAS) data to calculate AGB at the footprint level. We then mapped a 30 m resolution national forest AGB by relating the GLAS footprint AGB to various variables derived from Landsat images and Phased Array L-band Synthetic Aperture Radar (PALSAR) data. We estimated the average forest AGB to be 69.88 Mg/ha with a standard deviation of 35.38 Mg/ha and the total AGB carbon stock to be 5.44 PgC in China. Our AGB estimates corresponded reasonably well with AGB inventories from the top ten provinces in the forested area, and the coefficient of determination and root mean square error were 0.73 and 20.65 Mg/ha, respectively. We found that the main uncertainties for AGB estimation could be attributed to errors in allometric models and in height measurements by the GLAS. We also found that Landsat-derived variables outperform PALSAR-derived variables and that the textural features of PALSAR better support forest AGB estimates than backscattered intensity.

A Forest Attribute Mapping Framework: A Pilot Study in a Northern Boreal Forest, Northwest Territories, Canada

A methods framework is presented that utilizes field plots, airborne light detection and ranging (LiDAR), and spaceborne Geoscience Laser Altimeter System (GLAS) data to estimate forest attributes over a 20 Mha area in Northern Canada. The framework was implemented to scale up forest attribute models from field data to intersecting airborne LiDAR data, and then to GLAS footprints. GLAS data were sequentially filtered and submitted to the k-nearest neighbour (k-NN) imputation algorithm to yield regional estimates of stand height and crown closure at a 30 m resolution. Resulting outputs were assessed against independent airborne LiDAR data to evaluate regional estimates of stand height (mean difference = −1 m, RMSE = 5 m) and crown closure (mean difference = −5%, RMSE = 9%). Additional assessments were performed as a function of dominant vegetation type and ecoregion to further evaluate regional products. These attributes form the primary descriptive structure attributes that are typical of forest inventory mapping programs, and provide insight into how they can be derived in northern boreal regions where field information and physical access is often limited.

利用波形匹配实现卫星激光测高脚点精确定位的方法

The footprints of a satellite laser altimeter have an elevation accuracy of the decimeter order, which satisfies the elevation accuracy needs of ground control points (GCP) for mapping. However, the horizontal accuracy of footprints is only few tens of meters, and only the footprints illuminating on flat ground targets can be used as GCPs. In this paper, the waveform model was derived and used to develop a waveform simulator of laser altimeters. Compared with the current method of waveform matching, the new simulator considered more detailed effects arising from the device and target, e.g., time and spatial distribution of lasers, surface profile, and surface reflectivity. The airborne LiDAR point cloud and the Geoscience Laser Altimeter System(GLAS) data were involved to match the best-fit waveform by maximizing the correlation coefficients of the simulated waveform and GLAS waveform, and the precise location of every GLAS footprint could be acquired where the correlation coefficient was the maximum. The results show that, the mean of maximum correlation coefficients is more than 0.9 during the GLAS operating periods with normal received energies, and the horizontal accuracy of footprints is approximately 2 m after the waveform matching. The proposed method can be used to extract the laser GCPs on complex relief of the surface.

Continental Estimates of Canopy Gap Fraction by Active Remote Sensing

ABSTRACTThis study estimates canopy gap fraction (GF) from Geoscience Laser Altimeter System (GLAS) waveform data retrieved from a cross section of the Australian ecosystem, represented by 4 sites. However, as the literature supports, waveform profiles from systems such as GLAS require scaling in order to retrieve more accurate estimates of GF; the degree of scaling required is not a priori known, which has hindered large-scale estimates of GF by such systems to date. We employ spatially coincident GLAS footprints and airborne laser scanning (ALS) point cloud data from each site to scale initial estimates of GLAS GF to match ALS estimates, which allows the inference of scale factors to all footprints across the Australian forest landscape via the Random Forest (RF) technique. These refined estimates of GF are submitted to a second RF model to yield regional predictions of GF across the forest landscape. These predictions are validated against ALS data from 20 sites throughout Australia (R2 = 0.56, RMSE = 25.95%) and compared to MODIS, and also compared to ALS data (R2 = 0.13, RMSE = 66.46%). The developed method for producing large-scale estimates of GF shows promise for future studies and refinement, and is particularly pertinent as future spaceborne waveform systems come online (ICESat-2 and the Global Ecosystem Dynamics Investigation LiDAR).

Interest of Integrating Spaceborne LiDAR Data to Improve the Estimation of Biomass in High Biomass Forested Areas

Mapping forest AGB (Above Ground Biomass) is of crucial importance to estimate the carbon emissions associated with tropical deforestation. This study proposes a method to overcome the saturation at high AGB values of existing AGB map (Vieilledent’s AGB map) by using a map of correction factors generated from GLAS (Geoscience Laser Altimeter System) spaceborne LiDAR data. The Vieilledent’s AGB map of Madagascar was established using optical images, with parameters calculated from the SRTM Digital Elevation Model, climatic variables, and field inventories. In the present study, first, GLAS LiDAR data were used to obtain a spatially distributed (GLAS footprints geolocation) estimation of AGB (GLAS AGB) covering Madagascar forested areas, with a density of 0.52 footprint/km2. Second, the difference between the AGB from the Vieilledent’s AGB map and GLAS AGB at each GLAS footprint location was calculated, and additional spatially distributed correction factors were obtained. Third, an ordinary kriging interpolation was thus performed by taking into account the spatial structure of these additional correction factors to provide a continuous correction factor map. Finally, the existing and the correction factor maps were summed to improve the Vieilledent’s AGB map. The results showed that the integration of GLAS data improves the precision of Vieilledent’s AGB map by approximately 7 t/ha. By integrating GLAS data, the RMSE on AGB estimates decreases from 81 t/ha (R2 = 0.62) to 74.1 t/ha (R2 = 0.71). Most importantly, we showed that this approach using LiDAR data avoids underestimating high biomass values (new maximum AGB of 650 t/ha compared to 550 t/ha with the first approach).

Aboveground biomass mapping in French Guiana by combining remote sensing, forest inventories and environmental data

Mapping forest aboveground biomass (AGB) has become an important task, particularly for the reporting of carbon stocks and changes. AGB can be mapped using synthetic aperture radar data (SAR) or passive optical data. However, these data are insensitive to high AGB levels (>150Mg/ha, and >300Mg/ha for P-band), which are commonly found in tropical forests. Studies have mapped the rough variations in AGB by combining optical and environmental data at regional and global scales. Nevertheless, these maps cannot represent local variations in AGB in tropical forests. In this paper, we hypothesize that the problem of misrepresenting local variations in AGB and AGB estimation with good precision occurs because of both methodological limits (signal saturation or dilution bias) and a lack of adequate calibration data in this range of AGB values. We test this hypothesis by developing a calibrated regression model to predict variations in high AGB values (mean >300Mg/ha) in French Guiana by a methodological approach for spatial extrapolation with data from the optical geoscience laser altimeter system (GLAS), forest inventories, radar, optics, and environmental variables for spatial inter- and extrapolation. Given their higher point count, GLAS data allow a wider coverage of AGB values. We find that the metrics from GLAS footprints are correlated with field AGB estimations (R2=0.54, RMSE=48.3Mg/ha) with no bias for high values. First, predictive models, including remote-sensing, environmental variables and spatial correlation functions, allow us to obtain “wall-to-wall” AGB maps over French Guiana with an RMSE for the in situ AGB estimates of ∼50Mg/ha and R2=0.66 at a 1-km grid size. We conclude that a calibrated regression model based on GLAS with dependent environmental data can produce good AGB predictions even for high AGB values if the calibration data fit the AGB range. We also demonstrate that small temporal and spatial mismatches between field data and GLAS footprints are not a problem for regional and global calibrated regression models because field data aim to predict large and deep tendencies in AGB variations from environmental gradients and do not aim to represent high but stochastic and temporally limited variations from forest dynamics. Thus, we advocate including a greater variety of data, even if less precise and shifted, to better represent high AGB values in global models and to improve the fitting of these models for high values.

ESTIMATION OF REGIONAL FOREST ABOVEGROUND BIOMASS COMBINING ICESAT-GLAS WAVEFORMS AND HJ-1A/HSI HYPERSPECTRAL IMAGERIES

Estimation of forest aboveground biomass (AGB) is a critical challenge for understanding the global carbon cycle because it dominates the dynamics of the terrestrial carbon cycle. Light Detection and Ranging (LiDAR) system has a unique capability for estimating accurately forest canopy height, which has a direct relationship and can provide better understanding to the forest AGB. The Geoscience Laser Altimeter System (GLAS) onboard the Ice, Cloud, and land Elevation Satellite (ICESat) is the first polarorbiting LiDAR instrument for global observations of Earth, and it has been widely used for extracting forest AGB with footprints of nominally 70 m in diameter on the earth's surface. However, the GLAS footprints are discrete geographically, and thus it has been restricted to produce the regional full coverage of forest AGB. To overcome the limit of discontinuity, the Hyper Spectral Imager (HSI) of HJ-1A with 115 bands was combined with GLAS waveforms to predict the regional forest AGB in the study. Corresponding with the field investigation in Wangqing of Changbai Mountain, China, the GLAS waveform metrics were derived and employed to establish the AGB model, which was used further for estimating the AGB within GLAS footprints. For HSI imagery, the Minimum Noise Fraction (MNF) method was used to decrease noise and reduce the dimensionality of spectral bands, and consequently the first three of MNF were able to offer almost 98% spectral information and qualified to regress with the GLAS estimated AGB. Afterwards, the support vector regression (SVR) method was employed in the study to establish the relationship between GLAS estimated AGB and three of HSI MNF (i.e. <i>MNF1</i>, <i>MNF2</i> and <i>MNF3</i>), and accordingly the full covered regional forest AGB map was produced. The results showed that the adj.R<sup>2</sup> and RMSE of SVR-AGB models were 0.75 and 4.68 t hm<sup>−2</sup> for broadleaf forests, 0.73 and 5.39 t hm<sup>−2</sup> for coniferous forests and 0.71 and 6.15 t hm<sup>−2</sup> for mixed forests respectively. The full covered regional forest AGB map of the study area had 0.62 of accuracy and 11.11 t hm<sup>−2</sup> of RMSE. The study demonstrated that it holds great potential to achieve the full covered regional forest AGB distribution with higher accuracy by combing LiDAR data and hyperspectral imageries.

ESTIMATION OF REGIONAL FOREST ABOVEGROUND BIOMASS COMBINING ICESAT-GLAS WAVEFORMS AND HJ-1A/HSI HYPERSPECTRAL IMAGERIES

Abstract. Estimation of forest aboveground biomass (AGB) is a critical challenge for understanding the global carbon cycle because it dominates the dynamics of the terrestrial carbon cycle. Light Detection and Ranging (LiDAR) system has a unique capability for estimating accurately forest canopy height, which has a direct relationship and can provide better understanding to the forest AGB. The Geoscience Laser Altimeter System (GLAS) onboard the Ice, Cloud, and land Elevation Satellite (ICESat) is the first polarorbiting LiDAR instrument for global observations of Earth, and it has been widely used for extracting forest AGB with footprints of nominally 70 m in diameter on the earth's surface. However, the GLAS footprints are discrete geographically, and thus it has been restricted to produce the regional full coverage of forest AGB. To overcome the limit of discontinuity, the Hyper Spectral Imager (HSI) of HJ-1A with 115 bands was combined with GLAS waveforms to predict the regional forest AGB in the study. Corresponding with the field investigation in Wangqing of Changbai Mountain, China, the GLAS waveform metrics were derived and employed to establish the AGB model, which was used further for estimating the AGB within GLAS footprints. For HSI imagery, the Minimum Noise Fraction (MNF) method was used to decrease noise and reduce the dimensionality of spectral bands, and consequently the first three of MNF were able to offer almost 98% spectral information and qualified to regress with the GLAS estimated AGB. Afterwards, the support vector regression (SVR) method was employed in the study to establish the relationship between GLAS estimated AGB and three of HSI MNF (i.e. MNF1, MNF2 and MNF3), and accordingly the full covered regional forest AGB map was produced. The results showed that the adj.R2 and RMSE of SVR-AGB models were 0.75 and 4.68 t hm−2 for broadleaf forests, 0.73 and 5.39 t hm−2 for coniferous forests and 0.71 and 6.15 t hm−2 for mixed forests respectively. The full covered regional forest AGB map of the study area had 0.62 of accuracy and 11.11 t hm−2 of RMSE. The study demonstrated that it holds great potential to achieve the full covered regional forest AGB distribution with higher accuracy by combing LiDAR data and hyperspectral imageries.

Abiotic Controls on Macroscale Variations of Humid Tropical Forest Height

Spatial variation of tropical forest tree height is a key indicator of ecological processes associated with forest growth and carbon dynamics. Here we examine the macroscale variations of tree height of humid tropical forests across three continents and quantify the climate and edaphic controls on these variations. Forest tree heights are systematically sampled across global humid tropical forests with more than 2.5 million measurements from Geoscience Laser Altimeter System (GLAS) satellite observations (2004–2008). We used top canopy height (TCH) of GLAS footprints to grid the statistical mean and variance and the 90 percentile height of samples at 0.5 degrees to capture the regional variability of average and large trees globally. We used the spatial regression method (spatial eigenvector mapping-SEVM) to evaluate the contributions of climate, soil and topography in explaining and predicting the regional variations of forest height. Statistical models suggest that climate, soil, topography, and spatial contextual information together can explain more than 60% of the observed forest height variation, while climate and soil jointly explain 30% of the height variations. Soil basics, including physical compositions such as clay and sand contents, chemical properties such as PH values and cation-exchange capacity, as well as biological variables such as the depth of organic matter, all present independent but statistically significant relationships to forest height across three continents. We found significant relations between the precipitation and tree height with shorter trees on the average in areas of higher annual water stress, and large trees occurring in areas with low stress and higher annual precipitation but with significant differences across the continents. Our results confirm other landscape and regional studies by showing that soil fertility, topography and climate may jointly control a significant variation of forest height and influencing patterns of aboveground biomass stocks and dynamics. Other factors such as biotic and disturbance regimes, not included in this study, may have less influence on regional variations but strongly mediate landscape and small-scale forest structure and dynamics.

COMPARING ICESAT/GLAS BASED ELEVATION HEIGHTS WITH PHOTOGRAMMETRIC TERRAIN HEIGHTS FROM UAV-IMAGERY ON THE EAST TIBETAN PLATEAU

Abstract. Digital elevation models (DEMs) and height measurements are broadly used in environmental studies. Two common elevation sources are the Ice Cloud and land elevation Satellite (ICESat), which acquired laser range measurements with the Geoscience Laser Altimeter System (GLAS) across the globe and elevation data from the Shuttle Radar Topography Mission (SRTM). Current developments of small unmanned aerial vehicles (UAV) provide the opportunity to collect aerial images of remote areas at a high spatial resolution. These can be further processed to digital surface models by stereophotogrammetry and provide a reliable data source to evaluate coarse scale Digital Elevation Models (DEMs). This study compares ICESat/GLAS and SRTM90 elevation data against photogrammetric terrain heights within GLAS footprints on high altitudes on the East Tibetan Plateau. Without vegetation-bias, we were able to examine height differences under different topographic conditions and of different acquisition dates. Several resampling techniques were applied to SRTM90 data and averaged height within each footprint was calculated. ICESat/GLAS heights (n = 148) are most similar to UAV data based elevations with an averaged difference of −0.8m ±3.1m. Results furthermore indicate the validity of ICESat/GLAS heights, which are usually removed from analyses by applying different quality flags. Smallest difference of SRTM90 to UAV based heights could be observed by a natural neighbour resampling technique (averaged 3.6m ±14m), whereat other techniques achieved quite similar results. It can be confirmed that within a range of 3,800–4,200m above mean sea level the ICESat/GLAS heights are a precise source to determine elevation at footprint geolocation.

National Forest Aboveground Biomass Mapping from ICESat/GLAS Data and MODIS Imagery in China

Forest aboveground biomass (AGB) was mapped throughout China using large footprint LiDAR waveform data from the Geoscience Laser Altimeter System (GLAS) onboard NASA’s Ice, Cloud, and land Elevation Satellite (ICESat), Moderate Resolution Imaging Spectro-radiometer (MODIS) imagery and forest inventory data. The entire land of China was divided into seven zones according to the geographic characteristics of the forests. The forest AGB prediction models were separately developed for different forest types in each of the seven forest zones at GLAS footprint level from GLAS waveform parameters and biomass derived from height and diameter at breast height (DBH) field observation. Some waveform parameters used in the prediction models were able to reduce the effects of slope on biomass estimation. The models of GLAS-based biomass estimates were developed by using GLAS footprints with slopes less than 20° and slopes ≥ 20°, respectively. Then, all GLAS footprint biomass and MODIS data were used to establish Random Forest regression models for extrapolating footprint AGB to a nationwide scale. The total amount of estimated AGB in Chinese forests around 2006 was about 12,622 Mt vs. 12,617 Mt derived from the seventh national forest resource inventory data. Nearly half of all provinces showed a relative error (%) of less than 20%, and 80% of total provinces had relative errors less than 50%.

Estimates of forest structure parameters from GLAS data and multi-angle imaging spectrometer data

Quantitative estimates of forest vertical and spatial distribution using remote sensing technology play an important role in better understanding forest ecosystem function, forest carbon storage and the global carbon cycle. Although most remote sensing systems can provide horizontal distribution of canopies, information concerning the vertical distribution of canopies cannot be detected. Fortunately, laser radars have become available, such as GLAS (Geoscience laser altimeter system). Because laser radar can penetrate foliage, it is superior to other remote sensing technologies for detecting vertical forest structure and has higher accuracy. GLAS waveform data were used in this study to retrieve average tree height and biomass in a GLAS footprint area in Heilongjiang Province. However, GLAS data are not spatially continuous. To fill the gaps, MISR (multi- angle imaging spectrometer) spectral radiance was chosen to predict the regional continuous tree height by developing a multivariate linear regression model. We compared tree height estimated by the regression model and GLAS data. The results confirmed that estimates of tree height and biomass based on GLAS data are considerably more accurate than estimates based on traditional methods. The accuracy is approximately 90%. MISR can be used to estimate tree height in continuous areas with a robust regression model. The R2, precision and root mean square error of the regression model were 0.8, 83% and 1m, respectively. This study provides an important reference for mapping forest vertical parameters.

Sub-footprint analysis to uncover tree height variation using ICESat/GLAS

Detailed forest height data are an indispensable prerequisite for many forestry and earth science applications. Existing research of using Geoscience Laser Altimeter System (GLAS) data mainly focuses on deriving average or maximum tree heights within a GLAS footprint, i.e. an ellipse with a diameter of 65m. However, in most forests, it is likely that the tree heights within such ellipse are heterogeneous. Therefore, it is desired to uncover detailed tree height variation within a GLAS footprint. To the best of our knowledge, no such methods have been reported as of now. In this study, we aim to characterize tree heights’ variation within a GLAS footprint as different layers, each of which corresponds to trees with similar heights. As such, we developed a new method that embraces two steps: first, a refined Levenberg–Marquardt (LM) algorithm is proposed to decompose raw GLAS waveform into multiple Gaussian signals, within which it is hypothesized that each vegetation signal corresponds to a particular tree height layer. Second, for each layer, three parameters were first defined: Canopy Top Height (CTH), Crown Length (CL), and Cover Proportion (CP). Then we extracted the three parameters from each Gaussian signal through a defined model. In order to test our developed method, we set up a study site in Ejina, China where the dominant specie is Populus euphratica. Both simulated and field tree height data were adopted. With regard to the simulation data, results presented a very high agreement for the three predefined parameters between our results and simulation data. When our methods were applied to the field data, the respective R2 become 0.78 (CTH), CL (R2=0.76), CP (R2=0.74). Overall, our studies revealed that large footprint GLAS waveform data have the potentials for obtaining detailed forest height variation.

Forest canopy height mapping over China using GLAS and MODIS data

The Geoscience Laser Altimeter System (GLAS) accurately detects the vertical structural information of a target within its laser spot and is a promising system for the inversion of structural features and other biophysical parameters of forest ecosystems. Since the GLAS footprints are discontinuously distributed with a relativity low density, continuous vegetation height distributions cannot be mapped with a high accuracy using GLAS data alone. The MODIS BRDF product provides more forest structural information than other optical remote sensing data. This study aimed to map forest canopy heights over China from the GLAS and MODIS BRDF data. Firstly, the waveform characteristic parameters were extracted from the GLAS data by the method of wavelet analysis, and the terrain index was calculated using the ASTER GDEM data. Secondly, the model reducing the topographic influence was constructed from the waveform characteristic parameters and terrain index. Thirdly, the final canopy height estimation model was constructed from the neural network combining the canopy height estimated with the GLAS point and the MODIS BRDF data, and applied to get the continuous canopy height map over China. Finally, the map was validated by the measured data and the airborne LiDAR data, and the validation results indicated that forest canopy heights can be estimated with high accuracy from combined GLAS and MODIS data.

Large-scale retrieval of leaf area index and vertical foliage profile from the spaceborne waveform lidar (GLAS/ICESat)

Leaf area index (LAI) and canopy vertical profiles are important descriptors of ecosystem structure. The Geoscience Laser Altimeter System (GLAS) on board ICESat (Ice, Cloud, and land Elevation Satellite) provided three-dimensional observations that can be used to derive these canopy structure parameters globally. While several canopy height products have been produced globally from GLAS, no comparable data sets for LAI and canopy profiles exist across large areas. In this study we develop a physically based method of retrieving LAI and vertical foliage profiles (VFPs) from GLAS observations over the entire state of California, USA. This method refines lidar derived LAI and VFP through a recursive analysis of GLAS waveforms using ancillary data obtained from existing remote sensing products. Those supplemental inputs include canopy clumping index derived from POLDER, 500m land cover type and 1km LAI data derived from MODIS. Implementation of our method created state-level LAI and VFP data for the existing GLAS transects over California. We then analyzed the variability of LAI and VFP data sets across environmental gradients and as a function of land cover type and elevation. Both LAI and VFP showed strong variability across elevational gradients and among land cover types. We compared our results at the scale of GLAS footprints with an LAI map derived from Landsat (at 30m) and found moderate agreement (r2=0.34, bias=0.26, RMSD (Root Mean Square Difference)=1.85) between the two. In particular, Landsat LAI not only appeared to saturate relative to GLAS LAI at around LAI=5, but also showed an overestimation for LAI less than about 2. Best agreement between the two LAI data sets was shown to occur in areas with slope less than 20°. Results from our study suggest the possibility of retrieving global LAI and VFP data from GLAS data and the potential for synergetic observation of lidar and passive optical remote sensing data such as Landsat.