Integrating data to assess occupancy patterns of an endangered bumble bee

There is growing interest in integrating community science data with structured monitoring data to estimate changes in distribution patterns of imperiled species, including pollinators. However, significant challenges remain in determining how unstructured community science data should be incorporated into formal analyses of species distributions. We developed a dynamic framework for combining community science and structured monitoring data of bumble bees to estimate changes in occupancy of rusty-patched bumble bees (Bombus affinis), a federally endangered species in the United States. We applied traditional metapopulation theory and accounted for imperfect detection to estimate site-specific extirpation risk and colonization rates across the known distribution of B. affinis in the Upper Midwest (USA). Despite a 144% increase in presence-only detections from 2017 to 2022, occupancy probabilities and the estimated number of occupied sites remained static or declined slightly across a 4-state region during this period. Our results provide preliminary evidence that the probability of local extirpation risk of B. affinis increased in response to drought, but that effect was tempered with a high number of neighboring patches occupied by B. affinis (i.e., rescue effect). Our framework can be used by managers to track population recovery goals for B. affinis and other bumble bees of conservation concern. In addition, our study highlights the importance of accounting for imperfect detection and addressing spatial sampling biases in bumble bee monitoring efforts, particularly those for which a portion of the monitoring data are generated from community science projects.