Evaluating density-weighted connectivity of black bears (Ursus americanus) in Glacier National Park with spatial capture-recapture models.

BACKGROUND: Improved understanding of wildlife population connectivity among protected area networks can support effective planning for the persistence of wildlife populations in the face of land use and climate change. Common approaches to estimating connectivity often rely on small samples of individuals without considering the spatial structure of populations, leading to limited understanding of how individual movement links to demography and population connectivity. Recently developed spatial capture-recapture (SCR) models provide a framework to formally connect inference about individual movement, connectivity, and population density, but few studies have applied this approach to empirical data to support connectivity planning. METHODS: We used mark-recapture data collected from 924 genetic detections of 598 American black bears (Ursus americanus) in 2004 with SCR ecological distance models to simultaneously estimate density, landscape resistance to movement, and population connectivity in Glacier National Park northwest Montana, USA. We estimated density and movement parameters separately for males and females and used model estimates to calculate predicted density-weighted connectivity surfaces. RESULTS: Model results indicated that landscape structure influences black bear density and space use in Glacier. The mean density estimate was 16.08 bears/100 km2 (95% CI 12.52-20.6) for females and 9.27 bears/100 km2 (95% CI 7.70-11.14) for males. Density increased with forest cover for both sexes. For male black bears, density decreased at higher grizzly bear (Ursus arctos) densities. Drainages, valley bottoms, and riparian vegetation decreased estimates of landscape resistance to movement for male and female bears. For males, forest cover also decreased estimated resistance to movement, but a transportation corridor bisecting the study area strongly increased resistance to movement presenting a barrier to connectivity. CONCLUSIONS: Density-weighed connectivity surfaces highlighted areas important for population connectivity that were distinct from areas with high potential connectivity. For black bears in Glacier and surrounding landscapes, consideration of both vegetation and valley topography could inform the placement of underpasses along the transportation corridor in areas characterized by both high population density and potential connectivity. Our study demonstrates that the SCR ecological distance model can provide biologically realistic, spatially explicit predictions to support movement connectivity planning across large landscapes.

Bighorn sheep associations: understanding tradeoffs of sociality and implications for disease transmission.

Sociality directly influences mating success, survival rates, and disease, but ultimately likely evolved for its fitness benefits in a challenging environment. The tradeoffs between the costs and benefits of sociality can operate at multiple scales, resulting in different interpretations of animal behavior. We investigated the influence of intrinsic (e.g., relatedness, age) and extrinsic factors (e.g., land cover type, season) on direct contact (simultaneous GPS locations ≤ 25 m) rates of bighorn sheep (Ovis canadensis) at multiple scales near the Waterton-Glacier International Peace Park. During 2002-2012, male and female bighorn were equipped with GPS collars. Indirect contact (GPS locations ≤ 25 m regardless of time) networks identified two major breaks whereas direct contact networks identified an additional barrier in the population, all of which corresponded with prior disease exposure metrics. More direct contacts occurred between same-sex dyads than female-male dyads and between bighorn groups with overlapping summer home ranges. Direct contacts occurred most often during the winter-spring season when bighorn traveled at low speeds and when an adequate number of bighorn were collared in the area. Direct contact probabilities for all dyad types were inversely related to habitat quality, and differences in contact probability were driven by variables related to survival such as terrain ruggedness, distance to escape terrain, and canopy cover. We provide evidence that probabilities of association are higher when there is greater predation risk and that contact analysis provides valuable information for understanding fitness tradeoffs of sociality and disease transmission potential.

Recent and future declines of a historically widespread pollinator linked to climate, land cover, and pesticides.

