# Spatio-temporal optimization of sampling for bluetongue vectors (*Culicoides*) near grazing livestock

- Carsten Kirkeby
^{1}Email author, - Anders Stockmarr
^{1, 2}, - René Bødker
^{1}and - Peter Lind
^{1}

**6**:151

https://doi.org/10.1186/1756-3305-6-151

© Kirkeby et al.; licensee BioMed Central Ltd. 2013

**Received: **5 February 2013

**Accepted: **12 May 2013

**Published: **24 May 2013

## Abstract

### Background

Estimating the abundance of *Culicoides* using light traps is influenced by a large variation in abundance in time and place. This study investigates the optimal trapping strategy to estimate the abundance or presence/absence of *Culicoides* on a field with grazing animals. We used 45 light traps to sample specimens from the *Culicoides* obsoletus species complex on a 14 hectare field during 16 nights in 2009.

### Findings

The large number of traps and catch nights enabled us to simulate a series of samples consisting of different numbers of traps (1-15) on each night. We also varied the number of catch nights when simulating the sampling, and sampled with increasing minimum distances between traps. We used resampling to generate a distribution of different mean and median abundance in each sample. Finally, we used the hypergeometric distribution to estimate the probability of falsely detecting absence of vectors on the field. The variation in the estimated abundance decreased steeply when using up to six traps, and was less pronounced when using more traps, although no clear cutoff was found.

### Conclusions

Despite spatial clustering in vector abundance, we found no effect of increasing the distance between traps. We found that 18 traps were generally required to reach 90% probability of a true positive catch when sampling just one night. But when sampling over two nights the same probability level was obtained with just three traps per night. The results are useful for the design of vector monitoring programmes on fields with grazing animals.

## Keywords

*Culicoides*obsoletusSpatial variationLight trapsAbundanceBluetongue

## Findings

Estimates of the abundance and presence/absence of vectors are essential for understanding and modeling vector-borne diseases. Light trapping is the most widely used method to sample *Culicoides* and single light traps are often assumed to be representative for abundance in a large area (e.g. [1–3]). In a previous study the number of *Culicoides* trapped in 45 CDC light traps placed in a regular grid were found to vary dramatically [4]. Spatial clusters of higher abundance were found randomly on the field, causing up to 11 times higher abundance. Huge variations in abundance between catch nights was also observed [4]. This variation can potentially have a great impact on field studies where one or a few light traps are assumed to represent the mean vector abundance on a field with grazing animals such as sheep, cattle or horses which are not confined at night (e.g. [5–9]). In this paper we analyze the data from the study from a practical perspective, quantifying the impact of spatial and temporal variation and clustering on the uncertainty of the estimated vector abundance.

