Linking landscape characteristics to local grizzly bear abundance using multiple detection methods in a hierarchical model

Few studies link habitat to grizzly bear Ursus arctos abundance and these have not accounted for the variation in detection or spatial autocorrelation. We collected and genotyped bear hair in and around Glacier National Park in northwestern Montana during the summer of 2000. We developed a hierarchical Markov chain Monte Carlo model that extends the existing occupancy and count models by accounting for (1) spatially explicit variables that we hypothesized might influence abundance; (2) separate sub-models of detection probability for two distinct sampling methods (hair traps and rub trees) targeting different segments of the population; (3) covariates to explain variation in each sub-model of detection; (4) a conditional autoregressive term to account for spatial autocorrelation; (5) weights to identify most important variables. Road density and per cent mesic habitat best explained variation in female grizzly bear abundance; spatial autocorrelation was not supported. More female bears were predicted in places with lower road density and with more mesic habitat. Detection rates of females increased with rub tree sampling effort. Road density best explained variation in male grizzly bear abundance and spatial autocorrelation was supported. More male bears were predicted in areas of low road density. Detection rates of males increased with rub tree and hair trap sampling effort and decreased over the sampling period. We provide a new method to (1) incorporate multiple detection methods into hierarchical models of abundance; (2) determine whether spatial autocorrelation should be included in final models. Our results suggest that the influence of landscape variables is consistent between habitat selection and abundance in this system.

[1]  J. Andrew Royle.  Generalized estimators of avian abundance from count survey data , 2004 .

[2]  J. Nichols,et al.  Tigers and their prey: Predicting carnivore densities from prey abundance. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Kenneth H. Pollock,et al.  Separating Components of the Detection Process With Combined Methods: An Example With Northern Bobwhite , 2010 .

[4]  Katherine C. Kendall,et al.  Demography and Genetic Structure of a Recovering Grizzly Bear Population , 2009 .

[5]  A. Holmes,et al.  State-and-Transition Models for Assessing Grasshopper Sparrow Habitat Use , 2010 .

[6]  J. Andrew Royle,et al.  ESTIMATING ABUNDANCE FROM REPEATED PRESENCE–ABSENCE DATA OR POINT COUNTS , 2003 .

[7]  James D. Nichols.  The wildlife picture index, monitoring and conservation. , 2010 .

[8]  G. Stenhouse,et al.  Grizzly bears and forestry I: Selection of clearcuts by grizzly bears in west-central Alberta, Canada , 2004 .

[9]  J. Berger,et al.  Optimal predictive model selection , 2004, math/0406464.

[10]  J. Andrew Royle,et al.  HIERARCHICAL SPATIAL MODELS OF ABUNDANCE AND OCCURRENCE FROM IMPERFECT SURVEY DATA , 2007 .

[11]  Denis Hémon,et al.  On the variance of the sample correlation between two independent lattice processes , 1981, Journal of Applied Probability.

[12]  B. Mclellan,et al.  Natal dispersal of grizzly bears , 2001 .

[13]  E. Crist,et al.  Application of the Tasseled Cap concept to simulated thematic mapper data , 1984 .

[14]  W. Thuiller,et al.  Predicting species distribution: offering more than simple habitat models. , 2005, Ecology letters.

[15]  Katherine C. Kendall,et al.  A simulation test of the effectiveness of several methods for error‐checking non‐invasive genetic data , 2005 .

[16]  Mark S Boyce,et al.  Grizzly bear habitat selection is scale dependent. , 2007, Ecological applications : a publication of the Ecological Society of America.

[17]  T. L. Scheid-Cook,et al.  A Sand County Almanac , 1949 .

[18]  Marc Kéry,et al.  Hierarchical modelling and estimation of abundance and population trends in metapopulation designs. , 2010, The Journal of animal ecology.

[19]  David A. Roon,et al.  Multiple data sources improve DNA-based mark-recapture population estimates of grizzly bears. , 2008, Ecological applications : a publication of the Ecological Society of America.

[20]  James E. Hines,et al.  Approaches for the direct estimation of u , and demographic contributions to u , using capture-recapture data , 2002 .

[21]  Montana.,et al.  Grizzly bear recovery plan , 1982 .

[22]  J. Andrew Royle.  N‐Mixture Models for Estimating Population Size from Spatially Replicated Counts , 2004, Biometrics.

[23]  H. Pulliam,et al.  Sources, Sinks, and Population Regulation , 1988, The American Naturalist.

[24]  J. Andrew Royle,et al.  MODELING AVIAN ABUNDANCE FROM REPLICATED COUNTS USING BINOMIAL MIXTURE MODELS , 2005 .

[25]  William A Link,et al.  Model weights and the foundations of multimodel inference. , 2006, Ecology.

[26]  R. O’Hara,et al.  A review of Bayesian variable selection methods: what, how and which , 2009 .

[27]  J. S. Waller,et al.  Landscape Evaluation of Grizzly Bear Habitat in Western Montana , 1999 .

[28]  E. Beckwith,et al.  DROPOUT : a program to identify problem loci and samples for noninvasive genetic samples in a capture-mark-recapture framework , 2005 .

[29]  J. Battin.  When Good Animals Love Bad Habitats: Ecological Traps and the Conservation of Animal Populations , 2004 .

[30]  Kenneth H. Pollock,et al.  Bayesian spatial modeling of data from avian point count surveys , 2008 .

[31]  R. Dorazio.  On the choice of statistical models for estimating occurrence and extinction from animal surveys. , 2007, Ecology.

[32]  David A. Roon,et al.  Grizzly Bear Density in Glacier National Park, Montana , 2008 .

