A multidisciplinary approach to estimating wolf population size for long-term conservation.

The wolf (Canis lupus) is among the most controversial of wildlife species. Abundance estimates are required to inform public debate and policy decisions, but obtaining them at biologically relevant scales is challenging. We developed a system for comprehensive population estimation across the Italian alpine region (100,000 km2 ), involving 1513 trained operators representing 160 institutions. This extensive network allowed for coordinated genetic sample collection and landscape-level spatial capture-recapture analyses that transcended administrative boundaries to produce the first estimates of key parameters for wolf population status assessment. Wolf abundance was estimated at 952 individuals (95% credible interval 816-1120) and 135 reproductive units (i.e., packs) (95% credible interval 112-165). We also estimated that mature individuals accounted for 33-45% of the entire population. The monitoring effort was spatially estimated thereby overcoming an important limitation of citizen science data. This is an important approach for promoting wolf-human coexistence based on wolf abundance monitoring and an endorsement of large-scale harmonized conservation practices.

Una estrategia multidisciplinaria para la estimación del tamaño poblacional de los lobos para la conservación a largo plazo Resumen El lobo (Canis lupus) está entre las especies de fauna más controversiales. Se requieren estimaciones de abundancia para informar al debate público y las decisiones políticas, pero es un reto obtenerlos en escalas con relevancia biológica. Desarrollamos un sistema para la estimación completa de la población en la región alpina de Italia (100,000 km2 ), con la participación de 1,513 operadores entrenados que representan a 160 instituciones. Esta red extensa permitió una colecta coordinada de muestras genéticas y análisis de captura-recaptura espacial que trascendieron las fronteras administrativas para así producir las primeras estimaciones de los parámetros clave para la evaluación del estado de la población de los lobos. Se estimó la abundancia en 952 individuos (95% intervalo de confianza 816-1120) y 135 unidades reproductivas (es decir, manadas) (95% intervalo de confianza 112-165). También estimamos que los individuos maduros representaban el 33-45% de toda la población. El esfuerzo de monitoreo se estimó espacialmente, por lo que sobrepasó una limitación importante de la ciencia ciudadana. Esta estrategia es importante para promover la coexistencia entre lobos y humanos con base en el monitoreo de la abundancia y el apoyo a las prácticas armonizadas de conservación a gran escala.

Estimating and forecasting spatial population dynamics of apex predators using transnational genetic monitoring.

The ongoing recovery of terrestrial large carnivores in North America and Europe is accompanied by intense controversy. On the one hand, reestablishment of large carnivores entails a recovery of their most important ecological role, predation. On the other hand, societies are struggling to relearn how to live with apex predators that kill livestock, compete for game species, and occasionally injure or kill people. Those responsible for managing these species and mitigating conflict often lack fundamental information due to a long-standing challenge in ecology: How do we draw robust population-level inferences for elusive animals spread over immense areas? Here we showcase the application of an effective tool for spatially explicit tracking and forecasting of wildlife population dynamics at scales that are relevant to management and conservation. We analyzed the world's largest dataset on carnivores comprising more than 35,000 noninvasively obtained DNA samples from over 6,000 individual brown bears (Ursus arctos), gray wolves (Canis lupus), and wolverines (Gulo gulo). Our analyses took into account that not all individuals are detected and, even if detected, their fates are not always known. We show unequivocal quantitative evidence of large carnivore recovery in northern Europe, juxtaposed with the finding that humans are the single-most important factor driving the dynamics of these apex predators. We present maps and forecasts of the spatiotemporal dynamics of large carnivore populations, transcending national boundaries and management regimes.

Occupancy winners in tropical protected forests: a pantropical analysis.

