Evaluation and modelling of the Persian Leopard’s (Panthera pardus saxicolor) habitat suitability in the Gisakan area, Bushehr Province, Iran

Being at the top of the food pyramid, as a pioneer species, has led the Persian leopard (Panthera pardus saxicolor) to expose habitat degradation and prey reduction, in the recent years, in particular, necessitating the development of habitat conservation plans. In this study, habitat suitability modeling of the Persian Leopard was conducted using QGIS with only presence data as well as general linear model (GLM) with presence/absence data based on the binomial distribution, run in the R software. In QGIS, using topography, DEM, slope, and aspect maps as well as villages maps, the continuous population distribution map was delineated. The presence points, derived from field-working, were then inserted, and other mandatory layers were also collected and geo-rectified, where applicable. Afterward, each one of the layers was categorized and rasterized. Following the validation, the final suitability map of the Persian leopard was performed. In the GLM, simulations were generated using a full model with all variables and a leave-one-out procedure, gradually dropping the least contributing variable to form a minimal model. The results showed throughout the area, the most suitable habitats are scattered. Considering the substantial reduction in the preys population in Gisakan alongside the abundance of livestock, goats specifically, habitat suitability is largely increased near the rangelands and water resources. Furthermore, the habitats placed in northern aspects with > 1000 m height, slopes of 30-65%, 0 to 1000 meters distance to the roads, and 1000 to 2500 meters distance to the villages, showed higher suitability.