Next Article in Journal
Spatial Resolution Enhancement of Satellite Microwave Radiometer Data with Deep Residual Convolutional Neural Network
Previous Article in Journal
Ship Detection for PolSAR Images via Task-Driven Discriminative Dictionary Learning
Previous Article in Special Issue
Application of UAV-Based Methodology for Census of an Endangered Plant Species in a Fragile Habitat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Looking for Ticks from Space: Using Remotely Sensed Spectral Diversity to Assess Amblyomma and Hyalomma Tick Abundance

1
Earth and Life Institute, Université catholique de Louvain (UCLouvain), Place Louis Pasteur 3-L4.03.0, 1348 Louvain-la-Neuve, Belgium
2
Department of Life Sciences, University of Trieste, Via L. Giorgieri 10, 34127 Trieste, Italy
3
Department of Veterinary Tropical Diseases, Faculty of Veterinary Science, University of Pretoria, Onderstepoort, Pretoria 0110, South Africa
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(7), 770; https://doi.org/10.3390/rs11070770
Submission received: 12 February 2019 / Revised: 15 March 2019 / Accepted: 27 March 2019 / Published: 30 March 2019
(This article belongs to the Special Issue Remote Sensing for Biodiversity, Ecology and Conservation)

Abstract

:
Landscape heterogeneity, as measured by the spectral diversity of satellite imagery, has the potential to provide information on the resources available within the movement capacity range of arthropod vectors, and to help predict vector abundance. The Spectral Variation Hypothesis states that higher spectral diversity is positively related to a higher number of ecological niches present in the landscape, allowing more species to coexist regardless of the taxonomic group considered. Investigating the landscape heterogeneity as a proxy of the resources available to vectors may be relevant for complex and continuous agro-forest mosaics of small farmlands and degraded forests, where land cover classification is often imprecise. In this study, we hypothesized that larger spectral diversity would be associated with higher tick abundance due to the potentially higher number of hosts in heterogeneous landscapes. Specifically, we tested whether spectral diversity indices could represent heterogeneous landscapes, and if so, whether they explain Amblyomma and Hyalomma tick abundance in Benin and inform on their habitat preferences. Benin is a West-African country characterized by a mosaic landscape of farmland and degraded forests. Our results showed that both NDVI-derived and spectral predictors are highly collinear, with NDVI-derived predictors related to vegetated land cover classes and spectral predictors correlated to mosaic landscapes. Amblyomma abundance was not related to the predictors considered. Hyalomma abundance showed positive relationships to spectral diversity indices and negative relationships to NDVI-derived-ones. Though taxa dependent, our approach showed moderate performance in terms of goodness of fit (ca. 13–20% R2), which is a promising result considering the sampling and scale limitations. Spectral diversity indices coupled with classical SRS vegetation indices could be a complementary approach for providing further ecological aspects in the field of disease biogeography.

Graphical Abstract

1. Introduction

Environmental conditions and host distribution are important factors shaping disease vector distribution [1,2,3]. This is particularly true for ticks [4,5], which are the second most important arthropod disease vector after mosquitoes [6].
To investigate vector ecological requirements, satellite remote sensing (SRS) products are commonly used as proxies of environmental variables for modelling tick distribution and population dynamics at different spatial scales [2,7]. SRS is one of the most cost-effective, comprehensive approach for complete spatial coverage of the Earth’s surface under study over short (and repeated) periods of time [8]. However, since SRS cannot directly estimate the environmental variables important for tick populations (e.g., humidity), many predictive models make use of composite indices, such as vegetation indices, that correlate with climatic variables [6]. For instance, the Normalized Difference Vegetation Index (NDVI) is considered a milestone in epidemiological studies, since the pioneering studies of Rogers [9] and Randolph [10] on the biological link between disease vector distributions and NDVI. The use of SRS has often been assumed as capable of describing species’ habitats by characterizing them from single pixels (e.g., Ogden [11]), especially for species with a short movement range. However, disease vectors and hosts experience an assemblage of habitat resources that will often span across pixels. Moreover, epidemiological studies often only use vegetation indices, thereby exploiting only a portion of the ecological information derivable from the electromagnetic spectrum [12].
The Spectral Variation Hypothesis (SVH) states that spectral diversity, the variability in the spectral response of a remotely sensed image, is positively associated with a higher number of ecological niches, which can host more species [13]. Environments characterized across pixels and spectral diversity could provide information on the resources available within the movement capacity range of the organism considered. Inter-pixel spectral variance, which is directly related to landscape heterogeneity [14,15], together with theoretical approaches such as the Resource-Based Habitat Concept (RBHC) [16] could bring an interesting perspective to vector-borne disease modelling.
RBHC is strictly related to the Hutchinsonian niche concept [17], but RBHC stresses the functional aspect of the habitat, assuming that a functional habitat is a multi-dimensional entity resulting from the overlap or contiguity of ecological units (e.g., the ecological resources, such as reproduction sites and/or suitable microclimate) required to complete the life cycle of the organism [16,18]. Generalist disease vectors such as ticks can interact with different hosts during each stage of their life cycle, showing a wide range of habitat and climate preferences [19]. As heterogeneous landscapes potentially have a higher overlap of functional habitats, landscape heterogeneity, that is, the number and proportions of different cover types [20], may contribute to shape the distribution of disease vectors. Previous research showed that Ixodes spp. density tends to be higher in fragmented landscapes compared to more homogeneous and continuous ones [21,22,23]. Fragmented areas are usually considered more heterogeneous [24,25,26,27], harboring greater variability in terms of niche availability [28,29], feeding and reproduction sites, which in turn may promote a greater abundance of hosts [16]. Generally, a highly fragmented landscape characterized by a mosaic of crops and urban areas present higher spectral diversity values compared to more homogenous landscapes within the same study area [30]. The higher the spectral diversity, the higher the number of ecological niches available, which can allow more species to coexist regardless of the taxonomic group considered [13,14,31]. This assumption also bears interest for ticks. In light of RBHC, spectral diversity may be linked to disease vector resources availability through landscape heterogeneity [16]. This is particularly relevant for complex and continuous agro-forest mosaics of farmland and degraded forest, where land cover classification is often imprecise and subject to observer bias [32,33,34,35].
We hypothesized that larger spectral diversity may be associated with higher tick abundance due to the potentially higher number of hosts in heterogeneous landscapes. Specifically, we tested whether spectral diversity indices could represent heterogeneous landscapes, and if so, whether they explain tick abundance and inform on the habitat preferences of these ticks. In particular, we studied Amblyomma and Hyalomma in Benin, a West-African country characterized by a mosaic landscape of farmland and degraded forests.

2. Materials and Methods

2.1. Study Area

The Republic of Benin (13N°, 0–4°E) is characterized by a north-south aridity gradient, ranging from arid and desert regions in the north to sub-humid areas in the south [36]. The country has three major ecological regions with three distinct rainfall patterns [32,37]. The northern region has one short rainy season between August and September, while the central region features a long rainy season from the beginning of March to the end of October. The southern coastal region has a typical sub-equatorial climate with two rainy and two dry seasons per year. The main rainy season is from April to late July, with a shorter, less intense rainy period from late September to November. The main dry season is from December to April, with a short, cooler dry season from late July to early September. Bushy vegetation suitable for grazing covers about 65% of the whole country [38], which is characterized by a fragmented agricultural landscape [39,40,41].

2.2. Sampling Protocol