The structure of forest mammal communities appears surprisingly consistent across the continental tropics, presumably due to convergent evolution in similar environments. Whether such consistency extends to mammal occupancy, despite variation in species characteristics and context, remains unclear. Here we ask whether we can predict occupancy patterns and, if so, whether these relationships are consistent across biogeographic regions. Specifically, we assessed how mammal feeding guild, body mass and ecological specialization relate to occupancy in protected forests across the tropics. We used standardized camera-trap data (1002 camera-trap locations and 2-10 years of data) and a hierarchical Bayesian occupancy model. We found that occupancy varied by regions, and certain species characteristics explained much of this variation. Herbivores consistently had the highest occupancy. However, only in the Neotropics did we detect a significant effect of body mass on occupancy: large mammals had lowest occupancy. Importantly, habitat specialists generally had higher occupancy than generalists, though this was reversed in the Indo-Malayan sites. We conclude that habitat specialization is key for understanding variation in mammal occupancy across regions, and that habitat specialists often benefit more from protected areas, than do generalists. The contrasting examples seen in the Indo-Malayan region probably reflect distinct anthropogenic pressures.

Context-dependent fitness costs of reproduction despite stable body mass costs in an Arctic herbivore.

The cost of reproduction on demographic rates is often assumed to operate through changing body condition. Several studies have found that reproduction depresses body mass more if the current conditions are severe, such as high population densities or adverse weather, than under benign environmental conditions. However, few studies have investigated the association between the fitness components and body mass costs of reproduction. Using 25 years of individual-based capture-recapture data from Svalbard reindeer Rangifer tarandus platyrhynchus, we built a novel Bayesian state-space model that jointly estimated interannual change in mass, annual reproductive success and survival, while accounting for incomplete observations. The model allowed us to partition the differential effects of intrinsic and extrinsic factors on both non-reproductive mass change and the body mass cost of reproduction, and to quantify their consequences on demographic rates. Contrary to our expectation, the body mass cost of reproduction (mean = -5.8 kg) varied little between years (CV = 0.08), whereas the between-year variation in body mass changes, that were independent of the previous year's reproductive state, varied substantially (CV = 0.4) in relation to autumn temperature and the amount of rain-on-snow in winter. This body mass loss led to a cost of reproduction on the next reproduction, which was amplified by the same environmental covariates, from a 10% reduction in reproductive success in benign years, to a 50% reduction in harsh years. The reproductive mass loss also resulted in a small reduction in survival. Our results show how demographic costs of reproduction, driven by interannual fluctuations in individual body condition, result from the balance between body mass costs of reproduction and body mass changes that are independent of previous reproductive state. We illustrate how a strong context-dependent fitness cost of reproduction can occur, despite a relatively fixed body mass cost of reproduction. This suggests that female reindeer display a very conservative energy allocation strategy, either aborting their reproductive attempt at an early stage or weaning at a relatively constant cost. Such a strategy might be common in species living in a highly stochastic and food limited environment.

Estimating red fox density using non-invasive genetic sampling and spatial capture-recapture modelling.

Spatial capture-recapture modelling (SCR) is a powerful tool for estimating density, population size, and space use of elusive animals. Here, we applied SCR modelling to non-invasive genetic sampling (NGS) data to estimate red fox (Vulpes vulpes) densities in two areas of boreal forest in central (2016-2018) and southern Norway (2017-2018). Estimated densities were overall lower in the central study area (mean = 0.04 foxes per km2 in 2016, 0.10 in 2017, and 0.06 in 2018) compared to the southern study area (0.16 in 2017 and 0.09 in 2018). We found a positive effect of forest cover on density in the central, but not the southern study area. The absence of an effect in the southern area may reflect a paucity of evidence caused by low variation in forest cover. Estimated mean home-range size in the central study area was 45 km2 [95%CI 34-60] for females and 88 km2 [69-113] for males. Mean home-range sizes were smaller in the southern study area (26 km2 [16-42] for females and 56 km2 [35-91] for males). In both study areas, detection probability was session-dependent and affected by sampling effort. This study highlights how SCR modelling in combination with NGS can be used to efficiently monitor red fox populations, and simultaneously incorporate ecological factors and estimate their effects on population density and space use.

GPS collars have an apparent positive effect on the survival of a large carnivore.

