Multi-Resolution Mapping of Aboveground Biomass and Change in Puerto Rico’s Forests with Remote Sensing and Machine Learning

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.