Comparing airborne and spaceborne photon-counting LiDAR canopy structural estimates across different boreal forest types

Abstract

The monitoring of forested ecosystems relies on an accurate description of forest structure. The Ice, Cloud and land Elevation Satellite-2 (ICESat-2), launched in September 2018, carries the Advanced Topographic Laser Altimeter System (ATLAS), a Light Detection and Ranging (LiDAR) instrument capable of detecting individual photons reflected back from vegetation canopy. ICESat-2 data is delivering global estimates of forest structure; however, analysis of the performance of ICESat-2 on-orbit data across a range of forest conditions remains limited. This study derives structural estimates of (i) canopy height, (ii) canopy cover and (iii) canopy height variability from ICESat-2 data acquired in snow-free and low atmospheric scattering conditions over different boreal forest structural types in Ontario, Canada. ICESat-2 structural estimates were derived from the Global Geolocated Photon Data (ATL03) and Land and Vegetation Height (ATL08) data products and compared against single-photon detection airborne LiDAR (Leica SPL100). An extensive network of ground plots were used to stratify the study area into three distinct forest structural groups, each resulting from different stand development stages. ICESat-2 and SPL100 estimates of canopy height were compared at the ATL03 photon level, whereas estimates of height variability and canopy cover were compared for spatial analysis units (AU; mean size = 1287 m2). ICESat-2 photons returned from the top of the canopy underestimated canopy height relative to SPL100 by an average of 2.3 m overall and corresponded most strongly to the 90th percentile (P90) of coincident airborne SPL100 returns (root mean square difference (RMSD) = 2.9 m and correlation coefficient (r) = 0.84). The lowest average underestimation of SPL P90 was observed in homogeneous stands that were relatively simple, and single-layered with a single dominant species (RMSD = 2.5 m, r = 0.84). We observed the least agreement of ICESat-2 and SPL forest structural metrics in over-mature stands with complex structure and greater variability in canopy heights (RMSD = 3.5 m, r = 0.64). For the AUs, the strength of the relationship between SPL100 and ICESat-2 canopy height percentiles increased with increasing height percentiles (e.g. P25 RMSD% = 77.8%; P95 RMSD% = 23.7%). ICESat-2 generally underestimated canopy height variability relative to the SPL100 data, with both data having similar absolute variability (standard deviation of canopy heights RMSD% = 26.8%, r = 0.75), but lower agreement in relative variability (coefficient of variation of canopy heights RMSD% = 33.9%, r = 0.45). Herein we propose the use of the vegetation fill index as a method to estimate canopy cover with ICESat-2. Comparison of SPL100 and ICESat-2 vegetation fill indices at the AU level resulted in strong agreement overall (RMSD% = 19.7%; r = 0.57). These observations and results contribute to the overall objective of building a comprehensive understanding of the performance of ICESat-2 for characterizing vegetation structure in boreal forest environments.