Land use change and soil organic carbon stock (SOCS) depletion over time is one of the predominant world-wide environmental problems related to global warming and the need to secure food production for an increasing world population. In our research, satellite images from 1988 and 2018 were analyzed for a 177.48 km2 region in Kurdistan Province, Iran. Across the study area. 186 disturbed and undisturbed soil samples were collected at two depths (0–20 cm and 20–50 cm). Bulk density (BD), soil organic carbon (SOC), rock fragments (RockF) and SOCS were measured. Random forest was used to model the spatial variability of SOCS. Land use was mapped with supervised classification and maximum likelihood approaches. The Kappa index and overall accuracy of the supervised classification and maximum likelihood land use maps varied between 83 and 88% and 78 and 85%, respectively. The area of forest and high-quality rangeland covered 5286 ha in 1988 and decreased by almost 30% by 2018. Most of the decrease was due to the establishment of cropland and orchards, and due to overgrazing of high-quality rangeland. As expected, the results of the analysis of variance showed that mean values of SOCS for the high-quality rangeland and forest were significantly higher compared to other land use classes. Thus, transformation of land with natural vegetation like forest and high-quality rangeland led to a loss of 15,494 Mg C in the topsoil, 15,475 Mg C in the subsoil and 15,489 Mg C−1 in total. We concluded that the predominant causes of natural vegetation degradation in the study area were mostly due to the increasing need for food, anthropogenic activities such as cultivation and over grazing, lack of government landuse legislation and the results of this study are useful for land use monitoring, decision making, natural vegetation planning and other areas of research and development in Kurdistan province.