Due to the lack of an exhaustive list of farms, 104 sedentary cattle herds were selected following a spatially stratified sampling scheme, accounting for administrative unit areas and cattle presence (Figure 1). Locally, the sampling was opportunistic following the suggestions of local government veterinarians and willingness to participate. In each herd, ticks were sampled from at least two domestic bovine hosts. The sampling effort was mainly addressed to the invasive cattle tick Rhipicephalus microplus, however, during the sampling all ticks were collected, and other species were also found [37]. In warm regions, tick survival is largely dependent on humidity, and tick abundance starts to decline at the end of the rainy season. As the rainfall pattern features a north–south gradient, the field missions were organized from north to south accordingly, sampling from mid-September to mid-December 2011. For a detailed description of the sampling protocol, see De Clercq et al. [37]. In the present study, we focus on two native genera found during the field campaign, namely Amblyomma and Hyalomma. Tick identification was based on morphology using a stereoscope (80-fold magnification) and a microscope (100 to 200-fold magnification). Only adult specimens were identified down to species level, when possible, using taxonomic descriptions and morphological keys [42,43,44].

2.3. Tick Autecology

Only one species from the Amblyomma genus was found (A. variegatum), which is the only indigenous Amblyomma species in West Africa. A. variegatum is a vector of Ehrlichia ruminatum, the pathogenic agent of heartwater in cattle and goats, and Anaplasma phagocytophillum, agent of anaplasmosis. A. variegatum is a generalist species found in a wide range of ecological conditions, from rainforest to savanna and steppe [42]. Hyalomma ticks are abundant in arid areas such as deserts, steppe and savanna, and play an important role in the transmission of An. phagocytophillum, Crimean-Congo hemorrhagic fever virus and rickettsia [45]. The taxonomy of Hyalomma in West Africa is not well understood yet and recent analyses have showed that the species status of some lineages could be questioned [46]. For these reasons, we treated all Hyalomma specimens as representatives of the whole genus.

2.4. Satellite Images and Predictor Variables