[33]  D. MacKenzie.  Modeling the Probability of Resource Use: The Effect of, and Dealing with, Detecting a Species Imperfectly , 2006 .

[34]  D. P. A E T K A U,et al.  An empirical exploration of data quality in DNA-based population inventories , 2003 .

[35]  M. Bélisle,et al.  Effects of Landscape Structure on Male Density and Spacing Patterns in Wild Turkeys (Meleagris gallopavo) Depend on Winter Severity , 2009 .

[36]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[37]  Kevin S. McKelvey,et al.  Felid sex identification based on noninvasive genetic samples , 2005 .

[38]  I. Stirling,et al.  Microsatellite analysis of population structure in Canadian polar bears , 1995, Molecular ecology.

[39]  J. S. Waller,et al.  Grizzly Bear Habitat Selection in the Swan Mountains, Montana , 1997 .

[40]  Marc J. Mazerolle,et al.  LANDSCAPE CHARACTERISTICS INFLUENCE POND OCCUPANCY BY FROGS AFTER ACCOUNTING FOR DETECTABILITY , 2005 .

[41]  G. McDermid,et al.  Dynamic wildlife habitat models: Seasonal foods and mortality risk predict occupancy-abundance and habitat selection in grizzly bears , 2010 .

[42]  B. V. Horne,et al.  DENSITY AS A MISLEADING INDICATOR OF HABITAT QUALITY , 1983 .

[43]  J. S. Waller,et al.  EFFECTS OF TRANSPORTATION INFRASTRUCTURE ON GRIZZLY BEARS IN NORTHWESTERN MONTANA , 2005 .

[44]  M. Conroy,et al.  Analysis and Management of Animal Populations , 2002 .

[45]  Gordon B. Stenhouse,et al.  Grizzly bear movements relative to roads: application of step selection functions , 2010 .

[46]  J. Andrew Royle,et al.  Hierarchical Spatiotemporal Matrix Models for Characterizing Invasions , 2007, Biometrics.

[47]  J. Gibbs,et al.  EFFECTIVENESS OF CALL-BROADCAST SURVEYS FOR MONITORING MARSH BIRDS , 2005 .

[48]  Pierre Vernier,et al.  Grizzly bear use of open, closed, and restricted forestry roads , 2002 .

[49]  Scott E. Nielsen,et al.  Modelling the spatial distribution of human-caused grizzly bear mortalities in the Central Rockies ecosystem of Canada , 2004 .

[50]  Christopher Servheen.  Grizzly bear food habits, movements, and habitat selection in the Mission Mountains, Montana. , 1983 .

[51]  P. Taberlet,et al.  Sexing free‐ranging brown bears Ursus arctos using hairs found in the field , 1993, Molecular ecology.

[52]  Carl J. Schwarz,et al.  Discussion comments on: 'Approaches for the direct estimation of u and demographic contributions to u using capture-recapture data' , 2002 .

[53]  Raymond A. Webster,et al.  Bayesian Spatial Modeling of Data from Unit‐Count Surveys of Fish in Streams , 2008 .

[54]  P. Dixon,et al.  Accounting for Spatial Pattern When Modeling Organism- Environment Interactions , 2022 .

[55]  S. Ennis,et al.  A PCR-based sex-determination assay in cattle based on the bovine amelogenin locus. , 2009, Animal genetics.

[56]  Hans R. Zuuring,et al.  Relationships among grizzly bears, roads and habitat in the Swan Mountains, Montana , 1996 .

[57]  William A. Link,et al.  Bayesian Inference: With Ecological Applications , 2009 .

[58]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[59]  L. Maiorano,et al.  Assessing habitat quality for conservation using an integrated occurrence-mortality model , 2009 .

[60]  I. MacKenzie.  Occupancy estimation and modeling , 2013 .

[61]  R. Swihart,et al.  Absent or undetected? Effects of non-detection of species occurrence on wildlife-habitat models , 2004 .

[62]  M. Kittlein,et al.  Diet, habitat use, and relative abundance of pampas fox (Pseudalopex gymnocercus) in northern Patagonia, Argentina , 2005 .

[63]  G. Stenhouse,et al.  Grizzly Ursus arctos and black bear U. americanus densities in the interior mountains of North America , 2005 .

[64]  Andrew Thomas,et al.  WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility , 2000, Stat. Comput..

[65]  Brian P. Dreher,et al.  Noninvasive Estimation of Black Bear Abundance Incorporating Genotyping Errors and Harvested Bear , 2007 .

[66]  S. DeStefano,et al.  Scrub-Shrub Bird Habitat Associations at Multiple Spatial Scales in Beaver Meadows in Massachusetts , 2009 .

[67]  Gordon B. Stenhouse,et al.  A habitat-based framework for grizzly bear conservation in Alberta , 2006 .

[68]  J. Andrew Royle,et al.  Hierarchical Modeling and Inference in Ecology: The Analysis of Data from Populations, Metapopulations and Communities , 2008 .

[69]  Scott E. Nielsen,et al.  Can models of presence‐absence be used to scale abundance? Two case studies considering extremes in life history , 2005 .

[70]  W. Kendall,et al.  Making Great Leaps Forward: Accounting for Detectability in Herpetological Field Studies , 2007 .

[71]  M. McCarthy.  Bayesian Methods for Ecology: CASE STUDIES , 2007 .

[72]  M. Proctor,et al.  ESTIMATING GRIZZLY BEAR DISTRIBUTION AND ABUNDANCE RELATIVE TO HABITAT AND HUMAN INFLUENCE , 2004 .