Integrating human behavior and snake ecology with agent-based models to predict snakebite in high risk landscapes

Snakebite causes more than 1.8 million envenoming cases annually and is a major cause of death in the tropics especially for poor farmers. While both social and ecological factors influence the chance encounter between snakes and people, the spatio-temporal processes underlying snakebites remain poorly explored. Previous research has focused on statistical correlates between snakebites and ecological, sociological, or environmental factors, but the human and snake behavioral patterns that drive the spatio-temporal process have not yet been integrated into a single model. Here we use a bottom-up simulation approach using agent-based modelling (ABM) parameterized with datasets from Sri Lanka, a snakebite hotspot, to characterise the mechanisms of snakebite and identify risk factors. Spatio-temporal dynamics of snakebite risks are examined through the model incorporating six snake species and three farmer types (rice, tea, and rubber). We find that snakebites are mainly climatically driven, but the risks also depend on farmer types due to working schedules as well as species present in landscapes. Snake species are differentiated by both distribution and by habitat preference, and farmers are differentiated by working patterns that are climatically driven, and the combination of these factors leads to unique encounter rates for different landcover types as well as locations. Validation using epidemiological studies demonstrated that our model can explain observed patterns, including temporal patterns of snakebite incidence, and relative contribution of bites by each snake species. Our predictions can be used to generate hypotheses and inform future studies and decision makers. Additionally, our model is transferable to other locations with high snakebite burden as well.