Multi-Resolution Mapping of Aboveground Biomass and Change in Puerto Rico’s Forests with Remote Sensing and Machine Learning
| Authors: | Nafiseh Haghtalab, Tamara Heartsill Scalley, Tana E. Wood, J. Aaron Hogan, Humfredo Marcano-Vega, Thomas Brandeis, Thomas Ruzycki, Eileen Helmer |
| Year: | 2026 |
| Type: | Scientific Journal |
| Station: | International Institute of Tropical Forestry |
| DOI: | https://doi.org/10.3390/rs18081190 |
| Source: | Remote Sensing |
Abstract
Tropical forests are major contributors to the global carbon budget but are affected by disturbances such as hurricanes, which cause extensive yet spatially variable tree damage and mortality. High-resolution maps of forest aboveground biomass (AGB) and its temporal change aid in quantifying disturbance impacts, assessing resilience, and supporting forest management. This study presents wall-to-wall, high-resolution mapping of pre- and post-hurricane AGB and AGB change across Puerto Rico. The maps represent forest AGB measured 0–2 years before and after two major hurricanes (Irma and Maria), as well as longer-term conditions up to four years post-disturbance. AGB was modeled using Random Forest (RF) algorithms that integrated Forest Inventory and Analysis (FIA) plot data with canopy height and cover derived from discrete-return LiDAR, multi-temporal satellite imagery, and additional geospatial predictors. Model performance was evaluated using a 10% holdout dataset. Predicted versus observed regressions yielded, at 10 m and 90 m spatial resolutions, respectively, r = 0.75 and 0.79 with model residual mean standard deviation (RMSD) = 87.7 and 39.2 Mg ha−1 for pre-hurricane AGB, and r = 0.77 and 0.74 with RMSD = 69.7 and 58.1 Mg ha−1 for post-hurricane AGB. AGB change models at 10 m and 90 m resolutions yielded r = 0.58 and 0.73 with RMSD = 17.0 and 18.7 Mg ha−1, respectively. Ten-fold cross-validation produced stronger correlations and reduced RMSD values. Frequency distributions of mapped pixels of forest AGB and AGB change, in comparison with previously published maps and island-wide field-based estimates, indicate that, although hurricane-driven biomass reductions of up to 20% were recorded in field data, patterns consistent with longer-term recovery from historical deforestation are evident within four years after the hurricanes. The 10 m maps capture fine-scale heterogeneity in canopy damage and regrowth, whereas the 90 m maps emphasize broader regional patterns. This integrated framework provides a transferable approach for monitoring forest structure and biomass dynamics in disturbance-prone tropical ecosystems.