The acute decline in global biodiversity includes not only the loss of rare species, but also the rapid collapse of common species across many different taxa. The loss of pollinating insects is of particular concern because of the ecological and economic values these species provide. The western bumble bee (Bombus occidentalis) was once common in western North America, but this species has become increasingly rare through much of its range. To understand potential mechanisms driving these declines, we used Bayesian occupancy models to investigate the effects of climate and land cover from 1998 to 2020, pesticide use from 2008 to 2014, and projected expected occupancy under three future scenarios. Using 14,457 surveys across 2.8 million km2 in the western United States, we found strong negative relationships between increasing temperature and drought on occupancy and identified neonicotinoids as the pesticides of greatest negative influence across our study region. The mean predicted occupancy declined by 57% from 1998 to 2020, ranging from 15 to 83% declines across 16 ecoregions. Even under the most optimistic scenario, we found continued declines in nearly half of the ecoregions by the 2050s and mean declines of 93% under the most severe scenario across all ecoregions. This assessment underscores the tenuous future of B. occidentalis and demonstrates the scale of stressors likely contributing to rapid loss of related pollinator species throughout the globe. Scaled-up, international species-monitoring schemes and improved integration of data from formal surveys and community science will substantively improve the understanding of stressors and bumble bee population trends.

Precision and bias of spatial capture-recapture estimates: A multi-site, multi-year Utah black bear case study.

Spatial capture-recapture (SCR) models are powerful analytical tools that have become the standard for estimating abundance and density of wild animal populations. When sampling populations to implement SCR, the number of unique individuals detected, total recaptures, and unique spatial relocations can be highly variable. These sample sizes influence the precision and accuracy of model parameter estimates. Testing the performance of SCR models with sparse empirical data sets typical of low-density, wide-ranging species can inform the threshold at which a more integrated modeling approach with additional data sources or additional years of monitoring may be required to achieve reliable, precise parameter estimates. Using a multi-site, multi-year Utah black bear (Ursus americanus) capture-recapture data set, we evaluated factors influencing the uncertainty of SCR structural parameter estimates, specifically density, detection, and the spatial scale parameter, sigma. We also provided some of the first SCR density estimates for Utah black bear populations, which ranged from 3.85 to 74.33 bears/100 km2 . Increasing total detections decreased the uncertainty of density estimates, whereas an increasing number of total recaptures and individuals with recaptures decreased the uncertainty of detection and sigma estimates, respectively. In most cases, multiple years of data were required for precise density estimates (<0.2 coefficient of variation [CV]). Across study areas there was an average decline in CV of 0.07 with the addition of another year of data. One sampled population with very high estimated bear density had an atypically low number of spatial recaptures relative to total recaptures, apparently inflating density estimates. A complementary simulation study used to assess estimate bias suggested that when <30% of recaptured individuals were spatially recaptured, density estimates were unreliable and ranged widely, in some cases to >3 times the simulated density. Additional research could evaluate these requirements for other density scenarios. Large numbers of individuals detected, numbers of spatial recaptures, and precision alone may not be sufficient indicators of parameter estimate reliability. We provide an evaluation of simple summary statistics of capture-recapture data sets that can provide an early signal of the need to alter sampling design or collect auxiliary data before model implementation to improve estimate precision and accuracy.

Eyes on the herd: Quantifying ungulate density from satellite, unmanned aerial systems, and GPScollar data.