MOD09A1_V6 images as close as possible to the sampling period (taken on 11 December 2011) were downloaded from the MODIS website (https://ladsweb.modaps.eosdis.nasa.gov/search/, accessed February 2018). Except for a small portion of the southern part of the country, the images were cloud-free for the whole study area. The estimates between tick abundance and SRS predictors are not biased by clouds presence because clouds were not detected above the areas were the models were trained.
The MOD09A1_V6 product provides an atmospheric-corrected estimate of the surface spectral reflectance of Terra MODIS Bands 1–7 within an 8-day composite period, with a spatial resolution of 500 m. The spectral bands were used to calculate NDVI and the principal component analysis (PCA) was calculated on the seven spectral bands. The spectral landscape diversity measures were then computed on the NDVI, the first axis of the principal component analysis (PC1), and finally, on the full stack of spectral bands.
Vegetation indices highlight the properties of vegetation through a spectral transformation of two or more spectral bands [47]. The NDVI, computed at single-pixel level, is a measure of photosynthetic activity that has been widely used in spatial epidemiology as a proxy for forest cover or relative humidity within the vegetation layer [9,12,48,49].
Rocchini [30,50] summarized several approaches to measure spectral diversity and the associated landscape heterogeneity, mainly based on the spatial variation of pixel values. Image texture measurements can quantify the spatial variability in the reflectance of neighboring pixels and its spatial arrangement in a given area [51]. Textural measurements computed from SRS images can contribute to the investigation of landscape heterogeneity, adding relevant information that is not captured by single-pixel vegetation indices [52,53]. We used a grey-level co-occurrence matrix (GLCM) over all directions (0°, 45°, 90°, and 135°) as a texture measurement [54]. The GLCM used a moving window of 5 × 5 pixels, corresponding to the landscape grazed daily around each farm, within a distance of 2500 m. The GLCM counts the frequency of occurrence of two pixels of any combination of grey levels in a certain spatial relationship. A matrix is produced by considering any combination of grey levels in the image, in order to represent the spatial correlation of the grey levels, describing the image from the interval of adjacent pixels, direction and the extent of variation [55]. Haralick [54] proposed various measures of image texture based on the GLCM. Since texture measures are often collinear, following Wood [56,57], we computed three texture measures (variance, entropy and contrast) previously shown to be strongly related with foliage height diversity [56]. These three texture measures were computed on NDVI and on the scores of the first axis of PCA calculated on the seven MODIS spectral bands (PC1).
Together with texture measurements, we computed Rao’s Quadratic entropy (Rao’s Q), to measure spectral diversity [52]. This was computed both on NDVI and PC1 using a moving window of 5 × 5 pixels. Unlike indices based on entropy, such as Shannon’s H, Rao’s Q accounts for the numerical value of reflectance per se and the distances among reflectance values [52]. Rao’s Q index is computed as the sum of all the spectral distances between pairs of pixels, multiplied by the relative abundance of each type of pair of pixels in the image analyzed (see Rocchini [52] for a detailed explanation of the index).
Using the function available in Rocchini [52], we calculated multispectral Rao’s Q on the seven MODIS spectral bands within a moving window of 5 × 5 pixels. A flowchart of the analyses is shown in Figure 2.
Predictors were divided in two groups (spectral and NDVI-derived, respectively; Table 1). All herds sampled can be considered sedentary as none move across long distances or long periods. Because the majority of the grazing time is spent at less than 3 km from the farm, we arbitrarily defined a circular buffer of 2500 m radius around each sampling location as the landscape experienced by each herd, and we summarized the value of each variable in the buffer area.

2.5. Data Analysis

We computed a PCA and the Spearman’s rank correlation coefficient between the above-mentioned predictors and the GlobCover 2009 land cover map (http://due.esrin.esa.int/page_globcover.php, [58]) (Table 2), in order to assess the collinearity among the above-mentioned predictors and to investigate which land cover class is more related to spectral diversity indexes. Land cover classes are expressed as the proportions of each land cover class over the circular buffer of 2500 m radius around each sampling location.
We used a Zero-Inflated Poisson (ZIP) Generalized Linear Model (link = log) to test for association between tick abundance (response variable) and spectral diversity metrics. The choice of ZIP models was necessary to model the excess of zeros present in the response variable, which have been assumed to derive from a combination of independent processes [59]. The first process, the zero-component distribution, is modelled by a Bernoulli process where the probability of a true zero ( π i ) may depend upon environmental variables presumably leading to the difficulty in observing the phenomenon of interest, generating “excess” or “inflated” zeros observed in the data in addition to the zeros expected under the Poisson distribution. The second process is governed by a Poisson distribution that generates counts, of which some may be zero. The intensity of the Poisson process (mean number of individuals) at location i is μi. The ZIP model probability mass function was
Pr ( Y i = y ) = { π i + ( 1 π i ) exp ( μ i ) y = 0 ( 1 π i ) exp ( μ i ) μ i y / y ! y = 1 , 2 , 3 ,
Pr ( Y i = y ) = { π i + ( 1 π i ) exp ( μ i ) y = 0 ( 1 π i ) exp ( μ i ) μ i y / y ! y = 1 ,   2 ,   3 ,  
If the zero comes from the Bernoulli distribution, the observation is free from the probability of having a positive outcome. The coefficients estimated by the ZIP models were expressed as the exponential of the log coefficient. Thus, the coefficients of the zero-component distribution should be interpreted as an Odds Ratio (OR), while coefficients estimated by the count model as a Relative Risk (RR). Both OR and RR express the impact on the response variable y when increasing the independent variable x by one unit. Vuong’s non-nested hypothesis test [60] was used to compare the ZIP model with their non-zero-inflated analogs. A significant z-test indicates that the zero-inflated model should be preferred over a classical Generalized Linear Model (GLM).
Due to the presence of collinearity among the predictors (Spearman’s |rho| > 0.7; [61]), single-variable models were computed for each taxon. Models were compared using Akaike’s information criterion (AIC), where lower values indicate a better model. Statistical differences between models were based on ΔAIC scores larger than two units [62]. All the analyses were performed in R 3.4 [63] and the full code is available in the Supplementary Materials.

3. Results

The spatial pattern of abundance of the two genera differed (Figure 3). As expected, Hyalomma ssp. was found mainly in the drier areas of the north-east of the country, whereas A. variegatum was homogenously distributed over the whole country. A. variegatum was the most abundant species with an average abundance of 14 ± 23 individuals (maximum of 144 individuals), while the abundance of Hyalomma ssp. was on average 6 ± 11 individuals, with a maximum of 61.

3.1. Land Cover Analysis

The first two axes of the PCA on all variables explained 31.95% and 15.49% of the total variance, respectively (Figure 4). Hyalomma ssp. and A. variegatum had slightly different patterns and were associated with different sets of variables. Hyalomma ssp. was positively correlated to spectral predictors such as multispectral Rao’s Q and GLCM PC1 entropy, and to land cover classes 14 and 20 (Figure 4 and Figure S1). Additionally, Hyalomma spp. showed a negative correlation with GLCM PC1 variance being the spectral predictor more associated with vegetated surfaces and broadleaved forests. A. variegatum was not strongly associated to any of the variables considered and appeared orthogonal to most of them (Figure 4). A. variegatum seemed poorly correlated to the predictors considered (Spearman’s |rho| ~ 0.1, Figure S1), except for land cover class 60, which seems slightly negatively correlated. All spectral predictors except GLCM PC1 variance and NDVI-derived predictors were positively associated with anthropic mosaic landscapes (Spearman’s |rho| ~ 0.5, Figure S1), such as mosaic cropland and grassland/shrubland/forest (land cover class 20). Spectral predictors are highly correlated to each other and four main clusters are observable. Multispectral Rao’s Q, PC1 Rao’s Q, GLCM PC1 entropy and GLCM PC1 contrast are extremely correlated to each other and positively correlated with land cover classes 14 and 20, while, as expected, NDVI-derived predictors are positively correlated to two forest land cover classes (30 and 40) but they are negatively correlated to mosaic classes 20 and 14. Among them, NDVI and GLCM NDVI variance are extremely collinear (Spearman’s |rho| ~ 0.99, Figure S1).

3.2. Poisson Zero-Inflated Models

For both species, the Vuong test suggested the use of ZIP models instead of the Poisson GLMs (Table S1).
Table 3 and Table 4 report the model coefficients calculated for each taxon. Regarding A. variegatum, the ZIP models with lower AICs in each group of variables had multispectral Rao’s Q and NDVI Rao’s Q as predictors. The NDVI-based model had the overall lowest AIC score (Table 3). In the count part of the model, these two models have RR < 1, thus suggesting a negative relationship between the response variable and the predictors. However, in both cases, the binomial component of the ZIP model was not significant and the R2 calculated between the observed and predicted values was ~ 0.00, confirming that the predictors used did not catch the variability of the response variable.
Most of the ZIP models tested for Hyalomma ssp. had a significant predictor (Table 4). The models with lower AICs were the ones having GLCM PC1 variance as a predictor and those having NDVI and GLCM NDVI variance. None of the NDVI and GLCM NDVI variance models could be declared as the best model since the ΔAIC among them was lower than two. Comparing the best spectral group model with the two NDVI group models, the best spectral model has the overall lowest AIC. In the count part of the model, these three models showed a negative relationship between the response variable and the predictors (RR < 1, Table 4). In the binomial part of the model, all the estimate odds are positive, meaning that as the value of the predictors positively increases, it decreases the opportunity of having a value > 0. The R2 calculated between the observed and predicted values showed an R2 = 0.21 for the model using GLCM PC1 variance, while the two NDVI-based models have an R2 = 0.13.
The maps of the predicted value for the three best ZIP models are in Figure 5. Higher predicted Hyalomma ssp. abundance was mostly in the northern part of the country. However, the spatial pattern of predictions of the three models differed. While Figure 5a shows three main clusters of higher abundance in the uppermost part of the country and near Burkina Faso and Niger, Figure 5b, c shows lower and more homogeneously dispersed predicted abundance values. However, Figure 5b has a finer spatial pattern prediction compared to Figure 5c, which appears more smoothly distributed. False positives due to cloud presence were observed in the southern part of the map in Figure 5b,c but not in Figure 5a.

4. Discussion

In this study, we assessed the capacity of SRS spectral-based diversity metrics to represent heterogeneous landscapes and to explain tick abundance of two Western African genera of ticks, Amblyomma and Hyalomma. We framed our choice of indices within the RBHC, hypothesizing that greater landscape heterogeneity, as measured by spectral diversity, is associated to higher resource availability [16].
Spectral diversity tends to be higher in a heterogeneous landscape, such as a mosaic landscape [64,65]. According to Randolph [66], tick abundance is related to large spectral diversity among the band values, i.e., where the vegetation is very heterogeneous and diverse. The vegetation affects both the micro-climate to which ticks are exposed and the host communities that ticks feed on [19,29]. The larger the number of potential ecological niches present in heterogeneous landscapes, the higher species diversity there may be [50,67]). In mosaic croplands, land cover heterogeneity is expected to be higher, and thus likely associated with arthropod abundance and diversity [68,69,70,71]. Specifically referring to ticks, mosaic-shaped landscapes experienced higher tick densities than homogeneous forested landscapes [5]), which may be in part due to high host numbers or feeding and reproduction [16,19].
Unfortunately, the SRS predictors chosen were not able to explain the variability of A. variegatum abundance, as shown by both the PCA and modelling results, though the models computed always reported an AIC score lower than the null model.
The lack of explanatory and predictive power in the models was not totally unexpected. The use of spectral texture metrics to explain species abundance are highly taxa and scale dependent. Spectral texture metrics yielded good results for, e.g., avian abundances [56] or plants [51], but others studies have showed how this approach remains highly taxa dependent [53,72]. Our sampling may not have been synchronous everywhere with population dynamics, following more closely the dynamics of R. microplus, the target species of this sampling effort.
The best models for Hyalomma ssp. abundance highlighted GLCM PC1 variance, NDVI and GLCM NDVI variance as the best predictors. The latter two predictors being extremely correlated (Spearman’s rho ~ 0.99, Figure 2) suggests why these two NDVI derived models have a ΔAIC < 2. Overall, all three models showed a negative relationship between the response and the explanatory variables. Both the PCA and correlation matrix showed a negative correlation between these predictors and Hyalomma spp. abundance, which is indeed positively correlated to mosaic croplands land cover, and a positive correlation to more vegetated/forested land cover classes. Moreover, multispectral Rao’s Q, PC1 Rao’s Q and GLCM PC1 entropy are three predictors strongly correlated to mosaic cropland and grassland that also have significant positive ZIP estimates for Hyalomma ssp. abundances (Table 4).
Our findings support the general understanding about Hyalomma ssp. ecology. The genus is related to non-forested areas (negative relationship to forested areas) and it could be positively affected by heterogeneous landscape such as mosaic cropland and grassland, where the species may also find small mammals and birds as blood meals. Several studies pointed out how avian biodiversity and abundance is positively associated to spectral diversity in savanna landscapes (see for instance [56,57]). Interestingly, Hyalomma species [73,74], but also A. variegatum [75], have been shown to parasitize small mammals and birds in immature stages, while the adults parasitize mainly ungulates and cattle. The host preferences of both Amblyomma spp. and Hyalomma ssp. (among others) were investigated by Spengler and Estrada-Peña [76] using a network-based comparative analysis. Again, they highlighted how both genera feed on small vertebrates in the juvenile stages and mainly on large mammals as adults, however, Hyalomma ssp. feed on less phylogenetically diverse hosts than any other tick genus. It also feeds on small mammals during the adult stage. Findings by Spengler and Estrada-Peña [76] underline the importance of rural mosaic landscapes, which could provide more favorable host densities for this genus.
Though the spatial pattern of the predictions differs among the models (Figure 5), they highlight the northern drier part of the country as the areas with highest predicted abundance, according to the south-north aridity gradient and species ecology.

4.1. Ecological Perspectives of Spectral Diversity Measures

As shown in the PCA, vegetation indexes, textural indices and Rao’s Q are collinear, however, some of them showed correlation with the land cover classes considered. Ecologically speaking, NDVI remains the quickest and most easily interpretable solution to represent continuously photosynthetic active areas, without computing CPU intensive metrics such as the spectral diversity index, which may have a similar ecological meaning. However, in order to continuously represent heterogeneous landscapes, spectral diversity measures computed from the raw spectral bands present an interesting alternative to heterogeneity measures computed from class-based land cover maps. However, heterogeneous landscapes like those encountered in Benin suffer from oversimplification when investigated through land cover classes [35]. Several studies (e.g., [77]) have suggested that heterogeneous landscapes have a role in shaping tick distribution. However, mixed/mosaic environments represent habitat where various host species can live, having an ecological value that cannot be identified using a categorical or single pixel approach. Thus, continuous variables such as those derived from SRS could bring further information to monitor tick habitat (see for instance [78,79]).
The use of spectral diversity and the RBHC framework is attractive in order to better associate ecological meaning, in this case resource availability, to SRS products. Other studies (e.g., [80,81]) already acknowledged the importance of computing continuous spectral diversity measures from high spatial resolution spectral bands (e.g., Sentinel-2) in order to better understand and discriminate the various landscape components. To our knowledge, this is the first study that has explicitly investigated arthropod vectors abundance under the SVH. Though our approach showed modest results in term of goodness-of-fit between predicted and observed value, we consider the R2 = 13–20% a promising result for a single-variable model, taking into consideration the sampling and scale limitations. The spectral diversity index coupled with the classical SRS vegetation index approach could be a complementary approach to bring interesting ecological considerations into the field of disease biogeography.

4.2. Limits of the Study and Future Perspective

This study suffered from some limitations, which may have affected the results. Whereas measuring the occurrence of ticks presents us with the difficulty of distinguishing true from false absences, data on tick abundance come with additional biases. Within a herd, the two most infested animals are often selected for sampling. The number of ticks collected is thus not an indication of the mean tick load per animal in the herd, but rather an estimation of maximal tick burden.
Finer resolution of the spectral images may provide better results in predicting tick abundance. We used MODIS products (500 m spatial resolution) as finer spatial resolution products (≤30 m) were not available for 2011 (e.g., Sentinel-2) or the images were qualitatively poor in the area of interest and had a high cloud coverage (LANDSAT 5–7). Several studies have already showed how spectral diversity is sensor and scale dependent [82], thus, to overcome the scale issue, future applications of the methodology foresee the use of Sentinel-2 products and the analysis of the connectivity among pixels with higher spectral diversity [50]. As habitat diversity measures are scale-dependent, the usefulness of various SRS derived measures may vary with scale. Moreover, further applications of this methodology should harmonize the sampling period with the revisiting time of high spatial resolution satellites in order to increase the spatial and temporal consistency of the analysis.

5. Conclusions

Landscape heterogeneity is an important feature in tick ecology but it is notoriously difficult to measure in a coherent and measurable way since man-made, class-based landscapes can be irrelevant to organisms. Indeed, pixel-based remotely sensed spectral spatial diversity is a promising tool to investigate tick abundance considering a continuous measure of landscape heterogeneity.
Though taxa dependent, our approach showed a goodness-of-fit between predicted and observed values of 13–20% for Hyalomma ssp., which we consider an interesting and promising result since we modeled tick abundance using univariate models with a 500 m spatial resolution.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/11/7/770/s1, Figure S1: Spearman correlation matrix, Table S1: Vuong test results for the comparison between GLM and ZIP models. The ZIP models is preferred over the GLM, and only the p values are showed. Signif. codes: 0 ‘***’; 0.001 ‘**’; 0.01 ‘*’; 0.05 ‘.’; 0.1 ‘ ’ 1.

Author Contributions

All the authors of the present work contributed to the results and discussion of the results. D.D.R., R.R., E.M.D. and S.O.V. conceptualized the study; D.D.R., E.T., E.M.D. and S.O.V. conceived, designed, and performed the experiments; D.D.R. and E.T. analyzed the data; D.D.R. wrote the paper and all the other authors contributed to review and editing.

Funding

Tick sampling, tick identification and exploratory modelling was funded by the Belgian Science Policy Program (TickRisk project—Belspo, SR/00/144). D.D.R. is supported by the FRFS-WISD Walloon Institute for Sustainable Development PDR “Mapping livestock’s transition” (PDR-WISD X302317F).

Acknowledgments

Computational resources were provided by the supercomputing facilities of the Université catholique de Louvain (CISM/UCL) and the Consortium des Equipements de Calcul Intensif en Fédération Wallonie Bruxelles (CECI) funded by the Fond de la Recherche Scientifique de Belgique (FRS-FNRS). We acknowledge the ESA GlobCover 2009 Project for the land cover product provided.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zeman, P.; Lynen, G. Conditions for stable parapatric coexistence between Boophilusdecoloratus and B. microplus ticks: A simulation study using the competitive Lotka-Volterra model. Exp. Appl. Acarol. 2010, 52, 409–426. [Google Scholar] [CrossRef] [PubMed]
  2. Escobar, L.E.; Craft, M.E. Advances and limitations of disease biogeography using ecological niche modeling. Front. Microbiol. 2016, 7, 1174. [Google Scholar] [CrossRef] [PubMed]
  3. Medlock, J.M.; Hansford, K.M.; Bormane, A.; Derdakova, M.; Estrada-Peña, A.; George, J.C.; Golovljova, I.; Jaenson, T.G.T.; Jensen, J.K.; Jensen, P.M.; et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit. Vectors 2013, 6, 1. [Google Scholar] [CrossRef]
  4. Dautel, H.; Knülle, W. Cold hardiness, supercooling ability and causes of low-temperature mortality in the soft tick, Argas reflexus, and the hard tick, Ixodes ricinus (Acari: Ixodoidea) from Central Europe. J. Insect Physiol. 1997, 43, 843–854. [Google Scholar] [CrossRef]
  5. Pfäffle, M.; Littwin, N.; Muders, S.V.; Petney, T.N. The ecology of tick-borne diseases. Int. J. Parasitol. 2013, 43, 1059–1077. [Google Scholar] [CrossRef]
  6. Bowman, A.S.; Nuttall, P.A. Ticks: Biology, Disease and Control; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  7. Kalluri, S.; Gilruth, P.; Rogers, D.; Szczur, M. Surveillance of arthropod vector-borne infectious diseases using remote sensing techniques: A review. PLoS Pathog. 2007, 3, e116. [Google Scholar] [CrossRef]
  8. Rocchini, D.; Boyd, D.S.; Féret, J.B.; Foody, G.M.; He, K.S.; Lausch, A.; Nagendra, H.; Wegmann, M.; Pettorelli, N. Satellite remote sensing to monitor species diversity: Potential and pitfalls. Remote Sens. Ecolo. Conserv. 2016, 2, 25–36. [Google Scholar] [CrossRef]
  9. Rogers, D.J.; Hay, S.I.; & Packer, M.J. Predicting the distribution of tsetse flies in West Africa using temporal Fourier processed meteorological satellite data. Ann. Trop. Med. Parasitol. 1996, 90, 225–241. [Google Scholar] [CrossRef]
  10. Randolph, S.E. Abiotic and biotic determinants of the seasonal dynamics of the tick Rhipicephalus appendiculatus in South Africa. Med. Vet. Entomol. 1997, 11, 25–37. [Google Scholar] [CrossRef]
  11. Ogden, N.H.; Trudel, L.; Artsob, H.; Barker, I.K.; Beauchamp, G.; Charron, D.; Drebot, M.A.; Galloway, T.D.; O’Handley, R.; Thompson, R.A.; et al. Ixodes scapularis ticks collected by passive surveillance in Canada: Analysis of geographic distribution and infection with the Lyme borreliosis agent Borrelia burgdorferi. J. Med. Entomol. 2006, 43, 600–609. [Google Scholar] [CrossRef]
  12. Altobelli, A.; Boemo, B.; Mignozzi, K.; Bandi, M.; Floris, R.; Menardi, G.; Cinco, M. Spatial Lyme borreliosis risk assessment in north-eastern Italy. Int. J. Med. Microbiol. 2008, 298, 125–128. [Google Scholar] [CrossRef]
  13. Palmer, M.W.; Earls, P.G.; Hoagland, B.W.; White, P.S.; Wohlgemuth, T. Quantitative tools for perfecting species lists. Environmetrics 2002, 13, 121–137. [Google Scholar] [CrossRef]
  14. Gillespie, T.W.; Foody, G.M.; Rocchini, D.; Giorgi, A.P.; Saatchi, S. Measuring and modelling biodiversity from space. Prog. Phys. Geogr. 2008, 32, 203–221. [Google Scholar] [CrossRef] [Green Version]
  15. Skidmore, A.K.; Pettorelli, N. Agree on biodiversity metrics to track from space: Ecologists and space agencies must forge a global monitoring strategy. Nature 2015, 523, 403–406. [Google Scholar] [CrossRef] [PubMed]
  16. Hartemink, N.; Vanwambeke, S.O.; Purse, B.V.; Gilbert, M.; Van Dyck, H. Towards a resource-based habitat approach for spatial modelling of vector-borne disease risks. Biol. Rev. 2015, 90, 1151–1162. [Google Scholar] [CrossRef] [PubMed]
  17. Hutchinson, G.E. Concluding remarks. Cold Spring Harbor Symp. Quant. Biol. 1957, 22, 415–427. [Google Scholar] [CrossRef]
  18. Dennis, R.L.; Shreeve, T.G.; Van Dyck, H. Towards a functional resource-based concept for habitat: A butterfly biology viewpoint. Oikos 2003, 417–426. [Google Scholar]
  19. Jones, E.O.; Webb, S.D.; Ruiz-Fons, F.J.; Albon, S.; Gilbert, L. The effect of landscape heterogeneity and host movement on a tick-borne pathogen. Theor. Ecol. 2011, 4, 435–448. [Google Scholar] [CrossRef]
  20. Fahrig, L.; Baudry, J.; Brotons, L.; Burel, F.G.; Crist, T.O.; Fuller, R.J.; Sirami, C.; Siriwardena, G.M.; Martin, J.L. Functional landscape heterogeneity and animal biodiversity in agricultural landscapes. Ecol. Lett. 2011, 14, 101–112. [Google Scholar] [CrossRef]
  21. Li, S.; Hartemink, N.; Speybroeck, N.; Vanwambeke, S.O. Consequences of landscape fragmentation on Lyme disease risk: A cellular automata approach. PLoS ONE 2012, 7, e39612. [Google Scholar] [CrossRef]
  22. Li, S.; Heyman, P.; Cochez, C.; Simons, L.; Vanwambeke, S.O. A multi-level analysis of the relationship between environmental factors and questing Ixodes ricinus dynamics in Belgium. Parasites Vector 2012, 5, 149. [Google Scholar] [CrossRef]
  23. De Keukeleire, M.; Vanwambeke, S.O.; Somassè, E.; Kabamba, B.; Luyasu, V.; Robert, A. Scouts, forests, and ticks: Impact of landscapes on human-tick contacts. Ticks Tick Borne Dis. 2015, 6, 636–644. [Google Scholar] [CrossRef] [PubMed]
  24. Nupp, T.E.; Swihart, R.K. Effects of forest fragmentation on population attributes of white footed mice and eastern chipmunks. J. Mammol. 1998, 79, 1234–1243. [Google Scholar] [CrossRef]
  25. Allan, B.F.; Keesing, F.; Ostfeld, R.S. Effect of forest fragmentation on Lyme disease risk. Conserv. Biol. 2003, 17, 267–272. [Google Scholar] [CrossRef]
  26. Anderson, C.S.; Meikle, D.B. Annual changes in structural complexity of understory vegetation and relative abundance of Peromyscus leucopus in fragmented habitats. Acta Theriol. 2006, 51, 43–51. [Google Scholar] [CrossRef]
  27. Killilea, M.E.; Swei, A.; Lane, R.S.; Briggs, C.J.; Ostfeld, R.S. Spatial dynamics of Lyme disease: A review. EcoHealth 2008, 5, 167–195. [Google Scholar] [CrossRef]
  28. Gilbert, L.; Maffey, G.L.; Ramsay, S.L.; Hester, A.J. The effect of deer management on the abundance of Ixodes ricinus in Scotland. Ecol. Appl. 2012, 22, 658–667. [Google Scholar] [CrossRef]
  29. Gilbert, L. How landscapes shape Lyme borreliosis risk. In Ecology and Prevention of Lyme Borreliosis; Braks, M.A., van Wieren, S.E., Takken, W., Sprong, H., Eds.; Wageningen Academic Publishers: Wageningen, The Netherlands, 2016. [Google Scholar]
  30. Rocchini, D.; Bacaro, G.; Chirici, G.; Da Re, D.; Feilhauer, H.; Foody, G.M.; Gillespie, T.W.; He, K.S.; Kark, S.; Levin, N.; et al. Remotely sensed spatial heterogeneity as an exploratory tool for taxonomic and functional diversity study. Ecol. Indic. 2018, 85, 983–990. [Google Scholar] [CrossRef] [Green Version]
  31. Rocchini, D.; Delucchi, L.; Bacaro, G.; Cavallini, P.; Feilhauer, H.; Foody, G.M.; He, K.S.; Nagendra, H.; Porta, C.; Ricotta, C.; et al. Calculating landscape diversity with information-theory based indices: A GRASS GIS solution. Ecol. Inform. 2013, 17, 82–93. [Google Scholar] [CrossRef]
  32. De Clercq, E.M.; Vanwambeke, S.O. Classification of anthropogenic landscapes. In Analysis of Multitemporal Remote Sensing Images (MultiTemp). In Proceedings of the 9th International Workshop on the Analysis of Multitemporal Remote Sensing Images (MultiTemp), Brugge, Belgium, 27–29 June 2017. [Google Scholar] [CrossRef]
  33. Cherrill, A.; McClean, C. Between-observer variation in the application of a standard method of habitat mapping by environmental consultants in the UK. J. App. Ecol. 1999, 36, 989–1008. [Google Scholar] [CrossRef]
  34. Hearn, S.M.; Healey, J.R.; McDonald, M.A.; Turner, A.J.; Wong, J.L.G.; Stewart, G.B. The repeatability of vegetation classification and mapping. J. Environ. Manag. 2011, 92, 1174–1184. [Google Scholar] [CrossRef]
  35. Amici, V.; Filibeck, G.; Rocchini, D.; Geri, F.; Landi, S.; Giorgini, D.; Scoppola, A.; Chiarucci, A. Are CORINE land cover classes reliable proxies of plant species assemblages? A test in Mediterranean forest landscapes. Plant Biosyst. 2017, 152, 994–1001. [Google Scholar] [CrossRef]
  36. Guidigan, M.L.G.; Azihou, F.; Idohou, R.; Okhimamhe, A.A.; Fandohan, A.B.; Sinsin, B.; Adet, L. Modelling the current and future distribution of Kigelia africana under climate change in Benin, West Africa. Model. Earth Syst. Environ. 2018, 4, 1225–1238. [Google Scholar] [CrossRef]
  37. De Clercq, E.M.; Vanwambeke, S.O.; Sungirai, M.; Adehan, S.; Lokossou, R.; Madder, M. Geographic distribution of the invasive cattle tick Rhipicephalus microplus, a country-wide survey in Benin. Exp. Appl. Acarol. 2012, 58, 441–452. [Google Scholar] [CrossRef] [Green Version]
  38. De Clercq, E.M.; Leta, S.; Estrada-Peña, A.; Madder, M.; Adehan, S.; Vanwambeke, S.O. Species distribution modelling for Rhipicephalus microplus (Acari: Ixodidae) in Benin, West Africa: Comparing datasets and modelling algorithms. Prev. Vet. Med. 2015, 118, 8–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Vancutsem, C.; Marinho, E.; Kayitakire, F.; See, L.; Fritz, S. Harmonizing and combining existing land cover/land use datasets for cropland area monitoring at the African continental scale. Remote Sens. 2012, 5, 19–41. [Google Scholar] [CrossRef]
  40. Delrue, J.; Bydekerke, L.; Eerens, H.; Gilliams, S.; Piccard, I.; Swinnen, E. Crop mapping in countries with small-scale farming: A case study for West Shewa, Ethiopia. Int. J. Remote Sens. 2013, 34, 2566–2582. [Google Scholar] [CrossRef]
  41. Lambert, M.J.; Waldner, F.; Defourny, P. Cropland mapping over Sahelian and Sudanian agrosystems: A knowledge-based approach using PROBA-V time series at 100-m. Remote Sens. 2016, 8, 232. [Google Scholar] [CrossRef]
  42. Walker, A.R. Ticks of Domestic Animals in Africa: A Guide to Identification of Species; Bioscience Reports; Bioscience: Edinburgh, UK, 2003; pp. 3–210. [Google Scholar]
  43. Madder, M.; Adehan, S.; De Deken, R.; Adehan, R.; Lokossou, R. New foci of Rhipicephalus microplus in West Africa. Exp. Appl. Acarol. 2012, 56, 385–390. [Google Scholar] [CrossRef] [Green Version]
  44. Madder, M.; Walker, J.G.; Van Rooyen, J.; Knobel, D.; De Clercq, E.M.; Vanwambeke, S.O.; Berkvens, D. e-Surveillance in animal health: Use and evaluation of mobile tools. Parasitology 2012, 139, 1831–1842. [Google Scholar] [CrossRef]
  45. Wallménius, K.; Barboutis, C.; Fransson, T.; Jaenson, T.G.; Lindgren, P.E.; Nyström, F.; Olsen, B.; Salaneck, E.; Nilsson, K. Spotted fever Rickettsia species in Hyalomma and Ixodes ticks infesting migratory birds in the European Mediterranean area. Parasit. Vectors 2014, 7, 318. [Google Scholar] [CrossRef]
  46. Sands, A.F.; Apanaskevich, D.A.; Matthee, S.; Horak, I.G.; Harrison, A.; Karim, S.; Mohammad, M.K.; Mumcuouglu, K.Y.; Rajakaruna, R.S.; Santos-Silva, M.M.; et al. Effects of tectonics and large scale climatic changes on the evolutionary history of Hyalomma ticks. Mol. Phylogenet. Evol. 2017, 114, 153–165. [Google Scholar] [CrossRef]
  47. Hill, M.J. Vegetation index suites as indicators of vegetation state in grassland and savanna: An analysis with simulated SENTINEL 2 data for a North American transect. Remote Sens. Environ. 2013, 137, 94–111. [Google Scholar] [CrossRef]
  48. Estrada-Peña, A.; Gray, J.S.; Kahl, O.; Lane, R.S.; Nijhoff, A.M. Research on the ecology of ticks and tick-borne pathogens—methodological principles and caveats. Cell. Infect. Microbiol. 2013, 3, 29. [Google Scholar] [CrossRef]
  49. Estrada-Peña, A.; Alexander, N.; Wint, G.W. Perspectives on modelling the distribution of ticks for large areas: So far so good? Parasit. Vectors 2016, 9, 179. [Google Scholar] [CrossRef]
  50. Rocchini, D.; Balkenhol, N.; Carter, G.A.; Foody, G.M.; Gillespie, T.W.; He, K.S.; Kark, S.; Levin, N.; Lucas, K.; Levin, N.; et al. Remotely sensed spectral heterogeneity as a proxy of species diversity: Recent advances and open challenges. Ecol. Inform. 2010, 5, 318–329. [Google Scholar] [CrossRef]
  51. Hernández-Stefanoni, J.L.; Gallardo-Cruz, J.A.; Meave, J.A.; Rocchini, D.; Bello-Pineda, J.; López-Martínez, J.O. Modeling α-and β-diversity in a tropical forest from remotely sensed and spatial data. Int. J. Appl. Earth Obs. 2012, 19, 359–368. [Google Scholar] [CrossRef]
  52. Rocchini, D.; Marcantonio, M.; Ricotta, C. Measuring Rao’s Q diversity index from remote sensing: An open source solution. Ecol. Indic. 2017, 72, 234–238. [Google Scholar] [CrossRef]
  53. Wallis, C.I.; Brehm, G.; Donoso, D.A.; Fiedler, K.; Homeier, J.; Paulsch, D.; Süßenbach, D.; Tiede, Y.; Brandl, R.; Farwig, N.; et al. Remote sensing improves prediction of tropical montane species diversity but performance differs among taxa. Ecol. Indic. 2017, 83, 538–549. [Google Scholar] [CrossRef]
  54. Haralick, R.M.; Shanmugam, K. Textural features for image classification. IEEE Trans. Syst. Man. Cybern. 1973, 3, 610–621. [Google Scholar] [CrossRef]
  55. Zhang, X.; Cui, J.; Wang, W.; Lin, C. A study for texture feature extraction of high-resolution satellite images based on a direction measure and gray level co-occurrence matrix fusion algorithm. Sensors 2017, 17, 1474. [Google Scholar] [CrossRef]
  56. Wood, E.M.; Pidgeon, A.M.; Radeloff, V.C.; Keuler, N.S. Image texture as a remotely sensed measure of vegetation structure. Remote Sens. Environ. 2012, 121, 516–526. [Google Scholar] [CrossRef]
  57. Wood, E.M.; Pidgeon, A.M.; Radeloff, V.C.; Keuler, N.S. Image texture predicts avian density and species richness. PLoS ONE 2013, 8, e63211. [Google Scholar] [CrossRef] [PubMed]
  58. Arino, O.; Ramos Perez, J.J.; Kalogirou, V.; Bontemps, S.; Defourny, P.; Van Bogaert, E. Global land cover map for 2009 (GlobCover 2009). ESA & UCL. 2012. Available online: http://due.esrin.esa.int/page_globcover.php (accessed on February 2018).
  59. Lambert, D. Zero-inflated Poisson regression, with an application to defects in manufacturing. Technometrics 1992, 34, 1–14. [Google Scholar] [CrossRef]
  60. Vuong, Q.H. Likelihood ratio tests for model selection and non-nested hypotheses. Econometrica 1989, 307–333. [Google Scholar] [CrossRef]
  61. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; Garcia Marquéz, J.R.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 2013, 36, 27–46. [Google Scholar] [CrossRef]
  62. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach; Springer: New York, NY, USA, 2002. [Google Scholar]
  63. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: http://www.R-project.org/ (accessed on 1 February 2019).
  64. Rocchini, D.; Ricotta, C.; Chiarucci, A.; De Dominicis, V.; Cirillo, I.; Maccherini, S. Relating spectral and species diversity through rarefaction curves. Int. J. Remote Sens. 2009, 30, 2705–2711. [Google Scholar] [CrossRef]
  65. Rocchini, D.; Neteler, M. Spectral rank–abundance for measuring landscape diversity. Int. J. Remote Sens. 2012, 33, 4458–4470. [Google Scholar] [CrossRef] [Green Version]
  66. Randolph, S.E. Ticks and tick-borne disease systems in space and from space. Adv. Parasitol. 2000, 47, 217–243. [Google Scholar] [CrossRef]
  67. Warren, S.D.; Alt, M.; Olson, K.D.; Irl, S.D.; Steinbauer, M.J.; Jentsch, A. The relationship between the spectral diversity of satellite imagery, habitat heterogeneity, and plant species richness. Ecol. Inform. 2014, 24, 160–168. [Google Scholar] [CrossRef]
  68. Marino, P.C.; Landis, D.A. Effect of landscape structure onparasitoid diversity and parasitism in agroecosystems. Ecol. Appl. 1996, 6, 276–284. [Google Scholar] [CrossRef]
  69. Guerrini, L.; Bord, J.P.; Ducheyne, E.; Bouyer, J. Fragmentation analysis for prediction of suitable habitat for vectors: Example of riverine tsetse flies in Burkina Faso. J. Med. Entomol. 2008, 45, 1180–1186. [Google Scholar] [CrossRef] [PubMed]
  70. Gardiner, M.; Landis, D.; Gratton, C.; DiFonzo, C.; O’neal, M.; Chacon, J.M.; Wayo, M.T.; Schmidt, N.P.; Mueller, E.E.; Heimpel, G.E. Landscape diversity enhances biological control of an introduced crop pest in the north-central USA. Ecol. Appl. 2009, 19, 143–154. [Google Scholar] [CrossRef] [Green Version]
  71. Penn, H.J.; Athey, K.J.; Lee, B.D. Land cover diversity increases predator aggregation and consumption of prey. Ecol. Lett. 2017, 20, 609–618. [Google Scholar] [CrossRef]
  72. Hofmann, S.; Everaars, J.; Schweiger, O.; Frenzel, M.; Bannehr, L.; Cord, A.F. Modelling patterns of pollinator species richness and diversity using satellite image texture. PLoS ONE 2017, 12, e0185591. [Google Scholar] [CrossRef] [PubMed]
  73. Apanaskevich, D.A.; Schuster, A.L.; Horak, I.G. The genus Hyalomma: VII. Redescription of all parasitic stages of H.(Euhyalomma) dromedarii and H.(E.) schulzei (Acari: Ixodidae). J. Med. Entomol. 2008, 45, 817–831. [Google Scholar] [CrossRef]
  74. Guglielmone, A.A.; Robbins, R.G.; Apanaskevich, D.A.; Petney, T.N.; Estrada-Peña, A.; Horak, I.G. The Hard Ticks of the World; Springer: Dordrecht, The Netherlands, 2014; Volume 10, pp. 978–994. [Google Scholar]
  75. Miguel, E.; Boulinier, T.; de Garine-Wichatitsky, M.; Caron, A.; Fritz, H.; Grosbois, V. Characterising African tick communities at a wild–domestic interface using repeated sampling protocols and models. Acta Trop. 2014, 138, 5–14. [Google Scholar] [CrossRef] [Green Version]
  76. Spengler, J.R.; Estrada-Peña, A. Host preferences support the prominent role of Hyalomma ticks in the ecology of Crimean-Congo hemorrhagic fever. PLoS Negl. Trop. Dis. 2018, 12, e0006248. [Google Scholar] [CrossRef]
  77. Vanwambeke, S.O.; Li, S.; Hartemink, N.A. A resource-based habitat concept for tick-borne diseases. In Ecology and Prevention of Lyme Borreliosis; Braks, M.A., van Wieren, S.E., Takken, W., Sprong, H., Eds.; Wageningen Academic Publishers: Wageningen, The Netherlands, 2016; pp. 205–216. [Google Scholar]
  78. England, M. Understanding the Risks and Factors Associated with the Introduction of Crimean-Congo Haemorrahagic Fever Virus into Great Britain. Ph.D. Dissertation, University of Southampton, Southampton, UK, 2015. [Google Scholar]
  79. Gilot, B.; Guiguen, C.; Degeilh, B.; Doche, B.; Pichot, J.; Beaucournu, J.C. Phytoecological mapping of Ixodes ricinus as an approach to the distribution of Lyme borreliosis in France. In Lyme Borreliosis; Axford, J.S., Rees, D.H.E., Eds.; NATO ASI Series (Series A: Life Sciences); Springer: Boston, MA, USA, 1994; Volume 260, pp. 105–112. [Google Scholar]
  80. Karlson, M.; Ostwald, M.; Reese, H.; Sanou, J.; Tankoano, B.; Mattsson, E. Mapping TreeCanopy Cover and Aboveground Biomass in Sudano-Sahelian Woodlands Using Landsat 8 and Random Forest. Remote Sens. 2015, 7, 10017–10041. [Google Scholar] [CrossRef]
  81. Godinho, S.; Guiomar, N.; Gil, A. Estimating tree canopy cover percentage in a mediterranean silvopastoral systems using Sentinel-2A imagery and the stochastic gradient boosting algorithm. Int. J. Remote Sens. 2018, 39, 4640–4662. [Google Scholar] [CrossRef]
  82. Rocchini, D. Effects of spatial and spectral resolution in estimating ecosystem α-diversity by satellite imagery. Remot. Sen. Environ. 2007, 111, 423–434. [Google Scholar] [CrossRef]
Figure 1. Locations of the herds that were sampled for the country-wide survey of Benin.
Figure 1. Locations of the herds that were sampled for the country-wide survey of Benin.
Remotesensing 11 00770 g001
Figure 2. Spectral predictors computation flowchart.
Figure 2. Spectral predictors computation flowchart.
Remotesensing 11 00770 g002
Figure 3. Tick abundance maps: (a) Amblyomma variegatum; (b) Hyalomma ssp.
Figure 3. Tick abundance maps: (a) Amblyomma variegatum; (b) Hyalomma ssp.
Remotesensing 11 00770 g003
Figure 4. Principal component analysis (PCA) biplot. The variables are species abundance, satellite remote sensing (SRS) predictors and land use classes relative proportion.
Figure 4. Principal component analysis (PCA) biplot. The variables are species abundance, satellite remote sensing (SRS) predictors and land use classes relative proportion.
Remotesensing 11 00770 g004
Figure 5. Map of the predicted abundance value for the three best models for Hyalommma ssp. ZIP model using as predictor: (a) GLCM PC1 variance, (b) NDVI and (c) GLCM NDVI variance.
Figure 5. Map of the predicted abundance value for the three best models for Hyalommma ssp. ZIP model using as predictor: (a) GLCM PC1 variance, (b) NDVI and (c) GLCM NDVI variance.
Remotesensing 11 00770 g005
Table 1. List of predictors used to model tick abundance.
Table 1. List of predictors used to model tick abundance.
PredictorsAbbreviation
Spectral predictors
Multispectral Rao’s Qrao_ms
Rao’s Q PC1rao_pc1
GLCM PC1 variancepc1_var
GLCM PC1 entropypc1_entr
GLCM PC1 contrastpc1_contr
NDVI
NDVIndvi
Rao’s Q NDVIrao_ndvi
GLCM NDVI variancendvi_var
GLCM NDVI entropyndvi_entr
GLCM NDVI contrastndvi_contr
Table 2. GlobCover 2009 legend. The percentage is the proportion of land cover in the total area considered (2041 km2).
Table 2. GlobCover 2009 legend. The percentage is the proportion of land cover in the total area considered (2041 km2).
ValueLabel%
14Rainfed croplands2.59
20Mosaic cropland (50–70%)/vegetation (grassland/shrubland/forest) (20–50%)28.54
30Mosaic vegetation (grassland/shrubland/forest) (50–70%)/cropland (20–50%)10.95
40Closed to open (>15%) broadleaved evergreen or semi-deciduous forest (>5 m)1.44
60Open (15–40%) broadleaved deciduous forest/woodland (>5 m)9.67
110Mosaic forest or shrubland (50–70%)/grassland (20–50%)2.68
120Mosaic grassland (50–70%)/forest or shrubland (20–50%)0.29
130Closed to open (>15%) (broadleaved or needleleaved, evergreen or deciduous) shrubland (<5 m)43.71
140Closed to open (>15%) herbaceous vegetation (grassland, savannas or lichens/mosses)2.68
Table 3. Zero-inflated Poisson (ZIP) model estimates for A. variegatum. The estimates are expressed as the exponent of the log coefficient of the model. The count part of the model should be interpreted as a relative risk while the binomial part as an odds ratio. P < 0.001 ‘***’; P < 0.01 ‘**’; P < 0.05 ‘*’; ‘ns’ not significant. Multispectral Rao’s Q: rao_ms; Rao’s Q PC1: rao_pc1; GLCM PC1 variance: pc1_var; GLCM PC1 entropy: pc1_entr; GLCM PC1 contrast: pc1_contr; NDVI: ndvi; Rao’s Q NDVI: rao_ndvi; GLCM NDVI variance: ndvi_var; GLCM NDVI entropy: ndvi_entr; GLCM NDVI contrast: ndvi_contr.
Table 3. Zero-inflated Poisson (ZIP) model estimates for A. variegatum. The estimates are expressed as the exponent of the log coefficient of the model. The count part of the model should be interpreted as a relative risk while the binomial part as an odds ratio. P < 0.001 ‘***’; P < 0.01 ‘**’; P < 0.05 ‘*’; ‘ns’ not significant. Multispectral Rao’s Q: rao_ms; Rao’s Q PC1: rao_pc1; GLCM PC1 variance: pc1_var; GLCM PC1 entropy: pc1_entr; GLCM PC1 contrast: pc1_contr; NDVI: ndvi; Rao’s Q NDVI: rao_ndvi; GLCM NDVI variance: ndvi_var; GLCM NDVI entropy: ndvi_entr; GLCM NDVI contrast: ndvi_contr.
Spectral PredictorNDVI Predictors
PredictorsCount modelZI BinomialAICPredictorsCount modelZI BinomialAIC
Estimate 2.5% CI97.5% CIp.valEstimate 2.5% CI97.5% CIp.val Estimate 2.5% CI97.5% CIp.valEstimate 2.5% CI97.5% CIp.val
Intercept16.3215.5117.17***0.140.080.26***2454.71----------
Intercept16.1915.3817.04***0.130.070.25***2403.68Intercept16.3215.5117.17***0.140.080.25***2458.56
Multispectral Rao’s Q0.830.780.87***0.610.311.18nsNDVI1.000.951.05ns1.120.632.00ns
Intercept16.2815.4717.13***0.130.070.25***2432.81Intercept15.3614.5516.22***0.140.080.26***2303.15
Rao’s Q PC10.880.840.93***0.670.351.28nsRao’s Q NDVI0.680.630.72***0.950.521.74ns
Intercept16.3215.5117.17***0.120.060.24***2454.96Intercept16.3215.5117.17***0.140.080.26***2458.53
GLCM PC1 variance1.000.951.05ns1.950.934.10nsGLCM NDVI variance0.990.951.04ns1.100.621.97ns
Intercept16.3215.5117.17***0.140.080.25***2450.14Intercept15.9715.1616.82***0.140.080.25***2373.31
GLCM PC1 entropy0.930.890.98**0.740.421.30nsGLCM NDVI entropy0.780.730.82***0.780.421.43ns
Intercept16.2515.4417.10***0.130.070.25***2422.28Intercept15.9615.1516.82***0.140.080.25ns2382.06
GLCM PC1 contrast0.850.810.90***0.570.251.28nsGLCM NDVI contrast0.770.720.82***0.740.351.54***
Table 4. ZIP model estimates for Hyalomma ssp. The estimates are expressed as the exponent of the log coefficient of the model. The count part of the model should be interpreted as a relative risk while the binomial part as an odds ratio. Signif. codes: P < 0.001 ‘***’; P < 0.01 ‘**’; P < 0.05 ‘*’; ‘ns’ not significant. Multispectral Rao’s Q: rao_ms; Rao’s Q PC1: rao_pc1; GLCM PC1 variance: pc1_var; GLCM PC1 entropy: pc1_entr; GLCM PC1 contrast: pc1_contr; NDVI: ndvi; Rao’s Q NDVI: rao_ndvi; GLCM NDVI variance: ndvi_var; GLCM NDVI entropy: ndvi_entr; GLCM NDVI contrast: ndvi_contr.
Table 4. ZIP model estimates for Hyalomma ssp. The estimates are expressed as the exponent of the log coefficient of the model. The count part of the model should be interpreted as a relative risk while the binomial part as an odds ratio. Signif. codes: P < 0.001 ‘***’; P < 0.01 ‘**’; P < 0.05 ‘*’; ‘ns’ not significant. Multispectral Rao’s Q: rao_ms; Rao’s Q PC1: rao_pc1; GLCM PC1 variance: pc1_var; GLCM PC1 entropy: pc1_entr; GLCM PC1 contrast: pc1_contr; NDVI: ndvi; Rao’s Q NDVI: rao_ndvi; GLCM NDVI variance: ndvi_var; GLCM NDVI entropy: ndvi_entr; GLCM NDVI contrast: ndvi_contr.
Spectral PredictorsNDVI Predictors
PredictorsCount ModelZI BinomialAICPredictorsCount modelZI BinomialAIC
Estimate 2.5% CI97.5% CIp.valEstimate 2.5% CI97.5% CIp.val Estimate 2.5% CI97.5% CIp.valEstimate 2.5% CI97.5% CIp.val
Intercept14.9613.8716.13***1.310.891.93ns785.47----------
Intercept14.3913.2215.66***1.350.882.05ns767.84Intercept14.0212.7715.39***1.450.922.27ns757.86
Multispectral Rao’s Q1.081.011.16*0.390.240.65***NDVI0.890.820.98*3.271.915.61***
Intercept14.4113.2115.71***1.350.882.08ns765.24Intercept13.5912.4314.87***1.360.912.05ns760.05
Rao’s Q PC11.071.001.16ns0.350.210.59***Rao’s Q NDVI0.760.670.85***1.861.133.07*
Intercept12.7011.5014.02***1.330.842.10***724.21Intercept13.8912.6415.27***1.470.932.30ns756.63
GLCM PC1 variance0.790.730.86***4.102.247.50nsGLCM NDVI variance0.880.800.96**3.281.905.66***
Intercept13.7312.4715.11***1.460.942.28ns756.78Intercept14.7513.6015.99***1.340.902.01ns781.49
GLCM PC1 entropy1.171.061.29**0.320.190.55***GLCM NDVI entropy0.960.871.04ns1.731.132.66*
Intercept14.6913.5115.96***1.280.831.97ns767.07Intercept14.6513.5315.87***1.340.901.99ns782.91
GLCM PC1 contrast1.040.971.11ns0.330.190.58***GLCM NDVI contrast0.910.831.01ns1.510.942.42ns

Share and Cite

MDPI and ACS Style

Da Re, D.; De Clercq, E.M.; Tordoni, E.; Madder, M.; Rousseau, R.; Vanwambeke, S.O. Looking for Ticks from Space: Using Remotely Sensed Spectral Diversity to Assess Amblyomma and Hyalomma Tick Abundance. Remote Sens. 2019, 11, 770. https://doi.org/10.3390/rs11070770

AMA Style

Da Re D, De Clercq EM, Tordoni E, Madder M, Rousseau R, Vanwambeke SO. Looking for Ticks from Space: Using Remotely Sensed Spectral Diversity to Assess Amblyomma and Hyalomma Tick Abundance. Remote Sensing. 2019; 11(7):770. https://doi.org/10.3390/rs11070770

Chicago/Turabian Style

Da Re, Daniele, Eva M. De Clercq, Enrico Tordoni, Maxime Madder, Raphaël Rousseau, and Sophie O. Vanwambeke. 2019. "Looking for Ticks from Space: Using Remotely Sensed Spectral Diversity to Assess Amblyomma and Hyalomma Tick Abundance" Remote Sensing 11, no. 7: 770. https://doi.org/10.3390/rs11070770

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop