Snake co-occurrence patterns are best explained by habitat and hypothesized effects of interspecific interactions

Snakes often occur in species-rich assemblages, and sympatry is thought to be facilitated primarily by low diet overlap, not interspecific interactions. We selected, a priori, three species pairs consisting of species that are morphologically and taxonomically similar and may therefore be likely to engage in interspecific, consumptive competition. We then examined a large-scale database of snake detection/nondetection data and used occupancy modelling to determine whether these species occur together more or less frequently than expected by chance while accounting for variation in detection probability among species and incorporating important habitat categories in the models. For some snakes, we obtained evidence that the probabilities that habitat patches are used are influenced by the presence of potentially competing congeneric species. Specifically, timber rattlesnakes (Crotalus horridus) were less likely than expected by chance to use areas that also contained eastern diamond-backed rattlesnakes (Crotalus adamanteus) when the proportion of evergreen forest was relatively high. Otherwise, they occurred together more often than expected by chance. Complex relationships were revealed between habitat use, detection probabilities and occupancy probabilities of North American racers (Coluber constrictor) and coachwhips (Coluber flagellum) that indicated the probability of competitive exclusion increased with increasing area of grassland habitat, although there was some model uncertainty. Cornsnakes (Pantherophis guttatus or Pantherophis slowinskii) and ratsnakes (Pantherophis alleghaniensis, Pantherophis spiloides, or Pantherophis obsoletus) exhibited differences in habitat selection, but we obtained no evidence that patterns of use for this species pair were influenced by current interspecific interactions. Overall, our results are consistent with the hypothesis that competitive interactions influence snake assemblage composition; the strength of these effects was affected by landscape-scale habitat features. Furthermore, we suggest that current interspecific interactions may influence snake occupancy, challenging the paradigm that contemporary patterns of snake co-occurrence are largely a function of diet partitioning that arose over evolutionary time.