Novel approaches to quantifying density and distributions could help biologists adaptively manage wildlife populations, particularly if methods are accurate, consistent, cost-effective, rapid, and sensitive to change. Such approaches may also improve research on interactions between density and processes of interest, such as disease transmission across multiple populations. We assess how satellite imagery, unmanned aerial system (UAS) imagery, and Global Positioning System (GPS) collar data vary in characterizing elk density, distribution, and count patterns across times with and without supplemental feeding at the National Elk Refuge (NER) in the US state of Wyoming. We also present the first comparison of satellite imagery data with traditional counts for ungulates in a temperate system. We further evaluate seven different aggregation metrics to identify the most consistent and sensitive metrics for comparing density and distribution across time and populations. All three data sources detected higher densities and aggregation locations of elk during supplemental feeding than non-feeding at the NER. Kernel density estimates (KDEs), KDE polygon areas, and the first quantile of interelk distances detected differences with the highest sensitivity and were most highly correlated across data sources. Both UAS and satellite imagery provide snapshots of density and distribution patterns of most animals in the area at lower cost than GPS collars. While satellite-based counts were lower than traditional counts, aggregation metrics matched those from UAS and GPS data sources when animals appeared in high contrast to the landscape, including brown elk against new snow in open areas. UAS counts of elk were similar to traditional ground-based counts on feed grounds and are the best data source for assessing changes in small spatial extents. Satellite, UAS, or GPS data can provide appropriate data for assessing density and changes in density from adaptive management actions. For the NER, where high elk densities are beneath controlled airspace, GPS collar data will be most useful for evaluating how management actions, including changes in the dates of supplemental feeding, influence elk density and aggregation across large spatial extents. Using consistent and sensitive measures of density may improve research on the drivers and effects of density within and across a wide range of species.

The smell of success: Reproductive success related to rub behavior in brown bears.

Several species of bears are known to rub deliberately against trees and other objects, but little is known about why bears rub. Patterns in rubbing behavior of male and female brown bears (Ursus arctos) suggest that scent marking via rubbing functions to communicate among potential mates or competitors. Using DNA from bear hairs collected from rub objects in southwestern Alberta from 2011-2014 and existing DNA datasets from Montana and southeastern British Columbia, we determined sex and individual identity of each bear detected. Using these data, we completed a parentage analysis. From the parentage analysis and detection data, we determined the number of offspring, mates, unique rub objects where an individual was detected, and sampling occasions during which an individual was detected for each brown bear identified through our sampling methods. Using a Poisson regression, we found a positive relationship between bear rubbing behavior and reproductive success; both male and female bears with a greater number of mates and a greater number of offspring were detected at more rub objects and during more occasions. Our results suggest a fitness component to bear rubbing, indicate that rubbing is adaptive, and provide insight into a poorly understood behaviour.

Evaluating wildlife translocations using genomics: A bighorn sheep case study.

Wildlife restoration often involves translocation efforts to reintroduce species and supplement small, fragmented populations. We examined the genomic consequences of bighorn sheep (Ovis canadensis) translocations and population isolation to enhance understanding of evolutionary processes that affect population genetics and inform future restoration strategies. We conducted a population genomic analysis of 511 bighorn sheep from 17 areas, including native and reintroduced populations that received 0-10 translocations. Using the Illumina High Density Ovine array, we generated datasets of 6,155 to 33,289 single nucleotide polymorphisms and completed clustering, population tree, and kinship analyses. Our analyses determined that natural gene flow did not occur between most populations, including two pairs of native herds that had past connectivity. We synthesized genomic evidence across analyses to evaluate 24 different translocation events and detected eight successful reintroductions (i.e., lack of signal for recolonization from nearby populations) and five successful augmentations (i.e., reproductive success of translocated individuals) based on genetic similarity with the source populations. A single native population founded six of the reintroduced herds, suggesting that environmental conditions did not need to match for populations to persist following reintroduction. Augmentations consisting of 18-57 animals including males and females succeeded, whereas augmentations of two males did not result in a detectable genetic signature. Our results provide insight on genomic distinctiveness of native and reintroduced herds, information on the relative success of reintroduction and augmentation efforts and their associated attributes, and guidance to enhance genetic contribution of augmentations and reintroductions to aid in bighorn sheep restoration.

Using bear rub data and spatial capture-recapture models to estimate trend in a brown bear population.

