Species distribution modelling for Rhipicephalus microplus (Acari : Ixodidae) in Benin, West Africa : comparing datasets and modelling algorithms

Abstract

Rhipicephalus microplus is one of the most widely distributed and economically important ticks, transmitting Babesia bigemina, B. bovis and Anaplasma marginale. It was recently introduced to West Africa on live animals originating from Brazil. Knowing the precise environmental suitability for the tick would allow veterinary health officials to draft vector control strategies for different regions of the country. To test the performance of modelling algorithms and different sets of environmental explanatory variables, species distribution models for this tick species in Benin were developed using generalized linear models, lin-ear discriminant analysis and random forests. The training data for these models were a dataset containing reported absence or presence in 104 farms, randomly selected across Benin. These farms were sampled at the end of the rainy season, which corresponds with an annual peak in tick abundance. Two environmental datasets for the country of Benin were compared: one based on interpolated climate data (WorldClim) and one based on remotely sensed images (MODIS). The pixel size for both environmental datasets was 1 km. Highly suitable areas occurred mainly along the warmer and humid coast extending northwards to central Benin. The northern hot and drier areas were found to be unsuitable. The models developed and tested on data from the entire country were generally found to perform well, having an AUC value greater than 0.92. Although statistically significant, only small differences in accuracy measures were found between the modelling algorithms, or between the environmental datasets. The resulting risk maps differed nonetheless. Models based on interpolated climate suggested gradual variations in habitat suitability, while those based on remotely sensed data indicated a sharper contrast between suitable and unsuitable areas, and a patchy distribution of the suitable areas. Remotely sensed data yielded more spatial detail in the predictions. When computing accuracy measures on a subset of data along the invasion front, the modelling technique Random Forest outperformed the other modelling approaches, and results with MODIS-derived variables were better than those using WorldClim data.