Climate change diminishes the potential habitat of the bont tick (Amblyomma hebraeum): evidence from Mashonaland Central Province, Zimbabwe
Parasites & Vectors volume 15, Article number: 237 (2022)
Understanding the response of vector habitats to climate change is essential for vector management. Increasingly, there is fear that climate change may cause vectors to be more important for animal husbandry in the future. Therefore, knowledge about the current and future spatial distribution of vectors, including ticks (Ixodida), is progressively becoming more critical to animal disease control.
Our study produced present (2018) and future (2050) bont tick (Amblyomma hebraeum) niche models for Mashonaland Central Province, Zimbabwe. Specifically, our approach used the Ensemble algorithm in Biomod2 package in R 3.4.4 with a suite of physical and anthropogenic covariates against the tick’s presence-only location data obtained from cattle dipping facilities.
Our models showed that currently (the year 2018) the bont tick potentially occurs in 17,008 km2, which is 60% of Mashonaland Central Province. However, the models showed that in the future (the year 2050), the bont tick will occur in 13,323 km2, which is 47% of Mashonaland Central Province. Thus, the models predicted an ~ 13% reduction in the potential habitat, about 3685 km2 of the study area. Temperature, elevation and rainfall were the most important variables explaining the present and future potential habitat of the bont tick.
Results of our study are essential in informing programmes that seek to control the bont tick in Mashonaland Central Province, Zimbabwe and similar environments.
There is a notable increase in tick distribution and tick-borne diseases, particularly in the tropical and sub-tropical regions. Tick species of particular interest in the tropics include species in the genera Hyalomma, Amblyomma and Rhipicephalus [1, 2]. Some of these species, e.g., the bont tick (Amblyomma hebraeum), are challenging to control as they have very high reproductive rates and sometimes develop resistance to standard remedies. In addition, the bont tick has a wide spatial distribution because of its ability to survive in both dry and humid environments [3, 4]. Thus, the bont tick is a significant vector for various pathogens, including the obligate intracellular proteobacterium Ehrlichia ruminantium, which causes heartwater disease in ruminants .
Ticks spend much of their life away from their host. As such, ticks are prone to the effects of climatic and other environmental dynamics . The predicted change in climate, primarily to warmer conditions, motivates the modelling of future distribution of ticks using climate change scenarios . Those niche models use various environmental variables, including elevation temperature, rainfall and humidity . Recently, Land Use Land Cover Change (LULCC) has been a significant covariate in modelling tick species and tick host distribution and is often used in modelling overlaps in space and time between human activities and ticks habitats .
The present and future distribution of ticks must be known to support vector management programmes . However, while bont ticks are notorious for animal husbandry, their present and future spatial distribution and the drivers of that distribution are primarily unknown. This study used bont tick presence-only data and a suite of climate change scenarios data and environmental covariates to model the present and future distribution of one of the most notorious ticks in Mashonaland Central province, Zimbabwe. The objectives of this study were twofold: (1) identify the change in the potential habitat of bont tick between current time (2018) and future 2050 and (2) identify the drivers of the spatial distribution of bont tick. Results of this study are essential in informing programmes that seek to control the bont tick in Mashonaland Central Province, Zimbabwe, and similar environments.
The study area was Mashonaland Central Province, Zimbabwe, located between 30.04° E–32.75° E and 17.97° S–15.40° S, comprising 28,347 km2 of land (Fig. 1). According to the Zimbabwe National Census 2012 statistics, the province’s population was 1,152,520 inhabitants. A wide range of land tenure types coexists in the province, including communal areas, newly resettled small scale farming (A1), newly resettled large scale farming (A2), small-scale commercial farming, extensive-scale commercial farming and old resettlement areas . The province comprises regions of varying agricultural potential, including Agro-ecological Region 2, which receives as much as 1000 mm of annual rainfall, to Agro-ecological Region 5, which receives < 450 mm of annual rainfall . The area experiences a wet season (October to April) and a dry season (May to October).
Administratively, Mashonaland Province consists of eight districts: Bindura, Centenary, Mt Darwin, Guruve, Mbire, Rushinga, Shamva and Mazowe. Agriculture is the basis of the economy as the province is predominantly rural . As such, the majority of the inhabitants of the province are small-holder farmers whose livelihoods depend on crops (e.g. maize, wheat, tobacco) and livestock (e.g. cattle, sheep, goats and poultry). However, cattle farming experiences economic losses due to tick-borne diseases, typical to resource-poor cattle owners in Africa .
Like other provinces in Zimbabwe, Mashonaland Central Province’s tick control is based on dipping cattle in plunge dips (concrete plunge tanks) containing acaricides diluted in large volumes of water. Farmers have to dip their cattle weekly during the wet season and fortnightly during the dry season. It has been the responsibility of the government to supply acaricides. However, the communal farmer pays a reasonable fee to the Department of Veterinary Services so that their herds of cattle benefit from the tick control exercise . Dipping in acaricides has been a compulsory exercise since its legislation in 1914, and it has been one of the successful methods of tick control in livestock .
Tick species geographical data
Veterinary extension officers collected ticks at all the cattle dipping facilities for which location coordinates were known. The procedure of getting the ticks involves manually restraining every tick-infested animal at the dipping facility to allow the physical examination. The veterinary officers used the morphological identification method to identified the tick species following  plus guides including  and completed tick inspection forms using morphology. The veterinary officers are experts in tick species identification. Since cattle dipping is done weekly during the wet season and fortnightly during the dry season, veterinary officers inspected cattle for ticks and recorded tick presence data every month from 2010. The tick species collected for the 2018 year of study were Rhipicephalus decoloratus (African blue tick), Rhipicephalus sanguineus (s.l.) (brown dog tick), Rhipicephalus evertsi (red-legged tick) and Amblyomma hebraeum (bont tick).
We used bioclimatic data obtained from the WorldClim Version 2 Server for both present (2018) and future (2050) predictions . Specifically, we used the CCSM4 climate change scenario rcp45 with a spatial resolution of 2.5 min (Table 1).
We obtained the topographic and landcover variables from the ESA Climate Change Initiative  (Table 1). To model LULCC, we used the Land Change Modeler of Terrset Geospatial Monitoring and Modelling System because of its capability to incorporate Markov-Chain . To model future land use and cover change up to 2050, we used a first-order Markov-Chain, a stochastic model for quantitative land change detection . In addition, the Earth Explorer Aster Global DEM provided the elevation covariate  (Table 1).
All predictor variables were obtained in raster formats and adjusted to 1-km resolution and converted to ASCII format using ArcGIS 10.5 (Environmental Systems Research Institute, Redlands, CA, USA).
To deal with multicollinearity, we used the Variance Inflation Factor (VIF), where the covariate with a VIF > 10 was excluded in the modelling following Tagwireyi et al. . The VIF was calculated in R using USDM package . After excluding the collinear variables, we remained with ten for both the present and future modelling (Tables 2 and 3).
The Ensemble algorithm in Biomod2 package was used for the modelling in R 3.4.4 . The approach selected eight modelling techniques in fitting and averaging present and future species habitat predictions. These are generalized linear models (GLM), classification tree analysis (CTA) artificial neural networks (ANN), surface range envelope (SRE), generalized boosting model (GBM), Breiman and Cutler’s random forest for classification and regression (RF), mixture discriminant analysis (MDA) and multiple adaptive regression splines (MARS).
The choice of eight modelling techniques was based on each model's efficiency, especially in the evaluation of model response curves . The default settings in all chosen modelling techniques were used because default settings usually come enhanced for species distribution models (SDMs) . All models used a maximum of 100 iterations . Both present and future prediction models were created using 80% of occurrences and pseudo-absence data; the remaining 20% was used for model evaluation. Binary maps were then used to classify maps as suitable or unsuitable following Muhoyi et al. .
The partial receiving operating characteristic (partial ROC)  was for validating our Ensemble models by making use of area under curve (AUC) analysis and true skill statistics (TSS), which considers omission and commission errors . The TSS ranges from 0.0 to 1.0. A TSS value in range 0–0.5 was considered a poor model fit, 0.6 to 0.8 an acceptable fit and valid; any values > 0.8 were considered good to excellent . The AUC criteria considered a model with values from 0.0 to 0.6 as inferior, 0.6–0.9 as useful and > 0.9 as excellent following Ndaimani et al. . A TSS threshold of 0.6 was used to select models for the Ensemble model, following Gallien et al. .
The suitable area was calculated for both present and future predictions from binary maps using the raster zonal geometry calculation function in ArcGIS 10.5. To assess whether there is a reduction or increase in the suitable habitat, a two-sample test for equal proportions using a chi-square test was executed.
The Ensemble model achieved excellent success in predicting the current habitat of the bont tick (TSS = 0.900, ROC = 0.979). In addition, the approach was also excellent at predicting the future potential distribution of bont tick (TSS = 0.898, ROC = 0.976).
The variable importance analysis showed that elevation and maximum temperature of the warmest month, precipitation of the warmest month and precipitation of the wettest month are variables with high importance in explaining the occurrences of bont ticks (Tables 4 and 5). For example, the species habitat suitability is high in areas with moderate elevation ranging from 800 to 1000 m above sea level and high in areas with moderate temperatures ranging from 26 to 29 °C. Conversely, land cover, precipitation of coldest quarter and isothermality have the lowest contribution in determining bont tick habitat for present and future predictions.
Binary maps using the TSS = 0.6 thresholds showed the suitable habitat to be 17,008 km2 for the present model and 13,323 km2 for the future model (Fig. 2a, b), representing a decrease of bont tick potential habitat of ~ 13% (3685 km2) for Mashonaland Central Province, Zimbabwe.
The proportion of suitable habitat for the present model was 0.60, while that for the future model was 0.47. The chi-square test showed a significant difference between the present and future habitat sizes (χ2 = 3.40, df = 1, P = 2E−16).
We observed that three covariates, i.e., rainfall, elevation and temperature, are essential in explaining the predicted spatial distribution of the bont tick for both the present (2018) and future (2050). These three covariates directly affect the living conditions of tick species and tick hosts . The potential occurrences of bont ticks were high in moderate climatic conditions and moderate elevation. However, a combination of very high temperatures and high elevation and a combination of very low elevation and very low temperature did not favour the occurrence of the bont tick.
Based on the prediction results, there would be a significant ~ 13% reduction in bont tick potential habitat by 2050. The mechanisms explaining the decline include an increase in temperature, which leads to saturation deficit, a thermodynamic factor that has an inverse relationship with temperature . Ticks spend part of their life cycle stage of their time below the ground surface. Therefore, below ground conditions are critical for their survival . Bont ticks prefer moderate ground conditions, neither too cold nor too hot. As such, the ~ 13% reduction of bont tick habitat by 2050 reflects a strong association between climate change and tick phenology .
Current studies and literature suggest land cover-land use as one of the essential factors in predicting the potential distribution of tick species [1, 2]. However, in our models, land cover and land use were not essential in explaining the potential distribution of bont tick for present and future predictions. Therefore, we speculate that land use and land cover change in our study area is insignificant to influence the spatial distribution of the bont tick.
The response of probability of the presence of the bont tick against the covariates for the present and the future models had very few differences. The models showed an increase in temperature (BIO5) leading to a decrease in the probability of bont tick present, with high sensitivity at around 26 °C to 29 °C (Figs. 3 and 4). Analysis of the elevation response curve showed high sensitivity response at moderate elevation (Fig. 3), with increases in altitude leading to a low probability of bont tick presence. Since elevation and temperature are related directly, the exact mechanisms of the saturation deficit suggested by Herrmann and Gern  could apply.
The weakness of our study is onefold. Using only cattle controlling points for the tick presence data technically excludes ticks not attached to cattle in the study. Future researchers could improve the modelling process by using other data collection methods with more comprehensive sample coverage than our study. However, our results are useful because of two major strengths of our approach. The first strength is our use of authenticated bioclimatic variables that have global spatial coverage. Thus, the variables are available to other scientists to validate our findings . Second, our models maintain high validity and reliability by using the Ensemble model, which extracted the best of eight models . As such, the results of our study are essential in informing programmes that seek to control the bont tick in Mashonaland Central Province, Zimbabwe, and similar environments.
We successfully evaluated the association between environmental variables and bont tick occurrences in Mashonaland Central Province, Zimbabwe. We observed that the most critical environmental covariates which drive bont tick distribution are temperature, rainfall and elevation. We also successfully predicted both present (2018) and future (2050) bont tick potential habitat and observed that there would be a huge reduction of the potential habitat of the ticks by 2050. The study, however, could not include other tick host species in the modelling process; neither was it able to include tick samples outside tick control points. Therefore, we recommend considering tick data in all species that offer a host to bont tick as well as considering the abundance of those species in future work. This study provides baseline information for long-term tick eradication programmes.
Availability of data and materials
The Amblyomma hebraeum presence only raw data used in the study is published at https://figshare.com/articles/dataset/Climate_change_diminishes_the_potential_habitat_of_the_bont_tick_Amblyomma_hebraeum_evidence_from_Mashonaland_Central_Province_Zimbabwe_Raw_data/19816729.
Watt G, Parola P. Scrub typhus and tropical rickettsioses. Curr Opin Infect Dis. 2003;16:429–36.
Dantas-Torres F. Biology and ecology of the brown dog tick, Rhipicephalus sanguineus. Parasites Vectors. 2010;3:26.
Dabasa G, Zewdei W, Shanko T, Jilo K, Gurmesa G, Lolo G. Composition, prevalence and abundance of Ixodid cattle ticks at Ethio-Kenyan Border, Dillo district of Borana Zone, southern Ethiopia. J Vet Med Anim Health. 2017;9:204–12.
Sands AF, Apanaskevich DA, Matthee S, Horak IG, Matthee CA. The effect of host vicariance and parasite life history on the dispersal of the multi-host ectoparasite, Hyalomma truncatum. J Biogeogr. 2017;44:1124–36.
Ribeiro JM, Anderson JM, Manoukis NC, Meng Z, Francischetti IM. A further insight into the sialome of the tropical bont tick, Amblyomma variegatum. BMC Genom. 2011;12:1–11.
Randolph S. The impact of tick ecology on pathogen transmission dynamics. In: Bowman AS, Nuttall PA, editors. Ticks: biology, disease and control. Cambridge: Cambridge University Press; 2008. p. 40–72.
Kearney M, Phillips BL, Tracy CR, Christian KA, Betts G, Porter WP. Modelling species distributions without using species distributions: the cane toad in Australia under current and future climates. Ecography. 2008;314:423–34.
Otranto D, Wall R. New strategies for the control of arthropod vectors of disease in dogs and cats. MUVE. 2008;22:291–302.
Lambin EF, Tran A, Vanwambeke SO, Linard C, Soti V. Pathogenic landscapes: interactions between land, people, disease vectors, and their animal hosts. Int J Health Geogr. 2010;9:54.
Ehrmann S, Liira J, Gärtner S, Hansen K, Brunet J, Cousins SA, et al. Environmental drivers of Ixodes ricinus abundance in forest fragments of rural European landscapes. BMC Ecol. 2017;17:1–14.
Foti R, Nyakudya I, Moyo M, Chikuvire J, Mlambo N. Determinants of farmer demand for fee-for-service extension in Zimbabwe: the case of Mashonaland Central province. J Agric Educ Ext. 2007;14:95–104.
Vincent V, Thomas R. An agroecological survey of southern Rhodesia part 1: agro-ecological survey: government printers, Salisbury. 1960.
Njaya T, Mazuru N. Emerging new farming practices and their impact on the management of woodlots in A1 resettlement areas of Mashonaland Central Province in Zimbabwe. ADPR. 2014;2:1–19.
Sungirai M, Baron S, Moyo DZ, De Clercq P, Maritz-Olivier C, Madder M. Genotyping acaricide resistance profiles of Rhipicephalus microplus tick populations from communal land areas of Zimbabwe. Ticks Tick Borne Dis. 2018;9:2–9.
Norval R. Tick infestations and tick-borne diseases in Zimbabwe Rhodesia. J S Afr Vet Assoc. 1979;50:289–92.
Dantas-Torres F, Latrofa MS, Annoscia G, Giannelli A, Paris A, Otranto D. Morphological and genetic diversity of Rhipicephalus sanguineus sensu lato from the New and Old Worlds. Parasites Vectors. 2013;6:1–17.
Walke AR. Ticks of domestic animals in Africa: a guide to identification of species. Edinurgh: Bioscience Reports; 2003.
Mpakairi KS, Ndaimani H, Tagwireyi P, Zvidzai M, Madiri TH. Futuristic climate change scenario predicts a shrinking habitat for the African elephant Loxodonta africana: evidence from Hwange National Park, Zimbabwe. Eur J Wildl Res. 2020;66:1.
Blasi C, Capotorti G, Ortí MMA, Anzellotti I, Attorre F, Azzella MM, et al. Ecosystem mapping for the implementation of the European biodiversity strategy at the national level: the case of Italy. Environ Sci Policy. 2017;78:173–84.
Kumar KS, Bhaskar PU, Padmakumari K. Application of land change modeler for prediction of future land use land cover: a case study of Vijayawada City. IJATES. 2015;3:773–83.
Muller MR, Middleton J. A Markov model of land-use change dynamics in the Niagara Region, Ontario, Canada. J Landsc Ecol. 1994;9:151–7.
Faundeen JL, Kanengieter RL, Buswell MD. US geological survey spatial data access. J Geospat Eng. 2002;42:145–145.
Robinson TP, Wint GW, Conchedda G, Van Boeckel TP, Ercoli V, Palamara E, et al. Mapping the global distribution of livestock. PLoS ONE. 2014;9:e96084.
Gilbert M, Nicolas G, Cinardi G, Van Boeckel TP, Vanwambeke SO, Wint GW, et al. Global distribution data for cattle, buffaloes, horses, sheep, goats, pigs, chickens and ducks in 2010. Sci Data. 2018;5:1–11.
Tagwireyi P, Wenga T, Ndaimani H, Mpakairi KS. Environmental correlates of cheetah Acinonyx jubatus space-use in a savanna landscape. Afr J Wildl Res. 2020;50:157–66.
R Development Core Team. R: a language environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2010. http://www.r-project.org.
Thuiller W, Lafourcade B, Engler R, Araújo MB. BIOMOD—a platform for ensemble forecasting of species distributions. Ecography. 2009;32:369–73.
Ren Z, Peng H, Liu Z-W. The rapid climate change-caused dichotomy on subtropical evergreen broad-leaved forest in Yunnan: reduction in habitat diversity and increase in species diversity. Plant Divers. 2016;38:142–8.
Gomes VH, Ijff SD, Raes N, Amaral IL, Salomão RP, de Souza Coelho L, et al. Species distribution modelling: contrasting presence-only models with plot abundance data. Sci Rep. 2018;8:1–12.
Hardlife M, Henry N, Paradzayi T, Mpakairi KS, Eliah G. Predicting the invasion of a southern African savannah by the black wattle Acacia mearnsii. J For Res. 2019;31:1–9.
Peterson AT, Papeş M, Soberón J. Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecol Modell. 2008;213:63–72.
Allouche O, Tsoar A, Kadmon R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic TSS. J Appl Ecol. 2006;43:1223–32.
Morin X, Thuiller W. Comparing niche-and process-based models to reduce prediction uncertainty in species range shifts under climate change. J Ecol. 2009;90:1301–13.
Ndaimani H, Tagwireyi P, Sebele L, Madzikanda H. An ecological paradox: the African wild dog Lycaon pictus is not attracted to water points when water is scarce in Hwange National Park, Zimbabwe. PLoS ONE. 2016;11:e0146263.
Gallien L, Douzet R, Pratte S, Zimmermann NE, Thuiller W. Invasive species distribution models—how violating the equilibrium assumption can create new insights. Glob Ecol Biogeogr. 2012;21:1126–36.
Gazi U, Yapar D, Karasartova D, Gureser AS, Akdogan O, Unal O, et al. The role of Treg population in pathogenesis of Crimean Congo hemorrhagic fever. Virus Res. 2018;250:1–6.
Herrmann C, Gern L. Survival of Ixodes ricinus (Acari: Ixodidae) under challenging conditions of temperature and humidity is influenced by Borrelia burgdorferi sensu lato infection. J Med Entomol. 2014;47:1196–204.
Ostfeld RS, Brunner JL. Climate change and Ixodes tick-borne diseases of humans. Philos Trans R Soc Lond B Biol Sci. 2015;370:20140051.
Eisen RJ, Eisen L, Ogden NH, Beard CB. Linkages of weather and climate with Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae), enzootic transmission of Borrelia burgdorferi, and Lyme disease in North America. J Med Entomol. 2016;53:250–61.
We are grateful to the Veterinary extension officers who collected ticks at all the cattle dipping facilities in Mashonaland Central Province, Zimbabwe, and identified the ticks to species level. We are also thankful to the Department of Veterinary Services for sharing their data with the authors and allowing the data and the study to be published.
This study did not receive any funding.
Ethics approval and consent to participate
This study did not require ethical approval.
Consent for publication
All authors read and approved the manuscript and approved it for publication.
The three authors have no conflict of interest regarding this manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tagwireyi, P., Ndebele, M. & Chikurunhe, W. Climate change diminishes the potential habitat of the bont tick (Amblyomma hebraeum): evidence from Mashonaland Central Province, Zimbabwe. Parasites Vectors 15, 237 (2022). https://doi.org/10.1186/s13071-022-05346-z