Trends in population abundance can be challenging to quantify during range expansion and contraction, when there is spatial variation in trend, or the conservation area is large. We used genetic detection data from natural bear rubbing sites and spatial capture-recapture (SCR) modeling to estimate local density and population growth rates in a grizzly bear population in northwestern Montana, USA. We visited bear rubs to collect hair in 2004, 2009-2012 (3,579-4,802 rubs) and detected 249-355 individual bears each year. We estimated the finite annual population rate of change 2004-2012 was 1.043 (95% CI = 1.017-1.069). Population density shifted from being concentrated in the north in 2004 to a more even distribution across the ecosystem by 2012. Our genetic detection sampling approach coupled with SCR modeling allowed us to estimate spatially variable growth rates of an expanding grizzly bear population and provided insight into how those patterns developed. The ability of SCR to utilize unstructured data and produce spatially explicit maps that indicate where population change is occurring promises to facilitate the monitoring of difficult-to-study species across large spatial areas.

Circuit-theory applications to connectivity science and conservation.

Conservation practitioners have long recognized ecological connectivity as a global priority for preserving biodiversity and ecosystem function. In the early years of conservation science, ecologists extended principles of island biogeography to assess connectivity based on source patch proximity and other metrics derived from binary maps of habitat. From 2006 to 2008, the late Brad McRae introduced circuit theory as an alternative approach to model gene flow and the dispersal or movement routes of organisms. He posited concepts and metrics from electrical circuit theory as a robust way to quantify movement across multiple possible paths in a landscape, not just a single least-cost path or corridor. Circuit theory offers many theoretical, conceptual, and practical linkages to conservation science. We reviewed 459 recent studies citing circuit theory or the open-source software Circuitscape. We focused on applications of circuit theory to the science and practice of connectivity conservation, including topics in landscape and population genetics, movement and dispersal paths of organisms, anthropogenic barriers to connectivity, fire behavior, water flow, and ecosystem services. Circuit theory is likely to have an effect on conservation science and practitioners through improved insights into landscape dynamics, animal movement, and habitat-use studies and through the development of new software tools for data analysis and visualization. The influence of circuit theory on conservation comes from the theoretical basis and elegance of the approach and the powerful collaborations and active user community that have emerged. Circuit theory provides a springboard for ecological understanding and will remain an important conservation tool for researchers and practitioners around the globe.

Aplicaciones de la Teoría de Circuitos a la Conservación y a la Ciencia de la Conectividad Resumen Quienes practican la conservación han reconocido durante mucho tiempo que la conectividad ecológica es una prioridad mundial para la preservación de la biodiversidad y el funcionamiento del ecosistema. Durante los primeros años de la ciencia de la conservación los ecólogos difundieron los principios de la biografía de islas para evaluar la conectividad con base en la proximidad entre el origen y el fragmento, así como otras medidas derivadas de los mapas binarios de los hábitats. Entre 2006 y 2008 el fallecido Brad McRae introdujo la teoría de circuitos como una estrategia alternativa para modelar el flujo génico y la dispersión o las rutas de movimiento de los organismos. McRae propuso conceptos y medidas de la teoría de circuitos eléctricos como una manera robusta para cuantificar el movimiento a lo largo de múltiples caminos posibles en un paisaje, no solamente a lo largo de un camino o corredor de menor costo. La teoría de circuitos ofrece muchos enlaces teóricos, conceptuales y prácticos con la ciencia de la conservación. Revisamos 459 estudios recientes que citan la teoría de circuitos o el software de fuente abierta Circuitscape. Nos enfocamos en las aplicaciones de la teoría de circuitos a la ciencia y a la práctica de la conservación de la conectividad, incluyendo temas como la genética poblacional y del paisaje, movimiento y caminos de dispersión de los organismos, barreras antropogénicas de la conectividad, comportamiento ante incendios, flujo del agua, y servicios ambientales. La teoría de circuitos probablemente tenga un efecto sobre la ciencia de la conservación y quienes la practican por medio de una percepción mejorada de las dinámicas del paisaje, el movimiento animal, y los estudios de uso de hábitat, y por medio del desarrollo de nuevas herramientas de software para el análisis de datos y su visualización. La influencia de la teoría de circuitos sobre la conservación viene de la base teórica y la elegancia de la estrategia y de las colaboraciones fuertes y la comunidad activa de usuarios que han surgido recientemente. La teoría de circuitos proporciona un trampolín para el entendimiento ecológico y seguirá siendo una importante herramienta de conservación para los investigadores y practicantes en todo el mundo.

