Integrating sign surveys and telemetry data for estimating brown bear (Ursus arctos) density in the Romanian Carpathians

Accurate population size estimates are important information for sustainable wildlife management. The Romanian Carpathians harbor the largest brown bear (Ursus arctos) population in Europe, yet current management relies on estimates of density that lack statistical oversight and ignore uncertainty deriving from track surveys. In this study, we investigate an alternative approach to estimate brown bear density using sign surveys along transects within a novel integration of occupancy models and home range methods. We performed repeated surveys along 2‐km segments of forest roads during three distinct seasons: spring 2011, fall‐winter 2011, and spring 2012, within three game management units and a Natura 2000 site. We estimated bears abundances along transects using the number of unique tracks observed per survey occasion via N‐mixture hierarchical models, which account for imperfect detection. To obtain brown bear densities, we combined these abundances with the effective sampling area of the transects, that is, estimated as a function of the median (± bootstrapped SE) of the core home range (5.58 ± 1.08 km(2)) based on telemetry data from 17 bears tracked for 1‐month periods overlapping our surveys windows. Our analyses yielded average brown bear densities (and 95% confidence intervals) for the three seasons of: 11.5 (7.8–15.3), 11.3 (7.4–15.2), and 12.4 (8.6–16.3) individuals/100 km(2). Across game management units, mean densities ranged between 7.5 and 14.8 individuals/100 km(2). Our method incorporates multiple sources of uncertainty (e.g., effective sampling area, imperfect detection) to estimate brown bear density, but the inference fundamentally relies on unmarked individuals only. While useful as a temporary approach to monitor brown bears, we urge implementing DNA capture–recapture methods regionally to inform brown bear management and recommend increasing resources for GPS collars to improve estimates of effective sampling area.