Entangling snakebite dynamics: the spatiotemporal role of rainfall on snake envenoming in Colombia
The role of climate forcing on the population dynamics of infectious diseases has typically been addressed via retrospective analyses of aggregated incidence records over whole political regions. A central question in epidemiology has been whether seasonal and interannual cycles are driven by climate variation, or instead generated by other factors such as poverty or underreporting. Here, we use process-based models to determine the role of rainfall in the dynamics of snakebite, which is one of the main neglected tropical diseases around the world. We combined space-time datasets of snakebite incidence and rainfall for Colombia in combination with stochastic epidemiological models and iterated filtering methods to show the response to rainfall forcing, specifically, modulating the encounter frequency with venomous snakes. We identified six zones with different rainfall patterns to demonstrate that the relationship between rainfall and snakebite incidence was heterogeneous. Rainfall only drives snakebite incidence in regions with marked dry seasons, where rainfall becomes the limiting resource. In addition, the encounter frequency also differs between regions, and it is higher in regions where Bothrops atrox can be found. Our results show how the heterogeneous spatial distribution of snakebite risk seasonality in the country may be related to important traits of venomous snakes' natural history.