Density estimation is a key issue in wildlife management but is particularly challenging and labour-intensive for elusive species. Recently developed approaches based on remotely collected data and capture-recapture models, though representing a valid alternative to more traditional methods, have found little application to species with limited morphological variation. We implemented a camera trap capture-recapture study to survey wolf packs in a 560-km2 area of Central Italy. Individual recognition of focal animals (alpha) in the packs was possible by relying on morphological and behavioural traits and was validated by non-invasive genotyping and inter-observer agreement tests. Two types (Bayesian and likelihood-based) of spatially explicit capture-recapture (SCR) models were fitted on wolf pack capture histories, thus obtaining an estimation of pack density in the area. In two sessions of camera trapping surveys (2014 and 2015), we detected a maximum of 12 wolf packs. A Bayesian model implementing a half-normal detection function without a trap-specific response provided the most robust result, corresponding to a density of 1.21 ± 0.27 packs/100 km2 in 2015. Average pack size varied from 3.40 (summer 2014, excluding pups and lone-transient wolves) to 4.17 (late winter-spring 2015, excluding lone-transient wolves). We applied for the first time a camera-based SCR approach in wolves, providing the first robust estimate of wolf pack density for an area of Italy. We showed that this method is applicable to wolves under the following conditions: i) the existence of sufficient phenotypic/behavioural variation and the recognition of focal individuals (i.e. alpha, verified by non-invasive genotyping); ii) the investigated area is sufficiently large to include a minimum number of packs (ideally 10); iii) a pilot study is carried out to pursue an adequate sampling design and to train operators on individual wolf recognition. We believe that replicating this approach in other areas can allow for an assessment of density variation across the wolf range and would provide a reliable reference parameter for ecological studies.

We estimated the current size and dynamics of the wolf population in Tuscany and investigated the trends and demographic drivers of population changes. Estimates were obtained by two different approaches: (i) mixed-technique field monitoring (from 2014 to 2016) that found the minimum observed pack number and estimated population size, and (ii) an individual-based model (run by Vortex software v. 10.3.8.0) with demographic inputs derived from a local intensive study area and historic data on population size. Field monitoring showed a minimum population size of 558 wolves (SE = 12.005) in 2016, with a density of 2.74 individuals/100 km². The population model described an increasing trend with an average annual rate of increase λ = 1.075 (SE = 0.014), an estimated population size of about 882 individuals (SE = 9.397) in 2016, and a density of 4.29 wolves/100 km². Previously published estimates of wolf population were as low as 56.2% compared to our field monitoring estimation and 34.6% in comparison to our model estimation. We conducted sensitivity tests to analyze the key parameters driving population changes based on juvenile and adult mortality rates, female breeding success, and litter size. Mortality rates played a major role in determining intrinsic growth rate changes, with adult mortality accounting for 62.5% of the total variance explained by the four parameters. Juvenile mortality was responsible for 35.8% of the variance, while female breeding success and litter size had weak or negligible effects. We concluded that reliable estimates of population abundance and a deeper understanding of the role of different demographic parameters in determining population dynamics are crucial to define and carry out appropriate conservation and management strategies to address human–wildlife conflicts.

European wolf populations are expanding into human-dominated landscapes, triggering novel interactions with citizens and public concerns that may disrupt the traditional urban–rural divide in wolf attitudes and reshape conservation paradigms. We modelled the spatiotemporal distribution and valence of the wolf reports received through a dedicated phone service in Tuscany, Italy (2021–2024). Reports were significantly more common in: (i) late winter, aligning with the peak dispersal period and increased wolf movements. (ii) recently recolonized areas, suggesting a wolf-novelty effect, and (iii) urban areas, where negative valence was also more likely. Public concerns about wolves are increasingly emerging in urban areas, potentially disrupting the traditionally more supportive urban stance on wolf presence. In one of the European regions where wolf recovery began earlier and progressed further, our findings signal a broader shift in public attitudes that may weaken support for wolf conservation, potentially anticipating similar developments in areas of more recent recovery. Graphical Abstract