Capability of GLAS/ICESat Data to Estimate Forest Canopy Height and Volume in Mountainous Forests of Iran

Abstract

The importance of measuring forest biophysical properties for ecosystem health monitoring and forest management encourages researchers to find precise, yet low-cost methods especially in mountainous and large areas. In the present study, geoscience laser altimeter system (GLAS) on board Ice, Cloud, and land Elevation Satellite (ICESat) was used to estimate three biophysical characteristics of forests located in the north of Iran: 1) maximum canopy height ( ${\text {H}_{{\max}}}$ ); 2) Lorey’s height ( ${{\text H}_{\text{Lorey}}}$ ); and 3) forest volume (V). A large number of multiple linear regressions (MLR) and also random forest (RF) regressions were developed using different sets of variables including waveform metrics, principal components (PCs) produced from principal component analysis (PCA) and wavelet coefficients (WCs) generated from wavelet transformation (WT). To validate and compare models, statistical criteria were calculated based on a fivefold cross validation. Best model concerning the maximum height was an MLR ( ${\text {RMSE}} = {5.0}\;{\text m}$ ) which combined two metrics extracted from waveforms (waveform extent “ ${{\text W}_{\text{ext}}}$ ” and height at 50% of waveform energy “ ${{\text H}_{{50}}}$ ”), and one from digital elevation model (terrain index, TI). The mean absolute percentage error (MAPE) of maximum height estimates was 16.4%. For Lorey’s height, a simple MLR (including ${{\text W}_{\text{ext}}}$ and TI) represented the highest performance ( ${\text {RMSE}} = {5.1}\;{\text m}$ , ${\text {MAPE}} = {24.0}\% $ ). Generally, MLR models had a better performance when compared to the RF models. In addition, the accuracy of height estimations using waveform metrics was greater than those based on PCs or WCs. Concerning forest volume, regression models estimating volume directly from GLAS data led to a better result ( ${\text {RMSE}} = {128.8}\;{{\text m}^{3}}/\text {ha}$ ) rather than volume– ${{\text H}_{\text{Lorey}}}$ relationship ( ${\text {RMSE}} = {167.8}\;{{\text m}^{3}}/\text {ha}$ ).