Thematic mapping of complex landscapes, with various phenological patterns from satellite imagery, is a particularly challenging task. However, supplementary information, such as multitemporal data and/or land surface temperature (LST), has the potential to improve the land cover classification accuracy and efficiency. In this paper, in order to map land covers, we evaluated the potential of multitemporal Landsat 8's spectral and thermal imageries using a random forest (RF) classifier. We used a grid search approach based on the out-of-bag (OOB) estimate of error to optimize the RF parameters. Four different scenarios were considered in this research: (1) RF classification of multitemporal spectral images, (2) RF classification of multitemporal LST images, (3) RF classification of all multitemporal LST and spectral images, and (4) RF classification of selected important or optimum features. The study area in this research was Naghadeh city and its surrounding region, located in West Azerbaijan Province, northwest of Iran. The overall accuracies of first, second, third, and fourth scenarios were equal to 86.48, 82.26, 90.63, and 91.82%, respectively. The quantitative assessments of the results demonstrated that the most important or optimum features increase the class separability, while the spectral and thermal features produced a more moderate increase in the land cover mapping accuracy. In addition, the contribution of the multitemporal thermal information led to a considerable increase in the user and producer accuracies of classes with a rapid temporal change behavior, such as crops and vegetation.