Population and spatial dynamics of desert bighorn sheep in Grand Canyon during an outbreak of respiratory pneumonia

IntroductionTerrestrial species in riverine ecosystems face unique constraints leading to diverging patterns of population structure, connectivity, and disease dynamics. Desert bighorn sheep (Ovis canadensis nelsoni) in Grand Canyon National Park, a large native population in the southwestern USA, offer a unique opportunity to evaluate population patterns and processes in a remote riverine system with ongoing anthropogenic impacts. We integrated non-invasive, invasive, and citizen-science methods to address questions on abundance, distribution, disease status, genetic structure, and habitat fragmentation.MethodsWe compiled bighorn sightings collected during river trips by park staff, commercial guides, and private citizens from 2000–2018 and captured bighorn in 2010–2016 to deploy GPS collars and test for disease. From 2011–2015, we non-invasively collected fecal samples and genotyped them at 9–16 microsatellite loci for individual identification and genetic structure. We used assignment tests to evaluate genetic structure and identify subpopulations, then estimated gene flow and recent migration to evaluate fragmentation. We used spatial capture-recapture to estimate annual population size, distribution, and trends after accounting for spatial variation in detection with a resource selection function model.Results and discussionFrom 2010–2018, 3,176 sightings of bighorn were reported, with sightings of 56–145 bighorn annually on formal surveys. From 2012–2016, bighorn exhibiting signs of respiratory disease were observed along the river throughout the park. Of 25 captured individuals, 56% were infected by Mycoplasma ovipneumoniae, a key respiratory pathogen, and 81% were recently exposed. Pellet sampling for population estimation from 2011–2015 yielded 1,250 genotypes and 453 individuals. We detected 6 genetic clusters that exhibited mild to moderate genetic structure (FST 0.022–0.126). The river, distance, and likely topography restricted recent gene flow, but we detected cross-river movements in one section via genetic recaptures, no subpopulation appeared completely isolated, and genetic diversity was among the highest reported. Recolonization of one large stretch of currently empty habitat appears limited by the constrained topology of this system. Annual population estimates ranged 536–552 (95% CrI range 451–647), lamb:ewe ratios varied, and no significant population decline was detected. We provide a multi-method sampling framework useful for sampling other wildlife in remote riverine systems.