ABSTRACTMaximum seasonal thaw depth, referred to as active layer thickness (ALT), is one of the key parameters used to monitor permafrost conditions. ALT maps based on interpolation of point measurements or derived from coarse or moderate spatial resolution satellite data often hide small‐scale spatial variations in thaw depth resulting from differences in surface characteristics and microtopography. To model and predict changes in hydrological and biogeochemical processes in permafrost areas accurately, high‐resolution remote sensing‐based estimations of ALT are needed. Therefore, we applied random forest (RF) regression on a set of topographical and spectral vegetation indices derived from optical unoccupied aerial systems data, Landsat 8 land surface temperature (LST) data, and field measurements to estimate thaw depths in palsas at three mires in north‐west Finland. We also analyzed differences in thaw depths between mires located at different elevations, between dome and plateau‐shaped palsas, and between different vegetation and surface cover classes. The RF models resulted in root mean square errors from 2.4 to 5.7 cm between predicted and observed thaw depths and the R2 values of 0.57–0.96. Height from the surrounding fen surface and LST were the most important variables in thaw depth models, although high‐accuracy results were also achieved without LST. The mean thaw depths did not differ between the sites with lowest and highest elevation, whereas the thaw depths were significantly deeper in dome‐shaped palsas compared to plateaus. The thaw depths were significantly different between vegetation cover classes only on plateau‐shaped palsas. The results indicate the high impact of the topography on the palsa thaw depth, thus highlighting the importance of accurate elevation models in spatial modeling of palsa ALT. The methodology presented in this study can be applied to other permafrost regions where field measurements of ALT are accompanied with high‐resolution topographical and multispectral data.