Quantifying spatiotemporal occupancy dynamics and multi-year core-use areas at a species range boundary

• Published in: Diversity & distributions 2020-07, Vol.26 (7), p.795-805
• Main Authors: Hostetter, Nathan J., Ryan, Daniel, Grosshuesch, David, Catton, Timothy, Malick-Wahls, Sarah, Smith, Tamara A., Gardner, Beth
• Format: Article
• Language: English
• Subjects: Anthropogenic factors; BIODIVERSITY RESEARCH; Canada lynx; Climate change; Climate prediction; Colonization; Coniferous forests; core area; detection probability; dynamic occupancy; Endangered & extinct species; Endangered species; Forests; Habitats; Landscape; Lynx canadensis; Lynx lynx; Monitoring; Occupancy; range dynamics; site occupancy; Snowfall; turnover; Wildlife conservation
• ISSN: 1366-9516; 1472-4642
• DOI: 10.1111/ddi.13066

Aim Many species face large‐scale range contractions and predicted distributional shifts in response to climate change, shifting forest characteristics and anthropogenic disturbances. Canada lynx (Lynx canadensis) are listed as threatened under the U.S. Endangered Species Act and were recently recommended for delisting. Predicted climate‐driven losses in habitat quality and quantity may negatively affect the north‐eastern Minnesota lynx population, one of the six remaining resident populations in the contiguous United States. We develop a large‐scale monitoring protocol and dynamic occupancy modelling framework to identify multi‐year core‐use areas and track spatiotemporal occurrence at the southern periphery of the species range. Location North‐eastern Minnesota lynx geographic unit, Superior National Forest and designated critical habitat, Minnesota, USA. Methods Spatially and temporally replicated snow track surveys were used to collect lynx detection/non‐detection data across five winters (2014–15 to 2018–19) covering >17,000 km within the 22,100 km2 study area. We used a dynamic occupancy model to evaluate lynx occupancy, persistence, colonization and habitat covariates affecting these processes. Results Lynx occupancy probabilities displayed high spatial and temporal variability, with grid cell‐specific probabilities ranging from 0.0 in periphery regions to consistently near 1.0 in multi‐year core‐use areas, indicating low turnover rates in those areas. Lynx colonization and persistence increased in areas with more evergreen forest and greater average snowfall, while forest characteristics (3–5 and 10–30 m vegetation density) had mixed relationships with occupancy dynamics. We identified 55 grid cells classified as multi‐year core‐use areas across relatively contiguous regions of high average snowfall and per cent conifer forest. Main conclusions Our study demonstrates a landscape‐scale multi‐year monitoring programme assessing the effects of habitat characteristics and anthropogenic factors on species distributional changes and landscape‐level occupancy dynamics. Our framework incorporating landscape‐scale resource selection, core‐use area concepts and dynamic occupancy models provides a flexible approach to identify population‐level mechanisms driving species persistence and key areas for conservation protection.