State-wide forest canopy height and aboveground biomass map for New York with 10 m resolution, integrating GEDI, Sentinel-1, and Sentinel-2 data

Investigating the quantity of forest aboveground biomass (AGB) is crucial for understanding the role forests play in the global carbon cycle. The canopy height model (CHM) is a critical component in estimating AGB, as it provides a three-dimensional representation of the tree canopy. Traditional CHM estimation methods are time-consuming, labor-intensive, and expensive, particularly at large-scales. Remote sensing is a cost-effective and efficient alternative approach, providing valuable information over large areas in a timely manner. The Global Ecosystem Dynamics Investigation (GEDI) onboard the International Space Station is a space-based light detection and ranging (LiDAR) system designed to collect information on vertical structures of vegetation. One major problem with the collection of GEDI data is that it provides limited information over discrete ground samples, also known as footprints, and thus do not provide wall-to-wall gridded height products. The objective of this study was twofold: a) to integrate the GEDI LiDAR footprint heights with Sentinel-2 multispectral imagery to generate a 10 m wall-to-wall CHM map of New York State (NYS), USA for the year 2019 and b) to improve our previously generated AGB map (both accuracy and resolution) of NYS for the year 2019 by fusing Sentinel-2 multispectral, Sentinel-1 synthetic aperture radar (SAR), and the produced CHM. To generate the 10 m CHM map, the GEDI footprints height measurements were extrapolated using Sentinel-2 imagery and a random forest model. The CHM that was produced was assessed by using GEDI footprints that were not part of the training phase and were therefore independent (extrapolated). Comparing our 10 m CHM with the available global 30 m CHM map provided by Potapov et al. (2021) over NYS shows significant improvement not only in terms of spatial resolution, but also in terms of accuracy. The root mean square error (RMSE) of our 10 m CHM is 4.4 m while this value is 7.49 m for the 30 m CHM over NYS. Similarly, the R² value for the 10 m CHM map is 0.74, while that of the 30 m CHM is 0.46. Finally, the integration of produced 10 m CHM, Sentinel-1, and Sentinel-2 datasets were utilized to create a 10 m AGB map of NYS with the RMSE of 39.49 Mg/ha, and R² of 0.65. The results demonstrate the potential of integrating GEDI, Sentinel-1, and Sentinel-2 data for providing a valuable tool for large-scale mapping of forest canopy structure and biomass, which can help to inform forest management and carbon accounting efforts.