Abstract

Mapping the spatial distribution of soil organic carbon (SOC) in lands covered by tropical forests is important to understand the relationship and dynamics of SOC in this type of ecosystem. In this study, the Random Forest (RF) algorithm was used to map SOC stocks of topsoil (0–15 cm) in forest lands of the Dominican Republic. The methodology was developed using geospatial datasets available in the Google Earth Engine (GEE) platform combined with a set of 268 soil samples. Twenty environmental covariates were analyzed, including climate, topography, and vegetation. The results indicate that Model A (combining all 20 covariates) was only marginally better than Model B (combining topographic and climatic covariates), and Model C (only combining multispectral remote sensing data derived from Landsat 8 OLI images). Model A and Model B yielded SOC mean values of 110.35 and 110.87 Mg C ha−1, respectively. Model A reported the lowest prediction error and uncertainty with an R2 of 0.83, an RMSE of 35.02 Mg C ha−1. There was a strong dependence of SOC stocks on multispectral remote sensing data. Therefore, multispectral remote sensing proved accurate to map SOC stocks in forest ecosystems in the region.