Are instrumented animals representative of the population, given the potential bias caused by selective sampling and the influence of capture, handling and wearing bio-loggers? The answer is elusive owing to the challenges of obtaining comparable data from individuals with and without bio-loggers. Using non-invasive genetic data of a large carnivore, the wolverine (Gulo gulo) in Scandinavia, and an open-population spatial capture-recapture model, we found a 16 (credible interval: 4-30) percentage points lower mortality probability for GPS-collared individuals compared with individuals without GPS collars. While the risk of dying from legal culling was comparable for collared and non-collared wolverines, the former experienced lower probability of mortality due to causes other than legal culling. The aforementioned effect was pronounced despite a potentially lower age-and therefore likely higher natural mortality-of collared individuals. Reports of positive effects of bio-loggers on the survival of individuals are uncommon and we argue that GPS collars could shield animals from poaching. Our results highlight the challenges of drawing population-level inferences for populations subjected to poaching when using data from instrumented individuals.

Border Security Fencing and Wildlife: The End of the Transboundary Paradigm in Eurasia?

The ongoing refugee crisis in Europe has seen many countries rush to construct border security fencing to divert or control the flow of people. This follows a trend of border fence construction across Eurasia during the post-9/11 era. This development has gone largely unnoticed by conservation biologists during an era in which, ironically, transboundary cooperation has emerged as a conservation paradigm. These fences represent a major threat to wildlife because they can cause mortality, obstruct access to seasonally important resources, and reduce effective population size. We summarise the extent of the issue and propose concrete mitigation measures.

Sociodemographic factors modulate the spatial response of brown bears to vacancies created by hunting.

There is a growing recognition of the importance of indirect effects from hunting on wildlife populations, e.g. social and behavioural changes due to harvest, which occur after the initial offtake. Nonetheless, little is known about how the removal of members of a population influences the spatial configuration of the survivors. We studied how surviving brown bears (Ursus arctos) used former home ranges that had belonged to casualties of the annual bear hunting season in southcentral Sweden (2007-2015). We used resource selection functions to explore the effects of the casualty's and survivor's sex, age and their pairwise genetic relatedness, population density and hunting intensity on survivors' spatial responses to vacated home ranges. We tested the competitive release hypothesis, whereby survivors that increase their use of a killed bear's home range are presumed to have been released from intraspecific competition. We found strong support for this hypothesis, as survivors of the same sex as the casualty consistently increased their use of its vacant home range. Patterns were less pronounced or absent when the survivor and casualty were of opposite sex. Genetic relatedness between the survivor and the casualty emerged as the most important factor explaining increased use of vacated male home ranges by males, with a stronger response from survivors of lower relatedness. Relatedness was also important for females, but it did not influence use following removal; female survivors used home ranges of higher related female casualties more, both before and after death. Spatial responses by survivors were further influenced by bear age, population density and hunting intensity. We have shown that survivors exhibit a spatial response to vacated home ranges caused by hunting casualties, even in nonterritorial species such as the brown bear. This spatial reorganization can have unintended consequences for population dynamics and interfere with management goals. Altogether, our results underscore the need to better understand the short- and long-term indirect effects of hunting on animal social structure and their resulting distribution in space.

Evolutionary history of enigmatic bears in the Tibetan Plateau-Himalaya region and the identity of the yeti.

Although anecdotally associated with local bears (Ursus arctos and U. thibetanus), the exact identity of 'hominid'-like creatures important to folklore and mythology in the Tibetan Plateau-Himalaya region is still surrounded by mystery. Recently, two purported yeti samples from the Himalayas showed genetic affinity with an ancient polar bear, suggesting they may be from previously unrecognized, possibly hybrid, bear species, but this preliminary finding has been under question. We conducted a comprehensive genetic survey of field-collected and museum specimens to explore their identity and ultimately infer the evolutionary history of bears in the region. Phylogenetic analyses of mitochondrial DNA sequences determined clade affinities of the purported yeti samples in this study, strongly supporting the biological basis of the yeti legend to be local, extant bears. Complete mitochondrial genomes were assembled for Himalayan brown bear (U. a. isabellinus) and black bear (U. t. laniger) for the first time. Our results demonstrate that the Himalayan brown bear is one of the first-branching clades within the brown bear lineage, while Tibetan brown bears diverged much later. The estimated times of divergence of the Tibetan Plateau and Himalayan bear lineages overlap with Middle to Late Pleistocene glaciation events, suggesting that extant bears in the region are likely descendants of populations that survived in local refugia during the Pleistocene glaciations.

