Projecting Invasion Risk of Non-Native Watersnakes (Nerodia fasciata and Nerodia sipedon) in the Western United States

Species distribution models (SDMs) are increasingly used to project the potential distribution of introduced species outside their native range. Such studies rarely explicitly evaluate potential conflicts with native species should the range of introduced species expand. Two snake species native to eastern North America, Nerodia fasciata and Nerodia sipedon, have been introduced to California where they represent a new stressor to declining native amphibians, fish, and reptiles. To project the potential distributions of these non-native watersnakes in western North America, we built ensemble SDMs using MaxEnt, Boosted Regression Trees, and Random Forests and habitat and climatic variables. We then compared the overlap between the projected distribution of invasive watersnakes and the distributions of imperiled native amphibians, fish, and reptiles that can serve as prey or competitors for the invaders, to estimate the risk to native species posed by non-native watersnakes. Large areas of western North America were projected to be climatically suitable for both species of Nerodia according to our ensemble SDMs, including much of central California. The potential distributions of both N. fasciata and N. sipedon overlap extensively with the federally threatened Giant Gartersnake, Thamnophis gigas, which inhabits a similar ecological niche. N. fasciata also poses risk to the federally threatened California Tiger Salamander, Ambystoma californiense, whereas N. sipedon poses risk to some amphibians of conservation concern, including the Foothill Yellow-legged Frog, Rana boylii. We conclude that non-native watersnakes in California can likely inhabit ranges of several native species of conservation concern that are expected to suffer as prey or competing species for these invaders. Action should be taken now to eradicate or control these invasions before detrimental impacts on native species are widespread. Our methods can be applied broadly to quantify the risk posed by incipient invasions to native biodiversity.