Modeling relative habitat suitability of southern Florida for invasive Burmese pythons (Python molurus bivittatus)

Invasive Burmese pythons are altering the ecology of southern Florida and their distribution is expanding northward. Understanding their habitat use is an important step in understanding the pathways of the invasion. This study identifies key landscape variables in predicting relative habitat suitability for pythons at the present stage of invasion through presence-only ecological niche modeling using geographical sampling bias correction. We used 2014 presence-only observations from the EDDMapS database and three landscape variables to model habitat suitability: fine-scale land cover, home range-level land cover, and distance to open freshwater or wetland. Ten geographical sampling bias correction scenarios based on road presence and sampling effort were evaluated to improve the efficacy of modeling. The best performing models treated road presence as a binary factor rather than a continuous decrease in sampling effort with distance from roads. Home range-level cover contributed the most to the final prediction, followed by proximity to water and fine-scale land cover. Estuarine habitat and freshwater wetlands were the most important variables to contribute to python habitat suitability at both the home range-level and fine-scale. Suitability was highest within 30 m of open freshwater and wetlands. This study provides quantifiable, predictive relationships between habitat types and python presence at the current stage of invasion. This knowledge can elucidate future targeted studies of python habitat use and behavior and help inform management efforts. Furthermore, it illustrates how estimates of relative habitat suitability derived from MaxEnt can be improved by both multi-scale perspectives on habitat and consideration of a variety of bias correction scenarios for selecting background points.