Temporal and spatial stability of Anopheles gambiae larval habitat distribution in Western Kenya highlands
© Li et al. 2009
Received: 12 October 2009
Accepted: 18 December 2009
Published: 18 December 2009
Localized mosquito larval habitat management and the use of larvicides have been proposed as important control tools in integrated malaria vector management programs. In order to optimize the utility of these tools, detailed knowledge of the spatial distribution patterns of mosquito larval habitats is crucial. However, the spatial and temporal changes of habitat distribution patterns under different climatic conditions are rarely quantified and their implications to larval control are unknown.
Using larval habitat data collected in western Kenya highlands during both dry and rainy seasons of 2003-2005, this study analyzed the seasonal and inter-annual changes in the spatial patterns in mosquito larval habitat distributions. We found that the spatial patterns of larval habitats had significant temporal variability both seasonally and inter-annually.
The pattern of larval habitats is extremely important to the epidemiology of malaria because it results in spatial heterogeneity in the adult mosquito population and, subsequently, the spatial distribution of clinical malaria cases. Results from this study suggest that larval habitat management activities need to consider the dynamic nature of malaria vector habitats.
A series of malaria outbreaks with high case-fatality rates has occurred in the African highlands in the past two decades [1–3]. The severe malaria situation in these highlands coupled with the spread of drug resistant parasites calls for the rapid implementation of effective malaria control programs. Currently, the use of bednets treated with pyrethroids is the major malaria control approach in Africa [4, 5]. However, there are many limitations with this approach. First, bednets require a continuous retreatment with pyrethroids and regular replacement, which may be economically prohibitive for some Africans . Second, bednets only protect users while they are under it, thereby resulting in reduced, but not eliminated, human-vector contact rates . Third, bednet use imposes strong selection pressures for mosquito vectors to develop resistance to the insecticides . It would therefore be prudent to develop alternative and more sustainable vector control measures, such as environmental management of larval habitat resources. Environmentally-friendly biological larvicides have shown some success in reducing malaria transmission [9, 10]. Before this tool can become a widespread practice, however, it needs to be made more cost-effective. A key component of this would be improvements in the understanding of the spatial and temporal distribution of mosquito larval habitats in order to enhance our efforts to target control at key focal breeding sites [11, 12].
Our current knowledge on the spatial and temporal distribution of mosquito larval habitats can be improved in the following two ways. First, spatial and environmental patterns of larval habitats could be quantified to provide more accurate guidance to larval habitat management. Second, rather than basing larval habitat distribution projections on a snap-shot of data acquired under one particular climatic condition, the extent to which data can be extrapolated to other seasons or years with different climatic conditions need to be explored [13–15]. The objective of this study is therefore to quantify seasonal and inter-annual variations in larval habitat distribution patterns.
Mosquito larval data
Anopheles gambiae sensu stricto is the primary malaria vector species in western Kenya highlands, constituting > 95% of the indoor-resting adult vector population and about 60-70% of the larval population . The number of malaria cases increases dramatically from the dry season to the rainy season . February and May are the most representative months of dry and rainy seasons, respectively. Thorough searches of all identifiable aquatic habitats (e.g., animal footprints and water ponds) were conducted in the first two weeks of February and May in 2003, 2004, and 2005, generating a total of six point maps of aquatic habitat locations. Each identified aquatic habitat was dipped up to 20 times with a standard 350 ml dipper to collect water samples. When a site was too small to make 20 dips, water was dipped as many times as possible. The larval occurrence and species in the water samples was examined. Depending on whether a given aquatic habitat contained An. gambiae larvae or not, it was recorded as either an anopheline-positive or negative habitat. The location and elevation of each habitat was recorded using the global positioning system (GPS). This study focuses on the presence of An. gambiae larvae rather than their abundance in an aquatic environment because quantification of larval abundance is prone to sampling errors, particularly in large aquatic environments.
Precipitation has long been considered as the primary determinant of the spatial variation in larval habitats in western Kenya. The precipitation data of Kakamega for the period of 2003, 2004, and 2005 were obtained from the Kenya Meteorological Service. Although mosquito development is also affected by other climatic factors such as temperature, preliminary analysis based on our daily records revealed that inter-annual variations in temperature in the study area were statistically insignificant. We therefore excluded temperature as an explanatory factor in this study.
Analysis of the temporal variation in the larval habitat locations
The temporal variations in mosquito larval habitat locations were analyzed at two levels: (1) the recurrence of larval habitats in the same locations, and (2) occurrence of habitats at nearby locations. The analysis at these two levels can help answer two key questions: (1) how likely it is to find habitats at locations where habitats were previously found, and (2) to what are extent are new habitats spatially associated with habitats that have been found previously (where overlapping habitats are excluded)? These two questions are critical to mosquito habitat management efforts. If larval positive habitats tend to occur in the same or very proximal locations between seasons/years, resource management can be optimized much more easily. For the first level analysis, each of the six point maps was converted into a raster map, in which each habitat location is represented by a pixel and overlapping analysis was selected to determine the changes in habitat pixels between time points. Based on the field observations, the diameter of habitats ranged from 4 m to 15 m. To take into account effects of the varying size of habitats, the pixel sizes range from 4 m to 20 m (allowing for a 5 m GPS error), with a 2 m interval (half of the smallest diameter of habitats). The smallest and largest pixel sizes were chosen based on the observed diameters of the smallest and largest habitats in the field. Even though the pixel sizes were carefully selected, this overlapping approach still has a limitation: the overlapping area is expected to increase with increasing pixel size. To obtain a realistic estimation of the overlapping area, we plotted the percentages of overlapping aquatic habitats as a function of the size of pixels. For the rest of the analysis only the data on anopheline-positive habitats were used, since some stagnant aquatic habitats might not be suitable for mosquito reproduction and our focus is on the habitats that are productive.
For the second level analysis, a Nearest Neighborhood analysis was used to assess the distance between a habitat and its first and second nearest habitat locations. To determine the climatic impact, three types of comparison were performed: (1) the spatial overlapping and nearest neighbor of habitats between the February and the May within the same year; (2) the spatial overlapping and nearest neighbor of habitats in the February between different years and (3) the spatial overlapping and nearest neighbor of habitats in the May between different years. Even though the six datasets were collected separately in six different field surveys, the consistency in the sampling was maintained through the use of the same research group with the same training and sampling methods in the same study area.
Analysis of the temporal variation in the geographic extent of larval habitats
The spatial patterns of habitats were quantified using two statistics: the dispersion pattern analysis and the compact analysis. The dispersion pattern analysis measures how dispersed habitat locations are around the geographic center of the study area. It takes into account dominant landforms (e.g., a river or cliff) and reveals global trends in a point pattern, such as orientation. The compactness analysis complements the dispersion pattern analysis by quantifying the shape of the range of mosquito larval habitats. The dispersion and compactness of habitat distribution have a direct impact on habitat management: the more dispersed and compact the habitats are in an area, the more management efforts are required for these habitats.
To examine the dispersion pattern of anopheline-positive habitats, the six sets of point data were analyzed using a Standard Deviation Ellipse (SDE) test. SDE calculates an ellipse to describe the dispersion of points, and has been widely used in the geographic analysis of point patterns [21, 22]. In this study, four statistics were derived from the ellipse: the length of the major axis (the maximum dispersion distance), the length of the minor axis (the minimum dispersion distance), the product of the lengths of the major and the minor axes (an approximation of the area of an ellipse), and the orientation of the dispersion. The orientation of an ellipse is measured as the degree by which the major axis is rotated from the geographic east in the counterclockwise direction . The Standard Deviation Ellipse (SDE) was computed using CrimeStat 3.0 .
Because habitat geographies often have irregular shapes that can have important implications for habitat management , the shape of mosquito larval habitat ranges were explicitly analyzed using convex polygon analysis  combined with the Boyce-Clark index . These two methods are widely adopted measures for habitat compactness and shape . The convex polygon analysis identifies the smallest polygon that can enclose all habitat locations. For each set of habitat locations, a convex polygon was calculated and a total of six convex polygons was obtained. The Boyce-Clark index expresses how closely the shape of the habitat distribution (represented by the convex polygon) approximates that of a hypothetical circle. It is calculated by drawing a set of equally spaced radials from a location within a shape to its perimeter and measure. The variation in lengths of the radials is then calculated. In this study, 13 radials with 30 degree space were used (30 degree is selected based on the smallest angle between any two of the neighboring vertexes in each of the six convex polygons). This index ranges from 0 to 200. A value of 0 indicates the shape is close to a perfect circle. A value of 200 indicates the shape is close to a line.
After SDE analysis, convex polygon analysis and Boyce-Clark index were applied to each of the six sets of point data, pair-wise comparisons were carried out to determine the seasonal or inter-annual changes in larval distribution patterns. The three types of comparisons outlined in the first subsection of the statistical analysis were also performed: (1) the difference between the February and the May within the same year (three pairs of datasets); (2) the difference in the February between different years (three pairs of datasets); and (3) the difference in the May between different years (three pairs of datasets).
where mg is the mean value of an environmental gradient in the study area, me is the mean value of the environmental gradient in the locations occupied by habitats, and σg is the standard deviation of the environmental gradient in the study area.
where σg is the standard deviation of an environmental gradient in the study area and σe is the standard deviation of the environmental gradient in the areas occupied by habitats.
A habitat distribution was considered specialized if it has a narrow niche breadth (its specialization value is larger than one). If a species differs considerably in its marginality and specializations along an environmental gradient, this gradient is less likely to be adequately predictive for the distribution of this species [29, 30].
Mosquito larval habitats
Numbers of stagnant aquatic habitats and anopheline-positive habitats recorded during the six field surveys from 2003 to 2005
Number of aquatic habitats
Number of anopheline-positive habitats
Percentage of anopheline-positive habitats
Average number of habitats in February
Average number of habitats in May
Total number of habitats
Seasonal and inter-annual variation in habitat locations
Nearest Neighbor analysis in anopheline-positive habitat distribution in western Kenya highlands
Second nearest neighbor
Average first and second nearest neighbor distance between the rainy and dry seasons within a year (in meters)
Average first and second nearest neighbor distance between the same seasons in different years (in meters)
Seasonal and inter-annual variation in the geographic extent of larval habitats
Changes in dispersion pattern analysis in anopheline-positive habitat distribution across seasons and years.
Ratio of long to long axis
Ratio of short to short axis
Ratio of the approximate areas of the two ellipses
Difference in the angles of the two long axes
Variation in dispersion patterns between the rainy and dry seasons within a year
Inter-annual variation in dispersion patterns
Dry season habitats were more dispersed than the rainy season habitats in the same year in 2003 and 2005, but little seasonal difference in dispersion patterns was found in 2004, the year with low precipitation in the rainy season (as shown in Figure 2). However, the differences in the orientation of the ellipses between the dry and rainy seasons were greatest in 2004, measuring 77.5 degrees (Table 3).
Changes in Boyce-Clark index for anopheline-positive habitats across seasons and years.
Difference between the Februrary and May within a year
Difference between the same seasons in different years
Seasonal and inter-annual variation in niche patterns of habitats
Changes in marginality and specialization coefficients for anopheline-positive habitats across seasons and years.
Distance to streams
Distance to streams
Difference between the rainy and dry seasons within a year
Difference between the same seasons in different years
Among these four environmental variables, the differences in S values were also the largest for distance to streams. As shown in Table 5, for distance to streams, the anopheline-positive larval habitats in the dry season were more specialized than in the rainy season of the same year in 2004. The opposite association was found for habitats in 2003 and 2005. For this variable, the seasonal difference in S values was more apparent in 2005. The inter-annual changes in S values for dry season anopheline-positive larval habitats indicates that dry season habitats were most specialized for distance to streams in 2005. The inter-annual changes in S values for rainy season larval habitats indicate that rainy season anopheline-positive larval habitats were most specialized in 2004 (As shown in Figure 2, the precipitation level in the May of 2004 is lower than such levels in the May of other years). Other environmental variables are not discussed here, since the changes in their M and S values are not significant.
In this study, we systematically investigated the temporal and spatial variability in anopheline-positive larval habitats using several statistical methods that are rarely used in mosquito studies. The results extend the existing mosquito studies by providing detailed quantitative analysis of the spatial patterns of anopheline-positive larval habitats . We identified and quantified significant seasonal and inter-annual changes in anopheline-positive larval habitat distribution patterns. The quantification of these changes could provide important information for mosquito larval habitat management. Firstly, we found that percentages of overlapping anopheline-positive habitats across seasons or years were generally below 15% and more than half of the locations occupied by dry season habitats were not occupied by habitats in the May of the same year. If a map showing the distribution of the observed habitats at a particular time is to be used in larval habitat identification in a different time, it is necessary to search the area in a radius of 190 m or 280 m of each habitat location for the nearest one or two habitats respectively. Second, during a year with an irregular amount of precipitation, a habitat map generated from a previous time point could be particularly misleading. This conclusion is drawn from the comparison between habitat patterns in the May of 2004 and 2005. As shown in Figure 2, the precipitation in the May of 2004 was lower than the precipitation in the May of 2005. Our analysis shows that the spatial patterns of rainy season habitats in these two years were very different. A more spatially dispersed and environmentally marginal distribution was recorded during the rainy season of 2004 than during other rainy seasons. In this case, the total number of habitats was not extremely sensitive to the variation in the precipitation, since we observed a comparable number of anopheline-positive habitats in 2004 and 2005. This indicates that even there are little changes in total numbers of habitats in different time periods, habitat patterns in these time periods could be different. Finally, the inter-annual changes in habitat patterns between years with similar amount of precipitations could be smaller than the seasonal changes. This conclusion is based on the results from the nearest neighbor analysis and the niche analysis. As shown in Figure 2, the precipitation levels in the May of 2003 and 2005 are similar. The nearest distance analysis and niche analysis both revealed that the inter-annual differences in the patterns of rainy season habitats between these two years are smaller than the seasonal differences in these two years. Taken together, our results suggest that caution should be practiced in extrapolating potential focal points of mosquito habitats based on data collected at different time periods. It is believed that the pattern of larval habitats is extremely important to the epidemiology of malaria because it results in spatial heterogeneity in the adult mosquito population and, subsequently, the spatial distribution of clinical malaria cases [32, 33]. This study illustrates the extreme extent of inter-annual and seasonal variations in vector habitat distributions and the necessity of continually updating habitat maps to track this highly dynamic system.
Although the topography of the study area is representative of the western Kenya highlands, the extent to which our findings can be applied to lowland areas certainly necessitates further study. It is generally acknowledged that spatial analysis is sensitive to the geographic scale of a study. The size of our study area was limited in order to minimize expenses. An important next step would be to determine the extent to which our results scale up to larger geographical regions. Finally, this study utilized mosquito habitat data collected in six different time periods with apparent seasonal and inter-annual climatic variation. Using these data, this study explored and revealed the impact of climatic variability on the distribution of mosquito larval habitats. However, to project the long term changes in patterns of mosquito habitats in accordance with long term trends of climate, time series analysis based on more temporally extensive data is needed.
One potential vector control tool which has had considerable historical success in Africa is larval habitat modification [34–36]. With the advent of geographical information systems, the use of larval habitat maps has become increasingly popular in malaria control programs . There has been a recent call for vector management strategies integrating multiple tools. Our study provides a systematic analysis of the spatial and temporal patterns of mosquito habitats with serious implications with the ability of habitat maps to predict the future distribution of habitats. Significant inter-seasonal and -annual changes in the spatial location and distribution patterns of habitats were observed at a micro-spatial scale. Our results indicate that larval source reduction should be a dynamic and adaptive process and that integrated vector control approaches must take into account the temporal variability in larval habitat locations. Obviously, continuous monitoring of habitat location can help focus the larval control in the areas where habitats actually occur. However, it is financially prohibitive for most African countries to conduct larval distribution surveys continuously in a large geographical area. A possible solution is to identify the possible causes of habitat location changes and establish dynamic models to help the habitat source reduction activities. Munga et al (in press) investigated the relationship between habitat location and landuse changes in western Kenyan using habitat and landuse data collected from 2002 to 2005. In these four years, land cover changes demonstrated a great impact on the occurrence of anopheline larval habitats. It is possible that changes on habitat locations can be modelled by a combination of relevant environmental and anthropogenic factors. Future research on such models is especially needed in the areas where continuous survey of habitats is difficult.
We thank Steve Munga, Emmanuel Mushinzimana and Andrew K. Githeko for organizing the field research and providing archived data. This work was supported by NIH grant R01 AI050243.
- Malakooti MA, Biomndo K, Shanks GD: Reemergence of epidemic malaria in the highlands of western Kenya. Emerging Infectious Diseases 1998, 4:671–676.View ArticlePubMed
- Shanks GD, Biomndo K, Hay SI, Snow RW: Changing patterns of clinical malaria since 1965 among a tea estate population located in the Kenyan highlands. Transactions of the Royal Society of Tropical Medicine and Hygiene 2000, 94:253–255.View ArticlePubMed
- Akhwale WS, Lum JK, Kaneko A, Eto H, Obonyo C, Bjorkman A, Kobayakawa T: Anemia and malaria at different altitudes in the western highlands of Kenya. Acta Tropica 2004, 91:167–175.View ArticlePubMed
- Lengeler C: Insecticide treated bednets and curtains for malaria control (Cochrane Review). The Cochrane Library Oxford 1998.
- Lengeler C: Insecticide-treated nets for malaria control: real gains: Editorials Christian Lengeler. Bulletin of the World Health Organization 2004, 82:84.PubMed
- Schellenberg JA, Minja H, Mponda H, Kikumbih N, Mushi A, Nathan R, Abdulla S, Mukasa O, Marchant TJ, Tanner M, et al.: Re-treatment of mosquito nets with insecticide. Transactions of the Royal Society of Tropical Medicine and Hygiene 2002, 96:368–369.View ArticlePubMed
- Phillips-Howard PA, Nahlen BL, Kolczak MS, Hightower AW, Kuile FOT, Alaii JA, Gimnig JE, Arudo J, Vulule JM, Odhacha A, et al.: Efficacy of permethrin-treated bed nets in the prevention of mortality in young children in an area of high perennial malaria transmission in western Kenya. The American journal of tropical medicine and hygiene 2003, 68:23–29.PubMed
- Hemingway J, Field L, Vontas J: An overview of insecticide resistance. Science 2002, 298:96.View ArticlePubMed
- Fillinger U, Ndegwa B, Githeko A, Lindsay SW: Integrated malaria vector control with microbial larvicides and insecticide treated nets in the western Kenyan highlands: a controlled trial. Bull World Health Organ 2009, 87:655–665.View ArticlePubMed
- GeissbAhler Y, Kannady K, Chaki PP, Emidi B, Govella NJ, Mayagaya V, Kiama M, Mtasiwa D, Mshinda H, Lindsay SW: Microbial larvicide application by a large-scale, community-based program reduces malaria infection prevalence in urban Dar Es Salaam, Tanzania. PLoS ONE 2009, 4:345.
- Dale PE, Ritchie SA, Territo BM, Morris CD, Muhar A, Kay BH: An overview of remote sensing and GIS for surveillance of mosquito vector habitats and risk assessment. Journal of Vector Ecology 1998, 23:54–61.PubMed
- Mushinzimana E, Munga S, Minakawa N, Li L, Feng C, Bian L, Kitron U, Schmidt C, Beck L, Zhou G: Landscape determinants and remote sensing of anopheline mosquito larval habitats in the western Kenya highlands. Malaria Journal 2006, 5:13.View ArticlePubMed
- Gu W, Novak R: Habitat-based modeling of impacts of mosquito larval interventions on entomological inoculation rates, incidence, and prevalence of malaria. The American journal of Tropical Medicine and Hygiene 2005, 73:546–552.PubMed
- Rejmankova E, Savage HM, Rejmanek M, Arredondo-Jimenez JI, Roberts DR: Multivariate Analysis of Relationships Between Habitats, Environmental Factors and Occurrence of Anopheline Mosquito Larvae Anopheles albimanus and A. pseudopunctipennis in Southern Chiapas, Mexico. Journal of Applied Ecology 1991, 28:827–841.View Article
- Beck LR, Rodriguez MH, Dister SW, Rodriguez AD, Rejmankova E, Ulloa A, Meza RA, Roberts DR, Paris JF, Spanner MA: Remote sensing as a landscape epidemiologic tool to identify villages at high risk for malaria transmission. The American Journal of Tropical Medicine and Hygiene 1994, 51:271.PubMed
- Lindsay SW, Martens WJM: Malaria in the African highlands: past, present and future. Bulletin of the World Health Organization 1998, 76:33–46.PubMed
- Munyekenye OG, Githeko AK, Zhou G, Mushinzimana E, Minakawa N: Plasmodium falciparum spatial analysis, western Kenya highlands. Emerging Infectious Diseases 2005, 11:1571–1577.PubMed
- Minakawa N, Githure JI, Beier JC, Yan G: Anopheline Mosquito Survival Strategies During the Dry Period in Western Kenya. Journal of Medical Entomology 2001, 38:388–392.View ArticlePubMed
- Bodker R, Akida J, Shayo D, Kisinza W, Msangeni HA, Pedersen EM, Lindsay SW: Relationship between altitude and intensity of malaria transmission in the Usambara Mountains, Tanzania. Journal of Medical Entomology 2003, 40:706–717.View ArticlePubMed
- Minakawa N, Seda P, Yan G: Influence of host and larval habitat distribution on the abundance of African malaria vectors in western Kenya. American Journal of Tropical Medicine and Hygiene 2002, 67:32.PubMed
- Shannon GW, Spurlock CW: Urban ecological containers, environmental risk cells, and the use of medical services. Economic Geography 1976, 52:171–180.View Article
- Sherman JE, Spencer J, Preisser JS, Gesler WM, Arcury TA: A suite of methods for representing activity space in a healthcare accessibility study. International Journal of Health Geographics 2005, 4:24.View ArticlePubMed
- Levine N: CrimeStat: A spatial statistics program for the analysis of crime incident locations (v 2.0). Ned Levine & Associates, Houston, TX, and the National Institute of Justice, Washington, DC 2002.
- Franzreb KE: Implications of home-range estimation in the management of red-cockaded woodpeckers in South Carolina. Forest Ecology and Management 2006, 228:274–284.View Article
- Mohr CO: Table of Equivalent Populations of North American Small Mammals. The American Midland Naturalist 1947, 3:223.
- Boyce RR, Clark WAV: The Concept of Shape in Geography. Geographical Review 1964, 54:561–572.View Article
- Schoener T: The contribution of ecology to an understanding of the natural world, Ecological concepts Blackwell Scientific Publications 1989 chap 8:79–113.
- Hirzel AH, Hausser J, Chessel D, Perrin N: Ecological-niche factor analysis: how to compute habitat-suitability maps without absence data? Ecology 2002, 83:2027–2036.View Article
- Grinnell J: Field Tests of Theories Concerning Distributional Control. American Naturalist 1917, 51:115.View Article
- Colwell RK, Futuyma DJ: On the measurement of niche breadth and overlap. Ecology 1971, 3:567–576.View Article
- Munga S, Mushinzimana E, Minakawa N, Zhou G, Okeyo-Owuor J, Githeko A, Yan G: Land use and Land cover changes and spatiotemporal dynamics of Anopheline Larval habitats during a period of four years in Western Kenya Highland. American Journal of Tropical Medicine and Hygiene 2009, in press.
- Brooker S, Clarke S, Njagi JK, Polack S, Mugo B, Estambale B, Muchiri E, Magnussen P, Cox J: Spatial clustering of malaria and associated risk factors during an epidemic in a highland area of western Kenya. Tropical Medicine & International Health 2004, 9:757–766.View Article
- Smith D, Dushoff J, McKenzie F: The Risk of a Mosquito-Borne Infection a Heterogeneous Environment. PLoS Biol 2004, 2:368.View Article
- Gilroy AB: Malaria control by coastal swamp drainage in West Africa Ross Institute of Tropical Hygiene 1948.
- Utzinger J, Tozan Y, Singer BH: Efficacy and cost-effectiveness of environmental management for malaria control. Tropical Medicine & International Health 2001, 6:677–687.View Article
- Coetzee M, Craig M, Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitology Today 2000, 16:74–77.View ArticlePubMed
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.