The European rabbit (Oryctolagus cuniculus) presents a significant conservation and management challenge in Greece. While it has been listed in national biodiversity assessments, its population dynamics on the island of Lemnos demonstrate the characteristics of a highly adaptable and rapidly expanding species, exerting substantial ecological and economic impacts. Addressing this issue requires a spatially explicit understanding of its distribution patterns and habitat preferences, particularly given its extensive population growth over the past three decades. To this end, we conducted 40 field surveys across the island, documenting 1534 presence records of the species. We applied Kernel Density Estimation, Getis-Ord Gi *, and Anselin Local Moran’s I to identify the spatial distribution patterns and significant hotspots. A spatial lag model was used to quantify hotspot intensity and clustering dynamics, while abiotic, biotic, and anthropogenic factors were analyzed to assess habitat associations. Our results revealed that rabbit hotspots are predominantly concentrated in fertile lowland agroecosystems, with nearly 60% of high-density areas overlapping conservation zones. Soil and field conditions, grazing-supporting landscapes, and arable and subsidized agricultural areas emerged as significant predictors of O. cuniculus presence. The observed spatial dependencies indicated that while hotspot intensities and clustering dynamics are influenced by the conditions in neighboring areas, habitat characteristics remain fundamental in shaping their distribution, highlighting the broader landscape-scale spatial patterns affecting rabbit populations. These findings underscore the necessity of adopting spatially informed management strategies that mitigate agricultural impacts while accounting for interconnected spatial dynamics, providing a foundation for informed decision-making to manage rabbit populations while balancing conservation and agricultural priorities.

Understanding fine-scale spatial ecology is essential for defining effective conservation priorities, particularly at the range margins of vulnerable species. Here, we investigate the spatial ecology and habitat associations of the Persian squirrel (Sciurus anomalus) on Lesvos Island, Greece, representing the species’ westernmost distribution. Using a randomized grid-based survey, we recorded 424 presence records across the island and applied a suite of spatial analyses, including Kernel Density Estimation, Getis-Ord Gi*, and Anselin Local Moran’s I, to detect hotspots, coldspots, and spatial outliers. Binomial Logistic Regression, supported by Principal Component Analysis, identified key ecological drivers of habitat use, while spatial regression models (Spatial Lag and Spatial Error Models) quantified the influence of land-use characteristics and spatial dependencies on hotspot intensity and clustering dynamics. Our results showed that hotspots were primarily associated with olive-dominated and broadleaved landscapes, while coldspots and Low–Low clusters were concentrated in fragmented or degraded habitats, often outside protected areas. Spatial outliers revealed fine-scale deviations from broader patterns, indicating local habitat disruptions and emerging conservation risks not captured by existing Natura 2000 boundaries. Spatial regression confirmed that both hotspot intensity and clustering patterns were shaped by specific land-use features and spatially structured processes. Collectively, our findings underscore the fragmented nature of suitable habitats and the absence of cohesive population cores, reinforcing the need for connectivity-focused, landscape-scale conservation.