Examining the occupancy–density relationship for a low‐density carnivore

Summary The challenges associated with monitoring low-density carnivores across large landscapes have limited the ability to implement and evaluate conservation and management strategies for such species. Non-invasive sampling techniques and advanced statistical approaches have alleviated some of these challenges and can even allow for spatially explicit estimates of density, one of the most valuable wildlife monitoring tools. For some species, individual identification comes at no cost when unique attributes (e.g. pelage patterns) can be discerned with remote cameras, while other species require viable genetic material and expensive laboratory processing for individual assignment. Prohibitive costs may still force monitoring efforts to use species distribution or occupancy as a surrogate for density, which may not be appropriate under many conditions. Here, we used a large-scale monitoring study of fisher Pekania pennanti to evaluate the effectiveness of occupancy as an approximation to density, particularly for informing harvest management decisions. We combined remote cameras with baited hair snares during 2013–2015 to sample across a 70 096-km2 region of western New York, USA. We fit occupancy and Royle–Nichols models to species detection–non-detection data collected by cameras, and spatial capture–recapture (SCR) models to individual encounter data obtained by genotyped hair samples. Variation in the state variables within 15-km2 grid cells was modelled as a function of landscape attributes known to influence fisher distribution. We found a close relationship between grid cell estimates of fisher state variables from the models using detection–non-detection data and those from the SCR model, likely due to informative spatial covariates across a large landscape extent and a grid cell resolution that worked well with the movement ecology of the species. Fisher occupancy and density were both positively associated with the proportion of coniferous-mixed forest and negatively associated with road density. As a result, spatially explicit management recommendations for fisher were similar across models, though relative variation was dampened for the detection–non-detection data. Synthesis and applications. Our work provides empirical evidence that models using detection–non-detection data can make similar inferences regarding relative spatial variation of the focal population to models using more expensive individual encounters when the selected spatial grain approximates or is marginally smaller than home range size. When occupancy alone is chosen as a cost-effective state variable for monitoring, simulation and sensitivity analyses should be used to understand how inferences from detection–non-detection data will be affected by aspects of study design and species ecology.

[1]  J. Andrew Royle,et al.  Spatially explicit models for inference about density in unmarked or partially marked populations , 2011, 1112.3250.

[2]  M. Rossetto,et al.  Microsatellite markers for American mink ( Mustela vison ) and ermine ( Mustela erminea ) , 1999 .

[3]  Kevin J. Gaston,et al.  Abundance–occupancy relationships , 2000 .

[4]  C. Davis,et al.  Isolation, variability, and cross‐species amplification of polymorphic microsatellite loci in the family Mustelidae , 1998, Molecular ecology.

[5]  J. Hodges,et al.  Adding Spatially-Correlated Errors Can Mess Up the Fixed Effect You Love , 2010 .

[6]  Darryl I. MacKenzie,et al.  Designing occupancy studies: general advice and allocating survey effort , 2005 .

[7]  S. Piertney,et al.  Microsatellite primers for the Eurasian otter. , 1998, Molecular ecology.

[8]  William J. Zielinski,et al.  Carnivore Translocations and Conservation: Insights from Population Models and Field Data for Fishers (Martes pennanti) , 2012, PloS one.

[9]  D L Borchers,et al.  Spatially Explicit Maximum Likelihood Methods for Capture–Recapture Studies , 2008, Biometrics.

[10]  S. Creel,et al.  Population size estimation in Yellowstone wolves with error‐prone noninvasive microsatellite genotypes , 2003, Molecular ecology.

[11]  Christopher J. Kyle,et al.  Empirical comparison of density estimators for large carnivores , 2010 .

[12]  W. Krohn,et al.  Home Range Characteristics of Adult Fishers , 1989 .

[13]  Darryl I. MacKenzie,et al.  Assessing the fit of site-occupancy models , 2004 .

[14]  J. Nichols,et al.  Monitoring carnivore populations at the landscape scale: occupancy modelling of tigers from sign surveys , 2011 .

[15]  D. Macdonald,et al.  To bait or not to bait: A comparison of camera-trapping methods for estimating leopard Panthera pardus density , 2014 .

[16]  James H. Brown On the Relationship between Abundance and Distribution of Species , 1984, The American Naturalist.

[17]  Jacob S. Ivan,et al.  rSPACE: Spatially based power analysis for conservation and ecology , 2015 .

[18]  J. Andrew Royle,et al.  ESTIMATING SITE OCCUPANCY RATES WHEN DETECTION PROBABILITIES ARE LESS THAN ONE , 2002, Ecology.

[19]  D. Dawson,et al.  Occupancy in continuous habitat , 2012 .

[20]  James D. Nichols,et al.  Monitoring of biological diversity in space and time , 2001 .

[21]  Larissa L. Bailey,et al.  Inference for Occupancy and Occupancy Dynamics , 2011 .

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

[23]  P. Lancaster,et al.  Fishers, Farms, and Forests in Eastern North America , 2008, Environmental management.

[24]  Erin M. Bayne,et al.  REVIEW: Wildlife camera trapping: a review and recommendations for linking surveys to ecological processes , 2015 .

[25]  Jacob S. Ivan,et al.  Spatially Explicit Power Analyses for Occupancy‐Based Monitoring of Wolverine in the U.S. Rocky Mountains , 2014, Conservation biology : the journal of the Society for Conservation Biology.

[26]  Julia P. G. Jones Monitoring species abundance and distribution at the landscape scale , 2011 .

[27]  Mevin B. Hooten,et al.  Spatial occupancy models for large data sets , 2013 .

[28]  Michael K. Schwartz,et al.  Development of 22 new microsatellite loci for fishers (Martes pennanti) with variability results from across their range , 2007 .

[29]  Kenneth H. Pollock,et al.  Large scale wildlife monitoring studies: statistical methods for design and analysis , 2002 .

[30]  N. Pettorelli,et al.  Management by proxy?:the use of indices in applied ecology , 2015 .

[31]  P. Diggle,et al.  Negative Binomial Quadrat Counts and Point Processes , 2016 .

[32]  Richard B. Chandler,et al.  unmarked: An R Package for Fitting Hierarchical Models of Wildlife Occurrence and Abundance , 2011 .

[33]  Rahel Sollmann,et al.  Improving density estimates for elusive carnivores: Accounting for sex-specific detection and movements using spatial capture–recapture models for jaguars in central Brazil , 2011 .

[34]  Suresh A. Sethi,et al.  Guidelines for MSAT and SNP panels that lead to high-quality data for genetic mark–recapture studies , 2014 .

[35]  J. Andrew Royle,et al.  Management decision making for fisher populations informed by occupancy modeling , 2016 .

[36]  T. J. Roper,et al.  Reliable microsatellite genotyping of the Eurasian badger (Meles meles) using faecal DNA , 2003, Molecular ecology.

[37]  J. Wiens Spatial Scaling in Ecology , 1989 .

[38]  John D. J. Clare,et al.  Predicting bobcat abundance at a landscape scale and evaluating occupancy as a density index in central Wisconsin , 2015 .

[39]  Thomas L. Serfass,et al.  Ten new polymorphic microsatellite loci for North American river otters (Lontra canadensis) and their utility in related mustelids , 2005 .

[40]  Richard B. Chandler,et al.  Spatially explicit integrated population models , 2014 .

[41]  J. Andrew Royle,et al.  Likelihood analysis of spatial capture-recapture models for stratified or class structured populations , 2015 .

[42]  J Andrew Royle,et al.  A hierarchical model for spatial capture-recapture data. , 2008, Ecology.