High-resolution forest age mapping based on forest height maps derived from GEDI and ICESat-2 space-borne lidar data

Forest age is a key parameter for estimating forest growth and carbon uptake and for forest management. Remote sensing provides indirect but useful information for mapping forest age at large scales. However, existing regional and global forest age products were generated at low spatial resolutions (often 1000 m) and are not useful for most forest stands in China that are smaller than 1000 m. This study aims to map forest age at the 30 m resolution based on forest height maps mainly derived from the Global Ecosystem Dynamics Investigation (GEDI) and Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) data, and analyze the roles of auxiliary data including temperature, precipitation, slope, and aspect in forest age mapping. Forest age is defined as the average age of dominant tree species within a pixel. Five commonly-used stand growth equations and twelve machine learning methods were tested for their suitability for mapping forest age of different tree species. We found that the Logistic equation performed the best among the tested stand growth equations and the Random Forest (RF) was the best among the tested machine learning methods. According to RF, forest height contributed predominantly to the variance in forest age mapping, while temperature, precipitation, slope, and aspect also had an overall non-negligible and variable contribution among different tree species. By integrating the climate and topographical variables, RF was applicable for forest age mapping without classifying the tree species. These results show that forest height maps derived from space-borne lidar data such as GEDI and ICESat-2 data are highly useful for mapping forest stand age, and the methodology developed in this study highlights a perspective for generating national and global forest age products at a high spatial resolution.