Leave before it's too late: anthropogenic and environmental triggers of autumn migration in a hunted ungulate population.

Autumn has to a large extent been neglected in the climate effect literature, yet autumn events, e.g., plant senescence and animal migration, affect fitness of animals differently than spring events. Understanding how variables including plant phenology influence timing of autumn migrations is important to gain a comprehensive understanding of the full annual cycle of migratory species. Here we use 13 yr of data from 60 male and 168 female red deer (Cervus elaphus) to identify triggers of autumn migration. We relate the timing of autumn migration to environmental variables like snow fall, temperature, and plant phenology (NDVI), and to onset of hunting, sex, and migration distance. Severe weather has been suggested as the main trigger of autumn migration, but we found that the majority of the individuals had left the summer range well before snow fall (80.3%) and frost (70.5%), and also before the peak deterioration in forage quality (71.9%). Declining temperatures were associated with a higher daily migration potential. Onset of hunting showed the largest effect on migration potential, with a marked increase during the first days of hunting. Individuals still present in the summer range when snow fall, frost, or peak forage deterioration occurred showed a significantly higher migration potential around these events. Males were less responsive to environmental cues, suggesting rutting activity, starting earlier in males, initiate movement prior to such conditions. Also, individuals with longer migration distances had a higher migration potential late in the season than individuals with shorter migration distances. Our study shows that factors beyond weather and plant phenology, such as onset of hunting, may be important triggers of autumn migration. Severe weather and forage deterioration were important triggers for the individuals experiencing this, which suggests a hierarchical response to environmental cues. The trade-off between staying longer in the summer range and increased energy expenditures if surprised by severe weather is asymmetric, and leaving well in advance can be seen as a risk-averse tactic.

Leave before it's too late: anthropogenic and environmental triggers of autumn migration in a hunted ungulate population.

Autumn has to a large extent been neglected in the climate effect literature, yet autumn events, e.g., plant senescence and animal migration, affect fitness of animals differently than spring events. Understanding how variables including plant phenology influence timing of autumn migrations is important to gain a comprehensive understanding of the full annual cycle of migratory species. Here we use 13 yr of data from 60 male and 168 female red deer (Cervus elaphus) to identify triggers of autumn migration. We relate the timing of autumn migration to environmental variables like snow fall, temperature, and plant phenology (NDVI), and to onset of hunting, sex, and migration distance. Severe weather has been suggested as the main trigger of autumn migration, but we found that the majority of the individuals had left the summer range well before snow fall (80.3%) and frost (70.5%), and also before the peak deterioration in forage quality (71.9%). Declining temperatures were associated with a higher daily migration potential. Onset of hunting showed the largest effect on migration potential, with a marked increase during the first days of hunting. Individuals still present in the summer range when snow fall, frost, or peak forage deterioration occurred showed a significantly higher migration potential around these events. Males were less responsive to environmental cues, suggesting rutting activity, starting earlier in males, initiate movement prior to such conditions. Also, individuals with longer migration distances had a higher migration potential late in the season than individuals with shorter migration distances. Our study shows that factors beyond weather and plant phenology, such as onset of hunting, may be important triggers of autumn migration. Severe weather and forage deterioration were important triggers for the individuals experiencing this, which suggests a hierarchical response to environmental cues. The trade-off between staying longer in the summer range and increased energy expenditures if surprised by severe weather is asymmetric, and leaving well in advance can be seen as a risk-averse tactic.

Determinants of lifetime reproduction in female brown bears: early body mass, longevity, and hunting regulations.

In iteroparous mammals, conditions experienced early in life may have long-lasting effects on lifetime reproductive success. Human-induced mortality is also an important demographic factor in many populations of large mammals and may influence lifetime reproductive success. Here, we explore the effects of early development, population density, and human hunting on survival and lifetime reproductive success in brown bear (Ursus arctos) females, using a 25-year database of individually marked bears in two populations in Sweden. Survival of yearlings to 2 years was not affected by population density or body mass. Yearlings that remained with their mother had higher survival than independent yearlings, partly because regulations prohibit the harvest of bears in family groups. Although mass as a yearling did not affect juvenile survival, it was positively associated with measures of lifetime reproductive success and individual fitness. The majority of adult female brown bear mortality (72%) in our study was due to human causes, mainly hunting, and many females were killed before they reproduced. Therefore, factors allowing females to survive several hunting seasons had a strong positive effect on lifetime reproductive success. We suggest that, in many hunted populations of large mammals, sport harvest is an important influence on both population dynamics and life histories.

