Archive/Forest Cover Change in the Nevado de Colima Using Sentinel-2 and an Enriched Random Forest Classifier with Slope and Spectral Indices
Forest Cover Change in the Nevado de Colima Using Sentinel-2 and an Enriched Random Forest Classifier with Slope and Spectral Indices
Guilherme Amorim Homem de Abreu Loureiro, Víctor David Cibrián-Llanderal, David Cibrián-Tovar
May 25, 2026
en

Abstract

Methodological opacity and the omission of environmental variables in forest masks can generate biased estimates. The objective of this study was to validate a reproducible workflow for quantifying forest cover change in the area adjacent to Nevado de Colima over the 2019–2025 period, subdivided into nine assessment areas with standardized sampling based on 3 × 3 pixel kernels (900 m2). An enriched Random Forest model with slope and spectral indices (NDVI, NBR, NDWI-Gao, and BSI) classified six spectral combinations derived from Sentinel-2 L2A bands B2, B3, B4, B8, B11, and B12, together with a new index proposed in this study, Red-Enhanced Normalized Burn Ratio (RE-NBR), used as a conservative classifier and auxiliary classifier output in the probabilistic cross-check estimation. Validation employed thematic and areal metrics. All combinations reached OA values between 89.44% and 92.53% and Kappa values between 0.79 and 0.85, with Shortwave Infrared (B12, B8, B4) as the most consistent configuration across dates. Allocation disagreement systematically exceeded quantity disagreement on all dates. The Seasonal Stability Index increased from 0.73 in 2019 to 0.77 in 2025, with persistent positive asymmetry between February and April. The probabilistic cross-check adjustment produced an adjusted forest loss of 1594.74 ha and an adjusted gain of 802.65 ha over 120,289.70 ha. Within the protected natural areas, expected change was distributed unevenly among vegetation types, with pine–oak forest showing the highest total expected loss, whereas high-mountain meadow showed the highest expected gain and also remained among the covers with the highest expected loss, indicating active spatial reconfiguration in the upper ecological domain where Pinus hartwegii Lindl. is the dominant species, though no species-level classification was performed. These results provide spatial evidence to support field verification, forest-health monitoring, and management decisions in the protected high-mountain study area.

Keywords

forestcoverchangenevadocolimasentinel-2enrichedrandomclassifierslopespectralindicesforestsmethodologicalopacityomissionenvironmentalvariablesmasksgeneratebiasedestimatesobjectivevalidate
Reference this publication

€ 4.00