A Comparison of Machine Learning and Geostatistical Approaches for Mapping Forest Canopy Height over the Southeastern US Using ICESat-2

Abstract

The availability of canopy height information in the Ice, Cloud, and Land Elevation Satellite-2’s (ICESat-2’s) land and vegetation product, or ATL08, presents opportunities for developing full-coverage products over broad spatial scales. The primary goal of this study was to develop a 30-meter canopy height map over the southeastern US, for the Southeastern Plains ecoregion and the Middle Atlantic Coastal Plains ecoregion. More specifically, this work served to compare well-known modeling approaches for upscaling canopy information from ATL08 to develop a wall-to-wall product. Focusing on only strong beams from nighttime acquisitions, the h_canopy parameter was extracted from ATL08 data. Landsat-8 bands and derived vegetation indices (normalized difference vegetation index, enhanced vegetation index, and modified soil-adjusted vegetation index) along with National Land Cover Database’s canopy cover and digital elevation models were used to extrapolate ICESat-2 canopy height from tracks to the regional level. Two different modeling techniques, random forest (RF) and regression kriging (RK), were applied for estimating canopy height. The RF model estimated canopy height with a coefficient of determination (R2) value of 0.48, root-mean-square error (RMSE) of 4.58 m, mean absolute error (MAE) of 3.47 and bias of 0.23 for independent validation, and an R2 value of 0.38, RMSE of 6.39 m, MAE of 5.04 and bias of −1.39 when compared with airborne lidar-derived canopy heights. The RK model estimated canopy heights with an R2 value of 0.69, RMSE of 3.49 m, MAE of 2.61 and bias of 0.03 for independent validation, and an R value of 0.68, R2 value of 0.47, RMSE of 5.96m, MAE of 4.52 and bias of −1.81 when compared with airborne lidar-derived canopy heights. The results suggest feasibility for the implementation of the RK method over a larger spatial extent and potential for combining other remote sensing and satellite data for future monitoring of canopy height dynamics.