Estimating spatially variable and density-dependent survival using open-population spatial capture-recapture models.

Open-population spatial capture-recapture (OPSCR) models use the spatial information contained in individual detections collected over multiple consecutive occasions to estimate not only occasion-specific density, but also demographic parameters. OPSCR models can also estimate spatial variation in vital rates, but such models are neither widely used nor thoroughly tested. We developed a Bayesian OPSCR model that not only accounts for spatial variation in survival using spatial covariates but also estimates local density-dependent effects on survival within a unified framework. Using simulations, we show that OPSCR models provide sound inferences on the effect of spatial covariates on survival, including multiple competing sources of mortality, each with potentially different spatial determinants. Estimation of local density-dependent survival was possible but required more data due to the greater complexity of the model. Not accounting for spatial heterogeneity in survival led to up to 10% positive bias in abundance estimates. We provide an empirical demonstration of the model by estimating the effect of country and density on cause-specific mortality of female wolverines (Gulo gulo) in central Sweden and Norway. The ability to make population-level inferences on spatial variation in survival is an essential step toward a fully spatially explicit OPSCR model capable of disentangling the role of multiple spatial drivers of population dynamics.

Sexual segregation results in pronounced sex-specific density gradients in the mountain ungulate, Rupicapra rupicapra.

Sex-specific differences in habitat selection and space use are common in ungulates. Yet, it is largely unknown how this behavioral dimorphism, ultimately leading to sexual segregation, translates to population-level patterns and density gradients across landscapes. Alpine chamois (Rupicapra rupicapra r.) predominantly occupy habitat above tree line, yet especially males may also take advantage of forested habitats. To estimate male and female chamois density and determinants thereof, we applied Bayesian spatial capture-recapture (SCR) models in two contrasting study areas in the Alps, Germany, during autumn. We fitted SCR models to non-invasive individual encounter data derived from genotyped feces. Sex-specific densities were modeled as a function of terrain ruggedness, forest canopy cover, proportion of barren ground, and site severity. We detected pronounced differences in male and female density patterns, driven primarily by terrain ruggedness, rather than by sex-specific effects of canopy cover. The positive effect of ruggedness on density was weaker for males which translated into a higher proportion of males occupying less variable terrain, frequently located in forests, compared to females. By estimating sex-specific variation in both detection probabilities and density, we were able to quantify and map how individual behavioral differences scale up and shape spatial patterns in population density.

A flexible and efficient Bayesian implementation of point process models for spatial capture-recapture data.

Spatial capture-recapture (SCR) is now routinely used for estimating abundance and density of wildlife populations. A standard SCR model includes sub-models for the distribution of individual activity centers (ACs) and for individual detections conditional on the locations of these ACs. Both sub-models can be expressed as point processes taking place in continuous space, but there is a lack of accessible and efficient tools to fit such models in a Bayesian paradigm. Here, we describe a set of custom functions and distributions to achieve this. Our work allows for more efficient model fitting with spatial covariates on population density, offers the option to fit SCR models using the semi-complete data likelihood (SCDL) approach instead of data augmentation, and better reflects the spatially continuous detection process in SCR studies that use area searches. In addition, the SCDL approach is more efficient than data augmentation for simple SCR models while losing its advantages for more complicated models that account for spatial variation in either population density or detection. We present the model formulation, test it with simulations, quantify computational efficiency gains, and conclude with a real-life example using non-invasive genetic sampling data for an elusive large carnivore, the wolverine (Gulo gulo) in Norway.

Spatial variation in red deer density in a transboundary forest ecosystem.