Relationship between effective and demographic population size in continuously distributed populations.

Genetic monitoring of wild populations can offer insights into demographic and genetic information simultaneously. However, widespread application of genetic monitoring is hindered by large uncertainty in the estimation and interpretation of target metrics such as contemporary effective population size, Ne . We used four long-term genetic and demographic studies (≥9 years) to evaluate the temporal stability of the relationship between Ne and demographic population size (Nc ). These case studies focused on mammals that are continuously distributed, yet dispersal-limited within the spatial scale of the study. We estimated local, contemporary Ne with single-sample methods (LDNE, Heterozygosity Excess, and Molecular Ancestry) and demographic abundance with either mark-recapture estimates or catch-per-unit effort indices. Estimates of Ne varied widely within each case study suggesting interpretation of estimates is challenging. We found inconsistent correlations and trends both among estimates of Ne and between Ne and Nc suggesting the value of Ne as an indicator of Nc is limited in some cases. In the two case studies with consistent trends between Ne and Nc , FIS was more stable over time and lower, suggesting FIS may be a good indicator that the population was sampled at a spatial scale at which genetic structure is not biasing estimates of Ne . These results suggest that more empirical work on the estimation of Ne in continuous populations is needed to understand the appropriate context to use LDNe as a useful metric in a monitoring programme to detect temporal trends in either Ne or Nc .

Evaluating sample size to estimate genetic management metrics in the genomics era.

Inbreeding and relationship metrics among and within populations are useful measures for genetic management of wild populations, but accuracy and precision of estimates can be influenced by the number of individual genotypes analysed. Biologists are confronted with varied advice regarding the sample size necessary for reliable estimates when using genomic tools. We developed a simulation framework to identify the optimal sample size for three widely used metrics to enable quantification of expected variance and relative bias of estimates and a comparison of results among populations. We applied this approach to analyse empirical genomic data for 30 individuals from each of four different free-ranging Rocky Mountain bighorn sheep (Ovis canadensis canadensis) populations in Montana and Wyoming, USA, through cross-species application of an Ovine array and analysis of approximately 14,000 single nucleotide polymorphisms (SNPs) after filtering. We examined intra- and interpopulation relationships using kinship and identity by state metrics, as well as FST between populations. By evaluating our simulation results, we concluded that a sample size of 25 was adequate for assessing these metrics using the Ovine array to genotype Rocky Mountain bighorn sheep herds. However, we conclude that a universal sample size rule may not be able to sufficiently address the complexities that impact genomic kinship and inbreeding estimates. Thus, we recommend that a pilot study and sample size simulation using R code we developed that includes empirical genotypes from a subset of populations of interest would be an effective approach to ensure rigour in estimating genomic kinship and population differentiation.

Trade-offs and efficiencies in optimal budget-constrained multispecies corridor networks.