We used 45 battery-operated CDC 4 W light traps (http://www.johnwhock.com) on a field (length: 750 m, width: 250 m) near Vallø, Denmark, during the summer (July - September) of 2009. A detailed description of the study is given in Kirkeby et al. [4]. The traps were evenly dispersed throughout the field in 50 by 50 m grid points [4]. 260 sheep were confined to an enclosure in one end of the field during the study. Light traps turned automatically on at dusk and off at dawn and were emptied daily. Only females of the Obsoletus group were included in the dataset. For eight catch nights half of the traps were sorted and counted due to time constraints, resulting in a 100 by 100 m grid. We carried out a series of analyses to investigate the sampling variation of *Culicoides* abundance caused by light trapping. All statistical calculations were carried out in R (R 2.12.2).

A first analysis was carried out to quantify the sampling variation on a single catch night, and was repeated for all nights. For each catch night we generated 10,000 random samples using n (1-15) traps per sample, and we then calculated the mean number of vectors in each sample, resulting in a distribution of mean vector abundances for each sample size.

A second analysis was performed to investigate if a minimum distance between the traps in a sample would improve the mean abundance estimate. We resampled one to ten traps 10,000 times using two different minimum distances between the traps: 100 m and 150 m. We used the variation in mean abundance estimates to evaluate if larger distance between traps in a sample resulted in less variance. This was carried out on data from the nights of July 21^{
st
} and August 31^{
st
} where strong clusters were found on the field [4].

In a third analysis we quantified the probability of getting a false negative result (falsely detecting absence), using a given number of traps (1-20) per sample on the field. The proportion of negative traps per night, weighted with the number of analyzed traps per night, was used as the general probability for a negative trap during the whole study period. The hypergeometric distribution was then used to calculate the probability of detecting false absence in a sample consisting of n (1-20) traps. We also removed five catch nights with more than 90% negative samples and repeated the calculations with this modified dataset.

Finally, a fourth analysis was carried out to identify the number of traps necessary to detect presence of vectors when using from one to ten traps on one, two or three randomly selected nights. Using the hypergeometric distribution as above, we determined the number of traps required for detecting the presence of vectors with 90% or 95% certainty. The number of traps required to reach the certainty levels when sampling on two and three nights was calculated by exponentiating the probabilities for one catch night to the power of two and three. This procedure was also repeated with the modified dataset. For each catch night we also calculated the probability of falsely detecting absence of vectors using one to ten traps, by exponentiating the probabilities as above. No sheep were harmed in this study. Permission to move the sheep was obtained by the owners.

*C. obsoletus*group were counted (Table 1). Out of total 530 samples, 224 (42%) were negative for

*C. obsoletus*and 306 (58%) were positive.

**Descriptive statistics for each catch night: The number of Obsoletus group specimens caught, the mean catch per trap, the number of analyzed traps, the percentage of zero-catches, the minimum catch and the maximum catch**

Date | 20.07 | 21.07 | 27.07 | 03.08 | 04.08 | 06.08 | 17.08 | 18.08 | 21.08 | 24.08 | 25.08 | 27.08 | 28.08 | 31.08 | 03.09 | 04.09 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|

| 4 | 872 | 316 | 173 | 522 | 612 | 2 | 93 | 95 | 29 | 427 | 1086 | 1 | 253 | 2 | 1 |

| 0.08 | 19.38 | 14.36 | 6.92 | 11.86 | 13.91 | 0.08 | 2.02 | 4.52 | 1.31 | 18.56 | 72.40 | 0.04 | 5.62 | 0.04 | 0.02 |

| 45 | 45 | 22 | 25 | 44 | 44 | 23 | 46 | 21 | 22 | 23 | 15 | 23 | 45 | 42 | 45 |

| 91 | 0 | 5 | 20 | 0 | 0 | 91 | 35 | 14 | 73 | 0 | 0 | 96 | 33 | 95 | 98 |

| 0 | 2 | 0 | 0 | 1 | 2 | 0 | 0 | 0 | 0 | 2 | 24 | 0 | 0 | 0 | 0 |

| 1 | 79 | 68 | 106 | 79 | 48 | 1 | 20 | 18 | 12 | 58 | 176 | 1 | 44 | 1 | 1 |

The second analysis examined whether a more precise abundance estimate may be obtained with a larger distance between the traps in each sample. There was only a marginal decrease in the variance of the abundance estimate when sampling with one to ten traps using a minimum 100 and 150 m distance between traps (data not shown).

**Mean probabilities for a false negative result, depending of the number of traps**

Traps | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|

| 0.42 | 0.34 | 0.30 | 0.28 | 0.26 | 0.24 | 0.22 | 0.21 | 0.19 | 0.16 |

| 0.16 | 0.06 | 0.03 | 0.02 | 0.01 | 0.01 | 0.00 | 0.00 | 0.00 | 0.00 |

**Number of traps needed to reach a certainty of 90% or 95% of excluding a false negative result when sampling one, two or three nights**

Full datasetProbability / Nights | 1 | 2 | 3 |
---|---|---|---|

| 18 | 3 | 1 |

| 25 | 7 | 2 |

| |||

| 1 | 2 | 3 |

| 2 | 1 | 1 |

| 3 | 1 | 1 |

Light traps are the preferred tool for estimating *Culicoides* abundance, and often only a single trap is used. However, the large spatial variation and clustering observed in *Culicoides* abundance within a grazed field may produce large uncertainties when sampling with few traps [4]. In this study, we resampled a unique dataset to quantify this variation.

Unexpectedly, we found no effect of spatial separation of traps in a sample within the range investigated. However, the estimate of vector abundance was improved when separating the trap catches temporally, due to a large temporal variation. This yields a high probability of falsely detecting absence on some catch nights (Figure 2 left). We therefore suggest avoiding catch nights with low vector activity by using e.g. weather forecasts [10], which will increase the certainty of detecting vectors if they are present in an area.

The abundance pattern on the field was measured with CDC light traps, but the result is applicable to studies with other light traps, such as the often used Onderstepoort trap, because the mechanisms are the same. However, because the range of attraction for the Onderstepoort trap is larger than for the CDC trap, less traps may be needed to properly cover the field.

There are four main concerns about the dataset used. Firstly, there was spatial autocorrelation (clusters) in the data, and the aim of the present analysis was to address the practical implications of this effect. Secondly, we are resampling from a limited number of traps. This will cause the variation in samples with many traps to decrease more than if we were not restricted to a fixed number of traps. Thirdly, it has previously been shown that the maximum range of attraction for the CDC traps is 15.25 (12.7-18.3) m, and thus the catch areas did not overlap in this study [11]. A fourth concern is that the traps can compete with each other: specimens that are caught in one trap cannot be caught in another trap. Rigot et al. [12] investigated the competition between the more powerful Onderstepoort 8 W traps and found a statistically significant competition between the traps when placed 50 m apart, but not when placed 100 m apart. The CDC 4 W traps in this will compete less and therefore we do not consider competition a problem in the present study. Furthermore, if there was important competition in the present study, it is most likely that it would have caused a depletion of specimens in the middle of the field, which we did not find (results not shown).

To our knowledge, this is the first time that empirical data has been used to quantify the optimal sampling strategy for *Culicoides* vectors within a single field. The present study represents a worst case scenario for presence/absence studies, i.e. where vectors are present, but in low numbers and are therefore difficult to detect.

To conclude, we suggest using more than one trap for sampling a field of this size, and that the sampling is repeated over more than one catch night to obtain a more precise estimate of abundance.

## Declarations

### Acknowledgements

The authors would like to thank the shepherds Frank Hansen and Rune Ploug from ValløLam for their great help with the experimental setup. We also thank Simon Haarder and Peter Iversen for their help with microscopy, and Lasse Engbo Christiansen (DTU Compute) for comments. This study was partially funded by the Danish Ministry of Food, Agriculture and Fisheries and by EU grant GOCE-2003-010284 EDENext and is catalogued by the EDENext Steering Committee as EDENext 126 (http://www.edenext.eu). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission.

## Authors’ Affiliations

## References

- Meiswinkel R, Baldet T, De Deken R, Takken W, Delécolle J, Mellor P: The 2006 outbreak of bluetongue in northern Europe, the entomological perspective. Prev Vet Med. 2008, 87: 55-63. 10.1016/j.prevetmed.2008.06.005.View ArticlePubMedGoogle Scholar
- Viennet E, Garros C, Lancelot R, Allène X, Gardès L, Rakotoarivony I, Crochet D, Delécolle J, Moulia C, Baldet T: Assessment of vector/host contact: comparison of animal-baited traps and UV-light/suction trap for collecting
*Culicoides*biting midges (Diptera Ceratopogonidae), vectors of Orbiviruses. Parasites Vectors. 2011, 4: 1-12. 10.1186/1756-3305-4-1.View ArticleGoogle Scholar - Purse B, Falconer D, Sullivan M, Carpenter S, Mellor P, Piertney S, Albon S, Gunn G, Blackwell A: Impacts of climate, host and landscape factors on
*Culicoides*species in Scotland. Med Vet Entomol. 2011, 26 (2): 168-177.View ArticlePubMedGoogle Scholar - Kirkeby C, Bødker R, Stockmarr A, Lind P: Spatial abundance and clustering of
*Culicoides*(Diptera: Ceratopogonidae) on a local scale. Parasites Vectors. 2013, 6 (43): 1-14.Google Scholar - Mullens BA, Gerry AC, Monteys VSI, Pinna M, Gonzalez A: Field studies on
*Culicoides*(Diptera: Ceratopogonidae) activity and response to deltamethrin applications to sheep in Northeastern Spain. J Med Entomol. 2010, 47: 106-110. 10.1603/033.047.0115.View ArticlePubMedGoogle Scholar - Foxi C, Delrio G: Larval habitats and seasonal abundance of
*Culicoides*biting midges found in association with sheep in northern Sardinia, Italy. Med Vet Entomol. 2010, 24 (2): 199-209. 10.1111/j.1365-2915.2010.00861.x.View ArticlePubMedGoogle Scholar - Garcia-Saenz A, McCarter P, Baylis M: The influence of host number on the attraction of biting midges,
*Culicoides*spp., to light traps. Med Vet Entomol. 2011, 25: 113-5. 10.1111/j.1365-2915.2010.00904.x.View ArticlePubMedGoogle Scholar - Romón P, Higuera M, Delécolle J, Baldet T, Aduriz G, Goldarazena A: Phenology and attraction of potential
*Culicoides*vectors of bluetongue virus in Basque Country (northern Spain). Vet Parasitol. 2011, 186 (3): 415-424.PubMedGoogle Scholar - González M, Baldet T, Delécolle JC, López S, Romón P, Goldarazena A: Monitoring of
*Culicoides*Latreille (Diptera Ceratopogonidae) After Btv Outbreaks, in Sheep Farms and Natural Habitats from the Basque Country (Northern Spain). Proc Entomol Soc Washington. 2013, 115: 48-69. 10.4289/0013-8797.115.1.48.View ArticleGoogle Scholar - Sanders C, Shortall C, Gubbins S, Burgin L, Gloster J, Harrington R, Reynolds D, Mellor P, Carpenter S: Influence of season and meteorological parameters on flight activity of
*Culicoides*biting midges. J Appl Ecol. 2011, 48: 309-368.View ArticleGoogle Scholar - Kirkeby C, Græsbøll K, Stockmarr A, Christiansen L, Bødker R: The range of attraction for light traps catching
*Culicoides*biting midges (Diptera: Ceratopogonidae). Parasites Vectors. 2013, 6 (6): 1-11.Google Scholar - Rigot T, Gilbert M: Quantifying the spatial dependence of
*Culicoides*midge samples collected by Onderstepoort-type blacklight traps: an experimental approach to infer the range of attraction of light traps. Med Vet Entomol. 2011, 26 (2): 152-161.View ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.