Forests in Europe are exposed to increasingly frequent and severe disturbances. The resulting changes in the structure and composition of forests can have profound consequences for the wildlife inhabiting them. Moreover, wildlife populations in Europe are often subjected to differential management regimes as they regularly extend across multiple national and administrative borders. The red deer Cervus elaphus population in the Bohemian Forest Ecosystem, straddling the Czech-German border, has experienced forest disturbances, primarily caused by windfalls and bark beetle Ips typographus outbreaks during the past decades. To adapt local management strategies to the changing environmental conditions and to coordinate them across the international border, reliable estimates of red deer density and abundance are highly sought-after by policymakers, wildlife managers, and stakeholders. Covering a 1081-km2 study area, we conducted a transnational non-invasive DNA sampling study in 2018 that yielded 1578 genotyped DNA samples from 1120 individual red deer. Using spatial capture-recapture models, we estimated total and jurisdiction-specific abundance of red deer throughout the ecosystem and quantified the role of forest disturbance and differential management strategies in shaping spatial heterogeneity in red deer density. We hypothesised that (a) forest disturbances provide favourable habitat conditions (e.g., forage and cover), and (b) contrasting red deer management regimes in different jurisdictions create a differential risk landscape, ultimately shaping density distributions. Overall, we estimated that 2851 red deer (95% Credible Interval = 2609-3119) resided in the study area during the sampling period, with a relatively even overall sex ratio (1406 females, 95% CI = 1229-1612 and 1445 males, 95% CI = 1288-1626). The average red deer density was higher in Czechia (3.5 km-2, 95% CI = 1.2-12.3) compared to Germany (2 km-2, 95% CI = 0.2-11). The effect of forest disturbances on red deer density was context-dependent. Forest disturbances had a positive effect on red deer density at higher elevations and a negative effect at lower elevations, which could be explained by partial migration and its drivers in this population. Density of red deer was generally higher in management units where hunting is prohibited. In addition, we found that sex ratios differed between administrative units and were more balanced in the non-intervention zones. Our results show that the effect of forest disturbances on wild ungulates is modulated by additional factors, such as elevation and ungulate management practices. Overall density patterns and sex ratios suggested strong gradients in density between administrative units. With climate change increasing the severity and frequency of forest disturbances, population-level monitoring and management are becoming increasingly important, especially for wide-ranging species as both wildlife and global change transcend administrative boundaries.

Occurrence dynamics of mammals in protected tropical forests respond to human presence and activities.

Protected areas (PAs) play a vital role in wildlife conservation. Nonetheless there is concern and uncertainty regarding how and at what spatial scales anthropogenic stressors influence the occurrence dynamics of wildlife populations inside PAs. Here we assessed how anthropogenic stressors influence occurrence dynamics of 159 mammal species in 16 tropical PAs from three biogeographic regions. We quantified these relationships for species groups (habitat specialists and generalists) and individual species. We used long-term camera-trap data (1,002 sites) and fitted Bayesian dynamic multispecies occupancy models to estimate local colonization (the probability that a previously empty site is colonized) and local survival (the probability that an occupied site remains occupied). Multiple covariates at both the local scale and landscape scale influenced mammal occurrence dynamics, although responses differed among species groups. Colonization by specialists increased with local-scale forest cover when landscape-scale fragmentation was low. Survival probability of generalists was higher near the edge than in the core of the PA when landscape-scale human population density was low but the opposite occurred when population density was high. We conclude that mammal occurrence dynamics are impacted by anthropogenic stressors acting at multiple scales including outside the PA itself.

A migratory northern ungulate in the pursuit of spring: jumping or surfing the green wave?

The forage-maturation hypothesis (FMH) states that herbivores migrate along a phenological gradient of plant development in order to maximize energy intake. Despite strong support for the FMH, the actual relationship between plant phenology and ungulate movement has remained enigmatic. We linked plant phenology (MODIS-normalized difference vegetation index [NDVI] data) and space use of 167 migratory and 78 resident red deer (Cervus elaphus), using a space-time-time matrix of "springness," defined as the instantaneous rate of green-up. Consistent with the FMH, migrants experienced substantially greater access to early plant phenology than did residents. Deer were also more likely to migrate in areas where migration led to greater gains in springness. Rather than "surfing the green wave" during migration, migratory red deer moved rapidly from the winter to the summer range, thereby "jumping the green wave." However, migrants and, to a lesser degree, residents did track phenological green-up through parts of the growing season by making smaller-scale adjustments in habitat use. Despite pronounced differences in their life histories, we found only marginal differences between male and female red deer in this study. Those differences that we did detect pointed toward additional constraints on female space-use tactics, such as those posed by calving and caring for dependent offspring. We conclude that whereas in some systems migration itself is a way to surf the green wave, in others it may simply be a means to reconnect with phenological spring at the summer range. In the light of ubiquitous anthropogenic environmental change, understanding the relationship between the green wave and ungulate space use has important consequences for the management and conservation of migratory ungulates and the phenomenon of migration itself.