Conservation biologists recognize that a system of isolated protected areas will be necessary but insufficient to meet biodiversity objectives. Current approaches to connecting core conservation areas through corridors consider optimal corridor placement based on a single optimization goal: commonly, maximizing the movement for a target species across a network of protected areas. We show that designing corridors for single species based on purely ecological criteria leads to extremely expensive linkages that are suboptimal for multispecies connectivity objectives. Similarly, acquiring the least-expensive linkages leads to ecologically poor solutions. We developed algorithms for optimizing corridors for multispecies use given a specific budget. We applied our approach in western Montana to demonstrate how the solutions may be used to evaluate trade-offs in connectivity for 2 species with different habitat requirements, different core areas, and different conservation values under different budgets. We evaluated corridors that were optimal for each species individually and for both species jointly. Incorporating a budget constraint and jointly optimizing for both species resulted in corridors that were close to the individual species movement-potential optima but with substantial cost savings. Our approach produced corridors that were within 14% and 11% of the best possible corridor connectivity for grizzly bears (Ursus arctos) and wolverines (Gulo gulo), respectively, and saved 75% of the cost. Similarly, joint optimization under a combined budget resulted in improved connectivity for both species relative to splitting the budget in 2 to optimize for each species individually. Our results demonstrate economies of scale and complementarities conservation planners can achieve by optimizing corridor designs for financial costs and for multiple species connectivity jointly. We believe that our approach will facilitate corridor conservation by reducing acquisition costs and by allowing derived corridors to more closely reflect conservation priorities.

Nature vs. Nurture: Evidence for Social Learning of Conflict Behaviour in Grizzly Bears.

The propensity for a grizzly bear to develop conflict behaviours might be a result of social learning between mothers and cubs, genetic inheritance, or both learning and inheritance. Using non-invasive genetic sampling, we collected grizzly bear hair samples during 2011-2014 across southwestern Alberta, Canada. We targeted private agricultural lands for hair samples at grizzly bear incident sites, defining an incident as an occurrence in which the grizzly bear caused property damage, obtained anthropogenic food, or killed or attempted to kill livestock or pets. We genotyped 213 unique grizzly bears (118 M, 95 F) at 24 microsatellite loci, plus the amelogenin marker for sex. We used the program COLONY to assign parentage. We evaluated 76 mother-offspring relationships and 119 father-offspring relationships. We compared the frequency of problem and non-problem offspring from problem and non-problem parents, excluding dependent offspring from our analysis. Our results support the social learning hypothesis, but not the genetic inheritance hypothesis. Offspring of problem mothers are more likely to be involved in conflict behaviours, while offspring from non-problem mothers are not likely to be involved in incidents or human-bear conflicts themselves (Barnard's test, p = 0.05, 62.5% of offspring from problem mothers were problem bears). There was no evidence that offspring are more likely to be involved in conflict behaviour if their fathers had been problem bears (Barnard's test, p = 0.92, 29.6% of offspring from problem fathers were problem bears). For the mother-offspring relationships evaluated, 30.3% of offspring were identified as problem bears independent of their mother's conflict status. Similarly, 28.6% of offspring were identified as problem bears independent of their father's conflict status. Proactive mitigation to prevent female bears from becoming problem individuals likely will help prevent the perpetuation of conflicts through social learning.

Demographic mechanisms underpinning genetic assimilation of remnant groups of a large carnivore.

Current range expansions of large terrestrial carnivores are occurring following human-induced range contraction. Contractions are often incomplete, leaving small remnant groups in refugia throughout the former range. Little is known about the underlying ecological and evolutionary processes that influence how remnant groups are affected during range expansion. We used data from a spatially explicit, long-term genetic sampling effort of grizzly bears (Ursus arctos) in the Northern Continental Divide Ecosystem (NCDE), USA, to identify the demographic processes underlying spatial and temporal patterns of genetic diversity. We conducted parentage analysis to evaluate how reproductive success and dispersal contribute to spatio-temporal patterns of genetic diversity in remnant groups of grizzly bears existing in the southwestern (SW), southeastern (SE) and east-central (EC) regions of the NCDE. A few reproductively dominant individuals and local inbreeding caused low genetic diversity in peripheral regions that may have persisted for multiple generations before eroding rapidly (approx. one generation) during population expansion. Our results highlight that individual-level genetic and reproductive dynamics play critical roles during genetic assimilation, and show that spatial patterns of genetic diversity on the leading edge of an expansion may result from historical demographic patterns that are highly ephemeral.

Current approaches using genetic distances produce poor estimates of landscape resistance to interindividual dispersal.

Landscape resistance reflects how difficult it is for genes to move across an area with particular attributes (e.g. land cover, slope). An increasingly popular approach to estimate resistance uses Mantel and partial Mantel tests or causal modelling to relate observed genetic distances to effective distances under alternative sets of resistance parameters. Relatively few alternative sets of resistance parameters are tested, leading to relatively poor coverage of the parameter space. Although this approach does not explicitly model key stochastic processes of gene flow, including mating, dispersal, drift and inheritance, bias and precision of the resulting resistance parameters have not been assessed. We formally describe the most commonly used model as a set of equations and provide a formal approach for estimating resistance parameters. Our optimization finds the maximum Mantel r when an optimum exists and identifies the same resistance values as current approaches when the alternatives evaluated are near the optimum. Unfortunately, even where an optimum existed, estimates from the most commonly used model were imprecise and were typically much smaller than the simulated true resistance to dispersal. Causal modelling using Mantel significance tests also typically failed to support the true resistance to dispersal values. For a large range of scenarios, current approaches using a simple correlational model between genetic and effective distances do not yield accurate estimates of resistance to dispersal. We suggest that analysts consider the processes important to gene flow for their study species, model those processes explicitly and evaluate the quality of estimates resulting from their model.

Spatial capture--recapture models for jointly estimating population density and landscape connectivity.

Population size and landscape connectivity are key determinants of population viability, yet no methods exist for simultaneously estimating density and connectivity parameters. Recently developed spatial capture--recapture (SCR) models provide a framework for estimating density of animal populations but thus far have not been used to study connectivity. Rather, all applications of SCR models have used encounter probability models based on the Euclidean distance between traps and animal activity centers, which implies that home ranges are stationary, symmetric, and unaffected by landscape structure. In this paper we devise encounter probability models based on "ecological distance," i.e., the least-cost path between traps and activity centers, which is a function of both Euclidean distance and animal movement behavior in resistant landscapes. We integrate least-cost path models into a likelihood-based estimation scheme for spatial capture-recapture models in order to estimate population density and parameters of the least-cost encounter probability model. Therefore, it is possible to make explicit inferences about animal density, distribution, and landscape connectivity as it relates to animal movement from standard capture-recapture data. Furthermore, a simulation study demonstrated that ignoring landscape connectivity can result in negatively biased density estimators under the naive SCR model.

Balancing precision and risk: should multiple detection methods be analyzed separately in N-mixture models?

Using multiple detection methods can increase the number, kind, and distribution of individuals sampled, which may increase accuracy and precision and reduce cost of population abundance estimates. However, when variables influencing abundance are of interest, if individuals detected via different methods are influenced by the landscape differently, separate analysis of multiple detection methods may be more appropriate. We evaluated the effects of combining two detection methods on the identification of variables important to local abundance using detections of grizzly bears with hair traps (systematic) and bear rubs (opportunistic). We used hierarchical abundance models (N-mixture models) with separate model components for each detection method. If both methods sample the same population, the use of either data set alone should (1) lead to the selection of the same variables as important and (2) provide similar estimates of relative local abundance. We hypothesized that the inclusion of 2 detection methods versus either method alone should (3) yield more support for variables identified in single method analyses (i.e. fewer variables and models with greater weight), and (4) improve precision of covariate estimates for variables selected in both separate and combined analyses because sample size is larger. As expected, joint analysis of both methods increased precision as well as certainty in variable and model selection. However, the single-method analyses identified different variables and the resulting predicted abundances had different spatial distributions. We recommend comparing single-method and jointly modeled results to identify the presence of individual heterogeneity between detection methods in N-mixture models, along with consideration of detection probabilities, correlations among variables, and tolerance to risk of failing to identify variables important to a subset of the population. The benefits of increased precision should be weighed against those risks. The analysis framework presented here will be useful for other species exhibiting heterogeneity by detection method.