Linking noninvasive genetic sampling and traditional monitoring to aid management of a trans-border carnivore population.

Noninvasive genetic sampling has been embraced by wildlife managers and ecologists, especially those charged with monitoring rare and elusive species over large areas. Challenges arise when desired population measures are not directly attainable from genetic data and when monitoring targets trans-border populations. Norwegian management authorities count individual brown bears (Ursus arctos) using noninvasive genetic sampling but express management goals in the annual number of bear reproductions (females that produce cubs), a measure that is not directly available from genetic data. We combine noninvasive genetic sampling data with information obtained from a long-term intensive monitoring study in neighboring Sweden to estimate the number of annual reproductions by females detected within Norway. Most female brown bears in Norway occur near the border with neighboring countries (Sweden, Finland, and Russia) and their potential reproduction can therefore only partially be credited to Norway. Our model includes a simulation-based method that corrects census data to account for this. We estimated that 4.3 and 5.7 reproductions can be credited to females detected with noninvasive genetic sampling in Norway in 2008 and 2009, respectively. These numbers fall substantially short of the national target (15 annual reproductions). Ignoring the potential for home ranges to extend beyond Norway's borders leads to an increase in the estimate of the number of reproductions by -30%. Our study shows that combining noninvasive genetic sampling with information obtained from traditional intensive/invasive monitoring can help answer contemporary management questions in the currency desired by managers and policy makers. Furthermore, combining methodologies and thereby accounting for space use increases the accuracy of the information on which decisions are based. It is important that the information derived from multiple approaches is applicable to the same focal population and that predictions are cross-validated. When monitoring and management are constrained to administrative units, census data should be adjusted by discounting portions of individual space utilization that extend beyond the focal jurisdiction. Our simulation-based approach for making such an adjustment may be useful in other situations where management authorities target portions of trans-border populations.

Does the punishment fit the crime? Consequences and diagnosis of misspecified detection functions in Bayesian spatial capture-recapture modeling.

Spatial capture-recapture (SCR) analysis is now used routinely to inform wildlife management and conservation decisions. It is therefore imperative that we understand the implications of and can diagnose common SCR model misspecifications, as flawed inferences could propagate to policy and interventions. The detection function of an SCR model describes how an individual's detections are distributed in space. Despite the detection function's central role in SCR, little is known about the robustness of SCR-derived abundance estimates and home range size estimates to misspecifications. Here, we set out to (a) determine whether abundance estimates are robust to a wider range of misspecifications of the detection function than previously explored, (b) quantify the sensitivity of home range size estimates to the choice of detection function, and (c) evaluate commonly used Bayesian p-values for detecting misspecifications thereof. We simulated SCR data using different circular detection functions to emulate a wide range of space use patterns. We then fit Bayesian SCR models with three detection functions (half-normal, exponential, and half-normal plateau) to each simulated data set. While abundance estimates were very robust, estimates of home range size were sensitive to misspecifications of the detection function. When misspecified, SCR models with the half-normal plateau and exponential detection functions produced the most and least reliable home range size, respectively. Misspecifications with the strongest impact on parameter estimates were easily detected by Bayesian p-values. Practitioners using SCR exclusively for density estimation are unlikely to be impacted by misspecifications of the detection function. However, the choice of detection function can have substantial consequences for the reliability of inferences about space use. Although Bayesian p-values can aid the diagnosis of detection function misspecification under certain conditions, we urge the development of additional custom goodness-of-fit diagnostics for Bayesian SCR models to identify a wider range of model misspecifications.